Mycobacterium tuberculosis pyrazinamide resistance determinants: a multicenter study

Potential of nanopore sequencing for tuberculosis diagnosis and drug resistance detection

OBJECTIVES

This study evaluates the effectiveness of nanopore sequencing for accurate detection of Mycobacterium tuberculosis pathogens and drug resistance mutations in clinical specimens.

METHODS

A retrospective analysis of 2,421 specimens from suspected tuberculosis patients admitted to Xi'an Chest Hospital from 2022 to 2023 was conducted, with 131 specimens undergoing via real-time, fluorescence-based quantitative Polymerase Chain Reaction (qPCR), simultaneous amplification and testing RNA (RNA), Mycobacterium culture, Mycobacterium smear, and nanopore sequencing. Employing clinical tuberculosis diagnoses as the gold standard, sensitivity, specificity, positive predictive value, negative predictive value, concordance rate, and Kappa coefficient were measured for the five detection techniques. We compared nanopore sequencing with the Melting Curve method to detect drug-resistant gene mutations.

RESULTS

Nanopore sequencing has a significantly higher sensitivity (0.786) for tuberculosis diagnosis compared to qPCR (0.411), RNA (0.411), Mycobacterium culture (0.402), and Mycobacterium smear (0.241), against the gold-standard clinical diagnosis. It also exhibited a greater concordance rate (0.809) and Kappa coefficient (0.488), and outperformed the other methods in terms of the area under the ROC curve. Nanopore sequencing surpassed the Melting Curve method in identifying drug-resistant mutations.

CONCLUSION

Nanopore sequencing significantly enhances the detection of tuberculosis pathogens and drug-resistant genes.

Biphasic Medium Using Nicotinamide for Detection of Pyrazinamide Resistance in Mycobacterium tuberculosis

Reliable drug susceptibility testing of pyrazinamide (PZA) is technically difficult, since PZA activity is pH sensitive. The aim of this study was to evaluate a biphasic medium assay (BMA) for the reliable detection of PZA resistance in Mycobacterium tuberculosis (MTB) using nicotinamide (NIC) as a surrogate for PZA and identifying the appropriate cut-off value for the assay. The PZA susceptibility of 122 multidrug-resistant tuberculosis (MDR-TB) isolates and 39 drug-susceptible tuberculosis (DS-TB) isolates was examined using the BMA with NIC at four different concentrations (250, 500, 1000, and 2000 mg/L) and comparing the results with results from the BACTEC MGIT 960 reference method. Out of 122 MDR-TB isolates, 40 were identified as resistant by the BACTEC MGIT 960 system, of which 92.5% contained mutations within their pncA gene plus promoter region. A minimum inhibitory concentration of NIC ≥ 1000 mg/L was used as the cut-off concentration to define resistance in correlation with the MGIT 960 outcomes. NIC-BMA had a sensitivity of 90.91%, a specificity of 100%, and an accuracy of 97.52% compared with the MGIT 960 method. NIC-BMA is a promising assay to screen PZA resistance in microbiological laboratories without automation or advanced molecular instruments.

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.

Prediction of pyrazinamide resistance in Mycobacterium tuberculosis using structure-based machine-learning approaches

Background

Pyrazinamide is one of four first-line antibiotics used to treat tuberculosis; however, antibiotic susceptibility testing for pyrazinamide is challenging. Resistance to pyrazinamide is primarily driven by genetic variation in pncA, encoding an enzyme that converts pyrazinamide into its active form.

Methods

We curated a dataset of 664 non-redundant, missense amino acid mutations in PncA with associated high-confidence phenotypes from published studies and then trained three different machine-learning models to predict pyrazinamide resistance. All models had access to a range of protein structural-, chemical- and sequence-based features.

Results

The best model, a gradient-boosted decision tree, achieved a sensitivity of 80.2% and a specificity of 76.9% on the hold-out test dataset. The clinical performance of the models was then estimated by predicting the binary pyrazinamide resistance phenotype of 4027 samples harbouring 367 unique missense mutations in pncA derived from 24 231 clinical isolates.

Conclusions

This work demonstrates how machine learning can enhance the sensitivity/specificity of pyrazinamide resistance prediction in genetics-based clinical microbiology workflows, highlights novel mutations for future biochemical investigation, and is a proof of concept for using this approach in other drugs.

Prevalence, Transmission and Genetic Diversity of Pyrazinamide Resistance Among Multidrug-Resistant Mycobacterium tuberculosis Isolates in Hunan, China

Background

China is a country with a burden of high rates of both TB and multidrug-resistant TB (MDR-TB). However, published data on pyrazinamide (PZA) resistance are still limited in Hunan province, China. This study investigated the prevalence, transmission, and genetic diversity of PZA resistance among multidrug-resistant Mycobacterium tuberculosis isolates in Hunan province.

Methods

Drug susceptibility testing (DST) with the Bactec MGIT 960 PZA kit and pyrazinamidase (PZase) testing were conducted on all 298 MDR clinical isolates. Moreover, 24-locus MIRU-VNTR and DNA sequencing of pncA, rpsA, and panD genes were conducted on 180 PZA-resistant (PZA-R) isolates.

Results

The prevalence of PZA resistance among MDR-TB strains reached 60.4%. Newly diagnosed PZA-R TB patients and clustered isolates with identical pncA, rpsA, and panD mutations showed that transmission of PZA-R isolates played a significant role in the formation of PZA-R TB. Ninety-eight mutation patterns were observed in the pncA among 180 PZA-R isolates, and seventy-one (72.4%) were point mutations. Twenty-four of these mutations are new, including 2 base substitutions (V93G and T153S) and 22 nucleotide deletions or insertions. The W119C was found in PZA-S isolates, on the other hand, F94L and V155A mutations were found in both PZA resistant and susceptible isolates with positive PZase activity, indicating that they were not associated with PZA resistance. This is not entirely in line with the WHO catalogue. Ten novel rpsA mutations were found in 10 PZA-R isolates, which all combined with mutations in pncA. Thus, it is unpredictable whether these mutations in rpsA can impact PZA resistance. No panD mutation was found in all PZA-R isolates.

Conclusion

DNA sequencing of pncA and PZase activity testing have great potential in predicting PZA resistance.

Whole-genome sequencing-based analyses of drug-resistant Mycobacterium tuberculosis from Taiwan

Drug-resistant tuberculosis (DR-TB) posed challenges to global TB control. Whole-genome sequencing (WGS) is recommended for predicting drug resistance to guide DR-TB treatment and management. Nevertheless, data are lacking in Taiwan. Phenotypic drug susceptibility testing (DST) of 12 anti-TB drugs was performed for 200 Mycobacterium tuberculosis isolates. WGS was performed using the Illumina platform. Drug resistance profiles and lineages were predicted in silico using the Total Genotyping Solution for TB (TGS-TB). Using the phenotypic DST results as a reference, WGS-based prediction demonstrated high concordance rates of isoniazid (95.0%), rifampicin (RIF) (98.0%), pyrazinamide (98.5%) and fluoroquinolones (FQs) (99.5%) and 96.0% to 99.5% for second-line injectable drugs (SLIDs); whereas, lower concordance rates of ethambutol (87.5%), streptomycin (88.0%) and ethionamide (84.0%). Furthermore, minimum inhibitory concentrations confirmed that RIF rpoB S450L, FQs gyrA D94G and SLIDs rrs a1401g conferred high resistance levels. Besides, we identified lineage-associated mutations in lineage 1 (rpoB H445Y and fabG1 c-15t) and predominant lineage 2 (rpoB S450L and rpsL K43R). The WGS-based prediction of drug resistance is highly concordant with phenotypic DST results and can provide comprehensive genetic information to guide DR-TB precision therapies in Taiwan.

Mutations Associated with Pyrazinamide Resistance in Mycobacterium tuberculosis: A Review and Update

Pyrazinamide (PZA) has remained a keystone of tuberculosis (TB) therapy, and it possesses high imperative sterilizing action that can facilitate reduction in the present chemotherapy regimen. The combination of PZA works both with first- and second-line TB drugs, notably fluoroquinolones, clofazimine, bedaquiline, delamanid and pretomanid. Pyrazinamide inhibits various targets that are involved in different cellular processes like energy production (pncA), trans-translation (rpsA) and pantothenate/coenzyme A (panD) which are required for persistence of the pathogen. It is well known that pncA gene encoding pyrazinamidase is involved in the transition of PZA into the active form of pyrazinoic acid, which implies that mutation in the pncA gene can develop PZA resistance in Mycobacterium tuberculosis (M. tuberculosis) strain leading to a major clinical and public health concern. Therefore, it is very crucial to understand its resistance mechanism and to detect it precisely to help in the management of the disease. Scope of this review is to have a deep understanding of molecular mechanism of PZA resistance with its multiple targets which would help study the association of mutations and its resistance in M. tuberculosis. This will in turn help learn about the resistance of PZA and develop more accurate molecular diagnostic tool for drug-resistant TB in future TB therapy.

