Decoding genome-wide GadEWX-transcriptional regulatory networks reveals multifaceted cellular responses to acid stress in Escherichia coli

An acid-tolerance response system protecting exponentially growing Escherichia coli

The ability to grow at moderate acidic conditions (pH 4.0–5.0) is important to Escherichia coli colonization of the host’s intestine. Several regulatory systems are known to control acid resistance in E. coli, enabling the bacteria to survive under acidic conditions without growth. Here, we characterize an acid-tolerance response (ATR) system and its regulatory circuit, required for E. coli exponential growth at pH 4.2. A two-component system CpxRA directly senses acidification through protonation of CpxA periplasmic histidine residues, and upregulates the fabA and fabB genes, leading to increased production of unsaturated fatty acids. Changes in lipid composition decrease membrane fluidity, F₀F₁-ATPase activity, and improve intracellular pH homeostasis. The ATR system is important for E. coli survival in the mouse intestine and for production of higher level of 3-hydroxypropionate during fermentation. Furthermore, this ATR system appears to be conserved in other Gram-negative bacteria.

The ability to grow at acidic pH is crucial for E. coli colonization of the host’s intestine. Here, the authors identify an acid-tolerance response system that is important for E. coli exponential growth at pH 4.2, survival in the mouse intestine, and production of 3-hydroxypropionate during fermentation.

Effect of DR1558, a Deinococcus radiodurans response regulator, on the production of GABA in the recombinant Escherichia coli under low pH conditions

BACKGROUND

Gamma aminobutyric acid (GABA) is an important platform chemical, which has been used as a food additive and drug. Additionally, GABA is a precursor of 2-pyrrolidone, which is used in nylon synthesis. GABA is usually synthesized from glutamate in a reaction catalyzed by glutamate decarboxylase (GAD). Currently, there are several reports on GABA production from monosodium glutamate (MSG) or glucose using engineered microbes. However, the optimal pH for GAD activity is 4, which is the limiting factor for the efficient microbial fermentative production of GABA as fermentations are performed at pH 7. Recently, DR1558, a response regulator in the two-component signal transduction system was identified in Deinococcus radiodurans. DR1558 is reported to confer cellular robustness to cells by binding the promoter regions of genes via DNA-binding domains or by binding to the effector molecules, which enable the microorganisms to survive in various environmental stress conditions, such as oxidative stress, high osmotic shock, and low pH.

RESULTS

In this study, the effect of DR1558 in enhancing GABA production was examined using two different strategies: whole-cell bioconversion of GABA from MSG and direct fermentative production of GABA from glucose under acidic culture conditions. In the whole-cell bioconversion, GABA produced by E. coli expressing GadBC and DR1558 (6.52 g/L GABA from 13 g/L MSG·H₂O) in shake flask culture at pH 4.5 was 2.2-fold higher than that by E. coli expressing only GadBC (2.97 g/L of GABA from 13 g/L MSG·H₂O). In direct fermentative production of GABA from glucose, E. coli ∆gabT expressing isocitrate dehydrogenase (IcdA), glutamate dehydrogenase (GdhA), GadBC, and DR1558 produced 1.7-fold higher GABA (2.8 g/L of GABA from 30 g/L glucose) than E. coli ∆gabT expressing IcdA, GdhA, and GadBC (1.6 g/L of GABA from 30 g/L glucose) in shake flask culture at an initial pH 7.0. The transcriptional analysis of E. coli revealed that DR1558 conferred acid resistance to E. coli during GABA production. The fed-batch fermentation of E. coli expressing IcdA, GdhA, GadBC, and DR1558 performed at pH 5.0 resulted in the final GABA titer of 6.16 g/L by consuming 116.82 g/L of glucose in 38 h.

CONCLUSION

This is the first report to demonstrate GABA production by acidic fermentation and to provide an engineering strategy for conferring acid resistance to the recombinant E. coli for GABA production.

Genome Analysis and Physiological Characterization of Four Streptococcus thermophilus Strains Isolated From Chinese Traditional Fermented Milk

Streptococcus thermophilus plays important roles in the dairy industry and is widely used as a dairy starter in the production of fermented dairy products. The genomes of S. thermophilus strains CS5, CS9, CS18, and CS20 from fermented milk in China were sequenced and used for biodiversity analysis. In the present study, the phylogenetic analysis of all 34 S. thermophilus genomes publicly available including these four strains reveals that the phylogenetic reconstruction does not match geographic distribution as strains isolated from the same continent are not even clustered on the nearby branches. The core and variable genes were also identified, which vary among strains from 0 to 202. CS9 strain contained 127 unique genes from a variety of distantly related species. It was speculated that CS9 had undergone horizontal gene transfer (HGT) during the long evolutionary process. The safety evaluation of these four strains indicated that none of them contains antibiotic resistance genes and that they are all sensitive to multiple antibiotics. In addition, the strains do not contain any pathogenic virulence factors or plasmids and thus can be considered safe. Furthermore, these strains were investigated in terms of their technological properties including milk acidification, exopolysaccharide (EPS) and γ-aminobutyric acid (GABA) production, and in vitro survival capacity in the gastrointestinal tract. CS9 possesses a special eps gene cluster containing significant traces of HGT, while the eps gene clusters of CS5, CS18, and CS20 are almost the same. The monosaccharide compositional analysis indicated that crude EPS-CS5, EPS-CS9, EPS-CS18, and EPS-CS20 contain similar monosaccharide compositions with different ratios. Furthermore, CS9 was one of a few GABA-producing strains that could ferment glutamate to produce GABA, which is beneficial for improving the acid tolerance of the strain. CS18 has the most potential for the production of fermented food among these four strains because of its fast growth rate, rapid acidifying capacity, and stronger acid and bile salt resistance capacity. This study focused on the genome analysis of the four new S. thermophilus strains to investigate the diversity of strains and provides a reference for selecting excellent strains by use of the genome data.

