Fitness-compensatory mutations in rifampicin-resistant RNA polymerase

Strain-specific Growth Parameters are Important to Accurately Model Bacterial Growth on Baby Spinach in Simulation Models

Digital tools to predict produce shelf life have the potential to reduce food waste and improve consumer satisfaction. To address this need, we (i) performed an observational study on the microbial quality of baby spinach, (ii) completed growth experiments of bacteria that are representative of the baby spinach microbiota, and (iii) developed an initial simulation model of bacterial growth on baby spinach. Our observational data showed that the predominant genera found on baby spinach were Pseudomonas, Pantoea and Exiguobacterium. Rifampicin-resistant mutants (rifR mutants) of representative bacterial subtypes were subsequently generated to obtain strain-specific growth parameters on baby spinach. These experiments showed that: (i) it is difficult to select rifR mutants that do not have fitness costs affecting growth (9 of 15 rifR mutants showed substantial differences in growth, compared to their corresponding wild-type strain) and (ii) based on estimates from primary growth models, the mean (geometric) maximum population of rifR mutants on baby spinach (7.6 log10 CFU/g, at 6°C) appears lower than that of the spinach microbiota (9.6 log10 CFU/g, at 6°C), even if rifR mutants did not have substantial growth-related fitness costs. Thus, a simulation model, parameterized with the data obtained here as well as literature data on home refrigeration temperatures, underestimated bacterial growth on baby spinach. The root mean square error of the simulation's output, compared against data from the observational study, was 1.11 log10 CFU/g. Sensitivity analysis was used to identify key parameters (e.g., strain maximum population) that impact the simulation model's output, allowing for prioritization of future data collection to improve the simulation model. Overall, this study provides a roadmap for the development of models to predict bacterial growth on leafy vegetables with strain-specific parameters and suggests that additional data are required to improve these models.

Persistent Colonization of Ciprofloxacin-Resistant and Extended-Spectrum β-Lactamase (ESBL)-Producing Salmonella enterica Serovar Kentucky ST198 in a Patient with Inflammatory Bowel Disease

Objective

Salmonella enterica serovar Kentucky ST198 has emerged as a global threat to humans. In this study, we aimed to characterize the prolonged carriage of ciprofloxacin-resistant and extended-spectrum β-lactamase (ESBL)-producing S. Kentucky ST198 in a single patient with inflammatory bowel disease (IBD).

Methods

Three S. Kentucky strains were collected from a single patient with IBD on 11th January, 23rd January, and 8th February, 2022, respectively. Antimicrobial susceptibility testing, whole-genome sequencing, and phylogenetic analysis with 38 previously described Chinese S. Kentucky ST198 strains from patients and food were performed.

Results

All three S. Kentucky isolates belonged to ST198. They carried identical 16 resistance genes, such as blaCTX-M-55, tet(A), and qnrS1, and had identical mutations within gyrA (S83F and D87N) and parC (S80I). Therefore, they exhibited identical multidrug-resistant profiles, including the clinically important antibiotics cephalosporins (ceftazidime and cefepime), fluoroquinolones (ciprofloxacin and levofloxacin), and third-generation tetracycline (tigecycline). Our three S. Kentucky strains were classified into the subclade ST198.2-2, and were genetically identical (2-6 SNPs) to each other. They exhibited a close genetic similarity (15-20 SNPs) to the isolate NT-h3189 from a patient and AH19MCS1 from chicken meat in China, indicating a possible epidemiological link between these S. Kentucky ST198 isolates from the patients and chicken meat.

Conclusion

Long-term colonization of ciprofloxacin-resistant and ESBL-producing S. Kentucky ST198 in a single patient is a matter of concern. Due to the potential transfer of S. Kentucky ST198 from food sources to humans, ongoing surveillance of this particular clone in animals, animal-derived food products, and humans should be strengthened.

Bacteria can compensate the fitness costs of amplified resistance genes via a bypass mechanism

Antibiotic heteroresistance is a phenotype in which a susceptible bacterial population includes a small subpopulation of cells that are more resistant than the main population. Such resistance can arise by tandem amplification of DNA regions containing resistance genes that in single copy are not sufficient to confer resistance. However, tandem amplifications often carry fitness costs, manifested as reduced growth rates. Here, we investigated if and how these fitness costs can be genetically ameliorated. We evolved four clinical isolates of three bacterial species that show heteroresistance to tobramycin, gentamicin and tetracyclines at increasing antibiotic concentrations above the minimal inhibitory concentration (MIC) of the main susceptible population. This led to a rapid enrichment of resistant cells with up to an 80-fold increase in the resistance gene copy number, an increased MIC, and severely reduced growth rates. When further evolved in the presence of antibiotic, these strains acquired compensatory resistance mutations and showed a reduction in copy number while maintaining high-level resistance. A deterministic model indicated that the loss of amplified units was driven mainly by their fitness costs and that the compensatory mutations did not affect the loss rate of the gene amplifications. Our findings suggest that heteroresistance mediated by copy number changes can facilitate and precede the evolution towards stable resistance.

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.

Adaptive evolution in asymptomatic host confers MDR Salmonella with enhanced environmental persistence and virulence

