Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid.

OBJECTIVES: To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug. METHODS: The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid. RESULTS: We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes. CONCLUSIONS: In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST.

Deep amplicon sequencing for culture-free prediction of susceptibility or resistance to 13 anti-tuberculous drugs.

Conventional molecular tests for detecting Mycobacterium tuberculosis complex (MTBC) drug resistance on clinical samples cover a limited set of mutations. Whole-genome sequencing (WGS) typically requires culture.Here, we evaluated the Deeplex Myc-TB targeted deep-sequencing assay for prediction of resistance to 13 anti-tuberculous drugs/drug classes, directly applicable on sputum.With MTBC DNA tests, the limit of detection was 100-1000 genome copies for fixed resistance mutations. Deeplex Myc-TB captured in silico 97.1-99.3% of resistance phenotypes correctly predicted by WGS from 3651 MTBC genomes. On 429 isolates, the assay predicted 92.2% of 2369 first- and second-line phenotypes, with a sensitivity of 95.3% and a specificity of 97.4%. 56 out of 69 (81.2%) residual discrepancies with phenotypic results involved pyrazinamide, ethambutol and ethionamide, and low-level rifampicin or isoniazid resistance mutations, all notoriously prone to phenotypic testing variability. Only two out of 91 (2.2%) resistance phenotypes undetected by Deeplex Myc-TB had known resistance-associated mutations by WGS analysis outside Deeplex Myc-TB targets. Phenotype predictions from Deeplex Myc-TB analysis directly on 109 sputa from a Djibouti survey matched those of MTBSeq/PhyResSE/Mykrobe, fed with WGS data from subsequent cultures, with a sensitivity of 93.5/98.5/93.1% and a specificity of 98.5/97.2/95.3%, respectively. Most residual discordances involved gene deletions/indels and 3-12% heteroresistant calls undetected by WGS analysis or natural pyrazinamide resistance of globally rare "Mycobacterium canettii" strains then unreported by Deeplex Myc-TB. On 1494 arduous sputa from a Democratic Republic of the Congo survey, 14 902 out of 19 422 (76.7%) possible susceptible or resistance phenotypes could be predicted culture-free.Deeplex Myc-TB may enable fast, tailored tuberculosis treatment.

Revised Interpretation of the Hain Lifescience GenoType MTBC To Differentiate Mycobacterium canettii and Members of the Mycobacterium tuberculosis Complex.

Using 894 phylogenetically diverse genomes of the Mycobacterium tuberculosis complex (MTBC), we simulated in silico the ability of the Hain Lifescience GenoType MTBC assay to differentiate the causative agents of tuberculosis. Here, we propose a revised interpretation of this assay to reflect its strengths (e.g., it can distinguish some strains of Mycobacterium canettii and variants of Mycobacterium bovis that are not intrinsically resistant to pyrazinamide) and limitations (e.g., Mycobacterium orygis cannot be differentiated from Mycobacterium africanum).

Rapid Microarray-Based Detection of Rifampin, Isoniazid, and Fluoroquinolone Resistance in Mycobacterium tuberculosis by Use of a Single Cartridge.

The rapid and robust identification of mutations in Mycobacterium tuberculosis complex (MTBC) strains mediating multidrug-resistant (MDR) and extensively drug-resistant (XDR) phenotypes is crucial to combating the MDR tuberculosis (TB) epidemic. Currently available molecular anti-TB drug susceptibility tests either are restricted to a single target or drug (i.e., the Xpert MTB/RIF test) or present a risk of cross-contamination due to the design limitations of the open platform (i.e., line probe assays). With a good understanding of the technical and commercial boundaries, we designed a test cartridge based on an oligonucleotide array into which dried reagents are introduced and which has the ability to identify MTBC strains resistant to isoniazid, rifampin, and the fluoroquinolones. The melting curve assay interrogates 43 different mutations in the rifampin resistance-determining region (RRDR) of rpoB, rpoB codon 572, katG codon 315, the inhA promoter region, and the quinolone resistance-determining region (QRDR) of gyrA in a closed cartridge system within 90 min. Assay performance was evaluated with 265 clinical MTBC isolates, including MDR/XDR, non-MDR, and fully susceptible isolates, from a drug resistance survey performed in Swaziland in 2009 and 2010. In 99.5% of the cases, the results were consistent with data previously acquired utilizing Sanger sequencing. The assay, which uses a closed cartridge system in combination with a battery-powered Alere q analyzer and which has the potential to extend the current gene target panel, could serve as a rapid and robust point-of-care test in settings lacking a comprehensive molecular laboratory infrastructure to differentiate TB patients infected with MDR and non-MDR strains and to assist clinicians with their early treatment decisions.

