An improved catalogue for whole-genome sequencing prediction of bedaquiline resistance in M. tuberculosis using a reproduciblealgorithmic approach.

Abstract

Bedaquiline (BDQ) has only been approved for use for a little over a decade yet is a key drug for treating multi-drug resistant tuberculosis, however rising levels of resistance threaten to reduce its effectiveness. Catalogues of mutations associated with resistance to bedaquiline are key to detecting resistance genetically for either diagnosis or surveillance. At present building catalogues requires considerable expert knowledge, often requires the use of complex grading rules, and is an irreproducible process. We developed an automated method, catomatic, that associates genetic variants with resistance (or susceptibility) using a two-tailed binomial test with a stated background rate and applied it to a dataset of 11,867 Mycobacterium tuberculosis samples with whole genome and bedaquline susceptibility testing data. Using this framework we investigated how to best classify variants and the phenotypic significance of minor alleles. The genes mmpS5 and mmpL5 are not directly associated with bedaquline resistance, and our catalogue of Rv0678, atpE, and pepQ variants attains a cross-validated sensitivity and specificity of 79.4 {+/-} 1.8 % and 98.5 {+/-} 0.3%, respectively, for 94 {+/-} 0.4% of samples. Identifying samples with subpopulations containing Rv0678 variants improves sensitivity, and detection thresholds in bioinformatic pipelines should therefore be lowered. By using a more permissive and deterministic algorithm trained on a sufficient number of resistant samples we have reproducibly constructed an AMR catalogue for BDQ resistance-associated variants that is comprehensive and accurate. Impact StatementBedaquiline has recently received global endorsement for tuberculosis treatment, yet the genetic determinants of antimicrobial resistance remain incompletely understood. Existing gold-standard methods for building mutation catalogs lack public accessibility and reproducibility. We introduce catomatic, a reproducible and publicly available method that employs simpler statistics to increase sensitivity to resistance-associated variants. This approach has enabled investigations into mechanisms of resistance, the significance of genetic subpopulations, and key data attributes that influence the ease of classifying effects and the accuracy of BDQ resistance phenotype prediction in clinical samples. We strongly emphasise the utility in using reproducible statistics and sustainably developed software in genetics-focussed microbiology.