The role of bacterial cell envelope structures in acid stress resistance in E. coli

Acid resistance (AR) is an indispensable mechanism for the survival of neutralophilic bacteria, such as Escherichia coli (E. coli) strains that survive in the gastrointestinal tract. E. coli acid tolerance has been extensively studied during past decades, with most studies focused on gene regulation and mechanisms. However, the role of cell membrane structure in the context of acid stress resistance has not been discussed in depth. Here, we provide a comprehensive review of the roles and mechanisms of the E. coli cell envelope from different membrane components, such as membrane proteins, fatty acids, chaperones, and proton-consuming systems, and particularly focus on the innovative effects revealed by recent studies. We hope that the information guides us to understand the bacterial survival strategies under acid stress and to further explore the AR regulatory mechanisms to prevent or treat E. coli and other related Gram-negative bacteria infection, or to enhance the AR of engineering E. coli.

Plastic Circuits: Regulatory Flexibility in Fine Tuning Pathogen Success

Bacterial pathogens employ diverse fitness and virulence mechanisms to gain an advantage in competitive niches. These lifestyle-specific traits require integration into the regulatory network of the cell and are often controlled by pre-existing transcription factors. In this review, we highlight recent advances that have been made in characterizing this regulatory flexibility in prominent members of the Enterobacteriaceae. We focus on the direct global interactions between transcription factors and their target genes in pathogenic Escherichia coli and Salmonella revealed using chromatin immunoprecipitation coupled with next-generation sequencing. Furthermore, the implications and advantages of such regulatory adaptations in benefiting distinct pathogenic lifestyles are discussed.

Limits to a classic paradigm: most transcription factors in E. coli regulate genes involved in multiple biological processes

Transcription factors (TFs) are important drivers of cellular decision-making. When bacteria encounter a change in the environment, TFs alter the expression of a defined set of genes in order to adequately respond. It is commonly assumed that genes regulated by the same TF are involved in the same biological process. Examples of this are methods that rely on coregulation to infer function of not-yet-annotated genes. We have previously shown that only 21% of TFs involved in metabolism regulate functionally homogeneous genes, based on the proximity of the gene products’ catalyzed reactions in the metabolic network. Here, we provide more evidence to support the claim that a 1-TF/1-process relationship is not a general property. We show that the observed functional heterogeneity of regulons is not a result of the quality of the annotation of regulatory interactions, nor the absence of protein–metabolite interactions, and that it is also present when function is defined by Gene Ontology terms. Furthermore, the observed functional heterogeneity is different from the one expected by chance, supporting the notion that it is a biological property. To further explore the relationship between transcriptional regulation and metabolism, we analyzed five other types of regulatory groups and identified complex regulons (i.e. genes regulated by the same combination of TFs) as the most functionally homogeneous, and this is supported by coexpression data. Whether higher levels of related functions exist beyond metabolism and current functional annotations remains an open question.

RstA, a two-component response regulator, plays important roles in multiple virulence-associated processes in enterohemorrhagic Escherichia coli O157:H7

Background

Enterohemorrhagic Escherichia coli O157:H7 (EHEC O157) causes bloody diarrhea and hemolytic-uremic syndrome. EHEC O157 encounters varied microenvironments during infection, and can efficiently adapt to these using the two-component system (TCS). Recently, a functional TCS, RstAB, has been implicated in the regulation of virulence of several bacterial pathogens. However, the regulatory function of RstAB in EHEC O157 is poorly understood. This study aimed at providing insights into the global effects of RstA on gene expression in EHEC O157.

Results

In the present study, we analyzed gene expression differences between the EHEC O157 wild-type strain and a ΔrstA mutant using RNA-seq technology. Genes with differential expression in the ΔrstA mutant compared to that in the wild-type strain were identified and grouped into clusters of orthologous categories. RstA promoted EHEC O157 LEE gene expression, adhesion in vitro, and colonization in vivo by indirect regulation. We also found that RstA could bind directly to the promoter region of hdeA and yeaI to enhance acid tolerance and decrease biofilm formation by modulating the concentration of c-di-GMP.

Conclusions

In summary, the RstAB TCS in EHEC O157 plays a major role in the regulation of virulence, acid tolerance, and biofilm formation. We clarified the regulatory function of RstA, providing an insight into mechanisms that may be potential drug targets for treatment of EHEC O157-related infections.

Inverted Regulation of Multidrug Efflux Pumps, Acid Resistance, and Porins in Benzoate-Evolved Escherichia coli K-12

