Genome-scale reconstruction of the Lrp regulatory network in Escherichia coli

Beyond membrane permeability: A role for the small RNA MicF in regulation of chromosome replication and partitioning

Small regulatory RNAs (sRNA) have been shown to play a large role in the management of stress responses in Escherichia coli and other bacteria. sRNAs act post-transcriptionally on target mRNA through an imperfect base pairing mechanism to regulate downstream protein expression. The imperfect base pairing allows a single sRNA to bind and regulate a variety mRNA targets which can form intricate regulatory networks that connect different physiological processes for the cell's response. Upon exposure to antimicrobials and superoxide generating agents, the MicF sRNA in E. coli has been shown to regulate a small set of genes involved in the management of membrane permeability. Currently, it is unknown whether MicF acts on other processes to mediate the response to these agents. Using an sRNA interaction prediction tool, we identified genes in E. coli that are potentially regulated by MicF. Through subsequent analysis using a sfGFP-based reporter-gene fusion, we have validated two novel targets of MicF regulation: SeqA, a negative modulator of DNA replication, and ObgE, a GTPase crucial for chromosome partitioning. Importantly, the interaction between MicF and these target mRNAs is contingent upon the presence of the RNA chaperone protein, Hfq. Furthermore, our findings affirm the role of MicF's conserved 5' seed pairing region in initiating these regulatory interactions. Our study suggests that, beyond its established role in membrane permeability management, MicF exerts control over chromosome dynamics in response to distinct environmental cues, implicating a more multifaceted regulatory function in bacterial stress adaptation.

Inferred regulons are consistent with regulator binding sequences in E. coli

The transcriptional regulatory network (TRN) of E. coli consists of thousands of interactions between regulators and DNA sequences. Regulons are typically determined either from resource-intensive experimental measurement of functional binding sites, or inferred from analysis of high-throughput gene expression datasets. Recently, independent component analysis (ICA) of RNA-seq compendia has shown to be a powerful method for inferring bacterial regulons. However, it remains unclear to what extent regulons predicted by ICA structure have a biochemical basis in promoter sequences. Here, we address this question by developing machine learning models that predict inferred regulon structures in E. coli based on promoter sequence features. Models were constructed successfully (cross-validation AUROC > = 0.8) for 85% (40/47) of ICA-inferred E. coli regulons. We found that: 1) The presence of a high scoring regulator motif in the promoter region was sufficient to specify regulatory activity in 40% (19/47) of the regulons, 2) Additional features, such as DNA shape and extended motifs that can account for regulator multimeric binding, helped to specify regulon structure for the remaining 60% of regulons (28/47); 3) investigating regulons where initial machine learning models failed revealed new regulator-specific sequence features that improved model accuracy. Finally, we found that strong regulatory binding sequences underlie both the genes shared between ICA-inferred and experimental regulons as well as genes in the E. coli core pan-regulon of Fur. This work demonstrates that the structure of ICA-inferred regulons largely can be understood through the strength of regulator binding sites in promoter regions, reinforcing the utility of top-down inference for regulon discovery.

The roles of nucleoid-associated proteins and topoisomerases in chromosome structure, strand segregation, and the generation of phenotypic heterogeneity in bacteria

How to adapt to a changing environment is a fundamental, recurrent problem confronting cells. One solution is for cells to organize their constituents into a limited number of spatially extended, functionally relevant, macromolecular assemblies or hyperstructures, and then to segregate these hyperstructures asymmetrically into daughter cells. This asymmetric segregation becomes a particularly powerful way of generating a coherent phenotypic diversity when the segregation of certain hyperstructures is with only one of the parental DNA strands and when this pattern of segregation continues over successive generations. Candidate hyperstructures for such asymmetric segregation in prokaryotes include those containing the nucleoid-associated proteins (NAPs) and the topoisomerases. Another solution to the problem of creating a coherent phenotypic diversity is by creating a growth-environment-dependent gradient of supercoiling generated along the replication origin-to-terminus axis of the bacterial chromosome. This gradient is modulated by transcription, NAPs, and topoisomerases. Here, we focus primarily on two topoisomerases, TopoIV and DNA gyrase in Escherichia coli, on three of its NAPs (H-NS, HU, and IHF), and on the single-stranded binding protein, SSB. We propose that the combination of supercoiling-gradient-dependent and strand-segregation-dependent topoisomerase activities result in significant differences in the supercoiling of daughter chromosomes, and hence in the phenotypes of daughter cells.

Genomic footprinting uncovers global transcription factor responses to amino acids in Escherichia coli

Our knowledge of transcriptional responses to changes in nutrient availability comes primarily from few well-studied transcription factors (TFs), often lacking an unbiased genome-wide perspective. Leveraging recent advances allowing bacterial genomic footprinting, we comprehensively mapped the genome-wide regulatory responses of Escherichia coli to exogenous leucine, methionine, alanine, and lysine. The global TF Lrp was found to individually sense three amino acids and mount three different target gene responses. Overall, 531 genes had altered RNA polymerase occupancy, and 32 TFs responded directly or indirectly to the presence of amino acids, including regulators of membrane and osmotic pressure homeostasis. About 70% of the detected TF-DNA interactions had not been reported before. We thus identified 682 previously unknown TF-binding locations, for a subset of which the involved TFs were identified by affinity purification. This comprehensive map of amino acid regulation illustrates the incompleteness of the known transcriptional regulation network, even in E. coli.

Escherichia coli Leucine-Responsive Regulatory Protein Bridges DNA In Vivo and Tunably Dissociates in the Presence of Exogenous Leucine

Feast-famine response proteins are a widely conserved class of global regulators in prokaryotes, the most highly studied of which is the Escherichia coli leucine-responsive regulatory protein (Lrp). Lrp senses the environmental nutrition status and subsequently regulates up to one-third of the genes in E. coli, either directly or indirectly. Lrp exists predominantly as octamers and hexadecamers (16mers), where leucine is believed to shift the equilibrium toward the octameric state. In this study, we analyzed the effects of three oligomerization state mutants of Lrp in terms of their ability to bind to DNA and regulate gene expression in response to exogenous leucine. We find that oligomerization beyond dimers is required for Lrp's regulatory activity and that, contrary to previous speculation, exogenous leucine modulates Lrp activity at its target promoters exclusively by inhibiting Lrp binding to DNA. We also show evidence that Lrp binding bridges DNA over length scales of multiple kilobases, revealing a new range of mechanisms for Lrp-mediated transcriptional regulation. IMPORTANCE Leucine-responsive regulatory protein (Lrp) is one of the most impactful regulators in E. coli and other bacteria. Lrp senses nutrient conditions and responds by controlling strategies for virulence, cellular motility, and nutrient acquisition. Despite its importance and being evolutionarily highly conserved across bacteria and archaea, several mysteries remain regarding Lrp, including how it actually responds to leucine to change its regulation of targets. Previous studies have led to the hypothesis that Lrp switches between two states, an octamer (8 Lrp molecules together) and a hexadecamer (16 Lrp molecules together), upon exposure to leucine; these are referred to as different oligomerization states. Here, we show that contrary to previous expectations, it is Lrp's propensity to bind DNA, rather than its oligomerization state, that is directly affected by leucine in the cell's environment. Our new understanding of Lrp activity will aid in identifying and disrupting pathways used by bacteria to cause disease.

Hfq protein and GcvB small RNA tailoring of oppA target mRNA to levels allowing translation activation by MicF small RNA in Escherichia coli