As a common cause for food-borne diseases, the Salmonella spp. are generally prevalent among livestock, whereby they are likely to be transmitted to human via environmental contamination. To explore the potential mechanism for prevalence of MDR Salmonella and its risk for dissemination via contaminated environments, we profiled the colonization dynamics of MDR Salmonella in chicken, herein we found that an adaptive evolution, driven by mutagenesis in a small protein-encoding gene (STM14_1829), conferred the multidrug resistant (MDR) Salmonella with increased fitness in asymptomatic host. Then the mechanistic study demonstrated that only one amino acid substitution in small protein STM14_1829 rendered MDR Salmonella capable to better invade and persist in phagocytotic cells by modulating bacterial flagella overexpression. Concerningly, the evolved Salmonella was also more resilient to the potential stressors generally found in environments and food processing, including heat, cold, adverse pH and oxidations. It implied that the evolved subpopulations are plausibly more persistent in environments once they contaminated through animal manure or human excreta. Moreover, the evolution promoted the pathogenesis caused by MDR Salmonella in susceptible hosts, resulting in higher risk for dissemination of pathogens via contaminated environments. Together, our data provided the novel insights into that in vivo adaptive evolution benefits Salmonella colonization, persistence and pathogenesis, by promoting bacterial tolerance via modulating flagella expression. These findings may explain the rationale behind the increasing prevalence of certain MDR Salmonella clones in livestock and associated environment, and underscoring the need for advanced strategies to tackle the possible evolution of such zoonotic pathogens.

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.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs across hundreds of fluconazole-resistant yeast strains

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into 6 classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, and some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings empower multidrug strategies that leverage tradeoffs to combat resistance. Finally, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Mathematical Model for Growth and Rifampicin-Dependent Killing Kinetics of Escherichia coli Cells

Antibiotic resistance is a global health threat. We urgently need better strategies to improve antibiotic use to combat antibiotic resistance. Currently, there are a limited number of antibiotics in the treatment repertoire of existing bacterial infections. Among them, rifampicin is a broad-spectrum antibiotic against various bacterial pathogens. However, during rifampicin exposure, the appearance of persisters or resisters decreases its efficacy. Hence, to benefit more from rifampicin, its current standard dosage might be reconsidered and explored using both computational tools and experimental or clinical studies. In this study, we present the mathematical relationship between the concentration of rifampicin and the growth and killing kinetics of Escherichia coli cells. We generated time-killing curves of E. coli cells in the presence of 4, 16, and 32 μg/mL rifampicin exposures. We specifically focused on the oscillations with decreasing amplitude over time in the growth and killing kinetics of rifampicin-exposed E. coli cells. We propose the solution form of a second-order linear differential equation for a damped oscillator to represent the mathematical relationship. We applied a nonlinear curve fitting solver to time-killing curve data to obtain the model parameters. The results show a high fitting accuracy.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

Novel antibiotic susceptibility of an RNA polymerase α-subunit mutant in Pseudomonas aeruginosa

BACKGROUND

RNA polymerase (RNAP) is highly conserved and essential for prokaryotic housekeeping activities, making it an important target for the development of new antibiotics. The rpoB gene, encoding a β-subunit of bacterial RNAP, has a well-known association with rifampicin resistance. However, the roles of other RNAP component genes such as rpoA, encoding an α-subunit of RNAP, in antibiotic resistance remain unexplored.

OBJECTIVES

To characterize the antibiotic resistance-related role of RpoA.

METHODS

We measured the expression of the MexEF-OprN efflux pump in an RpoA mutant using a transcriptional reporter. The MICs of various antibiotics for this RpoA mutant were determined.

RESULTS

We uncover a novel role of antibiotic susceptibility for an RpoA mutant in Pseudomonas aeruginosa. We found that a single amino acid substitution in RpoA resulted in reduced activity of the MexEF-OprN efflux pump, which is responsible for the exportation of various antibiotics, including ciprofloxacin, chloramphenicol, ofloxacin and norfloxacin. This attenuated efflux pump activity, caused by the RpoA mutation, conferred the bacteria further susceptibility to antibiotics regulated by MexEF-OprN. Our work further revealed that certain clinical P. aeruginosa isolates also contained the same RpoA mutation, providing clinical relevance to our findings. Our results elucidate why this new antibiotic-susceptible function of RpoA mutants would have remained undetected in conventional screens for mutants involving antibiotic resistance.

CONCLUSIONS

The discovery of antibiotic susceptibility in an RpoA mutant implicates a new therapeutic approach for treating clinical isolates of P. aeruginosa with RpoA mutations, using specific antibiotics regulated by MexEF-OprN. More generally, our work suggests that RpoA could serve as a promising candidate target for anti-pathogen therapeutic purposes.

Attenuated but immunostimulatory Mycobacterium tuberculosis variant bovis strain Ravenel shows variation in T cell epitopes

Tuberculosis, caused by Mycobacterium tuberculosis complex (MTBC) organisms, affects a range of humans and animals globally. Mycobacterial pathogenesis involves manipulation of the host immune system, partially through antigen presentation. Epitope sequences across the MTBC are evolutionarily hyperconserved, suggesting their recognition is advantageous for the bacterium. Mycobacterium tuberculosis var. bovis (MBO) strain Ravenel is an isolate known to provoke a robust immune response in cattle, but typically fails to produce lesions and persist. Unlike attenuated MBO BCG strains that lack the critical RD1 genomic region, Ravenel is classic-type MBO structurally, suggesting genetic variation is responsible for defective pathogenesis. This work explores variation in epitope sequences in MBO Ravenel by whole genome sequencing, and contrasts such variation against a fully virulent clinical isolate, MBO strain 10-7428. Validated MTBC epitopes (n = 4818) from the Immune Epitope Database were compared to their sequences in MBO Ravenel and MBO 10-7428. Ravenel yielded 3 modified T cell epitopes, in genes rpfB, argC, and rpoA. These modifications were predicted to have little effect on protein stability. In contrast, T cells epitopes in 10-7428 were all WT. Considering T cell epitope hyperconservation across MTBC variants, these altered MBO Ravenel epitopes support their potential contribution to overall strain attenuation. The affected genes may provide clues on basic pathogenesis, and if so, be feasible targets for reverse vaccinology.

