Putative compensatory mutations in the rpoC gene of rifampin-resistant Mycobacterium tuberculosis are associated with ongoing transmission

Large-scale statistical analysis of Mycobacterium tuberculosis genome sequences identifies compensatory mutations associated with multi-drug resistance

Tuberculosis (TB), caused by Mycobacterium tuberculosis, has a significant impact on global health worldwide. The development of multi-drug resistant strains that are resistant to the first-line drugs isoniazid and rifampicin threatens public health security. Rifampicin and isoniazid resistance are largely underpinned by mutations in rpoB and katG respectively and are associated with fitness costs. Compensatory mutations are considered to alleviate these fitness costs and have been observed in rpoC/rpoA (rifampicin) and oxyR'-ahpC (isoniazid). We developed a framework (CompMut-TB) to detect compensatory mutations from whole genome sequences from a large dataset comprised of 18,396 M. tuberculosis samples. We performed association analysis (Fisher's exact tests) to identify pairs of mutations that are associated with drug-resistance, followed by mediation analysis to identify complementary or full mediators of drug-resistance. The analyses revealed several potential mutations in rpoC (N = 47), rpoA (N = 4), and oxyR'-ahpC (N = 7) that were considered either 'highly likely' or 'likely' to confer compensatory effects on drug-resistance, including mutations that have previously been reported and validated. Overall, we have developed the CompMut-TB framework which can assist with identifying compensatory mutations which is important for more precise genome-based profiling of drug-resistant TB strains and to further understanding of the evolutionary mechanisms that underpin drug-resistance.

Genetic diversity, evolution and drug resistance of Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis causes a chronic infectious disease called tuberculosis. Phylogenetic lineage 2 (L2) of M. tuberculosis, also known as the East Asian lineage, is associated with high virulence, increased transmissibility, and the spread of multidrug-resistant strains. This review article examines the genomic characteristics of the M. tuberculosis genome and M. tuberculosis lineage 2, such as the unique insertion sequence and spoligotype patterns, as well as MIRU-VNTR typing, and SNP-based barcoding. The review describes the geographical distribution of lineage 2 and its history of origin. In addition, the article discusses recent studies on drug resistance and compensatory mechanisms of M. tuberculosis lineage 2 and its impact on the pathogen's transmissibility and virulence. This review article discusses the importance of establishing a unified classification for lineage 2 to ensure consistency in terminology and criteria across different studies and settings.

A dose-response model for statistical analysis of chemical genetic interactions in CRISPRi screens

An important application of CRISPR interference (CRISPRi) technology is for identifying chemical-genetic interactions (CGIs). Discovery of genes that interact with exposure to antibiotics can yield insights to drug targets and mechanisms of action or resistance. The objective is to identify CRISPRi mutants whose relative abundance is suppressed (or enriched) in the presence of a drug when the target protein is depleted, reflecting synergistic behavior. Different sgRNAs for a given target can induce a wide range of protein depletion and differential effects on growth rate. The effect of sgRNA strength can be partially predicted based on sequence features. However, the actual growth phenotype depends on the sensitivity of cells to depletion of the target protein. For essential genes, sgRNA efficiency can be empirically measured by quantifying effects on growth rate. We observe that the most efficient sgRNAs are not always optimal for detecting synergies with drugs. sgRNA efficiency interacts in a non-linear way with drug sensitivity, producing an effect where the concentration-dependence is maximized for sgRNAs of intermediate strength (and less so for sgRNAs that induce too much or too little target depletion). To capture this interaction, we propose a novel statistical method called CRISPRi-DR (for Dose-Response model) that incorporates both sgRNA efficiencies and drug concentrations in a modified dose-response equation. We use CRISPRi-DR to re-analyze data from a recent CGI experiment in Mycobacterium tuberculosis to identify genes that interact with antibiotics. This approach can be generalized to non-CGI datasets, which we show via an CRISPRi dataset for E. coli growth on different carbon sources. The performance is competitive with the best of several related analytical methods. However, for noisier datasets, some of these methods generate far more significant interactions, likely including many false positives, whereas CRISPRi-DR maintains higher precision, which we observed in both empirical and simulated data.

Actionable mechanisms of drug tolerance and resistance in Mycobacterium tuberculosis

The emergence of antimicrobial resistance (AMR) across bacterial pathogens presents a serious threat to global health. This threat is further exacerbated in tuberculosis (TB), mainly due to a protracted treatment regimen involving a combination of drugs. A diversity of factors contributes to the emergence of drug resistance in TB, which is caused by the pathogen Mycobacterium tuberculosis (Mtb). While the traditional genetic mutation-driven drug resistance mechanisms operate in Mtb, there are also several additional unique features of drug resistance in this pathogen. Research in the past decade has enriched our understanding of such unconventional factors as efflux pumps, bacterial heterogeneity, metabolic states, and host microenvironment. Given that the discovery of new antibiotics is outpaced by the emergence of drug resistance patterns displayed by the pathogen, newer strategies for combating drug resistance are desperately needed. In the context of TB, such approaches include targeting the efflux capability of the pathogen, modulating the host environment to prevent bacterial drug tolerance, and activating the host anti-mycobacterial pathways. In this review, we discuss the traditional mechanisms of drug resistance in Mtb, newer understandings and the shaping of a set of unconventional approaches to target both the emergence and treatment of drug resistance in TB.

Compensatory mutations are associated with increased in vitro growth in resistant clinical samples of Mycobacterium tuberculosis

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

Laboratory Evolution of Antimicrobial Resistance in Bacteria to Develop Rational Treatment Strategies

Laboratory evolution studies, particularly with Escherichia coli, have yielded invaluable insights into the mechanisms of antimicrobial resistance (AMR). Recent investigations have illuminated that, with repetitive antibiotic exposures, bacterial populations will adapt and eventually become tolerant and resistant to the drugs. Through intensive analyses, these inquiries have unveiled instances of convergent evolution across diverse antibiotics, the pleiotropic effects of resistance mutations, and the role played by loss-of-function mutations in the evolutionary landscape. Moreover, a quantitative analysis of multidrug combinations has shed light on collateral sensitivity, revealing specific drug combinations capable of suppressing the acquisition of resistance. This review article introduces the methodologies employed in the laboratory evolution of AMR in bacteria and presents recent discoveries concerning AMR mechanisms derived from laboratory evolution. Additionally, the review outlines the application of laboratory evolution in endeavors to formulate rational treatment strategies.

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.

How do Mutations of Mycobacterium Genes Cause Drug Resistance in Tuberculosis?

A steady increase in the prevalence of drug-resistant tuberculosis (DR-TB) has already been reported in Pakistan. In addition, DR-TB is gradually changing from one-drug resistance to multi-drug resistance, which is a serious challenge for tuberculosis treatment. This review provides an overview of the anti-tuberculosis drugs and focuses on the molecular mechanisms of drug resistance in Mycobacterium tuberculosis, with the hope that it will contribute to the study of drug resistance in response to the emergence of multidrug-resistant tuberculosis.

Pan-genome de Bruijn graph using the bidirectional FM-index

BACKGROUND

Pan-genome graphs are gaining importance in the field of bioinformatics as data structures to represent and jointly analyze multiple genomes. Compacted de Bruijn graphs are inherently suited for this purpose, as their graph topology naturally reveals similarity and divergence within the pan-genome. Most state-of-the-art pan-genome graphs are represented explicitly in terms of nodes and edges. Recently, an alternative, implicit graph representation was proposed that builds directly upon the unidirectional FM-index. As such, a memory-efficient graph data structure is obtained that inherits the FM-index' backward search functionality. However, this representation suffers from a number of shortcomings in terms of functionality and algorithmic performance.

RESULTS

