Elucidating emergence and transmission of multidrug-resistant tuberculosis in treatment experienced patients by whole genome sequencing

Define SNP thresholds for delineation of tuberculosis transmissions using whole-genome sequencing

For facilitating tuberculosis (TB) control, we used a whole-genome sequencing (WGS)-based approach to delineate transmission networks in a country with an intermediate burden of TB. A cluster was defined as Mycobacterium tuberculosis isolates with identical genotypes, and an outbreak was defined as clustered cases with epidemiological links (epi-links). To refine a cluster predefined using space oligonucleotide typing and mycobacterial interspersed repetitive unit variable tandem repeat typing, we analyzed one pansusceptible TB (C1) and three multidrug-resistant (MDR)-TB (C2-C4) clusters from different scenarios. Pansusceptible TB cluster (C1) consisting of 28 cases had ≤5 single nucleotide polymorphisms (SNPs) difference between their isolates. C1 was a definite outbreak, with cases attending the same junior high school in 2012. Three MDR-TB clusters (C2-C4) with distinct genotypes were identified, each consisting of 12-22 cases. Some of the cases had either ≤5 or ≤15 SNPs difference with clear or probable epi-links. Of note, even though WGS could effectively assist TB contact tracing, we still observed missing epi-links in some cases within the same cluster. Our results showed that thresholds of ≤5 and ≤15 SNPs difference between isolates were used to categorize definite and probable TB transmission, respectively. Furthermore, a higher SNP threshold might be required to define an MDR-TB outbreak. WGS still needs to be combined with classical epidemiological methods for improving outbreak investigations. Importantly, different SNP thresholds have to be applied to define outbreaks.

IMPORTANCE

TB is a chronic disease. Depending on host factors and TB burden, clusters of cases may continue to increase for several years. Conventional genotyping methods overestimate TB transmission, hampering precise detection of outbreaks and comprehensive surveillance. WGS can be used to obtain SNP information of M. tuberculosis to improve discriminative limitations of conventional methods and to strengthen delineation of transmission networks. It is important to define the country-specific SNP thresholds for investigation of transmission. This study demonstrated the use of thresholds of ≤5 and ≤15 SNPs difference between isolates to categorize definite and probable transmission, respectively. Different SNP thresholds should be applied while a higher cutoff was required to define an MDR-TB outbreak. The utilization of SNP thresholds proves to be crucial for guiding public health interventions, eliminating the need for unnecessary public health actions, and potentially uncovering undisclosed TB transmissions.

GenoMycAnalyzer: a web-based tool for species and drug resistance prediction for Mycobacterium genomes

BACKGROUND

Drug-resistant tuberculosis (TB) is a major threat to global public health. Whole-genome sequencing (WGS) is a useful tool for species identification and drug resistance prediction, and many clinical laboratories are transitioning to WGS as a routine diagnostic tool. However, user-friendly and high-confidence automated bioinformatics tools are needed to rapidly identify M. tuberculosis complex (MTBC) and non-tuberculous mycobacteria (NTM), detect drug resistance, and further guide treatment options.

RESULTS

We developed GenoMycAnalyzer, a web-based software that integrates functions for identifying MTBC and NTM species, lineage and spoligotype prediction, variant calling, annotation, drug-resistance determination, and data visualization. The accuracy of GenoMycAnalyzer for genotypic drug susceptibility testing (gDST) was evaluated using 5,473 MTBC isolates that underwent phenotypic DST (pDST). The GenoMycAnalyzer database was built to predict the gDST for 15 antituberculosis drugs using the World Health Organization mutational catalogue. Compared to pDST, the sensitivity of drug susceptibilities by the GenoMycAnalyzer for first-line drugs ranged from 95.9% for rifampicin (95% CI 94.8-96.7%) to 79.6% for pyrazinamide (95% CI 76.9-82.2%), whereas those for second-line drugs ranged from 98.2% for levofloxacin (95% CI 90.1-100.0%) to 74.9% for capreomycin (95% CI 69.3-80.0%). Notably, the integration of large deletions of the four resistance-conferring genes increased gDST sensitivity. The specificity of drug susceptibilities by the GenoMycAnalyzer ranged from 98.7% for amikacin (95% CI 97.8-99.3%) to 79.5% for ethionamide (95% CI 76.4-82.3%). The incorporated Kraken2 software identified 1,284 mycobacterial species with an accuracy of 98.8%. GenoMycAnalyzer also perfectly predicted lineages for 1,935 MTBC and spoligotypes for 54 MTBC.

CONCLUSIONS

GenoMycAnalyzer offers both web-based and graphical user interfaces, which can help biologists with limited access to high-performance computing systems or limited bioinformatics skills. By streamlining the interpretation of WGS data, the GenoMycAnalyzer has the potential to significantly impact TB management and contribute to global efforts to combat this infectious disease. GenoMycAnalyzer is available at http://www.mycochase.org .

Evolution of Laboratory Diagnosis of Tuberculosis

Tuberculosis (TB) is an infectious disease of global public health importance caused by the Mycobacterium tuberculosis complex. Despite advances in diagnosis and treatment, this disease has worsened with the emergence of multidrug-resistant strains of tuberculosis. We aim to present and review the history, progress, and future directions in the diagnosis of tuberculosis by evaluating the current methods of laboratory diagnosis of tuberculosis, with a special emphasis on microscopic examination and cultivation on solid and liquid media, as well as an approach to molecular assays. The microscopic method, although widely used, has its limitations, and the use and evaluation of other techniques are essential for a complete and accurate diagnosis. Bacterial cultures, both in solid and liquid media, are essential methods in the diagnosis of TB. Culture on a solid medium provides specificity and accuracy, while culture on a liquid medium brings rapidity and increased sensitivity. Molecular tests such as LPA and Xpert MTB/RIF have been found to offer significant benefits in the rapid and accurate diagnosis of TB, including drug-resistant forms. These tests allow the identification of resistance mutations and provide essential information for choosing the right treatment. We conclude that combined diagnostic methods, using several techniques and approaches, provide the best result in the laboratory diagnosis of TB. Improving the quality and accessibility of tests, as well as the implementation of advanced technologies, is essential to help improve the sensitivity, efficiency, and accuracy of TB diagnosis.

Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019

Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.

Impact of Whole-Genome Sequencing of Mycobacterium tuberculosis on Treatment Outcomes for MDR-TB/XDR-TB: A Systematic Review

The emergence and persistence of drug-resistant tuberculosis is a major threat to global public health. Our objective was to assess the applicability of whole-genome sequencing (WGS) to detect genomic markers of drug resistance and explore their association with treatment outcomes for multidrug-resistant/extensively drug-resistant tuberculosis (MDR/XDR-TB).

METHODS

Five electronic databases were searched for studies published in English from the year 2000 onward. Two reviewers independently conducted the article screening, relevant data extraction, and quality assessment. The data of the included studies were synthesized with a narrative method and are presented in a tabular format.

RESULTS

The database search identified 949 published articles and 8 studies were included. An unfavorable treatment outcome was reported for 26.6% (488/1834) of TB cases, which ranged from 9.7 to 51.3%. Death was reported in 10.5% (194/1834) of total cases. High-level fluoroquinolone resistance (due to gyrA 94AAC and 94GGC mutations) was correlated as the cause of unfavorable treatment outcomes and reported in three studies. Other drug resistance mutations, like kanamycin high-level resistance mutations (rrs 1401G), rpoB Ile491Phe, and ethA mutations, conferring prothionamide resistance were also reported. The secondary findings from this systematic review involved laboratory aspects of WGS, including correlations with phenotypic DST, cost, and turnaround time, or the impact of WGS results on public health actions, such as determining transmission events within outbreaks.

CONCLUSIONS

WGS has a significant capacity to provide accurate and comprehensive drug resistance data for MDR/XDR-TB, which can inform personalized drug therapy to optimize treatment outcomes.

pncA Large Deletion is the Characteristic of Pyrazinamide-Resistant Mycobacterium tuberculosis belonging to the East Asian Lineage

BACKGROUND

Pyrazinamide (PZA) is often used as an add-on agent in the treatment of multidrug-resistant tuberculosis, regardless of phenotypic drug susceptibility testing (pDST) results. However, evaluating the effectiveness of PZA is challenging because of its low pH activity, which can result in unreliable pDST results. This study aimed to investigate the genomic characteristics associated with PZA resistance that can be used to develop genotypic DST.