Efficient, robust, and versatile fluctuation data analysis using MLE MUtation Rate calculator (mlemur)

The fluctuation assay remains an important tool for analyzing the levels of mutagenesis in microbial populations. The mutant counts originating from some average number of mutations are usually assumed to obey the Luria-Delbrück distribution. While several tools for estimating mutation rates are available, they sometimes lack accuracy or versatility under non-standard conditions. In this work, extensions to the Luria-Delbrück protocol to account for phenotypic lag and cellular death with either perfect or partial plating were developed. Hence, the novel MLE MUtation Rate calculator, or mlemur, is the first tool that provides a user-friendly graphical interface allowing the researchers to model their data with consideration for partial plating, differential growth of mutants and non-mutants, phenotypic lag, cellular death, variability of the final number of cells, post-exponential-phase mutations, and the size of the inoculum. Additionally, mlemur allows the users to incorporate most of these special conditions at the same time to obtain highly accurate estimates of mutation rates and P values, confidence intervals for an arbitrary function of data (such as fold), and perform power analysis and sample size determination for the likelihood ratio test. The accuracy of point and interval estimates produced by mlemur against historical and simulated fluctuation experiments are assessed. Both mlemur and the analyses in this work might be of great help when evaluating fluctuation experiments and increase the awareness of the limitations of the widely-used Lea-Coulson formulation of the Luria-Delbrück distribution in the more realistic biological contexts.

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.

The evolution of antibiotic resistance is associated with collateral drug phenotypes in Mycobacterium tuberculosis

The increasing incidence of drug resistance in Mycobacterium tuberculosis has diminished the efficacy of almost all available antibiotics, complicating efforts to combat the spread of this global health burden. Alongside the development of new drugs, optimised drug combinations are needed to improve treatment success and prevent the further spread of antibiotic resistance. Typically, antibiotic resistance leads to reduced sensitivity, yet in some cases the evolution of drug resistance can lead to enhanced sensitivity to unrelated drugs. This phenomenon of collateral sensitivity is largely unexplored in M. tuberculosis but has the potential to identify alternative therapeutic strategies to combat drug-resistant strains that are unresponsive to current treatments. Here, by using drug susceptibility profiling, genomics and evolutionary studies we provide evidence for the existence of collateral drug sensitivities in an isogenic collection M. tuberculosis drug-resistant strains. Furthermore, in proof-of-concept studies, we demonstrate how collateral drug phenotypes can be exploited to select against and prevent the emergence of drug-resistant strains. This study highlights that the evolution of drug resistance in M. tuberculosis leads to collateral drug responses that can be exploited to design improved drug regimens.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

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.

Mutations compensating for the fitness cost of rifampicin resistance in Escherichia coli exert pleiotropic effect on RNA polymerase catalysis

The spread of drug-resistant bacteria represents one of the most significant medical problems of our time. Bacterial fitness loss associated with drug resistance can be counteracted by acquisition of secondary mutations, thereby enhancing the virulence of such bacteria. Antibiotic rifampicin (Rif) targets cellular RNA polymerase (RNAP). It is potent broad spectrum drug used for treatment of bacterial infections. We have investigated the compensatory mechanism of the secondary mutations alleviating Rif resistance (Rifr) on biochemical, structural and fitness indices. We find that substitutions in RNAP genes compensating for the growth defect caused by βQ513P and βT563P Rifr mutations significantly enhanced bacterial relative growth rate. By assaying RNAP purified from these strains, we show that compensatory mutations directly stimulated basal transcriptional machinery (2-9-fold) significantly improving promoter clearance step of the transcription pathway as well as elongation rate. Molecular modeling suggests that compensatory mutations affect transcript retention, substrate loading, and nucleotidyl transfer catalysis. Strikingly, one of the identified compensatory substitutions represents mutation conferring rifampicin resistance on its own. This finding reveals an evolutionary process that creates more virulent species by simultaneously improving the fitness and augmenting bacterial drug resistance.

Revisiting Antibiotic Resistance: Mechanistic Foundations to Evolutionary Outlook

Antibiotics are the pivotal pillar of contemporary healthcare and have contributed towards its advancement over the decades. Antibiotic resistance emerged as a critical warning to public wellbeing because of unsuccessful management efforts. Resistance is a natural adaptive tool that offers selection pressure to bacteria, and hence cannot be stopped entirely but rather be slowed down. Antibiotic resistance mutations mostly diminish bacterial reproductive fitness in an environment without antibiotics; however, a fraction of resistant populations 'accidentally' emerge as the fittest and thrive in a specific environmental condition, thus favouring the origin of a successful resistant clone. Therefore, despite the time-to-time amendment of treatment regimens, antibiotic resistance has evolved relentlessly. According to the World Health Organization (WHO), we are rapidly approaching a 'post-antibiotic' era. The knowledge gap about antibiotic resistance and room for progress is evident and unified combating strategies to mitigate the inadvertent trends of resistance seem to be lacking. Hence, a comprehensive understanding of the genetic and evolutionary foundations of antibiotic resistance will be efficacious to implement policies to force-stop the emergence of resistant bacteria and treat already emerged ones. Prediction of possible evolutionary lineages of resistant bacteria could offer an unswerving impact in precision medicine. In this review, we will discuss the key molecular mechanisms of resistance development in clinical settings and their spontaneous evolution.

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.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the Allele Frequency Spectrum (AFS) of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a super-abundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of co-expressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Identification and Characterization of Pleiotropic High-Persistence Mutations in the Beta Subunit of the Bacterial RNA Polymerase