We present a data structure for a pan-genome, compacted de Bruijn graph that aims to address these shortcomings. It is built on the bidirectional FM-index, extending the ability of its unidirectional counterpart to navigate and search the graph in both directions. All basic graph navigation steps can be performed in constant time. Based on these features, we implement subgraph visualization as well as lossless approximate pattern matching to the graph using search schemes. We demonstrate that we can retrieve all occurrences corresponding to a read within a certain edit distance in a very efficient manner. Through a case study, we show the potential of exploiting the information embedded in the graph's topology through visualization and sequence alignment.

CONCLUSIONS

We propose a memory-efficient representation of the pan-genome graph that supports subgraph visualization and lossless approximate pattern matching of reads against the graph using search schemes. The C++ source code of our software, called Nexus, is available at https://github.com/biointec/nexus under AGPL-3.0 license.

Genomic Characterization of Drug-Resistant Mycobacterium tuberculosis L2/Beijing Isolates from Astana, Kazakhstan

Kazakhstan ranks among the countries with the highest number of MDR-TB patients per 100,000 population worldwide. The successful transmission of local MDR strains of Mycobacterium tuberculosis (Mtb) poses a significant threat to disease control. In this study, we employed whole-genome sequencing to examine drug resistance, compensatory mutations, population structure, and transmission patterns in a sample of 24 clinical isolates of L2/Beijing Mtb collected in Astana, Kazakhstan between 2021 and 2022. The genotypic prediction of Mtb susceptibility to anti-TB agents was consistent with the phenotypic susceptibility, except for bedaquiline. An analysis of resistance-associated genes characterized most of the isolates as pre-extensively drug-resistant tuberculosis (pre-XDR-TB) (n = 15; 62.5%). The phylogenetic analysis grouped the isolates into four transmission clusters; the dominant cluster was assigned to the "aggressive" Central Asia outbreak (CAO) clade of L2/Beijing (n = 15; 62.5%). Thirteen mutations with putative compensatory effects were observed exclusively in Mtb isolates containing the rpoB S450L mutation. The putative compensatory mutations had a stabilizing effect on RpoABC protein stability and dynamics. The high prevalence of the CAO clade in the population structure of Mtb may explain the rapid spread of MDR-TB in Kazakhstan.

In silico structural and mechanical insights into bedaquiline resistance associated with high-grade non-synonymous mutations in atpE, mmpR5, and pepQ

Clinical resistance against bedaquiline (BDQ) remains intractable to anti-tuberculosis therapies since its introduction to the market over a decade ago. Herein, we investigated the structural and mechanical aspects of BDQ resistance in AtpE, MmpR5, and PepQ. The known target-specific resistant single non-synonymous mutations were refined to high-grade candidates. Thus, 7 (AtpE), 5 (MmpR5), and 1 (PepQ) single nucleotide polymorphisms (SNPs) and one insertion frameshift mutation in MmpR5 were recreated at the molecular level, and these phenotypic models were then directed to stringent dynamics to define time-scaled changes. The AtpE variants destabilized the structure; mainly, L59V, E61D, and I66M were detrimental to the complex fitness, while L74V and L114P boosted the BDQ binding to MmpR5. The first three and last two alterations gave rise to loss- and gain-of-function to AtpE and MmpR5, respectively. Hence, these five mutants are functionally relevant and therapeutically targetable hotspots of BDQ resistance. There were no noticeable changes in PepQ data analysis. The present study revealed that MmpR5 mutations confer BDQ resistance, whereas AtpE and PepQ SNPs display low susceptibility. These results were tallied with the published findings, which testified to the pursued method's reliability and accuracy. We hope these data and inferences could be helpful for the futuristic design of novel TB drugs.Communicated by Ramaswamy H. Sarma.

Molecular Analysis of Anti-Tuberculosis Drug Resistance of Mycobacterium tuberculosis Isolated in the Republic of Korea

Rapid and accurate detection of tuberculosis (TB) drug resistance is critical for the successful treatment and control of TB. Here, we investigated resistance to anti-TB drugs and genetic variations in 215 drug-resistant Mycobacterium tuberculosis isolates in Korea. Genetic variations were observed in rpoB Ser531Leu, katG Ser315Thr, and gyrA Asp94Gly; however, the minimum inhibitory concentrations varied, which can be attributed to other resistance mechanisms. Examination of genetic relatedness among drug-resistant isolates revealed that the cluster size of resistant bacteria was less than six strains, suggesting no evidence of a large-scale epidemic caused by a specific strain. However, rpoC mutants of the rifampicin-resistant isolates were composed of five types of clusters, suggesting that these compensatory mutations advance propagation. In the present study, more than 90% of the resistance mechanisms to major anti-TB drugs were identified, and the effect of each mutation on drug resistance was estimated. With the clinical application of recent next-generation sequencing-based susceptibility testing, the present study is expected to improve the clinical utilization of genotype-based drug susceptibility testing for the diagnosis and treatment of patients with drug-resistant TB.

Effect of compensatory evolution in the emergence and transmission of rifampicin-resistant Mycobacterium tuberculosis in Cape Town, South Africa: a genomic epidemiology study

BACKGROUND

Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa.

METHODS

We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains.

FINDINGS

Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08-2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28-1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13-2·12), independent of other patient and bacterial factors.

INTERPRETATION

Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented.

FUNDING

Funding for this study was provided by a Swiss and South Africa joint research award (grant numbers 310030_188888, CRSII5_177163, and IZLSZ3_170834), the European Research Council (grant number 883582), and a Wellcome Trust fellowship (to HC; reference number 099818/Z/12/Z). ZS-D was funded through a PhD scholarship from the South African National Research Foundation and RMW was funded through the South African Medical Research Council.

The relative transmission fitness of multidrug-resistant Mycobacterium tuberculosis in a drug resistance hotspot

Multidrug-resistant tuberculosis (MDR-TB) is among the most frequent causes of death due to antimicrobial resistance. Although only 3% of global TB cases are MDR, geographical hotspots with up to 40% of MDR-TB have been observed in countries of the former Soviet Union. While the quality of TB control and patient-related factors are known contributors to such hotspots, the role of the pathogen remains unclear. Here we show that in the country of Georgia, a known hotspot of MDR-TB, MDR Mycobacterium tuberculosis strains of lineage 4 (L4) transmit less than their drug-susceptible counterparts, whereas most MDR strains of L2 suffer no such defect. Our findings further indicate that the high transmission fitness of these L2 strains results from epistatic interactions between the rifampicin resistance-conferring mutation RpoB S450L, compensatory mutations in the RNA polymerase, and other pre-existing genetic features of L2/Beijing clones that circulate in Georgia. We conclude that the transmission fitness of MDR M. tuberculosis strains is heterogeneous, but can be as high as drug-susceptible forms, and that such highly drug-resistant and transmissible strains contribute to the emergence and maintenance of hotspots of MDR-TB. As these strains successfully overcome the metabolic burden of drug resistance, and given the ongoing rollout of new treatment regimens against MDR-TB, proper surveillance should be implemented to prevent these strains from acquiring resistance to the additional drugs.

Comparison of Phenotype and Genotype Virulence and Antimicrobial Factors of Salmonella Typhimurium Isolated from Human Milk

Salmonella is a common foodborne infection. Many serovars belonging to Salmonella enterica subsp. enterica are present in the gut of various animal species. They can cause infection in human infants via breast milk or cross-contamination with powdered milk. In the present study, Salmonella BO was isolated from human milk in accordance with ISO 6579-1:2017 standards and sequenced using whole-genome sequencing (WGS), followed by serosequencing and genotyping. The results also allowed its pathogenicity to be predicted. The WGS results were compared with the bacterial phenotype. The isolated strain was found to be Salmonella enterica subsp. enterica serovar Typhimurium 4:i:1,2_69M (S. Typhimurium 69M); it showed a very close similarity to S. enterica subsp. enterica serovar Typhimurium LT2. Bioinformatics sequence analysis detected eleven SPIs (SPI-1, SPI-2, SPI-3, SPI-4, SPI-5, SPI-9, SPI-12, SPI-13, SPI-14, C63PI, CS54_island). Significant changes in gene sequences were noted, causing frameshift mutations in yeiG, rfbP, fumA, yeaL, ybeU (insertion) and lpfD, avrA, ratB, yacH (deletion). The sequences of several proteins were significantly different from those coded in the reference genome; their three-dimensional structure was predicted and compared with reference proteins. Our findings indicate the presence of a number of antimicrobial resistance genes that do not directly imply an antibiotic resistance phenotype.

