A standardised method for interpreting the association between mutations and phenotypic drug resistance in Mycobacterium tuberculosis

Diagnostic performance analysis of Xpert MTB/RIF in lymph node tuberculosis: A retrospective study

AIMS

To retrospectively analyze the diagnostic efficacy of Xpert MTB/RIF (Xpert) in lymph node tuberculosis (LNTB).

METHODS

Sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV) and area under the curve (AUC) of Xpert, pathological examination and culture for LNTB were calculated.

RESULTS

421 suspected LNTB cases were categorized into the LNTB group (377 cases) and non-LNTB group (44 cases). The sensitivities of Xpert, pathological examination, and culture were 72.15%, 20.69%, 30.24%, respectively, with NPVs of 29.53%, 12.83%, 14.33%. The AUC values were 0.861, 0.603, 0.651, respectively. The sensitivity of Xpert varied across sample types: tissue (64.73%), puncture fluid (74.42%), and pus (96.05%). For specific lymph node locations, the sensitivity was head-and-neck (72.51%), mediastinal (84.21%), and axillary (45.83%).

CONCLUSIONS

Xpert demonstrates high diagnostic value for LNTB, particularly in pus samples. It also performs better in mediastinal and head-and-neck lymph node samples compared to axillary lymph node samples.

Sequencing by binding rivals SMOR error-corrected sequencing by synthesis technology for accurate detection and quantification of minor (< 0.1%) subpopulation variants

BACKGROUND

Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications, including the detection of drug resistant pathogens and somatic variant detection in oncology. A recently available sequencing approach termed 'sequencing by binding (SBB)' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using SBB for the detection of ultra-rare drug resistant subpopulations in Mycobacterium tuberculosis (Mtb) using a targeted amplicon assay and compares the performance of SBB to single molecule overlapping reads (SMOR) error corrected sequencing by synthesis (SBS) data.

RESULTS

SBS displayed an elevated error rate when compared to SMOR error-corrected SBS and SBB techniques. SMOR error-corrected SBS and SBB technologies performed similarly within the linear range studies and error rate studies.

CONCLUSIONS

With lower sequencing error rates within SBB sequencing, this technique looks promising for both targeted and unbiased whole genome sequencing, leading to the identification of minor (< 1%) subpopulations without the need for error correction methods.

Accuracy of cobas MTB and MTB-RIF/INH for Detection of Mycobacterium tuberculosis and Drug Resistance

This study evaluated the performance of cobas MTB and cobas MTB-RIF/INH for the diagnosis of tuberculosis and detection of rifampicin (RIF) and isoniazid (INH) resistance. Adults presenting with pulmonary tuberculosis symptoms were recruited in South Africa, Moldova, and India. Performance of cobas MTB was assessed against culture, whereas cobas MTB-RIF/INH was assessed using phenotypic drug susceptibility testing and whole-genome sequencing as composite reference standards. Xpert MTB/RIF (Xpert) or Xpert MTB/RIF Ultra (Ultra) was used as a comparator. The overall sensitivity and specificity of cobas MTB were 95% (95% CI, 93%-96%) and 96% (95% CI, 95%-97%). Among smear-negatives, the sensitivity of cobas MTB was 75% (95% CI, 66%-83%). Among participants tested with both cobas MTB and Xpert, sensitivity was 96% (95% CI, 94%-97%) for cobas MTB and 95% (95% CI, 93%-97%) for Xpert. Among participants tested with both cobas MTB and Ultra, sensitivity was 88% (95% CI, 81%-92%) for cobas MTB and 89% (95% CI, 83%-93%) for Ultra. Sensitivity and specificity of cobas MTB-RIF/INH for RIF and INH detection were 90% (95% CI, 84%-94%) and 100% (95% CI, 99%-100%), and 89% (95% CI, 84%-93%) and 99.5% (95% CI, 98%-100%), respectively. The cobas MTB and cobas MTB-RIF/INH assays exhibited high performance in a diverse population and present a suitable option for molecular detection of tuberculosis and RIF and INH resistance.

Surveillance of fluoroquinolones resistance in rifampicin-susceptible tuberculosis in eastern China with whole-genome sequencing-based approach

Background

Leveraging well-established DNA-level drug resistance mechanisms, whole-genome sequencing (WGS) has emerged as a valuable methodology for predicting drug resistance. As the most effective second-line anti-tuberculosis (anti-TB) drugs, fluoroquinoloness (FQs) are generally used to treat multidrug-resistant tuberculosis (MDR-TB, defined as being resistant to resistant to rifampicin and isoniazid) or rifampicin-resistant tuberculosis (RR-TB). However, FQs are also commonly used in the management of other bacterial infections. There are few published data on the rates of FQs resistance among rifampicin-susceptible TB. The prevalence of FQs resistance among TB patients who are rifampicin-susceptible has not been studied in Zhejiang Province, China. The goal of this study was to provide a baseline characterization of the prevalence of FQs resistance, particularly among rifampicin-susceptible TB in Zhejiang Province, China.

Methods

Based on WGS, we have investigated the prevalence of FQs resistance among rifampicin-susceptible TB in Zhejiang Province. All pulmonary TB patients with positive cultures who were identified in Zhejiang area during TB drug resistance surveillance from 2018 to 2019 have enrolled in this population-based retrospective study.

Results

The rate of FQs resistance was 4.6% (32/698) among TB, 4.0% (27/676) among rifampicin-susceptible TB, and 22.7% (5/22) among RR-TB. According to WGS, strains that differ within 12 single-nucleotide polymorphisms (SNPs) were considered to be transmission of FQ-resistant strains. Specifically, 3.7% (1/27) of FQs resistance was caused by the transmission of FQs-resistant strains among the rifampicin-susceptible TB and 40.7% (11/27) of FQs resistance was identified as hetero-resistance.

Conclusion

The prevalence of FQs resistance among TB patients who were rifampicin-susceptible was severe in Zhejiang. The emergence of FQs resistance in TB isolates that are rifampicin-susceptible was mainly caused by the selection of drug-resistant strains. In order to prevent the emergence of FQs resistance, the WGS-based surveillance system for TB should be urgently established, and clinical awareness of the responsible use of FQs for respiratory infections should be enhanced.

Evaluation of MeltPro Assay in Identification of Second-Line Injectable Drug Resistance in Multidrug-Resistant Tuberculosis Isolates

Objective

We compared the MeltPro assay to whole-genome sequencing (WGS) to investigate the molecular characterization of second-line injectable drug (SLID) resistance in multidrug-resistant tuberculosis (MDR-TB) isolates in Chongqing, China.

Methods

A total of 122 MDR-TB patient isolates were collected between March 2019 and June 2020 from Chongqing Municipality, China. Conventional drug-susceptibility testing was performed using the proportion method, followed to generate minimum inhibitory concentrations (MICs) of SLIDs determined by microplate alamarblue assay. All strains were subjected to both MeltPro and WGS assays.

Results

Among 122 MDR-TB isolates, 30 (24.6%), 22 (18.0%), and 14 (11.5%) were resistant to kanamycin (KM), amikacin (AM), and capreomycin (CM), respectively. Of the 31 SLID-resistant isolates, 24 (77.4%, 24/31) isolates harbored mutations in the rrs gene, with the most prevalent mutations in rrs A1401G (22/24, 91.7%). Mutation in rrs A1401G was associated with high levels of resistance to KM (MIC, ≥40 μg/mL) and AM (MIC, ≥64 μg/mL), but disparities in CM-resistance levels. Using phenotypic drug-susceptibility testing as gold standard, we found that the overall sensitivity of MeltPro and WGS was 87.1% and 90.32% and specificity 100% and 97.8%, respectively. Seven isolates had discordant results between phenotypic and genotypic resistance of SLIDs.

Conclusion

MeltPro is a promising diagnostic tool for accurate identification of SLID-resistant MTB isolates with mutations in the rrs and eis genes. There was a disparity between MeltPro with WGS results in the proportion of heterogeneous drug-resistant bacteria with rrs mutation and limited probes. Resistance mechanisms other than genetic mutations will affect the consistency of MeltPro and WGS with phenotypic drug-susceptibility results.

Efflux pump gene single-nucleotide variants associated with resistance in Mycobacterium tuberculosis isolates with discrepant drug genotypes

INTRODUCTION

Single-nucleotide variants (SNVs) in Mycobacterium tuberculosis (M. tuberculosis) genomes can predict multidrug resistance (MDR) but not all phenotype-genotype correlations can be explained. We investigated SNVs in efflux pumps (EPs) in the context of M. tuberculosis drug resistance.

METHODS

We analysed 2221 M. tuberculosis genomes from 1432 susceptible and 200 MDR, 172 pre-extensively drug resistant (XDR) and 417 XDR isolates. Analysis of 47 EP genes was conducted using MTB-VCF, an in-house bioinformatics pipeline. SNVs were categorized according to their SIFT/Polyphen scores. Resistance genotypes were also called using the TB-Profiler tool.

RESULTS