Mutations conferring resistance to bactericidal antibiotics reduce the average susceptibility of mutant populations. It is unknown, however, how those mutations affect the survival of individual bacteria. Since surviving bacteria can be a reservoir for recurring infections, it is important to know how survival rates may be affected by resistance mutations and by the choice of antibiotics. Here, we present evidence that (i) Escherichia coli mutants with 100 to 1,000 times increased frequency of survival in ciprofloxacin, an archetypal fluoroquinolone antibiotic, can be readily obtained in a stepwise selection; (ii) the high survival frequency is conferred by mutations in the switch region of the beta subunit of the RNA polymerase; (iii) the switch-region mutations are (p)ppGpp mimics, partially analogous to rpoB stringent mutations; (iv) the stringent and switch region rpoB mutations frequently occur in clinical isolates of E. coli, Acinetobacter baumannii, Mycobacterium tuberculosis, and Staphylococcus aureus, and at least one of them, RpoB S488L, which is a common rifampicin resistance mutations, dramatically increases the survival of a clinical methicillin-resistant S. aureus (MRSA) strain in ampicillin; and (v) the RpoB-associated high-survival phenotype can be reversed by subinhibitory concentrations of chloramphenicol.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

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.

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.

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.

A Dual-Responsive Antibiotic-Loaded Nanoparticle Specifically Binds Pathogens and Overcomes Antimicrobial-Resistant Infections

Antimicrobial resistant (AMR) infections are a growing threat to public health and there is a general lack of development in new antibiotics. Here, a dextran-coated stimuli-responsive nanoparticle (NP) that encapsulates the hydrophobic antibiotic, rifampicin, and specifically binds bacteria to overcome AMR infections is reported. The NP shows a strong affinity with a variety of pathogens in vitro and effectively accumulates in the bacterial infected tissues. The NP is activated by either low pH or high reactive oxygen species in the infectious microenvironment, and releases both cationic polymer and rifampicin that display synergistic activity against AMR pathogens. The NP carrier also enables the antibiotic to penetrate both bacterial biofilms and mammalian cells, thus allowing the elimination of biofilm and intracellular infections. The NP formulation demonstrates both safety and efficacy in two animal infection models against either Gram-negative or Gram-positive AMR pathogens.

Fidaxomicin has high in vitro activity against Mycobacterium tuberculosis

This study aimed to evaluate whether the antibiotic fidaxomicin has in vitro activity against Mycobacterium tuberculosis (Mtb). 38 fully drug-sensitive Mtb strains and 34 multidrug-resistant tuberculosis (MDR-TB) strains were tested using the microplate alamar blue assay (MABA) method to determine the minimum inhibitory concentrations (MICs) for fidaxomicin and rifampicin. Fidaxomicin has high in vitro activity against Mtb and is a potential drug to treat Mtb, and MDR-TB infections in particular.

Dysbiosis individualizes the fitness effect of antibiotic resistance in the mammalian gut

In the absence of antibiotics, it is essential that antibiotic resistance has a fitness cost for microorganisms if suspending antibiotics treatment is to be a useful strategy for reducing antibiotic resistance. However, the cost of antibiotic resistance within the complex ecosystem of the mammalian gut is not well understood. Here, using mice, we show that the same antibiotic resistance mutation can reduce fitness in one host, while being neutral or even increasing fitness in other hosts. Such antagonistic pleiotropy is shaped by the microbiota because resistance in germ-free mice is consistently costly across all hosts, and the host-specific effect on antibiotic resistance is reduced in hosts with similar microbiotas. Using an eco-evolutionary model of competition for resources, we identify a general mechanism that underlies between-host variation and predicts that the dynamics of compensatory evolution of resistant bacteria should be host specific, a prediction that was supported by experimental evolution in vivo. The microbiome of each human is close to unique, and our results suggest that the short-term cost of resistances and their long-term within-host evolution are also highly personalized, a finding that may contribute to the observed variable outcome of withdrawing antibiotics to reduce resistance levels.

Using a Heat Diffusion Model to Detect Potential Drug Resistance Genes of Mycobacterium tuberculosis

BACKGROUND

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is one of the oldest known and most dangerous diseases. Although the spread of TB was controlled in the early 20th century using antibiotics and vaccines, TB has again become a threat because of increased drug resistance. There is still a lack of effective treatment regimens for a person who is already infected with multidrug-resistant Mtb (MDR-Mtb) or extensively drug-resistant Mtb (XDRMtb). In the past decades, many research groups have explored the drug resistance profiles of Mtb based on sequence data by GWAS, which identified some mutations that were significantly linked with drug resistance, and attempted to explain the resistance mechanisms. However, they mainly focused on several significant mutations in drug targets (e.g. rpoB, katG). Some genes which are potentially associated with drug resistance may be overlooked by the GWAS analysis.

OBJECTIVE

In this article, our motivation is to detect potential drug resistance genes of Mtb using a heat diffusion model.

METHODS