Benzoic acid, a partial uncoupler of the proton motive force (PMF), selects for sensitivity to chloramphenicol and tetracycline during the experimental evolution of Escherichia coli K-12. Transcriptomes of E. coli isolates evolved with benzoate showed the reversal of benzoate-dependent regulation, including the downregulation of multidrug efflux pump genes, the gene for the Gad acid resistance regulon, the nitrate reductase genes narHJ, and the gene for the acid-consuming hydrogenase Hyd-3. However, the benzoate-evolved strains had increased expression of OmpF and other large-hole porins that admit fermentable substrates and antibiotics. Candidate genes identified from benzoate-evolved strains were tested for their roles in benzoate tolerance and in chloramphenicol sensitivity. Benzoate or salicylate tolerance was increased by deletion of the Gad activator ariR or of the acid fitness island from slp to the end of the gadX gene encoding Gad regulators and the multidrug pump genes mdtEF Benzoate tolerance was also increased by deletion of multidrug component gene emrA, RpoS posttranscriptional regulator gene cspC, adenosine deaminase gene add, hydrogenase gene hyc (Hyd-3), and the RNA chaperone/DNA-binding regulator gene hfq Chloramphenicol resistance was decreased by mutations in genes for global regulators, such as RNA polymerase alpha subunit gene rpoA, the Mar activator gene rob, and hfq Deletion of lipopolysaccharide biosynthetic kinase gene rfaY decreased the rate of growth in chloramphenicol. Isolates from experimental evolution with benzoate had many mutations affecting aromatic biosynthesis and catabolism, such as aroF (encoding tyrosine biosynthesis) and apt (encoding adenine phosphoribosyltransferase). Overall, benzoate or salicylate exposure selects for the loss of multidrug efflux pumps and of hydrogenases that generate a futile cycle of PMF and upregulates porins that admit fermentable nutrients and antibiotics.IMPORTANCE Benzoic acid is a common food preservative, and salicylic acid (2-hydroxybenzoic acid) is the active form of aspirin. At high concentrations, benzoic acid conducts a proton across the membrane, depleting the proton motive force. In the absence of antibiotics, benzoate exposure selects against proton-driven multidrug efflux pumps and upregulates porins that admit fermentable substrates but that also allow the entry of antibiotics. Thus, evolution with benzoate and related molecules, such as salicylates, requires a trade-off for antibiotic sensitivity, a trade-off that could help define a stable gut microbiome. Benzoate and salicylate are naturally occurring plant signal molecules that may modulate the microbiomes of plants and animal digestive tracts so as to favor fermenters and exclude drug-resistant pathogens.

Experimental Evolution of Escherichia coli K-12 in the Presence of Proton Motive Force (PMF) Uncoupler Carbonyl Cyanide m-Chlorophenylhydrazone Selects for Mutations Affecting PMF-Driven Drug Efflux Pumps

Experimental evolution of Escherichia coli K-12 with benzoate, a partial uncoupler of the proton motive force (PMF), selects for mutations that decrease antibiotic resistance. We conducted experimental evolution in the presence of carbonyl cyanide m-chlorophenylhydrazone (CCCP), a strong uncoupler. Cultures were serially diluted daily 1:100 in LBK medium containing 20 to 150 µM CCCP buffered at pH 6.5 or at pH 8.0. After 1,000 generations, the populations tolerated up to 150 µM CCCP. Sequenced isolates had mutations in mprA (emrR), which downregulates the EmrAB-TolC pump that exports CCCP. A mprA::kanR deletion conferred growth at 60 μM CCCP, though not at the higher levels resisted by evolved strains (150 µM). Some mprA mutant strains also had point mutations affecting emrA, but deletion of emrA abolished the CCCP resistance. Thus, CCCP-evolved isolates contained additional adaptations. One isolate lacked emrA or mprA mutations but had mutations in cecR (ybiH), whose product upregulates drug pumps YbhG and YbhFSR, and in gadE, which upregulates the multidrug pump MdtEF. A cecR::kanR deletion conferred partial resistance to CCCP. Other multidrug efflux genes that had mutations included ybhR and acrAB The acrB isolate was sensitive to the AcrAB substrates chloramphenicol and tetracycline. Other mutant genes in CCCP-evolved strains include rng (RNase G) and cyaA (adenylate cyclase). Overall, experimental evolution revealed a CCCP-dependent fitness advantage for mutations increasing CCCP efflux via EmrA and for mutations that may deactivate proton-driven pumps for drugs not present (cecR, gadE, acrAB, and ybhR). These results are consistent with our previous report of drug sensitivity associated with evolved benzoate tolerance.IMPORTANCE The genetic responses of bacteria to depletion of proton motive force (PMF), and their effects on drug resistance, are poorly understood. PMF drives export of many antibiotics, but the energy cost may decrease fitness when antibiotics are absent. Our evolution experiment reveals genetic mechanisms of adaptation to the PMF uncoupler CCCP, including selection for increased CCCP efflux but also against the expression of PMF-driven pumps for drugs not present. The results have implications for our understanding of the gut microbiome, which experiences high levels of organic acids that decrease PMF.

Systematic discovery of uncharacterized transcription factors in Escherichia coli K-12 MG1655

Transcriptional regulation enables cells to respond to environmental changes. Of the estimated 304 candidate transcription factors (TFs) in Escherichia coli K-12 MG1655, 185 have been experimentally identified, but ChIP methods have been used to fully characterize only a few dozen. Identifying these remaining TFs is key to improving our knowledge of the E. coli transcriptional regulatory network (TRN). Here, we developed an integrated workflow for the computational prediction and comprehensive experimental validation of TFs using a suite of genome-wide experiments. We applied this workflow to (i) identify 16 candidate TFs from over a hundred uncharacterized genes; (ii) capture a total of 255 DNA binding peaks for ten candidate TFs resulting in six high-confidence binding motifs; (iii) reconstruct the regulons of these ten TFs by determining gene expression changes upon deletion of each TF and (iv) identify the regulatory roles of three TFs (YiaJ, YdcI, and YeiE) as regulators of l-ascorbate utilization, proton transfer and acetate metabolism, and iron homeostasis under iron-limited conditions, respectively. Together, these results demonstrate how this workflow can be used to discover, characterize, and elucidate regulatory functions of uncharacterized TFs in parallel.

