Pyrazinamide-resistant Tuberculosis Obscured From Common Targeted Molecular Diagnostics.

Here, we describe a clinical case of pyrazinamide-resistant (PZA-R) tuberculosis (TB) reported as PZA-susceptible (PZA-S) by common molecular diagnostics. Phenotypic susceptibility testing (pDST) indicated PZA-R TB. Targeted Sanger sequencing reported wild-type PncA, indicating PZA-S TB. Whole Genome Sequencing (WGS) by PacBio and IonTorrent both detected deletion of a large portion of pncA, indicating PZA-R. Importantly, both WGS methods showed deletion of part of the primer region targeted by Sanger sequencing. Repeating Sanger sequencing from a culture in presence of PZA returned no result, revealing that 1) two minority susceptible subpopulations had vanished, 2) the PZA-R majority subpopulation harboring the pncA deletion could not be amplified by Sanger primers, and was thus obscured by amplification process. This case demonstrates how a small susceptible subpopulation can entirely obscure majority resistant populations from targeted molecular diagnostics and falsely imply homogenous susceptibility, leading to incorrect diagnosis. To our knowledge, this is the first report of a minority susceptible subpopulation masking a majority resistant population, causing targeted molecular diagnostics to call false susceptibility. The consequence of such genomic events is not limited to PZA. This phenomenon can impact molecular diagnostics' sensitivity whenever the resistance-conferring mutation is not fully within primer-targeted regions. This can be caused by structural changes of genomic context with phenotypic consequence as we report here, or by uncommon mechanisms of resistance. Such false susceptibility calls promote suboptimal treatment and spread of strains that challenge targeted molecular diagnostics. This motivates development of molecular diagnostics unreliant on primer conservation, and impels frequent WGS surveillance for variants that evade prevailing molecular diagnostics.

Evaluating impacts of syntenic block detection strategies on rearrangement phylogeny using Mycobacterium tuberculosis isolates.

MOTIVATION: The phylogenetic signal of structural variation informs a more comprehensive understanding of evolution. As (near-)complete genome assembly becomes more commonplace, the next methodological challenge for inferring genome rearrangement trees is the identification of syntenic blocks of orthologous sequences. In this article, we studied 94 reference quality genomes of primarily Mycobacterium tuberculosis (Mtb) isolates as a benchmark to evaluate these methods. The clonal nature of Mtb evolution, the manageable genome sizes, along with substantial levels of structural variation make this an ideal benchmarking dataset. RESULTS: We tested several methods for detecting homology and obtaining syntenic blocks and two methods for inferring phylogenies from them, then compared the resulting trees to the standard method's tree, inferred from nucleotide substitutions. We found that, not only the choice of methods, but also their parameters can impact results, and that the tree inference method had less impact than the block determination method. Interestingly, a rearrangement tree based on blocks from the Cactus whole-genome aligner was fully compatible with the highly supported branches of the substitution-based tree, enabling the combination of the two into a high-resolution supertree. Overall, our results indicate that accurate trees can be inferred using genome rearrangements, but the choice of the methods for inferring homology requires care. AVAILABILITY AND IMPLEMENTATION: Analysis scripts and code written for this study are available at https://gitlab.com/LPCDRP/rearrangement-homology.pub and https://gitlab.com/LPCDRP/syntement. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel and reported compensatory mutations in rpoABC genes found in drug resistant tuberculosis outbreaks.

