Search

Mutations in Mycobacterium tuberculosis associated with resistance to antibiotics often come with a fitness cost for the bacteria. Resistance to the first-line drug rifampicin leads to lower competitive fitness of M. tuberculosis populations when compared to susceptible populations. This fitness cost, introduced by resistance mutations in the RNA polymerase, can be alleviated by compensatory mutations (CMs) in other regions of the affected protein. CMs are of particular interest clinically since they could lock in resistance mutations, encouraging the spread of resistant strains worldwide. Here, we report the statistical inference of a comprehensive set of CMs in the RNA polymerase of M. tuberculosis , using over 70 000 M . tuberculosis genomes that were collated as part of the CRyPTIC project. The unprecedented size of this data set gave the statistical tests more power to investigate the association of putative CMs with resistance-conferring mutations. Overall, we propose 51 high-confidence CMs by means of statistical association testing and suggest hypotheses for how they exert their compensatory mechanism by mapping them onto the protein structure. In addition, we were able to show an association of CMs with higher in vitro growth densities, and hence presumably with higher fitness, in resistant samples in the more virulent M. tuberculosis lineage 2. Our results suggest the association of CM presence with significantly higher in vitro growth than for wild-type samples, although this association is confounded with lineage and sub-lineage affiliation. Our findings emphasize the integral role of CMs and lineage affiliation in resistance spread and increases the urgency of antibiotic stewardship, which implies accurate, cheap and widely accessible diagnostics for M. tuberculosis infections to not only improve patient outcomes but also prevent the spread of resistant strains.

Background: Bedaquiline is a core drug for the treatment of multidrug-resistant tuberculosis; however, the understanding of resistance mechanisms is poor, which is hampering rapid molecular diagnostics. Some bedaquiline-resistant mutants are also cross-resistant to clofazimine. To decipher bedaquiline and clofazimine resistance determinants, we combined experimental evolution, protein modelling, genome sequencing, and phenotypic data., Methods: For this in-vitro and in-silico data analysis, we used a novel in-vitro evolutionary model using subinhibitory drug concentrations to select bedaquiline-resistant and clofazimine-resistant mutants. We determined bedaquiline and clofazimine minimum inhibitory concentrations and did Illumina and PacBio sequencing to characterise selected mutants and establish a mutation catalogue. This catalogue also includes phenotypic and genotypic data of a global collection of more than 14 000 clinical Mycobacterium tuberculosis complex isolates, and publicly available data. We investigated variants implicated in bedaquiline resistance by protein modelling and dynamic simulations., Findings: We discerned 265 genomic variants implicated in bedaquiline resistance, with 250 (94%) variants affecting the transcriptional repressor (Rv0678) of the MmpS5-MmpL5 efflux system. We identified 40 new variants in vitro, and a new bedaquiline resistance mechanism caused by a large-scale genomic rearrangement. Additionally, we identified in vitro 15 (7%) of 208 mutations found in clinical bedaquiline-resistant isolates. From our in-vitro work, we detected 14 (16%) of 88 mutations so far identified as being associated with clofazimine resistance and also seen in clinically resistant strains, and catalogued 35 new mutations. Structural modelling of Rv0678 showed four major mechanisms of bedaquiline resistance: impaired DNA binding, reduction in protein stability, disruption of protein dimerisation, and alteration in affinity for its fatty acid ligand., Interpretation: Our findings advance the understanding of drug resistance mechanisms in M tuberculosis complex strains. We have established an extended mutation catalogue, comprising variants implicated in resistance and susceptibility to bedaquiline and clofazimine. Our data emphasise that genotypic testing can delineate clinical isolates with borderline phenotypes, which is essential for the design of effective treatments., Funding: Leibniz ScienceCampus Evolutionary Medicine of the Lung, Deutsche Forschungsgemeinschaft, Research Training Group 2501 TransEvo, Rhodes Trust, Stanford University Medical Scientist Training Program, National Institute for Health and Care Research Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, Bill & Melinda Gates Foundation, Wellcome Trust, and Marie Skłodowska-Curie Actions., Competing Interests: Declaration of interests AS is the cofounder and owner of shares in BiKi Technology, which sells the BiKi Life Sciences software suite used in this study's analysis. PWF has received consulting fees from Global Pathogen Analysis Service. All other authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Drug resistant Mycobacterium tuberculosis, which mostly results from single nucleotide polymorphisms in antibiotic target genes, poses a major threat to tuberculosis treatment outcomes. Relative binding free energy (RBFE) calculations can rapidly predict the effects of mutations, but this approach has not been tested on large, complex proteins. We use RBFE calculations to predict the effects of M. tuberculosis RNA polymerase and DNA gyrase mutations on rifampicin and moxifloxacin susceptibility respectively. These mutations encompass a range of amino acid substitutions with known effects and include large steric perturbations and charged moieties. We find that moderate numbers (n = 3-15) of short RBFE calculations can predict resistance in cases where the mutation results in a large change in the binding free energy. We show that the method lacks discrimination in cases with either a small change in energy or that involve charged amino acids, and we investigate how these calculation errors may be decreased., (© 2022 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