All sequencing data, which contained 127 samples of Mtb, i.e. 34 ethambutol-, 65 isoniazid-, 53 rifampicin- and 45 streptomycin-resistant strains. The raw sequence data were preprocessed using Trimmomatic software and aligned to the Mtb H37Rv reference genome using Bowtie2. From the resulting alignments, SAMtools and VarScan were used to filter sequences and call SNPs. The GWAS was performed by the PLINK package to obtain the significant SNPs, which were mapped to genes. The P-values of genes calculated by GWAS were transferred into a heat vector. The heat vector and the Mtb protein-protein interactions (PPI) derived from the STRING database were inputted into the heat diffusion model to obtain significant subnetworks by HotNet2. Finally, the most significant (P < 0.05) subnetworks associated with different phenotypes were obtained. To verify the change of binding energy between the drug and target before and after mutation, the method of molecular dynamics simulation was performed using the AMBER software.

RESULTS

We identified significant subnetworks in rifampicin-resistant samples. Excitingly, we found rpoB and rpoC, which are drug targets of rifampicin. From the protein structure of rpoB, the mutation location was extremely close to the drug binding site, with a distance of only 3.97 Å. Molecular dynamics simulation revealed that the binding energy of rpoB and rifampicin decreased after D435V mutation. To a large extent, this mutation can influence the affinity of drug-target binding. In addition, topA and pyrG were reported to be linked with drug resistance, and might be new TB drug targets. Other genes that have not yet been reported are worth further study.

CONCLUSION

Using a heat diffusion model in combination with GWAS results and protein-protein interactions, the significantly mutated subnetworks in rifampicin-resistant samples were found. The subnetwork not only contained the known targets of rifampicin (rpoB, rpoC), but also included topA and pyrG, which are potentially associated with drug resistance. Together, these results offer deeper insights into drug resistance of Mtb, and provides potential drug targets for finding new antituberculosis drugs.

CRISPR/Cas9 recombineering-mediated deep mutational scanning of essential genes in Escherichia coli

Deep mutational scanning can provide significant insights into the function of essential genes in bacteria. Here, we developed a high-throughput method for mutating essential genes of Escherichia coli in their native genetic context. We used Cas9-mediated recombineering to introduce a library of mutations, created by error-prone PCR, within a gene fragment on the genome using a single gRNA pre-validated for high efficiency. Tracking mutation frequency through deep sequencing revealed biases in the position and the number of the introduced mutations. We overcame these biases by increasing the homology arm length and blocking mismatch repair to achieve a mutation efficiency of 85% for non-essential genes and 55% for essential genes. These experiments also improved our understanding of poorly characterized recombineering process using dsDNA donors with single nucleotide changes. Finally, we applied our technology to target rpoB, the beta subunit of RNA polymerase, to study resistance against rifampicin. In a single experiment, we validate multiple biochemical and clinical observations made in the previous decades and provide insights into resistance compensation with the study of double mutants.

Differential Expression of Resistant and Efflux Pump Genes in MDR-TB Isolates

Background:

Numerous investigations demonstrate efflux as a worldwide bacterial mode of action which contributes to the resistance of drugs. The activity of antibiotics, which subjects to efflux, can be improved by the combined usage of efflux inhibitors. However, the efflux role to the overall levels of antibiotic resistance of clinical M. tuberculosis isolates is inadequately comprehended and is still disregarded by many.

Method:

Here, we assessed the contribution of resistant genes associated with isoniazid (INH) and rifampin (R) resistance to the levels of drug resistance in the (27) clinical isolates of MDR-TB. Additionally, the role of the resistance for six putative drug efflux pump genes to the antibiotics was investigated. The level of katG expression was down-regulated in 24/27 (88.88%) of MDR-TB isolates. Of the 27 MDR-TB isolates, inhA, oxyR-ahpC, and rpoB showed either overexpression or up-regulation in 8 (29.62%), 4 (14.81 %), and 24 (88.88%), respectively. Moreover, the efflux pump genes drrA, drrB, efpA, Rv2459, Rv1634, and Rv1250 were overexpressed under INH/RIF plus fresh pomegranate juice (FPJ) stress signifying the efflux pumps contribution to the overall levels of the resistance of MDR-TB isolates.

Conclusion:

These results displayed that the levels of drug resistance of MDR-TB clinical isolates are due to combination among drug efflux pump and the presence of mutations in target genes, a truth which is often ignored by the specialists of tuberculosis in favour of the almost undoubted significance of drug target- gene mutations for the resistance in M. tuberculosis.

Genomic Analysis of Mycobacterium tuberculosis Isolates and Construction of a Beijing Lineage Reference Genome

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, kills over 1 million people worldwide annually. Development of drug resistance (DR) in the pathogen is a major challenge for TB control. We conducted whole-genome analysis of seven Taiwan M. tuberculosis isolates: One drug susceptible (DS) and five DR Beijing lineage isolates and one DR Euro-American lineage isolate. Developing a new method for DR mutation identification and applying it to the next-generation sequencing (NGS) data from the 6 Beijing lineage isolates, we identified 13 known and 6 candidate DR mutations and provided experimental support for 4 of them. We assembled the genomes of one DS and two DR Beijing lineage isolates and the Euro-American lineage isolate using NGS data. Moreover, using both PacBio and NGS sequencing data, we obtained a high-quality assembly of an extensive DR Beijing lineage isolate. Comparative analysis of these five newly assembled genomes and two published complete genomes revealed a large number of genetic changes, including gene gains and losses, indels and translocations, suggesting rapid evolution of M. tuberculosis. We found the MazEF toxin–antitoxin system in all the seven isolates studied and several interesting mutations in MazEF proteins. Finally, we used the four assembled Beijing lineage genomes to construct a high-quality Beijing lineage reference genome that is DS and contains all the genes in the four genomes. It contains 212 genes not found in the standard reference H37Rv, which is Euro-American. It is therefore a better reference than H37Rv for the Beijing lineage, the predominant lineage in Asia.