MATERIALS AND METHODS

A publicly available whole genome sequencing (WGS) dataset of 10,725 Mycobacterium tuberculosis complex genomes (3,326 phenotypically PZA-resistant and 7,399 phenotypically PZA-susceptible isolates) were analyzed.

RESULTS

In total, 2,934 pncA non-silent mutations were identified in 2,880 isolates (26.9%). Detected mutations were found throughout the entire coding region of pncA in a scattered pattern, of which the most frequent mutation was p.Q10P (n = 278), followed by p.H57D (n = 167) and c.-11A>G (n = 122). The sensitivity and specificity of the group 1 or 2 mutations reported by the World Health Organization (WHO) mutational catalogue were 73.0% and 98.9%, respectively. We further identified 18 novel pncA mutations that were significantly associated with phenotypically PZA-resistant. In addition to these mutations, we identified 102 large deletions in the pncA gene, and all but two isolates were phenotypically resistant to PZA isolates. Notably, pncA deletions were mutually exclusive to pncA mutations, and more than half of the isolates with pncA large deletions belonged to the East Asian lineage (67.6%). The sensitivity, specificity, positive predictive value, and negative predictive value of the pooled variants (group 1 or 2 mutations, novel resistance-associated mutations, and large deletions of the pncA gene) were 79.0%, 98.9%, 97.0%, and 91.3%, respectively. The area under the curve (AUC) value for the pooled variants was significantly higher than the AUC value for the group 1 or 2 mutations (P <0.001), indicating that the pooled variants have a better discriminative ability for predicting PZA resistance.

CONCLUSION

Using WGS, we found that the pncA mutations are scattered without specific mutational hotspots, and large deletions associated with PZA resistance are more common in the East Asian lineage of M. tuberculosis isolates. Our data also demonstrated the reliability of group 1 or 2 mutations presented in the WHO mutation catalogue and the need for further investigation on group 3 mutations, contributing to the evaluation of the current knowledge base on mutations associated with the PZA-resistant M. tuberculosis complex.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

GWAS and functional studies suggest a role for altered DNA repair in the evolution of drug resistance in Mycobacterium tuberculosis

The emergence of drug resistance in Mycobacterium tuberculosis (Mtb) is alarming and demands in-depth knowledge for timely diagnosis. We performed genome-wide association analysis using 2237 clinical strains of Mtb to identify novel genetic factors that evoke drug resistance. In addition to the known direct targets, we identified for the first time, a strong association between mutations in DNA repair genes and the multidrug-resistant phenotype. To evaluate the impact of variants identified in the clinical samples in the evolution of drug resistance, we utilized knockouts and complemented strains in Mycobacterium smegmatis and Mtb. Results show that variant mutations compromised the functions of MutY and UvrB. MutY variant showed enhanced survival compared with wild-type (Rv) when the Mtb strains were subjected to multiple rounds of ex vivo antibiotic stress. In an in vivo guinea pig infection model, the MutY variant outcompeted the wild-type strain. We show that novel variant mutations in the DNA repair genes collectively compromise their functions and contribute to better survival under antibiotic/host stress conditions.

A convolutional neural network highlights mutations relevant to antimicrobial resistance in Mycobacterium tuberculosis

Long diagnostic wait times hinder international efforts to address antibiotic resistance in M. tuberculosis. Pathogen whole genome sequencing, coupled with statistical and machine learning models, offers a promising solution. However, generalizability and clinical adoption have been limited by a lack of interpretability, especially in deep learning methods. Here, we present two deep convolutional neural networks that predict antibiotic resistance phenotypes of M. tuberculosis isolates: a multi-drug CNN (MD-CNN), that predicts resistance to 13 antibiotics based on 18 genomic loci, with AUCs 82.6-99.5% and higher sensitivity than state-of-the-art methods; and a set of 13 single-drug CNNs (SD-CNN) with AUCs 80.1-97.1% and higher specificity than the previous state-of-the-art. Using saliency methods to evaluate the contribution of input sequence features to the SD-CNN predictions, we identify 18 sites in the genome not previously associated with resistance. The CNN models permit functional variant discovery, biologically meaningful interpretation, and clinical applicability.

Whole-Genome Sequencing for Resistance Level Prediction in Multidrug-Resistant Tuberculosis

Defining the precise relationship between resistance mutations and quantitative phenotypic drug susceptibility testing will increase the value of whole-genome sequencing (WGS) for predicting tuberculosis drug resistance. However, a large number of WGS data sets currently lack corresponding quantitative phenotypic data—the MICs. Using MYCOTBI plates, we determined the MICs to nine antituberculosis drugs for 154 clinical multidrug-resistant tuberculosis isolates from the Shenzhen Center for Chronic Disease Control in Shenzhen, China. Comparing MICs with predicted drug-resistance profiles inferred by WGS showed that WGS could predict the levels of resistance to isoniazid, rifampicin, streptomycin, fluoroquinolones, and aminoglycosides. We also found some mutations that may not be associated with drug resistance, such as EmbB D328G, mutations in the gid gene, and C−12T in the eis promoter. However, some strains carrying the same mutations showed different levels of resistance to the corresponding drugs. The MICs of different strains with the RpsL K88R, fabG1 C−15T mutations and some with mutations in embB and rpoB, had MICs to the corresponding drugs that varied by 8-fold or more. This variation is unexplained but could be influenced by the bacterial genetic background. Additionally, we found that 32.3% of rifampicin-resistant isolates were rifabutin-susceptible, particularly those with rpoB mutations H445D, H445L, H445S, D435V, D435F, L452P, S441Q, and S441V. Studying the influence of bacterial genetic background on the MIC and the relationship between rifampicin-resistant mutations and rifabutin resistance levels should improve the ability of WGS to guide the selection of medical treatment regimens.

IMPORTANCE Whole-genome sequencing (WGS) has excellent potential in drug-resistance prediction. The MICs are essential indications of adding a particular antituberculosis drug dosage or changing the entire treatment regimen. However, the relationship between many known drug-resistant mutations and MICs is unclear, especially for rarer ones. The results showed that WGS could predict resistance levels to isoniazid, rifampicin, streptomycin, fluoroquinolones, and aminoglycosides. However, some mutations may not be associated with drug resistance, and some others may confer various MICs to strains carrying them. Also, 32.3% of rifampicin (RIF)-resistant strains were classified as sensitive to rifabutin (RFB), and some mutations in the rpoB gene may be associated with this phenotype. Our data on the MIC distribution of strains with some rarer mutations add to the accumulated data on the resistance level associated with such mutations to help guide further research and draw meaningful conclusions.

A Whole-Genome Sequencing-Based Study to Delineate the Risk and Characteristics of Tuberculosis Transmission in an Insular Population Over 10 Years in Shanghai

Background

Tuberculosis (TB) has remained a tough problem in China. This study aims to identify the risk of tuberculosis transmission and to assess its characteristics.

Methods

We performed a molecular epidemiological study for patients with culture-positive Mycobacterium tuberculosis (M. tuberculosis) in Shanghai, from 2009 to 2018. Demographic information was obtained from the Tuberculosis Information Management System. Whole-genome sequencing (WGS) was conducted with a threshold of 12 single-nucleotide polymorphisms (SNPs) to distinguish the genomic cluster. To analyze the characteristics of TB transmission, the contact investigation for clustered cases was performed.

Results

In total, 94 (27.25%) of the 345 enrolled patients were grouped into 42 genomic clusters, indicating local transmission of M. tuberculosis strains. Compared to a health system delay <14 days, patients with a health system delay ≥14 days [adjusted odds ratios (AOR) = 2.57, 95% confidence interval (CI): 1.34–4.95] were more likely to be clustered. Patients under 65 years old (AOR = 3.11, 95% CI: 1.76–5.49), residents (AOR = 2.43, 95% CI: 1.18–4.99), and Beijing genotype strains (AOR = 3.35, 95% CI: 1.32–8.53) were associated with increased risk of clustering. Interestingly, patients with resistance to isoniazid (AOR = 2.36, 95% CI: 1.15–4.88) had a higher risk of transmission. Sixteen confirmed/probable epidemiological links were identified. Local transmission of imported cases and household transmission were prominent.

Conclusion

Health system delay is a crucial factor for TB transmission. Patients with resistance to isoniazid should be priority targets for contact investigation to reduce transmission.