Genome comparisons between susceptible and drug resistant (DR) isolates identified 418 unique SNVs in EP of which; 53.5% were in MDR, 68.9% in pre-XDR and 61.3% in XDR isolates. Twenty EPs had unique SNVs with a high SIFT/PolyPhen score, comprising 38 unique SNVs. Sixteen SNVs across 12 EP genes were significantly associated with drug resistance and enriched in pre-XDR and XDR strains. These comprised 12 previously reported SNVs (in Rv0191, Rv0507, Rv0676, Rv1217, Rv1218, Rv1273, Rv1458, Rv1819, and Rv2688) and 4 novel SNVs (in Rv1877 and Rv2333). We investigated their presence in genomes of 52 MDR isolates with phenotype-genotype discrepancies to rifampicin (RIF), isoniazid (INH), or fluoroquinolones. SNVs associated with RIF and INH (Rv1217_1218, Rv1819, Rv0450, Rv1458, Rv3827, Rv0507, Rv0676, Rv1273, and Rv2333), and with fluoroquinolone (Rv2688) resistance were present in these discrepant strains.

CONCLUSIONS

Considering SNVs in EPs as part of M. tuberculosis genome-based resistance interpretation may add value, especially in evaluation of XDR resistance in strains with phenotype-genotype discrepancies.

Sequencing by Binding rivals error-corrected Sequencing by Synthesis technology for accurate detection and quantification of minor (<0.1%) subpopulation variants

Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications including detection of drug resistant pathogens and somatic variant detection in oncology. To enable these applications, wet lab enhancements and bioinformatic error correction methods have been developed for 'sequencing by synthesis' technology to reduce its inherent sequencing error rate. A recently available sequencing approach termed 'sequencing by binding' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using 'sequencing by binding' for the detection of ultra-rare subpopulations down to 0.001%.

Xpert MTB/XDR assay: rapid TB drug resistance detection

PURPOSE

To assess the Xpert MTB/XDR assay's efficiency in promptly detecting resistance to isoniazid, fluoroquinolones, ethionamide, and second-line injectable drugs among tuberculosis (TB) patients.

METHODS

From August 2020 to July 2021, TB suspected patient samples were enrolled at a tertiary care center for our study. We conducted mycobacterial culture, phenotypic DST using proportion method in liquid culture at WHO-recommended concentrations, and the line probe assay (LPA). Simultaneously, the Index test, Xpert MTB/XDR, was performed following the manufacturer's instructions.

RESULTS

Among 360 samples, 107 were excluded due to incomplete information. Resistance to isoniazid, levofloxacin and moxifloxacin was found in 45/251, 21/251 and 20/251 samples, respectively by phenotypic DST. The diagnostic accuracy of Index test, taking phenotypic DST as a reference standard, was 95.8%, 99.04%, and 99.05% for isoniazid, levofloxacin, and moxifloxacin, respectively. The Index test assay demonstrated a specificity of 99.1% for detecting SLID resistance, yielding a diagnostic accuracy of 99.2. Comparing the Index test with LPA revealed a significant enhancement in sensitivity for detecting isoniazid resistance (86.7% vs. 82.2%).

CONCLUSIONS

The Index test exhibited promising outcomes in identifying resistance to isoniazid and fluoroquinolones, surpassing the performance of the LPA. This could be valuable for promptly initiating treatment in cases of drug-resistant tuberculosis.

A robust computational quest: Discovering potential hits to improve the treatment of pyrazinamide-resistant Mycobacterium tuberculosis

The rise of pyrazinamide (PZA)-resistant strains of Mycobacterium tuberculosis (MTB) poses a major challenge to conventional tuberculosis (TB) treatments. PZA, a cornerstone of TB therapy, must be activated by the mycobacterial enzyme pyrazinamidase (PZase) to convert its active form, pyrazinoic acid, which targets the ribosomal protein S1. Resistance, often associated with mutations in the RpsA protein, complicates treatment and highlights a critical gap in the understanding of structural dynamics and mechanisms of resistance, particularly in the context of the G97D mutation. This study utilizes a novel integration of computational techniques, including multiscale biomolecular and molecular dynamics simulations, physicochemical and medicinal chemistry predictions, quantum computations and virtual screening from the ZINC and Chembridge databases, to elucidate the resistance mechanism and identify lead compounds that have the potential to improve treatment outcomes for PZA-resistant MTB, namely ZINC15913786, ZINC20735155, Chem10269711, Chem10279789 and Chem10295790. These computational methods offer a cost-effective, rapid alternative to traditional drug trials by bypassing the need for organic subjects while providing highly accurate insight into the binding sites and efficacy of new drug candidates. The need for rapid and appropriate drug development emphasizes the need for robust computational analysis to justify further validation through in vitro and in vivo experiments.

Performance of the BD MAX MDR-TB assay in a clinical setting and its impact on the clinical course of patients with pulmonary tuberculosis: a retrospective before-after study

BACKGROUND

Missing isoniazid (INH) resistance during tuberculosis (TB) diagnosis can worsen the outcomes of INH-resistant TB. The BD MAX MDR-TB assay (BD MAX) facilitates the rapid detection of TB and INH and rifampin (RIF) resistance; however, data related to its performance in clinical setting remain limited. Moreover, its effect on treatment outcomes has not yet been studied.

METHODS

We compared the performance of BD MAX for the detection of INH/RIF resistances to that of the line probe assay (LPA) in patients with pulmonary TB (PTB), using the results of a phenotypic drug sensitivity test as a reference standard. The treatment outcomes of patients who used BD MAX were compared with those of patients who did not.

RESULTS

Of the 83 patients included in the study, the BD MAX was used for an initial PTB diagnosis in 39 patients. The sensitivity of BD MAX for detecting PTB was 79.5%. The sensitivity and specificity of BD MAX for INH resistance were both 100%, whereas these were 50.0% and 95.8%, respectively, for RIF resistance. The sensitivity and specificity of BD MAX were comparable to those of LPA. The BD MAX group had a shorter time interval from specimen request to the initiation of anti-TB drugs (2.0 days vs. 5.5 days, p=0.001).

CONCLUSION

BD MAX showed comparable performance to conventional tests for detecting PTB and INH/RIF resistances. The implementation of BD MAX as a diagnostic tool for PTB resulted in a shorter turnaround time for the initiation of PTB treatment.

Discordance in Phenotypic and Genotypic Susceptibility Testing for Streptomycin Due to Nonsynonymous/Nonsense/Deletion Frame-Shift Mutations in Gidb Among Clinical Mycobacterium tuberculosis Isolates in Kuwait

OBJECTIVE

Increasing reports of resistance to newer anti-tuberculosis drugs have prompted the search for other alternative drugs. Streptomycin could be used for the treatment of drug-resistant tuberculosis if susceptibility of Mycobacterium tuberculosis isolate to streptomycin could be accurately detected. We performed phenotypic and genotypic drug susceptibility testing (DST) of 118 M. tuberculosis isolates for streptomycin.

MATERIALS AND METHODS

Fifty pansusceptible and 68 multidrug-resistant M. tuberculosis (MDR-TB) isolates were used. Phenotypic DST for streptomycin, rifampicin, isoniazid and ethambutol was performed by mycobacteria growth indicator tube (MGIT) 960 System. Genotypic DST was done by GenoTypeMTBDRplus assay for rifampicin and isoniazid and by PCR-sequencing of rpsL, rrs and gidB genes for streptomycin. MDR-TB isolates were genotyped by spoligotyping.

RESULTS

Phenotypic DST identified 50 isolates susceptible to all four drugs (pansusceptible). Sixty-one of 68 MDR-TB isolates were resistant to streptomycin. Genotypic testing for rifampicin and isoniazid yielded expected results. Fifty pansusceptible and 7 streptomycin-susceptible MDR-TB isolates contained no mutation in rpsL or rrs, while 47, 2 and 1 STR-resistant isolate contained rpsL, rrs and rpsL + rrs mutations, respectively. Of the remaining 11 STR-resistant MDR-TB, 9 isolates contained deletion frame-shift/nonsynonymous mutations in gidB. Surprisingly, 13 pansusceptible isolates also contained deletion frame-shift/nonsense/nonsynonymous mutations in gidB. Also, 30 of 68 MDR-TB but only 2 of 50 pansusceptible isolates belonged to the Beijing genotype.

CONCLUSIONS

Our data show that, like ifampicin, ethambutol and pyrazinamide, streptomycin also exhibits discordant phenotypic and genotypic DST results for some M. tuberculosis isolates. Hence, streptomycin should be included in therapy regimens only if both phenotypic and genotypic resistance testing indicate susceptibility to avoid amplification of resistance and drug toxicity.

Drug-resistant tuberculosis: a persistent global health concern

Drug-resistant tuberculosis (TB) is estimated to cause 13% of all antimicrobial resistance-attributable deaths worldwide and is driven by both ongoing resistance acquisition and person-to-person transmission. Poor outcomes are exacerbated by late diagnosis and inadequate access to effective treatment. Advances in rapid molecular testing have recently improved the diagnosis of TB and drug resistance. Next-generation sequencing of Mycobacterium tuberculosis has increased our understanding of genetic resistance mechanisms and can now detect mutations associated with resistance phenotypes. All-oral, shorter drug regimens that can achieve high cure rates of drug-resistant TB within 6-9 months are now available and recommended but have yet to be scaled to global clinical use. Promising regimens for the prevention of drug-resistant TB among high-risk contacts are supported by early clinical trial data but final results are pending. A person-centred approach is crucial in managing drug-resistant TB to reduce the risk of poor treatment outcomes, side effects, stigma and mental health burden associated with the diagnosis. In this Review, we describe current surveillance of drug-resistant TB and the causes, risk factors and determinants of drug resistance as well as the stigma and mental health considerations associated with it. We discuss recent advances in diagnostics and drug-susceptibility testing and outline the progress in developing better treatment and preventive therapies.