Pyrazinamide resistance in rifampicin discordant tuberculosis

Introduction

Mycobacterium tuberculosis strains with phenotypically susceptible rpoB mutations (rifampicin discordant) have emerged following implementation of rapid molecular drug resistance testing for tuberculosis. Whilst rifampicin resistance is known to be associated with resistance to other rifamycins (rifapentine and rifabutin) as well as isoniazid and pyrazinamide, rifampicin discordant strains have shown high rates of susceptibility to isoniazid and rifabutin. However, pyrazinamide susceptibly testing results have not been reported.

Materials and methods

We evaluated pyrazinamide resistance in 80 rifampicin discordant and 25 rifampicin and isoniazid susceptible isolates from KwaZulu-Natal in South Africa using Mycobacteria Growth Indicator Tube method and sequencing of the pncA. We also compared susceptibility of pyrazinamide with that of isoniazid.

Results

Pyrazinamide resistance was found in 6/80 (7.5%) rifampicin discordant isolates. All pyrazinamide resistant isolates were also resistant to isoniazid and pyrazinamide resistance was found to be associated with isoniazid resistance. No pyrazinamide resistance was found among the isoniazid susceptible isolates.

Conclusion

Given the low prevalence of pyrazinamide resistance in rifampicin discordant TB, this anti-TB drug still has a significant role in the treatment of these patients. Performing pyrazinamide susceptibility testing remains a challenge, our findings show that isoniazid susceptible isolates are unlikely to be resistant to pyrazinamide among the discordant TB isolates.

Detection of Mutations in pncA in Mycobacterium tuberculosis Clinical Isolates from Nepal in Association with Pyrazinamide Resistance

Without the proper information on pyrazinamide (PZA) susceptibility of Mycobacterium tuberculosis (MTB), PZA is inappropriately recommended for the treatment of both susceptible and multidrug-resistant tuberculosis (MDR-TB) in Nepal. This study aimed to collect information regarding PZA susceptibility in MTB isolates from Nepal by analyzing pncA and its upstream regulatory region (URR). A total of 211 MTB isolates were included in this study. Sequence analysis of pncA and its URR was performed to assess PZA resistance. First-line drug susceptibility testing, spoligotyping, and sequence analysis of rpoB, katG, the inhA regulatory region, gyrA, gyrB, and rrs were performed to assess their association with pncA mutation. Sequencing results reveal that 125 (59.2%) isolates harbored alterations in pncA and its URR. A total of 57 different mutation types (46 reported and 11 novel) were scattered throughout the whole length of the pncA gene. Eighty-seven isolates (41.2%) harbored mutations in pncA, causing PZA resistance in MTB. There was a more significant association of pncA alterations in MDR/pre-extensively drug-resistant (Pre-XDR) TB than in mono-resistant/pan-susceptible TB (p < 0.005). This first report on the increasing level of PZA resistance in DR-TB in Nepal highlights the importance of PZA susceptibility testing before DR-TB treatment.

Prediction of drug resistance by Sanger sequencing of Mycobacterium tuberculosis complex strains isolated from multidrug resistant tuberculosis suspect patients in Ethiopia

Background

Ethiopia is one of the high multidrug-resistant tuberculosis (MDR-TB) burden countries. However, phenotypic drug susceptibility testing can take several weeks due to the slow growth of Mycobacterium tuberculosis complex (MTBC) strains. In this study, we assessed the performance of a Sanger sequencing approach to predict resistance against five anti-tuberculosis drugs and the pattern of resistance mediating mutations.

Methods

We enrolled 226 MTBC culture-positive MDR-TB suspects and collected sputum specimens and socio-demographic and TB related data from each suspect between June 2015 and December 2016 in Addis Ababa, Ethiopia. Phenotypic drug susceptibility testing (pDST) for rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin using BACTEC MGIT 960 was compared with the results of a Sanger sequencing analysis of seven resistance determining regions in the genes rpoB, katG, fabG-inhA, pncA, embB, rpsL, and rrs.

Result

DNA isolation for Sanger sequencing was successfully extracted from 92.5% (209/226) of the MTBC positive cultures, and the remaining 7.5% (17/226) strains were excluded from the final analysis. Based on pDST results, drug resistance proportions were as follows: isoniazid: 109/209 (52.2%), streptomycin: 93/209 (44.5%), rifampicin: 88/209 (42.1%), ethambutol: 74/209 (35.4%), and pyrazinamide: 69/209 (33.0%). Resistance against isoniazid was mainly mediated by the mutation katG S315T (97/209, 46.4%) and resistance against rifampicin by rpoB S531L (58/209, 27.8%). The dominating resistance-conferring mutations for ethambutol, streptomycin, and pyrazinamide affected codon 306 in embB (48/209, 21.1%), codon 88 in rpsL (43/209, 20.6%), and codon 65 in pncA (19/209, 9.1%), respectively. We observed a high agreement between phenotypic and genotypic DST, such as 89.9% (at 95% confidence interval [CI], 84.2%–95.8%) for isoniazid, 95.5% (95% CI, 91.2%–99.8%) for rifampicin, 98.6% (95% CI, 95.9–100%) for ethambutol, 91.3% (95% CI, 84.6–98.1%) for pyrazinamide and 57.0% (95% CI, 46.9%–67.1%) for streptomycin.

Conclusion

We detected canonical mutations implicated in resistance to rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin. High agreement with phenotypic DST results for all drugs renders Sanger sequencing promising to be performed as a complementary measure to routine phenotypic DST in Ethiopia. Sanger sequencing directly from sputum may accelerate accurate clinical decision-making in the future.

Drug resistance prediction and resistance genes identification in Mycobacterium tuberculosis based on a hierarchical attentive neural network utilizing genome-wide variants

Prediction of antimicrobial resistance based on whole-genome sequencing data has attracted greater attention due to its rapidity and convenience. Numerous machine learning-based studies have used genetic variants to predict drug resistance in Mycobacterium tuberculosis (MTB), assuming that variants are homogeneous, and most of these studies, however, have ignored the essential correlation between variants and corresponding genes when encoding variants, and used a limited number of variants as prediction input. In this study, taking advantage of genome-wide variants for drug-resistance prediction and inspired by natural language processing, we summarize drug resistance prediction into document classification, in which variants are considered as words, mutated genes in an isolate as sentences, and an isolate as a document. We propose a novel hierarchical attentive neural network model (HANN) that helps discover drug resistance-related genes and variants and acquire more interpretable biological results. It captures the interaction among variants in a mutated gene as well as among mutated genes in an isolate. Our results show that for the four first-line drugs of isoniazid (INH), rifampicin (RIF), ethambutol (EMB) and pyrazinamide (PZA), the HANN achieves the optimal area under the ROC curve of 97.90, 99.05, 96.44 and 95.14% and the optimal sensitivity of 94.63, 96.31, 92.56 and 87.05%, respectively. In addition, without any domain knowledge, the model identifies drug resistance-related genes and variants consistent with those confirmed by previous studies, and more importantly, it discovers one more potential drug-resistance-related gene.

A Comprehensive Evaluation of GeneLEAD VIII DNA Platform Combined to Deeplex Myc-TB® Assay to Detect in 8 Days Drug Resistance to 13 Antituberculous Drugs and Transmission of Mycobacterium tuberculosis Complex Directly From Clinical Samples

The GeneLEAD VIII (Diagenode, Belgium) is a new, fully automated, sample-to-result precision instrument for the extraction of DNA and PCR detection of Mycobacterium tuberculosis complex (MTBC) directly from clinical samples. The Deeplex Myc-TB^® assay (Genoscreen, France) is a diagnostic kit based on the deep sequencing of a 24-plexed amplicon mix allowing simultaneously the detection of resistance to 13 antituberculous (antiTB) drugs and the determination of spoligotype. We evaluated the performance of a strategy combining the both mentioned tools to detect directly from clinical samples, in 8 days, MTBC and its resistance to 13 antiTB drugs, and identify potential transmission of strains from patient-to-patient. Using this approach, we screened 112 clinical samples (65 smear-negative) and 94 MTBC cultured strains. The sensitivity and the specificity of the GeneLEAD/Deeplex Myc-TB approach for MTBC detection were 79.3% and 100%, respectively. One hundred forty successful Deeplex Myc-TB results were obtained for 46 clinical samples and 94 strains, a total of 85.4% of which had a Deeplex Myc-TB susceptibility and resistance prediction consistent with phenotypic drug susceptibility testing (DST). Importantly, the Deeplex Myc-TB assay was able to detect 100% of the multidrug-resistant (MDR) MTBC tested. The lowest concordance rates were for pyrazinamide, ethambutol, streptomycin, and ethionamide (84.5%, 81.5%, 73%, and 55%, respectively) for which the determination of susceptibility or resistance is generally difficult with current tools. One of the main difficulties of Deeplex Myc-TB is to interpret the non-synonymous uncharacterized variants that can represent up to 30% of the detected single nucleotide variants. We observed a good level of concordance between Deeplex Myc-TB-spoligotyping and MIRU-VNTR despite a lower discriminatory power for spoligotyping. The median time to obtain complete results from clinical samples was 8 days (IQR 7–13) provided a high-throughput NGS sequencing platform was available. Our results highlight that the GeneLEAD/Deeplex Myc-TB approach could be a breakthrough in rapid diagnosis of MDR TB in routine practice.