Elucidating the mechanism of antimicrobial resistance in Mycobacterium tuberculosis using gene interaction networks

Antimicrobial resistance (AMR) in microorganisms is an urgent global health threat. AMR of Mycobacterium tuberculosis is associated with significant morbidity and mortality. It is of great importance to underpin the resistance pathways involved in the mechanisms of AMR and identify the genes that are directly involved in AMR. The focus of the current study was the bacteria M. tuberculosis, which carries AMR genes that give resistance that lead to multidrug resistance. We, therefore, built a network of 43 genes and examined for potential gene-gene interactions. Then we performed a clustering analysis and identified three closely related clusters that could be involved in multidrug resistance mechanisms. Through the bioinformatics pipeline, we consistently identified six-hub genes (dnaN, polA, ftsZ, alr, ftsQ, and murC) that demonstrated the highest number of interactions within the clustering analysis. This study sheds light on the multidrug resistance of MTB and provides a protocol for discovering genes that might be involved in multidrug resistance, which will improve the treatment of resistant strains of TB.

Genome-Wide Association with Uncertainty in the Genetic Similarity Matrix

Genome-wide association studies (GWASs) are often confounded by population stratification and structure. Linear mixed models (LMMs) are a powerful class of methods for uncovering genetic effects, while controlling for such confounding. LMMs include random effects for a genetic similarity matrix, and they assume that a true genetic similarity matrix is known. However, uncertainty about the phylogenetic structure of a study population may degrade the quality of LMM results. This may happen in bacterial studies in which the number of samples or loci is small, or in studies with low-quality genotyping. In this study, we develop methods for linear mixed models in which the genetic similarity matrix is unknown and is derived from Markov chain Monte Carlo estimates of the phylogeny. We apply our model to a GWAS of multidrug resistance in tuberculosis, and illustrate our methods on simulated data.

The development of variation-based rifampicin resistance in Staphylococcus aureus deciphered through genomic and transcriptomic study

Rifampicin (RIF) resistance imposes a challenge on the antimicrobial treatment of pathogen infections. Figuring out the development mechanism of RIF resistance is critical to improving antimicrobial therapy strategy in clinics and biological treatment strategy of RIF polluted sewage in environmental engineering. The RIF resistance development of Staphylococcus aureus (S. aureus) with exposure to RIF at sub-inhibitory concentrations was comprehensively investigated via genomic and transcriptomic approaches in this study. RIF minimal inhibitory concentration (MIC) for S. aureus rapidly increased from 0.032 to 256 mg/L. Membrane permeability decrease, biofilm formation enhancement, and ROS production increase associated with RIF resistance were observed in RIF-induced strains. Through comparative genomic analysis, mutations in rpoB and rpoC were considered to be associated with RIF resistance in S. aureus mutants. Pan-genome-wide single-nucleotide variant analysis indicated that mutations at rpoB-1412, rpoB-1451, and rpoB-1457 were prevalent in 13849 public genomes of S. aureus, while mutations at rpoB-2256, and rpoC-3092 were first discovered in this study. The panorama of adaptative alteration of cellular physiological processes was observed via transcriptomic analysis. The oxidation pressure responses, metabolism, transporters, virulence factors, and multiple steps of DNA and RNA machinery were found to be perturbed by RIF in S. aureus.

Molecular mechanism of antibiotic resistance induced by mono- and twin-chained quaternary ammonium compounds

The usage of quaternary ammonium compounds (QACs) as disinfectants has increased dramatically since the outbreak of COVID-19 pandemic, leading to potentially accelerated emergence of antibiotic resistance. Long-term exposure to subinhibitory level QACs can lead to multidrug resistance, but the contribution of mutagenesis to resistance evolution is obscure. In this study, we subcultured E. coli K-12 under subinhibitory (0.25 × and 0.5 × Minimum Inhibitory Concentration, MIC) or inhibitory (1 × and 2 × MIC) concentrations of benzalkonium chloride (BAC, mono-chained) or didecyldimethylammonium chloride (DDAC, twin-chained) for 60 days. The sensitivity of QAC-adapted cells to five typical antibiotics decreased significantly, and in particular, the MIC of rifampicin increased by 85 times. E. coli adapted faster to BAC but developed 20-167% higher antibiotic resistance with 56% more mutations under DDAC exposure. The broader mutations induced by QACs, including negative regulators (acrR, marR, soxR, and crp), outer membrane proteins and transporters (mipA and sbmA), and RNA polymerase (rpoB and rpoC), potentially contributed to the high multi-drug resistance. After QACs stresses were removed, the phenotypic resistance induced by subinhibitory concentrations of QACs was reversible, whereas that induced by inhibitory concentrations of QACs was irreversible. The different patterns and molecular mechanism of antibiotic resistance induced by BAC and DDAC is informative to estimating the risks of broader QACs present at varied concentrations in the environment.

Detection of Isoniazid and Rifampin Resistance in Mycobacterium tuberculosis Clinical Isolates from Sputum Samples by High-Resolution Melting Analysis

The effective management of multidrug-resistant tuberculosis (MDR-TB) and the need for rapid and accurate screening of rifampin (RIF) and isoniazid (INH)-resistant Mycobacterium tuberculosis (Mtb) isolates are the most fundamental and difficult challenges facing the global TB control. The present study aimed to compare the diagnostic accuracy of high-resolution melting-curve analysis (HRMA) in comparison to multiplex allele-specific PCR (MAS-PCR) and xpert MTB/RIF as well as the conventional drug-susceptibility test (DST) and gene sequencing for the detection of INH and RIF resistance in the Mtb isolates. In the present study, a total of 431 Mtb isolates including 11 MDR (%2.55), 7 INH resistance (%1.62), two RIF resistance (%0.46), and 411 sensitive isolates were phenotypically confirmed. HRMA assay identified katG gene mutations and the mabA-inhA promoter region in 15 of 18 INH-resistant samples and rpoB gene mutations were successfully evaluated in 11 out of 13 RIF-resistant samples. The sensitivity and specificity of the HRMA method were 83.3% and 98.8% for INH and 84.6% and 99% for RIF, respectively. The most common mutation in RIF-resistance-determining region (RRDR) occurred at codon 531 (TCG → TTG)(84.6%) and then at codon 513 (CAA → GTA)(7.6%) and 526 (CAC → TAC) (7.6%), which resulted in the amino-acid changes. Also, 88.8% of INH-resistant samples had mutations in the katG gene and the mabA-inhA promoter region, of which the highest mutation occurred at codon 315 (AGC → ACC) of the katG gene. In conclusion, all these results indicated that the sensitivity and specificity of the HRM method were increased when the katG gene and the mabA-inhA promoter region were used as a target.

Whole genome sequencing of multidrug-resistant Mycobacterium tuberculosis isolates collected in the Czech Republic, 2005-2020