Whole genome sequencing for the prediction of resistant tuberculosis strains from northern India

PURPOSE

Tuberculosis is an important public health problem among infectious diseases. The problem becomes more concerning with the emergence of MDR-TB and pre-XDR-TB. Whole genome sequencing (WGS) detection of resistance has recently gained popularity as it has advantages over other commercial techniques.

METHODS

We performed in-house WGS followed by detailed analysis by an in-house pipeline to identify the resistance markers. This was accompanied by Phenotypic DST, and Sanger sequencing on all the 12 XDR, 06 pre-XDR, and 06 susceptible M. tb isolates. These results were collated with online M. tb WGS pipelines (TB profiler, PhyResSE, Mykrobe predictor) for comparative analysis.

RESULTS

Following our in-house analysis, we observed 64 non-synonymous SNPs, fifteen synonymous SNPs, and five INDELs in 25 drug resistance-associated genes/intergenic regions (IGRs) in M. tb isolates. Sensitivity for detecting XDR is 33%, 58%, 83%, and 83%, respectively, using Mykrobe predictor, PhyResSE, TB-profiler, and in-house pipeline for WGS analysis, respectively. TB-profiler detected a rare mutation H70R in the gyrA gene in one pre-XDR isolate. Lineage 2.2.1 East-Asian (Beijing sublineage type) predominated (60%) in WGS data analysis of the XDR isolates.

CONCLUSIONS

Our findings suggest that in-house analysis of WGS data and TB-profiler sensitivity was better for the detection of second-line resistance as compared to other automated tested tools. Frequent upgradation of newer mutations associated with resistance needs to be updated, as it potentiates tailored treatment for patients.

Opportunities and limitations of genomics for diagnosing bedaquiline-resistant tuberculosis: a systematic review and individual isolate meta-analysis

BACKGROUND

Clinical bedaquiline resistance predominantly involves mutations in mmpR5 (Rv0678). However, mmpR5 resistance-associated variants (RAVs) have a variable relationship with phenotypic Mycobacterium tuberculosis resistance. We did a systematic review to assess the maximal sensitivity of sequencing bedaquiline resistance-associated genes and evaluate the association between RAVs and phenotypic resistance, using traditional and machine-based learning techniques.

METHODS

We screened public databases for articles published from database inception until Oct 31, 2022. Eligible studies performed sequencing of at least mmpR5 and atpE on clinically sourced M tuberculosis isolates and measured bedaquiline minimum inhibitory concentrations (MICs). A bias risk scoring tool was used to identify bias. Individual genetic mutations and corresponding MICs were aggregated, and odds ratios calculated to determine association of mutations with resistance. Machine-based learning methods were used to define test characteristics of parsimonious sets of diagnostic RAVs, and mmpR5 mutations were mapped to the protein structure to highlight mechanisms of resistance. This study was registered in the PROSPERO database (CRD42022346547).

FINDINGS

18 eligible studies were identified, comprising 975 M tuberculosis isolates containing at least one potential RAV (mutation in mmpR5, atpE, atpB, or pepQ), with 201 (20·6%) showing phenotypic bedaquiline resistance. 84 (29·5%) of 285 resistant isolates had no candidate gene mutation. Sensitivity and positive predictive value of taking an any mutation approach was 69% and 14%, respectively. 13 mutations, all in mmpR5, had a significant association with a resistant MIC (adjusted p<0·05). Gradient-boosted machine classifier models for predicting intermediate or resistant and resistant phenotypes both had receiver operator characteristic c statistic of 0·73 (95% CI 0·70-0·76). Frameshift mutations clustered in the α1 helix DNA-binding domain, and substitutions in the α2 and α3 helix hinge region and in the α4 helix-binding domain.

INTERPRETATION

Sequencing candidate genes is insufficiently sensitive to diagnose clinical bedaquiline resistance, but where identified, some mutations should be assumed to be associated with resistance. Genomic tools are most likely to be effective in combination with rapid phenotypic diagnostics. This study was limited by selective sampling in contributing studies and only considering single genetic loci as causative of resistance.

FUNDING

Francis Crick Institute and National Institute of Allergy and Infectious Diseases at the National Institutes of Health.

Phenotype versus genotype discordant rifampicin susceptibility testing in tuberculosis: implications for a diagnostic accuracy

IMPORTANCE

An accurate diagnosis of drug resistance in clinical isolates is an important step for better treatment outcomes. The current study observed a higher discordance rate of rifampicin resistance on Mycobacteria Growth Indicator Tube (MGIT) drug susceptibility testing (DST) than Lowenstein-Jenson (LJ) DST when compared with the rpoB sequencing. We detected a few novel mutations and their combination in rifampicin resistance isolates that were missed by MGIT DST and may be useful for the better management of tuberculosis (TB) treatment outcomes. Few novel deletions in clinical isolates necessitate the importance of rpoB sequencing in large data sets in geographic-specific locations, especially high-burden countries. We explored the discordance rate on MGIT and LJ, which is important for the clinical management of rifampicin resistance to avoid the mistreatment of drug-resistant TB. Furthermore, MGIT-sensitive isolates may be subjected to molecular methods of diagnosis for further confirmation and treatment options.

Role of the first WHO mutation catalogue in the diagnosis of antibiotic resistance in Mycobacterium tuberculosis in the Valencia Region, Spain: a retrospective genomic analysis

BACKGROUND

In June, 2021, WHO published the most complete catalogue to date of resistance-conferring mutations in Mycobacterium tuberculosis. Here, we aimed to assess the performance of genome-based antimicrobial resistance prediction using the catalogue and its potential for improving diagnostics in a real low-burden setting.

METHODS

In this retrospective population-based genomic study M tuberculosis isolates were collected from 25 clinical laboratories in the low-burden setting of the Valencia Region, Spain. Culture-positive tuberculosis cases reported by regional public health authorities between Jan 1, 2014, and Dec 31, 2016, were included. The drug resistance profiles of these isolates were predicted by the genomic identification, via whole-genome sequencing (WGS), of the high-confidence resistance-causing variants included in the catalogue and compared with the phenotype. We determined the minimum inhibitory concentration (MIC) of the isolates with discordant resistance profiles using the resazurin microtitre assay.

FINDINGS

WGS was performed on 785 M tuberculosis complex culture-positive isolates, and the WGS resistance prediction sensitivities were: 85·4% (95% CI 70·8-94·4) for isoniazid, 73·3% (44·9-92·2) for rifampicin, 50·0% (21·1-78·9) for ethambutol, and 57·1% (34·0-78·2) for pyrazinamide; all specificities were more than 99·6%. Sensitivity values were lower than previously reported, but the overall pan-susceptibility accuracy was 96·4%. Genotypic analysis revealed that four phenotypically susceptible isolates carried mutations (rpoB Leu430Pro and rpoB Ile491Phe for rifampicin and fabG1 Leu203Leu for isoniazid) known to give borderline resistance in standard phenotypic tests. Additionally, we identified three putative resistance-associated mutations (inhA Ser94Ala, katG Leu48Pro, and katG Gly273Arg for isoniazid) in samples with substantially higher MICs than those of susceptible isolates. Combining both genomic and phenotypic data, in accordance with the WHO diagnostic guidelines, we could detect two new multidrug-resistant cases. Additionally, we detected 11 (1·6%) of 706 isolates to be monoresistant to fluoroquinolone, which had been previously undetected.

INTERPRETATION

We showed that the WHO catalogue enables the detection of resistant cases missed in phenotypic testing in a low-burden region, thus allowing for better patient-tailored treatment. We also identified mutations not included in the catalogue, relevant at the local level. Evidence from this study, together with future updates of the catalogue, will probably lead in the future to the partial replacement of culture testing with WGS-based drug susceptibility testing in our setting.

FUNDING

European Research Council and the Spanish Ministerio de Ciencia.

Use of Whole Genome Sequencing for Mycobacterium tuberculosis Complex Antimicrobial Susceptibility Testing: From Sequence Data to Resistance Profiles

Whole genome sequencing of Mycobacterium tuberculosis complex (MTBC) isolates has been shown to provide accurate predictions for resistance and susceptibility for many first- and second-line anti-tuberculosis drugs. However, bioinformatic pipelines and mutation catalogs to predict antimicrobial resistances in MTBC isolates are often customized and detailed protocols are difficult to access. Here, we provide a step-by-step workflow for the processing and interpretation of short-read sequencing data and give an overview of available analysis pipelines.

Isoniazid resistance-conferring mutations are associated with highly variable phenotypic resistance

Background

High-dose isoniazid is recommended in the 9-12 months short-course regimen for multidrug-resistant tuberculosis with inhA mutation. However, there is insufficient evidence to support the assumption of genotypic-phenotypic concordance. This study aimed to identify the genetic mutations associated with high-level phenotypic isoniazid resistance.

Methods

Clinical isolates from patients with drug-resistant tuberculosis were profiled by whole-genome sequencing and subjected to minimum inhibitory concentration (MIC) testing using MGIT based-method. MICs were performed in concentration ranges based on the mutation present: isolates with no isoniazid resistance-conferring mutations and H37Rv, 0.016-0.256 µg/ml; inhA, 0.256-4.0 µg/ml, katG 1.0-16.0 µg/ml; and inhA + katG, 4.0-64.0 µg/ml. Isolates demonstrating resistance at the upper limit of the concentration range were tested up to the maximum of 64.0 µg/ml. Bootstrap of the mean MICs was performed to increase the robustness of the estimates and an overlap index was used to compare the distributions of the MICs for each mutation profile.