Application of Computational Methods in Understanding Mutations in Mycobacterium tuberculosis Drug Resistance

The emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb) impedes the End TB Strategy by the World Health Organization aiming for zero deaths, disease, and suffering at the hands of tuberculosis (TB). Mutations within anti-TB drug targets play a major role in conferring drug resistance within Mtb; hence, computational methods and tools are being used to understand the mechanisms by which they facilitate drug resistance. In this article, computational techniques such as molecular docking and molecular dynamics are applied to explore point mutations and their roles in affecting binding affinities for anti-TB drugs, often times lowering the protein’s affinity for the drug. Advances and adoption of computational techniques, chemoinformatics, and bioinformatics in molecular biosciences and resources supporting machine learning techniques are in abundance, and this has seen a spike in its use to predict mutations in Mtb. This article highlights the importance of molecular modeling in deducing how point mutations in proteins confer resistance through destabilizing binding sites of drugs and effectively inhibiting the drug action.

Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages

Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.

Overcoming the Challenges of Pyrazinamide Susceptibility Testing in Clinical Mycobacterium tuberculosis Isolates

Pyrazinamide (PZA) is one of the first-line agents used for the treatment of tuberculosis. However, current phenotypic PZA susceptibility testing in the Bactec MGIT 960 system is unreliable, and false resistance is well documented. Rapid identification of resistance-associated mutations can confirm the phenotypic result. This study aimed to investigate the use of genotypic methods in combination with phenotypic susceptibility testing for confirmation of PZA-resistant Mycobacterium tuberculosis isolates. Sanger sequencing and/or whole-genome sequencing were performed to detect mutations in pncA, rpsA, panD, and clpC1. Isolates were screened for heteroresistance, and PZA susceptibility testing was performed using the Bactec MGIT 960 system using a reduced inoculum to investigate false resistance. Overall, 40 phenotypically PZA-resistant isolates were identified. Of these, PZA resistance was confirmed in 22/40 (55%) isolates by detecting mutations in the pncA, rpsA, and panD genes. Of the 40 isolates, 16 (40%) were found to be susceptible using the reduced inoculum method (i.e., false resistance). No mutations were detected in two PZA-resistant isolates. False resistance was observed in isolates with MICs close to the critical concentration. In particular, East African Indian strains (lineage 1) appeared to have an elevated MIC that is close to the critical concentration. While this study illustrates the complexity and challenges associated with PZA susceptibility testing of M. tuberculosis, we conclude that a combination of genotypic and phenotypic drug susceptibility testing methods is required for accurate detection of PZA resistance.

Phenotypic and molecular characterization of pyrazinamide resistance among multidrug-resistant Mycobacterium tuberculosis isolates in Ningbo, China

Background

Detection of pyrazinamide (PZA) resistance in Mycobacterium tuberculosis (TB) patients is critical, especially in dealing with multidrug-resistant Mycobacterium tuberculosis (MDR-TB) case. Up to date, PZA drug susceptibility testing (DST) has not been regularly performed in China. The prevalence and molecular characteristics of PZA resistance in M.tuberculosis isolates, especially MDR-TB have not been studied in Ningbo, China. This study aimed to analyze the phenotypic and molecular characterization of PZA resistance among MDR-TB isolates in Ningbo.

Methods

A total of 110 MDR-TB isolates were collected from the TB patients who were recorded at local TB dispensaries in Ningbo. All clinical isolates were examined by drug susceptibility testing and genotyping. DNA sequencing was used to detect mutations in the pncA gene associated with PZA resistance.

Results

The prevalence of PZA resistance among MDR-TB strains in Ningbo was 59.1%. With regard to the history and the outcome of treatments among MDR-TB cases, the percentages of re-treated MDR-TB patients in the PZA-resistant group and of successful patients in PZA-susceptible group were significantly higher than the ones in the PZA-susceptible group and in the PZA-resistant group, respectively (P = 0.027, P = 0.020). The results showed that the resistance of streptomycin (67.7% vs 46.7%, P = 0.027), ethambutol (56.9% vs 33.3%, P = 0.015), ofloxacin (43.1% vs 11.1%, P = 0.000), levofloxacin (43.1% vs 11.1%, P = 0.000), pre-XDR (pre-Xtensively Drug Resistance) (38.5% vs 15.6%, P = 0.009), were more frequently adverted among PZA-resistant isolates compared with PZA-susceptible isolates. In addition, 110 MDR-TB was composed of 87 (PZA resistant, 78.5%) Beijing strains and 23 (PZA resistant, 21.5%) non-Beijing strains. Fifty-four out of 65 (83.1%) PZA-resistant MDR strains harbored a mutation located in the pncA gene and the majority (90.7%) were point mutations. Compared with the phenotypic characterization, DNA sequencing of pncA has sensitivity and specificity of 83.1 and 95.6%.

Conclusion

The mutations within pncA gene was the primary mechanism of PZA resistance among MDR-TB and DNA sequencing of pncA gene could provide a rapid detection evidence in PZA drug resistance of MDR-TB in Ningbo.

Facile synthesis and antimycobacterial activity of isoniazid, pyrazinamide and ciprofloxacin derivatives

Several rationally designed isoniazid (INH), pyrazinamide (PZA) and ciprofloxacin (CPF) derivatives were conveniently synthesized and evaluated in vitro against H37Rv Mycobacterium tuberculosis (M. tb) strain. CPF derivative 16 displayed a modest activity (MIC = 16 µg/ml) and was docked into the M. tb DNA gyrase. Isoniazid-pyrazinoic acid (INH-POA) hybrid 21a showed the highest potency in our study (MIC = 2 µg/ml). It also retained its high activity against the other tested M. tb drug-sensitive strain (DS) V4207 (MIC = 4 µg/ml) and demonstrated negligible cytotoxicity against Vero cells (IC50  ≥ 64 µg/ml). Four tested drug-resistant (DR) M. tb strains were refractory to 21a, similar to INH, whilst being sensitive to CPF. Compound 21a was also inactive against two non-tuberculous mycobacterial (NTM) strains, suggesting its selective activity against M. tb. The noteworthy activity of 21a against DS strains and its low cytotoxicity highlight its potential to treat DS M. tb.

Quantifying transmission fitness costs of multi-drug resistant tuberculosis

As multi-drug resistant tuberculosis (MDR-TB) continues to spread, investigating the transmission potential of different drug-resistant strains becomes an ever more pressing topic in public health. While phylogenetic and transmission tree inferences provide valuable insight into possible transmission chains, phylodynamic inference combines evolutionary and epidemiological analyses to estimate the parameters of the underlying epidemiological processes, allowing us to describe the overall dynamics of disease spread in the population. In this study, we introduce an approach to Mycobacterium tuberculosis (M. tuberculosis) phylodynamic analysis employing an existing computationally efficient model to quantify the transmission fitness costs of drug resistance with respect to drug-sensitive strains. To determine the accuracy and precision of our approach, we first perform a simulation study, mimicking the simultaneous spread of drug-sensitive and drug-resistant tuberculosis (TB) strains. We analyse the simulated transmission trees using the phylodynamic multi-type birth-death model (MTBD, (Kühnert et al., 2016)) within the BEAST2 framework and show that this model can estimate the parameters of the epidemic well, despite the simplifying assumptions that MTBD makes compared to the complex TB transmission dynamics used for simulation. We then apply the MTBD model to an M. tuberculosis lineage 4 dataset that primarily consists of MDR sequences. Some of the MDR strains additionally exhibit resistance to pyrazinamide - an important first-line anti-tuberculosis drug. Our results support the previously proposed hypothesis that pyrazinamide resistance confers a transmission fitness cost to the bacterium, which we quantify for the given dataset. Importantly, our sensitivity analyses show that the estimates are robust to different prior distributions on the resistance acquisition rate, but are affected by the size of the dataset - i.e. we estimate a higher fitness cost when using fewer sequences for analysis. Overall, we propose that MTBD can be used to quantify the transmission fitness cost for a wide range of pathogens where the strains can be appropriately divided into two or more categories with distinct properties.

Insights Into Mutations Induced Conformational Changes and Rearrangement of Fe2+ Ion in pncA Gene of Mycobacterium tuberculosis to Decipher the Mechanism of Resistance to Pyrazinamide