Potential contribution of HIV during first-line tuberculosis treatment to subsequent rifampicin-monoresistant tuberculosis and acquired tuberculosis drug resistance in South Africa: a retrospective molecular epidemiology study

Background

South Africa has a high burden of rifampicin-resistant tuberculosis (including multidrug-resistant [MDR] tuberculosis), with increasing rifampicin-monoresistant (RMR) tuberculosis over time. Resistance acquisition during first-line tuberculosis treatment could be a key contributor to this burden, and HIV might increase the risk of acquiring rifampicin resistance. We assessed whether HIV during previous treatment was associated with RMR tuberculosis and resistance acquisition among a retrospective cohort of patients with MDR or rifampicin-resistant tuberculosis.

Methods

In this retrospective cohort study, we included all patients routinely diagnosed with MDR or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa, between Jan 1, 2008, and Dec 31, 2017. Patient-level data were obtained from a prospective database, complemented by data on previous tuberculosis treatment and HIV from a provincial health data exchange. Stored MDR or rifampicin-resistant tuberculosis isolates from patients underwent whole-genome sequencing (WGS). WGS data were used to infer resistance acquisition versus transmission, by identifying genomically unique isolates (single nucleotide polymorphism threshold of five). Logistic regression analyses were used to assess factors associated with RMR tuberculosis and genomic uniqueness.

Findings

The cohort included 2041 patients diagnosed with MDR or rifampicin-resistant tuberculosis between Jan 1, 2008, and Dec 31, 2017; of those, 463 (22·7%) with RMR tuberculosis and 1354 (66·3%) with previous tuberculosis treatment. In previously treated patients, HIV positivity during previous tuberculosis treatment versus HIV negativity (adjusted odds ratio [OR] 2·07, 95% CI 1·35–3·18), and three or more previous tuberculosis treatment episodes versus one (1·96, 1·21–3·17) were associated with RMR tuberculosis. WGS data showing MDR or rifampicin-resistant tuberculosis were available for 1169 patients; 360 (30·8%) isolates were identified as unique. In previously treated patients, RMR tuberculosis versus MDR tuberculosis (adjusted OR 4·96, 3·40–7·23), HIV positivity during previous tuberculosis treatment (1·71, 1·03–2·84), and diagnosis in 2013–17 (1·42, 1·02–1·99) versus 2008–12, were associated with uniqueness. In previously treated patients with RMR tuberculosis, HIV positivity during previous treatment (adjusted OR 5·13, 1·61–16·32) was associated with uniqueness as was female sex (2·50 [1·18–5·26]).

Interpretation

These data suggest that HIV contributes to rifampicin-resistance acquisition during first-line tuberculosis treatment and that this might be driving increasing RMR tuberculosis over time. Large-scale prospective cohort studies are required to further quantify this risk.

Funding

Swiss National Science Foundation, South African National Research Foundation, and Wellcome Trust.

Whole-genome sequencing as a tool for studying the microevolution of drug-resistant serial Mycobacterium tuberculosis isolates

Treatment of drug-resistant tuberculosis requires extended use of more toxic and less effective drugs and may result in retreatment cases due to failure, abandonment or disease recurrence. It is therefore important to understand the evolutionary process of drug resistance in Mycobacterium tuberculosis. We here in describe the microevolution of drug resistance in serial isolates from six previously treated patients. Drug resistance was initially investigated through phenotypic methods, followed by genotypic approaches. The use of whole-genome sequencing allowed the identification of mutations in the katG, rpsL and rpoB genes associated with drug resistance, including the detection of rare mutations in katG and mixed populations of strains. Molecular docking simulation studies of the impact of observed mutations on isoniazid binding were also performed. Whole-genome sequencing detected 266 single nucleotide polymorphisms between two isolates obtained from one patient, suggesting a case of exogenous reinfection. In conclusion, sequencing technologies can detect rare mutations related to drug resistance, identify subpopulations of resistant strains, and identify diverse populations of strains due to exogenous reinfection, thus improving tuberculosis control by guiding early implementation of appropriate clinical and therapeutic interventions.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

Common Adaptive Strategies Underlie Within-Host Evolution of Bacterial Pathogens

Within-host adaptation is a hallmark of chronic bacterial infections, involving substantial genomic changes. Recent large-scale genomic data from prolonged infections allow the examination of adaptive strategies employed by different pathogens and open the door to investigate whether they converge toward similar strategies. Here, we compiled extensive data of whole-genome sequences of bacterial isolates belonging to miscellaneous species sampled at sequential time points during clinical infections. Analysis of these data revealed that different species share some common adaptive strategies, achieved by mutating various genes. Although the same genes were often mutated in several strains within a species, different genes related to the same pathway, structure, or function were changed in other species utilizing the same adaptive strategy (e.g., mutating flagellar genes). Strategies exploited by various bacterial species were often predicted to be driven by the host immune system, a powerful selective pressure that is not species specific. Remarkably, we find adaptive strategies identified previously within single species to be ubiquitous. Two striking examples are shifts from siderophore-based to heme-based iron scavenging (previously shown for Pseudomonas aeruginosa) and changes in glycerol-phosphate metabolism (previously shown to decrease sensitivity to antibiotics in Mycobacterium tuberculosis). Virulence factors were often adaptively affected in different species, indicating shifts from acute to chronic virulence and virulence attenuation during infection. Our study presents a global view on common within-host adaptive strategies employed by different bacterial species and provides a rich resource for further studying these processes.

Emergence of additional drug resistance during treatment of multidrug-resistant tuberculosis in China: a prospective cohort study

OBJECTIVES

Little is known about how additional second-line drug resistance emerges during multidrug-resistant tuberculosis (MDR-TB) treatment. The present study aimed to investigate the influence of microevolution, exogenous reinfection and mixed infection on second-line drug resistance during the recommended 2-year MDR-TB treatment.

METHODS

Individuals with MDR-TB were enrolled between 2013 and 2016 in a multicentre prospective observational cohort study and were followed up for 2 years until treatment completion. Whole-genome sequencing (WGS) was applied for serial Mycobacterium tuberculosis isolates from study participants throughout the treatment, to study the role of microevolution, exogenous reinfection and mixed infection in the development of second-line drug resistance.

RESULTS

Of the 286 enrolled patients with MDR-TB, 63 (22.0%) M. tuberculosis isolates developed additional drug resistance during the MDR-TB treatment, including 5 that fulfilled the criteria of extensively drug-resistant TB. By comparing WGS data of serial isolates retrieved from the patients throughout treatment, 41 (65.1%) of the cases of additional second-line drug resistance were the result of exogenous reinfection, 18 (28.6%) were caused by acquired drug resistance, i.e. microevolution, while the remaining 4 (6.3%) were caused by mixed infections with drug-resistant and drug-susceptible strains. In multivariate analysis, previous TB treatment (adjusted hazard ratio (aHR) 2.51, 95% CI 1.51-4.18), extensive disease on chest X-ray (aHR 3.39, 95% CI 2.03-5.66) and type 2 diabetes mellitus (aHR 4.00, 95% CI 2.22-7.21) were independent risk factors associated with the development of additional second-line drug resistance.

CONCLUSIONS

A large proportion of additional second-line drug resistance emerging during MDR-TB treatment was attributed to exogenous reinfection, indicating the urgency of infection control in health facilities as well as the need for repeated drug susceptibility testing throughout MDR-TB treatment.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Approaches for characterizing and tracking hospital-associated multidrug-resistant bacteria

Hospital-associated infections are a major concern for global public health. Infections with antibiotic-resistant pathogens can cause empiric treatment failure, and for infections with multidrug-resistant bacteria which can overcome antibiotics of "last resort" there exists no alternative treatments. Despite extensive sanitization protocols, the hospital environment is a potent reservoir and vector of antibiotic-resistant organisms. Pathogens can persist on hospital surfaces and plumbing for months to years, acquire new antibiotic resistance genes by horizontal gene transfer, and initiate outbreaks of hospital-associated infections by spreading to patients via healthcare workers and visitors. Advancements in next-generation sequencing of bacterial genomes and metagenomes have expanded our ability to (1) identify species and track distinct strains, (2) comprehensively profile antibiotic resistance genes, and (3) resolve the mobile elements that facilitate intra- and intercellular gene transfer. This information can, in turn, be used to characterize the population dynamics of hospital-associated microbiota, track outbreaks to their environmental reservoirs, and inform future interventions. This review provides a detailed overview of the approaches and bioinformatic tools available to study isolates and metagenomes of hospital-associated bacteria, and their multi-layered networks of transmission.