Results

A total of 52 clinical isolates were included in this analysis. Bootstrap MIC means for inhA, katG and inhA + katG were 33.64 (95% CI, 9.47, 56.90), 6.79 (4.45, 9.70) and 52.34 (42.750, 61.66) µg/ml, respectively. There was high overlap between inhA and inhA + katG mutations (eta = 0.45) but not with inhA and katG (eta = 0.19). Furthermore, katG showed poor overlap with inhA + katG mutations (eta = 0.09). Unexpectedly, 4/8 (50.0%) of all InhA mutants demonstrated high-level resistance, while 20/24 (83.3%) of katG mutants demonstrated moderate-level resistance.

Conclusions

InhA mutations demonstrated unexpectedly high MICs and showed high overlap with inhA + katG. Contrary to the common belief that katG mutants are associated with high-level resistance, this mutation primarily showed moderate-level resistance.

Implementation of national whole-genome sequencing of Mycobacterium tuberculosis, National Public Health Laboratory, Singapore, 2019-2022

The National Tuberculosis Programme (NTBP) monitors the occurrence and spread of tuberculosis (TB) and multidrug-resistant TB (MDR-TB) in Singapore. Since 2020, whole-genome sequencing (WGS) of Mycobacterium tuberculosis isolates has been performed at the National Public Health Laboratory (NPHL) for genomic surveillance, replacing spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeats analysis (MIRU-VNTR). Four thousand three hundred and seven samples were sequenced from 2014 to January 2023, initially as research projects and later developed into a comprehensive public health surveillance programme. Currently, all newly diagnosed culture-positive cases of TB in Singapore are prospectively sent for WGS, which is used to perform lineage classification, predict drug resistance profiles and infer genetic relationships between TB isolates. This paper describes NPHL's operational and technical experiences with implementing the WGS service in an urban TB-endemic setting, focusing on cluster detection and genomic drug susceptibility testing (DST). Cluster detection: WGS has been used to guide contact tracing by detecting clusters and discovering unknown transmission networks. Examples have been clusters in a daycare centre, housing apartment blocks and a horse-racing betting centre. Genomic DST: genomic DST prediction (gDST) identifies mutations in core genes known to be associated with TB drug resistance catalogued in the TBProfiler drug resistance mutation database. Mutations are reported with confidence scores according to a standardized approach referencing NPHL's internal gDST confidence database, which is adapted from the World Health Organization (WHO) TB drug mutation catalogue. Phenotypic-genomic concordance was observed for the first-line drugs ranging from 2959/2998 (98.7 %) (ethambutol) to 2983/2996 (99.6 %) (rifampicin). Aspects of internal database management, reporting standards and caveats in results interpretation are discussed.

Trends in the Antibiotic Resistance of Non-Tuberculous Mycobacteria in Iran: A Systematic Review and Meta-Analysis

Background

Non-tuberculous mycobacteria (NTM) infections have been continuously increasing as major concerns of public health in Iran. Because innate resistance of NTM species, the treatment of these infections is difficult task, but until now resistance pattern of NTM and suitable regimens are not determined.

Methods

We systematically searched the relevant studies in PubMed, Scopus, and Embase (Until Dec 2022). All statistical analyses were carried out using the statistical package R.

Results

Eleven studies included in the analysis were performed in 6 provinces and investigated 1223 NTM clinical species. The majority of the studies originated in Tehran. Among the first-line anti-TB drugs, almost all NTM species were highly resistant to first-line anti-TB drugs. No significant difference in the isoniazid resistance rate was found in the slow or rapid-growing species and Runyon's classification of NTM isolates. A decreased in the prevalence of ciprofloxacin, clarithromycin, and moxifloxacin resistance were showed in during 2013-2022 years.

Conclusion

Most investigated antibiotics have a minor effect on NTM species and a steady increase of resistance has been seen in last few years then, need more-effective alternative regimens is clear.

Rapid Detection of Extensive Drug Resistance by Xpert MTB/XDR Optimizes Therapeutic Decision-Making in Rifampin-Resistant Tuberculosis Patients

The Xpert MTB/XDR assay met the critical need for etiologic diagnosis of tuberculosis and rifampin resistance in previous studies. However, its benefits in tailoring the treatment regimen and improving the outcome for patients with rifampin-resistant tuberculosis (RR-TB) require further investigation. In this study, the Xpert MTB/XDR assay was used to determine the resistance profile of second-line drugs for RR-TB patients in two registered multicenter clinical trials, TB-TRUST (NCT03867136) and TB-TRUST-plus (NCT04717908), with the aim of testing the efficacy of all-oral shorter regimens in RR-TB patients in China. Patients would receive the fluoroquinolone-based all-oral shorter regimen, the injectable-containing regimen, or the bedaquiline-based regimen depending on fluoroquinolone susceptibility by using Xpert MTB/XDR. Among the 497 patients performed with Xpert MTB/XDR, 128 (25.8%) had infections resistant to fluoroquinolones and/or second-line injectable drugs (SLIDs). A total of 371 participants were recruited for the trials, and whole-genome sequencing (WGS) was performed on all corresponding culture-positive baseline strains. Taking the WGS results as the standard, the accuracy of the Xpert MTB/XDR assay in terms of resistance detection was 95.2% to 99.0% for all drugs. A total of 33 cases had inconsistent results, 9 of which were due to resistance heterogeneity. Most of the patients (241/281, 85.8%) had sputum culture conversion at 2 months. In conclusion, the Xpert MTB/XDR assay has the potential to serve as a quick reflex test in patients with RR-TB, as detected via Xpert MTB/RIF, to provide a reliable drug susceptibility profile of the infecting Mycobacterium tuberculosis strain and to initiate optimized treatment promptly.

A Bayesian approach to estimate the probability of resistance to bedaquiline in the presence of a genomic variant

BACKGROUND

Bedaquiline is a core drug for treatment of rifampicin-resistant tuberculosis. Few genomic variants have been statistically associated with bedaquiline resistance. Alternative approaches for determining the genotypic-phenotypic association are needed to guide clinical care.

METHODS

Using published phenotype data for variants in Rv0678, atpE, pepQ and Rv1979c genes in 756 Mycobacterium tuberculosis isolates and survey data of the opinion of 33 experts, we applied Bayesian methods to estimate the posterior probability of bedaquiline resistance and corresponding 95% credible intervals.

RESULTS

Experts agreed on the role of Rv0678, and atpE, were uncertain about the role of pepQ and Rv1979c variants and overestimated the probability of bedaquiline resistance for most variant types, resulting in lower posterior probabilities compared to prior estimates. The posterior median probability of bedaquiline resistance was low for synonymous mutations in atpE (0.1%) and Rv0678 (3.3%), high for missense mutations in atpE (60.8%) and nonsense mutations in Rv0678 (55.1%), relatively low for missense (31.5%) mutations and frameshift (30.0%) in Rv0678 and low for missense mutations in pepQ (2.6%) and Rv1979c (2.9%), but 95% credible intervals were wide.

CONCLUSIONS

Bayesian probability estimates of bedaquiline resistance given the presence of a specific mutation could be useful for clinical decision-making as it presents interpretable probabilities compared to standard odds ratios. For a newly emerging variant, the probability of resistance for the variant type and gene can still be used to guide clinical decision-making. Future studies should investigate the feasibility of using Bayesian probabilities for bedaquiline resistance in clinical practice.

Opportunities and limitations of genomics for diagnosing bedaquiline-resistant tuberculosis: an individual isolate metaanalysis

Background

Clinical bedaquiline resistance predominantly involves mutations in mmpR5 (Rv0678). However, mmpR5 resistance-associated variants (RAVs) have a variable relationship with phenotypic M. tuberculosis resistance. We performed a systematic review to (1) assess the maximal sensitivity of sequencing bedaquiline resistance-associated genes and (2) evaluate the association between RAVs and phenotypic resistance, using traditional and machine-based learning techniques.

Methods

We screened public databases for articles published until October 2022. Eligible studies performed sequencing of at least mmpR5 and atpE on clinically-sourced M. tuberculosis isolates and measured bedaquiline minimum inhibitory concentrations (MICs). We performed genetic analysis for identification of phenotypic resistance and determined the association of RAVs with resistance. Machine-based learning methods were employed to define test characteristics of optimised sets of RAVs, and mmpR5 mutations were mapped to the protein structure to highlight mechanisms of resistance.

Results

Eighteen eligible studies were identified, comprising 975 M. tuberculosis isolates containing ≥1 potential RAV (mutation in mmpR5, atpE, atpB or pepQ), with 201 (20.6%) demonstrating phenotypic bedaquiline resistance. 84/285 (29.5%) resistant isolates had no candidate gene mutation. Sensitivity and positive predictive value of taking an 'any mutation' approach was 69% and 14% respectively. Thirteen mutations, all in mmpR5, had a significant association with a resistant MIC (adjusted p<0.05). Gradient-boosted machine classifier models for predicting intermediate/resistant and resistant phenotypes both had receiver operator characteristic c-statistics of 0.73. Frameshift mutations clustered in the alpha 1 helix DNA binding domain, and substitutions in the alpha 2 and 3 helix hinge region and in the alpha 4 helix binding domain.

Discussion

