Whole-genome sequencing for prediction of Mycobacterium tuberculosis drug susceptibility and resistance: a retrospective cohort study

Summary Background Diagnosing drug-resistance remains an obstacle to the elimination of tuberculosis. Phenotypic drug-susceptibility testing is slow and expensive, and commercial genotypic assays screen only common resistance-determining mutations. We used whole-genome sequencing to characterise common and rare mutations predicting drug resistance, or consistency with susceptibility, for all first-line and second-line drugs for tuberculosis. Methods Between Sept 1, 2010, and Dec 1, 2013, we sequenced a training set of 2099 Mycobacterium tuberculosis genomes. For 23 candidate genes identified from the drug-resistance scientific literature, we algorithmically characterised genetic mutations as not conferring resistance (benign), resistance determinants, or uncharacterised. We then assessed the ability of these characterisations to predict phenotypic drug-susceptibility testing for an independent validation set of 1552 genomes. We sought mutations under similar selection pressure to those characterised as resistance determinants outside candidate genes to account for residual phenotypic resistance. Findings We characterised 120 training-set mutations as resistance determining, and 772 as benign. With these mutations, we could predict 89·2% of the validation-set phenotypes with a mean 92·3% sensitivity (95% CI 90·7–93·7) and 98·4% specificity (98·1–98·7). 10·8% of validation-set phenotypes could not be predicted because uncharacterised mutations were present. With an in-silico comparison, characterised resistance determinants had higher sensitivity than the mutations from three line-probe assays (85·1% vs 81·6%). No additional resistance determinants were identified among mutations under selection pressure in non-candidate genes. Interpretation A broad catalogue of genetic mutations enable data from whole-genome sequencing to be used clinically to predict drug resistance, drug susceptibility, or to identify drug phenotypes that cannot yet be genetically predicted. This approach could be integrated into routine diagnostic workflows, phasing out phenotypic drug-susceptibility testing while reporting drug resistance early. Funding Wellcome Trust, National Institute of Health Research, Medical Research Council, and the European Union.

[1]  Samuel A. Assefa,et al.  Elucidating Emergence and Transmission of Multidrug-Resistant Tuberculosis in Treatment Experienced Patients by Whole Genome Sequencing , 2013, PloS one.

[2]  Tim E A Peto,et al.  Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007–12, with whole pathogen genome sequences: an observational study , 2014, The Lancet. Respiratory medicine.

[3]  T. Dallman,et al.  Performance comparison of benchtop high-throughput sequencing platforms , 2012, Nature Biotechnology.

[4]  K. Kam,et al.  Ethambutol Resistance as Determined by Broth Dilution Method Correlates Better than Sequencing Results with embB Mutations in Multidrug-Resistant Mycobacterium tuberculosis Isolates , 2013, Journal of Clinical Microbiology.

[5]  E. Böttger The ins and outs of Mycobacterium tuberculosis drug susceptibility testing. , 2011, Clinical microbiology and infection : the official publication of the European Society of Clinical Microbiology and Infectious Diseases.

[6]  Jukka Corander,et al.  Evolution and transmission of drug resistant tuberculosis in a Russian population , 2014, Nature Genetics.

[7]  K. Holt,et al.  Out-of-Africa migration and Neolithic co-expansion of Mycobacterium tuberculosis with modern humans , 2013, Nature Genetics.

[8]  Y. Balabanova,et al.  Rapid diagnostics of tuberculosis and drug resistance in the industrialized world: clinical and public health benefits and barriers to implementation , 2013, BMC Medicine.

[9]  A. Crook,et al.  Four-month moxifloxacin-based regimens for drug-sensitive tuberculosis. , 2014, The New England journal of medicine.

[10]  Jianming Wang,et al.  Rapid Diagnosis of Drug Resistance to Fluoroquinolones, Amikacin, Capreomycin, Kanamycin and Ethambutol Using Genotype MTBDRsl Assay: A Meta-Analysis , 2013, PloS one.

[11]  M. Pallen,et al.  Culture-independent detection and characterisation of Mycobacterium tuberculosis and M. africanum in sputum samples using shotgun metagenomics on a benchtop sequencer , 2014, PeerJ.