A tuberculin skin test survey among healthcare workers in two public tertiary care hospitals in Bangladesh

In Bangladesh, there is currently no data on the burden of latent TB infection (LTBI) amongst hospital healthcare workers (HCWs). This study aimed to determine the prevalence of LTBI and compare the prevalence among HCWs in two public tertiary care hospitals. Between September 2018 and August 2019, we conducted a cross-sectional study in two public tertiary care general hospitals. Using a survey and tuberculin skin test (TST), we assessed risk factors for LTBI, adjusting for known and plausible confounders. In addition, a facility assessment was undertaken to understand the implementation of relevant IPC measures. The prevalence of LTBI among HCWs was 42%. HCWs spent a median of 6 hours (SD = 1.76, IQR 2.00) per day and attended an average of 1.87 pulmonary TB patients per week. HCWs did not receive any TB IPC training, the wards lacked a symptom checklist to screen patients for TB, and no masks were available for coughing patients. Seventy-seven percent reportedly did not use any facial protection (masks or respirators) while caring for patients. In the multivariable model adjusting for hospital level clustering effect, TST positivity was significantly higher among HCWs aged 35-45 years (aOR1.36, 95% CI: 1.06-1.73) and with >3 years of service (aOR 1.67, 95% CI: 1.62-1.72). HCWs working in the medicine ward had 3.65 (95% CI: 2.20-6.05) times, and HCWs in the gynecology and obstetrics ward had 2.46 (95% CI: 1.42-4.27) times higher odds of TST positivity compared to HCWs working in administrative areas. This study identified high prevalence of LTBI among HCWs. This may be due to the level of exposure to pulmonary TB patients, and/or limited use of personal protective equipment along with poor implementation of TB IPC in the hospitals. Considering the high prevalence of LTBI, we recommend the national TB program consider providing preventative therapy to the HCWs as the high-risk group, and implement TB IPC in the hospitals.

Mycobacterium bovis: From Genotyping to Genome Sequencing

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

Differential Expression of Resistant and Efflux Pump Genes in MDR-TB Isolates

Background:

Numerous investigations demonstrate efflux as a worldwide bacterial mode of action which contributes to the resistance of drugs. The activity of antibiotics, which subjects to efflux, can be improved by the combined usage of efflux inhibitors. However, the efflux role to the overall levels of antibiotic resistance of clinical M. tuberculosis isolates is inadequately comprehended and is still disregarded by many.

Method:

Here, we assessed the contribution of resistant genes associated with isoniazid (INH) and rifampin (R) resistance to the levels of drug resistance in the (27) clinical isolates of MDR-TB. Additionally, the role of the resistance for six putative drug efflux pump genes to the antibiotics was investigated. The level of katG expression was down-regulated in 24/27 (88.88%) of MDR-TB isolates. Of the 27 MDR-TB isolates, inhA, oxyR-ahpC, and rpoB showed either overexpression or up-regulation in 8 (29.62%), 4 (14.81 %), and 24 (88.88%), respectively. Moreover, the efflux pump genes drrA, drrB, efpA, Rv2459, Rv1634, and Rv1250 were overexpressed under INH/RIF plus fresh pomegranate juice (FPJ) stress signifying the efflux pumps contribution to the overall levels of the resistance of MDR-TB isolates.

Conclusion:

These results displayed that the levels of drug resistance of MDR-TB clinical isolates are due to combination among drug efflux pump and the presence of mutations in target genes, a truth which is often ignored by the specialists of tuberculosis in favour of the almost undoubted significance of drug target- gene mutations for the resistance in M. tuberculosis.

QuantTB - a method to classify mixed Mycobacterium tuberculosis infections within whole genome sequencing data

BACKGROUND

Mixed infections of Mycobacterium tuberculosis and antibiotic heteroresistance continue to complicate tuberculosis (TB) diagnosis and treatment. Detection of mixed infections has been limited to molecular genotyping techniques, which lack the sensitivity and resolution to accurately estimate the multiplicity of TB infections. In contrast, whole genome sequencing offers sensitive views of the genetic differences between strains of M. tuberculosis within a sample. Although metagenomic tools exist to classify strains in a metagenomic sample, most tools have been developed for more divergent species, and therefore cannot provide the sensitivity required to disentangle strains within closely related bacterial species such as M. tuberculosis. Here we present QuantTB, a method to identify and quantify individual M. tuberculosis strains in whole genome sequencing data. QuantTB uses SNP markers to determine the combination of strains that best explain the allelic variation observed in a sample. QuantTB outputs a list of identified strains, their corresponding relative abundances, and a list of drugs for which resistance-conferring mutations (or heteroresistance) have been predicted within the sample.

RESULTS

We show that QuantTB has a high degree of resolution and is capable of differentiating communities differing by less than 25 SNPs and identifying strains down to 1× coverage. Using simulated data, we found QuantTB outperformed other metagenomic strain identification tools at detecting strains and quantifying strain multiplicity. In a real-world scenario, using a dataset of 50 paired clinical isolates from a study of patients with either reinfections or relapses, we found that QuantTB could detect mixed infections and reinfections at rates concordant with a manually curated approach.

CONCLUSION

QuantTB can determine infection multiplicity, identify hetero-resistance patterns, enable differentiation between relapse and re-infection, and clarify transmission events across seemingly unrelated patients - even in low-coverage (1×) samples. QuantTB outperforms existing tools and promises to serve as a valuable resource for both clinicians and researchers working with clinical TB samples.

MycoResistance: a curated resource of drug resistance molecules in Mycobacteria

The emergence and spread of drug-resistant Mycobacterium tuberculosis is of global concern. To improve the understanding of drug resistance in Mycobacteria, numerous studies have been performed to discover diagnostic markers and genetic determinants associated with resistance to anti-tuberculosis drug. However, the related information is scattered in a massive body of literature, which is inconvenient for researchers to investigate the molecular mechanism of drug resistance. Therefore, we manually collected 1707 curated associations between 73 compounds and 132 molecules (including coding genes and non-coding RNAs) in 6 mycobacterial species from 465 studies. The experimental details of molecular epidemiology and mechanism exploration research were also summarized and recorded in our work. In addition, multidrug resistance and extensively drug resistance molecules were also extracted to interpret the molecular mechanisms that are responsible for cross resistance among anti-tuberculosis drugs. Finally, we constructed an omnibus repository named MycoResistance, a user friendly interface to conveniently browse, search and download all related entries. We hope that this elaborate database will serve as a beneficial resource for mechanism explanations, precise diagnosis and effective treatment of drug-resistant mycobacterial strains.

Dispersal of Mycobacterium tuberculosis Driven by Historical European Trade in the South Pacific

Mycobacterium tuberculosis (Mtb) is a globally distributed bacterial pathogen whose population structure has largely been shaped by the activities of its obligate human host. Oceania was the last major global region to be reached by Europeans and is the last region for which the dispersal and evolution of Mtb remains largely unexplored. Here, we investigated the evolutionary history of the Euro-American L4.4 sublineage and its dispersal to the South Pacific. Using a phylodynamics approach and a dataset of 236 global Mtb L4.4 genomes we have traced the origins and dispersal of L4.4 strains to New Zealand. These strains are predominantly found in indigenous Māori and Pacific people and we identify a clade of European, likely French, origin that is prevalent in indigenous populations in both New Zealand and Canada. Molecular dating suggests the expansion of European trade networks in the early 19th century drove the dispersal of this clade to the South Pacific. We also identify historical and social factors within the region that have contributed to the local spread and expansion of these strains, including recent Pacific migrations to New Zealand and the rapid urbanization of Māori in the 20th century. Our results offer new insight into the expansion and dispersal of Mtb in the South Pacific and provide a striking example of the role of historical European migrations in the global dispersal of Mtb.

Beyond the SNP Threshold: Identifying Outbreak Clusters Using Inferred Transmissions