Background: Rifampicin (RIF) is a key first-line drug used to treat tuberculosis, a primarily pulmonary disease caused by Mycobacterium tuberculosis. RIF resistance is caused by mutations in rpoB, at the cost of slower growth and reduced transcription efficiency. Antibiotic resistance to RIF is prevalent despite this fitness cost. Compensatory mutations in rpoABC genes have been shown to alleviate the fitness cost of rpoB:S450L, explaining how RIF resistant strains harbor this mutation can spread so rapidly. Unfortunately, the full set of RIF compensatory mutations is still unknown, particularly those compensating for rarer RIF resistance mutations. Objectives: We performed an association study on a globally representative set of 4,309 whole genome sequenced clinical M. tuberculosis isolates to identify novel putative compensatory mutations, determine the prevalence of known and previously reported putative compensatory mutations, and determine which RIF resistance markers associate with these compensatory mutations. Results and conclusions: Of the 1,079 RIF resistant isolates, 638 carried previously reported putative and high-probability compensatory mutations. Our strict criteria identified 46 additional mutations in rpoABC for which no strong prior evidence of their compensatory role exists. Of these, 35 have previously been reported. As such, our independent corroboration adds to the mounting evidence that these 35 also carry a compensatory role. The remaining 11 are novel putative compensatory markers, reported here for the first time. Six of these 11 novel putative compensatory mutations had two or more mutation events. Most compensatory mutations appear to be specifically compensating for the fitness loss due to rpoB:S450L. However, an outbreak of 22 closely related isolates each carried three rpoB mutations, the rare RIFR markers D435G and L452P and the putative compensatory mutation I1106T. This suggests compensation may require specific combinations of rpoABC mutations. Here, we report only mutations that met our very strict criteria. It is highly likely that many additional rpoABC mutations compensate for rare resistance-causing mutations and therefore did not carry the statistical power to be reported here. These findings aid in the identification of RIF resistant M. tuberculosis strains with restored fitness, which pose a greater risk of causing resistant outbreaks.

Distribution of Common and Rare Genetic Markers of Second-Line-Injectable-Drug Resistance in Mycobacterium tuberculosis Revealed by a Genome-Wide Association Study.

Point mutations in the rrs gene and the eis promoter are known to confer resistance to the second-line injectable drugs (SLIDs) amikacin (AMK), capreomycin (CAP), and kanamycin (KAN). While mutations in these canonical genes confer the majority of SLID resistance, alternative mechanisms of resistance are not uncommon and threaten effective treatment decisions when using conventional molecular diagnostics. In total, 1,184 clinical Mycobacterium tuberculosis isolates from 7 countries were studied for genomic markers associated with phenotypic resistance. The markers rrs:A1401G and rrs:G1484T were associated with resistance to all three SLIDs, and three known markers in the eis promoter (eis:G-10A, eis:C-12T, and eis:C-14T) were similarly associated with kanamycin resistance (KAN-R). Among 325, 324, and 270 AMK-R, CAP-R, and KAN-R isolates, 274 (84.3%), 250 (77.2%), and 249 (92.3%) harbored canonical mutations, respectively. Thirteen isolates harbored more than one canonical mutation. Canonical mutations did not account for 103 of the phenotypically resistant isolates. A genome-wide association study identified three genes and promoters with mutations that, on aggregate, were associated with unexplained resistance to at least one SLID. Our analysis associated whiB7 5'-untranslated-region mutations with KAN resistance, supporting clinical relevance for this previously demonstrated mechanism of KAN resistance. We also provide evidence for the novel association of CAP resistance with the promoter of the Rv2680-Rv2681 operon, which encodes an exoribonuclease that may influence the binding of CAP to the ribosome. Aggregating mutations by gene can provide additional insight and therefore is recommended for identifying rare mechanisms of resistance when individual mutations carry insufficient statistical power.

Structure-Aware Mycobacterium tuberculosis Functional Annotation Uncloaks Resistance, Metabolic, and Virulence Genes.