Sequencing candidate genes is insufficiently sensitive to diagnose clinical bedaquiline resistance, but where identified a limited number of mutations should be assumed to be associated with resistance. Genomic tools are most likely to be effective in combination with rapid phenotypic diagnostics.

Whole-genome sequencing and transcriptome-characterized in vitro evolution of aminoglycoside resistance in Mycobacterium tuberculosis

Antibiotic resistance of Mycobacterium tuberculosis (Mtb) is a major public health concern worldwide. Therefore, it is of great significance to characterize the mutational pathways by which susceptible Mtb evolves into drug resistance. In this study, we used laboratory evolution to explore the mutational pathways of aminoglycoside resistance. The level of resistance in amikacin inducing Mtb was also associated with changes in susceptibility to other anti-tuberculosis drugs such as isoniazid, levofloxacin and capreomycin. Whole-genome sequencing (WGS) revealed that the induced resistant Mtb strains had accumulated diverse mutations. We found that rrs A1401G was the predominant mutation in aminoglycoside-resistant clinical Mtb isolates from Guangdong. In addition, this study provided global insight into the characteristics of the transcriptome in four representative induced strains and revealed that rrs mutated and unmutated aminoglycoside-resistant Mtb strains have different transcriptional profiles. WGS analysis and transcriptional profiling of Mtb strains during evolution revealed that Mtb strains harbouring rrs A1401G have an evolutionary advantage over other drug-resistant strains under the pressure of aminoglycosides because of their ultra-high resistance level and low physiological impact on the strain. The results of this study should advance our understanding of aminoglycoside resistance mechanisms.

Clinical implications of molecular drug resistance testing for Mycobacterium tuberculosis: a 2023 TBnet/RESIST-TB consensus statement

Drug-resistant tuberculosis is a substantial health-care concern worldwide. Despite culture-based methods being considered the gold standard for drug susceptibility testing, molecular methods provide rapid information about the Mycobacterium tuberculosis mutations associated with resistance to anti-tuberculosis drugs. This consensus document was developed on the basis of a comprehensive literature search, by the TBnet and RESIST-TB networks, about reporting standards for the clinical use of molecular drug susceptibility testing. Review and the search for evidence included hand-searching journals and searching electronic databases. The panel identified studies that linked mutations in genomic regions of M tuberculosis with treatment outcome data. Implementation of molecular testing for the prediction of drug resistance in M tuberculosis is key. Detection of mutations in clinical isolates has implications for the clinical management of patients with multidrug-resistant or rifampicin-resistant tuberculosis, especially in situations when phenotypic drug susceptibility testing is not available. A multidisciplinary team including clinicians, microbiologists, and laboratory scientists reached a consensus on key questions relevant to molecular prediction of drug susceptibility or resistance to M tuberculosis, and their implications for clinical practice. This consensus document should help clinicians in the management of patients with tuberculosis, providing guidance for the design of treatment regimens and optimising outcomes.

Drug-Resistant Tuberculosis on the Balkan Peninsula: Determination of Drug Resistance Mechanisms with Xpert MTB/XDR and Whole-Genome Sequencing Analysis

The new molecular assay Xpert MTB/XDR (Cepheid, Sunnyvale, CA, USA) was launched in 2021 to detect Mycobacterium tuberculosis (MT) complex with mutations conferring resistance to isoniazid (INH), ethionamide (ETH), fluoroquinolone (FQ), and second-line injectable drugs (SLIDs). The aim of our study was to evaluate the performance of the Xpert MTB/XDR rapid molecular assay on rifampicin-resistant, multidrug-resistant, and pre-extensively resistant tuberculosis (TB) isolates in a clinical laboratory in the Balkan Peninsula compared to a phenotypic drug susceptibility test (pDST). Xpert MTB/XDR was used to test positive Bactec MGIT 960 (Becton, Dickinson and Co., Franklin Lakes, NJ, USA) cultures or DNA isolates. In the case of discrepant results between Xpert MTB/XDR and pDST, the usefulness of whole-genome sequencing (WGS) was emphasized. In our study, 80 MT isolates from different Balkan countries were selectively chosen from the National Mycobacterial Strain Collection in Golnik, Slovenia. Isolates were tested with the Xpert MTB/XDR assay, conventional pDST, and WGS. Xpert MTB/XDR showed high sensitivities of 91.9%, 100%, and 100% for detecting INH, FQ, and SLID resistance, respectively, compared to pDST. In contrast, low sensitivity (51.9%) for ETH resistance was achieved because isolates harbored widespread mutations across the ethA gene. The specificity of Xpert MTB/XDR was 100% for all drugs except for INH (66.7%). Further investigation with WGS revealed -57c→t mutations in the oxyR-ahpC region marked with uncertain significance, which caused the low specificity for detecting INH resistance with the new assay. Xpert MTB/XDR can be used in clinical laboratories for the rapid detection of INH, FQ, and SLID resistance. Moreover, it can be used to rule in resistance to ETH. Additional use of WGS is recommended in cases of discrepant results between pDST and Xpert MTB/XDR. Future improvements of Xpert MTB/XDR with the inclusion of additional genes may increase the usefulness of the assay. IMPORTANCE The Xpert MTB/XDR was tested on drug-resistant Mycobacterium tuberculosis complex isolates from the Balkan Peninsula. Positive Bactec MGIT 960 cultures or DNA isolates were tested as starting material. According to the results of our study with Xpert MTB/XDR, sensitivities for the detection of SLID, FQ, and INH resistance were sufficient (>90%) for the assay to be implemented into diagnostic algorithms. In our study, WGS revealed lesser-known mutations in genes conferring INH and ETH resistance, and their impact on resistance is still unknown. Mutations in the ethA gene causing resistance to ETH were scattered along structural gene without high-confidence markers for resistance. Therefore, resistance to ETH should be reported based on a combination of methods. Because the Xpert MTB/XDR assay was found to have good performance, we propose that it should be the method of choice for confirming resistance to INH, FQ, and SLID and conditionally for resistance to ETH.

Experiences from 4 Years of Organization of an External Quality Assessment for Mycobacterium tuberculosis Whole-Genome Sequencing in the European Union/European Economic Area

Here, we report the development and key features of the first external quality assessment (EQA) scheme for Mycobacterium tuberculosis whole-genome sequencing (WGS). The results of four rounds (2017 to 2020) of implementation within the European tuberculosis reference laboratories network (ERLTB-Net-2) are presented and discussed. EQA panels comprising 10 genomic DNAs were distributed to ERLTB-Net 2 laboratories volunteering to participate in this exercise. Since 2018, five FASTQ files were added to better assess the dry WGS processes, and in 2020, three of the five files were replaced by synthetic files (providing additional flexibility for the mutations included in the panels). Ten National tuberculosis reference laboratories participated in all four EQA rounds, and seven participated in at least one. High-confidence resistance mutations were correctly identified by all laboratories, but challenges remained with respect to the identification of mixed loci and interpretation of rare mutations. M. tuberculosis genotyping and clustering analysis was >90% accurate for pure samples with the main challenges being related to the analysis of mixed genotypes and DNA FASTQ files. The development and implementation of this WGS EQA scheme has contributed to the continuous improvement in performance of participating laboratories in M. tuberculosis WGS and data analysis. This scheme can serve as a model of comprehensive quality assessment for M. tuberculosis WGS that can be replicated in different settings worldwide. IMPORTANCE The wider availability of whole-genome sequencing (WGS) coupled to new developments in bioinformatic tools and databases to interpret Mycobacterium tuberculosis complex WGS data has accelerated the adoption of this method for the routine prediction of antimycobacterial drug resistance and genotyping, thus necessitating the establishment of a comprehensive external quality control system. Here, we report 4 years of development and results from such a panel.

Investigation of Multi-Subunit Mycobacterium tuberculosis DNA-Directed RNA Polymerase and Its Rifampicin Resistant Mutants

Emerging Mycobacterium tuberculosis (Mtb) resistant strains have continued to limit the efficacies of existing antitubercular therapies. More specifically, mutations in the RNA replicative machinery of Mtb, RNA polymerase (RNAP), have been widely linked to rifampicin (RIF) resistance, which has led to therapeutic failures in many clinical cases. Moreover, elusive details on the underlying mechanisms of RIF-resistance caused by Mtb-RNAP mutations have hampered the development of new and efficient drugs that are able to overcome this challenge. Therefore, in this study we attempt to resolve the molecular and structural events associated with RIF-resistance in nine clinically reported missense Mtb RNAP mutations. Our study, for the first time, investigated the multi-subunit Mtb RNAP complex and findings revealed that the mutations commonly disrupted structural-dynamical attributes that may be essential for the protein's catalytic functions, particularly at the βfork loop 2, β'zinc-binding domain, the β' trigger loop and β'jaw, which in line with previous experimental reports, are essential for RNAP processivity. Complementarily, the mutations considerably perturbed the RIF-BP, which led to alterations in the active orientation of RIF needed to obstruct RNA extension. Consequentially, essential interactions with RIF were lost due to the mutation-induced repositioning with corresponding reductions in the binding affinity of the drug observed in majority of the mutants. We believe these findings will significantly aid future efforts in the discovery of new treatment options with the potential to overcome antitubercular resistance.

Omics analysis of Mycobacterium tuberculosis isolates uncovers Rv3094c, an ethionamide metabolism-associated gene