Traffic of molecules across the bacterial membrane mainly relies on porins and transporters, whose expression must adapt to environmental conditions. To ensure bacterial fitness, synthesis and assembly of functional porins and transporters are regulated through a plethora of mechanisms. Among them, small regulatory RNAs (sRNAs) are known to be powerful post-transcriptional regulators. In Escherichia coli, the MicF sRNA is known to regulate only four targets, a very narrow targetome for a sRNA responding to various stresses, such as membrane stress, osmotic shock, or thermal shock. Using an in vivo pull-down assay combined with high-throughput RNA sequencing, we sought to identify new targets of MicF to better understand its role in the maintenance of cellular homoeostasis. Here, we report the first positively regulated target of MicF, the oppA mRNA. The OppA protein is the periplasmic component of the Opp ATP-binding cassette (ABC) oligopeptide transporter and regulates the import of short peptides, some of them bactericides. Mechanistic studies suggest that oppA translation is activated by MicF through a mechanism of action involving facilitated access to a translation-enhancing region in oppA 5'UTR. Intriguingly, MicF activation of oppA translation depends on cross-regulation by negative trans-acting effectors, the GcvB sRNA and the RNA chaperone protein Hfq.

The DNA relaxation-dependent OFF-to-ON biasing of the type 1 fimbrial genetic switch requires the Fis nucleoid-associated protein

The structural genes expressing type 1 fimbriae in Escherichia coli alternate between expressed (phase ON) and non-expressed (phase OFF) states due to inversion of the 314 bp fimS genetic switch. The FimB tyrosine integrase inverts fimS by site-specific recombination, alternately connecting and disconnecting the fim operon, encoding the fimbrial subunit protein and its associated secretion and adhesin factors, to and from its transcriptional promoter within fimS. Site-specific recombination by the FimB recombinase becomes biased towards phase ON as DNA supercoiling is relaxed, a condition that occurs when bacteria approach the stationary phase of the growth cycle. This effect can be mimicked in exponential phase cultures by inhibiting the negative DNA supercoiling activity of DNA gyrase. We report that this bias towards phase ON depends on the presence of the Fis nucleoid-associated protein. We mapped the Fis binding to a site within the invertible fimS switch by DNase I footprinting. Disruption of this binding site by base substitution mutagenesis abolishes both Fis binding and the ability of the mutated switch to sustain its phase ON bias when DNA is relaxed, even in bacteria that produce the Fis protein. In addition, the Fis binding site overlaps one of the sites used by the Lrp protein, a known directionality determinant of fimS inversion that also contributes to phase ON bias. The Fis–Lrp relationship at fimS is reminiscent of that between Fis and Xis when promoting DNA relaxation-dependent excision of bacteriophage λ from the E. coli chromosome. However, unlike the co-binding mechanism used by Fis and Xis at λ attR, the Fis–Lrp relationship at fimS involves competitive binding. We discuss these findings in the context of the link between fimS inversion biasing and the physiological state of the bacterium.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, E_max, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with E_max but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and E_max values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Gene Networks and Pathways Involved in Escherichia coli Response to Multiple Stressors

Stress response helps microorganisms survive extreme environmental conditions and host immunity, making them more virulent or drug resistant. Although both reductionist approaches investigating specific genes and systems approaches analyzing individual stress conditions are being used, less is known about gene networks involved in multiple stress responses. Here, using a systems biology approach, we mined hundreds of transcriptomic data sets for key genes and pathways involved in the tolerance of the model microorganism Escherichia coli to multiple stressors. Specifically, we investigated the E. coli K-12 MG1655 transcriptome under five stresses: heat, cold, oxidative stress, nitrosative stress, and antibiotic treatment. Overlaps of transcriptional changes between studies of each stress factor and between different stressors were determined: energy-requiring metabolic pathways, transport, and motility are typically downregulated to conserve energy, while genes related to survival, bona fide stress response, biofilm formation, and DNA repair are mainly upregulated. The transcription of 15 genes with uncharacterized functions is higher in response to multiple stressors, which suggests they may play pivotal roles in stress response. In conclusion, using rank normalization of transcriptomic data, we identified a set of E. coli stress response genes and pathways, which could be potential targets to overcome antibiotic tolerance or multidrug resistance.

Auxotrophic and prototrophic conditional genetic networks reveal the rewiring of transcription factors in Escherichia coli

Bacterial transcription factors (TFs) are widely studied in Escherichia coli. Yet it remains unclear how individual genes in the underlying pathways of TF machinery operate together during environmental challenge. Here, we address this by applying an unbiased, quantitative synthetic genetic interaction (GI) approach to measure pairwise GIs among all TF genes in E. coli under auxotrophic (rich medium) and prototrophic (minimal medium) static growth conditions. The resulting static and differential GI networks reveal condition-dependent GIs, widespread changes among TF genes in metabolism, and new roles for uncharacterized TFs (yjdC, yneJ, ydiP) as regulators of cell division, putrescine utilization pathway, and cold shock adaptation. Pan-bacterial conservation suggests TF genes with GIs are co-conserved in evolution. Together, our results illuminate the global organization of E. coli TFs, and remodeling of genetic backup systems for TFs under environmental change, which is essential for controlling the bacterial transcriptional regulatory circuits.

The bacterium E. coli has around 300 transcriptional factors, but the functions of many of them, and the interactions between their respective regulatory networks, are unclear. Here, the authors study genetic interactions among all transcription factor genes in E. coli, revealing condition-dependent interactions and roles for uncharacterized transcription factors.

The Escherichia coli Amino Acid Uptake Protein CycA: Regulation of Its Synthesis and Practical Application in l-Isoleucine Production

Amino acid transport systems perform important physiological functions; their role should certainly be considered in microbial production of amino acids. Typically, in the context of metabolic engineering, efforts are focused on the search for and application of specific amino acid efflux pumps. However, in addition, importers can also be used to improve the industrial process as a whole. In this study, the protein CycA, which is known for uptake of nonpolar amino acids, was characterized from the viewpoint of regulating its expression and range of substrates. We prepared a cycA-overexpressing strain and found that it exhibited high sensitivity to branched-chain amino acids and their structural analogues, with relatively increased consumption of these amino acids, suggesting that they are imported by CycA. The expression of cycA was found to be dependent on the extracellular concentrations of substrate amino acids. The role of some transcription factors in cycA expression, including of Lrp and Crp, was studied using a reporter gene construct. Evidence for the direct binding of Crp to the cycA regulatory region was obtained using a gel-retardation assay. The enhanced import of named amino acids due to cycA overexpression in the l-isoleucine-producing strain resulted in a significant reduction in the generation of undesirable impurities. This work demonstrates the importance of uptake systems with respect to their application in metabolic engineering.

What Flips the Switch? Signals and Stress Regulating Extraintestinal Pathogenic Escherichia coli Type 1 Fimbriae (Pili)

Pathogens are exposed to a multitude of harmful conditions imposed by the environment of the host. Bacterial responses against these stresses are pivotal for successful host colonization and pathogenesis. In the case of many E. coli strains, type 1 fimbriae (pili) are an important colonization factor that can contribute to diseases such as urinary tract infections and neonatal meningitis. Production of type 1 fimbriae in E. coli is dependent on an invertible promoter element, fimS, which serves as a phase variation switch determining whether or not a bacterial cell will produce type 1 fimbriae. In this review, we present aspects of signaling and stress involved in mediating regulation of type 1 fimbriae in extraintestinal E. coli; in particular, how certain regulatory mechanisms, some of which are linked to stress response, can influence production of fimbriae and influence bacterial colonization and infection. We suggest that regulation of type 1 fimbriae is potentially linked to environmental stress responses, providing a perspective for how environmental cues in the host and bacterial stress response during infection both play an important role in regulating extraintestinal pathogenic E. coli colonization and virulence.