Whole-genome sequencing (WGS) is increasingly used to aid the understanding of pathogen transmission. A first step in analyzing WGS data is usually to define "transmission clusters," sets of cases that are potentially linked by direct transmission. This is often done by including two cases in the same cluster if they are separated by fewer single-nucleotide polymorphisms (SNPs) than a specified threshold. However, there is little agreement as to what an appropriate threshold should be. We propose a probabilistic alternative, suggesting that the key inferential target for transmission clusters is the number of transmissions separating cases. We characterize this by combining the number of SNP differences and the length of time over which those differences have accumulated, using information about case timing, molecular clock, and transmission processes. Our framework has the advantage of allowing for variable mutation rates across the genome and can incorporate other epidemiological data. We use two tuberculosis studies to illustrate the impact of our approach: with British Columbia data by using spatial divisions; with Republic of Moldova data by incorporating antibiotic resistance. Simulation results indicate that our transmission-based method is better in identifying direct transmissions than a SNP threshold, with dissimilarity between clusterings of on average 0.27 bits compared with 0.37 bits for the SNP-threshold method and 0.84 bits for randomly permuted data. These results show that it is likely to outperform the SNP-threshold method where clock rates are variable and sample collection times are spread out. We implement the method in the R package transcluster.

Extensive global movement of multidrug-resistantM. tuberculosis strains revealed by whole-genome analysis

Background

While the international spread of multidrug-resistant (MDR) Mycobacterium tuberculosis strains is an acknowledged public health threat, a broad and more comprehensive examination of the global spread of MDR-tuberculosis (TB) using whole-genome sequencing has not yet been performed.

Methods

In a global dataset of 5310 M. tuberculosis whole-genome sequences isolated from five continents, we performed a phylogenetic analysis to identify and characterise clades of MDR-TB with respect to geographic dispersion.

Results

Extensive international dissemination of MDR-TB was observed, with identification of 32 migrant MDR-TB clades with descendants isolated in 17 unique countries. Relatively recent movement of strains from both Beijing and non-Beijing lineages indicated successful global spread of varied genetic backgrounds. Migrant MDR-TB clade members shared relatively recent common ancestry, with a median estimate of divergence of 13–27 years. Migrant extensively drug-resistant (XDR)-TB clades were not observed, although development of XDR-TB within migratory MDR-TB clades was common.

Conclusions

Application of genomic techniques to investigate global MDR migration patterns revealed extensive global spread of MDR clades between countries of varying TB burden. Further expansion of genomic studies to incorporate isolates from diverse global settings into a single analysis, as well as data sharing platforms that facilitate genomic data sharing across country lines, may allow for future epidemiological analyses to monitor for international transmission of MDR-TB. In addition, efforts to perform routine whole-genome sequencing on all newly identified M. tuberculosis, like in England, will serve to better our understanding of the transmission dynamics of MDR-TB globally.

Whole genome sequencing identifies bacterial factors affecting transmission of multidrug-resistant tuberculosis in a high-prevalence setting

Whole genome sequencing (WGS) can elucidate Mycobacterium tuberculosis (Mtb) transmission patterns but more data is needed to guide its use in high-burden settings. In a household-based TB transmissibility study in Peru, we identified a large MIRU-VNTR Mtb cluster (148 isolates) with a range of resistance phenotypes, and studied host and bacterial factors contributing to its spread. WGS was performed on 61 of the 148 isolates. We compared transmission link inference using epidemiological or genomic data and estimated the dates of emergence of the cluster and antimicrobial drug resistance (DR) acquisition events by generating a time-calibrated phylogeny. Using a set of 12,032 public Mtb genomes, we determined bacterial factors characterizing this cluster and under positive selection in other Mtb lineages. Four of the 61 isolates were distantly related and the remaining 57 isolates diverged ca. 1968 (95%HPD: 1945-1985). Isoniazid resistance arose once and rifampin resistance emerged subsequently at least three times. Emergence of other DR types occurred as recently as within the last year of sampling. We identified five cluster-defining SNPs potentially contributing to transmissibility. In conclusion, clusters (as defined by MIRU-VNTR typing) may be circulating for decades in a high-burden setting. WGS allows for an enhanced understanding of transmission, drug resistance, and bacterial fitness factors.

Genome-wide analysis of Mycobacterium tuberculosis polymorphisms reveals lineage-specific associations with drug resistance

Background

Continuing evolution of the Mycobacterium tuberculosis (Mtb) complex genomes associated with resistance to anti-tuberculosis drugs is threatening tuberculosis disease control efforts. Both multi- and extensively drug resistant Mtb (MDR and XDR, respectively) are increasing in prevalence, but the full set of Mtb genes involved are not known. There is a need for increased sensitivity of genome-wide approaches in order to elucidate the genetic basis of anti-microbial drug resistance and gain a more detailed understanding of Mtb genome evolution in a context of widespread antimicrobial therapy. Population structure within the Mtb complex, due to clonal expansion, lack of lateral gene transfer and low levels of recombination between lineages, may be reducing statistical power to detect drug resistance associated variants.

Results

To investigate the effect of lineage-specific effects on the identification of drug resistance associations, we applied genome-wide association study (GWAS) and convergence-based (PhyC) methods to multiple drug resistance phenotypes of a global dataset of Mtb lineages 2 and 4, using both lineage-wise and combined approaches. We identify both well-established drug resistance variants and novel associations; uniquely identifying associations for both lineage-specific and -combined GWAS analyses. We report 17 potential novel associations between antimicrobial resistance phenotypes and Mtb genomic variants.

Conclusions

For GWAS, both lineage-specific and -combined analyses are useful, whereas PhyC may perform better in contexts of greater diversity. Unique associations with XDR in lineage-specific analyses provide evidence of diverging evolutionary trajectories between lineages 2 and 4 in response to antimicrobial drug therapy.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5615-3) contains supplementary material, which is available to authorized users.

Deciphering Within-Host Microevolution of Mycobacterium tuberculosis through Whole-Genome Sequencing: the Phenotypic Impact and Way Forward

The Mycobacterium tuberculosis genome is more heterogenous and less genetically stable within the host than previously thought. Currently, only limited data exist on the within-host microevolution, diversity, and genetic stability of M. tuberculosis As a direct consequence, our ability to infer M. tuberculosis transmission chains and to understand the full complexity of drug resistance profiles in individual patients is limited. Furthermore, apart from the acquisition of certain drug resistance-conferring mutations, our knowledge on the function of genetic variants that emerge within a host and their phenotypic impact remains scarce. We performed a systematic literature review of whole-genome sequencing studies of serial and parallel isolates to summarize the knowledge on genetic diversity and within-host microevolution of M. tuberculosis We identified genomic loci of within-host emerged variants found across multiple studies and determined their functional relevance. We discuss important remaining knowledge gaps and finally make suggestions on the way forward.

Few eligible for the newly recommended short course MDR-TB regimen at a large Mumbai private clinic

BACKGROUND

India has the world's highest tuberculosis burden, and Mumbai is particularly affected by multidrug resistant tuberculosis (MDR-TB). WHO recommends short, intensive treatment ("Short Course") for previously untreated pulmonary MDR-TB patients but does not require universal drug susceptibility testing (DST) before Short Course. DST would likely screen out many MDR-TB patients in places like Mumbai with significant drug resistance.

METHODS

MDR-TB patients at a private clinic were recruited for a prospective observational cohort. Short Course eligibility was evaluated by clinical criteria and DST results. Eligibility by DST was classified as rifampin monoresistance (as tested by Xpert MTB/RIF), rifampin, fluoroquinolones, and 2nd-line injectable drugs resistance (as tested by line probe assays) and resistance to other drugs.

RESULTS

Of 559 participants with MDR-TB, 33% met clinical eligibility for Short Course. DST for rifampin, fluoroquinolones, and 2nd-line injectable drugs excluded 74.7% of participants. Complete phenotypic DST excluded 96.6% of participants. Prior treatment with either 1st or 2nd-line drugs did not significantly affect eligibility.

CONCLUSIONS

In a global MDR-TB hotspot, < 5% of participants with MDR-TB were appropriate for Short Course by clinical characteristics and DST results. Rapid molecular testing would not sufficiently identify drug resistance in this population. Eligibility rates were not significantly reduced by prior TB treatment.

The use of whole-genome sequencing in cluster investigation of a multidrug-resistant tuberculosis outbreak