Accurate and timely functional genome annotation is essential for translating basic pathogen research into clinically impactful advances. Here, through literature curation and structure-function inference, we systematically update the functional genome annotation of Mycobacterium tuberculosis virulent type strain H37Rv. First, we systematically curated annotations for 589 genes from 662 publications, including 282 gene products absent from leading databases. Second, we modeled 1,711 underannotated proteins and developed a semiautomated pipeline that captured shared function between 400 protein models and structural matches of known function on Protein Data Bank, including drug efflux proteins, metabolic enzymes, and virulence factors. In aggregate, these structure- and literature-derived annotations update 940/1,725 underannotated H37Rv genes and generate hundreds of functional hypotheses. Retrospectively applying the annotation to a recent whole-genome transposon mutant screen provided missing function for 48% (13/27) of underannotated genes altering antibiotic efficacy and 33% (23/69) required for persistence during mouse tuberculosis (TB) infection. Prospective application of the protein models enabled us to functionally interpret novel laboratory generated pyrazinamide (PZA)-resistant mutants of unknown function, which implicated the emerging coenzyme A depletion model of PZA action in the mutants' PZA resistance. Our findings demonstrate the functional insight gained by integrating structural modeling and systematic literature curation, even for widely studied microorganisms. Functional annotations and protein structure models are available at https://tuberculosis.sdsu.edu/H37Rv in human- and machine-readable formats. IMPORTANCE Mycobacterium tuberculosis, the primary causative agent of tuberculosis, kills more humans than any other infectious bacterium. Yet 40% of its genome is functionally uncharacterized, leaving much about the genetic basis of its resistance to antibiotics, capacity to withstand host immunity, and basic metabolism yet undiscovered. Irregular literature curation for functional annotation contributes to this gap. We systematically curated functions from literature and structural similarity for over half of poorly characterized genes, expanding the functionally annotated Mycobacterium tuberculosis proteome. Applying this updated annotation to recent in vivo functional screens added functional information to dozens of clinically pertinent proteins described as having unknown function. Integrating the annotations with a prospective functional screen identified new mutants resistant to a first-line TB drug, supporting an emerging hypothesis for its mode of action. These improvements in functional interpretation of clinically informative studies underscore the translational value of this functional knowledge. Structure-derived annotations identify hundreds of high-confidence candidates for mechanisms of antibiotic resistance, virulence factors, and basic metabolism and other functions key in clinical and basic tuberculosis research. More broadly, they provide a systematic framework for improving prokaryotic reference annotations.

Atypical Genetic Basis of Pyrazinamide Resistance in Monoresistant Mycobacterium tuberculosis.

Pyrazinamide (PZA) is a widely used antitubercular chemotherapeutic. Typically, PZA resistance (PZA-R) emerges in Mycobacterium tuberculosis strains with existing resistance to isoniazid and rifampin (i.e., multidrug resistance [MDR]) and is conferred by loss-of-function pncA mutations that inhibit conversion to its active form, pyrazinoic acid (POA). PZA-R departing from this canonical scenario is poorly understood. Here, we genotyped pncA and purported alternative PZA-R genes (panD, rpsA, and clpC1) with long-read sequencing of 19 phenotypically PZA-monoresistant isolates collected in Sweden and compared their phylogenetic and genomic characteristics to a large set of MDR PZA-R (MDRPZA-R) isolates. We report the first association of ClpC1 mutations with PZA-R in clinical isolates, in the ClpC1 promoter (clpC1p-138) and the N terminus of ClpC1 (ClpC1Val63Ala). Mutations have emerged in both these regions under POA selection in vitro, and the N-terminal region of ClpC1 has been implicated further, through its POA-dependent efficacy in PanD proteolysis. ClpC1Val63Ala mutants spanned 4 Indo-Oceanic sublineages. Indo-Oceanic isolates invariably harbored ClpC1Val63Ala and were starkly overrepresented (odds ratio [OR] = 22.2, P < 0.00001) among PZA-monoresistant isolates (11/19) compared to MDRPZA-R isolates (5/80). The genetic basis of Indo-Oceanic isolates' overrepresentation in PZA-monoresistant tuberculosis (TB) remains undetermined, but substantial circumstantial evidence suggests that ClpC1Val63Ala confers low-level PZA resistance. Our findings highlight ClpC1 as potentially clinically relevant for PZA-R and reinforce the importance of genetic background in the trajectory of resistance development.

Exact mapping of Illumina blind spots in the Mycobacterium tuberculosis genome reveals platform-wide and workflow-specific biases.

Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina 'blind spots' in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the 'Illumina-sequenceable' genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.

Pyrazinamide Susceptibility Is Driven by Activation of the SigE-Dependent Cell Envelope Stress Response in Mycobacterium tuberculosis.