[12]  S. Borrell,et al.  KvarQ: targeted and direct variant calling from fastq reads of bacterial genomes , 2014, BMC Genomics.

[13]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[14]  G. Fischer,et al.  Next-Generation Ion Torrent Sequencing of Drug Resistance Mutations in Mycobacterium tuberculosis Strains , 2012, Journal of Clinical Microbiology.

[15]  B. D. de Jong,et al.  A four-month gatifloxacin-containing regimen for treating tuberculosis. , 2015, The New England journal of medicine.

[16]  M. Pallen,et al.  Whole-Genome Sequencing for Rapid Susceptibility Testing of M. tuberculosis , 2013 .

[17]  Linda C. van der Gaag,et al.  Probabilistic Graphical Models , 2014, Lecture Notes in Computer Science.

[18]  B. Barrell,et al.  Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence , 1998, Nature.

[19]  Phelim Bradley,et al.  Rapid antibiotic resistance predictions from genome sequence data for S. aureus and M. tuberculosis , 2015, bioRxiv.

[20]  Martin Goodson,et al.  Stampy: a statistical algorithm for sensitive and fast mapping of Illumina sequence reads. , 2011, Genome research.

[21]  V. Jarlier,et al.  Performance of MTBDR plus for detecting high/low levels of Mycobacterium tuberculosis resistance to isoniazid. , 2009, The international journal of tuberculosis and lung disease : the official journal of the International Union against Tuberculosis and Lung Disease.

[22]  George M Church,et al.  Tuberculosis Drug Resistance Mutation Database , 2009, PLoS medicine.

[23]  J. Musser,et al.  Mutations associated with pyrazinamide resistance in pncA of Mycobacterium tuberculosis complex organisms , 1997, Antimicrobial agents and chemotherapy.

[24]  V. Jarlier,et al.  Detection by GenoType MTBDRsl Test of Complex Mechanisms of Resistance to Second-Line Drugs and Ethambutol in Multidrug-Resistant Mycobacterium tuberculosis Complex Isolates , 2010, Journal of Clinical Microbiology.

[25]  F. Dafae,et al.  Sequence analysis for detection of first-line drug resistance in Mycobacterium tuberculosis strains from a high-incidence setting , 2012, BMC Microbiology.

[26]  Daniel J. Wilson,et al.  Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study , 2013, The Lancet. Infectious diseases.

[27]  Michael Eisenstein,et al.  Oxford Nanopore announcement sets sequencing sector abuzz , 2012, Nature Biotechnology.

[28]  L. Gabbasova,et al.  Global tuberculosis report (2014) , 2014 .

[29]  W. Sougakoff,et al.  Crystal Structure of the Pyrazinamidase of Mycobacterium tuberculosis: Insights into Natural and Acquired Resistance to Pyrazinamide , 2011, PloS one.

[30]  Phelim Bradley,et al.  Rapid antibiotic-resistance predictions from genome sequence data for Staphylococcus aureus and Mycobacterium tuberculosis , 2015, Nature Communications.

[31]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[32]  G. Kahlmeter,et al.  Challenging a dogma: antimicrobial susceptibility testing breakpoints for Mycobacterium tuberculosis. , 2012, Bulletin of the World Health Organization.

[33]  S. Niemann,et al.  Mycobacterium tuberculosis Pyrazinamide Resistance Determinants: a Multicenter Study , 2014, mBio.

[34]  W. Hanage,et al.  Comprehensive Identification of Single Nucleotide Polymorphisms Associated with Beta-lactam Resistance within Pneumococcal Mosaic Genes , 2014, PLoS genetics.

[35]  G. McVean,et al.  De novo assembly and genotyping of variants using colored de Bruijn graphs , 2011, Nature Genetics.

[36]  Razvan Sultana,et al.  Genomic Analysis Identifies Targets of Convergent Positive Selection in Drug Resistant Mycobacterium tuberculosis , 2013, Nature Genetics.

[37]  K. Floyd,et al.  Rapid molecular TB diagnosis: evidence, policy-making and global implementation of Xpert®MTB/RIF Weyer, , 2013 .

[38]  Daniel J. Wilson,et al.  ClonalFrameML: Efficient Inference of Recombination in Whole Bacterial Genomes , 2015, PLoS Comput. Biol..