Synthetic biosensor accelerates evolution by rewiring carbon metabolism toward a specific metabolite

Proper carbon flux distribution between cell growth and production of a target compound is important for biochemical production because improper flux reallocation inhibits cell growth, thus adversely affecting production yield. Here, using a synthetic biosensor to couple production of a specific metabolite with cell growth, we spontaneously evolve cells under the selective condition toward the acquisition of genotypes that optimally reallocate cellular resources. Using 3-hydroxypropionic acid (3-HP) production from glycerol in Escherichia coli as a model system, we determine that mutations in the conserved regions of proteins involved in global transcriptional regulation alter the expression of several genes associated with central carbon metabolism. These changes rewire central carbon flux toward the 3-HP production pathway, increasing 3-HP yield and reducing acetate accumulation by alleviating overflow metabolism. Our study provides a perspective on adaptive laboratory evolution (ALE) using synthetic biosensors, thereby supporting future efforts in metabolic pathway optimization.

Leucine-Responsive Regulatory Protein in Acetic Acid Bacteria Is Stable and Functions at a Wide Range of Intracellular pH Levels

Acetic acid bacteria grow while producing acetic acid, resulting in acidification of the culture. Limited reports elucidate the effect of changes in intracellular pH on transcriptional factors. In the present study, the intracellular pH of Komagataeibacter europaeus was monitored with a pH-sensitive green fluorescent protein, showing that the intracellular pH decreased from 6.3 to 4.7 accompanied by acetic acid production during cell growth. The leucine-responsive regulatory protein of K. europaeus (KeLrp) was used as a model to examine pH-dependent effects, and its properties were compared with those of the Escherichia coli ortholog (EcLrp) at different pH levels. The DNA-binding activities of EcLrp and KeLrp with the target DNA (Ec-ilvI and Ke-ilvI) were examined by gel mobility shift assays under various pH conditions. EcLrp showed the highest affinity with the target at pH 8.0 (K_d [dissociation constant], 0.7 μM), decreasing to a minimum of 3.4 μM at pH 4.0. Conversely, KeLrp did not show significant differences in binding affinity between pH 4 and 7 (K_d, 1.0 to 1.5 μM), and the highest affinity was at pH 5.0 (K_d, 1.0 μM). Circular dichroism spectroscopy revealed that the α-helical content of KeLrp was the highest at pH 5.0 (49%) and was almost unchanged while being maintained at >45% over a range of pH levels examined, while that of EcLrp decreased from its maximum (49% at pH 7.0) to its minimum (36% at pH 4.0). These data indicate that KeLrp is stable and functions over a wide range of intracellular pH levels. IMPORTANCE Lrp is a highly conserved transcriptional regulator found in bacteria and archaea and regulates transcriptions of various genes. The intracellular pH of acetic acid bacteria (AAB) changes accompanied by acetic acid production during cell growth. The Lrp of AAB K. europaeus (KeLrp) was structurally stable over a wide range of pH and maintained DNA-binding activity even at low pH compared with Lrp from E. coli living in a neutral environment. An in vitro experiment showed DNA-binding activity of KeLrp to the target varied with changes in pH. In AAB, change of the intracellular pH during a cell growth would be an important trigger in controlling the activity of Lrp in vivo.

The EcoCyc Database in 2021

The EcoCyc model-organism database collects and summarizes experimental data for Escherichia coli K-12. EcoCyc is regularly updated by the manual curation of individual database entries, such as genes, proteins, and metabolic pathways, and by the programmatic addition of results from select high-throughput analyses. Updates to the Pathway Tools software that supports EcoCyc and to the web interface that enables user access have continuously improved its usability and expanded its functionality. This article highlights recent improvements to the curated data in the areas of metabolism, transport, DNA repair, and regulation of gene expression. New and revised data analysis and visualization tools include an interactive metabolic network explorer, a circular genome viewer, and various improvements to the speed and usability of existing tools.

The leucine-responsive regulatory proteins/feast-famine regulatory proteins: an ancient and complex class of transcriptional regulators in bacteria and archaea

Since the discovery of the Escherichia coli leucine-responsive regulatory protein (Lrp) almost 50 years ago, hundreds of Lrp homologs have been discovered, occurring in 45% of sequenced bacteria and almost all sequenced archaea. Lrp-like proteins are often referred to as the feast/famine regulatory proteins (FFRPs), reflecting their common regulatory roles. Acting as either global or local transcriptional regulators, FFRPs detect the environmental nutritional status by sensing small effector molecules (usually amino acids) and regulate the expression of genes involved in metabolism, virulence, motility, nutrient transport, stress tolerance, and antibiotic resistance to implement appropriate behaviors for the specific ecological niche of each organism. Despite FFRPs' complexity, a significant role in gene regulation, and prevalence throughout prokaryotes, the last comprehensive review on this family of proteins was published about a decade ago. In this review, we integrate recent notable findings regarding E. coli Lrp and other FFRPs across bacteria and archaea with previous observations to synthesize a more complete view on the mechanistic details and biological roles of this ancient class of transcription factors.

A Nitric Oxide-Responsive Transcriptional Regulator NsrR Cooperates With Lrp and CRP to Tightly Control the hmpA Gene in Vibrio vulnificus

Nitric oxide (NO) is an important antimicrobial effector produced by the host innate immune system to counteract invading pathogens. To survive and establish a successful infection, a fulminating human pathogen Vibrio vulnificus expresses the hmpA gene encoding an NO dioxygenase in an NO-responsive manner. In this study, we identified an Rrf2-family transcriptional regulator NsrR that is predicted to contain the Fe-S cluster coordinated by three cysteine residues. Transcriptome analysis showed that NsrR controls the expression of multiple genes potentially involved in nitrosative stress responses. Particularly, NsrR acts as a strong repressor of hmpA transcription and relieves the repression of hmpA upon exposure to NO. Notably, nsrR and hmpA are transcribed divergently, and their promoter regions overlap with each other. Molecular biological analyses revealed that NsrR directly binds to this overlapping promoter region, which is alleviated by loss of the Fe-S cluster, leading to the subsequent derepression of hmpA under nitrosative stress. We further found that a leucine-responsive regulatory protein (Lrp) negatively regulates hmpA in an NsrR-dependent manner by directly binding to the promoter region, presumably resulting in a DNA conformation change to support the repression by NsrR. Meanwhile, a cyclic AMP receptor protein (CRP) positively regulates hmpA probably through repression of nsrR and lrp by directly binding to each promoter region in a sequential cascade. Altogether, this collaborative regulation of NsrR along with Lrp and CRP enables an elaborate control of hmpA transcription, contributing to survival under host-derived nitrosative stress and thereby the pathogenesis of V. vulnificus.

Microbial methionine transporters and biotechnological applications

Methionine (Met) is an essential amino acid with commercial value in animal feed, human nutrition, and as a chemical precursor. Microbial production of Met has seen intensive investigation towards a more sustainable alternative to the chemical synthesis that currently meets the global Met demand. Indeed, efficient Met biosynthesis has been achieved in genetically modified bacteria that harbor engineered enzymes and streamlined metabolic pathways. Very recently, the export of Met as the final step during its fermentative production has been studied and optimized, primarily through identification and expression of microbial Met efflux transporters. In this mini-review, we summarize the current knowledge on four families of Met export and import transporters that have been harnessed for the production of Met and other valuable biomolecules. These families are discussed with respect to their function, gene regulation, and biotechnological applications. We cover methods for identification and characterization of Met transporters as the basis for the further engineering of these proteins and for exploration of other solute carrier families. The available arsenal of Met transporters from different species and protein families provides blueprints not only for fermentative production but also synthetic biology systems, such as molecular sensors and cell-cell communication systems. KEY POINTS: • Sustainable production of methionine (Met) using microbes is actively explored. • Met transporters of four families increase production yield and specificity. • Further applications include other biosynthetic pathways and synthetic biology.