Characteristics of compensatory mutations in the rpoC gene and their association with compensated transmission of Mycobacterium tuberculosis

The aim of this study was to characterize rpoC gene mutations in Mycobacterium tuberculosis (MTB) and investigate the factors associated with rpoC mutations and the relation between rpoC mutations and tuberculosis (TB) transmission. A total of 245 MTB clinical isolates from patients with TB in six provinces and two municipalities in China were characterized based on gene mutations through DNA sequencing of rpoC and rpoB genes, phenotyping via standard drug susceptibility testing, and genotypic profiling by mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing. Approximately 36.4% of the rifampin-resistant isolates harbored nonsynonymous mutations in the rpoC gene. Twenty-nine nonsynonymous single mutations and three double mutations were identified. The rpoC mutations at locus 483 (11.3%) were predominant, and the mutations at V483G, W484G, I491V, L516P, L566R, N698K, and A788E accounted for 54.5% of the total detected mutations. Fifteen new mutations in the rpoC gene were identified. Rifampin resistance and rpoB mutations at locus 531 were significantly associated with rpoC mutations. MIRU-VNTR genotype results indicated that 18.4% of the studied isolates were clustered, and the rpoC mutations were not significantly associated with MIRU-VNTR clusters. A large proportion of rpoC mutation was observed in the rifampicin-resistant MTB isolates. However, the findings of this study do not support the association of rpoC mutation with compensated transmissibility.

Primary drug resistance of mycobacterium tuberculosis in Shandong, China, 2004-2018

BACKGROUND

Primary drug-resistant tuberculosis (DR-TB) has contributed to a significant health and economic burden on a global scale, especially in China. we sought to estimate epidemiological characteristics of primary DR-TB in China from 2004 to 2018.

METHODS

Eleven thousand four hundred sixty-seven newly diagnosed and 1981 retreated TB cases with drug susceptibility data were included. Chi-Square test for trends, linear regression, a joinpoint regression model and temporal trend in proportions of the different resistance patterns were carried out.

RESULTS

The proportion of primary DR-TB and mono-resistant TB (MR-TB) in China had reduced by more than 12% since 2004, and were 21.38%, 13.35% in 2018 respectively. Among primary DR-TB cases (2173,18.95%), the percentage of multiresistant TB (MDR-TB, from 5.41 to 17.46%), male (from 77.03 to 84.13%), cavity (from 13.51 to 43.92%), rifampicin(RFP)-resistant TB (from 8.11 to 26.98%), streptomycin(SM)-resistant TB (from 50.00 to 71.43%) increased significantly (P < 0.05). On the contrary, the proportion of female, non-cavity, isoniazide(INH)-resistant TB (from 55.41 to 48.15%) and MR-TB (from 82.43 to 62.43%) decreased significant (P < 0.05). The primary drug resistance rate among female, cavity, smoking, drinking, 15 to 44 year-old TB subgroups increased by 0.16, 6.24, 20.95, 158.85, 31.49%, respectively. The percentage of primary DR-TB, RFP-resistant TB dropped significantly during 2004-2007 in Joinpoint regression model.

CONCLUSION

The total rate of drug resistance among new TB cases showed a downward trend in Shandong, China, from 2004 to 2018. Primary drug resistance patterns were shifting from female, non-cavity, INH-resistant TB, and MR-TB groups to male, cavity, RFP/SM-resistant TB, and MDR-TB groups. Considering the rising drug resistance rate among some special population, future control of primary DR-TB in China may require an increased focus on female, cavity, smoking, drinking, or 15 to 44 year-old TB subgroups.

Rapid decline of bacterial drug-resistance in an antibiotic-free environment through phenotypic reversion

Antibiotic resistance typically induces a fitness cost that shapes the fate of antibiotic-resistant bacterial populations. However, the cost of resistance can be mitigated by compensatory mutations elsewhere in the genome, and therefore the loss of resistance may proceed too slowly to be of practical importance. We present our study on the efficacy and phenotypic impact of compensatory evolution in Escherichia coli strains carrying multiple resistance mutations. We have demonstrated that drug-resistance frequently declines within 480 generations during exposure to an antibiotic-free environment. The extent of resistance loss was found to be generally antibiotic-specific, driven by mutations that reduce both resistance level and fitness costs of antibiotic-resistance mutations. We conclude that phenotypic reversion to the antibiotic-sensitive state can be mediated by the acquisition of additional mutations, while maintaining the original resistance mutations. Our study indicates that restricting antimicrobial usage could be a useful policy, but for certain antibiotics only.

In vivo exposure of insect AMP resistant Staphylococcus aureus to an insect immune system

Antimicrobial peptides (AMPs) are important immune effectors in insects. Bacteria have a limited number of ways to resist AMPs, and AMP-resistance is often costly. Recently, it has become clear that AMP activities in vitro and in vivo differ. Although some studies have followed the in vivo survival of AMP resistant pathogens, studying a pathogen resistant to the AMPs of that particular host has never been reported. Here, we infected the mealworm beetle Tenebrio molitor with Staphylococcus aureus strains that were evolved in vitro in the presence of one or two antimicrobial peptides from T. molitor. We found that the Tenebrio immune system could clear mutant Tenecin resistant strains at least as efficiently as sensitive controls. The bacterial load of Tenecin resistant S. aureus segregated by mutation. Strains with mutations in both the pmt and rpo operons showed the highest in vivo survival and therefore showed the lowest fitness cost amongst the evolved resistance mutations. In contrast, Tenecin resistant strains with mutations in the nsa and rpo operons showed much lower survival within the hosts. Our study shows that Tenecin resistant strains are phagocytosed at a lower rate. The nsa/rpo mutants were phagocytosed at a higher rate than other Tenecin resistant S. aureus strains. The differences in resistance against AMPs and phagocytosis did not translate into changes in virulence. AMP resistance, while a prerequisite for an infection in vertebrates, does not provide a survival advantage to S. aureus in a host environment that is dominated by AMPs.