We used whole-genome sequencing (WGS) to delineate transmission networks and investigate the benefits of WGS during cluster investigation.We included clustered cases of multidrug-resistant (MDR) tuberculosis (TB)/extensively drug-resistant (XDR) TB linked by mycobacterial interspersed repetitive unit variable tandem repeat (MIRU-VNTR) strain typing or epidemiological information in the national cluster B1006, notified between 2007 and 2013 in the UK. We excluded from further investigation cases whose isolates differed by greater than 12 single nucleotide polymorphisms (SNPs). Data relating to patients' social networks were collected.27 cases were investigated and 22 had WGS, eight of which (36%) were excluded as their isolates differed by more than 12 SNPs to other cases. 18 cases were ruled into the transmission network based on genomic and epidemiological information. Evidence of transmission was inconclusive in seven out of 18 cases (39%) in the transmission network following WGS and epidemiological investigation.This investigation of a drug-resistant TB cluster illustrates the opportunities and limitations of WGS in understanding transmission in a setting with a high proportion of migrant cases. The use of WGS should be combined with classical epidemiological methods. However, not every cluster will be solvable, regardless of the quality of genomic data.

Prevalence and patterns of rifampicin and isoniazid resistance conferring mutations in Mycobacterium tuberculosis isolates from Uganda

BACKGROUND

Accurate diagnosis of tuberculosis, especially by using rapid molecular assays, can reduce transmission of drug resistant tuberculosis in communities. However, the frequency of resistance conferring mutations varies with geographic location of Mycobacterium tuberculosis, and this affects the efficiency of rapid molecular assays in detecting resistance. This has created need for characterizing drug resistant isolates from different settings to investigate frequencies of resistance conferring mutations. Here, we describe the prevalence and patterns of rifampicin- and isoniazid- resistance conferring mutations in isolates from Uganda, which could be useful in the management of MDR-TB patients in Uganda and other countries in sub-Saharan Africa.

RESULTS

Ninety seven M. tuberculosis isolates were characterized, of which 38 were MDR, seven rifampicin-resistant, 12 isoniazid-mono-resistant, and 40 susceptible to rifampicin and isoniazid. Sequence analysis of the rpoB rifampicin-resistance determining region (rpoB/RRDR) revealed mutations in six codons: 588, 531, 526, 516, 513, and 511, of which Ser531Leu was the most frequent (40%, 18/45). Overall, the three mutations (Ser531Leu, His526Tyr, Asp516Tyr) frequently associated with rifampicin-resistance occurred in 76% of the rifampicin resistant isolates while 18% (8/45) of the rifampicin-resistant isolates lacked mutations in rpoB/RRDR. Furthermore, sequence analysis of katG and inhA gene promoter revealed mainly the Ser315Thr (76%, 38/50) and C(-15)T (8%, 4/50) mutations, respectively. These two mutations combined, which are frequently associated with isoniazid-resistance, occurred in 88% of the isoniazid resistant isolates. However, 20% (10/50) of the isoniazid-resistant isolates lacked mutations both in katG and inhA gene promoter. The sensitivity of sequence analysis of rpoB/RRDR for rifampicin-resistance via detection of high confidence mutations (Ser531Leu, His526Tyr, Asp516Tyr) was 81%, while it was 77% for analysis of katG and inhA gene promoter to detect isoniazid-resistance via detection of high confidence mutations (Ser315Thr, C(-15)T, T(-8)C). Furthermore, considering the circulating TB genotypes in Uganda, the isoniazid-resistance conferring mutations were more frequent in M. tuberculosis lineage 4/sub-lineage Uganda, perhaps explaining why this genotype is weakly associated with MDR-TB.

CONCLUSION

Sequence analysis of rpoB/RRDR, katG and inhA gene promoter is useful in detecting rifampicin/isoniazid resistant M. tuberculosis isolates in Uganda however, about ≤20% of the resistant isolates lack known resistance-conferring mutations hence rapid molecular assays may not detect them as resistant.

Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis

To characterize the genetic determinants of resistance to antituberculosis drugs, we performed a genome-wide association study (GWAS) of 6,465 Mycobacterium tuberculosis clinical isolates from more than 30 countries. A GWAS approach within a mixed-regression framework was followed by a phylogenetics-based test for independent mutations. In addition to mutations in established and recently described resistance-associated genes, novel mutations were discovered for resistance to cycloserine, ethionamide and para-aminosalicylic acid. The capacity to detect mutations associated with resistance to ethionamide, pyrazinamide, capreomycin, cycloserine and para-aminosalicylic acid was enhanced by inclusion of insertions and deletions. Odds ratios for mutations within candidate genes were found to reflect levels of resistance. New epistatic relationships between candidate drug-resistance-associated genes were identified. Findings also suggest the involvement of efflux pumps (drrA and Rv2688c) in the emergence of resistance. This study will inform the design of new diagnostic tests and expedite the investigation of resistance and compensatory epistatic mechanisms.

Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples

Routine full characterization of Mycobacterium tuberculosis is culture based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by using the Illumina MiSeq sequencer (40 smear-positive respiratory samples obtained after routine clinical testing and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis DNA in direct samples.

Genomic analysis of globally diverse Mycobacterium tuberculosis strains provides insights into the emergence and spread of multidrug resistance

Multidrug-resistant tuberculosis (MDR-TB), caused by drug resistant strains of Mycobacterium tuberculosis, is an increasingly serious problem worldwide. In this study, we examined a dataset of 5,310 M. tuberculosis whole genome sequences from five continents. Despite great diversity with respect to geographic point of isolation, genetic background and drug resistance, patterns of drug resistance emergence were conserved globally. We have identified harbinger mutations that often precede MDR. In particular, the katG S315T mutation, conferring resistance to isoniazid, overwhelmingly arose before rifampicin resistance across all lineages, geographic regions, and time periods. Molecular diagnostics that include markers for rifampicin resistance alone will be insufficient to identify pre-MDR strains. Incorporating knowledge of pre-MDR polymorphisms, particularly katG S315, into molecular diagnostics will enable targeted treatment of patients with pre-MDR-TB to prevent further development of MDR-TB.

A brief primer on genomic epidemiology: lessons learned from Mycobacterium tuberculosis

Genomics is now firmly established as a technique for the investigation and reconstruction of communicable disease outbreaks, with many genomic epidemiology studies focusing on revealing transmission routes of Mycobacterium tuberculosis. In this primer, we introduce the basic techniques underlying transmission inference from genomic data, using illustrative examples from M. tuberculosis and other pathogens routinely sequenced by public health agencies. We describe the laboratory and epidemiological scenarios under which genomics may or may not be used, provide an introduction to sequencing technologies and bioinformatics approaches to identifying transmission-informative variation and resistance-associated mutations, and discuss how variation must be considered in the light of available clinical and epidemiological information to infer transmission.

The variability and reproducibility of whole genome sequencing technology for detecting resistance to anti-tuberculous drugs

BACKGROUND

The emergence of resistance to anti-tuberculosis drugs is a serious and growing threat to public health. Next-generation sequencing is rapidly gaining traction as a diagnostic tool for investigating drug resistance in Mycobacterium tuberculosis to aid treatment decisions. However, there are few little data regarding the precision of such sequencing for assigning resistance profiles.

METHODS

We investigated two sequencing platforms (Illumina MiSeq, Ion Torrent PGM™) and two rapid analytic pipelines (TBProfiler, Mykrobe predictor) using a well characterised reference strain (H37Rv) and clinical isolates from patients with tuberculosis resistant to up to 13 drugs. Results were compared to phenotypic drug susceptibility testing. To assess analytical robustness individual DNA samples were subjected to repeated sequencing.

RESULTS

The MiSeq and Ion PGM systems accurately predicted drug-resistance profiles and there was high reproducibility between biological and technical sample replicates. Estimated variant error rates were low (MiSeq 1 per 77 kbp, Ion PGM 1 per 41 kbp) and genomic coverage high (MiSeq 51-fold, Ion PGM 53-fold). MiSeq provided superior coverage in GC-rich regions, which translated into incremental detection of putative genotypic drug-specific resistance, including for resistance to para-aminosalicylic acid and pyrazinamide. The TBProfiler bioinformatics pipeline was concordant with reported phenotypic susceptibility for all drugs tested except pyrazinamide and para-aminosalicylic acid, with an overall concordance of 95.3%. When using the Mykrobe predictor concordance with phenotypic testing was 73.6%.

CONCLUSIONS

We have demonstrated high comparative reproducibility of two sequencing platforms, and high predictive ability of the TBProfiler mutation library and analytical pipeline, when profiling resistance to first- and second-line anti-tuberculosis drugs. However, platform-specific variability in coverage of some genome regions may have implications for predicting resistance to specific drugs. These findings may have implications for future clinical practice and thus deserve further scrutiny, set within larger studies and using updated mutation libraries.

The role of whole genome sequencing in antimicrobial susceptibility testing of bacteria: report from the EUCAST Subcommittee

Whole genome sequencing (WGS) offers the potential to predict antimicrobial susceptibility from a single assay. The European Committee on Antimicrobial Susceptibility Testing established a subcommittee to review the current development status of WGS for bacterial antimicrobial susceptibility testing (AST). The published evidence for using WGS as a tool to infer antimicrobial susceptibility accurately is currently either poor or non-existent and the evidence / knowledge base requires significant expansion. The primary comparators for assessing genotypic-phenotypic concordance from WGS data should be changed to epidemiological cut-off values in order to improve differentiation of wild-type from non-wild-type isolates (harbouring an acquired resistance). Clinical breakpoints should be a secondary comparator. This assessment will reveal whether genetic predictions could also be used to guide clinical decision making. Internationally agreed principles and quality control (QC) metrics will facilitate early harmonization of analytical approaches and interpretive criteria for WGS-based predictive AST. Only data sets that pass agreed QC metrics should be used in AST predictions. Minimum performance standards should exist and comparative accuracies across different WGS laboratories and processes should be measured. To facilitate comparisons, a single public database of all known resistance loci should be established, regularly updated and strictly curated using minimum standards for the inclusion of resistance loci. For most bacterial species the major limitations to widespread adoption for WGS-based AST in clinical laboratories remain the current high-cost and limited speed of inferring antimicrobial susceptibility from WGS data as well as the dependency on previous culture because analysis directly on specimens remains challenging. For most bacterial species there is currently insufficient evidence to support the use of WGS-inferred AST to guide clinical decision making. WGS-AST should be a funding priority if it is to become a rival to phenotypic AST. This report will be updated as the available evidence increases.

In silico region of difference (RD) analysis of Mycobacterium tuberculosis complex from sequence reads using RD-Analyzer

BACKGROUND

Whole-genome sequencing is increasingly used in clinical diagnosis of tuberculosis and study of Mycobacterium tuberculosis complex (MTC). MTC consists of several genetically homogenous mycobacteria species which can cause tuberculosis in humans and animals. Regions of difference (RDs) are commonly regarded as gold standard genetic markers for MTC classification.

RESULTS

We develop RD-Analyzer, a tool that can accurately infer the species and lineage of MTC isolates from sequence reads based on the presence and absence of a set of 31 RDs. Applied on a publicly available diverse set of 377 sequenced MTC isolates from known major species and lineages, RD-Analyzer achieved an accuracy of 98.14 % (370/377) in species prediction and a concordance of 98.47 % (257/261) in Mycobacterium tuberculosis lineage prediction compared to predictions based on single nucleotide polymorphism markers. By comparing respective sequencing read depths on each genomic position between isolates of different sublineages, we were able to identify the known RD markers in different sublineages of Lineage 4 and provide support for six potential delineating markers having high sensitivities and specificities for sublineage prediction. An extended version of RD-Analyzer was thus developed to allow user-defined RDs for lineage prediction.

CONCLUSIONS

RD-Analyzer is a useful and accurate tool for species, lineage and sublineage prediction using known RDs of MTC from sequence reads and is extendable to accepting user-defined RDs for analysis. RD-Analyzer is written in Python and is freely available at https://github.com/xiaeryu/RD-Analyzer .

Mycobacterium tuberculosis lineage 4 comprises globally distributed and geographically restricted sublineages

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration.

Whole genome sequencing of Mycobacterium tuberculosis for detection of drug resistance: a systematic review

OBJECTIVES

We conducted a systematic review to determine the diagnostic accuracy of whole genome sequencing (WGS) of Mycobacterium tuberculosis for the detection of resistance to first- and second-line anti-tuberculosis (TB) drugs.

METHODS

The study was conducted according to the criteria of the Preferred Reporting Items for Systematic Reviews group. A total of 20 publications were included. The sensitivity, specificity, positive-predictive value and negative-predictive value of WGS using phenotypic drug susceptibility testing methods as a reference standard were determined.

RESULTS

Anti-TB agents tested included all first-line drugs, a variety of reserve drugs, as well as new drugs. Polymorphisms in a total of 53 genes were tested for associations with drug resistance. Pooled sensitivity and specificity values for detection of resistance to selected first-line drugs were 0.98 (95% CI 0.93-0.98) and 0.98 (95% CI 0.98-1.00) for rifampicin and 0.97 (95% CI 0.94-0.99) and 0.93 (95% CI 0.91-0.96) for isoniazid, respectively. Due to high heterogeneity in study designs, lack of data, knowledge of resistance mechanisms and clarity on exclusion of phylogenetic markers, there was a significant variation in analytical performance of WGS for the remaining first-line, reserved drugs and new drugs.

CONCLUSIONS

Whole genome sequencing could be considered a promising alternative to existing phenotypic and molecular drug susceptibility testing methods for rifampicin and isoniazid pending standardization of analytical pipelines. To ensure clinical relevance of WGS for detection of M. tuberculosis complex drug resistance, future studies should include information on clinical outcomes.

Preventing the spread of multidrug-resistant tuberculosis and protecting contacts of infectious cases

Prevention of multidrug-resistant and extensively drug-resistant tuberculosis (MDR/XDR-TB) is a top priority for global TB control, given the need to limit epidemic spread and considering the high cost, toxicity and poor treatment outcomes with available therapies. We performed a systematic literature review to evaluate the evidence for strategies to reduce MDR/XDR-TB transmission and disease progression. Rapid detection and timely initiation of effective treatment is critical to rendering MDR/XDR-TB cases non-infectious. The scale-up of rapid molecular testing has transformed the capacity of high-incidence settings to identify and treat patients with MDR/XDR-TB. Optimized infection control measures in hospitals and clinics are critical to protect other patients and healthcare workers, whereas creative measures to reduce transmission within community hotspots require consideration. Targeted screening of high-risk communities may enhance early case-detection and limit the spread of MDR/XDR-TB. Among infected contacts, preventive therapy promises to reduce the risk of disease progression. This is supported by observational cohort studies, but randomized trials are urgently needed to confirm these observations and guide policy formulation. Substantial investment in MDR/XDR-TB prevention and care will be critical if the ambitious global goal of TB elimination is to be realized.

Evaluation of resistance acquisition during tuberculosis treatment using whole genome sequencing

Tuberculosis (TB) is still considered a major global public health problem in the world and there is a concern about the worldwide increase of drug-resistance (DR). This paper describes the analysis of three Mycobacterium tuberculosis isolates from a single patient collected over a long treatment period of time. DR was initially investigated through phenotypic testing, followed by line probe assays (LPAs) and whole genome sequencing (WGS). It presents an intriguing situation where a multidrug-resistant (MDR-) TB case was diagnosed and treated based only on late phenotypic drug susceptibility testing of isolate 1. During the treatment, another two isolates were cultivated: isolate 2, nine months after starting MDR-TB treatment; and isolate 3, cultivated five months later, during regular use of anti-TB drugs. These two isolates were evaluated using molecular LPA and WGS, retrospectively. All mutations detected by LPA were also detected in the WGS, including conversion from fluoroquinolones susceptibility to resistance from isolate 2 to isolate 3. WGS showed additional mutations, including some which may confer resistance to other drugs not tested (terizidone/cycloserine) and mutations with no correspondent resistance in drug susceptibility testing (streptomycin and second-line injectable drugs).

Mycobacterium tuberculosis whole genome sequencing and protein structure modelling provides insights into anti-tuberculosis drug resistance

BACKGROUND

Combating the spread of drug resistant tuberculosis is a global health priority. Whole genome association studies are being applied to identify genetic determinants of resistance to anti-tuberculosis drugs. Protein structure and interaction modelling are used to understand the functional effects of putative mutations and provide insight into the molecular mechanisms leading to resistance.

METHODS

To investigate the potential utility of these approaches, we analysed the genomes of 144 Mycobacterium tuberculosis clinical isolates from The Special Programme for Research and Training in Tropical Diseases (TDR) collection sourced from 20 countries in four continents. A genome-wide approach was applied to 127 isolates to identify polymorphisms associated with minimum inhibitory concentrations for first-line anti-tuberculosis drugs. In addition, the effect of identified candidate mutations on protein stability and interactions was assessed quantitatively with well-established computational methods.

RESULTS

The analysis revealed that mutations in the genes rpoB (rifampicin), katG (isoniazid), inhA-promoter (isoniazid), rpsL (streptomycin) and embB (ethambutol) were responsible for the majority of resistance observed. A subset of the mutations identified in rpoB and katG were predicted to affect protein stability. Further, a strong direct correlation was observed between the minimum inhibitory concentration values and the distance of the mutated residues in the three-dimensional structures of rpoB and katG to their respective drugs binding sites.

CONCLUSIONS

Using the TDR resource, we demonstrate the usefulness of whole genome association and convergent evolution approaches to detect known and potentially novel mutations associated with drug resistance. Further, protein structural modelling could provide a means of predicting the impact of polymorphisms on drug efficacy in the absence of phenotypic data. These approaches could ultimately lead to novel resistance mutations to improve the design of tuberculosis control measures, such as diagnostics, and inform patient management.

Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review

BACKGROUND

Whole genome sequencing (WGS) is becoming an important part of epidemiological investigations of infectious diseases due to greater resolution and cost reductions compared to traditional typing approaches. Many public health and clinical teams will increasingly use WGS to investigate clusters of potential pathogen transmission, making it crucial to understand the benefits and assumptions of the analytical methods for investigating the data. We aimed to understand how different approaches affect inferences of transmission dynamics and outline limitations of the methods.

METHODS

We comprehensively searched electronic databases for studies that presented methods used to interpret WGS data for investigating tuberculosis (TB) transmission. Two authors independently selected studies for inclusion and extracted data. Due to considerable methodological heterogeneity between studies, we present summary data with accompanying narrative synthesis rather than pooled analyses.

RESULTS

Twenty-five studies met our inclusion criteria. Despite the range of interpretation tools, the usefulness of WGS data in understanding TB transmission often depends on the amount of genetic diversity in the setting. Where diversity is small, distinguishing re-infections from relapses may be impossible; interpretation may be aided by the use of epidemiological data, examining minor variants and deep sequencing. Conversely, when within-host diversity is large, due to genetic hitchhiking or co-infection of two dissimilar strains, it is critical to understand how it arose. Greater understanding of microevolution and mixed infection will enhance interpretation of WGS data.

CONCLUSIONS

As sequencing studies have sampled more intensely and integrated multiple sources of information, the understanding of TB transmission and diversity has grown, but there is still much to be learnt about the origins of diversity that will affect inferences from these data. Public health teams and researchers should combine epidemiological, clinical and WGS data to strengthen investigations of transmission.

Rapid, comprehensive, and affordable mycobacterial diagnosis with whole-genome sequencing: a prospective study

Background

Slow and cumbersome laboratory diagnostics for Mycobacterium tuberculosis complex (MTBC) risk delayed treatment and poor patient outcomes. Whole-genome sequencing (WGS) could potentially provide a rapid and comprehensive diagnostic solution. In this prospective study, we compare real-time WGS with routine MTBC diagnostic workflows.

Methods

We compared sequencing mycobacteria from all newly positive liquid cultures with routine laboratory diagnostic workflows across eight laboratories in Europe and North America for diagnostic accuracy, processing times, and cost between Sept 6, 2013, and April 14, 2014. We sequenced specimens once using local Illumina MiSeq platforms and processed data centrally using a semi-automated bioinformatics pipeline. We identified species or complex using gene presence or absence, predicted drug susceptibilities from resistance-conferring mutations identified from reference-mapped MTBC genomes, and calculated genetic distance to previously sequenced UK MTBC isolates to detect outbreaks. WGS data processing and analysis was done by staff masked to routine reference laboratory and clinical results. We also did a microcosting analysis to assess the financial viability of WGS-based diagnostics.

Findings

Compared with routine results, WGS predicted species with 93% (95% CI 90–96; 322 of 345 specimens; 356 mycobacteria specimens submitted) accuracy and drug susceptibility also with 93% (91–95; 628 of 672 specimens; 168 MTBC specimens identified) accuracy, with one sequencing attempt. WGS linked 15 (16% [95% CI 10–26]) of 91 UK patients to an outbreak. WGS diagnosed a case of multidrug-resistant tuberculosis before routine diagnosis was completed and discovered a new multidrug-resistant tuberculosis cluster. Full WGS diagnostics could be generated in a median of 9 days (IQR 6–10), a median of 21 days (IQR 14–32) faster than final reference laboratory reports were produced (median of 31 days [IQR 21–44]), at a cost of £481 per culture-positive specimen, whereas routine diagnosis costs £518, equating to a WGS-based diagnosis cost that is 7% cheaper annually than are present diagnostic workflows.

Interpretation

We have shown that WGS has a scalable, rapid turnaround, and is a financially feasible method for full MTBC diagnostics. Continued improvements to mycobacterial processing, bioinformatics, and analysis will improve the accuracy, speed, and scope of WGS-based diagnosis.

Funding

National Institute for Health Research, Department of Health, Wellcome Trust, British Colombia Centre for Disease Control Foundation for Population and Public Health, Department of Clinical Microbiology, Trinity College Dublin.

Burden of transmitted multidrug resistance in epidemics of tuberculosis: a transmission modelling analysis

Background

Multidrug-resistant tuberculosis (MDR-TB) can be acquired through de novo mutation during TB treatment or through transmission from other individuals with active MDR-TB. Understanding the balance between these two mechanisms is essential when allocating resources for MDR-TB.

Methods

We constructed a dynamic transmission model of an MDR-TB epidemic, allowing for both treatment-related acquisition and person-to-person transmission of resistance. We used national TB notification data to inform Bayesian estimates of the fraction of each country’s 2013 MDR-TB incidence that resulted from MDR transmission rather than treatment-related MDR acquisition.

Findings

Global estimates of 3·5% MDR-TB prevalence among new TB notifications and 20·5% among retreatment notifications translate into an estimate that resistance transmission rather than acquisition accounts for a median 96% (95% UR: 68–100%) of all incident MDR-TB, and 61% (16–95%) of incident MDR-TB in previously-treated individuals. The estimated percentage of MDR-TB resulting from transmission varied substantially with different countries’ notification data; for example, we estimated this percentage at 48% (30–75%) of MDR-TB in Bangladesh, versus 99% (91–100%) in Uzbekistan. Estimates were most sensitive to estimates of the transmissibility of MDR strains, the probability of acquiring MDR during tuberculosis treatment, and the responsiveness of MDR TB to first-line treatment.

Interpretation

Notifications of MDR prevalence from most high-burden settings are most consistent with the vast majority of incident MDR-TB resulting from transmission rather than new treatment-related acquisition of resistance. Merely improving the treatment of drug-susceptible TB is unlikely to greatly reduce future MDR-TB incidence. Improved diagnosis and treatment of MDR-TB – including new tests and drug regimens – should be highly prioritized.

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread

The transcontinental spread of multidrug-resistant (MDR) tuberculosis is poorly characterized in molecular epidemiologic studies. We used genomic sequencing to understand the establishment and dispersion of MDR Mycobacterium tuberculosis within a group of immigrants to the United States. We used a genomic epidemiology approach to study a genotypically matched (by spoligotype, IS6110 restriction fragment length polymorphism, and mycobacterial interspersed repetitive units-variable number of tandem repeat signature) lineage 2/Beijing MDR strain implicated in an outbreak of tuberculosis among refugees in Thailand and consecutive cases within California. All 46 MDR M. tuberculosis genomes from both Thailand and California were highly related, with a median difference of 10 single-nucleotide polymorphisms (SNPs). The Wat Tham Krabok (WTK) strain is a new sequence type distinguished from all known Beijing strains by 55 SNPs and a genomic deletion (Rv1267c) associated with increased fitness. Sequence data revealed a highly prevalent MDR strain that included several closely related but distinct allelic variants within Thailand, rather than the occurrence of a single outbreak. In California, sequencing data supported multiple independent introductions of WTK with subsequent transmission and reactivation within the state, as well as a potential super spreader with a prolonged infectious period. Twenty-seven drug resistance-conferring mutations and 4 putative compensatory mutations were found within WTK strains. Genomic sequencing has substantial epidemiologic value in both low- and high-burden settings in understanding transmission chains of highly prevalent MDR strains.