Global Transcriptional Regulators Fine-Tune the Translational and Metabolic Efficiency for Optimal Growth of Escherichia coli

Global transcriptional regulators coordinate complex genetic interactions that bestow better adaptability for an organism against external and internal perturbations. These transcriptional regulators are known to control an enormous array of genes with diverse functionalities. However, regulator-driven molecular mechanisms that underpin precisely tuned translational and metabolic processes conducive for rapid exponential growth remain obscure. Here, we comprehensively reveal the fundamental role of global transcriptional regulators FNR, ArcA, and IHF in sustaining translational and metabolic efficiency under glucose fermentative conditions in Escherichia coli By integrating high-throughput gene expression profiles and absolute intracellular metabolite concentrations, we illustrate that these regulators are crucial in maintaining nitrogen homeostasis, govern expression of otherwise unnecessary or hedging genes, and exert tight control on metabolic bottleneck steps. Furthermore, we characterize changes in expression and activity profiles of other coregulators associated with these dysregulated metabolic pathways, determining the regulatory interactions within the transcriptional regulatory network. Such systematic findings emphasize their importance in optimizing the proteome allocation toward metabolic enzymes as well as ribosomes, facilitating condition-specific phenotypic outcomes. Consequentially, we reveal that disruption of this inherent trade-off between ribosome and metabolic proteome economy due to the loss of regulators resulted in lowered growth rates. Moreover, our findings reinforce that the accumulations of intracellular metabolites in the event of proteome repartitions negatively affects the glucose uptake. Overall, by extending the three-partition proteome allocation theory concordant with multi-omics measurements, we elucidate the physiological consequences of loss of global regulators on central carbon metabolism restraining the organism to attain maximal biomass synthesis.IMPORTANCE Cellular proteome allocation in response to environmental or internal perturbations governs their eventual phenotypic outcome. This entails striking an effective balance between amino acid biosynthesis by metabolic proteins and its consumption by ribosomes. However, the global transcriptional regulator-driven molecular mechanisms that underpin their coordination remains unexplored. Here, we emphasize that global transcriptional regulators, known to control expression of a myriad of genes, are fundamental for precisely tuning the translational and metabolic efficiencies that define the growth optimality. Towards this, we systematically characterized the single deletion effect of FNR, ArcA, and IHF regulators of Escherichia coli on exponential growth under anaerobic glucose fermentative conditions. Their absence disrupts the stringency of proteome allocation, which manifests as impairment in key metabolic processes and an accumulation of intracellular metabolites. Furthermore, by incorporating an extension to the empirical growth laws, we quantitatively demonstrate the general design principles underlying the existence of these regulators in E. coli.

A Thermosensitive, Phase-Variable Epigenetic Switch: pap Revisited

It has been more than a decade since the last comprehensive review of the phase-variable uropathogen-associated pyelonephritis-associated pilus (pap) genetic switch. Since then, important data have come to light, including additional factors that regulate pap expression, better characterization of H-NS regulation, the structure of the Lrp octamer in complex with pap regulatory DNA, the temperature-insensitive phenotype of a mutant lacking the acetyltransferase RimJ, evidence that key components of the regulatory machinery are acetylated, and new insights into the role of DNA binding by key regulators in shaping both the physical structure and regulatory state of the papI and papBA promoters. This review revisits pap, integrating these newer observations with older ones to produce a new model for the concerted behavior of this virulence-regulatory region.

A MARTX Toxin rtxA Gene Is Controlled by Host Environmental Signals through a CRP-Coordinated Regulatory Network in Vibrio vulnificus

A MARTX toxin, RtxA, is an essential virulence factor of many pathogens, including Vibrio species. H-NS and HlyU repress and derepress, respectively, rtxA expression of a life-threatening pathogen, Vibrio vulnificus. We found that Lrp directly activates rtxA independently of H-NS and HlyU, and leucine inhibits the Lrp-mediated activation of rtxA. Furthermore, we demonstrated that CRP represses rtxA but derepresses in the presence of exogenous glucose. CRP represses rtxA not only directly by binding to upstream of rtxA but also indirectly by repressing lrp and hlyU. This is the first report of a regulatory network comprising CRP, Lrp, H-NS, and HlyU, which coordinates the rtxA expression in response to environmental signals such as leucine and glucose during infection. This elaborate regulatory network will enhance the fitness of V. vulnificus and contribute to its successful infection within the host.

A multifunctional autoprocessing repeats-in-toxin (MARTX) toxin plays an essential role in the virulence of many pathogens, including a fulminating human pathogen Vibrio vulnificus. H-NS and HlyU repress and derepress expression of the MARTX toxin gene rtxA in V. vulnificus, respectively. However, little is known about other regulatory proteins and environmental signals involved in rtxA regulation. In this study, we found that a leucine-responsive regulatory protein (Lrp) activates rtxA by binding directly and specifically to the rtxA promoter, P_rtxA. Phased hypersensitivity resulting from DNase I cleavage of the P_rtxA regulatory region suggests that Lrp probably induces DNA bending in P_rtxA. Lrp activates P_rtxA independently of H-NS and HlyU, and leucine inhibits Lrp binding to P_rtxA and reduces the Lrp-mediated activation. Furthermore, a cyclic AMP receptor protein (CRP) represses P_rtxA, and exogenous glucose relieves the CRP-mediated repression. Biochemical and mutational analyses demonstrated that CRP binds directly and specifically to the upstream region of P_rtxA, which presumably alters the DNA conformation in P_rtxA and thus represses rtxA. Moreover, CRP represses expression of lrp and hlyU by binding directly to their upstream regions, forming coherent feed-forward loops with Lrp and HlyU. In conclusion, expression of rtxA is controlled by a regulatory network comprising CRP, Lrp, H-NS, and HlyU in response to changes in host environmental signals such as leucine and glucose. This collaborative regulation enables the elaborate expression of rtxA, thereby enhancing the fitness and pathogenesis of V. vulnificus during the course of infection.

Covert Cross-Feeding Revealed by Genome-Wide Analysis of Fitness Determinants in a Synthetic Bacterial Mutualism

Microbial interactions abound in natural ecosystems and shape community structure and function. Substantial attention has been given to cataloging mechanisms by which microbes interact, but there is a limited understanding of the genetic landscapes that promote or hinder microbial interactions. We previously developed a mutualistic coculture pairing Escherichia coli and Rhodopseudomonas palustris, wherein E. coli provides carbon to R. palustris in the form of glucose fermentation products and R. palustris fixes N₂ gas and provides nitrogen to E. coli in the form of NH₄ ⁺ The stable coexistence and reproducible trends exhibited by this coculture make it ideal for interrogating the genetic underpinnings of a cross-feeding mutualism. Here, we used random barcode transposon sequencing (RB-TnSeq) to conduct a genome-wide search for E. coli genes that influence fitness during cooperative growth with R. palustris RB-TnSeq revealed hundreds of genes that increased or decreased E. coli fitness in a mutualism-dependent manner. Some identified genes were involved in nitrogen sensing and assimilation, as expected given the coculture design. The other identified genes were involved in diverse cellular processes, including energy production and cell wall and membrane biogenesis. In addition, we discovered unexpected purine cross-feeding from R. palustris to E. coli, with coculture rescuing growth of an E. coli purine auxotroph. Our data provide insight into the genes and gene networks that can influence a cross-feeding mutualism and underscore that microbial interactions are not necessarily predictable a priori IMPORTANCE Microbial communities impact life on Earth in profound ways, including driving global nutrient cycles and influencing human health and disease. These community functions depend on the interactions that resident microbes have with the environment and each other. Thus, identifying genes that influence these interactions will aid the management of natural communities and the use of microbial consortia as biotechnology. Here, we identified genes that influenced Escherichia coli fitness during cooperative growth with a mutualistic partner, Rhodopseudomonas palustris Although this mutualism centers on the bidirectional exchange of essential carbon and nitrogen, E. coli fitness was positively and negatively affected by genes involved in diverse cellular processes. Furthermore, we discovered an unexpected purine cross-feeding interaction. These results contribute knowledge on the genetic foundation of a microbial cross-feeding interaction and highlight that unanticipated interactions can occur even within engineered microbial communities.