The emergence and spread of resistant tuberculosis (TB) pose a threat to public health, so it is necessary to diagnose the drug-resistant forms in a clinically short time frame and closely monitor their transmission. In this study, we carried out a first whole genome sequencing (WGS)-based analysis of multidrug resistant (MDR) M. tuberculosis strains to explore the phylogenetic lineages diversity, drug resistance mechanisms, and ongoing transmission chains within the country. In total, 65 isolates phenotypically resistant to at least rifampicin and isoniazid collected in the Czech Republic in 2005-2020 were enrolled for further analysis. The agreement of the results obtained by WGS with phenotypic drug susceptibility testing (pDST) in the determination of resistance to isoniazid, rifampicin, pyrazinamide, streptomycin, second-line injectables and fluoroquinolones was more than 80%. Phylogenetic analysis of WGS data revealed that the majority of MDR M. tuberculosis isolates were the Beijing lineage 2.2.1 (n = 46/65; 70.8%), while the remaining strains belonged to Euro-American lineage. Cluster analysis with a predefined cut-off distance of less than 12 single nucleotide polymorphisms between isolates showed 19 isolates in 6 clusters (clustering rate 29.2%), located mainly in the region of the capital city of Prague. This study highlights the utility of WGS as a high-resolution approach in the diagnosis, characterization of resistance patterns, and molecular-epidemiological analysis of resistant TB in the country.

Multidrug-resistant tuberculosis control in Rwanda overcomes a successful clone that causes most disease over a quarter century

SUMMARY BACKGROUND

Multidrug-resistant (MDR) tuberculosis (TB) poses an important challenge in TB management and control. Rifampicin resistance (RR) is a solid surrogate marker of MDR-TB. We investigated the RR-TB clustering rates, bacterial population dynamics to infer transmission dynamics, and the impact of changes to patient management on these dynamics over 27 years in Rwanda.

METHODS

We analysed whole genome sequences of a longitudinal collection of nationwide RR-TB isolates. The collection covered three important periods: before programmatic management of MDR-TB (PMDT; 1991-2005), the early PMDT phase (2006-2013), in which rifampicin drug-susceptibility testing (DST) was offered to retreatment patients only, and the consolidated phase (2014-2018), in which all bacteriologically confirmed TB patients had rifampicin DST done mostly via Xpert MTB/RIF assay. We constructed clusters based on a 5 SNP cut-off and resistance conferring SNPs. We used Bayesian modelling for dating and population size estimations, TransPhylo to estimate the number of secondary cases infected by each patient, and multivariable logistic regression to assess predictors of being infected by the dominant clone.

RESULTS

Of 308 baseline RR-TB isolates considered for transmission analysis, the clustering analysis grouped 259 (84.1%) isolates into 13 clusters. Within these clusters, a single dominant clone was discovered containing 213 isolates (82.2% of clustered and 69.1% of all RR-TB), which we named the "Rwanda Rifampicin-Resistant clone" (R3clone). R3clone isolates belonged to Ugandan sub-lineage 4.6.1.2 and its rifampicin and isoniazid resistance were conferred by the Ser450Leu mutation in rpoB and Ser315Thr in katG genes, respectively. All R3clone isolates had Pro481Thr, a putative compensatory mutation in the rpoC gene that likely restored its fitness. The R3clone was estimated to first arise in 1987 and its population size increased exponentially through the 1990s', reaching maximum size (∼84%) in early 2000 s', with a declining trend since 2014. Indeed, the highest proportion of R3clone (129/157; 82·2%, 95%CI: 75·3-87·8%) occurred between 2000 and 13, declining to 64·4% (95%CI: 55·1-73·0%) from 2014 onward. We showed that patients with R3clone detected after an unsuccessful category 2 treatment were more likely to generate secondary cases than patients with R3clone detected after an unsuccessful category 1 treatment regimen.

CONCLUSIONS

RR-TB in Rwanda is largely transmitted. Xpert MTB/RIF assay as first diagnostic test avoids unnecessary rounds of rifampicin-based TB treatment, thus preventing ongoing transmission of the dominant R3clone. As PMDT was intensified and all TB patients accessed rifampicin-resistance testing, the nationwide R3clone burden declined. To our knowledge, our findings provide the first evidence supporting the impact of universal DST on the transmission of RR-TB.

Transmission, distribution and drug resistance-conferring mutations of extensively drug-resistant tuberculosis in the Western Cape Province, South Africa

Extensively drug-resistant tuberculosis (XDR-TB), defined as resistance to at least isoniazid (INH), rifampicin (RIF), a fluoroquinolone (FQ) and a second-line injectable drug (SLID), is difficult to treat and poses a major threat to TB control. The transmission dynamics and distribution of XDR Mycobacterium tuberculosis (Mtb) strains have not been thoroughly investigated. Using whole genome sequencing data on 461 XDR-Mtb strains, we aimed to investigate the geographical distribution of XDR-Mtb strains in the Western Cape Province of South Africa over a 10 year period (2006–2017) and assess the association between Mtb sub-lineage, age, gender, geographical patient location and membership or size of XDR-TB clusters. First, we identified transmission clusters by excluding drug resistance-conferring mutations and using the 5 SNP cutoff, followed by merging clusters based on their most recent common ancestor. We then consecutively included variants conferring resistance to INH, RIF, ethambutol (EMB), pyrazinamide (PZA), SLIDs and FQs in the cluster definition. Cluster sizes were classified as small (2–4 isolates), medium (5–20 isolates), large (21–100 isolates) or very large (>100 isolates) to reflect the success of individual strains. We found that most XDR-TB strains were clustered and that including variants conferring resistance to INH, RIF, EMB, PZA and SLIDs in the cluster definition did not significantly reduce the proportion of clustered isolates (85.5–82.2 %) but increased the number of patients belonging to small clusters (4.3–12.4 %, P=0.56). Inclusion of FQ resistance-conferring variants had the greatest effect, with 11 clustered isolates reclassified as unique while the number of clusters increased from 17 to 37. Lineage 2 strains (lineage 2.2.1 typical Beijing or lineage 2.2.2 atypical Beijing) showed the large clusters which were spread across all health districts of the Western Cape Province. We identified a significant association between residence in the Cape Town metropole and cluster membership (P=0.016) but no association between gender, age and cluster membership or cluster size (P=0.39). Our data suggest that the XDR-TB epidemic in South Africa probably has its origin in the endemic spread of MDR Mtb and pre-XDR Mtb strains followed by acquisition of FQ resistance, with more limited transmission of XDR Mtb strains. This only became apparent with the inclusion of drug resistance-conferring variants in the definition of a cluster. In addition to the prevention of amplification of resistance, rapid diagnosis of MDR, pre-XDR and XDR-TB and timely initiation of appropriate treatment is needed to reduce transmission of difficult-to-treat TB.

Mycobacterium tuberculosis Affects Protein and Lipid Content of Circulating Exosomes in Infected Patients Depending on Tuberculosis Disease State

Tuberculosis (TB), which is caused by the bacterium Mycobacterium tuberculosis (Mtb), is still one of the deadliest infectious diseases. Understanding how the host and pathogen interact in active TB will have a significant impact on global TB control efforts. Exosomes are increasingly recognized as a means of cell-to-cell contact and exchange of soluble mediators. In the case of TB, exosomes are released from the bacillus and infected cells. In the present study, a comprehensive lipidomics and proteomics analysis of size exclusion chromatography-isolated plasma-derived exosomes from patients with TB lymphadenitis (TBL) and treated as well as untreated pulmonary TB (PTB) was performed to elucidate the possibility to utilize exosomes in diagnostics and knowledge building. According to our findings, exosome-derived lipids and proteins originate from both the host and Mtb in the plasma of active TB patients. Exosomes from all patients are mostly composed of sphingomyelins (SM), phosphatidylcholines, phosphatidylinositols, free fatty acids, triacylglycerols (TAG), and cholesterylesters. Relative proportions of, e.g., SMs and TAGs, vary depending on the disease or treatment state and could be linked to Mtb pathogenesis and dormancy. We identified three proteins of Mtb origin: DNA-directed RNA polymerase subunit beta (RpoC), Diacyglycerol O-acyltransferase (Rv2285), and Formate hydrogenase (HycE), the latter of which was discovered to be differently expressed in TBL patients. Furthermore, we discovered that Mtb infection alters the host protein composition of circulating exosomes, significantly affecting a total of 37 proteins. All TB patients had low levels of apolipoproteins, as well as the antibacterial proteins cathelicidin, Scavenger Receptor Cysteine Rich Family Member (SSC5D), and Ficolin 3 (FCN3). When compared to healthy controls, the protein profiles of PTB and TBL were substantially linked, with 14 proteins being co-regulated. However, adhesion proteins (integrins, Intercellular adhesion molecule 2 (ICAM2), CD151, Proteoglycan 4 (PRG4)) were shown to be more prevalent in PTB patients, while immunoglobulins, Complement component 1r (C1R), and Glutamate receptor-interacting protein 1 (GRIP1) were found to be more abundant in TBL patients, respectively. This study could confirm findings from previous reports and uncover novel molecular profiles not previously in focus of TB research. However, we applied a minimally invasive sampling and analysis of circulating exosomes in TB patients. Based on the findings given here, future studies into host–pathogen interactions could pave the way for the development of new vaccines and therapies.