Modeling the multi-scale mechanisms of macromolecular resource allocation

As microbes face changing environments, they dynamically allocate macromolecular resources to produce a particular phenotypic state. Broad 'omics' data sets have revealed several interesting phenomena regarding how the proteome is allocated under differing conditions, but the functional consequences of these states and how they are achieved remain open questions. Various types of multi-scale mathematical models have been used to elucidate the genetic basis for systems-level adaptations. In this review, we outline several different strategies by which microbes accomplish resource allocation and detail how mathematical models have aided in our understanding of these processes. Ultimately, such modeling efforts have helped elucidate the principles of proteome allocation and hold promise for further discovery.

Deciphering bacterial xylose metabolism and metabolic engineering of industrial microorganisms for use as efficient microbial cell factories

The goal of sustainable production of biochemicals and biofuels has driven the engineering of microbial cell as factories that convert low-value substrates to high-value products. Xylose is the second most abundant sugar substrate in lignocellulosic hydrolysates. We analyzed the mechanisms of xylose metabolism using genome sequencing data of 492 industrially relevant bacterial species in the mini-review. The analysis revealed the xylose isomerase and Weimberg pathways as the major routes across diverse routes of bacterial xylose metabolism. In addition, we discuss recent developments in metabolic engineering of xylose metabolism in industrial microorganisms. Genome-scale analyses have revealed xylose pathway-specific flux landscapes. Overall, a comprehensive understanding of bacterial xylose metabolism could be useful for the feasible development of microbial cell factories.

Evolution of gene knockout strains of E. coli reveal regulatory architectures governed by metabolism

Biological regulatory network architectures are multi-scale in their function and can adaptively acquire new functions. Gene knockout (KO) experiments provide an established experimental approach not just for studying gene function, but also for unraveling regulatory networks in which a gene and its gene product are involved. Here we study the regulatory architecture of Escherichia coli K-12 MG1655 by applying adaptive laboratory evolution (ALE) to metabolic gene KO strains. Multi-omic analysis reveal a common overall schema describing the process of adaptation whereby perturbations in metabolite concentrations lead regulatory networks to produce suboptimal states, whose function is subsequently altered and re-optimized through acquisition of mutations during ALE. These results indicate that metabolite levels, through metabolite-transcription factor interactions, have a dominant role in determining the function of a multi-scale regulatory architecture that has been molded by evolution.

The function of metabolic genes in the context of regulatory networks is not well understood. Here, the authors investigate the adaptive responses of E. coli after knockout of metabolic genes and highlight the influence of metabolite levels in the evolution of regulatory function.

Systems assessment of transcriptional regulation on central carbon metabolism by Cra and CRP

Two major transcriptional regulators of carbon metabolism in bacteria are Cra and CRP. CRP is considered to be the main mediator of catabolite repression. Unlike for CRP, in vivo DNA binding information of Cra is scarce. Here we generate and integrate ChIP-exo and RNA-seq data to identify 39 binding sites for Cra and 97 regulon genes that are regulated by Cra in Escherichia coli. An integrated metabolic-regulatory network was formed by including experimentally-derived regulatory information and a genome-scale metabolic network reconstruction. Applying analysis methods of systems biology to this integrated network showed that Cra enables optimal bacterial growth on poor carbon sources by redirecting and repressing glycolysis flux, by activating the glyoxylate shunt pathway, and by activating the respiratory pathway. In these regulatory mechanisms, the overriding regulatory activity of Cra over CRP is fundamental. Thus, elucidation of interacting transcriptional regulation of core carbon metabolism in bacteria by two key transcription factors was possible by combining genome-wide experimental measurement and simulation with a genome-scale metabolic model.

Functional Prediction of Hypothetical Transcription Factors of Escherichia coli K-12 Based on Expression Data

The repertoire of 304 DNA-binding transcription factors (TFs) in Escherichia coli K-12 has been described recently, with 196 TFs experimentally characterized and 108 proteins predicted by sequence comparisons. Based on 303 expression profile patterns retrieved from the Colombos database 12 clusters were identified, including hypothetical and experimentally characterized TFs, using a spectral clustering algorithm based on a 3NN graph built using 14 principal components that represent 65% of the variance of the expression data. In a posterior step, clusters were characterized in terms of their associated overrepresented functions, based on KEGG, Supfam annotations and Pfam assignments among other functional categories using an enrichment test, reinforcing the notion that the identified clusters are functionally similar among them. Based on these data, the we identified 12 clusters in which hypothetical and known TFs share similar regulatory and physiological functions, such as module associations of toxin-antitoxin (TA) systems with DNA repair mechanisms, amino acid biosynthesis, and carbon metabolism/transport, among others. This analysis has increased our knowledge about gene regulation in E. coli K-12 and can be further expanded to other organisms.

Experimental evolution of diverse Escherichia coli metabolic mutants identifies genetic loci for convergent adaptation of growth rate