Pyrazinamide (PZA) is the first-line drug commonly used in treating Mycobacterium tuberculosis (Mtb) infections and reduces treatment time by 33%. This prodrug is activated and converted to an active form, Pyrazinoic acid (POA), by Pyrazinamidase (PZase) enzyme. Mtb resistance to PZA is the outcome of mutations frequently reported in pncA, rpsA, and panD genes. Among the mentioned genes, pncA mutations contribute to 72-99% of the total resistance to PZA. Thus, considering the vital importance of this gene in PZA resistance, its frequent mutations (D49N, Y64S, W68G, and F94A) were investigated through in-depth computational techniques to put conclusions that might be useful for new scaffolds design or structure optimization to improve the efficacy of the available drugs. Mutants and wild type PZase were used in extensive and long-run molecular dynamics simulations in triplicate to disclose the resistance mechanism induced by the above-mentioned point mutations. Our analysis suggests that these mutations alter the internal dynamics of PZase and hinder the correct orientation of PZA to the enzyme. Consequently, the PZA has a low binding energy score with the mutants compared with the wild type PZase. These mutations were also reported to affect the binding of Fe2+ ion and its coordinated residues. Conformational dynamics also revealed that β-strand two is flipped, which is significant in Fe2+ binding. MM-GBSA analysis confirmed that these mutations significantly decreased the binding of PZA. In conclusion, these mutations cause conformation alterations and deformities that lead to PZA resistance.

Mechanisms and detection methods of Mycobacterium tuberculosis rifampicin resistance: The phenomenon of drug resistance is complex

Tuberculosis (TB) is an infectious disease that poses a serious threat to human health. Rifampin (RIF) is an important first-line anti-TB drug, and rifampin resistance (RIF-R) is a key factor in formulating treatment regimen and evaluating the prognosis of TB. Compared with other drugs resistance, the RIF-R mechanism of Mycobacterium tuberculosis (M. tuberculosis) is one of the clearest, which is mainly caused by RIF resistance-related mutations in the rpoB gene. This provides a convenient condition for developing rapid detection methods, and also an ideal object for studying the general drug resistance mechanisms of M. tuberculosis. This review focuses on the mechanisms that influence the RIF resistance of M. tuberculosis and related detection methods. Besides the mutations in rpoB, M. tuberculosis can decrease the amount of drugs entering the cells, enhance the drugs efflux, and be heterogeneous RIF susceptibility to resist drug pressure. Based on the results of current researches, many genes participate in influencing the susceptibility to RIF, which indicates the phenomenon of M. tuberculosis drug resistance is very complex.

Comparative Utility of Genetic Determinants of Drug Resistance and Phenotypic Drug Susceptibility Profiling in Predicting Clinical Outcomes in Patients With Multidrug-Resistant Mycobacterium tuberculosis

Setting: Programmatic management of drug-resistant tuberculosis in Ningbo, China.

Objective: To assess whether data-driven genetic determinants of drug resistance patterns could outperform phenotypic drug susceptibility testing in predicting clinical meaningful outcomes among patients with multidrug-resistant tuberculosis (MDR-TB).

Design: We conducted a prospective cohort study of 104 MDR-TB patients. All MDR-TB isolates underwent drug susceptibility testing and genotyping for mutations that could cause drug resistance. Study outcomes were time to sputum smear conversion and probability of treatment success, as well as time to culture conversion within 6 months. Data were analyzed using latent class analysis, Kaplan–Meier curves, and Cox regression models.

Results: We report that latent class analysis of data identified two latent classes that predicted sputum smear conversion with P = 0.001 and area under receiver-operating characteristic curve of 0.73. The predicted latent class memberships were associated with superior capability in predicting sputum culture conversion at 6 months and overall treatment success compared to phenotypic drug susceptibility profiling using boosted logistic regression models.

Conclusion: These results suggest that genetic determinants of drug resistance in combination with phenotypic drug-resistant tests could serve as useful biomarkers in predicting treatment prognosis in MDR-TB.

Genome-wide analysis in Escherichia coli unravels a high level of genetic homoplasy associated with cefotaxime resistance

Cefotaxime (CTX) is a third-generation cephalosporin (3GC) commonly used to treat infections caused by Escherichia coli. Two genetic mechanisms have been associated with 3GC resistance in E. coli. The first is the conjugative transfer of a plasmid harbouring antibiotic-resistance genes. The second is the introduction of mutations in the promoter region of the ampC β-lactamase gene that cause chromosome-encoded β-lactamase hyperproduction. A wide variety of promoter mutations related to AmpC hyperproduction have been described. However, their link to CTX resistance has not been reported. We recultured 172 cefoxitin-resistant E. coli isolates with known CTX minimum inhibitory concentrations and performed genome-wide analysis of homoplastic mutations associated with CTX resistance by comparing Illumina whole-genome sequencing data of all isolates to a PacBio sequenced reference chromosome. We mapped the mutations on the reference chromosome and determined their occurrence in the phylogeny, revealing extreme homoplasy at the -42 position of the ampC promoter. The 24 occurrences of a T at the -42 position rather than the wild-type C, resulted from 18 independent C>T mutations in five phylogroups. The -42 C>T mutation was only observed in E. coli lacking a plasmid-encoded ampC gene. The association of the -42 C>T mutation with CTX resistance was confirmed to be significant (false discovery rate <0.05). To conclude, genome-wide analysis of homoplasy in combination with CTX resistance identifies the -42 C>T mutation of the ampC promotor as significantly associated with CTX resistance and underlines the role of recurrent mutations in the spread of antibiotic resistance.

Deep amplicon sequencing for culture-free prediction of susceptibility or resistance to 13 anti-tuberculous drugs

Conventional molecular tests for detecting Mycobacterium tuberculosis complex (MTBC) drug resistance on clinical samples cover a limited set of mutations. Whole-genome sequencing (WGS) typically requires culture.

Here, we evaluated the Deeplex Myc-TB targeted deep-sequencing assay for prediction of resistance to 13 anti-tuberculous drugs/drug classes, directly applicable on sputum.

With MTBC DNA tests, the limit of detection was 100–1000 genome copies for fixed resistance mutations. Deeplex Myc-TB captured in silico 97.1–99.3% of resistance phenotypes correctly predicted by WGS from 3651 MTBC genomes. On 429 isolates, the assay predicted 92.2% of 2369 first- and second-line phenotypes, with a sensitivity of 95.3% and a specificity of 97.4%. 56 out of 69 (81.2%) residual discrepancies with phenotypic results involved pyrazinamide, ethambutol and ethionamide, and low-level rifampicin or isoniazid resistance mutations, all notoriously prone to phenotypic testing variability. Only two out of 91 (2.2%) resistance phenotypes undetected by Deeplex Myc-TB had known resistance-associated mutations by WGS analysis outside Deeplex Myc-TB targets. Phenotype predictions from Deeplex Myc-TB analysis directly on 109 sputa from a Djibouti survey matched those of MTBSeq/PhyResSE/Mykrobe, fed with WGS data from subsequent cultures, with a sensitivity of 93.5/98.5/93.1% and a specificity of 98.5/97.2/95.3%, respectively. Most residual discordances involved gene deletions/indels and 3–12% heteroresistant calls undetected by WGS analysis or natural pyrazinamide resistance of globally rare “Mycobacterium canettii” strains then unreported by Deeplex Myc-TB. On 1494 arduous sputa from a Democratic Republic of the Congo survey, 14 902 out of 19 422 (76.7%) possible susceptible or resistance phenotypes could be predicted culture-free.

Deeplex Myc-TB may enable fast, tailored tuberculosis treatment.

The novel Deeplex Myc-TB molecular assay shows a high degree of accuracy for extensive prediction of susceptibility and resistance to 13 anti-tuberculous drugs, directly achievable without culture, which may enable fast, tailored tuberculosis treatmenthttps://bit.ly/3bAvcAt

Characterization of pncA Mutations and Prediction of PZA Resistance in Mycobacterium tuberculosis Clinical Isolates From Chongqing, China

Pyrazinamide (PZA) is widely used to treat drug-sensitive or multidrug resistance tuberculosis. However, conventional PZA susceptibility tests of clinical isolates are rather difficult because of the requirement of acid pH. Since resistance to pyrazinamide is primary mediated by mutation of pncA, an alternative way of PZA susceptibility test is to analyze the pyrazinamidase activities of Mycobacterium tuberculosis clinical isolates. Therefore, a database containing the full spectrum of pncA mutations along with pyrazinamidase activities will be beneficial. To characterize mutations of pncA in M. tuberculosis from Chongqing, China, the pncA gene was sequenced and analyzed in 465 clinical isolates. A total of 124 types of mutations were identified in 424 drug-resistant isolates, while no mutation was identified in the 31 pan-susceptible isolates. Ninety-four of the 124 mutations had previously been reported, and 30 new mutations were identified. Based on reported literatures, 294 isolates could be predicted resistant to pyrazinamide. Furthermore, pyrazinamidase activities of the 30 new mutations were tested using the Escherichia coli pncA gene knockout strain. The results showed that 24 of these new mutations (28 isolates) led to loss of pyrazinamidase activity and six (8 isolates) of them did not. Taken together, 322 isolates with pncA mutations could be predicted to be PZA resistant among the 424 drug-resistant isolates tested. Analysis of pncA mutations and their effects on pyrazinamidase activity will not only enrich our knowledge of comprehensive pncA mutations related with PZA resistance but also facilitate rapid molecular diagnosis of pyrazinamide resistance in M. tuberculosis.