Tuberculosis is a respiratory disease that is treatable with antibiotics. An increasing prevalence of resistance means that to ensure a good treatment outcome it is desirable to test the susceptibility of each infection to different antibiotics. Conventionally, this is done by culturing a clinical sample and then exposing aliquots to a panel of antibiotics, each being present at a pre-determined concentration, thereby determining if the sample isresistant or susceptible to each sample. The minimum inhibitory concentration (MIC) of a drug is the lowestconcentration that inhibits growth and is a more useful quantity but requires each sample to be tested at a range ofconcentrations for each drug. Using 96-well broth micro dilution plates with each well containing a lyophilised pre-determined amount of an antibiotic is a convenient and cost-effective way to measure the MICs of several drugs at once for a clinical sample. Although accurate, this is still an expensive and slow process that requires highly-skilled and experienced laboratory scientists. Here we show that, through the BashTheBug project hosted on the Zooniverse citizen science platform, a crowd of volunteers can reproducibly and accurately determine the MICs for 13 drugs and that simply taking the median or mode of 11-17 independent classifications is sufficient. There is therefore a potential role for crowds to support (but not supplant) the role of experts in antibiotic susceptibility testing., Competing Interests: PF, CW, HS, TZ, EB, SH, AG, AR, SK, TW, TP, GM, CL, DC, DC, AW No competing interests declared, (© 2022, Fowler et al.)