Fitness of Spontaneous Rifampicin-Resistant Staphylococcus aureus Isolates in a Biofilm Environment

Biofilms of S. aureus accumulate cells resistant to the antibiotic rifampicin. We show here that the accumulation of rifampicin resistant mutants (RifR) in biofilms is not equable but rather is a local event, suggesting that the growth of a few locally emerged mutants is responsible for this. Competition assays demonstrated that, compared to wild-type bacteria, the isolated RifR mutants have a growth advantage in biofilms, but not in planktonic culture. To gain insight into the mechanism of the growth advantage, we tested the involvement of the two-component systems (TCS) that sense and respond to environmental changes. We found that a deletion of SrrAB or NreBC has a drastic effect on the growth advantage of RifR mutants, suggesting the importance of oxygen/respiration responses. All six of the RifR isolates tested showed increased resistance to at least one of the common stresses found in the biofilm environment (i.e., oxidative, nitric acid, and organic acid stress). The RifR mutants also had a growth advantage in a biofilm flow model, which highlights the physiological relevance of our findings.

Deciphering the essentiality and function of the anti-σM factors in Bacillus subtilis

Bacteria use alternative sigma factors to adapt to different growth and stress conditions. The Bacillus subtilis extracytoplasmic function sigma factor SigM regulates genes for cell wall synthesis and is crucial for maintaining cell wall homeostasis under stress conditions. The activity of SigM is regulated by its anti-sigma factor, YhdL, and the accessory protein YhdK. Here, we show that dysregulation of SigM caused by the absence of either component of the anti-sigma factor complex leads to toxic levels of SigM and severe growth defects. High SigM activity results from a dysregulated positive feedback loop, and can be suppressed by overexpression of the housekeeping sigma, SigA. Using a sigM merodiploid strain, we selected for suppressor mutations that allow survival of yhdL depletion strain. The recovered suppressor mutations map to the beta and beta-prime subunits of RNA polymerase core enzyme and selectively reduce SigM activity, and in some cases increase the activity of other alternative sigma factors. This work highlights the ability of mutations in RNA polymerase that remodel the sigma-core interface to differentially affect sigma factor activity, and thereby alter the transcriptional landscape of the cell.

Alterations in Life-History Associated With Non-target-site Herbicide Resistance in Alopecurus myosuroides

The evolution of resistance to herbicides is a classic example of rapid contemporary adaptation in the face of a novel environmental stress. Evolutionary theory predicts that selection for resistance will be accompanied by fitness trade-offs in environments where the stress is absent. Alopecurus myosuroides, an autumn-germinating grass weed of cereal crops in North-West Europe, has evolved resistance to seven herbicide modes-of-action, making this an ideal species to examine the presence and magnitudes of such fitness costs. Here, we use two contrasting A. myosuroides phenotypes derived from a common genetic background, one with enhanced metabolism resistance to a commercial formulation of the sulfonylurea (ALS) actives mesosulfuron and iodosulfuron, and the other with susceptibility to these actives (S). Comparisons of plant establishment, growth, and reproductive potential were made under conditions of intraspecific competition, interspecific competition with wheat, and over a gradient of nitrogen deprivation. Herbicide dose response assays confirmed that the two lines had contrasting resistance phenotypes, with a 20-fold difference in resistance between them. Pleiotropic effects of resistance were observed during plant development, with R plants having a greater intraspecific competitive effect and longer tiller lengths than S plants during vegetative growth, but with S plants allocating proportionally more biomass to reproductive tissues during flowering. Direct evidence of a reproductive cost of resistance was evident in the nitrogen deprivation experiment with R plants producing 27% fewer seed heads per plant, and a corresponding 23% reduction in total seed head length. However, these direct effects of resistance on fecundity were not consistent across experiments. Our results demonstrate that a resistance phenotype based on enhanced herbicide metabolism has pleiotropic impacts on plant growth, development and resource partitioning but does not support the hypothesis that resistance is associated with a consistent reproductive fitness cost in this species. Given the continued difficulties associated with unequivocally detecting costs of herbicide resistance, we advocate future studies that adopt classical evolutionary quantitative genetics approaches to determine genetic correlations between resistance and fitness-related plant life history traits.

Analysis of immune response in the sea urchin larva

Sea urchin larvae deploy a complex immune system in the context of relatively simple morphology. Several types of phagocytic or granular immune cells respond rapidly to microbes and microbial components within the body cavity. Many of these cells also respond to microbial disturbances in the gut lumen. In the course of immune response, hundreds of genes are up- and downregulated, many of which have homologs involved in immunity in other species. Thus, the larval sea urchin provides an experimentally advantageous model for investigating the response to immune challenge at the level of cell behavior and gene regulatory networks. Importantly, the morphological simplicity and optical clarity of these larvae allow studies to be carried out within the intact animal. Here, we outline techniques to probe and visualize the immune system of the feeding sea urchin larva, particularly for quantifying gene expression and cell migration as the animal responds to both pathogens and symbionts. Techniques addressed in this chapter include (1) exposure of larvae to microbes and microbial products in sea water and by blastocoelar microinjection, (2) time-lapse imaging of immune response, (3) isolation of culturable bacteria associated with feeding larvae, (4) quantification of larval associations with isolated bacterial strains and (5) preparation of secreted products from isolated bacteria for testing in larval culture.