Global control of the tuberculosis epidemic is threatened by increasing prevalence of drug resistant M. tuberculosis isolates. Many genome-wide studies focus on SNP-associated drug resistance mechanisms, but drug resistance in 5-30% of M. tuberculosis isolates (varying with antibiotic) appears unrelated to reported SNPs, and alternative drug resistance mechanisms involving variation in gene/protein expression are not well-studied. Here, using an omics approach, we identify 388 genes with lineage-related differential expression and 68 candidate drug resistance-associated gene pairs/clusters in 11 M. tuberculosis isolates (variable lineage/drug resistance profiles). Structural, mutagenesis, biochemical and bioinformatic studies on Rv3094c from the Rv3093c-Rv3095 gene cluster, a gene cluster selected for further investigation as it contains a putative monooxygenase/repressor pair and is associated with ethionamide resistance, provide insights on its involvement in ethionamide sulfoxidation, the initial step in its activation. Analysis of the structure of Rv3094c and its complex with ethionamide and flavin mononucleotide, to the best of our knowledge the first structures of an enzyme involved in ethionamide activation, identify key residues in the flavin mononucleotide and ethionamide binding pockets of Rv3094c, and F221, a gate between flavin mononucleotide and ethionamide allowing their interaction to complete the sulfoxidation reaction. Our work broadens understanding of both lineage- and drug resistance-associated gene/protein expression perturbations and identifies another player in mycobacterial ethionamide metabolism.

Comparative genome analysis reveals high-level drug resistance markers in a clinical isolate of Mycobacterium fortuitum subsp. fortuitum MF GZ001

Introduction

Infections caused by non-tuberculosis mycobacteria are significantly worsening across the globe. M. fortuitum complex is a rapidly growing pathogenic species that is of clinical relevance to both humans and animals. This pathogen has the potential to create adverse effects on human healthcare.

Methods

The MF GZ001 clinical strain was collected from the sputum of a 45-year-old male patient with a pulmonary infection. The morphological studies, comparative genomic analysis, and drug resistance profiles along with variants detection were performed in this study. In addition, comparative analysis of virulence genes led us to understand the pathogenicity of this organism.

Results