Cell growth is determined by substrate availability and the cell’s metabolic capacity to assimilate substrates into building blocks. Metabolic genes that determine growth rate may interact synergistically or antagonistically, and can accelerate or slow growth, depending on genetic background and environmental conditions. We evolved a diverse set of Escherichia coli single-gene deletion mutants with a spectrum of growth rates and identified mutations that generally increase growth rate. Despite the metabolic differences between parent strains, mutations that enhanced growth largely mapped to core transcription machinery, including the β and β’ subunits of RNA polymerase (RNAP) and the transcription elongation factor, NusA. The structural segments of RNAP that determine enhanced growth have been previously implicated in antibiotic resistance and in the control of transcription elongation and pausing. We further developed a computational framework to characterize how the transcriptional changes that occur upon acquisition of these mutations affect growth rate across strains. Our experimental and computational results provide evidence for cases in which RNAP mutations shift the competitive balance between active transcription and gene silencing. This study demonstrates that mutations in specific regions of RNAP are a convergent adaptive solution that can enhance the growth rate of cells from distinct metabolic states.

The loss of a metabolic function caused by gene deletion can be compensated, in certain cases, by the concurrent mutation of a second gene. Whether such gene pairs share a local chemical or regulatory relationship or interact via a non-local mechanism has implications for the co-evolution of genetic changes, development of alternatives to gene therapy, and the design of combination antimicrobial therapies that select against resistance. Yet, we lack a comprehensive knowledge of adaptive responses to metabolic mutations, and our understanding of the mechanisms underlying genetic rescue remains limited. We present results of a laboratory evolution approach that has the potential to address both challenges, showing that mutations in specific regions of RNA polymerase enhance growth rates of distinct mutant strains of Escherichia coli with a spectrum of growth defects. Several of these adaptive mutations are deleterious when engineered directly into the original wild-type strain under alternative cultivation conditions, and thus have epistatic rescue properties when paired with the corresponding primary metabolic gene deletions. Our combination of adaptive evolution, directed genetic engineering, and mathematical analysis of transcription and growth rate distinguishes between rescue interactions that are specific or non-specific to a particular deletion. Our study further supports a model for RNA polymerase as a locus of convergent adaptive evolution from different sub-optimal metabolic starting points.

Global transcriptional regulatory network for Escherichia coli robustly connects gene expression to transcription factor activities

Transcriptional regulatory networks (TRNs) have been studied intensely for >25 y. Yet, even for the Escherichia coli TRN-probably the best characterized TRN-several questions remain. Here, we address three questions: (i) How complete is our knowledge of the E. coli TRN; (ii) how well can we predict gene expression using this TRN; and (iii) how robust is our understanding of the TRN? First, we reconstructed a high-confidence TRN (hiTRN) consisting of 147 transcription factors (TFs) regulating 1,538 transcription units (TUs) encoding 1,764 genes. The 3,797 high-confidence regulatory interactions were collected from published, validated chromatin immunoprecipitation (ChIP) data and RegulonDB. For 21 different TF knockouts, up to 63% of the differentially expressed genes in the hiTRN were traced to the knocked-out TF through regulatory cascades. Second, we trained supervised machine learning algorithms to predict the expression of 1,364 TUs given TF activities using 441 samples. The algorithms accurately predicted condition-specific expression for 86% (1,174 of 1,364) of the TUs, while 193 TUs (14%) were predicted better than random TRNs. Third, we identified 10 regulatory modules whose definitions were robust against changes to the TRN or expression compendium. Using surrogate variable analysis, we also identified three unmodeled factors that systematically influenced gene expression. Our computational workflow comprehensively characterizes the predictive capabilities and systems-level functions of an organism's TRN from disparate data types.

Interplay between tolerance mechanisms to copper and acid stress in Escherichia coli

Copper (Cu) is a key antibacterial component of the host innate immune system and almost all bacterial species possess systems that defend against the toxic effects of excess Cu. The Cu tolerance system in Gram-negative bacteria is composed minimally of a Cu sensor (CueR) and a Cu export pump (CopA). The cueR and copA genes are encoded on the chromosome typically as a divergent but contiguous operon. In Escherichia coli, cueR and copA are separated by two additional genes, ybaS and ybaT, which confer glutamine (Gln)-dependent acid tolerance and contribute to the glutamate (Glu)-dependent acid resistance system in this organism. Here we show that Cu strongly inhibits growth of a ∆copA mutant strain in acidic cultures. We further demonstrate that Cu stress impairs the pathway for Glu biosynthesis via glutamate synthase, leading to decreased intracellular levels of Glu. Addition of exogenous Glu rescues the ∆copA mutant from Cu stress in acidic conditions. Gln is also protective but this relies on the activities of YbaS and YbaT. Notably, expression of both enzymes is up-regulated during Cu stress. These results demonstrate a link between Cu stress, acid stress, and Glu/Gln metabolism, establish a role for YbaS and YbaT in Cu tolerance, and suggest that subtle changes in core metabolic pathways may contribute to overcoming host-imposed copper toxicity.

Revealing genome-scale transcriptional regulatory landscape of OmpR highlights its expanded regulatory roles under osmotic stress in Escherichia coli K-12 MG1655