Targeted next generation sequencing directly from sputum for comprehensive genetic information on drug resistant Mycobacterium tuberculosis

BACKGROUND

Timely drug resistance detection is essential to global tuberculosis management. Unfortunately, rapid molecular tests assess resistance to only a few drugs, with culture required for comprehensive susceptibility test results.

METHODS

We evaluated targeted next generation sequencing (tNGS) for tuberculosis on 40 uncultured sputum samples. Resistance profiles from tNGS were compared with profiles from Xpert MTB/RIF, line probe assay (LPA), pyrosequencing (PSQ), and phenotypic testing. Concordance, sensitivity, specificity, and overall test agreement were compared across assays.

RESULTS

tNGS provided results for 39 of 40 samples (97.5%) with faster turnaround than phenotypic testing (median 3 vs. 21 days, p = 0.0068). Most samples were isoniazid and rifampin resistant (N = 31, 79.5%), 21 (53.8%) were fluoroquinolone resistant, and 3 (7.7%) were also resistant to Kanamycin. Half were of the Beijing lineage (N = 20, 51.3%). tNGS from uncultured sputum identified all resistance to isoniazid, rifampin, fluoroquinolones, and second-line injectable drugs that was identified by other methods. Agreement between tNGS and existing assays was excellent for isoniazid, rifampin, and SLDs, very good for levofloxacin, and good for moxifloxacin.

CONCLUSION

tNGS can rapidly identify tuberculosis, lineage, and drug resistance with faster turnaround than phenotypic testing. tNGS is a potential alternative to phenotypic testing in high-burden settings.

Immunological detection of pyrazine-2-carboxylic acid for the detection of pyrazinamide resistance in Mycobacterium tuberculosis

Pyrazinamide (PZA) susceptibility testing in Mycobacterium tuberculosis (Mtb) is a current area of development and PZA-resistant strains are increasingly prevalent. Previous studies have demonstrated that the detection of pyrazinoic acid (POA), the metabolite produced by the deamidation of PZA, is a good predictor for PZA resistance since a resistant strain would not convert PZA into POA at a critical required rate, whereas a susceptible strain will do, expelling POA to the extracellular environment at a certain rate, and allowing for quantification of this accumulated analyte. In order to quantify POA, an indirect competitive ELISA (icELISA) test using hyperimmune polyclonal rabbit serum against POA was developed: for this purpose, pure POA was first covalently linked to the highly immunogenic Keyhole Limpet Hemocyanine, and inoculated in rabbits. A construct made of bovine serum albumin (BSA) linked to pure POA and fixed at the bottom of wells was used as a competitor against spiked samples and liquid Mtb culture supernatants. When spiked samples (commercial POA alone) were analyzed, the half maximal inhibitory concentration (IC50) was 1.16 mg/mL, the limit of detection 200 μg/mL and the assay was specific (it did not detect PZA, IC50 > 20 mg/mL). However, culture supernatants (7H9-OADC-PANTA medium) disrupted the competition and a proper icELISA curve was not obtainable. We consider that, although we have shown that it is feasible to induce antibodies against POA, matrix effects could damage its analytical usefulness; multiple, upcoming ways to solve this obstacle are suggested.

Pyrazinamide resistance and pncA mutations in drug resistant Mycobacterium tuberculosis clinical isolates from Myanmar

Pyrazinamide (PZA) is an important anti-tuberculosis drug, which is active against semi-dormant bacilli and used as a component of first-line drugs and drug-resistant tuberculosis regimens. Mutations in pncA and its promoter region are main cause of PZA resistance. There are limited PZA susceptibility data as there is no routine drug susceptibility testing (DST) for PZA. This study was aimed to determine the proportion of PZA resistance among rifampicin-resistant tuberculosis patients and to identify mutations which are responsible for PZA resistance in pncA and its promoter region. Liquid-based DST was performed to detect PZA susceptibility on 192 culture positive rifampicin-resistant isolates collected from National Tuberculosis Reference Laboratory. Sequencing on pncA including its promoter region was performed and analysis was done on 157 isolates. Phenotypic PZA resistance was detected in 58.9% of isolates. Sixty-five different mutations were distributed in pncA or promoter region of 82 isolates. Sensitivity and specificity of pncA sequencing in detection of PZA resistance showed 89.8% and 95.6% respectively. High proportion of PZA resistance among rifampicin-resistant cases highlighted the need for effective treatment regimen development for PZA-resistant MDR-TB. It is also suggested that routine PZA susceptibility test should be incorporated to treatment monitoring regimen and National Drug Resistance surveys.

Characterization of Genomic Variants Associated with Resistance to Bedaquiline and Delamanid in Naive Mycobacterium tuberculosis Clinical Strains

The role of mutations in genes associated with phenotypic resistance to bedaquiline (BDQ) and delamanid (DLM) in Mycobacterium tuberculosis complex (MTBc) strains is poorly characterized. A clear understanding of the genetic variants' role is crucial to guide the development of molecular-based drug susceptibility testing (DST). In this work, we analyzed all mutations in candidate genomic regions associated with BDQ- and DLM-resistant phenotypes using a whole-genome sequencing (WGS) data set from a collection of 4,795 MTBc clinical isolates from six countries with a high burden of tuberculosis (TB). From WGS analysis, we identified 61 and 163 unique mutations in genomic regions potentially involved in BDQ- and DLM-resistant phenotypes, respectively. Importantly, all strains were isolated from patients who likely have never been exposed to these medicines. To characterize the role of mutations, we calculated the free energy variation upon mutations in the available protein structures of Ddn (DLM), Fgd1 (DLM), and Rv0678 (BDQ) and performed MIC assays on a subset of MTBc strains carrying mutations to assess their phenotypic effect. The combination of structural and phenotypic data allowed for cataloguing the mutations clearly associated with resistance to BDQ (n = 4) and DLM (n = 35), only two of which were previously described, as well as about a hundred genetic variants without any correlation with resistance. Significantly, these results show that both BDQ and DLM resistance-related mutations are diverse and distributed across the entire region of each gene target, which is of critical importance for the development of comprehensive molecular diagnostic tools.

Genotypic resistance of pyrazinamide but not MIC is associated with longer time to sputum culture conversion in patients with multidrug-resistant tuberculosis

BACKGROUND

PZA resistance in multidrug-resistant tuberculosis (MDR-TB) is common and it is not clear how it affects interim and treatment outcomes. Although rarely performed, phenotypic drug susceptibility testing (pDST) is used to define PZA resistance but genotypic DST (gDST) and minimum inhibitory concentration (MIC) could be beneficial. We aimed to assess the impact of PZA gDST and MIC on time to sputum culture conversion (SCC) and treatment outcome in patients with MDR-TB.

METHODS

Clinical, microbiological and treatment data was collected in this cohort study for all patients diagnosed with MDR-TB in Sweden 1992-2014. MIC, pDST and whole genome sequencing of the pncA, rpsA and panD genes were used to define PZA resistance. A Cox regression model was used for statistical analyses.

RESULTS

Of 157 patients with MDR-TB, 56.1% (n=88) had PZA resistant strains and 49.7% (n=78) were treated with PZA. In crude and adjusted analyses, PZA gDST resistance was associated with a 29-day longer time to SCC (hazard ratio [HR] 0.57, 95% confidence interval [CI] 0.36-0.89, p=0.013 and HR 0.49, 95% CI 0.29-0.82, p=0.007, respectively). A two-fold decrease in dilutions of PZA MIC for PZA susceptible strains showed no association with SCC in crude or adjusted analyses (HR 0.98, 95% CI 0.73-1.31, p=0.89). Genotypic DST and MIC for PZA were not associated with treatment outcome.

CONCLUSION

In patients with MDR-TB, gDST PZA resistance was associated with a longer time to SCC. Rapid PZA gDST is important to identify patients who may benefit from PZA treatment.

How To Optimally Combine Genotypic and Phenotypic Drug Susceptibility Testing Methods for Pyrazinamide

False-susceptible phenotypic drug-susceptibility testing (DST) results for pyrazinamide due to mutations with MICs close to the critical concentration (CC) confound the classification of pncA resistance mutations, leading to an underestimate of the specificity of genotypic DST. This could be minimized by basing treatment decisions on well-understood mutations and by adopting an area of technical uncertainty for phenotypic DST rather than only testing the CC, as is current practice for the Mycobacterium tuberculosis complex.

Determining the unbinding events and conserved motions associated with the pyrazinamide release due to resistance mutations of Mycobacterium tuberculosis pyrazinamidase