Bacterial growth kinetics and morphology confirmed that MF GZ001 is a rapidly growing species with a rough morphotype. The MF GZ001 contains 6413573 bp genome size with 66.18 % high G+C content. MF GZ001 possesses a larger genome than other related mycobacteria and included 6156 protein-coding genes. Molecular phylogenetic tree, collinearity, and comparative genomic analysis suggested that MF GZ001 is a novel member of the M. fortuitum complex. We carried out the drug resistance profile analysis and found single nucleotide polymorphism (SNP) mutations in key drug resistance genes such as rpoB, katG, AAC(2')-Ib, gyrA, gyrB, embB, pncA, blaF, thyA, embC, embR, and iniA. In addition, the MF GZ001strain contains mutations in iniA, iniC, pncA, and ribD which conferred resistance to isoniazid, ethambutol, pyrazinamide, and para-aminosalicylic acid respectively, which are not frequently observed in rapidly growing mycobacteria. A wide variety of predicted putative potential virulence genes were found in MF GZ001, most of which are shared with well-recognized mycobacterial species with high pathogenic profiles such as M. tuberculosis and M. abscessus.

Discussion

Our identified novel features of a pathogenic member of the M. fortuitum complex will provide the foundation for further investigation of mycobacterial pathogenicity and effective treatment.

Clinical utility of target-based next-generation sequencing for drug-resistant TB

BACKGROUND: In high TB burden countries, access to drug susceptibility testing is a major bottleneck. Targeted next-generation sequencing (tNGS) is a promising technology for rapid resistance detection. This study assessed the role of tNGS for the diagnosis of drug-resistant TB (DR-TB).METHODS: A total of 161 samples from bacteriologically confirmed TB cases were subjected to tNGS using the Deeplex^® Myc-TB kit and sequenced using the MiSeq platform. These samples were also processed for conventional phenotypic DST (pDST) using 13 drugs on Mycobacteria Growth Indicator Tube and line-probe assays (MTBDRplus and MTBDRsl).RESULTS: There were 146 DR-TB and 15 drug-susceptible TB (DS-TB) samples. About 70% of patients with DR-TB had no previous TB treatment history. Overall, 88.2% had rifampicin-resistant/multidrug-resistant TB (RR/MDR-TB), 58.5% pre-extensively drug-resistant TB (pre-XDR-TB) and 9.2% had XDR-TB as defined by the WHO (2020). Around 8% (n = 13) of samples were non-culturable; however, identified 8 were resistant to first and second-line drugs using tNGS. Resistance frequency was similar across methods, with discordance in drugs less reliable using pDST or with limited mutational representation within databases. Sensitivities were aligned with literature reports for most drugs. We observed 10% heteroresistance, while 75% of strains were of Lineages 2 and 3.CONCLUSIONS: Programme data supported tNGS in the diagnosis of DR-TB for early treatment using individualised regimens.

Genomic analysis of Mycobacterium tuberculosis variant bovis strains isolated from bovine in the state of Mato Grosso, Brazil

The species Mycobacterium tuberculosis variant bovis (M. tuberculosis var. bovis) is associated with tuberculosis, mainly in cattle and buffaloes. This pathogen has the potential to infect other mammals, including humans. Tuberculosis caused by M. tuberculosis var. bovis is a zoonosis clinically identical to tuberculosis caused by Mycobacterium tuberculosis, and the recommended treatment in humans results in the use of antibiotics. In this study, we used the whole genome sequencing (WGS) methodology Illumina NovaSeq 6000 System platform to characterize the genome of M. tuberculosis var. bovis in cattle circulating in Mato Grosso, identify mutations related to drug resistance genes, compare with other strains of M. tuberculosis var. bovis brazilian and assess potential drug resistance. Four isolates of M. tuberculosis var. bovis of cattle origin representing the main livestock circuits, which had been more prevalent in previous studies in the state of Mato Grosso, were selected for the genomic study. The genome sizes of the sequenced strains ranged from 4,306,423 to 4,332,964 bp, and the GC content was 65.6%. The four strains from Mato Grosso presented resistance genes to pncA (pyrazinamide), characterized as drug-resistant strains. In addition to verifying several point mutations in the pncA, rpsA, rpsL, gid, rpoB, katG, gyrB, gyrA, tlyA, embA, embB, embC, fgd, fbiB, and fbiC genes, these genes were similar to antibiotic resistance in more than 92% of the Brazilian strains. Therefore, our results indicated a high genetic diversity between our isolates and other M. tuberculosis var. bovis isolated in Brazil. Thus, multiple transmission routes of this pathogen may be present in the production chain. So, to achieve a bovine tuberculosis-free health status, the use of the WGS as a control and monitoring tool will be crucial to determine these transmission routes.

Regional distribution of Mycobacterium tuberculosis infection and resistance to rifampicin and isoniazid as determined by high-resolution melt analysis

Background

Identifying the transmission mode and resistance mechanism of Mycobacterium tuberculosis (MTB) is key to prevent disease transmission. However, there is a lack of regional data. Therefore, the aim of this study was to identify risk factors associated with the transmission of MTB and regional patterns of resistance to isoniazid (INH) and rifampicin (RFP), as well as the prevalence of multidrug-resistant tuberculosis (MDR-TB).

Methods

High-resolution melt (HRM) analysis was conducted using sputum, alveolar lavage fluid, and pleural fluid samples collected from 17,515 patients with suspected or confirmed MTB infection in the downtown area and nine counties of Luoyang City from 2019 to 2021.

Results

Of the 17,515 patients, 82.6% resided in rural areas, and 96.0% appeared for an initial screening. The HRM positivity rate was 16.8%, with a higher rate in males than females (18.0% vs. 14.1%, p < 0.001). As expected, a positive sputum smear was correlated with a positive result for HRM analysis. By age, the highest rates of MTB infection occurred in males (22.9%) aged 26–30 years and females (28.1%) aged 21–25. The rates of resistance to RFP and INH and the incidence of MDR were higher in males than females (20.5% vs. 16.1%, p < 0.001, 15.9% vs. 12.0%, p < 0.001 and 12.9% vs. 10.2%, p < 0.001, respectively). The HRM positivity rate was much higher in previously treated patients than those newly diagnosed for MTB infection. Notably, males at the initial screening had significantly higher rates of HRM positive, INH resistance, RFP resistance, and MDR-TB than females (all, p < 0.05), but not those previously treated for MTB infection. The HRM positivity and drug resistance rates were much higher in the urban vs. rural population. By multivariate analyses, previous treatment, age < 51 years, residing in an urban area, and male sex were significantly and positively associated with drug resistance after adjusting for smear results and year of testing.

Conclusion

Males were at higher risks for MTB infection and drug resistance, while a younger age was associated with MTB infection, resistance to INH and RFP, and MDR-TB. Further comprehensive monitoring of resistance patterns is needed to control the spread of MTB infection and manage drug resistance locally.

Rapid Identification of Drug Resistance and Phylogeny in M. tuberculosis, Directly from Sputum Samples

Tuberculosis (TB) remains one of the most important infectious diseases globally. Establishing a resistance profile from the initial TB diagnosis is a priority. Rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, and Whole Genome Sequencing (WGS) needs culture prior to next-generation sequencing (NGS), limiting their clinical value. Targeted sequencing (TS) from clinical samples avoids these drawbacks, providing a signature of genetic markers that can be associated with drug resistance and phylogeny. In this study, a proof-of-concept protocol was developed for detecting genomic variants associated with drug resistance and for the phylogenetic classification of Mycobacterium Tuberculosis (Mtb) in sputum samples. Initially, a set of Mtb reference strains from the WHO were sequenced (WGS and TS). The results from the protocol agreed >95% with WHO reported data and phenotypic drug susceptibility testing (pDST). Lineage genetics results were 100% concordant with those derived from WGS. After that, the TS protocol was applied to sputum samples from TB patients to detect resistance to first- and second-line drugs and derive phylogeny. The accuracy was >90% for all evaluated drugs, except Eto/Pto (77.8%), and 100% were phylogenetically classified. The results indicate that the described protocol, which affords the complete drug resistance profile and phylogeny of Mtb from sputum, could be useful in the clinical area, advancing toward more personalized and more effective treatments in the near future. IMPORTANCE The COVID-19 pandemic negatively affected the progress in accessing essential Tuberculosis (TB) services and reducing the burden of TB disease, resulting in a decreased detection of new cases and increased deaths. Generating molecular diagnostic tests with faster results without losing reliability is considered a priority. Specifically, developing an antimicrobial resistance profile from the initial stages of TB diagnosis is essential to ensure appropriate treatment. Currently available rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, limiting their clinical value. In this work, targeted sequencing on sputum samples from TB patients was used to identify Mycobacterium tuberculosis mutations in genes associated with drug resistance and to derive a phylogeny of the infecting strain. This protocol constitutes a proof-of-concept toward the goal of helping clinicians select a timely and appropriate treatment by providing them with actionable information beyond current molecular approaches.

In silico evaluation of WHO-endorsed molecular methods to detect drug resistant tuberculosis

Universal drug susceptibility testing (DST) for tuberculosis is a major goal of the END TB strategy. PCR-based molecular diagnostic tests have been instrumental in increasing DST globally and several assays have now been endorsed by the World Health Organization (WHO) for use in the diagnosis of drug resistance. These endorsed assays, however, each interrogate a limited number of mutations associated with resistance, potentially limiting their sensitivity compared to sequencing-based methods. We applied an in silico method to compare the sensitivity and specificity of WHO-endorsed molecular based diagnostics to the mutation set identified by the WHO mutations catalogue using phenotypic DST as the reference. We found that, in silico, the mutation sets used by probe-based molecular diagnostic tests to identify rifampicin, isoniazid, pyrazinamide, levofloxacin, moxifloxacin, amikacin, capreomycin and kanamycin resistance produced similar sensitivities and specificities to the WHO mutation catalogue. PCR-based diagnostic tests were most sensitive for drugs where mechanisms of resistance are well established and localised to small genetic regions or a few prevalent mutations. Approaches using sequencing technologies can provide advantages for drugs where our knowledge of resistance is limited, or where complex resistance signatures exist.

“Upcycling” known molecules and targets for drug-resistant TB

Despite reinvigorated efforts in Tuberculosis (TB) drug discovery over the past 20 years, relatively few new drugs and candidates have emerged with clear utility against drug resistant TB. Over the same period, significant technological advances and learnings around target value have taken place. This has offered opportunities to re-assess the potential for optimization of previously discovered chemical matter against Mycobacterium tuberculosis (M.tb) and for reconsideration of clinically validated targets encumbered by drug resistance. A re-assessment of discarded compounds and programs from the “golden age of antibiotics” has yielded new scaffolds and targets against TB and uncovered classes, for example beta-lactams, with previously unappreciated utility for TB. Leveraging validated classes and targets has also met with success: booster technologies and efforts to thwart efflux have improved the potential of ethionamide and spectinomycin classes. Multiple programs to rescue high value targets while avoiding cross-resistance are making progress. These attempts to make the most of known classes, drugs and targets complement efforts to discover new chemical matter against novel targets, enhancing the chances of success of discovering effective novel regimens against drug-resistant TB.

Predicting antibiotic resistance in complex protein targets using alchemical free energy methods

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.

Epidemiological cut-off values for a 96-well broth microdilution plate for high-throughput research antibiotic susceptibility testing of M. tuberculosis

Drug susceptibility testing of M. tuberculosis is rooted in a binary susceptible/resistant paradigm. While there are considerable advantages in measuring the minimum inhibitory concentrations (MICs) of a panel of drugs for an isolate, it is necessary to measure the epidemiological cut-off values (ECOFF/ECVs) to permit comparison with qualitative data. Here we present ECOFF/ECVs for 13 anti-tuberculosis compounds, including bedaquiline and delamanid, derived from 20 637 clinical isolates collected by 14 laboratories based in 11 countries on five continents. Each isolate was incubated for 14 days on a dry 96-well broth microdilution plate and then read. Resistance to most of the drugs due to prior exposure is expected and the MIC distributions for many of the compounds are complex, and therefore a phenotypically wild-type population could not be defined. Since a majority of samples also underwent genetic sequencing, we defined a genotypically wild-type population and measured the MIC of the 99th percentile by direct measurement and via fitting a Gaussian using interval regression. The proposed ECOFF/ECVs were then validated by comparing with the MIC distributions of high-confidence genetic variants that confer resistance and with qualitative drug susceptibility tests obtained via the Mycobacterial Growth Indicator Tube (MGIT) system or Microscopic-Observation Drug Susceptibility (MODS) assay. These ECOFF/ECVs will inform and encourage the more widespread adoption of broth microdilution: this is a cheap culture-based method that tests the susceptibility of 12-14 antibiotics on a single 96-well plate and so could help personalise the treatment of tuberculosis.

Evaluation of xpert Mycobacterium tuberculosis rifampicin for tuberculosis diagnosis in a reference laboratory

Background

The Xpert Mycobacterium tuberculosis Rifampicin (MTB-RIF) is a technological innovation that presents precision and speed in the diagnosis of tuberculosis (TB). The study aimed to evaluate the performance of Xpert MTB/RIF in the TB diagnosis and compare its results with those of culture, species identification, Antimicrobial Susceptibility Testing (AST), and rpoB gene sequencing of discordant results in AST, used for the diagnosis of TB in a reference laboratory.

Methods

Cross-sectional descriptive study of pulmonary and extrapulmonary samples requesting Xpert MTB/RIF and culture for TB diagnosis from 2015 to 2019 at Adolfo Lutz Institute-São Paulo/Brazil. The analysis was performed with Epi-Info 7.2.1, presenting the distribution of frequencies and, for comparative analyses, Pearson's Chi-square test and Fisher's exact test were used, considering P ≤ 0.05 as statistically significant. For variables agreement, the Kappa method was used.

Results

A total of 1575 pulmonary and extrapulmonary samples were analyzed using Xpert MTB/RIF and culture, of which 293 were positive for the MTB Complex in both methodologies with a sensitivity of 94.55%; specificity of 95.97%; accuracy 95.69%; positive predictive value 85.53%; negative predictive value 98.59%, substantial agreement by Kappa 0.87, and detection sensitivity even at lower levels of bacillary load (P < 0.05). The Xpert MTB/RIF and AST showed concordant results (P < 0.05).

Conclusion

The study brings forward that the Xpert MTB/RIF shows substantial agreement with the results of culture and AST, contributing to the diagnosis of TB and the rapid detection of resistance. The sequencing of resistant cultures in Xpert MTB/RIF and sensitivity in AST identified the H526N mutation, characterized by a low level of resistance to RIF.

Investigating resistance in clinical Mycobacterium tuberculosis complex isolates with genomic and phenotypic antimicrobial susceptibility testing: a multicentre observational study

BACKGROUND

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis complex has become an important tool in diagnosis and management of drug-resistant tuberculosis. However, data correlating resistance genotype with quantitative phenotypic antimicrobial susceptibility testing (AST) are scarce.

METHODS

In a prospective multicentre observational study, 900 clinical M tuberculosis complex isolates were collected from adults with drug-resistant tuberculosis in five high-endemic tuberculosis settings around the world (Georgia, Moldova, Peru, South Africa, and Viet Nam) between Dec 5, 2014, and Dec 12, 2017. Minimum inhibitory concentrations (MICs) and resulting binary phenotypic AST results for up to nine antituberculosis drugs were determined and correlated with resistance-conferring mutations identified by WGS.

FINDINGS

Considering WHO-endorsed critical concentrations as reference, WGS had high accuracy for prediction of resistance to isoniazid (sensitivity 98·8% [95% CI 98·5-99·0]; specificity 96·6% [95% CI 95·2-97·9]), levofloxacin (sensitivity 94·8% [93·3-97·6]; specificity 97·1% [96·7-97·6]), kanamycin (sensitivity 96·1% [95·4-96·8]; specificity 95·0% [94·4-95·7]), amikacin (sensitivity 97·2% [96·4-98·1]; specificity 98·6% [98·3-98·9]), and capreomycin (sensitivity 93·1% [90·0-96·3]; specificity 98·3% [98·0-98·7]). For rifampicin, pyrazinamide, and ethambutol, the specificity of resistance prediction was suboptimal (64·0% [61·0-67·1], 83·8% [81·0-86·5], and 40·1% [37·4-42·9], respectively). Specificity for rifampicin increased to 83·9% when borderline mutations with MICs overlapping with the critical concentration were excluded. Consequently, we highlighted mutations in M tuberculosis complex isolates that are often falsely identified as susceptible by phenotypic AST, and we identified potential novel resistance-conferring mutations.

INTERPRETATION

The combined analysis of mutations and quantitative phenotypes shows the potential of WGS to produce a refined interpretation of resistance, which is needed for individualised therapy, and eventually could allow differential drug dosing. However, variability of MIC data for some M tuberculosis complex isolates carrying identical mutations also reveals limitations of our understanding of the genotype and phenotype relationships (eg, including epistasis and strain genetic background).

FUNDING

Bill & Melinda Gates Foundation, German Centre for Infection Research, German Research Foundation, Excellence Cluster Precision Medicine of Inflammation (EXC 2167), and Leibniz ScienceCampus EvoLUNG.

Evaluation of whole-genome sequence to predict drug resistance of nine anti-tuberculosis drugs and characterize resistance genes in clinical rifampicin-resistant Mycobacterium tuberculosis isolates from Ningbo, China

Setting

Controlling drug-resistant tuberculosis in Ningbo, China.

Objective

Whole-genome sequencing (WGS) has not been employed to comprehensively study Mycobacterium tuberculosis isolates, especially rifampicin-resistant tuberculosis, in Ningbo, China. Here, we aim to characterize genes involved in drug resistance in RR-TB and create a prognostic tool for successfully predicting drug resistance in patients with TB.

Design

Drug resistance was predicted by WGS in a "TB-Profiler" web service after phenotypic drug susceptibility tests (DSTs) against nine anti-TB drugs among 59 clinical isolates. A comparison of consistency, sensitivity, specificity, and positive and negative predictive values between WGS and DST were carried out for each drug.

Results

The sensitivities and specificities for WGS were 95.92 and 90% for isoniazid (INH), 100 and 64.1% for ethambutol (EMB), 97.37 and 100% for streptomycin (SM), 75 and 100% for amikacin (AM), 80 and 96.3%for capreomycin (CAP), 100 and 97.22% for levofloxacin (LFX), 93.33 and 90.91% for prothionamide (PTO), and 70 and 97.96% for para-aminosalicylic acid (PAS). Around 53 (89.83%) and 6 (10.17%) of the isolates belonged to lineage two (East-Asian) and lineage four (Euro-American), respectively.

Conclusion

Whole-genome sequencing is a reliable method for predicting resistance to INH, RIF, EMB, SM, AM, CAP, LFX, PTO, and PAS with high consistency, sensitivity, and specificity. There was no transmission that occurred among the patients with RR-TB in Ningbo, China.

A data compendium associating the genomes of 12,289 Mycobacterium tuberculosis isolates with quantitative resistance phenotypes to 13 antibiotics

The Comprehensive Resistance Prediction for Tuberculosis: an International Consortium (CRyPTIC) presents here a data compendium of 12,289 Mycobacterium tuberculosis global clinical isolates, all of which have undergone whole-genome sequencing and have had their minimum inhibitory concentrations to 13 antitubercular drugs measured in a single assay. It is the largest matched phenotypic and genotypic dataset for M. tuberculosis to date. Here, we provide a summary detailing the breadth of data collected, along with a description of how the isolates were selected, collected, and uniformly processed in CRyPTIC partner laboratories across 23 countries. The compendium contains 6,814 isolates resistant to at least 1 drug, including 2,129 samples that fully satisfy the clinical definitions of rifampicin resistant (RR), multidrug resistant (MDR), pre-extensively drug resistant (pre-XDR), or extensively drug resistant (XDR). The data are enriched for rare resistance-associated variants, and the current limits of genotypic prediction of resistance status (sensitive/resistant) are presented by using a genetic mutation catalogue, along with the presence of suspected resistance-conferring mutations for isolates resistant to the newly introduced drugs bedaquiline, clofazimine, delamanid, and linezolid. Finally, a case study of rifampicin monoresistance demonstrates how this compendium could be used to advance our genetic understanding of rare resistance phenotypes. The data compendium is fully open source and it is hoped that it will facilitate and inspire future research for years to come.

Genome-wide association studies of global Mycobacterium tuberculosis resistance to 13 antimicrobials in 10,228 genomes identify new resistance mechanisms

The emergence of drug-resistant tuberculosis is a major global public health concern that threatens the ability to control the disease. Whole-genome sequencing as a tool to rapidly diagnose resistant infections can transform patient treatment and clinical practice. While resistance mechanisms are well understood for some drugs, there are likely many mechanisms yet to be uncovered, particularly for new and repurposed drugs. We sequenced 10,228 Mycobacterium tuberculosis (MTB) isolates worldwide and determined the minimum inhibitory concentration (MIC) on a grid of 2-fold concentration dilutions for 13 antimicrobials using quantitative microtiter plate assays. We performed oligopeptide- and oligonucleotide-based genome-wide association studies using linear mixed models to discover resistance-conferring mechanisms not currently catalogued. Use of MIC over binary resistance phenotypes increased sample heritability for the new and repurposed drugs by 26% to 37%, increasing our ability to detect novel associations. For all drugs, we discovered uncatalogued variants associated with MIC, including in the Rv1218c promoter binding site of the transcriptional repressor Rv1219c (isoniazid), upstream of the vapBC20 operon that cleaves 23S rRNA (linezolid) and in the region encoding an α-helix lining the active site of Cyp142 (clofazimine, all p < 10-7.7). We observed that artefactual signals of cross-resistance could be unravelled based on the relative effect size on MIC. Our study demonstrates the ability of very large-scale studies to substantially improve our knowledge of genetic variants associated with antimicrobial resistance in M. tuberculosis.

Minos: variant adjudication and joint genotyping of cohorts of bacterial genomes

There are many short-read variant-calling tools, with different strengths and weaknesses. We present a tool, Minos, which combines outputs from arbitrary variant callers, increasing recall without loss of precision. We benchmark on 62 samples from three bacterial species and an outbreak of 385 Mycobacterium tuberculosis samples. Minos also enables joint genotyping; we demonstrate on a large (N=13k) M. tuberculosis cohort, building a map of non-synonymous SNPs and indels in a region where all such variants are assumed to cause rifampicin resistance. We quantify the correlation with phenotypic resistance and then replicate in a second cohort (N=10k).

Intra-host genetic population diversity: Role in emergence and persistence of drug resistance among Mycobacterium tuberculosis complex minor variants

Drug resistant tuberculosis (DR-TB) is an important public health issue in different parts of the world. Mycobacterium tuberculosis complex variants (MTBC vars) preferentially infect certain hosts, limiting their distribution to different ecosystems. However, MTBC vars can infect other hosts beyond their preferred target potentially contributing to persistence of drug resistance (DR) in other niches. Here, we performed a comprehensive intra-host genetic analysis for the identification of DR-related mutations among all MTBC minor vars whole genome sequences (8,095 strains) publicly available worldwide. High confidence drug-resistance mutations in katG (isoniazid), rpsL (streptomycin), pncA (pyrazinamide), rpoB (rifampicin) and gyrA (fluoroquinolones) genes were identified among intrahost minor sub-populations in 197 different strains (2.43%) belonging to vars africanum, bovis, caprae, microti, orygis and pinnipedii. In addition, a three-dimensional structure modeling analysis to assess the role of novel mutations was also performed. Our findings highlight the importance of detecting discrete intra-host populations carrying DR mutations.

Mycobacterium tuberculosis Lineages Associated with Mutations and Drug Resistance in Isolates from India

Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. In this study, we performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand the genomic diversity. A retrospective collection of 498 M. tuberculosis isolates submitted to the National Institute for Research in Tuberculosis for phenotypic susceptibility testing between 2014 to 2016 were sequenced. Genotypic resistance prediction was performed using known resistance-conferring determinants. Genotypic and phenotypic results for 12 antituberculosis drugs were compared, and sequence data were explored to characterize lineages and their association with drug resistance. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction. IMPORTANCE Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. We performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand genomic diversity. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction.

Whole genome sequencing analysis to evaluate the influence of T2DM on polymorphisms associated with drug resistance in M. tuberculosis

Background

Type 2 diabetes mellitus (T2DM) has been associated with treatment failure, and the development of drug resistance in tuberculosis (TB). Also, whole-genome sequencing has provided a better understanding and allowed the growth of knowledge about polymorphisms in genes associated with drug resistance. Considering the above, this study analyzes genome sequences to evaluate the influence of type 2 diabetes mellitus in the development of mutations related to tuberculosis drug resistance. M. tuberculosis isolates from individuals with (n = 74), and without (n = 74) type 2 diabetes mellitus was recovered from online repositories, and further analyzed.

Results

The results showed the presence of 431 SNPs with similar proportions between diabetics, and non-diabetics individuals (48% vs. 52%), but with no significant relationship. A greater number of mutations associated with rifampicin resistance was observed in the T2DM-TB individuals (23.2% vs. 16%), and the exclusive presence of rpoBQ432L, rpoBQ432P, rpoBS441L, and rpoBH445L variants. While these variants are not private to T2DM-TB cases they are globally rare highlighting a potential role of T2DM. The phylogenetic analysis showed 12 sublineages, being 4.1.1.3, and 4.1.2.1 the most prevalent in T2DM-TB individuals but not differing from those most prevalent in their geographic location. Four clonal complexes were found, however, no significant relationship with T2DM was observed. Samples size and potential sampling biases prevented us to look for significant associations.

Conclusions

The occurrence of globally rare rifampicin variants identified only in isolates from individuals with T2DM could be due to the hyperglycemic environment within the host. Therefore, further studies about the dynamics of SNPs’ generation associated with antibiotic resistance in patients with diabetes mellitus are necessary.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08709-z.