Diversity, versatility and complexity of bacterial gene regulation mechanisms: opportunities and drawbacks for applications in synthetic biology

Gene expression occurs in two essential steps: transcription and translation. In bacteria, the two processes are tightly coupled in time and space, and highly regulated. Tight regulation of gene expression is crucial. It limits wasteful consumption of resources and energy, prevents accumulation of potentially growth inhibiting reaction intermediates, and sustains the fitness and potential virulence of the organism in a fluctuating, competitive and frequently stressful environment. Since the onset of studies on regulation of enzyme synthesis, numerous distinct regulatory mechanisms modulating transcription and/or translation have been discovered. Mostly, various regulatory mechanisms operating at different levels in the flow of genetic information are used in combination to control and modulate the expression of a single gene or operon. Here, we provide an extensive overview of the very diverse and versatile bacterial gene regulatory mechanisms with major emphasis on their combined occurrence, intricate intertwinement and versatility. Furthermore, we discuss the potential of well-characterized basal expression and regulatory elements in synthetic biology applications, where they may ensure orthogonal, predictable and tunable expression of (heterologous) target genes and pathways, aiming at a minimal burden for the host.

Competitive Repression of the artPIQM Operon for Arginine and Ornithine Transport by Arginine Repressor and Leucine-Responsive Regulatory Protein in Escherichia coli

Two out of the three major uptake systems for arginine in Escherichia coli are encoded by the artJ-artPIQM gene cluster. ArtJ is the high-affinity periplasmic arginine-specific binding protein (ArgBP-I), whereas artI encodes the arginine and ornithine periplasmic binding protein (AO). Both ArtJ and ArtI are supposed to combine with the inner membrane-associated ArtQMP₂ transport complex of the ATP-binding cassette-type (ABC). Transcription of artJ is repressed by arginine repressor (ArgR) and the artPIQM operon is regulated by the transcriptional regulators ArgR and Leucine-responsive regulatory protein (Lrp). Whereas repression by ArgR requires arginine as corepressor, repression of P_artP by Lrp is partially counteracted by leucine, its major effector molecule. We demonstrate that binding of dimeric Lrp to the artP control region generates four complexes with a distinct migration velocity, and that leucine has an effect on both global binding affinity and cooperativity in the binding. We identify the binding sites for Lrp in the artP control region, reveal interferences in the binding of ArgR and Lrp in vitro and demonstrate that the two transcription factors act as competitive repressors in vivo, each one being a more potent regulator in the absence of the other. This competitive behavior may be explained by the partial steric overlap of their respective binding sites. Furthermore, we demonstrate ArgR binding to an unusual position in the control region of the lrp gene, downstream of the transcription initiation site. From this unusual position for an ArgR-specific operator, ArgR has little direct effect on lrp expression, but interferes with the negative leucine-sensitive autoregulation exerted by Lrp. Direct arginine and ArgR-dependent repression of lrp could be observed with a 25-bp deletion mutant, in which the ArgR binding site was artificially moved to a position immediately downstream of the lrp transcription initiation site. This finding is reminiscent of a previous observation made for the carAB operon encoding carbamoylphosphate synthase, where ArgR bound in overlap with the downstream promoter P2 does not block transcription initiated 67 bp upstream at the P1 promoter, and further supports the hypothesis that ArgR does not act as an efficient roadblock.

Escherichia coli Lrp Regulates One-Third of the Genome via Direct, Cooperative, and Indirect Routes

The global regulator Lrp plays a crucial role in regulating metabolism, virulence, and motility in response to environmental conditions. Lrp has previously been shown to activate or repress approximately 10% of the genes in Escherichia coli However, the full spectrum of targets, and how Lrp acts to regulate them, have stymied earlier study. We have combined matched chromatin-immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) under nine physiological conditions to comprehensively map the binding and regulatory activity of Lrp as it directs responses to nutrient abundance. In addition to identifying hundreds of novel Lrp targets, we observe two new global trends, as follows: first, that Lrp will often bind to promoters in a poised position under conditions when it has no regulatory activity to enable combinatorial interactions with other regulators, and second, that nutrient levels induce a global shift in the equilibrium between less-sequence-specific and more-sequence-specific DNA binding. The overall regulatory behavior of Lrp, which as we now show extends to 38% of E. coli genes directly or indirectly under at least one condition, thus arises from the interaction between changes in Lrp binding specificity and cooperative action with other regulators.IMPORTANCE To survive, bacteria such as E. coli must rapidly respond to changing environmental conditions, including nutrient levels. A decrease in nutrient availability causes bacteria to stop rapid replication and enter stationary phase, where they perform limited to no cell division. The E. coli global regulatory protein Lrp has been previously implicated in modulating the expression of genes particularly important at this transition from rapid to slowed growth. Here, we monitor Lrp's DNA binding locations and effect on gene expression under three different nutrient conditions across three growth stages. We find that Lrp's role is even broader than previously suspected and that it appears to interact with many other bacterial regulators to perform its function in a condition-specific manner.

Allosteric Feedback Inhibition Enables Robust Amino Acid Biosynthesis in E. coli by Enforcing Enzyme Overabundance

Microbes must ensure robust amino acid metabolism in the face of external and internal perturbations. This robustness is thought to emerge from regulatory interactions in metabolic and genetic networks. Here, we explored the consequences of removing allosteric feedback inhibition in seven amino acid biosynthesis pathways in Escherichia coli (arginine, histidine, tryptophan, leucine, isoleucine, threonine, and proline). Proteome data revealed that enzyme levels decreased in five of the seven dysregulated pathways. Despite that, flux through the dysregulated pathways was not limited, indicating that enzyme levels are higher than absolutely needed in wild-type cells. We showed that such enzyme overabundance renders the arginine, histidine, and tryptophan pathways robust against perturbations of gene expression, using a metabolic model and CRISPR interference experiments. The results suggested a sensitive interaction between allosteric feedback inhibition and enzyme-level regulation that ensures robust yet efficient biosynthesis of histidine, arginine, and tryptophan in E. coli.

•

Amino acid biosynthesis enzymes do not normally operate at maximum capacity

•

Allosteric feedback inhibition ensures that enzymes are overabundant

•

Enzyme overabundance provides robustness against decreases in gene expression

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.

Newly Identified Nucleoid-Associated-Like Protein YlxR Regulates Metabolic Gene Expression in Bacillus subtilis

Expression of genes encoding NAPs is often temporally regulated. According to results from single-cell analysis, the ylxR gene is induced by glucose and expressed in a bistable mode. These characteristics have not previously been reported for NAP gene expression. Transcriptional profiling of the ylxR disruptant revealed a change in the expression levels of approximately 400 genes, including genes for synthesis of 12 amino acids and 4 nucleotides, in addition to the SigX/M regulons. Thus, YlxR is a critical regulator of glucose response in B. subtilis.