Pyrazinamide (PZA) plays a crucial role in first-line tuberculosis drug therapy. Unlike other antimicrobial agents, PZA is active against Mycobacterium tuberculosis only at low pH. The basis for this conditional drug susceptibility remains undefined. In this study, we utilized a genome-wide approach to interrogate potentiation of PZA action. We found that mutations in numerous genes involved in central metabolism as well as cell envelope maintenance and stress response are associated with PZA resistance. Further, we demonstrate that constitutive activation of the cell envelope stress response can drive PZA susceptibility independent of environmental pH. Consequently, exposure to peptidoglycan synthesis inhibitors, such as beta-lactams and d-cycloserine, potentiate PZA action through triggering this response. These findings illuminate a regulatory mechanism for conditional PZA susceptibility and reveal new avenues for enhancing potency of this important drug through targeting activation of the cell envelope stress response. IMPORTANCE For decades, pyrazinamide has served as a cornerstone of tuberculosis therapy. Unlike any other antitubercular drug, pyrazinamide requires an acidic environment to exert its action. Despite its importance, the driver of this conditional susceptibility has remained unknown. In this study, a genome-wide approach revealed that pyrazinamide action is governed by the cell envelope stress response. This observation was validated by orthologous approaches that demonstrate that a central player of this response, SigE, is both necessary and sufficient for potentiation of pyrazinamide action. Moreover, constitutive activation of this response through deletion of the anti-sigma factor gene rseA or exposure of bacilli to drugs that target the cell wall was found to potently drive pyrazinamide susceptibility independent of environmental pH. These findings force a paradigm shift in our understanding of pyrazinamide action and open new avenues for improving diagnostic and therapeutic tools for tuberculosis.

Drivers and sites of diversity in the DNA adenine methylomes of 93 Mycobacterium tuberculosis complex clinical isolates.

This study assembles DNA adenine methylomes for 93 Mycobacterium tuberculosis complex (MTBC) isolates from seven lineages paired with fully-annotated, finished, de novo assembled genomes. Integrative analysis yielded four key results. First, methyltransferase allele-methylome mapping corrected methyltransferase variant effects previously obscured by reference-based variant calling. Second, heterogeneity analysis of partially active methyltransferase alleles revealed that intracellular stochastic methylation generates a mosaic of methylomes within isogenic cultures, which we formalize as 'intercellular mosaic methylation' (IMM). Mutation-driven IMM was nearly ubiquitous in the globally prominent Beijing sublineage. Third, promoter methylation is widespread and associated with differential expression in the ΔhsdM transcriptome, suggesting promoter HsdM-methylation directly influences transcription. Finally, comparative and functional analyses identified 351 sites hypervariable across isolates and numerous putative regulatory interactions. This multi-omic integration revealed features of methylomic variability in clinical isolates and provides a rational basis for hypothesizing the functions of DNA adenine methylation in MTBC physiology and adaptive evolution.

Publisher Correction: Genomic analysis of the emergence of drug-resistant strains of Mycobacterium tuberculosis in the Middle East.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Functional Connectivities Are More Informative Than Anatomical Variables in Diagnostic Classification of Autism.

Machine learning techniques have been implemented to reveal brain features that distinguish people with autism spectrum disorders (ASDs) from typically developing (TD) peers. However, it remains unknown whether different neuroimaging modalities are equally informative for diagnostic classification. We combined anatomical magnetic resonance imaging (aMRI), diffusion weighted imaging (DWI), and functional connectivity MRI (fcMRI) using conditional random forest (CRF) for supervised learning to compare how informative each modality was in diagnostic classification. In-house data (N = 93) included 47 TD and 46 ASD participants, matched on age, motion, and nonverbal IQ. Four main analyses consistently indicated that fcMRI variables were significantly more informative than anatomical variables from aMRI and DWI. This was found (1) when the top 100 variables from CRF (run separately in each modality) were combined for multimodal CRF; (2) when only 19 top variables reaching >67% accuracy in each modality were combined in multimodal CRF; and (3) when the large number of initial variables (before dimension reduction) potentially biasing comparisons in favor of fcMRI was reduced using a less granular region of interest scheme. Consistent superiority of fcMRI was even found (4) when 100 variables per modality were randomly selected, removing any such potential bias. Greater informative value of functional than anatomical modalities may relate to the nature of fcMRI data, reflecting more closely behavioral condition, which is also the basis of diagnosis, whereas brain anatomy may be more reflective of neurodevelopmental history.