Relationship between Resistance to Ethambutol and Rifampin and Clinical Outcomes in Mycobacterium avium Complex Pulmonary Disease

We evaluated the associations between the in vitro activities of ethambutol and rifampin and clinical outcomes of Mycobacterium avium complex (MAC) pulmonary disease (PD). Among 158 patients with MAC-PD, there was no relationship between high MICs for ethambutol and/or rifampin and treatment failure for MAC-PD. Ethambutol and rifampin resistance was common among MAC isolates (rates of 87% and 59%, respectively), but mutations in embB, rpoB, and rpoC were rare, with detection in only 4% of the drug-resistant MAC isolates.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

Heterogeneous fitness landscape cues, pknG low expression, and phthiocerol dimycocerosate low production of Mycobacterium tuberculosis ATCC25618 rpoB S450L in enriched broth

Multidrug-resistant tuberculosis (isoniazid/rifampin[RIF]-resistant TB) ravages developing countries. Fitness is critical in clinical outcomes. Previous studies on RIF-resistant TB (RR-TB) showed competitive fitness gains and losses, with rpoB-S450L as the most isolated/fit mutation. This study measured virulence/resistance genes, phthiocerol dimycocerosate (PDIM) levels and their relationship with rpoB S450L ATCC25618 RR-TB strain fitness. After obtaining 10 different RR-TB GenoType MTBDRplus 2.0-genotyped isolates (with nontyped, S441, H445 and S450 positions), only one S450L isolate (R9, rpoB-S450L ATCC 25618, RR 1 μg/mL) was observed, with H445Y being the most common. A competitive fitness in vitro assay with wild-type (wt) ATCC 25618: R9 1:1 in 50 mL Middlebrook 7H9/OADC was performed, and generation time (G) in vitro and relative fitness were obtained. mRNA and PDIM were extracted on log and stationary phases. Fitness decreased in rpoB S450L and H445Y strains, with heterogeneous fitness cues in three biological replicas of rpoB-S450L: one high and two low fitness replicas. S450L strain had significant pknG increase. Compared with S450L, wt-rpoB showed increased polyketide synthase ppsA expression and high PDIM peak measured by HPLC-MS in log phase compared to S450L. This contrasts with previously increased PDIM in other RR-TB isolates.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Expression Dysregulation as a Mediator of Fitness Costs in Antibiotic Resistance

Antimicrobial resistance (AMR) poses a threat to global health and the economy. Rifampicin-resistant Mycobacterium tuberculosis accounts for a third of the global AMR burden. Gaining the upper hand on AMR requires a deeper understanding of the physiology of resistance. AMR often results in a fitness cost in the absence of drug. Identifying the molecular mechanisms underpinning this cost could help strengthen future treatment regimens. Here, we used a collection of M. tuberculosis strains that provide an evolutionary and phylogenetic snapshot of rifampicin resistance and subjected them to genome-wide transcriptomic and proteomic profiling to identify key perturbations of normal physiology. We found that the clinically most common rifampicin resistance-conferring mutation, RpoB Ser450Leu, imparts considerable gene expression changes, many of which are mitigated by the compensatory mutation in RpoC Leu516Pro. However, our data also provide evidence for pervasive epistasis—the same resistance mutation imposed a different fitness cost and functionally distinct changes to gene expression in genetically unrelated clinical strains. Finally, we report a likely posttranscriptional modulation of gene expression that is shared in most of the tested strains carrying RpoB Ser450Leu, resulting in an increased abundance of proteins involved in central carbon metabolism. These changes contribute to a more general trend in which the disruption of the composition of the proteome correlates with the fitness cost of the RpoB Ser450Leu mutation in different strains.

Compensatory effects of M. tuberculosis rpoB mutations outside the rifampicin resistance-determining region

Mycobacterium tuberculosis has been observed to develop resistance to the frontline anti-tuberculosis drug rifampicin, primarily through mutations in the rifampicin resistance-determining region (RRDR) of rpoB. While these mutations have been determined to confer a fitness cost, compensatory mutations in rpoA and rpoC that may enhance the fitness of resistant strains have been demonstrated. Recent genomic studies identified several rpoB non-RRDR mutations that co-occurred with RRDR mutations in clinical isolates without rpoA/rpoC mutations and may confer fitness compensation. In this study, we identified 33 evolutionarily convergent rpoB non-RRDR mutations through phylogenomic analysis of public genomic data for clinical M. tuberculosis isolates. We found that none of these mutations, except V170F and I491F, can cause rifampin resistance in Mycolicibacterium smegmatis. The compensatory effects of five representative mutations across rpoB were evaluated by an in vitro competition assay, through which we observed that each of these mutations can significantly improve the relative fitness of the initial S450L mutant (0.97–1.08 vs 0.87). Furthermore, we observed that the decreased RNAP transcription efficiency introduced by S450L was significantly alleviated by each of the five mutations. Structural analysis indicated that the fitness compensation observed for the non-RRDR mutations might be achieved by modification of the RpoB active centre or by changes in interactions between RNAP subunits. Our results provide experimental evidence supporting that compensatory effects are exerted by several rpoB non-RRDR mutations, which could be utilized as additional molecular markers for predicting the fitness of clinical rifampin-resistant M. tuberculosis strains.

Inhibition of RNA Polymerase by Rifampicin and Rifamycin-Like Molecules

RNA polymerases (RNAPs) accomplish the first step of gene expression in all living organisms. However, the sequence divergence between bacterial and human RNAPs makes the bacterial RNAP a promising target for antibiotic development. The most clinically important and extensively studied class of antibiotics known to inhibit bacterial RNAP are the rifamycins. For example, rifamycins are a vital element of the current combination therapy for treatment of tuberculosis. Here, we provide an overview of the history of the discovery of rifamycins, their mechanisms of action, the mechanisms of bacterial resistance against them, and progress in their further development.

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.

Characterization of rifampicin-resistant Mycobacterium tuberculosis in Khyber Pakhtunkhwa, Pakistan

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is endemic in Pakistan. Resistance to both firstline rifampicin and isoniazid drugs (multidrug-resistant TB; MDR-TB) is hampering disease control. Rifampicin resistance is attributed to rpoB gene mutations, but rpoA and rpoC loci may also be involved. To characterise underlying rifampicin resistance mutations in the TB endemic province of Khyber Pakhtunkhwa, we sequenced 51 M. tuberculosis isolates collected between 2016 and 2019; predominantly, MDR-TB (n = 44; 86.3%) and lineage 3 (n = 30, 58.8%) strains. We found that known mutations in rpoB (e.g. S405L), katG (e.g. S315T), or inhA promoter loci explain the MDR-TB. There were 24 unique mutations in rpoA, rpoB, and rpoC genes, including four previously unreported. Five instances of within-host resistance diversity were observed, where two were a mixture of MDR-TB strains containing mutations in rpoB, katG, and the inhA promoter region, as well as compensatory mutations in rpoC. Heteroresistance was observed in two isolates with a single lineage. Such complexity may reflect the high transmission nature of the Khyber Pakhtunkhwa setting. Our study reinforces the need to apply sequencing approaches to capture the full-extent of MDR-TB genetic diversity, to understand transmission, and to inform TB control activities in the highly endemic setting of Pakistan.