Glucose is the most favorable carbon source for the majority of bacteria, which have several glucose-responsive gene networks. Recently, we found that in Bacillus subtilis, glucose induces expression of the extracellular sigma factor genes sigX/M. To explore the factors affecting this phenomenon, we performed a transposon mutagenesis screen for mutants with no glucose induction (GI) of sigX-lacZ and identified ylxR. YlxR is widely conserved in eubacteria. Further analysis revealed that ylxR is induced by glucose addition. In vitro DNA-binding and cytological studies suggested that YlxR is a nucleoid-associated protein (NAP) in B. subtilis. In many cases, NAPs influence transcription, recombination, and genome stability. Thus, we performed transcriptome sequencing (RNA-Seq) analysis to evaluate the impact of ylxR disruption on the transcriptome in the presence of glucose and observed that YlxR has a profound impact on metabolic gene expression in addition to that of four sigma factor genes. The wide fluctuations of gene expression may result in abolition of GI of sigX/M in the ylxR disruptant.

IMPORTANCE Expression of genes encoding NAPs is often temporally regulated. According to results from single-cell analysis, the ylxR gene is induced by glucose and expressed in a bistable mode. These characteristics have not previously been reported for NAP gene expression. Transcriptional profiling of the ylxR disruptant revealed a change in the expression levels of approximately 400 genes, including genes for synthesis of 12 amino acids and 4 nucleotides, in addition to the SigX/M regulons. Thus, YlxR is a critical regulator of glucose response in B. subtilis.

Control of proline utilization by the Lrp-like regulator PutR in Caulobacter crescentus

Cellular metabolism recently emerged as a central player modulating the bacterial cell cycle. The Alphaproteobacterium Caulobacter crescentus appears as one of the best models to study these connections, but its metabolism is still poorly characterized. Considering that it lives in oligotrophic environments, its capacity to use amino-acids is often critical for its growth. Here, we characterized the C. crescentus PutA bi-functional enzyme and showed that it is required for the utilization of proline as a carbon source. We also found that putA transcription and proline utilization by PutA are strictly dependent on the Lrp-like PutR activator. The activation of putA by PutR needs proline, which most likely acts as an effector molecule for PutR. Surprisingly, we also observed that an over-production of PutR leads to cell elongation in liquid medium containing proline, while it inhibits colony formation even in the absence of proline on solid medium. These cell division and growth defects were equally pronounced in a ΔputA mutant background, indicating that PutR can play other roles beyond the control of proline catabolism. Altogether, these findings suggest that PutR might connect central metabolism with cell cycle processes.

A Requirement for Global Transcription Factor Lrp in Licensing Replication of Vibrio cholerae Chromosome 2

The human pathogen, Vibrio cholerae, belongs to the 10% of bacteria in which the genome is divided. Each of its two chromosomes, like bacterial chromosomes in general, replicates from a unique origin at fixed times in the cell cycle. Chr1 initiates first, and upon duplication of a site in Chr1, crtS, Chr2 replication initiates. Recent in vivo experiments demonstrate that crtS binds the Chr2-specific initiator RctB and promotes its initiator activity by remodeling it. Compared to the well-defined RctB binding sites in the Chr2 origin, crtS is an order of magnitude longer, suggesting that other factors can bind to it. We developed an in vivo screen to identify additional crtS-binding proteins and identified the global transcription factor, Lrp, as one such protein. Studies in vivo and in vitro indicate that Lrp binds to crtS and facilitates RctB binding to crtS. Chr2 replication is severely defective in the absence of Lrp, indicative of a critical role of the transcription factor in licensing Chr2 replication. Since Lrp responds to stresses such as nutrient limitation, its interaction with RctB presumably sensitizes Chr2 replication to the physiological state of the cell.

A unified resource for transcriptional regulation in Escherichia coli K-12 incorporating high-throughput-generated binding data into RegulonDB version 10.0

BACKGROUND

Our understanding of the regulation of gene expression has benefited from the availability of high-throughput technologies that interrogate the whole genome for the binding of specific transcription factors and gene expression profiles. In the case of widely used model organisms, such as Escherichia coli K-12, the new knowledge gained from these approaches needs to be integrated with the legacy of accumulated knowledge from genetic and molecular biology experiments conducted in the pre-genomic era in order to attain the deepest level of understanding possible based on the available data.

RESULTS

In this paper, we describe an expansion of RegulonDB, the database containing the rich legacy of decades of classic molecular biology experiments supporting what we know about gene regulation and operon organization in E. coli K-12, to include the genome-wide dataset collections from 32 ChIP and 19 gSELEX publications, in addition to around 60 genome-wide expression profiles relevant to the functional significance of these datasets and used in their curation. Three essential features for the integration of this information coming from different methodological approaches are: first, a controlled vocabulary within an ontology for precisely defining growth conditions; second, the criteria to separate elements with enough evidence to consider them involved in gene regulation from isolated transcription factor binding sites without such support; and third, an expanded computational model supporting this knowledge. Altogether, this constitutes the basis for adequately gathering and enabling the comparisons and integration needed to manage and access such wealth of knowledge.

CONCLUSIONS

This version 10.0 of RegulonDB is a first step toward what should become the unifying access point for current and future knowledge on gene regulation in E. coli K-12. Furthermore, this model platform and associated methodologies and criteria can be emulated for gathering knowledge on other microbial organisms.

Regulation of stringent factor by branched-chain amino acids

When faced with amino acid starvation, prokaryotic cells induce a stringent response that modulates their physiology. The stringent response is manifested by production of signaling molecules guanosine 5'-diphosphate,3'-diphosphate (ppGpp) and guanosine 5'-triphosphate,3'-diphosphate (pppGpp) that are also called alarmones. In many species, alarmone levels are regulated by a multidomain bifunctional alarmone synthetase/hydrolase called Rel. In this enzyme, there is an ACT domain at the carboxyl region that has an unknown function; however, similar ACT domains are present in other enzymes that have roles in controlling amino acid metabolism. In many cases, these other ACT domains have been shown to allosterically regulate enzyme activity through the binding of amino acids. Here, we show that the ACT domain present in the Rhodobacter capsulatus Rel alarmone synthetase/hydrolase binds branched-chain amino acids valine and isoleucine. We further show that the binding of these amino acids stimulates alarmone hydrolase activity both in vitro and in vivo. Furthermore, we found that the ACT domain present in Rel proteins from many diverse species also binds branched-chain amino acids. These results indicate that the cellular concentration of amino acids can directly affect Rel alarmone synthetase/hydrolase activity, thus adding another layer of control to current models of cellular control of the stringent response.

ChIP-exo interrogation of Crp, DNA, and RNAP holoenzyme interactions

Numerous in vitro studies have yielded a refined picture of the structural and molecular associations between Cyclic-AMP receptor protein (Crp), the DNA motif, and RNA polymerase (RNAP) holoenzyme. In this study, high-resolution ChIP-exonuclease (ChIP-exo) was applied to study Crp binding in vivo and at genome-scale. Surprisingly, Crp was found to provide little to no protection of the DNA motif under activating conditions. Instead, Crp demonstrated binding patterns that closely resembled those generated by σ⁷⁰. The binding patterns of both Crp and σ⁷⁰ are indicative of RNAP holoenzyme DNA footprinting profiles associated with stages during transcription initiation that occur post-recruitment. This is marked by a pronounced advancement of the template strand footprint profile to the +20 position relative to the transcription start site and a multimodal distribution on the nontemplate strand. This trend was also observed in the familial transcription factor, Fnr, but full protection of the motif was seen in the repressor ArcA. Given the time-scale of ChIP studies and that the rate-limiting step in transcription initiation is typically post recruitment, we propose a hypothesis where Crp is absent from the DNA motif but remains associated with RNAP holoenzyme post-recruitment during transcription initiation. The release of Crp from the DNA motif may be a result of energetic changes that occur as RNAP holoenzyme traverses the various stable intermediates towards elongation complex formation.