A transcription factor (TF), OmpR, plays a critical role in transcriptional regulation of the osmotic stress response in bacteria. Here, we reveal a genome-scale OmpR regulon in Escherichia coli K-12 MG1655. Integrative data analysis reveals that a total of 37 genes in 24 transcription units (TUs) belong to OmpR regulon. Among them, 26 genes show more than two-fold changes in expression level in an OmpR knock-out strain. Specifically, we find that: 1) OmpR regulates mostly membrane-located gene products involved in diverse fundamental biological processes, such as narU (encoding nitrate/nitrite transporter), ompX (encoding outer membrane protein X), and nuoN (encoding NADH:ubiquinone oxidoreductase); 2) by investigating co-regulation of entire sets of genes regulated by other stress-response TFs, stresses are surprisingly independently regulated among each other; and, 3) a detailed investigation of the physiological roles of the newly discovered OmpR regulon genes reveals that activation of narU represents a novel strategy to significantly improve osmotic stress tolerance of E. coli. Thus, the genome-scale approach to elucidating regulons comprehensively identifies regulated genes and leads to fundamental discoveries related to stress responses.

Noisy Response to Antibiotic Stress Predicts Subsequent Single-Cell Survival in an Acidic Environment

Antibiotics elicit drastic changes in microbial gene expression, including the induction of stress response genes. While certain stress responses are known to "cross-protect" bacteria from other stressors, it is unclear whether cellular responses to antibiotics have a similar protective role. By measuring the genome-wide transcriptional response dynamics of Escherichia coli to four antibiotics, we found that trimethoprim induces a rapid acid stress response that protects bacteria from subsequent exposure to acid. Combining microfluidics with time-lapse imaging to monitor survival and acid stress response in single cells revealed that the noisy expression of the acid resistance operon gadBC correlates with single-cell survival. Cells with higher gadBC expression following trimethoprim maintain higher intracellular pH and survive the acid stress longer. The seemingly random single-cell survival under acid stress can therefore be predicted from gadBC expression and rationalized in terms of GadB/C molecular function. Overall, we provide a roadmap for identifying the molecular mechanisms of single-cell cross-protection between antibiotics and other stressors.

Genome-Wide Transcriptional Response to Varying RpoS Levels in Escherichia coli K-12

The alternative sigma factor RpoS is a central regulator of many stress responses in Escherichia coli The level of functional RpoS differs depending on the stress. The effect of these differing concentrations of RpoS on global transcriptional responses remains unclear. We investigated the effect of RpoS concentration on the transcriptome during stationary phase in rich media. We found that 23% of genes in the E. coli genome are regulated by RpoS, and we identified many RpoS-transcribed genes and promoters. We observed three distinct classes of response to RpoS by genes in the regulon: genes whose expression changes linearly with increasing RpoS level, genes whose expression changes dramatically with the production of only a little RpoS ("sensitive" genes), and genes whose expression changes very little with the production of a little RpoS ("insensitive"). We show that sequences outside the core promoter region determine whether an RpoS-regulated gene is sensitive or insensitive. Moreover, we show that sensitive and insensitive genes are enriched for specific functional classes and that the sensitivity of a gene to RpoS corresponds to the timing of induction as cells enter stationary phase. Thus, promoter sensitivity to RpoS is a mechanism to coordinate specific cellular processes with growth phase and may also contribute to the diversity of stress responses directed by RpoS.IMPORTANCE The sigma factor RpoS is a global regulator that controls the response to many stresses in Escherichia coli Different stresses result in different levels of RpoS production, but the consequences of this variation are unknown. We describe how changing the level of RpoS does not influence all RpoS-regulated genes equally. The cause of this variation is likely the action of transcription factors that bind the promoters of the genes. We show that the sensitivity of a gene to RpoS levels explains the timing of expression as cells enter stationary phase and that genes with different RpoS sensitivities are enriched for specific functional groups. Thus, promoter sensitivity to RpoS is a mechanism that coordinates specific cellular processes in response to stresses.

Common Distribution of gad Operon in Lactobacillus brevis and its GadA Contributes to Efficient GABA Synthesis toward Cytosolic Near-Neutral pH

Many strains of lactic acid bacteria (LAB) and bifidobacteria have exhibited strain-specific capacity to produce γ-aminobutyric acid (GABA) via their glutamic acid decarboxylase (GAD) system, which is one of amino acid-dependent acid resistance (AR) systems in bacteria. However, the linkage between bacterial AR and GABA production capacity has not been well established. Meanwhile, limited evidence has been provided to the global diversity of GABA-producing LAB and bifidobacteria, and their mechanisms of efficient GABA synthesis. In this study, genomic survey identified common distribution of gad operon-encoded GAD system in Lactobacillus brevis for its GABA production among varying species of LAB and bifidobacteria. Importantly, among four commonly distributed amino acid-dependent AR systems in Lb. brevis, its GAD system was a major contributor to maintain cytosolic pH homeostasis by consuming protons via GABA synthesis. This highlights that Lb. brevis applies GAD system as the main strategy against extracellular and intracellular acidification demonstrating its high capacity of GABA production. In addition, the abundant GadA retained its activity toward near-neutral pH (pH 5.5–6.5) of cytosolic acidity thus contributing to efficient GABA synthesis in Lb. brevis. This is the first global report illustrating species-specific characteristic and mechanism of efficient GABA synthesis in Lb. brevis.

Coordinated regulation of acid resistance in Escherichia coli

Background

Enteric Escherichia coli survives the highly acidic environment of the stomach through multiple acid resistance (AR) mechanisms. The most effective system, AR2, decarboxylates externally-derived glutamate to remove cytoplasmic protons and excrete GABA. The first described system, AR1, does not require an external amino acid. Its mechanism has not been determined. The regulation of the multiple AR systems and their coordination with broader cellular metabolism has not been fully explored.