Pyrazinamide (PZA) is the only first-line antitubercular drug active against latent Mycobacterium tuberculosis (Mtb). It is activated to pyrazinoic acid by the pncA-encoded pyrazinamidase enzyme (PZase). Despite the emergence of PZA drug resistance, the underlying mechanisms of resistance remain unclear. This study investigated part of these mechanisms by modelling a PZA-bound wild type and 82 mutant PZase structures before applying molecular dynamics (MD) with an accurate Fe2+ cofactor coordination geometry. After observing nanosecond-scale PZA unbinding from several PZase mutants, an algorithm was developed to systematically detect ligand release via centre of mass distances (COM) and ligand average speed calculations, before applying the statistically guided network analysis (SGNA) method to investigate conserved protein motions associated with ligand unbinding. Ligand and cofactor perspectives were also investigated. A conserved pair of lid-destabilising motions was found. These consisted of (1) antiparallel lid and side flap motions; (2) the contractions of a flanking region within the same flap and residue 74 towards the core. Mutations affecting the hinge residues (H51 and H71), nearby residues or L19 were found to destabilise the lid. Additionally, other metal binding site (MBS) mutations delocalised the Fe2+ cofactor, also facilitating lid opening. In the early stages of unbinding, a wider variety of PZA poses were observed, suggesting multiple exit pathways. These findings provide insights into the late events preceding PZA unbinding, which we found to occur in some resistant PZase mutants. Further, the algorithm developed here to identify unbinding events coupled with SGNA can be applicable to other similar problems.

Targeted-Sequencing Workflows for Comprehensive Drug Resistance Profiling of Mycobacterium tuberculosis Cultures Using Two Commercial Sequencing Platforms: Comparison of Analytical and Diagnostic Performance, Turnaround Time, and Cost

Background

The emergence of Mycobacterium tuberculosis with complex drug resistance profiles necessitates a rapid and comprehensive drug susceptibility test for guidance of patient treatment. We developed two targeted-sequencing workflows based on Illumina MiSeq and Nanopore MinION for the prediction of drug resistance in M. tuberculosis toward 12 antibiotics.

Methods

A total of 163 M. tuberculosis isolates collected from Hong Kong and Ethiopia were subjected to a multiplex PCR for simultaneous amplification of 19 drug resistance-associated genetic regions. The amplicons were then barcoded and sequenced in parallel on MiSeq and MinION in respective batch sizes of 24 and 12 samples. A web-based bioinformatics pipeline, BacterioChek-TB, was developed to translate the raw datasets into clinician-friendly reports.

Results

Both platforms successfully sequenced all samples with mean read depths of 1,127× and 1,649×, respectively. The variant calling by MiSeq and MinION could achieve 100% agreement if variants with an allele frequency of &lt;40% reported by MinION were excluded. Both workflows achieved a mean clinical sensitivity of 94.8% and clinical specificity of 98.0% when compared with phenotypic drug susceptibility test (pDST). Turnaround times for the MiSeq and MinION workflows were 38 and 15 h, facilitating the delivery of treatment guidance at least 17–18 days earlier than pDST, respectively. The higher cost per sample on the MinION platform ($71.56) versus the MiSeq platform ($67.83) was attributed to differences in batching capabilities.

Conclusion

Our study demonstrates the interchangeability of MiSeq and MinION platforms for generation of accurate and actionable results for the treatment of tuberculosis.

Direct Determination of Pyrazinamide (PZA) Susceptibility by Sputum Microscopic Observation Drug Susceptibility (MODS) Culture at Neutral pH: the MODS-PZA Assay

Pyrazinamide (PZA) is considered the pivot drug in all tuberculosis treatment regimens due to its particular action on the persistent forms of Mycobacterium tuberculosis. However, no drug susceptibility test (DST) is considered sufficiently reliable for routine application. Although molecular tests are endorsed, their application is limited to known PZA resistance associated mutations. Microbiological DSTs for PZA have been restricted by technical limitations, especially the necessity for an acidic pH.

Pyrazinamide (PZA) is considered the pivot drug in all tuberculosis treatment regimens due to its particular action on the persistent forms of Mycobacterium tuberculosis. However, no drug susceptibility test (DST) is considered sufficiently reliable for routine application. Although molecular tests are endorsed, their application is limited to known PZA resistance associated mutations. Microbiological DSTs for PZA have been restricted by technical limitations, especially the necessity for an acidic pH. Here, for the first time, MODS culture at neutral pH was evaluated using high PZA concentrations (400 and 800 μg/ml) to determine PZA susceptibility directly from sputum samples. Sputum samples were cultured with PZA for up to 21 days at 37°C. Plate reading was performed at two time points: R1 (mean, 10 days) and R2 (mean, 13 days) for each PZA concentration. A consensus reference test, composed of MGIT-PZA, pncA sequencing, and the classic Wayne test, was used. A total of 182 samples were evaluated. The sensitivity and specificity for 400 μg/ml ranged from 76.9 to 89.7 and from 93.0 to 97.9%, respectively, and for 800 μg/ml ranged from 71.8 to 82.1 and from 95.8 to 98.6%, respectively. Compared to MGIT-PZA, our test showed a similar turnaround time (medians of 10 and 12 days for PZA-sensitive and -resistant isolates, respectively). In conclusion, MODS-PZA is presented as a fast, simple, and low-cost DST that could complement the MODS assay to evaluate resistance to the principal first-line antituberculosis drugs. Further optimization of test conditions would be useful in order to increase its performance.

The Bewildering Antitubercular Action of Pyrazinamide

Pyrazinamide (PZA) is a cornerstone antimicrobial drug used exclusively for the treatment of tuberculosis (TB). Due to its ability to shorten drug therapy by 3 months and reduce disease relapse rates, PZA is considered an irreplaceable component of standard first-line short-course therapy for drug-susceptible TB and second-line treatment regimens for multidrug-resistant TB. Despite over 60 years of research on PZA and its crucial role in current and future TB treatment regimens, the mode of action of this unique drug remains unclear. Defining the mode of action for PZA will open new avenues for rational design of novel therapeutic approaches for the treatment of TB. In this review, we discuss the four prevailing models for PZA action, recent developments in modulation of PZA susceptibility and resistance, and outlooks for future research and drug development.

Structure guided prediction of Pyrazinamide resistance mutations in pncA

Pyrazinamide plays an important role in tuberculosis treatment; however, its use is complicated by side-effects and challenges with reliable drug susceptibility testing. Resistance to pyrazinamide is largely driven by mutations in pyrazinamidase (pncA), responsible for drug activation, but genetic heterogeneity has hindered development of a molecular diagnostic test. We proposed to use information on how variants were likely to affect the 3D structure of pncA to identify variants likely to lead to pyrazinamide resistance. We curated 610 pncA mutations with high confidence experimental and clinical information on pyrazinamide susceptibility. The molecular consequences of each mutation on protein stability, conformation, and interactions were computationally assessed using our comprehensive suite of graph-based signature methods, mCSM. The molecular consequences of the variants were used to train a classifier with an accuracy of 80%. Our model was tested against internationally curated clinical datasets, achieving up to 85% accuracy. Screening of 600 Victorian clinical isolates identified a set of previously unreported variants, which our model had a 71% agreement with drug susceptibility testing. Here, we have shown the 3D structure of pncA can be used to accurately identify pyrazinamide resistance mutations. SUSPECT-PZA is freely available at: http://biosig.unimelb.edu.au/suspect_pza/.

Value of pyrazinamide for composition of new treatment regimens for multidrug-resistant Mycobacterium tuberculosis in China

Background

Pyrazinamide still may be a useful drug for treatment of rifampin-resistant (RR-TB) or multidrug-resistant tuberculosis (MDR-TB) in China while awaiting scale up of new drugs and regimens including bedaquiline and linezolid. The level of pyrazinamide resistance among MDR-TB patients in China is not well established. Therefore, we assessed pyrazinamide resistance in a representative sample and explored determinants and patterns of pncA mutations.

Methods

MDR-TB isolates from the 2007 national drug resistance survey of China were sub-cultured and examined for pyrazinamide susceptibility by BACTEC MGIT 960 method. pncA mutations were identified by sequencing. Characteristics associated with pyrazinamide resistance were analyzed using univariable and multivariable log-binominal regression.

Results

Of 401 MDR-TB isolates, 324 were successfully sub-cultured and underwent drug susceptibility testing. Pyrazinamide resistance was prevalent in 40.7% of samples, similarly among new and previously treated MDR-TB patients. Pyrazinamide resistance in MDR-TB patients was associated with lower age (adjusted OR 0.54; 95% CI, 0.34–0.87 for those aged ≧60 years compared to < 40 years). Pyrazinamide resistance was not associated with gender, residential area, previous treatment history and Beijing genotype. Of 132 patients with pyrazinamide resistant MDR-TB, 97 (73.5%) had a mutation in the pncA gene; with 61 different point mutations causing amino acid change, and 11 frameshifts in the pncA gene. The mutations were scattered throughout the whole pncA gene and no hot spot region was identified.