Mutant RNA polymerase can reduce susceptibility to antibiotics via ppGpp-independent induction of a stringent-like response

Background

Mutations in RNA polymerase (RNAP) can reduce susceptibility to ciprofloxacin in Escherichia coli, but the mechanism of transcriptional reprogramming responsible is unknown. Strains carrying ciprofloxacin-resistant (Cip^R) rpoB mutations have reduced growth fitness and their impact on clinical resistance development is unclear.

Objectives

To assess the potential for Cip^RrpoB mutations to contribute to resistance development by estimating the number of distinct alleles. To identify fitness-compensatory mutations that ameliorate the fitness costs of Cip^RrpoB mutations. To understand how Cip^RrpoB mutations reprogramme RNAP.

Methods

E. coli strains carrying five different Cip^RrpoB alleles were evolved with selection for improved fitness and characterized for acquired mutations, relative fitness and MIC_Cip. The effects of dksA mutations and a ppGpp⁰ background on growth and susceptibility phenotypes associated with Cip^RrpoB alleles were determined.

Results

The number of distinct Cip^RrpoB mutations was estimated to be >100. Mutations in RNAP genes and in dksA can compensate for the fitness cost of Cip^RrpoB mutations. Deletion of dksA reduced the MIC_Cip for strains carrying Cip^RrpoB alleles. A ppGpp⁰ phenotype had no effect on drug susceptibility.

Conclusions

Cip^RrpoB mutations induce an ppGpp-independent stringent-like response. Approximately half of the reduction in ciprofloxacin susceptibility is caused by an increased affinity of RNAP to DksA while the other half is independent of DksA. Stringent-like response activating mutations might be the most diverse class of mutations reducing susceptibility to antibiotics.

Prisons as ecological drivers of fitness-compensated multidrug-resistant Mycobacterium tuberculosis

Multidrug-resistant tuberculosis (MDR-TB) accounts for one third of the annual deaths due to antimicrobial resistance¹. Drug resistance-conferring mutations frequently cause fitness costs in bacteria^2-5. Experimental work indicates that these drug resistance-related fitness costs might be mitigated by compensatory mutations^6-10. However, the clinical relevance of compensatory evolution remains poorly understood. Here we show that, in the country of Georgia, during a 6-year nationwide study, 63% of MDR-TB was due to patient-to-patient transmission. Compensatory mutations and patient incarceration were independently associated with transmission. Furthermore, compensatory mutations were overrepresented among isolates from incarcerated individuals that also frequently spilled over into the non-incarcerated population. As a result, up to 31% of MDR-TB in Georgia was directly or indirectly linked to prisons. We conclude that prisons fuel the epidemic of MDR-TB in Georgia by acting as ecological drivers of fitness-compensated strains with high transmission potential.

Generation of Tetracycline and Rifamycin Resistant Chlamydia Suis Recombinants

The Chlamydiaceae are a family of obligate intracellular, gram-negative bacteria known to readily exchange DNA by homologous recombination upon co-culture in vitro, allowing the transfer of antibiotic resistance residing on the chlamydial chromosome. Among all the obligate intracellular bacteria, only Chlamydia (C.) suis naturally integrated a tetracycline resistance gene into its chromosome. Therefore, in order to further investigate the readiness of Chlamydia to exchange DNA and especially antibiotic resistance, C. suis is an excellent model to advance existing co-culture protocols allowing the identification of factors crucial to promote homologous recombination in vitro. With this strategy, we co-cultured tetracycline-resistant with rifamycin group-resistant C. suis, which resulted in an allover recombination efficiency of 28%. We found that simultaneous selection is crucial to increase the number of recombinants, that sub-inhibitory concentrations of tetracycline inhibit rather than promote the selection of double-resistant recombinants, and identified a recombination-deficient C. suis field isolate, strain SWA-110 (1-28b). While tetracycline resistance was detected in field isolates, rifampicin/rifamycin resistance (RifR) had to be induced in vitro. Here, we describe the protocol with which RifR C. suis strains were generated and confirmed. Subsequent whole-genome sequencing then revealed that G530E and D461A mutations in rpoB, a gene encoding for the β-subunit of the bacterial RNA polymerase (RNAP), was likely responsible for rifampicin and rifamycin resistance, respectively. Finally, whole-genome sequencing of recombinants obtained by co-culture revealed that recombinants picked from the same plate may be sibling clones and confirmed C. suis genome plasticity by revealing variable, apparently non-specific areas of recombination.

Development of New Therapeutics to Meet the Current Challenge of Drug Resistant Tuberculosis

Tuberculosis (TB) is a prominent infective disease and a major reason of mortality/ morbidity globally. Mycobacterium tuberculosis causes a long-lasting latent infection in a significant proportion of human population. The increasing burden of tuberculosis is mainly caused due to multi drug-resistance. The failure of conventional treatment has been observed in large number of cases. Drugs that are used to treat extensively drug-resistant tuberculosis are expensive, have limited efficacy, and have more side effects for a longer duration of time and are often associated with poor prognosis. To regulate the emergence of multidrug resistant tuberculosis, extensively drug-resistant tuberculosis and totally drug resistant tuberculosis, efforts are being made to understand the genetic/molecular basis of target drug delivery and mechanisms of drug resistance. Understanding the molecular approaches and pathology of Mycobacterium tuberculosis through whole genome sequencing may further help in the improvement of new therapeutics to meet the current challenge of global health. Understanding cellular mechanisms that trigger resistance to Mycobacterium tuberculosis infection may expose immune associates of protection, which could be an important way for vaccine development, diagnostics, and novel host-directed therapeutic strategies. The recent development of new drugs and combinational therapies for drug-resistant tuberculosis through major collaboration between industry, donors, and academia gives an improved hope to overcome the challenges in tuberculosis treatment. In this review article, an attempt was made to highlight the new developments of drug resistance to the conventional drugs and the recent progress in the development of new therapeutics for the treatment of drugresistant and non-resistant cases.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Drug resistance, fitness and compensatory mutations in Mycobacterium tuberculosis

For tuberculosis to be eradicated, the transmission of Multi-Drug-Resistant and eXtensively Drug Resistant strains of Mycobacterium tuberculosis (MDR and XDR-TB) must be considerably reduced. Drug resistant strains were initially thought to have reduced fitness, and the majority of resistant strains may actually have compromised fitness because they are found in only one or a few patients. In contrast, some MDR/XDR-TB strains are highly transmitted and cause large outbreaks. Most antibiotics target essential bacterial functions and the mutations that confer resistance to anti-TB drugs can incur fitness costs manifested as slower growth and reduced viability. The fitness costs vary with different resistance mutations and the bacilli can also accumulate secondary mutations that compensate for the compromised functions and partially or fully restore lost fitness. The compensatory mutations (CM) are different for each antibiotic, as they mitigate the deleterious effects of the specific functions compromised by the resistance mutations. CM are generally more common in strains with resistance mutations incurring the greatest fitness costs, but for RIF resistance, CM are most frequent in strains with the mutation carrying the least fitness cost, Ser450Leu. Here, we review what is known about fitness costs, CM and mechanisms of resistance to the drugs that define a strain as MDR or XDR-TB. The relative fitness costs of the resistance mutations and the mitigating effects of CM largely explain why certain mutations are frequently found in highly transmitted clusters while others are less frequently, rarely or never found in clinical isolates. The CM illustrate how drug resistance affects bacteria and how bacteria evolve to overcome the effects of the antibiotics, and thus a paradigm for how mycobacteria can evolve in response to stress.