Results

We utilized a combination of ChIP-Seq and gene expression analysis to experimentally map the regulatory interactions of four TFs: nac, ntrC, ompR, and csiR. Our data identified all previously in vivo confirmed direct interactions and revealed several others previously inferred from gene expression data. Our data demonstrate that nac and csiR directly modulate AR, and leads to a regulatory network model in which all four TFs participate in coordinating acid resistance, glutamate metabolism, and nitrogen metabolism. This model predicts a novel mechanism for AR1 by which the decarboxylation enzymes of AR2 are used with internally derived glutamate. This hypothesis makes several testable predictions that we confirmed experimentally.

Conclusions

Our data suggest that the regulatory network underlying AR is complex and deeply interconnected with the regulation of GABA and glutamate metabolism, nitrogen metabolism. These connections underlie and experimentally validated model of AR1 in which the decarboxylation enzymes of AR2 are used with internally derived glutamate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0376-y) contains supplementary material, which is available to authorized users.

Synthetic redesign of Escherichia coli for cadaverine production from galactose

BACKGROUND

With increasing concerns over the environment, biological production of cadaverine has been suggested as an alternative route to replace polyamides generated from the petroleum-based process. For an ideal bioprocess, cadaverine should be produced with high yield and productivity from various sugars abundant in biomass. However, most microorganisms are not able to efficiently metabolize other biomass-derived sugars as fast as glucose. This results in reduced growth rate and low carbon flux toward the production of desired bio-chemicals. Thus, redesign of microorganisms is necessary for utilizing those carbon sources with enhanced carbon flux and product formation.

RESULTS