Genomic analysis of the emergence of drug-resistant strains of Mycobacterium tuberculosis in the Middle East.

Tuberculosis (TB) represents a significant challenge to public health authorities, especially with the emergence of drug-resistant (DR) and multidrug-resistant (MDR) isolates of Mycobacterium tuberculosis. We sought to examine the genomic variations among recently isolated strains of M. tuberculosis in two closely related countries with different population demography in the Middle East. Clinical isolates of M. tuberculosis from both Egypt and Saudi Arabia were subjected to phenotypic and genotypic analysis on gene and genome-wide levels. Isolates with MDR phenotypes were highly prevalent in Egypt (up to 35%) despite its relatively stable population structure (sympatric pattern). MDR-TB isolates were not identified in the isolates from Saudi Arabia despite its active guest worker program (allopatric pattern). However, tuberculosis isolates from Saudi Arabia, where lineage 4 was more prevalent (>65%), showed more diversity than isolates from Egypt, where lineage 3 was the most prevalent (>75%). Phylogenetic and molecular dating analyses indicated that lineages from Egypt were recently diverged (~78 years), whereas those from Saudi Arabia were diverged by over 200 years. Interestingly, DR isolates did not appear to cluster together or spread more widely than drug-sensitive isolates, suggesting poor treatment as the main cause for emergence of drug resistance rather than more virulence or more capacity to persist.

Modeling the structural origins of drug resistance to isoniazid via key mutations in Mycobacterium tuberculosis catalase-peroxidase, KatG.

WHO reported 10.4 million new tuberculosis (TB) cases and 1.8 million deaths in 2015, making M. tuberculosis the most successful human pathogen with highest mortality among infectious diseases [1,2]. Drug-resistant TB is a major threat to global TB control [2,3]. Recently Torres et al. [4] identified 14 novel substitutions in M. tuberculosis-KatG (the enzyme associated with resistance to isoniazid-an important first-line anti-TB drug) and demonstrated that 12 of the 14 can cause INH-resistance in M. smegmatis. This study presents an in silico structure-based analysis of these 14 amino acid substitutions using homology models and x-ray crystal structures (when available) in M. tuberculosis. Our models demonstrate that several of these mutations cluster around three openings in the KatG tertiary structure which appear to initiate channels to the heme group at the catalytic center of the enzyme. We studied the effects of these mutations on the tertiary structure of KatG, focusing on conformational changes in the three channels in the protein structure. Our results suggest that the 14 novel mutations sufficiently restrict one or more of these access channels, thus potentially preventing INH from reaching the catalytic heme. These observations provide valuable insights into the structure-based origins of INH resistance and provide testable hypotheses for future experimental studies.

SMRT genome assembly corrects reference errors, resolving the genetic basis of virulence in Mycobacterium tuberculosis.

BACKGROUND: The genetic basis of virulence in Mycobacterium tuberculosis has been investigated through genome comparisons of virulent (H37Rv) and attenuated (H37Ra) sister strains. Such analysis, however, relies heavily on the accuracy of the sequences. While the H37Rv reference genome has had several corrections to date, that of H37Ra is unmodified since its original publication. RESULTS: Here, we report the assembly and finishing of the H37Ra genome from single-molecule, real-time (SMRT) sequencing. Our assembly reveals that the number of H37Ra-specific variants is less than half of what the Sanger-based H37Ra reference sequence indicates, undermining and, in some cases, invalidating the conclusions of several studies. PE_PPE family genes, which are intractable to commonly-used sequencing platforms because of their repetitive and GC-rich nature, are overrepresented in the set of genes in which all reported H37Ra-specific variants are contradicted. Further, one of the sequencing errors in H37Ra masks a true variant in common with the clinical strain CDC1551 which, when considered in the context of previous work, corresponds to a sequencing error in the H37Rv reference genome. CONCLUSIONS: Our results constrain the set of genomic differences possibly affecting virulence by more than half, which focuses laboratory investigation on pertinent targets and demonstrates the power of SMRT sequencing for producing high-quality reference genomes.