A Bioinformatics Whole-Genome Sequencing Workflow for Clinical Mycobacterium tuberculosis Complex Isolate Analysis, Validated Using a Reference Collection Extensively Characterized with Conventional Methods and In Silico Approaches

The use of whole-genome sequencing (WGS) for routine typing of bacterial isolates has increased substantially in recent years. For Mycobacterium tuberculosis (MTB), in particular, WGS has the benefit of drastically reducing the time required to generate results compared to most conventional phenotypic methods. Consequently, a multitude of solutions for analyzing WGS MTB data have been developed, but their successful integration in clinical and national reference laboratories is hindered by the requirement for their validation, for which a consensus framework is still largely absent. We developed a bioinformatics workflow for (Illumina) WGS-based routine typing of MTB complex (MTBC) member isolates allowing complete characterization, including (sub)species confirmation and identification (16S, csb/RD, hsp65), single nucleotide polymorphism (SNP)-based antimicrobial resistance (AMR) prediction, and pathogen typing (spoligotyping, SNP barcoding, and core genome multilocus sequence typing). Workflow performance was validated on a per-assay basis using a collection of 238 in-house-sequenced MTBC isolates, extensively characterized with conventional molecular biology-based approaches supplemented with public data. For SNP-based AMR prediction, results from molecular genotyping methods were supplemented with in silico modified data sets, allowing us to greatly increase the set of evaluated mutations. The workflow demonstrated very high performance with performance metrics of >99% for all assays, except for spoligotyping, where sensitivity dropped to ∼90%. The validation framework for our WGS-based bioinformatics workflow can aid in the standardization of bioinformatics tools by the MTB community and other SNP-based applications regardless of the targeted pathogen(s). The bioinformatics workflow is available for academic and nonprofit use through the Galaxy instance of our institute at https://galaxy.sciensano.be.

The DnaK Chaperone System Buffers the Fitness Cost of Antibiotic Resistance Mutations in Mycobacteria

Chaperones aid in protein folding and maintenance of protein integrity. In doing so, they have the unique ability to directly stabilize resistance-conferring amino acid substitutions in drug targets and to counter the stress imparted by these substitutions, thus supporting heritable antimicrobial resistance (AMR). We asked whether chaperones support AMR in Mycobacterium smegmatis, a saprophytic model of Mycobacterium tuberculosis, the causative agent of tuberculosis (TB). We show that DnaK associates with many drug targets and that DnaK associates more with AMR-conferring mutant RNA polymerase (RNAP) than with wild-type RNAP. In addition, frequency-of-resistance (FOR) and fitness studies reveal that the DnaK system of chaperones supports AMR in antimicrobial targets in mycobacteria, including RNAP and the ribosome. These findings highlight chaperones as potential targets for drugs to overcome AMR in mycobacteria, including M. tuberculosis, as well as in other pathogens.IMPORTANCE AMR is a global problem, especially for TB. Here, we show that mycobacterial chaperones support AMR in M. smegmatis, a nonpathogenic model of M. tuberculosis, the causative agent of TB. In particular, the mycobacterial DnaK system of chaperones supports AMR in the antimicrobial targets RNA polymerase and the ribosome. This is the first report showing a role for protein chaperones in mediating AMR in mycobacteria. Given the widespread role of protein chaperones in enabling genomic diversity, we anticipate that our findings can be extended to other microbes.

Low-Level Rifampin Resistance and rpoB Mutations in Mycobacterium tuberculosis: an Analysis of Whole-Genome Sequencing and Drug Susceptibility Test Data in New York

Rapid and reliable detection of rifampin (RIF) resistance is critical for the diagnosis and treatment of drug-resistant and multidrug-resistant (MDR) tuberculosis. Discordant RIF phenotype/genotype susceptibility results remain a challenge due to the presence of rpoB mutations that do not confer high levels of RIF resistance, as have been exhibited in strains with mutations such as Ser450Leu. These strains, termed low-level RIF resistant, exhibit elevated RIF MICs compared to fully susceptible strains but remain phenotypically susceptible by mycobacterial growth indicator tube (MGIT) testing and have been associated with poor patient outcomes. Here, we assess RIF resistance prediction by whole-genome sequencing (WGS) among a set of 1,779 prospectively tested strains by both prevalence of rpoB gene mutation and phenotype as part of routine clinical testing during a 2.5-year period. During this time, 139 strains were found to have nonsynonymous rpoB mutations, 53 of which were associated with RIF resistance, including both low-level and high-level resistance. Resistance to RIF (1.0 μg/ml in MGIT) was identified in 43 (81.1%) isolates. The remaining 10 (18.9%) strains were susceptible by MGIT but were confirmed to be low-level RIF resistant by MIC testing. Full rpoB gene sequencing overcame the limitations of critical concentration phenotyping, probe-based genotyping, and partial gene sequencing methods. Universal clinical WGS with concurrent phenotypic testing provided a more complete understanding of the prevalence and type of rpoB mutations and their association with RIF resistance in New York.

Comprehensive review on mechanism of action, resistance and evolution of antimycobacterial drugs

Tuberculosis is one of the deadliest infectious diseases existing in the world since ancient times and still possesses serious threat across the globe. Each year the number of cases increases due to high drug resistance shown by Mycobacterium tuberculosis (Mtb). Available antimycobacterial drugs have been classified as First line, Second line and Third line antibiotics depending on the time of their discoveries and their effectiveness in the treatment. These antibiotics have a broad range of targets ranging from cell wall to metabolic processes and their non-judicious and uncontrolled usage in the treatment for years has created a significant problem called multi-drug resistant (MDR) tuberculosis. In this review, we have summarized the mechanism of action of all the classified antibiotics currently in use along with the resistance mechanisms acquired by Mtb. We have focused on the new drug candidates/repurposed drugs, and drug in combinations, which are in clinical trials for either treating the MDR tuberculosis more effectively or involved in reducing the time required for the chemotherapy of drug sensitive TB. This information is not discussed very adequately on a single platform. Additionally, we have discussed the recent technologies that are being used to discover novel resistance mechanisms acquired by Mtb and for exploring novel drugs. The story of intrinsic resistance mechanisms and evolution in Mtb is far from complete. Therefore, we have also discussed intrinsic resistance mechanisms of Mtb and their evolution with time, emphasizing the hope for the development of novel antimycobacterial drugs for effective therapy of tuberculosis.

Evolution of Drug-Resistant Mycobacterium tuberculosis Strains and Their Adaptation to the Human Lung Environment

In the last two decades, multi (MDR), extensively (XDR), extremely (XXDR) and total (TDR) drug-resistant Mycobacterium tuberculosis (M.tb) strains have emerged as a threat to public health worldwide, stressing the need to develop new tuberculosis (TB) prevention and treatment strategies. It is estimated that in the next 35 years, drug-resistant TB will kill around 75 million people and cost the global economy $16.7 trillion. Indeed, the COVID-19 pandemic alone may contribute with the development of 6.3 million new TB cases due to lack of resources and enforced confinement in TB endemic areas. Evolution of drug-resistant M.tb depends on numerous factors, such as bacterial fitness, strain's genetic background and its capacity to adapt to the surrounding environment, as well as host-specific and environmental factors. Whole-genome transcriptomics and genome-wide association studies in recent years have shed some insights into the complexity of M.tb drug resistance and have provided a better understanding of its underlying molecular mechanisms. In this review, we will discuss M.tb phenotypic and genotypic changes driving resistance, including changes in cell envelope components, as well as recently described intrinsic and extrinsic factors promoting resistance emergence and transmission. We will further explore how drug-resistant M.tb adapts differently than drug-susceptible strains to the lung environment at the cellular level, modulating M.tb-host interactions and disease outcome, and novel next generation sequencing (NGS) strategies to study drug-resistant TB.

Early Drug Development and Evaluation of Putative Antitubercular Compounds in the -Omics Era

