BackgroundThe continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.Methods and FindingsWe performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937–1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974–1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988–1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.ConclusionsIn the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic e...

A more complete understanding of the genetic basis of drug resistance in Mycobacterium tuberculosis is critical for prompt diagnosis and optimal treatment, particularly for toxic second-line drugs like D-cycloserine. Here, we used whole-genome sequences from 498 strains of M. tuberculosis to identify novel resistance-conferring genotypes. By combining association and correlated evolution tests with strategies for amplifying signal from rare variants, we found that loss-of-function mutations in ald (Rv2780), encoding L-alanine dehydrogenase, were associated with unexplained drug resistance. Convergent evolution of this loss-of-function was observed exclusively among multidrug-resistant strains. Drug susceptibility testing established that ald loss-of-function conferred resistance to D-cycloserine, and susceptibility to the drug was partially restored by complementation of ald. Clinical strains with mutations in ald and alr exhibited increased resistance to D-cycloserine when cultured in vitro. Incorporation of D-cycloserine resistance in novel molecular diagnostics could allow for targeted utilization of this toxic drug among patients with susceptible infections.

Tuberculosis chemotherapy is dependent on the use of the antibiotic pyrazinamide, which is being threatened by emerging drug resistance. Resistance is mediated through mutations in the bacterial gene pncA. Methods for testing pyrazinamide susceptibility are difficult and rarely performed, and this means that the full spectrum of pncA alleles that confer clinical resistance to pyrazinamide is unknown. Here, we performed in vitro saturating mutagenesis of pncA to generate a comprehensive library of PncA polymorphisms resultant from a single-nucleotide polymorphism. We then screened it for pyrazinamide resistance both in vitro and in an infected animal model. We identify over 300 resistance-conferring substitutions. Strikingly, these mutations map throughout the PncA structure and result in either loss of enzymatic activity and/or decrease in protein abundance. Our comprehensive mutational and screening approach should stand as a paradigm for determining resistance mutations and their mechanisms of action.

Genetics-based drug susceptibility testing has improved the diagnosis of drug-resistant tuberculosis but is limited by our lack of knowledge of all resistance mechanisms. Next-generation sequencing has assisted in identifying the principal genetic mechanisms of resistance for many drugs, but a significant proportion of phenotypic drug resistance is unexplained genetically. Few studies have formally compared the transcriptomes of susceptible and resistant Mycobacterium tuberculosis strains. We carried out comparative whole-genome transcriptomics of extensively drug-resistant (XDR) clinical isolates using RNA sequencing (RNA-seq) to find novel transcription-mediated mechanisms of resistance. We identified a promoter mutation (t to c) at position −11 (t−11c) relative to the start codon of ethA that reduces the expression of a monooxygenase (EthA) that activates ethionamide. (In this article, nucleotide changes are lowercase and amino acid substitutions are uppercase.) Using a flow cytometry-based reporter assay, we show that the reduced transcription of ethA is not due to transcriptional repression by ethR. Clinical strains harboring this mutation were resistant to ethionamide. Other ethA promoter mutations were identified in a global genomic survey of resistant M. tuberculosis strains. These results demonstrate a new mechanism of ethionamide resistance that can cause high-level resistance when it is combined with other ethionamide resistance-conferring mutations. Our study revealed many other genes which were highly up- or downregulated in XDR strains, including a toxin-antitoxin module (mazF5 mazE5) and tRNAs (leuX and thrU). This suggests that global transcriptional modifications could contribute to resistance or the maintenance of bacterial fitness have also occurred in XDR strains.

We developed a pyrazinamidase gene DNA-sequencing method to rapidly identify pyrazinamide resistance-causing mutations in GenoLyse-treated, smear-positive sputum specimens. The sensitivity and specificity were 90.9 and 100%, respectively, compared to those of MGIT drug susceptibility testing, after the exclusion of synonymous mutations and nonsynonymous mutations previously associated with susceptibility to pyrazinamide.