Structure-based analysis of Bacilli and plasmid dihydrofolate reductase evolution.

Dihydrofolate reductase (DHFR), a key enzyme in tetrahydrofolate-mediated biosynthetic pathways, has a structural motif known to be highly conserved over a wide range of organisms. Given its critical role in purine and amino acid synthesis, DHFR is a well established therapeutic target for treating a wide range of prokaryotic and eukaryotic infections as well as certain types of cancer. Here we present a structural-based computer analysis of bacterial (Bacilli) and plasmid DHFR evolution. We generated a structure-based sequence alignment using 7 wild-type DHFR x-ray crystal structures obtained from the RCSB Protein Data Bank and 350 chromosomal and plasmid homology models we generated from sequences obtained from the NCBI Protein Database. We used these alignments to compare active site and non-active site conservation in terms of amino acid residues, secondary structure and amino acid residue class. With respect to amino acid sequences and residue classes, active-site positions in both plasmid and chromosomal DHFR are significantly more conserved than non-active site positions. Secondary structure conservation was similar for active site and non-active site positions. Plasmid-encoded DHFR proteins have greater degree of sequence and residue class conservation, particularly in sequence positions associated with a network of concerted protein motions, than chromosomal-encoded DHFR proteins. These structure-based were used to build DHFR specific phylogenetic trees from which evidence for horizontal gene transfer was identified.

Performance Comparison of Three Rapid Tests for the Diagnosis of Drug-Resistant Tuberculosis.

BACKGROUND: The aim of this study was to compare the performance of several recently developed assays for the detection of multi- and extensively drug-resistant tuberculosis (M/XDR-TB) in a large, multinational field trial. METHODS: Samples from 1,128 M/XDR-TB suspects were examined by Line Probe Assay (LPA), Pyrosequencing (PSQ), and Microscopic Observation of Drug Susceptibility (MODS) and compared to the BACTEC MGIT960 reference standard to detect M/XDR-TB directly from patient sputum samples collected at TB clinics in India, Moldova, and South Africa. RESULTS: Specificity for all three assays was excellent: 97-100% for isoniazid (INH), rifampin (RIF), moxifloxacin (MOX) and ofloxacin (OFX) and 99-100% for amikacin (AMK), capreomycin (CAP) and kanamycin (KAN) resistance. Sensitivities were lower, but still very good: 94-100% for INH, RIF, MOX and OFX, and 84-90% for AMK and CAP, but only 48-62% for KAN. In terms of agreement, statistically significant differences were only found for detection of RIF (MODS outperformed PSQ) and KAN (MODS outperformed LPA and PSQ) resistance. Mean time-to-result was 1.1 days for LPA and PSQ, 14.3 days for MODS, and 24.7 days for MGIT. CONCLUSIONS: All three rapid assays evaluated provide clinicians with timely detection of resistance to the drugs tested; with molecular results available one day following laboratory receipt of samples. In particular, the very high specificity seen for detection of drug resistance means that clinicians can use the results of these rapid tests to avoid the use of toxic drugs to which the infecting organism is resistant and develop treatment regiments that have a higher likelihood of yielding a successful outcome.

Novel katG mutations causing isoniazid resistance in clinical M. tuberculosis isolates.