In this study, we engineered Escherichia coli to produce cadaverine with rapid assimilation of galactose, a promising future feedstock. To achieve this, genes related to the metabolic pathway were maximally expressed to amplify the flux toward cadaverine production via synthetic expression cassettes consisting of predictive and quantitative genetic parts (promoters, 5'-untranslated regions, and terminators). Furthermore, the feedback inhibition of metabolic enzymes and degradation/re-uptake pathways was inactivated to robustly produce cadaverine. Finally, the resultant strain, DHK4, produced 8.80 g/L cadaverine with high yield (0.170 g/g) and productivity (0.293 g/L/h) during fed-batch fermentation, which was similar to or better than the previous glucose fermentation.

CONCLUSIONS

Taken together, synthetic redesign of a microorganism with predictive and quantitative genetic parts is a prerequisite for converting sugars from abundant biomass into desired platform chemicals. This is the first report to produce cadaverine from galactose. Moreover, the yield (0.170 g/g) was the highest among engineered E. coli systems.

Principles of proteome allocation are revealed using proteomic data and genome-scale models

Integrating omics data to refine or make context-specific models is an active field of constraint-based modeling. Proteomics now cover over 95% of the Escherichia coli proteome by mass. Genome-scale models of Metabolism and macromolecular Expression (ME) compute proteome allocation linked to metabolism and fitness. Using proteomics data, we formulated allocation constraints for key proteome sectors in the ME model. The resulting calibrated model effectively computed the "generalist" (wild-type) E. coli proteome and phenotype across diverse growth environments. Across 15 growth conditions, prediction errors for growth rate and metabolic fluxes were 69% and 14% lower, respectively. The sector-constrained ME model thus represents a generalist ME model reflecting both growth rate maximization and "hedging" against uncertain environments and stresses, as indicated by significant enrichment of these sectors for the general stress response sigma factor σS. Finally, the sector constraints represent a general formalism for integrating omics data from any experimental condition into constraint-based ME models. The constraints can be fine-grained (individual proteins) or coarse-grained (functionally-related protein groups) as demonstrated here. This flexible formalism provides an accessible approach for narrowing the gap between the complexity captured by omics data and governing principles of proteome allocation described by systems-level models.

Genome-Wide Mapping of Binding Sites Reveals Multiple Biological Functions of the Transcription Factor Cst6p in Saccharomyces cerevisiae

In the model eukaryote Saccharomyces cerevisiae, the transcription factor Cst6p has been reported to play important roles in several biological processes. However, the genome-wide targets of Cst6p and its physiological functions remain unknown. Here, we mapped the genome-wide binding sites of Cst6p at high resolution. Cst6p binds to the promoter regions of 59 genes with various biological functions when cells are grown on ethanol but hardly binds to the promoter at any gene when cells are grown on glucose. The retarded growth of the CST6 deletion mutant on ethanol is attributed to the markedly decreased expression of NCE103, encoding a carbonic anhydrase, which is a direct target of Cst6p. The target genes of Cst6p have a large overlap with those of stress-responsive transcription factors, such as Sko1p and Skn7p. In addition, a CST6 deletion mutant growing on ethanol shows hypersensitivity to oxidative stress and ethanol stress, assigning Cst6p as a new member of the stress-responsive transcriptional regulatory network. These results show that mapping of genome-wide binding sites can provide new insights into the function of transcription factors and highlight the highly connected and condition-dependent nature of the transcriptional regulatory network in S. cerevisiae.

Transcription factors regulate the activity of various biological processes through binding to specific DNA sequences. Therefore, the determination of binding positions is important for the understanding of the regulatory effects of transcription factors. In the model eukaryote Saccharomyces cerevisiae, the transcription factor Cst6p has been reported to regulate several biological processes, while its genome-wide targets remain unknown. Here, we mapped the genome-wide binding sites of Cst6p at high resolution. We show that the binding of Cst6p to its target promoters is condition dependent and explain the mechanism for the retarded growth of the CST6 deletion mutant on ethanol. Furthermore, we demonstrate that Cst6p is a new member of a stress-responsive transcriptional regulatory network. These results provide deeper understanding of the function of the dynamic transcriptional regulatory network in S. cerevisiae.

AutA and AutR, Two Novel Global Transcriptional Regulators, Facilitate Avian Pathogenic Escherichia coli Infection

Bacteria can change its lifestyle during inhabiting in host niches where they survive and replicate by rapidly altering gene expression pattern to accommodate the new environment. In this study, two novel regulators in avian pathogenic Escherichia coli (APEC) were identified and designated as AutA and AutR. RT-PCR and β-galactosidase assay results showed that AutA and AutR co-regulated the expression of adhesin UpaB in APEC strain DE205B. Electrophoretic mobility shift assay showed that AutA and AutR could directly bind the upaB promoter DNA. In vitro transcription assay indicated that AutA could activate the upaB transcription, while AutR inhibited the upaB transcription due to directly suppressing the activating effect of AutA on UpaB expression. Transcriptome analysis showed that AutA and AutR coherently affected the expression of hundreds of genes. Our study confirmed that AutA and AutR co-regulated the expression of DE205B K1 capsule and acid resistance systems in E. coli acid fitness island (AFI). Moreover, phenotypic heterogeneity in expression of K1 capsule and acid resistance systems in AFI during host–pathogen interaction was associated with the regulation of AutA and AutR. Collectively speaking, our studies presented that AutA and AutR are involved in APEC adaptive lifestyle change to facilitate its infection.

Transcription profile of Escherichia coli: genomic SELEX search for regulatory targets of transcription factors

Bacterial genomes are transcribed by DNA-dependent RNA polymerase (RNAP), which achieves gene selectivity through interaction with sigma factors that recognize promoters, and transcription factors (TFs) that control the activity and specificity of RNAP holoenzyme. To understand the molecular mechanisms of transcriptional regulation, the identification of regulatory targets is needed for all these factors. We then performed genomic SELEX screenings of targets under the control of each sigma factor and each TF. Here we describe the assembly of 156 SELEX patterns of a total of 116 TFs performed in the presence and absence of effector ligands. The results reveal several novel concepts: (i) each TF regulates more targets than hitherto recognized; (ii) each promoter is regulated by more TFs than hitherto recognized; and (iii) the binding sites of some TFs are located within operons and even inside open reading frames. The binding sites of a set of global regulators, including cAMP receptor protein, LeuO and Lrp, overlap with those of the silencer H-NS, suggesting that certain global regulators play an anti-silencing role. To facilitate sharing of these accumulated SELEX datasets with the research community, we compiled a database, ‘Transcription Profile of Escherichia coli’ (www.shigen.nig.ac.jp/ecoli/tec/).

RegulonDB version 9.0: high-level integration of gene regulation, coexpression, motif clustering and beyond

RegulonDB (http://regulondb.ccg.unam.mx) is one of the most useful and important resources on bacterial gene regulation,as it integrates the scattered scientific knowledge of the best-characterized organism, Escherichia coli K-12, in a database that organizes large amounts of data. Its electronic format enables researchers to compare their results with the legacy of previous knowledge and supports bioinformatics tools and model building. Here, we summarize our progress with RegulonDB since our last Nucleic Acids Research publication describing RegulonDB, in 2013. In addition to maintaining curation up-to-date, we report a collection of 232 interactions with small RNAs affecting 192 genes, and the complete repertoire of 189 Elementary Genetic Sensory-Response units (GENSOR units), integrating the signal, regulatory interactions, and metabolic pathways they govern. These additions represent major progress to a higher level of understanding of regulated processes. We have updated the computationally predicted transcription factors, which total 304 (184 with experimental evidence and 120 from computational predictions); we updated our position-weight matrices and have included tools for clustering them in evolutionary families. We describe our semiautomatic strategy to accelerate curation, including datasets from high-throughput experiments, a novel coexpression distance to search for ‘neighborhood’ genes to known operons and regulons, and computational developments.

Genome-wide Reconstruction of OxyR and SoxRS Transcriptional Regulatory Networks under Oxidative Stress in Escherichia coli K-12 MG1655

Three transcription factors (TFs), OxyR, SoxR, and SoxS, play a critical role in transcriptional regulation of the defense system for oxidative stress in bacteria. However, their full genome-wide regulatory potential is unknown. Here, we perform a genome-scale reconstruction of the OxyR, SoxR, and SoxS regulons in Escherichia coli K-12 MG1655. Integrative data analysis reveals that a total of 68 genes in 51 transcription units (TUs) belong to these regulons. Among them, 48 genes showed more than 2-fold changes in expression level under single-TF-knockout conditions. This reconstruction expands the genome-wide roles of these factors to include direct activation of genes related to amino acid biosynthesis (methionine and aromatic amino acids), cell wall synthesis (lipid A biosynthesis and peptidoglycan growth), and divalent metal ion transport (Mn(2+), Zn(2+), and Mg(2+)). Investigating the co-regulation of these genes with other stress-response TFs reveals that they are independently regulated by stress-specific TFs.