Conclusions

Pyrazinamide resistance among MDR-TB patients in China is common, although less so in elderly patients. Therefore, pyrazinamide should only be used for treatment of RR/MDR-TB in China if susceptibility is confirmed. Molecular testing for detection of pyrazinamide resistance only based on pncA mutations has certain value for the rapid detection of pyrazinamide resistance in MDR-TB strains but other gene mutations conferring to pyrazinamide resistance still need to be explored to increase its predictive ability .

Molecular characterization of multidrug-resistant Mycobacterium tuberculosis (MDR-TB) isolates identifies local transmission of infection in Kuwait, a country with a low incidence of TB and MDR-TB

Background

Increasing incidence of multidrug-resistant Mycobacterium tuberculosis infections is hampering global tuberculosis control efforts. Kuwait is a low-tuberculosis-incidence country, and ~ 1% of M. tuberculosis strains are resistant to rifampicin and isoniazid (MDR-TB). This study detected mutations in seven genes predicting resistance to rifampicin, isoniazid, pyrazinamide, ethambutol and streptomycin in MDR-TB strains. Sequence data were combined with spoligotypes for detecting local transmission of MDR-TB in Kuwait.

Methods

Ninety-three MDR-TB strains isolated from 12 Kuwaiti and 81 expatriate patients and 50 pansusceptible strains were used. Phenotypic drug susceptibility was determined by MGIT 460 TB/960 system. Mutations conferring resistance to rifampicin, isoniazid, pyrazinamide, ethambutol and streptomycin were detected by genotype MTBDRplus assay and/or PCR sequencing of three rpoB regions, katG codon 315 (katG315) + inhA regulatory region, pncA, three embB regions and rpsL + rrs-500–900 regions. Spoligotyping kit was used, spoligotypes were identified by SITVIT2, and phylogenetic tree was constructed by using MIRU-VNTRplus software. Phylogenetic tree was also constructed from concatenated sequences by MEGA7 software. Additional PCR sequencing of gidB and rpsA was performed for cluster isolates.

Results

Pansusceptible isolates contained wild-type sequences. Mutations in rpoB and katG and/or inhA were detected in 93/93 and 92/93 MDR-TB strains, respectively. Mutations were also detected for pyrazinamide resistance, ethambutol resistance and streptomycin resistance in MDR-TB isolates in pncA, embB and rpsL + rrs, respectively. Spoligotyping identified 35 patterns with 18 isolates exhibiting unique patterns while 75 isolates grouped in 17 patterns. Beijing genotype was most common (32/93), and 11 isolates showed nine orphan patterns. Phylogenetic analysis of concatenated sequences showed unique patterns for 51 isolates while 42 isolates grouped in 16 clusters. Interestingly, 22 isolates in eight clusters by both methods were isolated from TB patients typically within a span of 2 years. Five of eight clusters were confirmed by additional gidB and rpsA sequence data.

Conclusions

Our study provides the first insight into molecular epidemiology of MDR-TB in Kuwait and identified several potential clusters of local transmission of MDR-TB involving 2–6 subjects which had escaped detection by routine surveillance studies. Prospective detection of resistance-conferring mutations can identify possible cases of local transmission of MDR-TB in low MDR-TB settings.

Pharmacological and Molecular Mechanisms Behind the Sterilizing Activity of Pyrazinamide

Inclusion of pyrazinamide (PZA) in the tuberculosis (TB) drug regimen during the 1970s enabled a reduction in treatment duration from 12 to 6 months. PZA has this remarkable effect in patients despite displaying poor potency against Mycobacterium tuberculosis (Mtb) in vitro. The pharmacological basis for the in vivo sterilizing activity of the drug has remained obscure and its bacterial target controversial. Recently it was shown that PZA penetrates necrotic caseous TB lung lesions and kills nongrowing, drug-tolerant bacilli. Furthermore, it was uncovered that PZA inhibits bacterial Coenzyme A biosynthesis. It may block this pathway by triggering degradation of its target, aspartate decarboxylase. The elucidation of the pharmacological and molecular mechanisms of PZA provides the basis for the rational discovery of the next-generation PZA with improved in vitro potency while maintaining attractive pharmacological properties.

Tuberculosis, HIV, and viral hepatitis diagnostics in eastern Europe and central Asia: high time for integrated and people-centred services

Globally, high rates (and in the WHO European region an increasing prevalence) of co-infection with tuberculosis and HIV and HIV and hepatitis C virus exist. In eastern European and central Asian countries, the tuberculosis, HIV, and viral hepatitis programmes, including diagnostic services, are separate vertical structures. In this Personal View, we consider underlying reasons for the poor integration for these diseases, particularly in the WHO European region, and how to address this with an initial focus on diagnostic services. In part, this low integration has reflected different diagnostic development histories, global funding sources, and sample types used for diagnosis (eg, typically sputum for tuberculosis and blood for HIV and hepatitis C). Cooperation between services improved as patients with tuberculosis needed routine testing for HIV and vice versa, but financial, infection control, and logistical barriers remain. Multidisease diagnostic platforms exist, but to be used optimally, appropriate staff training and sensible understanding of different laboratory and infection control risks needs rapid implementation. Technically these ideas are all feasible. Poor coordination between these vertical systems remains unhelpful. There is a need to increase political and operational integration of diagnostic and treatment services and bring them closer to patients.

Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of Mycobacterium tuberculosis carrying katG-S315T in Hanoi, Vietnam

Drug-resistant tuberculosis (TB) is a serious global problem, and pathogen factors involved in the transmission of isoniazid (INH)-resistant TB have not been fully investigated. We performed whole genome sequencing of 332 clinical Mycobacterium tuberculosis (Mtb) isolates collected from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam. Using a bacterial genome-wide approach based on linear mixed models, we investigated the associations between 31-bp k-mers and clustered strains harboring katG-S315T, a major INH-resistance mutation in the present cohort and in the second panel previously published in South Africa. Five statistically significant genes, namely, PPE18/19, gid, emrB, Rv1588c, and pncA, were shared by the two panels. We further identified variants of the genes responsible for these k-mers, which are relevant to the spread of INH-resistant strains. Phylogenetic convergence test showed that variants relevant to PPE46/47-like chimeric genes were significantly associated with the same phenotype in Hanoi. The associations were further confirmed after adjustment for the confounders. These findings suggest that genomic variations of the pathogen facilitate the expansion of INH-resistance TB, at least in part, and our study provides a new insight into the mechanisms by which drug-resistant Mtb maintains fitness and spreads in Asia and Africa.

'Happy the man, who, studying nature's laws, Thro' known effects can trace the secret cause.' Do we have enough pieces to solve the pyrazinamide puzzle?

A low pH was assumed to be required for the activity of pyrazinoic acid (the active form of pyrazinamide) against Mycobacterium tuberculosis, but recently activity has been demonstrated at neutral pH. Renewed interest in pyrazinamide has led to an increasing number of potential targets and the suspicion that pyrazinamide is a 'dirty drug'. However, it is our opinion that the recent demonstration that pyrazinoic acid is active against PanD provides an alternative explanation for the secret of pyrazinamide's unusual activity. In this article we propose that PanD is the primary target of pyrazinoic acid but expression of pyrazinoic acid susceptibility requires an intact stress response. As the mycobacterial stress response requires the interaction of a number of genes, disruption of any could result in an inability to enter the susceptible phenotype. We believe this model can explain most of the recent observations of the seemingly diverse spectrum of activity of pyrazinamide.

Resistance To First-Line Antituberculosis Drugs And Prevalence Of pncA Mutations In Clinical Isolates Of Mycobacterium tuberculosis From Zunyi, Guizhou Province Of China

Background

China is one of the high-burden countries for multidrug-resistant tuberculosis (MDR-TB), and pyrazinamide is one of the anti-TB drugs used for the shorter MDR-TB treatment regimen. The aim of this study was to determine the correlation between pncA gene mutations and resistance to four first-line anti-TB drugs as well as treatment history in clinical isolates of Mycobacterium tuberculosis.

Patients and methods

M. tuberculosis clinical isolates were collected from 318 in-patients with smear-positive TB between October 2008 and September 2016 at a major hospital in Zunyi, Guizhou Province of China, and used for drug susceptibility testing against four first-line anti-TB drugs. Genomic DNA extracted from clinical isolates was used for PCR amplification and DNA sequencing of the pncA gene.

Results

Among 318 clinical isolates, 129 (40.6%), 170 (53.5%), 66 (20.8%) and 109 (34.3%) were resistant to rifampicin, isoniazid, ethambutol and streptomycin respectively. In addition, 124 clinical isolates were MDR-TB and 71.8% of them were previously treated cases. Sequencing results showed that 46.8% of MDR-TB and 2.2% of drug susceptible isolates harbored a pncA mutation, and 52 types of pncA mutations were detected from 64 isolates. The prevalence of pncA mutations in isolates resistant to first-line anti-TB drugs and previously treated TB cases was significantly higher than that in drug-susceptible isolates and new cases of TB.