Mechanisms of fitness cost reduction for rifampicin-resistant strains with deletion or duplication mutations in rpoB

Rifampicin resistance (Rif^R) is caused by mutations in rpoB, encoding the β-subunit of RNA polymerase. Rif^R mutations generally incur a fitness cost and in resistant isolates are frequently accompanied by compensatory mutations in rpoA, rpoB or rpoC. Previous studies of fitness compensation focused on Rif^R caused by amino acid substitutions within rpoB. Rif^R is also caused by deletion and duplication mutations in rpoB but it is not known whether or how such mutants can ameliorate their fitness costs. Using experimental evolution of Salmonella carrying Rif^R deletion or duplication mutations we identified compensatory amino acid substitution mutations within rpoA, rpoB or rpoC in 16 of 21 evolved lineages. Additionally, we found one lineage where a large deletion was compensated by duplication of adjacent amino acids (possibly to fill the gap within the protein structure), two lineages where mutations occurred outside of rpoABC, and two lineages where a duplication mutant reverted to the wild-type sequence. All but the two revertant mutants maintained the Rif^R phenotype. These data suggest that amino acid substitution mutations are the major compensatory mechanism regardless of the nature of the primary Rif^R mutation.

Reduction of the fitness cost of antibiotic resistance caused by chromosomal mutations under poor nutrient conditions

The prevalence of antibiotic resistance in drinking water system is pressing public health risk. Antibiotic resistance conferred by chromosomal mutations often produces fitness cost, which may affect its spread and persistence. In this study, the rifampin-resistant strains were competed with their wild-type counterparts at different nutrient levels. It was observed that the ratio of the absolute number between resistant and wild-type cells quickly decreased under rich nutrient conditions, but it slowly reduced or remained stable in the poor nutrient medium. This finding suggested that poor nutrient conditions resulted in the reduction of fitness cost of antibiotic resistance, i.e. the resistant bacteria became more competitive. Implying mechanisms analysis found that the differences of metabolic activity between wild-type and rifampin-resistant strains was significant smaller (P < 0.05) at low nutrient levels. Additionally, distinguishable large colony size rifampin-resistant strains were observed during competition assay. DNA sequencing of RNA polymerase subunit genes further revealed that these colonies could be adaptive mutants from wild-type strain in rpoB gene. To our knowledge, this is the first study to reveal that the oligotrophic conditions facilitate the persistence of antibiotic resistance in drinking water by reducing the fitness cost of the resistant strains.

Insights into the processes that drive the evolution of drug resistance in Mycobacterium tuberculosis

At present, the successful transmission of drug-resistant Mycobacterium tuberculosis, including multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains, in human populations, threatens tuberculosis control worldwide. Differently from many other bacteria, M. tuberculosis drug resistance is acquired mainly through mutations in specific drug resistance-associated genes. The panel of mutations is highly diverse, but depends on the affected gene and M. tuberculosis genetic background. The variety of genetic profiles observed in drug-resistant clinical isolates underlines different evolutionary trajectories towards multiple drug resistance, although some mutation patterns are prominent. This review discusses the intrinsic processes that may influence drug resistance evolution in M. tuberculosis, such as mutation rate, drug resistance-associated mutations, fitness cost, compensatory mutations and epistasis. This knowledge should help to better predict the risk of emergence of highly resistant M. tuberculosis strains and to develop new tools and strategies to limit the development and spread of MDR and XDR strains.

Genetics and roadblocks of drug resistant tuberculosis

Considering the extensive evolutionary history of Mycobacterium tuberculosis, anti-Tuberculosis (TB) drug therapy exerts a recent selective pressure. However, in a microorganism devoid of horizontal gene transfer and with a strictly clonal populational structure such as M. tuberculosis the usual, but not sole, path to overcome drug susceptibility is through de novo mutations on a relatively strict set of genes. The possible allelic diversity that can be associated with drug resistance through several mechanisms such as target alteration or target overexpression, will dictate how these genes can become associated with drug resistance. The success demonstrated by this pathogenic microbe in this latter process and its ability to spread is currently one of the major obstacles to an effective TB elimination. This article reviews the action mechanism of the more important anti-TB drugs, including bedaquiline and delamanid, along with new findings on specific resistance mechanisms. With the development, validation and endorsement of new in vitro molecular tests for drug resistance, knowledge on these resistance mechanisms and microevolutionary dynamics leading to the emergence and fixation of drug resistance mutations within the host is highly important. Additionally, the fitness toll imposed by resistance development is also herein discussed together with known compensatory mechanisms. By elucidating the possible mechanisms that enable one strain to reacquire the original fitness levels, it will be theoretically possible to make more informed decisions and develop novel strategies that can force M. tuberculosis microevolutionary trajectory down through a path of decreasing fitness levels.

The Evolution of Gene Expression in cis and trans

There is abundant variation in gene expression between individuals, populations, and species. The evolution of gene regulation and expression within and between species is thought to frequently contribute to adaptation. Yet considerable evidence suggests that the primary evolutionary force acting on variation in gene expression is stabilizing selection. We review here the results of recent studies characterizing the evolution of gene expression occurring in cis (via linked polymorphisms) or in trans (through diffusible products of other genes) and their contribution to adaptation and response to the environment. We review the evidence for buffering of variation in gene expression at the level of both transcription and translation, and the possible mechanisms for this buffering. Lastly, we summarize unresolved questions about the evolution of gene regulation.