Adaptation and compensation in a bacterial gene regulatory network evolving under antibiotic selection
Abstract
Gene regulatory networks allow organisms to generate coordinated responses to environmental challenges. In bacteria, regulatory networks are re-wired and re-purposed during evolution, though the relationship between selection pressures and evolutionary change is poorly understood. In this study, we discover that the early evolutionary response of Escherichia coli to the antibiotic trimethoprim involves derepression of PhoPQ signaling, an Mg2+-sensitive two-component system, by inactivation of the MgrB feedback-regulatory protein. We report that derepression of PhoPQ confers trimethoprim-tolerance to E. coli by hitherto unrecognized transcriptional upregulation of dihydrofolate reductase (DHFR), target of trimethoprim. As a result, mutations in mgrB precede and facilitate the evolution of drug resistance. Using laboratory evolution, genome sequencing, and mutation re-construction, we show that populations of E. coli challenged with trimethoprim are faced with the evolutionary ‘choice’ of transitioning from tolerant to resistant by mutations in DHFR, or compensating for the fitness costs of PhoPQ derepression by inactivating the RpoS sigma factor, itself a PhoPQ-target. Outcomes at this evolutionary branch-point are determined by the strength of antibiotic selection, such that high pressures favor resistance, while low pressures favor cost compensation. Our results relate evolutionary changes in bacterial gene regulatory networks to strength of selection and provide mechanistic evidence to substantiate this link.
Introduction
The relationship between genotype and phenotype is complex and much of modern genetics aims at understanding how it is established and maintained. In cellular organisms, this relationship is dictated in large part by the regulation of gene expression in response to the internal and external environment (Hill et al., 2021; Kemble et al., 2019; Kim and Przytycka, 2012). By regulating gene expression, organisms select which traits are expressed, enabling a wide palette of phenotypes to be generated from the same genotype. Not surprisingly, the functioning and evolution of mechanisms that regulate gene expression are central questions at the cusp of genetics and evolutionary biology with implications for diverse phenomena such as development (Levine and Davidson, 2005), differentiation (Levine and Davidson, 2005), and disease (Singh et al., 2018).
In bacteria, like in more complex organisms, gene regulatory proteins, such as transcription factors and signal transducers are organized into networks (Shen-Orr et al., 2002). The predominant sensory modules in bacteria are two-component signaling systems, which are made up of membrane-bound receptors and cognate cytosolic response regulators (Perraud et al., 1999). Upon activation of the receptor by an appropriate ligand, a phospho-relay between the two components is initiated, which changes the phosphorylation state of the response regulator. Response regulators are often transcription factors themselves, which bind to DNA upon phosphorylation and produce changes in gene expression (Perraud et al., 1999). Individual two-component pathways are integrated into networks by ‘connector’ proteins, and also by cross-activation (Mitrophanov and Groisman, 2008). Further, two-component systems are integrated with other signaling molecules in bacteria such as regulatory RNAs (Valverde and Haas, 2008) and sigma factors (Hengge, 2008).
Gene regulatory networks orchestrate a number of physiological functions in bacteria such as sporulation (Stephenson and Lewis, 2005), virulence (Beier and Gross, 2006), dormancy (Sivaramakrishnan and de Montellano, 2013), and tolerance to stresses like acid (Xu et al., 2020), heavy metal (Perron et al., 2004), and osmolarity (Yuan et al., 2017). More recently, the role of bacterial gene regulatory networks in antimicrobial resistance is beginning to emerge (Tierney and Rather, 2019). VanRS (vancomycin resistance, Staphylococcus aureus) (Jansen et al., 2007), BaeRS (ciprofloxacin resistance, Salmonella typhimurium) (Guerrero et al., 2013), PhoPQ (Bader et al., 2005; Yadavalli et al., 2016), PmrAB (Gunn, 2008) (antimicrobial peptide resistance, Salmonella and Escherichia coli), MtrAB (multidrug resistance, Mycobacterium tuberculosis) (Nguyen et al., 2010), and CroRS (β-lactam resistance, Enterococcus faecalis) (Comenge et al., 2003) are some of the two-component systems that mediate drug resistance in bacteria. Two-component systems modulate antibiotic resistance by regulating the expression of downstream effectors, such as efflux pumps, porins, or enzymes that modify the cell envelope (Tierney and Rather, 2019). Not surprisingly, two-component sensor kinases have also been proposed as potential targets for the design of novel anti-microbials (Bem et al., 2015).
The repertoire of two-component systems that bacteria possess depends on environment and life history. For instance, Myxococcus xanthus that faces rapidly changing, diverse environments codes for more than 250 different two-component systems (Shi et al., 2008). Further, rewiring of two-component systems during evolution is also seen among bacteria. For example, the regulons of Mg2+-sensitive PhoPQ systems from closely related members of Enterobacteriaceae show several species-specific features. Most notably, Salmonella and Yersinia have only a marginal overlap in PhoP-regulons, even though they share a common ancestor (Perez et al., 2009). Similarly, co-option of existing regulatory modules for the expression of genes coded by horizontally acquired pathogenicity islands in Salmonella has also been described (Palmer et al., 2019). Thus, two-component systems and their associated transcriptional networks are likely to be ready substrates for evolution to suit the unique demands of different growth environments.
Current understanding of the evolution of two-component signaling systems and their associated gene regulatory networks comes from two main approaches. First of these relies on the comparison of orthologous genes from related bacterial species (Perez et al., 2009; Monsieurs et al., 2005). The second uses genetically engineered bacterial strains to understand general principles governing the specificity of receptor-regulator pairs and the potential for cross-activation in two-component systems (Schmidl et al., 2019). While invaluable insight has emerged from both approaches, the direct link between evolutionary changes in gene regulatory networks and selection pressures remains elusive. In this study, we use laboratory evolution to identify a specific environmental factor, that is, the antibiotic trimethoprim, as a selective pressure that leads to derepression of the PhoPQ signaling pathway in E. coli. This adaptation is mediated by loss-of-function mutations at the mgrB locus, which codes for a feedback attenuator of PhoPQ signaling. Intrinsic resistance to trimethoprim in E. coli is thought to be acquired primarily by mutations in the drug-target dihydrofolate reductase (DHFR), coded by the folA gene (Toprak et al., 2011; Watson et al., 2007). We find that mutations in DHFR are typically preceded by mutations at the mgrB locus and together, they have a synergistic effect on the level of trimethoprim resistance. We use this system to ask how the costs of adaptation by perturbation of gene regulatory networks are compensated over long-term evolution and identify the key role that strength of selection plays in determining the adaptive trajectory that bacteria take upon antibiotic challenge.
Results
Mutation at the mgrB locus rather than folA is an early evolutionary response of E. coli challenged with trimethoprim
We had earlier evolved trimethoprim-resistance in E. coli K-12 MG1655 by serially passaging bacteria in drug-supplemented medium for ~25 generations (Matange et al., 2018). From these experiments, 10 independently evolved trimethoprim-resistant isolates (designated TMPR1-10), with drug IC50 ranging from 6- to 300-fold over wild-type were obtained (Figure 1A). Unexpectedly, only 3 of the 10 trimethoprim-resistant isolates (TMPR 1, 3, and 5) had mutations in the coding region of the folA gene, which suggested that early adaptation to trimethoprim did not necessarily involve mutations in the drug target itself (Figure 1A). We, therefore, sequenced the genomes of 5 of the 10 isolates (TMPR1–5) to identify other genes that may contribute to trimethoprim resistance. Only one locus, coding for the mgrB feedback regulator protein, consistently harbored mutations in all five isolates (Table 1). Four of the five isolates (TMPR1–4) harbored IS-element insertions in the promoter of mgrB between the PhoP-binding site and the translational start site of the gene (Figure 1A, Table 1). Similar insertions of IS-elements in the promoter of mgrB attenuate its expression levels in other bacteria such as Klebsiella (Pitt et al., 2018; Uz Zaman et al., 2018). The fifth isolate (TMPR5) harbored a single nucleotide deletion in the coding region of mgrB that resulted in a frame-shift and altered sequence of the C-terminal domain of the protein (Figure 1A, Table 1). The C-terminus of MgrB is necessary for its inhibitory activity against PhoQ and mutations in this region alter the ability of MgrB to bind to PhoQ (Yadavalli et al., 2020). Sanger sequencing of the mgrB gene and its promoter from isolates TMPR6–10 confirmed the presence of similar mutations in these bacteria as well (Figure 1A), suggesting that lower MgrB expression or activity may be advantageous for E. coli in the presence of trimethoprim.
We tested this hypothesis using an independently generated ΔmgrB strain and found that loss of MgrB enhanced trimethoprim IC50 by ~3-fold (Figure 1B and C). Thus, mutation-driven loss of mgrB, rather than mutations in the folA gene, was the predominant early adaptive change in E. coli challenged with trimethoprim.
Loss of mgrB confers trimethoprim tolerance by PhoP-dependent upregulation of DHFR expression
To better understand how loss of mgrB reduced the effectiveness of trimethoprim against E. coli, we further characterized the mgrB knock-out strain. E. coli ΔmgrB had lower relative fitness than wild-type in antibiotic-free media (wΔmgrB=0.86±0.12). However, sub-inhibitory concentration (300 ng/ml) trimethoprim in growth media enhanced relative fitness of the mutant (wΔmgrB=1.4±0.08). Concomitantly, E. coli ΔmgrB was selected over wild-type in co-cultures when trimethoprim was added to growth media. The competitive advantage of E. coli ΔmgrB was apparent even when mixing ratios were biased in favor of wild-type (Figure 1D).
Interestingly, there was no detectable difference in trimethoprim minimum inhibitory concentration (MIC) between wild-type and mutant in broth dilution assay or E-test (Figure 1B and E). However, mgrB-deficient E. coli had higher colony-forming efficiency than wild-type on trimethoprim-supplemented solid media even at concentrations close to MIC (Figure 1E). Similarly, relatively fewer cells of E. coli ΔmgrB were required to colonize liquid growth media containing inhibitory concentrations of trimethoprim (Figure 1—figure supplement 1). Thus, loss of mgrB allowed better survival of E. coli over a wide range of trimethoprim concentrations, without a detectable increase in MIC. Therefore, we concluded that mgrB-deficiency conferred trimethoprim tolerance, rather than resistance, to E. coli. Tolerance has been understood recently as a change in the death rate of bacteria by bactericidal antibiotics without change in MIC. However, since trimethoprim is a bacteriostatic antibiotic, we define ‘tolerance’ for trimethoprim as reduced susceptibility to the antibiotic, without substantial change in MIC.
The primary role of the MgrB protein in E. coli is to attenuate PhoPQ signaling through negative feedback (Lippa and Goulian, 2009; Salazar et al., 2016). Thus, PhoPQ signaling is derepressed in an mgrB knockout strain (Lippa and Goulian, 2009; Salazar et al., 2016). In line with this role, deletion of phoP re-sensitized E. coli ΔmgrB to trimethoprim, while in a wild-type background had no detectable effect (Figure 1C, D and E, Figure 1—figure supplement 1). PhoP directly activates transcription of close to 50 genes in E. coli with diverse functions (Minagawa et al., 2003). Among these is iraM, which protects the stress-responsive sigma factor RpoS from degradation (Xu et al., 2019). Arguing that trimethoprim tolerance in the mgrB-knockout could be a result of enhanced general stress-response pathway, we generated knockout strains of iraM and rpoS in a ΔmgrB background. However, neither deletion altered trimethoprim susceptibility of ΔmgrB, ruling out this possibility (Figure 1C, D and E, Figure 1—figure supplement 1).
To identify the target of PhoP that led to trimethoprim-tolerance, we turned to the transcriptomics data published by Xu et al., 2019 that compared gene expression profiles of wild-type and mgrB-deficient E. coli. Among over 500 genes upregulated in mgrB-deficient E. coli, 97 are known targets of RpoS and hence could be excluded (Figure 2—figure supplement 1). Curiously, folA transcript levels were reported to be enhanced by ~3-fold in an mgrB-deficient strain (Xu et al., 2019). The folA gene codes for DHFR, target of trimethoprim, and its overexpression results in trimethoprim resistance (Flensburg and Sköld, 1987). We found that DHFR protein levels were indeed elevated in an mgrB knockout strain (Figure 2A). Further, DHFR levels could be lowered in the mgrB-knockout by deletion of phoP or phoQ, but not rpoS or iraM (Figure 2A). The increase in DHFR levels was due to greater activity of the folA promoter in E. coli ΔmgrB compared to wild-type (Figure 2B and C), and this too was phoP-dependent but rpoS-independent (Figure 2C). Importantly, the folA promoter retained trimethoprim-stimulation in the mgrB-knockout (Figure 2C). The folA promoter is not a known target for PhoP, and does not harbor sequences similar to the PhoP binding site. Consequently, we were unable to detect direct binding of PhoP to the folA promoter in vitro by gel shift assays (Figure 2—figure supplement 2). Thus, the enhancement of folA expression by PhoP was most probably an indirect effect.
In line with higher expression of DHFR in mgrB-deficient E. coli, 9 of the 10 trimethoprim-resistant isolates (TMPR1–2, 4–10) showed elevated DHFR protein levels (Figure 2D). The only isolate that did not (TMPR3), harbored the W30G mutation in folA, which we have shown earlier to result in lower steady-state levels of DHFR (Matange et al., 2018; Matange, 2020). Further, deletion of phoP resulted in lower expression of DHFR in TMPR1–5 (Figure 2E). These data demonstrated that loss of mgrB conferred trimethoprim tolerance in E. coli by enhancing DHFR protein levels, through PhoP-dependent, RpoS-independent upregulation of the folA promoter.
Derepression of PhoPQ facilitates resistance evolution by altering the fitness landscape of folA mutations
In addition to mgrB mutations, isolates TMPR1, 3, and 5 harbored mutations in folA (Table 1). Likewise, TMPR4 harbored a large genomic duplication encompassing the folA gene along with other genes that may confer trimethoprim resistance (Table 1). To assess the contribution of loss of mgrB to the phenotypes of these isolates, we tested trimethoprim MIC and IC50 of ΔphoP derivatives of TMPR1–5. Remarkably, loss of phoP greatly reduced the MIC of trimethoprim for all five isolates (Figure 3A). TMPR2 and 4 that did not harbor mutations in folA were completely re-sensitized to trimethoprim by phoP deletion (Figure 3A and B). TMPR1, 3, and 5, which harbored mutations in folA retained some trimethoprim-resistance after deletion of phoP, but with close to tenfold reduction in drug IC50 (Figure 3B). Thus, even though loss of mgrB itself did not appreciably alter drug MIC, it significantly potentiated the phenotypes of folA mutations in trimethoprim-resistant E. coli.
The above result prompted us to ask whether mutation at mgrB influenced the selection dynamics of trimethoprim resistance in E. coli. The genetic constitution of antibiotic-resistant bacteria is known to be influenced by drug concentration (Matange et al., 2019; Wistrand-Yuen et al., 2018; Zhou et al., 2000). Hence, we mapped the fitness landscape of drug-tolerant (E. coli ΔmgrB) and drug-resistant (TMPR1, 3, and 5) strains across different concentrations of trimethoprim (Figure 3C). Each of the above strains was allowed to compete against a GFP-expressing derivative of E. coli wild-type (E. coli-GFP). Since GFP fluorescence of the mixed culture would indicate the fractional abundance of E. coli-GFP in the population, we used it as a read-out of relative fitness of the test strains (designated as competitive quotient; a value of 0 indicated no change in fitness relative to wild-type, >0 signified higher relative fitness, <0 signified reduced relative fitness). Over the entire range of trimethoprim concentrations used by us, TMPR1, 3, and 5 performed better than ΔmgrB alone. This difference was most pronounced at concentrations that approached the MIC value for the wild-type. Next, we asked how ΔphoP derivatives of TMPR1, 3, and 5 would perform in a similar assay (Figure 3D). Interestingly, ΔmgrB was fitter than ΔphoP derivatives of TMPR1, 3, and 5 across all concentrations of trimethoprim tested by us. These results explained why mgrB mutations were the predominant early adaptive event in our selection experiments and indicated that resistance-conferring DHFR mutations were likely to be selected only after PhoPQ signaling had been derepressed.
The above experiments showed that mutations in folA and mgrB had a synergistic effect on the phenotype of trimethoprim-resistant E. coli. We had shown earlier that several mutations in DHFR, notably those at Trp30, are detrimental for its stability and result in aggregation, proteolysis, and reduced levels (Matange et al., 2018). Therefore, we hypothesized that enhanced expression due to PhoPQ derepression could compensate for loss of mutant DHFRs due to misfolding. To test this, we expressed wild-type DHFR or its W30G/W30R mutant alleles from an IPTG-inducible promoter and checked the impact that expression level had on drug IC50 (Figure 3E). While trimethoprim IC50 was enhanced by ~3-fold upon overproduction of wild-type DHFR, it was potentiated by ×16- and ×8-fold for W30G and W30R DHFR alleles, respectively. We performed the same assay with W30F and W30L mutant alleles of DHFR. These alleles have comparable in vivo sensitivity to trimethoprim, but only W30L destabilizes the protein (Matange et al., 2018). Overexpression affected the phenotype of W30F to a similar extent as wild-type DHFR, while the IC50 of trimethoprim for E. coli expressing W30L DHFR was stimulated ×8-fold upon overexpression (Figure 3E). Thus, mutant DHFR alleles resulted in higher gains in drug IC50 upon overexpression, explaining the synergy between mgrB and folA mutations, and this was related to their destabilizing effects on DHFR.
Cost compensation drives evolution of the PhoPQ transcriptional network during long-term antibiotic exposure
Having established a new mechanism for trimethoprim tolerance, and its effect as a potentiator of trimethoprim resistance, we next sought to understand the long-term evolutionary consequences of this adaptation. We established three lineages of E. coli evolving in high concentration of trimethoprim (300 ng/ml, ×0.5-fold MIC) for ~350 generations (designated WTMP300 A, B, and C). In agreement with results from short-term trimethoprim exposures (Figure 1A, Table 1), mutations in the mgrB gene or its promoter swept through all three lineages within the first ~20 generations of antibiotic exposure (Supplementary file 1). Subsequently, at this high drug pressure, there was rapid enrichment of resistant bacteria that were eventually fixed in all three lineages (Figure 4A).
Next, we tracked adaptation of the evolving lineages to trimethoprim by competing them against E. coli-GFP. For all three lineages, competitive quotient rose steeply initially and then plateaued (Figure 4B). The initial rapid phase of adaptation was concordant with enrichment of trimethoprim-resistant bacteria (Figure 4A and B). Interestingly, IC50 of trimethoprim for resistant isolates from all three lineages plateaued early during evolution and did not change significantly over time (Figure 4—figure supplement 1). Thus, once resistance had evolved, its fixation in the population did not require further gains in drug MIC/IC50. On the other hand, competitive quotients of evolving populations in trimethoprim-free media initially decreased (Figure 4B), most likely due to the costs of mutations at the mgrB locus (Figure 1D). However, all three lineages recovered over time such that at 350 generations of evolution their relative fitness was similar to or slightly better than the ancestor (Figure 4B). Taken together, we interpreted these results to mean that trimethoprim-resistant bacteria enhanced their fitness over long-term evolution by amelioration of fitness costs, rather than enhancement in drug IC50.
To understand the genetic basis of these adaptations, we sequenced the genomes of up to five resistant isolates from lineage A at three different time points, that is, 105, 147, and 350 generations of evolution (Supplementary file 1, Figure 4C). This lineage was chosen for genetic analyses, since it showed the slowest accumulation of resistant bacteria, providing an opportunity to capture the sequence of mutations accumulated during resistance evolution. Sequenced isolates from all three time points had mutations at the mgrB locus, re-iterating its role as a driver of early adaptation to trimethoprim (Supplementary file 1, Figure 4C). Mutations at the folA locus were also common among isolates at all three time points, explaining the high IC50 of trimethoprim for these isolates (Supplementary file 1, Figure 4C). In addition, some isolates harbored mutations at the baeS locus, a known regulator of efflux pump expression (Nagakubo et al., 2002), and safA, a regulator of PhoPQ signaling (Eguchi et al., 2011; Eguchi et al., 2007; Supplementary file 1, Figure 4C). Strikingly, all sequenced isolates from 350 generations, but not earlier time points, harbored substitution of Ile128 to Ser or Asn in RpoS (Supplementary file 1, Figure 4C). The Ile128 residue of the RpoS sigma factor is critical for its activity, and hydrophilic substitutions at this site inactivate it (Iwase et al., 2017). In order to confirm loss of RpoS activity in these lineages, we exploited the differential staining of RpoS-deficient and RpoS-expressing E. coli by Congo red. Congo red stains curli fibers and cellulose, which are produced only by bacteria with active RpoS (Smith et al., 2017). In agreement with results from genome-sequencing, all three lineages showed progressively poorer staining with Congo red (Figure 4D). We also PCR amplified and Sanger sequenced the rpoS gene from three trimethoprim isolates each derived from lineages B and C at 350 generations. Five of the six isolates sequenced harbored mutations at the rpoS locus (Figure 4E). Thus, while resistance-associated genetic changes were observed early during evolution, loss-of-function mutations in rpoS were associated with adaptation during long-term exposure to trimethoprim.
Loss of RpoS, we had found, did not itself confer trimethoprim resistance (Figure 1C). Further, since RpoS is stimulated by PhoPQ signaling, it was reasonable to expect that loss of RpoS activity compensated for the costs of mgrB-mutations. In support of this suggestion, RpoS-deficient bacteria were selected over RpoS-expressing bacteria only in a ΔmgrB background (Figure 5A). Moreover, this effect was relatively independent of the presence of trimethoprim (Figure 5A). RpoS activates the transcription of a number of stress-responsive and stationary phase genes. Several of these genes are upregulated in an mgrB-knockout (Xu et al., 2019). Could the costs of mgrB deletion be related to unrequired expression of stress and stationary phase genes? We tested this hypothesis by generating knockouts of four RpoS-regulated genes, that is, gadX, cbpA, fic, and sra that are upregulated in a mgrB-deficient background (Xu et al., 2019). Knockouts of three of these, that is, gadX, cbpA, and fic were all selected over E. coli ΔmgrB independent of the presence of trimethoprim (Figure 5A). These data confirmed that loss of RpoS compensated for the fitness costs of PhoPQ derepression by ameliorating the expression of unwanted stationary-phase and stress-related genes (Figure 5B).
Strength of antibiotic selection determines the choice of adaptive strategy during the evolution of the PhoPQ-folA-RpoS transcriptional network
Our experiments so far identified two possible strategies available to E. coli for enhancing fitness subsequent to evolving tolerance, that is, resistance and cost compensation. We next asked what motivated the "choice" of strategy that an evolving bacterial population would adopt. Cost compensation is known to allow resistant bacteria to enhance their fitness in drug-free or low antibiotic conditions. Therefore, we argued that antibiotic concentration may determine which adaptive strategy was favored. In order to test this, we mixed tolerant (E. coli ΔmgrB), cost-compensated tolerant (E. coli ΔmgrBΔrpoS), and resistant (TMPR1) strains (initial mixing ratio 98:1:1) and allowed them to compete for ~9 generations in different concentrations of trimethoprim. In the absence of trimethoprim or in low concentrations of the antibiotic, ΔmgrBΔrpoS was strongly enriched over ΔmgrB and TMPR1 (Figure 6A). However, at higher drug concentrations, TMPR1 was favored over its competitors (Figure 6A). This clear dependency on trimethoprim concentration suggested that at low drug pressures, evolutionary outcomes were likely to be biased in favor of cost compensation rather than resistance-evolution.
We empirically tested the above prediction in three replicate lineages of E. coli K-12 MG1655 evolving under low trimethoprim pressure (50 ng/ml, 0.08×MIC, designated WTMP50 A, B, and C). In support of our hypothesis, resistant bacteria did not exceed more than 1% of the population in these lineages for the duration of the experiment (Figure 6B). Indeed, in WTMP50 lineage C, the fraction of resistant bacteria stayed below 0.001% of the total population throughout the experiment (Figure 6B). Gain in relative fitness of WTMP50 lineages over the ancestor was lower than WTMP300, and showed very poor concordance with enrichment of trimethoprim-resistant bacteria over time (Figure 6C). Importantly, like WTMP300, all three WTMP50 lineages got fitter in antibiotic-free medium over time (Figure 6C) and lost RpoS activity (Figure 6D). These data suggested that there were at least two parallel evolving adaptive trajectories in the WTMP50 lines, that is, resistant bacteria, that made up a minority of the population and cost-compensated, tolerant bacteria that made up the majority of the population. Genome sequencing of resistant and tolerant isolates from 147 to 350 generations under selection from WTMP50 lineage A, confirmed that this was indeed the case (Supplementary file 2, Figure 6E). The genetic trajectory of resistant bacteria in WTMP50 was very similar to the WTMP300, in that mutations in folA were acquired before cost-compensatory mutations in rpoS (Supplementary file 2, Figure 6E). On the other hand, tolerant bacteria acquired rpoS mutations earlier during evolution and did not acquire mutations in folA even after long-term antibiotic exposure (Supplementary file 2, Figure 6E). Thus, cost compensation was indeed the preferred strategy for adaptation at lower antibiotic concentrations.
Drug concentrations are generally equated with strength of selection for antibiotic resistance (Andersson and Hughes, 2014; Gullberg et al., 2011). However, selection strength is determined by intrinsic factors as well, such as the genetic background (Matange, 2020; Matange et al., 2019). Therefore, we tested whether the link between selection strength and evolutionary trajectory would hold true in another genetic background. We have shown earlier that Lon protease-deficient E. coli are trimethoprim-tolerant due to a folA-independent, acrAB efflux pump-dependent mechanism (Matange, 2020). For E. coli Δlon, it is expected that 300 ng/ml of trimethoprim would represent a milder selection strength than wild-type. Therefore, we established three lineages of E. coli Δlon that experienced trimethoprim (300 ng/ml, 0.5×MIC, designated LTMP300 A, B, and C) for 350 generations and examined their adaptation over time. At the phenotypic level, LTMP300 lineages resembled WTMP300 in terms of fixation of resistant-bacteria (Figure 7A). However, at the genotypic level there was substantial similarity between LTMP300 and WTMP50 lineages. We found no folA mutations in resistant isolates from the LTMP300 lineage A (Supplementary file 3, Figure 7B). Instead, trimethoprim-resistant bacteria derived from an early time point (147 generations) had already accumulated inactivating mutations in rpoS in addition to mutations in mgrB (Supplementary file 3, Figure 7B). Thus, while the LTMP300 lineages had the genetic signature of evolution at low selection strength, they phenotypically resembled high selection strength. This discrepancy was explained by the fact that Δlon and ΔmgrB showed an additive effect on trimethoprim IC50 (Figure 7C). While loss of mgrB alone resulted in drug tolerance, together with Lon-deficiency it resulted in drug resistance. Importantly, once mgrB mutations were acquired, LTMP300 lineages evolved by compensatory mutations rather than resistance-conferring mutations in folA, similar to WTMP50. These results, therefore, re-iterated the role of selection strength in determining the genetic trajectory followed by E. coli during adaptation to trimethoprim (Figure 7D).
Discussion
The PhoPQ signaling cascade is a well-studied two-component system, not just in Enterobacteria but also other species. Traditionally thought to be regulated by Mg2+ (Groisman, 2001), this system also responds to low pH (Eguchi et al., 2011) and periplasmic redox state (Lippa and Goulian, 2012). Mounting evidence in recent years has thrown light on its role in resistance to antibiotics, particularly antimicrobial peptides such as polymyxins (Mmatli et al., 2020). The clinical relevance of these findings was brought to the fore when a number of studies identified inactivating mutations in mgrB in carbapenem/colistin-resistant Klebsiella pneumoniae (Cannatelli et al., 2013; Haeili et al., 2017; Silva et al., 2021). Derepression of PhoPQ in Klebsiella results in alterations in the structure of Lipid A, conferring resistance to colistin (Kidd et al., 2017). Inactivation of mgrB also collaterally enhances virulence in Klebsiella (Kidd et al., 2017), a particularly worrying observation given that polymyxins represent the last line of defense against hospital-borne Gram-negative infections at the moment. Similarly in E. coli, C18G, another antimicrobial peptide, leads to growth arrest and filamentation by activating PhoPQ signaling (Yadavalli et al., 2016).
In the context of trimethoprim resistance, mutations in mgrB were reported earlier by Baym et al., 2016 from drug-resistant E. coli evolving on mega petri-plates (Baym et al., 2016), though a mechanistic basis for this change was not known. Transcriptional upregulation of folA upon loss of mgrB, a main result from our study, fills this lacuna. Additionally, by showing the synergistic interaction between mutations in mgrB and folA, our study also establishes how loss of mgrB facilitates resistance evolution. These results are in line with other studies that have identified the facilitatory role that drug tolerance plays in the evolution of resistance (Levin-Reisman et al., 2017).
Early inactivation of mgrB in E. coli upon trimethoprim exposure raises two questions. First, how does PhoPQ enhance the expression of folA? The only known transcriptional regulator of folA in E. coli is the TyrR transcription factor (Yang et al., 2007). However, we were unable to detect any change in the expression level of DHFR protein in a tyrR deletion strain (Figure 2—figure supplement 3). Given the widespread gene regulatory impact of PhoPQ derepression in E. coli, including its effects on other transcription factors, its effects on folA levels are likely to be secondary or tertiary. These indirect regulatory effects of PhoPQ derepression were sufficient to drive the selection of mgrB mutations. Further investigation into the mechanism of folA promoter upregulation in mgrB-deficient E. coli is likely to reveal additional regulators for this housekeeping gene. It is important to note that a larger number of loss-of-function mutations at mgrB are possible than gain-of-function mutations at the folA locus. This may also be an important factor explaining why mgrB mutations occurred earlier in our experiments than folA mutations.
Second, is inactivation of MgrB relevant in clinical or environmental strains of E. coli? We detected four kinds of mutations at the mgrB locus in our experiments, the most frequent being disruption of the promoter of mgrB. Additionally, we found three different single nucleotide insertions or deletions within the coding region. These mutations are expected to result in a longer or shorter protein with markedly altered C-terminal sequence as a result of translational frame-shift (Figure 7—figure supplement 1). We searched for protein sequences similar to these frame-shifted alleles of MgrB in the non-redundant protein sequence database. BLAST analysis revealed nine clinical and environmental E. coli strains from four different pathogen surveillance and sequencing projects that had altered C-terminal sequences similar to the mutant alleles found by us (Figure 7—figure supplement 1). Indeed, one of these strains, E. coli NCTC9075 from the collection of Public Health England, harbored exactly the same deletion in mgrB (+100 ΔA) that we identified in some of our laboratory-evolved strains (Figure 7—figure supplement 1). We also identified mutant MgrB sequences from three isolates of Salmonella enterica and one isolate of Citrobacter freundii (Figure 7—figure supplement 2). Thus, the relevance of MgrB in trimethoprim tolerance, or indeed in other phenotypes of pathogenic Enterobacteria, is likely to be far more significant that currently understood.
Though we have focussed on the compensatory effects of RpoS mutations in this study, genomics analyses revealed a number of other loci that accumulated mutations in our laboratory evolution experiments. For instance, mutations in the s. For instance, mutations in the deoA locus, predicted to cause premature termination of the DeoA protein, were detected in resistant isolates from all selection schemes employed by us (Table 1, Supplementary files 1–3). DeoA, a thymidine phosphorylase, is part of the pyrimidine catabolism pathway in E. coli (Razzell and Casshyap, 1964). Deactivation of this pathway is known to be beneficial for E. coli harboring hypomorphic DHFR alleles by reducing the demand for de novo pyrimidine biosynthesis (Rodrigues and Shakhnovich, 2019). Interestingly, mutations in some of the other loci from our study showed strict dependency on drug concentration used during selection. Most notably, mutations in baeS, ecpC, and dnaA genes, which code for two-component kinase, fimbrial chaperone, and replication initiator proteins respectively, were found in a majority of resistant isolates from the WTMP300 lineages, but not in WTMP50 lineages (Supplementary files 1 and 2). BaeS, sensor kinase of the BaeRS two-component system, has been associated with multi-drug resistance in a number of bacteria including E. coli due to its effect on the expression of efflux pumps (Nagakubo et al., 2002). The fimbrial chaperone EcpC, on the other hand, has never been associated with resistance to antibiotics, and is potentially a novel locus worth investigating. Point mutations in DnaA (D269G and D280G) that were detected in trimethoprim-resistant isolates from WTMP300 lineages are expected to compromise its ATP-binding (Erzberger et al., 2006). Indeed, D269 lies in the ‘Sensor I’ motif of DnaA and its mutation to Ala reduces ATPase activity of the protein in vitro (Kawakami et al., 2006). We envisage two possible reasons for these mutations in DnaA to be selected during adaptation to trimethoprim. First, thymine-less death, induced by thymine starvation or trimethoprim treatment, can be prevented in E. coli by lowered rates of replication initiation (Guzmán and Martín, 2015). Mutations in DnaA that compromise its ATPase activity would help achieve this. Second, DnaA mutations in clinical and laboratory strains of M. tuberculosis have been associated with isoniazid resistance due to altered gene expression profiles (Hicks et al., 2020), and a similar mechanism may be operative in the case of E. coli and trimethoprim. The appearance of mutations at these loci only at high antibiotic pressure is likely to be a result of the higher pressure for resistance.
We also observed mutations specific to genetic background. Most notable of these were mutations in 5 S and 16 S ribosomal RNAs, promoter of sulA and the ArcB transcriptional regulator in trimethoprim-resistant isolates evolved from the Δlon strain. The Lon protease degrades a number of housekeeping proteins, but can also degrade misfolded mutant proteins (Tsilibaris et al., 2006). We have shown in the past that a Lon-deficient background permits access to a larger repertoire of mutations to adapting bacterial populations (Matange, 2020). This phenomenon may very well explain why certain mutations were detected only in the Lon-deficient background. Further, mutations in the sulA gene promoter are expected to prevent its expression (Cole, 1983), thus alleviating the accumulation of this cell division inhibitor in the Lon mutant. Loss of sulA is advantageous for E. coli Δlon as it rescues the hyper-filamentation defect of this strain (Matange, 2020; Gottesman et al., 1981), which may explain why mutations in the sulA promoter were identified only in the LTMPR300 lineages.
To our knowledge, the present study is one of a kind in being able to trace the evolutionary trajectory of a gene regulatory network from its initial perturbation in response to selection, to eventual evolutionary rescue by compensatory mutations. Three main observations from this study significantly contribute to a general understanding of how gene regulatory networks change in response to selection. First, changes in gene regulatory networks can alter the fitness landscape of adaptive mutations. The evolutionary potential of gene regulatory changes is well-recognized particularly in the context of vertebrate development (Spitz and Duboule, 2008). Similarly, gene regulatory changes as mediators of evolution have been suggested in a series of synthetic biology experiments in prokaryotes (Bayer, 2010). Our study demonstrates that gene regulatory mutations can also act as facilitators of subsequent evolutionary change by shifting the fitness landscape of mutant alleles. This finding adds a new dimension to understanding the role of gene regulatory networks in mediating the evolution of new functions. Further, mutant alleles of enzymes, like in the case of DHFR, are known to often display trade-offs with stability (Matange et al., 2018; Rodrigues et al., 2016). The role of gene regulatory networks in negotiating these trade-offs needs to be factored in to understanding the genetics of evolutionary adaptation.
Second, cost compensation emerged as an important force driving the evolution of the PhoPQ-folA-RpoS network. Evolution of gene regulatory networks is thought to be primarily under stabilizing selection, in light of empirically observed properties such as developmental canalization and mutational robustness (Gilad et al., 2006). Compensatory evolution in regulatory networks has been understood mainly in the context of cis or trans-regulatory changes that minimize variation in gene expression (Signor and Nuzhdin, 2018; Thompson et al., 2015). Our results demonstrate that the fitness effects of adaptive mutations can be significant. As a result, compensation of fitness costs can also drive additional changes in regulatory networks. Importantly, loss of RpoS activity, a key cost-compensatory change in our experiments, is a commonly encountered trait in Enterobacteria (Dong et al., 2009; Kotewicz et al., 2003; Martinez-Garcia et al., 2003). RpoS mutants are associated with fitness advantage in aging bacterial colonies (Wrande et al., 2008), and natural isolates of E. coli often have inactive alleles of this sigma factor (Dong et al., 2009; Kotewicz et al., 2003; Martinez-Garcia et al., 2003; Chiang et al., 2011; Liu et al., 2017). The driving force behind loss of RpoS in Enterobacteria is thought to be selected for enhanced growth rate (Dong et al., 2009; Liu et al., 2017). Since RpoS competes with housekeeping sigma factor RpoD for RNA polymerase occupancy, loss of RpoS drives growth and replication, and compromises stress tolerance as a trade-off (Dong et al., 2009; Liu et al., 2017). The results of our study support this thesis. However, in our study, selection of RpoS mutants was driven in response to gene regulatory perturbation and not by environmental factors. Given the high prevalence of RpoS-mutations in natural isolates of Enterobacteria, it is very likely that they may serve as compensatory changes for other gene regulatory perturbations as well.
Finally, the dependency of evolutionary strategy on selection pressure is a key finding from this study (Figure 6D). Over the last decade, there has been significant attention given to how sub-lethal drug concentrations, such as those found in natural reservoirs of bacteria, alter the evolution of drug resistance (Matange et al., 2019; Wistrand-Yuen et al., 2018; Zhou et al., 2000; Andersson and Hughes, 2014; Baquero, 2001; Gutierrez et al., 2013; Hiltunen et al., 2017). The present study provides direct evidence linking drug concentration to the genetics of adaptation to antibiotics in bacteria. We also provide a mechanistic understanding for this link. Cost compensation by second-site mutations is well-known for antimicrobial resistance and is indeed a widely-accepted phenomenon (Björkman et al., 2000; Handel et al., 2006; Melnyk et al., 2015; Qi et al., 2016; Schulz zur Wiesch et al., 2010). It is thought to occur under relaxed selection, that is, when drug is not present in the environment. Our study demonstrates that compensation of cost and resistance evolution can also be thought of alternative strategies for enhancing fitness in the presence of the antibiotic. The relative gain in fitness from each of these strategies is determined not just by drug concentration but also by genetic background, and both these factors have vital influence on the evolutionary outcomes of bacteria under antibiotic selection.
Materials and methods
Strains, plasmids, and culture conditions
Request a detailed protocolE. coli K-12 MG1655 or its derivates were cultured in Luria-Bertani broth supplemented with appropriate antibiotics at 37°C with shaking (180–200 rpm). The various strains and plasmids used in this study and their sources are listed in Table 2. Single gene knockouts were obtained from the Keio collection (Baba et al., 2006; Yamamoto et al., 2009) (from the NBRP Resource, National Institute of Genetics, Japan) and moved into the MG1655 background using P1 transduction. For double gene knockouts, antibiotic marker cassette was first removed using Flp recombinase (pCP20-flp plasmid) (Datsenko and Wanner, 2000), and then subjected to P1 transduction.
Antibiotics used in this study were purchased from Himedia (India) or Sigma-Merck (USA, Germany). Kanamycin (30 µg/ml), chloramphenicol (25 µg/ml), or ampicillin (100 µg/ml) were added to liquid or semi-solid media just before inoculation as required. For isolating trimethoprim-resistant mutants from laboratory evolution experiments Luria-Bertani agar was supplemented with 1 µg/ml trimethoprim.
Antibiotic sensitivity
Request a detailed protocolTrimethoprim sensitivity of gene knockouts or evolved clones was tested using the following methods:
Minimum inhibitory concentration (MIC): MIC determination was performed using E-test. Trimethoprim strips (x–x µg/ml) were procured from Himedia (India) and used as per the manufacturer’s instruction.
Inhibitory concentration-50 (IC50): IC50 of trimethoprim was determined using a broth dilution assay as described in Matange et al., 2019 and Matange et al., 2018.
Colony formation assay: Colony-forming efficiency of wild-type or mutant E. coli was calculated by spotting 10 µl of neat or serially diluted (tenfold dilution series) late log-phase cultures on Luria-Bertani agar supplemented with 0, 50, 100, 300, 500, and 1000 ng/ml trimethoprim. Plates were incubated for 18 hr at 37°C before counting colonies.
Colonization efficiency: Wild-type or mutant E. coli cultures were grown to saturation and then serially diluted. Increasing bacterial CFUs (~104, 105,106, and 107) were added to 1 ml of Luria-Bertani broth supplemented with 0, 500, 1000, or 5000 ng/ml trimethoprim. Cultures were grown for 18 hr at 37°C and optical density (at 600 nm) was measured in a spectrophotometer.
Relative fitness and competitive growth assays
Request a detailed protocolRelative fitness of gene knockouts or evolved clones was estimated using competition assays. Test strains were mixed with ancestral E. coli from saturated mono-cultures (1:1, unless otherwise mentioned) in 1–3 ml Luria-Bertani broth supplemented with trimethoprim at the appropriate concentration. Strains were allowed to compete for ~9 generations.
For 2-strain or 3-strain competition, CFU/ml for each strain was determined by plating serially diluted cultures on antibiotic-free (total CFU) and kanamycin (knock-out strains) or trimethoprim (trimethoprim resistant isolates) supplemented Luria-Bertani agar. Relative fitness (w) or selection coefficient (s) was calculated using the following formulae:
where Tf and Ti are the final and initial CFU/ml of the test strain and Rf and Ri are the final and initial CFU/ml of the reference strain, respectively.
For competitive quotients, GFP-expressing derivative of E. coli wild-type (E. coli-GFP) was used as the reference strain. At the end of the competition experiment, 200 µl of each mixed culture was aliquoted into a 96-well plate and GFP fluorescence and optical density (at 600 nm) were estimated in a multi-mode plate reader (Varioskan-Thermo Scientific or Enlight-PerkinElmer). For each individual experiment, two controls were maintained, that is, E. coli -GFP alone and control competition between E. coli-GFP and wild-type. GFP fluorescence was normalized to optical density (nGFP) and used to calculate competitive quotient using the following formula:
where nGFP(c), nGFP(T), and nGFP(wt) are the normalized GFP fluorescence values of E. coli-GFP alone, mixed culture of E. coli -GFP, and test strain, and mixed culture of E. coli-GFP and wild-type, respectively.
Immunoblotting
Request a detailed protocolImmunoblotting was used to determine the expression level of DHFR in mutants or resistant isolates as described in Matange, 2020. FtsZ was used as a loading control. Anti-DHFR IgG was used at a concentration of 100 ng/ml. Anti-FtsZ (kind gift from Prof. Manjula Reddy, CCMB, India) polyclonal antiserum was used at a dilution of 1:50,000. Band intensities were quantitated in ImageJ software.
Activity of folA gene promoter
Request a detailed protocolThe activity of the promoter of folA was measured in wild-type or mutant E. coli using the pUA66-PfolA-GFP reporter plasmid (kind gift from Dr. Sanchari Bhattacharya, Department of Chemistry & Chemical Biology, Harvard University). Plasmid reporter was transformed into appropriate strains and GFP fluorescence was measured over growth in a multi-mode plate reader (Varioskan-Thermo Scientific). For comparison across strains, GFP fluorescence was normalized to optical density.
Purification of His-tagged PhoP
Request a detailed protocolHexa-His tagged PhoP was overexpressed in E. coli K-12 MG1655 and purified using Ni-NTA affinity chromatography. Bacterial strain harboring pCA24N-PhoP plasmid was grown in 50–100 ml Luria-Bertani broth till an optical density of ~0.8. Protein production was induced using IPTG (500 µM). Induction was carried out for 3 hr at 37°C, following which bacterial cultures were centrifuged at 3000×g for 10 min. Bacterial pellets were stored at –80°C until needed. Cells were lysed by sonication in Lysis buffer (Tris-HCl pH 8.0 100 mM, NaCl 100 mM, β-mercaptoethanol 1 mM, and glycerol 10%) and loaded onto a Ni-NTA-agarose resin (Bio-Rad, USA). The resin was washed using Wash buffer (Tris-HCl pH 8.0 100 mM, NaCl 300 mM, β-mercaptoethanol 1 mM, and imidazole 20 mM) and eluted using Elution buffer (Tris-HCl pH 8.0 100 mM, NaCl 100 mM, β-mercaptoethanol 1 mM, imidazole 300 mM, and glycerol 10%). Eluted protein was desalted using PD-10 desalting columns (Sigma-Merck, USA/Germany). Concentration of the desalted protein was estimated using Bradford’s method and the protein was stored at –80°C until needed.
Electrophoretic mobility shift assay
Request a detailed protocolFor electrophoretic mobility shift assay (EMSA), biotinylated probes shown in Figure 2—figure supplement 2 were synthesized (Sigma-Aldrich, USA) and annealed to a final concentration of 25 pmol/µl. They were further diluted to 25 fmol/µl in nuclease-free water. Recombinant purified PhoP protein was phosphorylated in vitro by incubating with 20 mM acetyl-phosphate and 10 mM MgCl2 for 45 min at 37°C. The phosphorylated protein was used immediately for further experiments. For binding and visualization, LightShift Chemilunescent EMSA Kit (Thermo Fisher Scientific, USA) was used as per the manufacturer’s instructions.
Laboratory evolution of trimethoprim resistance
Request a detailed protocolFor long-term laboratory evolution of trimethoprim resistance, replicate lineages (three for each drug concentration/genotype) were established using serial passaging in trimethoprim-supplemented Luria-Bertani broth (2 ml). Bacteria (1 %) were transferred to fresh media every 10–14 hr such that each growth cycle had between six and seven generations. Laboratory evolution was continued for ~350 generations (i.e., 50 growth cycles). Periodically an aliquot of culture from each lineage was mixed with an equal volume of sterile glycerol (50%) and frozen at –80°C for analysis.
The titer of trimethoprim-resistant bacteria from each of the frozen stocks was determined by spotting 10 µl of serially diluted culture on Luria-Bertani agar supplemented with trimethoprim (1 µg/ml). Resistant colonies (3–6) were selected at random for MIC determination or sequencing. Isolation of tolerant bacteria from WTMP50 lineages was done by spotting serially diluted cultures on drug-free medium. Tolerant isolates (3–6) were re-inoculated and their MIC was determined to confirm that they were not resistant to trimethoprim.
Congo red staining of bacterial colonies
Request a detailed protocolCongo red staining was used to determine RpoS activity in evolved lineages. Saturated cultures (10 µl) of appropriate strains were spotted on petri-plates containing YESCA medium (0.1% w/v yeast extract, 1% w/v casein hydrolysate, and 2% w/v agar) supplemented with Congo red (50 µg/ml) and Coomassie Brilliant Blue G-250 (10 µg/ml). Plates were incubated at 28°C to induce Curli fibers and photographed after ~24 hr of growth.
Sequencing and genome analysis
Request a detailed protocolGenomic DNA was extracted from evolved resistant/tolerant isolates using phenol:chloroform:isoamyl alcohol extraction as described in Matange et al., 2019. Extracted gDNA was then cleaned up using QiaAMP DNA purification spin columns (QIAGEN, Germany) and quantitated spectrophotometrically. Integrity of the prepared gDNA was verified by gel electrophoresis before further applications.
The folA locus (encompassing the coding region and 195 bp upstream promoter) was PCR amplified from gDNA of trimethoprim-resistant isolates as described in Matange et al., 2018 and Sanger sequenced. The mgrB and rpoS loci were similarly PCR amplified using the mgrB_seq_fwd (5′- CAACCTCTTCTCTTTTTATGTTCGC-3′) and mgrB_seq_rev (5′-CAACCAAAGACGCAATGTTCATCACG-3′) primers or rpoS_seq_fwd (5′-CGGAACCAGGCTTTTGCTTGAATG-3′) and rpoS_seq_rev (5′-GGCCAGCCTCGCTTGAGACTGGCC-3′) primers respectively and Sanger sequenced.
For whole-genome resequencing, paired-end whole-genome next-generation sequencing (NGS) was performed on a MiSeq system (Illumina, USA) with read lengths of 150 bps. Processed reads were aligned to the reference genome E. coli K-12 MG1655 (NZ_CP025268) using bowtie2. A cutoff of a minimum of 20× coverage was employed, which made >97% of the reads for each of the samples permissible for analysis. Variant calling and prediction of new junctions in the genome to identify large structural mutations, such as insertions and deletions, were done using Breseq using default settings (Deatherage and Barrick, 2014). All variants and new junctions that were already present in the wild-type ancestor were excluded.
Sanger sequencing services were provided by FirstBase (Malaysia). Library preparation and NGS services were provided by Bencos (India).
Data availability
Raw data from whole genome sequencing experiments used for this study have been deposited to GenBank under Project ID PRJNA741586.
-
NCBI BioProjectID PRJNA741586. Trimethoprim resistance in E. coli.
References
-
Microbiological effects of sublethal levels of antibioticsNature Reviews. Microbiology 12:465–478.https://doi.org/10.1038/nrmicro3270
-
Construction of Escherichia coli k-12 in-frame, single-gene knockout mutants: The Keio CollectionMolecular Systems Biology 2:2006.0008.https://doi.org/10.1038/msb4100050
-
Low-level antibacterial resistance: A gateway to clinical resistanceDrug Resistance Updates 4:93–105.https://doi.org/10.1054/drup.2001.0196
-
Regulation of bacterial virulence by two-component systemsCurrent Opinion in Microbiology 9:143–152.https://doi.org/10.1016/j.mib.2006.01.005
-
Bacterial histidine kinases as novel antibacterial drug targetsACS Chemical Biology 10:213–224.https://doi.org/10.1021/cb5007135
-
Phenotypic diversity caused by differential rpos activity among environmental Escherichia coli isolatesApplied and Environmental Microbiology 77:7915–7923.https://doi.org/10.1128/AEM.05274-11
-
Characterisation of the promoter for the LEXA regulated sula gene of Escherichia coliMolecular & General Genetics 189:400–404.https://doi.org/10.1007/BF00325901
-
Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using BRESEQMethods in Molecular Biology 1151:165–188.https://doi.org/10.1007/978-1-4939-0554-6_12
-
Regulation of acid resistance by connectors of two-component signal transduction systems in Escherichia coliJournal of Bacteriology 193:1222–1228.https://doi.org/10.1128/JB.01124-10
-
Structural basis for atp-dependent DNAA assembly and replication-origin remodelingNature Structural & Molecular Biology 13:676–683.https://doi.org/10.1038/nsmb1115
-
Massive overproduction of dihydrofolate reductase in bacteria as a response to the use of trimethoprimEuropean Journal of Biochemistry 162:473–476.https://doi.org/10.1111/j.1432-1033.1987.tb10664.x
-
Natural selection on gene expressionTrends in Genetics 22:456–461.https://doi.org/10.1016/j.tig.2006.06.002
-
Role of sulA and sulB in filamentation by lon mutants of Escherichia coli K-12Journal of Bacteriology 148:265–273.https://doi.org/10.1128/jb.148.1.265-273.1981
-
The pleiotropic two-component regulatory system PHOP-PHOQJournal of Bacteriology 183:1835–1842.https://doi.org/10.1128/JB.183.6.1835-1842.2001
-
Thymineless death, at the originFrontiers in Microbiology 6:499.https://doi.org/10.3389/fmicb.2015.00499
-
Mgrb alterations mediate colistin resistance in klebsiella pneumoniae isolates from IranFrontiers in Microbiology 8:2470.https://doi.org/10.3389/fmicb.2017.02470
-
The role of compensatory mutations in the emergence of drug resistancePLOS Computational Biology 2:e137.https://doi.org/10.1371/journal.pcbi.0020137
-
The two-component network and the general stress Sigma Factor rpos (Sigma s) in Escherichia coliAdvances in Experimental Medicine and Biology 631:40–53.https://doi.org/10.1007/978-0-387-78885-2_4
-
Molecular and evolutionary processes generating variation in gene expressionNature Reviews. Genetics 22:203–215.https://doi.org/10.1038/s41576-020-00304-w
-
Role of insertion elements and yycfg in the development of decreased susceptibility to vancomycin in Staphylococcus aureusInternational Journal of Medical Microbiology 297:205–215.https://doi.org/10.1016/j.ijmm.2007.02.002
-
Recent insights into the genotype-phenotype relationship from massively parallel genetic assaysEvolutionary Applications 12:1721–1742.https://doi.org/10.1111/eva.12846
-
Perturbation of the oxidizing environment of the periplasm stimulates the phoq/phop system in Escherichia coliJournal of Bacteriology 194:1457–1463.https://doi.org/10.1128/JB.06055-11
-
Polymorphism in the Yclc-rpos region of EnterobacteriaCurrent Microbiology 46:365–370.https://doi.org/10.1007/s00284-002-3814-6
-
Highly contingent phenotypes of Lon protease deficiency in Escherichia coli upon antibiotic challengeJournal of Bacteriology 202:e00561-19.https://doi.org/10.1128/JB.00561-19
-
The fitness costs of antibiotic resistance mutationsEvolutionary Applications 8:273–283.https://doi.org/10.1111/eva.12196
-
Identification and molecular characterization of the mg2+ stimulon of Escherichia coliJournal of Bacteriology 185:3696–3702.https://doi.org/10.1128/JB.185.13.3696-3702.2003
-
Signal integration in bacterial two-component regulatory systemsGenes & Development 22:2601–2611.https://doi.org/10.1101/gad.1700308
-
Comparison of the PhoPQ regulon in Escherichia coli and Salmonella typhimuriumJournal of Molecular Evolution 60:462–474.https://doi.org/10.1007/s00239-004-0212-7
-
Phop-mediated repression of the spi1 type 3 secretion system in Salmonella enterica SEROVAR typhimuriumJournal of Bacteriology 201:e00264-19.https://doi.org/10.1128/JB.00264-19
-
Signalling pathways in two-component phosphorelay systemsTrends in Microbiology 7:115–120.https://doi.org/10.1016/s0966-842x(99)01458-4
-
Czcr-czcs, a two-component system involved in heavy metal and carbapenem resistance in Pseudomonas aeruginosaThe Journal of Biological Chemistry 279:8761–8768.https://doi.org/10.1074/jbc.M312080200
-
The genomic basis of adaptation to the fitness cost of rifampicin resistance in Pseudomonas aeruginosa proceedingsBiological Sciences 283:1822.
-
Substrate Specificity and Induction of Thymidine Phosphorylase in Escherichia coliThe Journal of Biological Chemistry 239:1789–1793.https://doi.org/10.1016/S0021-9258(18)91259-7
-
Rewiring bacterial two-component systems by modular dna-binding domain swappingNature Chemical Biology 15:690–698.https://doi.org/10.1038/s41589-019-0286-6
-
Compensation of fitness costs and reversibility of antibiotic resistance mutationsAntimicrobial Agents and Chemotherapy 54:2085–2095.https://doi.org/10.1128/AAC.01460-09
-
Bioinformatics and experimental analysis of proteins of two-component systems in Myxococcus xanthusJournal of Bacteriology 190:613–624.https://doi.org/10.1128/JB.01502-07
-
The Evolution of Gene Expression in cis and transTrends in Genetics 34:532–544.https://doi.org/10.1016/j.tig.2018.03.007
-
Insertion sequences disrupting mgrB in carbapenem-resistant Klebsiella pneumoniae strains in BrazilJournal of Global Antimicrobial Resistance 24:53–57.https://doi.org/10.1016/j.jgar.2020.11.003
-
Differential gene regulatory networks in development and diseaseCellular and Molecular Life Sciences 75:1013–1025.https://doi.org/10.1007/s00018-017-2679-6
-
Comparative analysis of gene regulatory networks: From network reconstruction to evolutionAnnual Review of Cell and Developmental Biology 31:399–428.https://doi.org/10.1146/annurev-cellbio-100913-012908
-
Roles of two-component regulatory systems in antibiotic resistanceFuture Microbiology 14:533–552.https://doi.org/10.2217/fmb-2019-0002
-
Biological roles of the Lon atp-dependent proteaseResearch in Microbiology 157:701–713.https://doi.org/10.1016/j.resmic.2006.05.004
-
Diverse fluoroquinolone resistance plasmids from retail meat E. coli in the United StatesFrontiers in Microbiology 10:2826.https://doi.org/10.3389/fmicb.2019.02826
-
Small RNAs controlled by two-component systemsAdvances in Experimental Medicine and Biology 631:54–79.https://doi.org/10.1007/978-0-387-78885-2_5
-
Directed evolution of trimethoprim resistance in Escherichia coliThe FEBS Journal 274:2661–2671.https://doi.org/10.1111/j.1742-4658.2007.05801.x
-
MGRB affects the acid stress response of Escherichia coli by modulating the expression of iramFEMS Microbiology Letters 366:fnz123.https://doi.org/10.1093/femsle/fnz123
-
Functional determinants of a small protein controlling a broadly conserved bacterial sensor kinaseJournal of Bacteriology 10:JB.00305-20.https://doi.org/10.1128/JB.00305-20
-
Update on the Keio collection of Escherichia coli single-gene deletion mutantsMolecular Systems Biology 5:335.https://doi.org/10.1038/msb.2009.92
-
FOLA, a new member of the tyrr regulon in Escherichia coli k-12Journal of Bacteriology 189:6080–6084.https://doi.org/10.1128/JB.00482-07
-
Selection of antibiotic-resistant bacterial mutants: Allelic diversity among fluoroquinolone-resistant mutationsThe Journal of Infectious Diseases 182:517–525.https://doi.org/10.1086/315708
Article and author information
Author details
Funding
Department of Science and Technology, Ministry of Science and Technology, India
- Nishad Matange
Science and Engineering Research Board
- Nishad Matange
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors would like to acknowledge Dr. Amrita Hazra (IISER, Pune, India) and Dr. Manjula Reddy (CCMB, Hyderabad, India) for strains and reagents. The authors would like to acknowledge National Bioresource, NIG, Japan, for E. coli mutants. This project was funded by the INSPIRE program (Department of Science and Technology, Govt. of India) and IISER-Pune. NM is a recipient of the SERB-Research Scientist Award (Science Education and Research Board, Govt. of India) and the DBT-Wellcome India Alliance Intermediate Fellowship.
Copyright
© 2021, Patel and Matange
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,747
- views
-
- 380
- downloads
-
- 19
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling due to steric similarities to endogenous hormones. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) and differentially expressed genes (DEGs) that predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many EDCs persists over multiple generations. Many studies have described direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDC-induced epimutations are introduced or subsequently propagated, whether there are cell type-specific susceptibilities to the same EDC, and whether this correlates with differential expression of relevant hormone receptors. We exposed cultured pluripotent (iPS), somatic (Sertoli and granulosa), and primordial germ cell-like (PGCLC) cells to BPS and found that differential incidences of BPS-induced epimutations and DEGs correlated with differential expression of relevant hormone receptors inducing epimutations near relevant hormone response elements in somatic and pluripotent, but not germ cell types. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and DEGs was largely retained, however, >90% of the specific epimutations and DEGs were replaced by novel epimutations and DEGs. These results suggest a unique mechanism by which an EDC-induced epimutated state may be propagated transgenerationally.
-
- Genetics and Genomics
Microbial secondary metabolites are a rich source for pharmaceutical discoveries and play crucial ecological functions. While tools exist to identify secondary metabolite clusters in genomes, precise sequence-to-function mapping remains challenging because neither function nor substrate specificity of biosynthesis enzymes can accurately be predicted. Here, we developed a knowledge-guided bioinformatic pipeline to solve these issues. We analyzed 1928 genomes of Pseudomonas bacteria and focused on iron-scavenging pyoverdines as model metabolites. Our pipeline predicted 188 chemically different pyoverdines with nearly 100% structural accuracy and the presence of 94 distinct receptor groups required for the uptake of iron-loaded pyoverdines. Our pipeline unveils an enormous yet overlooked diversity of siderophores (151 new structures) and receptors (91 new groups). Our approach, combining feature sequence with phylogenetic approaches, is extendable to other metabolites and microbial genera, and thus emerges as powerful tool to reconstruct bacterial secondary metabolism pathways based on sequence data.