Conclusion

High prevalence of pncA mutations in clinical isolates of M. tuberculosis from Zunyi, Guizhou Province of China, is correlated with resistance to four first-line anti-TB drugs, MDR-TB and previously treated TB cases.

Of Testing and Treatment: Implications of Implementing New Regimens for Multidrug-Resistant Tuberculosis

A novel, shorter-course regimen for treating multidrug-resistant (MDR) tuberculosis was recently recommended by the World Health Organization. However, the most appropriate use of drug susceptibility testing (DST) to support this regimen is less clear. Implementing countries must therefore often choose between using a standardized regimen despite high levels of underlying drug resistance or require more stringent DST prior to treatment initiation. The former carries a high likelihood of exposing patients to de facto monotherapy with a critical drug class (fluoroquinolones), whereas the latter could exclude large groups of patients from their most effective treatment option. We discuss the implications of this dilemma and argue for an approach that will integrate DST into the delivery of any novel antimicrobial regimen, without excessively stringent requirements. Such guidance could make the novel MDR tuberculosis regimen available to most patients while reducing the risk of generating additional drug resistance.

Pyrazinamide resistance and mutations L19R, R140H, and E144K in Pyrazinamidase of Mycobacterium tuberculosis

Pyrazinamide (PZA) is an important component of first-line antituberculosis drugs activated by Mycobacterium tuberculosis pyrazinamidase (PZase) into its active form pyrazinoic acid. Mutations in the pncA gene have been recognized as the major cause of PZA resistance. We detected some novel mutations, Leucine19Arginine (L19R), Arginine140Histidine (R140H), and Glutamic acid144 Lysine (E144K), in the pncA gene of PZA-resistant isolates in our wet lab PZA drug susceptibility testing and sequencing. As the molecular mechanism of resistance of these variants has not been reported earlier, we have performed multiple analyses to unveil different mechanisms of resistance because of PZase mutations L19R, R140H, and E144K. The mutants and native PZase structures were subjected to comprehensive computational molecular dynamics (MD) simulations at 100 nanoseconds in apo and drug-bound form. Mutants and native PZase binding pocket were compared to observe the consequence of mutations on the binding pocket size. Hydrogen bonding, Gibbs free energy, and natural ligand Fe ⁺² effect were also analyzed between native and mutants. A significant variation between native and mutant PZase structure activity was observed. The native PZase protein docking score was found to be the maximum, showing strong binding affinity in comparison with mutants. MD simulations explored the effect of the variants on the biological function of PZase. Hydrogen bonding, metal ion Fe ⁺² deviation, and fluctuation also seemed to be affected because of the mutations L19R, R140H, and E144K. The variants L19R, R140H, and E144K play a significant role in PZA resistance, altering the overall activity of native PZase, including metal ion Fe ⁺² displacement and free energy. This study offers valuable evidence for better management of drug-resistant tuberculosis.

Introducing molecular testing of pyrazinamide susceptibility improves multidrug-resistant tuberculosis treatment outcomes: a prospective cohort study

The current treatment for multidrug-resistant tuberculosis (MDR-TB) takes a lengthy period of 18-24 months and has a poor cure rate of 50-60%. A multicenter, prospective cohort study was conducted to assess the role of testing for molecular susceptibility to pyrazinamide (PZA) in optimising treatment for MDR-TB.We assigned 76 patients to an optimised molecular susceptibility group and 159 patients to a regular treatment group where PZA susceptibility was not determined. Of these patients, 152 were matched after propensity score matching (76 in the optimised group and 76 in the regular group). Treatment success rate was measured in the propensity-matched cohort as the primary outcome.Patients in the optimised group achieved a higher treatment success rate than those in the regular group (76.3% versus 55.3%, p=0.006). Of 51 patients with isolates that were susceptible to PZA and who were receiving a 12-month regimen, 42 (82.4%) were treated successfully. The optimised group showed faster culture conversion than the regular group (p=0.024). After exclusion of pre-extensively drug-resistant TB (pre-XDR-TB), the treatment outcome in the optimised group was still better than the regular group (83.1% versus 62.1%, p=0.009).Introducing molecular susceptibility testing for PZA improved the treatment outcomes for MDR-TB without the use of new drugs. Introducing PZA for patients with PZA-susceptible (PZA-S) MDR-TB allows the current regimen to be shortened to 12 months with comparable success rates to the World Health Organization (WHO) recommended shorter regimen.

Overcoming the pitfalls of automatic interpretation of whole genome sequencing data by online tools for the prediction of pyrazinamide resistance in Mycobacterium tuberculosis

OBJECTIVES

Automated online software tools that analyse whole genome sequencing (WGS) data without the need for bioinformatics expertise can motivate the implementation of WGS-based molecular drug susceptibility testing (DST) in routine diagnostic settings for tuberculosis (TB). Pyrazinamide (PZA) is a key drug for current and future TB treatment regimens; however, it was reported that predictive power for PZA resistance by the available tools is low. Therefore, this low predictive power may make users hesitant to use the tools. This study aimed to elucidate why and to uncover the real performance of the tools when taking into account their variation calling lists (manual inspection), not just their automated reporting system (default setting) that was evaluated by previous studies.

METHODS

WGS data from 191 datasets comprising 108 PZA-resistant and 83 susceptible strains were used to evaluate the potential performance of the available online tools (TB Profiler, TGS-TB, PhyResSE, and CASTB) for predicting phenotypic PZA resistance.

RESULTS

When taking into consideration the variation calling lists, 73 variants in total (47 non-synonymous mutations and 26 indels) in pncA were detected by TGS-TB and PhyResSE, covering all mutations for the 108 PZA-resistant strains. The 73 variants were confirmed by Sanger sequencing. TB Profiler also detected all but three complete loss, two large deletion at the 3'-end, and one relatively large insertion of pncA. On the other hand, many of the 73 variants were lacking in the automated reporting systems except by TGS-TB; of these variants, CASTB detected only 20. By applying the 'non-wild type sequence' approach for predicting PZA resistance, accuracy of the results significantly improved compared with that of the automated results obtained by each tool.

CONCLUSION

Users can obtain more accurate predictions for PZA resistance than previously reported by manually checking the results and applying the 'non-wild type sequence' approach.

Underestimated pyrazinamide resistance may compromise outcomes of pyrazinamide containing regimens for treatment of drug susceptible and multi-drug-resistant tuberculosis in Tanzania

BACKGROUND

Tuberculosis (TB) is the leading cause of death from an infectious disease and the roll-out of rapid molecular diagnostics for rifampin resistance has resulted in a steady rise in the number of patients with multidrug-resistant (MDR)-TB referred for treatment. Pyrazinamide is used in susceptible TB treatment for 6 months when used in combination with rifampin, isoniazid and ethambutol and is an important companion drug in novel MDR-TB trials. This study was undertaken to determine the prevalence of pyrazinamide resistance by either phenotypic or pncA testing among patients admitted to a referral hospital in Tanzania for drug-susceptible and MDR-TB treatment.

METHODS

Surveillance sputa were sent among subjects beginning TB therapy at the national MDR-TB referral hospital during a 6 month period in 2013-2014. Mycobacterial cultures of pretreatment sputa were performed at the Kilimanjaro Clinical Research Institute (KCRI) in the BACTEC mycobacterial growth indicator tubes (MGIT) 960 system. Speciation of M. tuberculosis complex was confirmed by MTBc assay. Isolates were sub-cultured on to Lowenstein-Jensen (LJ) slants. Phenotypic resistance to pyrazinamide was performed in the MGIT system while a real-time PCR with High Resolution Melt (HRM) technique was used to determine mutation in the pncA gene from the same pure subculture. Sputa were then collected monthly to determine the time to culture negativity. Final treatment outcome was determined.

RESULTS

Ninety-one M. tuberculosis isolates from individual patients were available for analysis of which 30 (32.9%) had MDR-TB, the mean (±SD) age was 33 ± 10 years, and the majority 23 (76.7%) were males. Of the 30 MDR-TB patients, 15(50%) had isolates with pyrazinamide resistance by conventional MGIT testing. This proportion expectedly exceeded the number with pyrazinamide resistance in the 61 patients without MDR-TB, 13 (21.3%) (p = 0.008). Six (20%) of MDR-TB patients had a poor outcome including treatment failure. Among patients with treatment failure, 5 (83%) had pyrazinamide resistance compared to only 10 (41.6%) with treatment success (p = 0.08). Two patients died, and both had isolates with pyrazinamide resistance. No other pretreatment characteristic was associated with treatment outcome.

CONCLUSION

Pyrazinamide susceptibility appears to be important in clinical outcomes for MDR-TB patients, and susceptibility testing appears to be a critical adjunct to TB care. The high proportion of PZA resistance in non-MDR TB cases calls for further local investigation.