The Global Acetylome of the Human Pathogen Vibrio cholerae V52 Reveals Lysine Acetylation of Major Transcriptional Regulators

Protein lysine acetylation is recognized as an important reversible post translational modification in all domains of life. While its primary roles appear to reside in metabolic processes, lysine acetylation has also been implicated in regulating pathogenesis in bacteria. Several global lysine acetylome analyses have been carried out in various bacteria, but thus far there have been no reports of lysine acetylation taking place in the important human pathogen Vibrio cholerae. In this study, we analyzed the lysine acetylproteome of the human pathogen V. cholerae V52. By applying a combination of immuno-enrichment of acetylated peptides and high resolution mass spectrometry, we identified 3,402 acetylation sites on 1,240 proteins. Of the acetylated proteins, more than half were acetylated on two or more sites. As reported for other bacteria, we observed that many of the acetylated proteins were involved in metabolic and cellular processes and there was an over-representation of acetylated proteins involved in protein synthesis. Of interest, we demonstrated that many global transcription factors such as CRP, H-NS, IHF, Lrp and RpoN as well as transcription factors AphB, TcpP, and PhoB involved in direct regulation of virulence in V. cholerae were acetylated. In conclusion, this is the first global protein lysine acetylome analysis of V. cholerae and should constitute a valuable resource for in-depth studies of the impact of lysine acetylation in pathogenesis and other cellular processes.

Dissecting the genetic and metabolic mechanisms of adaptation to the knockout of a major metabolic enzyme in Escherichia coli

Unraveling the mechanisms of microbial adaptive evolution following genetic or environmental challenges is of fundamental interest in biological science and engineering. When the challenge is the loss of a metabolic enzyme, adaptive responses can also shed significant insight into metabolic robustness, regulation, and areas of kinetic limitation. In this study, whole-genome sequencing and high-resolution ¹³C-metabolic flux analysis were performed on 10 adaptively evolved pgi knockouts of Escherichia coliPgi catalyzes the first reaction in glycolysis, and its loss results in major physiological and carbon catabolism pathway changes, including an 80% reduction in growth rate. Following adaptive laboratory evolution (ALE), the knockouts increase their growth rate by up to 3.6-fold. Through combined genomic-fluxomic analysis, we characterized the mutations and resulting metabolic fluxes that enabled this fitness recovery. Large increases in pyridine cofactor transhydrogenase flux, correcting imbalanced production of NADPH and NADH, were enabled by direct mutations to the transhydrogenase genes sthA and pntAB The phosphotransferase system component crr was also found to be frequently mutated, which corresponded to elevated flux from pyruvate to phosphoenolpyruvate. The overall energy metabolism was found to be strikingly robust, and what have been previously described as latently activated Entner-Doudoroff and glyoxylate shunt pathways are shown here to represent no real increases in absolute flux relative to the wild type. These results indicate that the dominant mechanism of adaptation was to relieve the rate-limiting steps in cofactor metabolism and substrate uptake and to modulate global transcriptional regulation from stress response to catabolism.

Methanol assimilation in Escherichia coli is improved by co-utilization of threonine and deletion of leucine-responsive regulatory protein

Methane, the main component of natural gas, can be used to produce methanol which can be further converted to other valuable products. There is increasing interest in using biological systems for the production of fuels and chemicals from methanol, termed methylotrophy. In this work, we have examined methanol assimilation metabolism in a synthetic methylotrophic E. coli strain. Specifically, we applied ¹³C-tracers and evaluated 25 different co-substrates for methanol assimilation, including amino acids, sugars and organic acids. In particular, co-utilization of threonine significantly enhanced methylotrophy. Through our investigations, we proposed specific metabolic pathways that, when activated, correlated with increased methanol assimilation. These pathways are normally repressed by the leucine-responsive regulatory protein (lrp), a global regulator of metabolism associated with the feast-or-famine response in E. coli. By deleting lrp, we were able to further enhance the methylotrophic ability of our synthetic strain, as demonstrated through increased incorporation of ¹³C carbon from ¹³C-methanol into biomass.

Stringent Response Regulators Contribute to Recovery from Glucose Phosphate Stress in Escherichia coli

In enteric bacteria such as Escherichia coli, the transcription factor SgrR and the small RNA SgrS regulate the response to glucose phosphate stress, a metabolic dysfunction that results in growth inhibition and stems from the intracellular accumulation of sugar phosphates. SgrR activates the transcription of sgrS, and SgrS helps to rescue cells from stress in part by inhibiting the uptake of stressor sugar phosphates. While the regulatory targets of this stress response are well described, less is known about how the SgrR-SgrS response itself is regulated. To further characterize the regulation of the glucose phosphate stress response, we screened global regulator gene mutants for growth changes during glucose phosphate stress. We found that deleting dksA, which encodes a regulator of the stringent response to nutrient starvation, decreases growth under glucose phosphate stress conditions. The stringent response alarmone regulator ppGpp (synthesized by RelA and SpoT) also contributes to recovery from glucose phosphate stress: as with dksA, mutating relA and spoT worsens the growth defect of an sgrS mutant during stress, although the sgrS relA spoT mutant defect was only detectable under lower stress levels. In addition, mutating dksA or relA and spoT lowers sgrS expression (as measured with a P _sgrS -lacZ fusion), suggesting that the observed growth defects may be due to decreased induction of the glucose phosphate stress response or related targets. This regulatory effect could occur through altered sgrR transcription, as dksA and relA spoT mutants also exhibit decreased expression of a P _sgrR -lacZ fusion. Taken together, this work supports a role for stringent response regulators in aiding the recovery from glucose phosphate stress.IMPORTANCE Glucose phosphate stress leads to growth inhibition in bacteria such as Escherichia coli when certain sugar phosphates accumulate in the cell. The transcription factor SgrR and the small RNA SgrS alleviate this stress in part by preventing further sugar phosphate transport. While the regulatory mechanisms of this response have been characterized, the regulation of the SgrR-SgrS response itself is not as well understood. Here, we describe a role for stringent response regulators DksA and ppGpp in the response to glucose phosphate stress. sgrS dksA and sgrS relA spoT mutants exhibit growth defects under glucose phosphate stress conditions. These defects may be due to a decrease in stress response induction, as deleting dksA or relA and spoT also results in decreased expression of sgrS and sgrR This research presents one of the first regulatory effects on the glucose phosphate stress response outside SgrR and SgrS and depicts a novel connection between these two metabolic stress responses.

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.

Lrp, a global regulator, regulates the virulence of Vibrio vulnificus

Background

An attenuated mutant (designated NY303) of Vibrio vulnificus, which causes serious wound infection and septicemia in humans, was isolated fortuitously from a clinical strain YJ016. This mutant was defective in cytotoxicity, migration on soft agar and virulence in the mouse. The purpose of this study was to map the mutation in this attenuated mutant and further explore how the gene thus identified is involved in virulence.

Methods

The whole genome sequence of mutant NY303 determined by next-generation sequencing was compared with that of strain YJ016 to map the mutations. By isolating and characterizing the specific gene-knockout mutants, the gene associated with the phenotype of mutant NY303 was identified. This gene encodes a global regulator, Lrp. A mutant, YH01, deficient in Lrp was isolated and examined in vitro, in vivo and ex vivo to find the affected virulence mechanisms. The target genes of Lrp were further identified by comparing the transcriptomes, which were determined by RNA-seq, of strain YJ016 and mutant YH01. The promoters bound by Lrp were identified by genome footprinting-sequencing, and those related with virulence were further examined by electrophoretic mobility shift assay.