Tuberculosis (TB) is an infectious disease caused by the bacterium Mycobacterium tuberculosis. According to the WHO, the disease is one of the top 10 causes of death of people worldwide. Mycobacterium tuberculosis is an intracellular pathogen with an unusually thick, waxy cell wall and a complex life cycle. These factors, combined with M. tuberculosis ability to enter prolonged periods of latency, make the bacterium very difficult to eradicate. The standard treatment of TB requires 6-20months, depending on the drug susceptibility of the infecting strain. The need to take cocktails of antibiotics to treat tuberculosis effectively and the emergence of drug-resistant strains prompts the need to search for new antitubercular compounds. This review provides a perspective on how modern -omic technologies facilitate the drug discovery process for tuberculosis treatment. We discuss how methods of DNA and RNA sequencing, proteomics, and genetic manipulation of organisms increase our understanding of mechanisms of action of antibiotics and allow the evaluation of drugs. We explore the utility of mathematical modeling and modern computational analysis for the drug discovery process. Finally, we summarize how -omic technologies contribute to our understanding of the emergence of drug resistance.

Abattoir survey of bovine tuberculosis in tanta, centre of the Nile delta, with in silico analysis of gene mutations and protein-protein interactions of the involved mycobacteria

Bovine tuberculosis is a transboundary disease of high economic and public health burden worldwide. In this study, post-mortem examination of 750 cattle and buffalo in Tanta abattoir, Centre of the Nile Delta, revealed visible TB in 4% of animals and a true prevalence of 6.85% (95% CI: 5.3%-8.9%). Mycobacterial culture, histopathology and RT-PCR targeting all members of M. tuberculosis complex were performed, upon which 85%, 80% and 100% of each tested lesions were confirmed as TB, respectively. Mpb70-targeting PCR was conducted on ten RT-PCR positive samples for sequencing and identified nine Mycobacterium (M.) bovis strains and, interestingly, one M. tuberculosis (Mtb) strain from a buffalo. Bioinformatics tools were used for prediction of mutations, nucleotide polymorphisms, lineages, drug resistance and protein-protein interactions (PPI) of the sequenced strains. The Mtb strain was resistant to rifampicin, isoniazid and streptomycin, and to the best of our knowledge, this is the first report of multidrug resistant (MDR)-Mtb originating from buffaloes. Seven M. bovis strains were resistant to ethambutol and ethionamide. Such resistances were associated with KatG, rpoB, rpsL, embB and ethA genes mutations. Other mutations and nucleotide polymorphisms were also predicted, some are reported for the first time and require experimental work for validation. PPI revealed more interactions than what would be expected for a random set of proteins of similar size and had dense interactions between nodes that are biologically connected, as a group. Two M. bovis strains belonged to BOV AFRI lineage (Spoligotypes BOV 1; BOV 2) and eight strains belonged to East-Asian (Beijing) lineage. In conclusion, visible TB was prevalent in the study area, RT-PCR is the best to confirm the disease, MDR-Mtb is associated with buffalo TB, and mycobacteria of different lineages carry many resistance genes to chemotherapeutic agents used in treatment of human TB constituting a major public health risk.

Clinical Interpretation of Drug Susceptibility Tests in Tuberculosis

:

Prompt and accurate diagnosis of drug resistance is essential for optimal treatment of drug-resistant tuberculosis. However, only 20% of the more than half a million patients eligible for the treatment of MDR-TB/RR-TB receive an appropriate drug regimen. Drug-resistant TB regimens must include a sufficient number of effective medications, a significant challenge for clinicians worldwide, as most are forced to make therapeutic decisions without any, or very little information on drug susceptibility testing. Although phenotypic DST is still commonly regarded as the gold standard for determining M. tuberculosis susceptibility to antituberculosis drugs, it has several limitations, mainly its prolonged turnaround time. Molecular methods based on M. tuberculosis genomic DNA sequencing have been developed during the past two decades, to identify the most common mutations involved in drug resistance. The Xpert ® MTB/RIF is a real-time polymerase chain reaction that offers results in less than two hours and has an overall sensitivity for rifampin resistance of 96% and 98% specificity. Line probe assays (LPAs) are commercial DNA strip-based tests for detecting the most frequent mutations responsible for resistance to rifampin, isoniazid, fluoroquinolones, and second-line injectable drugs.

:

Discrepancies between phenotypic and genotyping methods are a problem that the clinician will face in everyday practice. However, any resistance result (with any type of test) in a person with risk factors for harboring resistant microorganisms should be considered appropriate while the results of complementary tests are available.

Role of post-translational modifications in the acquisition of drug resistance in Mycobacterium tuberculosis

Tuberculosis (TB) is one of the primary causes of deaths due to infectious diseases. The current TB regimen is long and complex, failing of which leads to relapse and/or the emergence of drug resistance. There is a critical need to understand the mechanisms of resistance development. With increasing drug pressure, Mycobacterium tuberculosis (Mtb) activates various pathways to counter drug-related toxicity. Signaling modules steer the evolution of Mtb to a variant that can survive, persist, adapt, and emerge as a form that is resistant to one or more drugs. Recent studies reveal that about 1/3rd of the annotated Mtb proteome is modified post-translationally, with a large number of these proteins being essential for mycobacterial survival. Post-translational modifications (PTMs) such as phosphorylation, acetylation, and pupylation play a salient role in mycobacterial virulence, pathogenesis, and metabolism. The role of many other PTMs is still emerging. Understanding the signaling pathways and PTMs may assist clinical strategies and drug development for Mtb. In this review, we explore the contribution of PTMs to mycobacterial physiology, describe the related cellular processes, and discuss how these processes are linked to drug resistance. A significant number of drug targets, InhA, RpoB, EmbR, and KatG, are modified at multiple residues via PTMs. A better understanding of drug-resistance regulons and associated PTMs will aid in developing effective drugs against TB.

Characterization of Mycobacterium tuberculosis genotypes and their correlation to multidrug resistance in Lusaka, Zambia

OBJECTIVES

The burden of multidrug-resistant tuberculosis (MDR-TB) has been reported to be increasing in Zambia. The reasons for the increase are still unclear. This study determined the diversity of Mycobacterium tuberculosis genotypes among isolates in Lusaka, the capital city, and investigated their association with MDR-TB.

METHODS

Spoligotyping, large sequence polymorphism (LSP) analysis, and sequencing of MDR associated genes were performed on a total of 274 M. tuberculosis clinical isolates stored at the University Teaching Hospital from 2013 to 2017. Of these, 134 were MDR-TB while 126 were pan-susceptible.

RESULTS

Spoligotyping showed the LAM family as the most predominant genotype (149/274, 54.4%) followed by the CAS family (44/274, 16.1%), T family (39/274, 14.2%), and minor proportions of X, S, Harleem, EAI and Beijing spoligofamilies were identified. Three M. bovis isolates were also observed. Among those, CAS1-Kili (SIT 21) and LAM1 (SIT 20) subfamilies showed a propensity for MDR-TB with p = 0.0001 and p = 0.001, respectively.

CONCLUSIONS

This phenomenon might explain the future increase in the MDR-TB burden caused by specific lineages in Zambia. Therefore, it is recommended that the National TB control program in the country complements conventional control strategies with molecular analysis for monitoring and surveillance of MDR-TB epidemiology.

HIV Coinfection Is Associated with Low-Fitness rpoB Variants in Rifampicin-Resistant Mycobacterium tuberculosis

We analyzed 312 drug-resistant genomes of Mycobacterium tuberculosis isolates collected from HIV-coinfected and HIV-negative TB patients from nine countries with a high tuberculosis burden. We found that rifampicin-resistant M. tuberculosis strains isolated from HIV-coinfected patients carried disproportionally more resistance-conferring mutations in rpoB that are associated with a low fitness in the absence of the drug, suggesting these low-fitness rpoB variants can thrive in the context of reduced host immunity.