Background: Molecular diagnostics are considered the most promising route to achieving rapid, universal drug susceptibility testing for Mycobacterium tuberculosis complex (MTBC). We aimed to generate a WHO endorsed catalogue of mutations to serve as a global standard for interpreting molecular information for drug resistance prediction., Methods: A candidate gene approach was used to identify mutations as associated with resistance, or consistent with susceptibility, for 13 WHO endorsed anti-tuberculosis drugs. 38,215 MTBC isolates with paired whole-genome sequencing and phenotypic drug susceptibility testing data were amassed from 45 countries. For each mutation, a contingency table of binary phenotypes and presence or absence of the mutation computed positive predictive value, and Fisher's exact tests generated odds ratios and Benjamini-Hochberg corrected p-values. Mutations were graded as Associated with Resistance if present in at least 5 isolates, if the odds ratio was >1 with a statistically significant corrected p-value, and if the lower bound of the 95% confidence interval on the positive predictive value for phenotypic resistance was >25%. A series of expert rules were applied for final confidence grading of each mutation., Findings: 15,667 associations were computed for 13,211 unique mutations linked to one or more drugs. 1,149/15,667 (7·3%) mutations were classified as associated with phenotypic resistance and 107/15,667 (0·7%) were deemed consistent with susceptibility. For rifampicin, isoniazid, ethambutol, fluoroquinolones, and streptomycin, the mutations' pooled sensitivity was >80%. Specificity was over 95% for all drugs except ethionamide (91·4%), moxifloxacin (91·6%) and ethambutol (93·3%). Only two resistance mutations were classified for bedaquiline, delamanid, clofazimine, and linezolid as prevalence of phenotypic resistance was low for these drugs., Interpretation: This first WHO endorsed catalogue of molecular targets for MTBC drug susceptibility testing provides a global standard for resistance interpretation. Its existence should encourage the implementation of molecular diagnostics by National Tuberculosis Programmes., Funding: UNITAID, Wellcome, MRC, BMGF., Competing Interests: Conflicts of interest C.U.K. is a consultant Becton Dickinson, the Foundation for Innovative New Diagnostics and the TB Alliance. C.U.K. is collaborating with Janssen, PZA Innovation and Thermo Fisher Scientific. C.U.K. worked as a consultant for QuantuMDx, the Stop TB Partnership, the World Health Organization (WHO) Global TB Programme and the WHO Regional Office for Europe. C.U.K. gave a paid educational talk for Oxford Immunotec. Hain Lifescience covered C.U.K.’s and accommodation to present at a meeting. C.U.K. is an unpaid advisor to BioVersys and GenoScreen. E.R. is employed by Public Health England and holds an honorary contract with Imperial College London. I.F.L. is Director of the Scottish Mycobacteria Reference Laboratory. S.N. receives funding from German Center for Infection Research, Excellenz Cluster Precision Medicine in Chronic Inflammation, Leibniz Science Campus Evolutionary Medicine of the LUNG (EvoLUNG)tion EXC 2167. P.S. is a consultant at Genoscreen. T.R. is funded by NIH and DoD and receives salary support from the non-profit organization FIND. T.R. is a co-founder, board member and shareholder of Verus Diagnostics Inc, a company that was founded with the intent of developing diagnostic assays. Verus Diagnostics was not involved in any way with data collection, analysis or publication of the results. T.R. has not received any financial support from Verus Diagnostics. UCSD Conflict of Interest office has reviewed and approved T.R.’s role in Verus Diagnostics Inc. T.R. is a co-inventor of a provisional patent for a TB diagnostic assay (provisional patent #: 63/048.989). T.R. is a co-inventor on a patent associated with the processing of TB sequencing data (European Patent Application No. 14840432.0 & USSN 14/912,918). T.R. has agreed to “donate all present and future interest in and rights to royalties from this patent” to UCSD to ensure that he does not receive any financial benefits from this patent. S.S. is working and holding ESOPs at HaystackAnalytics Pvt. Ltd. (Product: Using whole genome sequencing for drug susceptibility testing for Mycobacterium tuberculosis). G.F.G. is listed as an inventor on patent applications for RBD-dimer-based CoV vaccines. The patents for RBD-dimers as protein subunit vaccines for SARS-CoV-2 have been licensed to Anhui Zhifei Longcom Biopharmaceutical Co. Ltd, China. No other authors declare a conflict of interest. I.C. is a consultant for the Foundation for Innovative New Diagnostics. DAC reports funding from GlaxoSmithKline and consultancy fees from Biobeats, Oxford University Innovation, Sensyne Health.C.C. reports funding from FIND to his institution (Pathology Queensland, Queensland Department of Health) for his laboratory to perform molecular analytic studies (limits of detection) for new molecular platforms manufactured by Cepheid and Bioneer.

Background: A comprehensive understanding of the pre-existing genetic variation in genes associated with antibiotic resistance in the Mycobacterium tuberculosis complex (MTBC) is needed to accurately interpret whole-genome sequencing data for genotypic drug susceptibility testing (DST)., Methods: We investigated mutations in 92 genes implicated in resistance to 21 anti-tuberculosis drugs using the genomes of 405 phylogenetically diverse MTBC strains. The role of phylogenetically informative mutations was assessed by routine phenotypic DST data for the first-line drugs isoniazid, rifampicin, ethambutol, and pyrazinamide from a separate collection of over 7000 clinical strains. Selected mutations/strains were further investigated by minimum inhibitory concentration (MIC) testing., Results: Out of 547 phylogenetically informative mutations identified, 138 were classified as not correlating with resistance to first-line drugs. MIC testing did not reveal a discernible impact of a Rv1979c deletion shared by M. africanum lineage 5 strains on resistance to clofazimine. Finally, we found molecular evidence that some MTBC subgroups may be hyper-susceptible to bedaquiline and clofazimine by different loss-of-function mutations affecting a drug efflux pump subunit (MmpL5)., Conclusions: Our findings underline that the genetic diversity in MTBC has to be studied more systematically to inform the design of clinical trials and to define sound epidemiologic cut-off values (ECOFFs) for new and repurposed anti-tuberculosis drugs. In that regard, our comprehensive variant catalogue provides a solid basis for the interpretation of mutations in genotypic as well as in phenotypic DST assays.

M. tuberculosis grows slowly and is challenging to work with experimentally compared with many other bacteria. Although microtitre plates have the potential to enable high-throughput phenotypic testing of M. tuberculosis, they can be difficult to read and interpret. Here we present a software package, the Automated Mycobacterial Growth Detection Algorithm (AMyGDA), that measures how much M. tuberculosis is growing in each well of a 96-well microtitre plate. The plate used here has serial dilutions of 14 anti-tuberculosis drugs, thereby permitting the MICs to be elucidated. The three participating laboratories each inoculated 38 96-well plates with 15 known M. tuberculosis strains (including the standard H37Rv reference strain) and, after 2 weeks' incubation, measured the MICs for all 14 drugs on each plate and took a photograph. By analysing the images, we demonstrate that AMyGDA is reproducible, and that the MICs measured are comparable to those measured by a laboratory scientist. The AMyGDA software will be used by the Comprehensive Resistance Prediction for Tuberculosis: an International Consortium (CRyPTIC) to measure the drug susceptibility profile of a large number (>30000) of samples of M. tuberculosis from patients over the next few years.

The UKMYC5 plate is a 96-well microtiter plate designed by the CRyPTIC Consortium (Comprehensive Resistance Prediction for Tuberculosis: an International Consortium) to enable the measurement of MICs of 14 different antituberculosis (anti-TB) compounds for >30,000 clinical Mycobacterium tuberculosis isolates. Unlike the MYCOTB plate, on which the UKMYC5 plate is based, the UKMYC5 plate includes two new (bedaquiline and delamanid) and two repurposed (clofazimine and linezolid) compounds. UKMYC5 plates were tested by seven laboratories on four continents by use of a panel of 19 external quality assessment (EQA) strains, including H37Rv. To assess the optimal combination of reading method and incubation time, MICs were measured from each plate by two readers, using three methods (mirrored box, microscope, and Vizion digital viewing system), after 7, 10, 14, and 21 days of incubation. In addition, all EQA strains were subjected to whole-genome sequencing and phenotypically characterized by the 7H10/7H11 agar proportion method (APM) and by use of MGIT960 mycobacterial growth indicator tubes. We concluded that the UKMYC5 plate is optimally read using the Vizion system after 14 days of incubation, achieving an interreader agreement of 97.9% and intra- and interlaboratory reproducibility rates of 95.6% and 93.1%, respectively. The mirrored box had a similar reproducibility. Strains classified as resistant by APM, MGIT960, or the presence of mutations known to confer resistance consistently showed elevated MICs compared to those for strains classified as susceptible. Finally, the UKMYC5 plate records intermediate MICs for one strain for which the APM measured MICs close to the applied critical concentration, providing early evidence that the UKMYC5 plate can quantitatively measure the magnitude of resistance to anti-TB compounds that is due to specific genetic variation., (Copyright © 2018 Rancoita et al.)