Results

A mutation in lrp was shown to be associated with the reduced cytotoxicity, chemotaxis and virulence of mutant NY303. Mutant YH01 exhibited a phenotype resembling that of mutant NY303, and was defective in colonization in the mouse and growth in mouse serum, but not the antiphagocytosis ability. 596 and 95 genes were down- and up-regulated, respectively, in mutant YH01. Many of the genes involved in secretion of the MARTX cytotoxin, chemotaxis and iron-acquisition were down-regulated in mutant YH01. The lrp gene, which was shown to be negatively autoregulated, and 7 down-regulated virulence-associated genes were bound by Lrp in their promoters. A 14-bp consensus sequence, mkCrTTkwAyTsTG, putatively recognized by Lrp was identified in the promoters of these genes.

Conclusions

Lrp is a global regulator involved in regulation of cytotoxicity, chemotaxis and iron-acquisition in V. vulnificus. Down-regulation of many of the genes associated with these properties may be responsible, at least partly, for loss of virulence in mutant NY303.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-017-0361-9) contains supplementary material, which is available to authorized users.

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.

Nucleoid-Associated Proteins: Genome Level Occupancy and Expression Analysis

The advent of Chromatin Immunoprecipitation sequencing (ChIP-Seq) has allowed the identification of genomic regions bound by a DNA binding protein in-vivo on a genome-wide scale. The impact of the DNA binding protein on gene expression can be addressed using transcriptome experiments in appropriate genetic settings. Overlaying the above two sources of data enables us to dissect the direct and indirect effects of a DNA binding protein on gene expression. Application of these techniques to Nucleoid Associated Proteins (NAPs) and Global Transcription Factors (GTFs) has underscored the complex relationship between DNA-protein interactions and gene expression change, highlighting the role of combinatorial control. Here, we demonstrate the usage of ChIP-Seq to infer binding properties and transcriptional effects of NAPs such as Fis and HNS, and the GTF CRP in the model organism Escherichia coli K12 MG1655 (E. coli).

The Bacterial Two-Hybrid System Uncovers the Involvement of Acetylation in Regulating of Lrp Activity in Salmonella Typhimurium

N^𝜀-lysine acetylation is an abundant and important Post-translational modification in bacteria. We used the bacterial two-hybrid system to screen the genome library of the Salmonella Typhimurium to identify potential proteins involved in acetyltransferase Pat - or deacetylase CobB-mediated acetylation. Then, the in vitro (de)acetylation assays were used to validate the potential targets, such as STM14_1074, NrdF, RhaR. Lrp, a leucine-responsive regulatory protein and global regulator, was shown to interact with Pat. We further demonstrate that Lrp could be acetylated by Pat and deacetylated by NAD⁺-dependent CobB in vitro. Specifically, the conserved lysine residue 36 (K36) in helix-turn-helix (HTH) DNA-binding domain of Lrp was acetylated. Acetylation of K36 impaired the function of Lrp through altering the affinity with the target promoter. The mutation of K36 in chromosome mimicking acetylation enhanced the transcriptional level of itself and attenuated the mRNA levels of Lrp-regulated genes including fimA, which was confirmed by yeast agglutination assay. These findings demonstrate that the acetylation regulates the DNA-binding activity of Lrp, suggesting that acetylation modification of transcription factors is a conserved regulatory manner to modulate gene expression in bacteria and eukaryotes.

Distinct transcriptional regulation of the two Escherichia coli transhydrogenases PntAB and UdhA

Transhydrogenases catalyse interconversion of the redox cofactors NADH and NADPH, thereby conveying metabolic flexibility to balance catabolic NADPH formation with anabolic or stress-based consumption of NADPH. Escherichia coli is one of the very few microbes that possesses two isoforms: the membrane-bound, proton-translocating transhydrogenase PntAB and the cytosolic, energy-independent transhydrogenase UdhA. Despite their physiological relevance, we have only fragmented information on their regulation and the signals coordinating their counteracting activities. Here we investigated PntAB and UdhA regulation by studying transcriptional responses to environmental and genetic perturbations. By testing pntAB and udhA GFP reporter constructs in the background of WT E. coli and 62 transcription factor mutants during growth on different carbon sources, we show distinct transcriptional regulation of the two transhydrogenase promoters. Surprisingly, transhydrogenase regulation was independent of the actual catabolic overproduction or underproduction of NADPH but responded to nutrient levels and growth rate in a fashion that matches the cellular need for the redox cofactors NADPH and/or NADH. Specifically, the identified transcription factors Lrp, ArgP and Crp link transhydrogenase expression to particular amino acids and intracellular concentrations of cAMP. The overall identified set of regulators establishes a primarily biosynthetic role for PntAB and link UdhA to respiration.

The Landscape of Evolution: Reconciling Structural and Dynamic Properties of Metabolic Networks in Adaptive Diversifications

The network of the interactions among genes, proteins, and metabolites delineates a range of potential phenotypic diversifications in a lineage, and realized phenotypic changes are the result of differences in the dynamics of the expression of the elements and interactions in this deterministic network. Regulatory mechanisms, such as hormones, mediate the relationship between the structural and dynamic properties of networks by determining how and when the elements are expressed and form a functional unit or state. Changes in regulatory mechanisms lead to variable expression of functional states of a network within and among generations. Functional properties of network elements, and the magnitude and direction of evolutionary change they determine, depend on their location within a network. Here, we examine the relationship between network structure and the dynamic mechanisms that regulate flux through a metabolic network. We review the mechanisms that control metabolic flux in enzymatic reactions and examine structural properties of the network locations that are targets of flux control. We aim to establish a predictive framework to test the contributions of structural and dynamic properties of deterministic networks to evolutionary diversifications.

Transcriptome Analysis of Escherichia coli during dGTP Starvation

Our laboratory recently discovered that Escherichia coli cells starved for the DNA precursor dGTP are killed efficiently (dGTP starvation) in a manner similar to that described for thymineless death (TLD). Conditions for specific dGTP starvation can be achieved by depriving an E. coli optA1 gpt strain of the purine nucleotide precursor hypoxanthine (Hx). To gain insight into the mechanisms underlying dGTP starvation, we conducted genome-wide gene expression analyses of actively growing optA1 gpt cells subjected to hypoxanthine deprivation for increasing periods. The data show that upon Hx withdrawal, the optA1 gpt strain displays a diminished ability to derepress the de novo purine biosynthesis genes, likely due to internal guanine accumulation. The impairment in fully inducing the purR regulon may be a contributing factor to the lethality of dGTP starvation. At later time points, and coinciding with cell lethality, strong induction of the SOS response was observed, supporting the concept of replication stress as a final cause of death. No evidence was observed in the starved cells for the participation of other stress responses, including the rpoS-mediated global stress response, reinforcing the lack of feedback of replication stress to the global metabolism of the cell. The genome-wide expression data also provide direct evidence for increased genome complexity during dGTP starvation, as a markedly increased gradient was observed for expression of genes located near the replication origin relative to those located toward the replication terminus.

IMPORTANCE

Control of the supply of the building blocks (deoxynucleoside triphosphates [dNTPs]) for DNA replication is important for ensuring genome integrity and cell viability. When cells are starved specifically for one of the four dNTPs, dGTP, the process of DNA replication is disturbed in a manner that can lead to eventual death. In the present study, we investigated the transcriptional changes in the bacterium E. coli during dGTP starvation. The results show increasing DNA replication stress with an increased time of starvation, as evidenced by induction of the bacterial SOS system, as well as a notable lack of induction of other stress responses that could have saved the cells from cell death by slowing down cell growth.