We report the discovery and confirmation of 23 novel mutations with previously undocumented role in isoniazid (INH) drug resistance, in catalase-peroxidase (katG) gene of Mycobacterium tuberculosis (Mtb) isolates. With these mutations, a synonymous mutation in fabG1 (g609a), and two canonical mutations, we were able to explain 98% of the phenotypic resistance observed in 366 clinical Mtb isolates collected from four high tuberculosis (TB)-burden countries: India, Moldova, Philippines, and South Africa. We conducted overlapping targeted and whole-genome sequencing for variant discovery in all clinical isolates with a variety of INH-resistant phenotypes. Our analysis showed that just two canonical mutations (katG 315AGC-ACC and inhA promoter-15C-T) identified 89.5% of resistance phenotypes in our collection. Inclusion of the 23 novel mutations reported here, and the previously documented point mutation in fabG1, increased the sensitivity of these mutations as markers of INH resistance to 98%. Only six (2%) of the 332 resistant isolates in our collection did not harbor one or more of these mutations. The third most prevalent substitution, at inhA promoter position -8, present in 39 resistant isolates, was of no diagnostic significance since it always co-occurred with katG 315. 79% of our isolates harboring novel mutations belong to genetic group 1 indicating a higher tendency for this group to go down an uncommon evolutionary path and evade molecular diagnostics. The results of this study contribute to our understanding of the mechanisms of INH resistance in Mtb isolates that lack the canonical mutations and could improve the sensitivity of next generation molecular diagnostics.

Systematic review of mutations in pyrazinamidase associated with pyrazinamide resistance in Mycobacterium tuberculosis clinical isolates.

Pyrazinamide (PZA) is an important first-line drug in the treatment of tuberculosis (TB) and of significant interest to the HIV-infected community due to the prevalence of TB-HIV coinfection in some regions of the world. The mechanism of resistance to PZA is unlike that of any other anti-TB drug. The gene pncA, encoding pyrazinamidase (PZase), is associated with resistance to PZA. However, because single mutations in PZase have a low prevalence, the individual sensitivities are low. Hundreds of distinct mutations in the enzyme have been associated with resistance, while some only appear in susceptible isolates. This makes interpretation of molecular testing difficult and often leads to the simplification that any PZase mutation causes resistance. This systematic review reports a comprehensive global list of mutations observed in PZase and its promoter region in clinical strains, their phenotypic association, their global frequencies and diversity, the method of phenotypic determination, their MIC values when given, and the method of MIC determination and assesses the strength of the association between mutations and phenotypic resistance to PZA. In this systematic review, we report global statistics for 641 mutations in 171 (of 187) codons from 2,760 resistant strains and 96 mutations from 3,329 susceptible strains reported in 61 studies. For diagnostics, individual mutations (or any subset) were not sufficiently sensitive. Assuming similar error profiles of the 5 phenotyping platforms included in this study, the entire enzyme and its promoter provide a combined estimated sensitivity of 83%. This review highlights the need for identification of an alternative mechanism(s) of resistance, at least for the unexplained 17% of cases.

Detection of Low-Level Mixed-Population Drug Resistance in Mycobacterium tuberculosis Using High Fidelity Amplicon Sequencing.

Undetected and untreated, low-levels of drug resistant (DR) subpopulations in clinical Mycobacterium tuberculosis (Mtb) infections may lead to development of DR-tuberculosis, potentially resulting in treatment failure. Current phenotypic DR susceptibility testing has a theoretical potential for 1% sensitivity, is not quantitative, and requires several weeks to complete. The use of "single molecule-overlapping reads" (SMOR) analysis with next generation DNA sequencing for determination of ultra-rare target alleles in complex mixtures provides increased sensitivity over standard DNA sequencing. Ligation free amplicon sequencing with SMOR analysis enables the detection of resistant allele subpopulations at ≥0.1% of the total Mtb population in near real-time analysis. We describe the method using standardized mixtures of DNA from resistant and susceptible Mtb isolates and the assay's performance for detecting ultra-rare DR subpopulations in DNA extracted directly from clinical sputum samples. SMOR analysis enables rapid near real-time detection and tracking of previously undetectable DR sub-populations in clinical samples allowing for the evaluation of the clinical relevance of low-level DR subpopulations. This will provide insights into interventions aimed at suppressing minor DR subpopulations before they become clinically significant.