Lab-on-Chip-Based Platform for Fast Molecular Diagnosis of Multidrug-Resistant Tuberculosis.

We evaluated the performance of the molecular lab-on-chip-based VerePLEX Biosystem for detection of multidrug-resistant tuberculosis (MDR-TB), obtaining a diagnostic accuracy of more than 97.8% compared to sequencing and MTBDRplus assay for Mycobacterium tuberculosis complex and rifampin and isoniazid resistance detection on clinical isolates and smear-positive specimens. The speed, user-friendly interface, and versatility make it suitable for routine laboratory use.

PhyResSE: a Web Tool Delineating Mycobacterium tuberculosis Antibiotic Resistance and Lineage from Whole-Genome Sequencing Data.

Antibiotic-resistant tuberculosis poses a global threat, causing the deaths of hundreds of thousands of people annually. While whole-genome sequencing (WGS), with its unprecedented level of detail, promises to play an increasingly important role in diagnosis, data analysis is a daunting challenge. Here, we present a simple-to-use web service (free for academic use at http://phyresse.org). Delineating both lineage and resistance, it provides state-of-the-art methodology to life scientists and physicians untrained in bioinformatics. It combines elaborate data processing and quality control, as befits human diagnostics, with a treasure trove of validated resistance data collected from well-characterized samples in-house and worldwide.

Phylogenetic polymorphisms in antibiotic resistance genes of the Mycobacterium tuberculosis complex.

OBJECTIVES: Sequence analysis of known antibiotic resistance genes of the Mycobacterium tuberculosis complex (MTBC) is increasingly being used to infer phenotypic resistance to a variety of antibiotics. However, a clear understanding of the genotype-phenotype relationship is required to interpret genotypic susceptibility results accurately. In this context, it is particularly important to distinguish phylogenetically informative neutral polymorphisms from true resistance-conferring mutations. METHODS: Using a collection of 71 strains that encompasses all major MTBC genotypes, we mapped the genetic diversity in 18 genes that are known to be involved or were previously implicated in antibiotic resistance to eight current as well as two novel antibiotics. This included bedaquiline, capreomycin, ethambutol, fluoroquinolones, isoniazid, PA-824, para-aminosalicylic acid, prothionamide, rifampicin and streptomycin. Moreover, we included data from one of our prior studies that focused on two of the three known pyrazinamide resistance genes. RESULTS: We found 58 phylogenetic polymorphisms that were markers for the genotypes M. tuberculosis Beijing, Haarlem, Latin American-Mediterranean (LAM), East African Indian (EAI), Delhi/Central Asian (CAS), Ghana, Turkey (Tur), Uganda I and II, Ural and X-type, as well as for Mycobacterium africanum genotypes West African I (WA I) and II (WA II), Mycobacterium bovis, Mycobacterium caprae, Mycobacterium pinnipedii, Mycobacterium microti and Mycobacterium canettii. CONCLUSIONS: This study represents one of the most extensive overviews of phylogenetically informative polymorphisms in known resistance genes to date, and will serve as a resource for the design and interpretation of genotypic susceptibility assays.

rplC T460C identified as a dominant mutation in linezolid-resistant Mycobacterium tuberculosis strains.

The ribosomal L3 protein was identified as a novel target in linezolid (LZD)-resistant Mycobacterium tuberculosis strains. Next-generation sequencing confirmed rplC T460C as the sole mutation in an LZD-resistant M. tuberculosis H37Rv strain selected in vitro. Sequencing analysis revealed the rplC T460C mutation in eight further LZD-resistant isolates (three in vitro-selected mutants and five patient isolates, including isolates from three different patients that developed LZD resistance during treatment) but in none of the susceptible control strains (n = 84).

Importance of the genetic diversity within the Mycobacterium tuberculosis complex for the development of novel antibiotics and diagnostic tests of drug resistance.

Despite being genetically monomorphic, the limited genetic diversity within the Mycobacterium tuberculosis complex (MTBC) has practical consequences for molecular methods for drug susceptibility testing and for the use of current antibiotics and those in clinical trials. It renders some representatives of MTBC intrinsically resistant against one or multiple antibiotics and affects the spectrum and consequences of resistance mutations selected for during treatment. Moreover, neutral or silent changes within genes responsible for drug resistance can cause false-positive results with hybridization-based assays, which have been recently introduced to replace slower phenotypic methods. We discuss the consequences of these findings and propose concrete steps to rigorously assess the genetic diversity of MTBC to support ongoing clinical trials.

Sequence analysis for detection of first-line drug resistance in Mycobacterium tuberculosis strains from a high-incidence setting.

BACKGROUND: Drug resistance displays a problem for the therapy of Mycobacterium tuberculosis infections. For molecular resistance testing, it is essential to have precise knowledge on genomic variations involved in resistance development. However, data from high-incidence settings are only sparely available. Therefore we performed a systematic approach and analyzed a total of 97 M. tuberculosis strains from previously treated patients in Sierra Leone for mutations in katG, rpoB, rrs, rpsL, gidB, embB, pncA and where applicable in inhA and ahpC. Of the strains investigated 50 were either mono- or poly-resistant to isoniazid, rifampin, streptomycin, ethambutol and pyrazinamide or MDR and 47 fully susceptible strains served as controls. RESULTS: The majority of isoniazid and rifampin resistant strains had mutations in katG315 (71.9%) and rpoB531 (50%). However, rpoB mutations in codons 511, 516 and 533 were also detected in five rifampin susceptible strains. MIC determinations revealed low-level rifampin resistance for those strains. Thus, the sensitivity and specificity of sequencing of katG for detection of drug resistance were 86.7% and 100% and for sequencing of rpoB 100% and 93.8%, respectively.Strikingly, none of the streptomycin resistant strains had mutations in rrs, but 47.5% harboured mutations in rpsL. Further changes were detected in gidB. Among ethambutol resistant strains 46.7% had mutations at embB306. Pyrazinamide resistant strains displayed a variety of mutations throughout pncA. The specificities of sequencing of rpsL, embB and pncA for resistance detection were high (96-100%), whereas sensitivities were lower (48.8%, 73.3%, 70%). CONCLUSIONS: Our study reveals a good correlation between data from molecular and phenotypic resistance testing in this high-incidence setting. However, the fact that particular mutations in rpoB are not linked to high-level resistance is challenging and demonstrates that careful interpretation of molecular resistance assays is mandatory. In addition, certain variations, especially in gidB, appear to be phylogenetically informative polymorphisms rather than markers for drug resistance.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis.

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a "perfect storm" that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

Impact of Fgd1 and ddn diversity in Mycobacterium tuberculosis complex on in vitro susceptibility to PA-824.

PA-824 is a promising drug candidate for the treatment of tuberculosis (TB). It is in phase II clinical trials as part of the first newly designed regimen containing multiple novel antituberculosis drugs (PA-824 in combination with moxifloxacin and pyrazinamide). However, given that the genes involved in resistance against PA-824 are not fully conserved in the Mycobacterium tuberculosis complex (MTBC), this regimen might not be equally effective against different MTBC genotypes. To investigate this question, we sequenced two PA-824 resistance genes (fgd1 [Rv0407] and ddn [Rv3547]) in 65 MTBC strains representing major phylogenetic lineages. The MICs of representative strains were determined using the modified proportion method in the Bactec MGIT 960 system. Our analysis revealed single-nucleotide polymorphisms in both genes that were specific either for several genotypes or for individual strains, yet none of these mutations significantly affected the PA-824 MICs (≤ 0.25 µg/ml). These results were supported by in silico modeling of the mutations identified in Fgd1. In contrast, "Mycobacterium canettii" strains displayed a higher MIC of 8 µg/ml. In conclusion, we found a large genetic diversity in PA-824 resistance genes that did not lead to elevated PA-824 MICs. In contrast, M. canettii strains had MICs that were above the plasma concentrations of PA-824 documented so far in clinical trials. As M. canettii is also intrinsically resistant against pyrazinamide, new regimens containing PA-824 and pyrazinamide might not be effective in treating M. canettii infections. This finding has implications for the design of multiple ongoing clinical trials.

Thr202Ala in thyA is a marker for the Latin American Mediterranean lineage of the Mycobacterium tuberculosis complex rather than para-aminosalicylic acid resistance.

Single nucleotide polymorphisms (SNPs) involved in the development of resistance represent powerful markers for the rapid detection of first- and second-line resistance in clinical Mycobacterium tuberculosis complex (MTBC) isolates. However, the association between particular mutations and phenotypic resistance is not always clear-cut, and phylogenetic SNPs have been misclassified as resistance markers in the past. In the present study, we investigated the utility of a specific polymorphism in thyA (Thr202Ala) as a marker for resistance to para-aminosalicyclic acid (PAS). Sixty-three PAS-susceptible MTBC strains comprising all major phylogenetic lineages, reference strain H37Rv, and 135 multidrug-resistant (MDR) strains from Germany (comprising 8 PAS-resistant isolates) were investigated for the presence of Thr202Ala. In both strain collections, the Thr202Ala SNP was found exclusively in strains of the Latin American Mediterranean (LAM) lineage irrespective of PAS resistance. Furthermore, PAS MICs (0.5 mg/liter) for selected LAM strains (all containing the SNP) and non-LAM strains (not containing the SNP), as well as the results of growth curve analyses performed in liquid 7H9 medium in the presence of increasing PAS concentrations (0 to 2.0 mg/liter), were identical. In conclusion, our data demonstrate that the Thr202Ala polymorphism in thyA is not a valid marker for PAS resistance but, instead, represents a phylogenetic marker for the LAM lineage of the M. tuberculosis complex. These findings challenge some of the previous understanding of PAS resistance and, as a consequence, warrant further in-depth investigations of the genetic variation in PAS-resistant clinical isolates and spontaneous mutants.

embCAB sequence variation among ethambutol-resistant Mycobacterium tuberculosis isolates without embB306 mutation.

OBJECTIVES: Mechanisms of resistance to ethambutol in Mycobacterium tuberculosis remain inadequately described. Although there is mounting evidence that mutations of codon 306 in embB play a key role, a significant number of phenotypically ethambutol-resistant strains do not carry mutations in this codon. Here, other mutations in the embCAB operon are suggested to be involved in resistance development. METHODS: The entire embCAB operon ( approximately 10 kb) was analysed in 34 phenotypically ethambutol-resistant M. tuberculosis strains without mutations in embB306 and in 12 ethambutol-susceptible strains. Furthermore, 106 control strains were investigated for the presence of particular mutations only. RESULTS: Overall, 18 non-synonymous mutations in 15 distinct codons of the embCAB operon were identified in ethambutol-resistant strains but not in ethambutol-susceptible isolates. The majority occurred in the embB gene (10 distinct codons), in a 570 bp region also encompassing embB306. Mutations in embC and embA were found rarely and in most cases in combination with polymorphisms in embB. One synonymous mutation (embA 228 bp) and two non-synonymous mutations (embCVal981Leu and embCArg738Gln) were found in ethambutol-susceptible strains as well as resistant strains and were confirmed to represent phylogenetic markers for strains of the Beijing, Haarlem and Delhi/CAS genotypes, respectively. CONCLUSIONS: Besides mutations in embB306, mutations in embB406 and embB497 were confirmed as hot spots for genomic variation in ethambutol-resistant clinical isolates. Of all resistant strains 70.6% carry a mutation in a relatively short region in embB, which therefore represents a promising target for inclusion in molecular assays for rapid detection of ethambutol resistance.

Detection of low-frequency resistance-mediating SNPs in next-generation sequencing data of Mycobacterium tuberculosis complex strains with binoSNP.

Accurate drug resistance detection is key for guiding effective tuberculosis treatment. While genotypic resistance can be rapidly detected by molecular methods, their application is challenged by mixed mycobacterial populations comprising both susceptible and resistant cells (heteroresistance). For this, next-generation sequencing (NGS) based approaches promise the determination of variants even at low frequencies. However, accurate methods for a valid detection of low-frequency variants in NGS data are currently lacking. To tackle this problem, we developed the variant detection tool binoSNP which allows the determination of low-frequency single nucleotide polymorphisms (SNPs) in NGS datasets from Mycobacterium tuberculosis complex (MTBC) strains. By taking a reference-mapped file as input, binoSNP evaluates each genomic position of interest using a binomial test procedure. binoSNP was validated using in-silico, in-vitro, and serial patient isolates datasets comprising varying genomic coverage depths (100-500×) and SNP allele frequencies (1-30%). Overall, the detection limit for low-frequency SNPs depends on the combination of coverage depth and allele frequency of the resistance-associated mutation. binoSNP allows for valid detection of resistance associated SNPs at a 1% frequency with a coverage ≥400×. In conclusion, binoSNP provides a valid approach to detect low-frequency resistance-mediating SNPs in NGS data from clinical MTBC strains. It can be implemented in automated, end-user friendly analysis tools for NGS data and is a step forward towards individualized TB therapy.

Sequence analyses of just four genes to detect extensively drug-resistant Mycobacterium tuberculosis strains in multidrug-resistant tuberculosis patients undergoing treatment.

The rapid detection of Mycobacterium tuberculosis isolates resistant to second-line drugs is crucial for the institution of appropriate treatment regimens as early as possible. Although molecular methods have successfully been used for the rapid detection of resistance to first-line drugs, there are limited data on mutations that confer resistance to second-line drugs. To address this question, we analyzed Mycobacterium tuberculosis strains resistant to ofloxacin (n = 26) and to capreomycin and/or amikacin (n = 48) from Uzbekistan for variations in target genes (gyrA, gyrB, rrs, and tlyA). Strains susceptible to ofloxacin (n = 49) and capreomycin and/or amikacin (n = 39) were included as controls. Mutations in gyrA or gyrB were found in 96% (25/26 strains) of the ofloxacin-resistant strains, while none of the susceptible strains displayed mutations in those two genes. The most common mutation occurred in gyrA at codon 94 (17/26 strains [65.4%]), followed by mutations at codons 90 and 91. Two strains showed a mutation in gyrB, at codons 485 and 543, respectively; both mutations have not been reported previously. The most frequent mutation in strains resistant to both amikacin and capreomycin was A1401G in rrs (34/40 strains [85.0%]). Three strains had mutations in tlyA, of which two (at codons 18 and 118) were associated with resistance to capreomycin alone. Overall, none of the 10 resistant strains (5 amikacin-resistant and capreomycin-susceptible strains) and none of the 39 susceptible control strains had mutations in the genes investigated. Our results clearly demonstrate the potential of sequence analyses of short regions of relatively few target genes for the rapid detection of resistance to second-line drugs among strains isolated from patients undergoing treatment for multidrug-resistant tuberculosis. The mechanisms that confer amikacin resistance in this setting remain unclear.

Rapid microarray-based assay for detection of pyrazinamide resistant Mycobacterium tuberculosis.

Pyrazinamide (PZA) is a key antibiotic for the treatment of drug susceptible tuberculosis. PZA-resistance is mainly mediated by mutations in the pncA gene; however the current gold standard is a phenotypic drug susceptibility test requiring a well-adjusted pH-value for reliable results. Our melting curve assay detects a non-wild type genotype in selected pncA regions in at least 3750 gene copies/mL within 2.5 hours. The prototype assay was further evaluated by analyzing 271 Mycobacterium tuberculosis complex isolates from Swaziland originating from a previously published drug resistance survey and including 118 isolates with pncA mutations. Sensitivity was 83% (95% CI 75-89%) and specificity was 100% (95% CI 98-100%). Under consideration of further improvements with regard to the target range our melting curve assay has the potential as a rapid rule-in test for PZA susceptibility (wild type pncA), however false resistant results (mutant pncA, but PZA susceptible) cannot be ruled out completely.