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

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.

sRNA roles in regulating transcriptional regulators: Lrp and SoxS regulation by sRNAs

Post-transcriptional regulation of transcription factors contributes to regulatory circuits. We created translational reporter fusions for multiple central regulators in Escherichia coli and examined the effect of Hfq-dependent non-coding RNAs on these fusions. This approach yields an 'RNA landscape,' identifying Hfq-dependent sRNAs that regulate a given fusion. No significant sRNA regulation of crp or fnr was detected. hns was regulated only by DsrA, as previously reported. Lrp and SoxS were both found to be regulated post-transcriptionally. Lrp, ' L: eucine-responsive R: egulatory P: rotein,' regulates genes involved in amino acid biosynthesis and catabolism and other cellular functions. sRNAs DsrA, MicF and GcvB each independently downregulate the lrp translational fusion, confirming previous reports for MicF and GcvB. MicF and DsrA interact with an overlapping site early in the lrp ORF, while GcvB acts upstream at two independent sites in the long lrp leader. Surprisingly, GcvB was found to be responsible for significant downregulation of lrp after oxidative stress; MicF also contributed. SoxS, an activator of genes used to combat oxidative stress, is negatively regulated by sRNA MgrR. This study demonstrates that while not all global regulators are subject to sRNA regulation, post-transcriptional control by sRNAs allows multiple environmental signals to affect synthesis of the transcriptional regulator.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNA synthetases (aaRSs) are modular enzymes globally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation. Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g., in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show huge structural plasticity related to function and limited idiosyncrasies that are kingdom or even species specific (e.g., the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS). Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably between distant groups such as Gram-positive and Gram-negative Bacteria. The review focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation, and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulated in last two decades is reviewed, showing how the field moved from essentially reductionist biology towards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRS paralogs (e.g., during cell wall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointed throughout the review and distinctive characteristics of bacterium-like synthetases from organelles are outlined.

Regulation of Expression of Uropathogenic Escherichia coli Nonfimbrial Adhesin TosA by PapB Homolog TosR in Conjunction with H-NS and Lrp

Urinary tract infections (UTIs) are a major burden to human health. The overwhelming majority of UTIs are caused by uropathogenic Escherichia coli (UPEC) strains. Unlike some pathogens, UPEC strains do not have a fixed core set of virulence and fitness factors but do have a variety of adhesins and regulatory pathways. One such UPEC adhesin is the nonfimbrial adhesin TosA, which mediates adherence to the epithelium of the upper urinary tract. The tos operon is AT rich, resides on pathogenicity island aspV, and is not expressed under laboratory conditions. Because of this, we hypothesized that tosA expression is silenced by H-NS. Lrp, based on its prominent function in the regulation of other adhesins, is also hypothesized to contribute to tos operon regulation. Using a variety of in vitro techniques, we mapped both the tos operon promoter and TosR binding sites. We have now identified TosR as a dual regulator of the tos operon, which could control the tos operon in association with H-NS and Lrp. H-NS is a negative regulator of the tos operon, and Lrp positively regulates the tos operon. Exogenous leucine also inhibits Lrp-mediated tos operon positive regulation. In addition, TosR binds to the pap operon, which encodes another important UPEC adhesin, P fimbria. Induction of TosR synthesis reduces production of P fimbria. These studies advance our knowledge of regulation of adhesin expression associated with uropathogen colonization of a host.

Immunization with individual proteins of the Lrp/AsnC family induces protection against Brucella melitensis 16M challenges in mice

Brucellosis is one of the most common zoonoses worldwide. Subunit vaccines are promising for the prevention of human brucellosis. In our previous protective antigen screening studies, we identified a new protective antigen, BMEI0357, which belongs to the Lrp/asnC protein family, a conserved transcriptional regulator in bacteria that is absent in eukaryotes. In the present study, the Brucella genome annotation was screened and a total of six proteins were identified as members of the Lrp/AsnC family. Lrp/AsnC proteins have two domains that are conserved among the family members. However, sequence similarities between these proteins ranged from 9 to 50%, indicating high sequence heterogeneity. To test whether proteins of this family have similar characteristics, all six proteins were cloned and expressed in Escherichia coli. The recombinant proteins were purified and their protective efficacy was evaluated in BALB/c mice challenged with Brucella melitensis 16M. The results show that all six Lrp/AsnC proteins could induce a protective immune response against Brucella melitensis 16M. Antibodies against the Lrp/AsnC proteins were detected in the immunized mice. However, levels of antibodies against these proteins were relatively variable in human brucellosis sera. Taken together, our results show that these six proteins of the Lrp/AsnC family in Brucella could induce protective immune responses in mice.

Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide

Biochemical network reconstructions have become popular tools in systems biology. Metabolicnetwork reconstructions are biochemically, genetically, and genomically (BiGG) structured databases of biochemical reactions and metabolites. They contain information such as exact reaction stoichiometry, reaction reversibility, and the relationships between genes, proteins, and reactions. Network reconstructions have been used extensively to study the phenotypic behavior of wild-type and mutant stains under a variety of conditions, linking genotypes with phenotypes. Such phenotypic simulations have allowed for the prediction of growth after genetic manipulations, prediction of growth phenotypes after adaptive evolution, and prediction of essential genes. Additionally, because network reconstructions are organism specific, they can be used to understand differences between organisms of species in a functional context.There are different types of reconstructions representing various types of biological networks (metabolic, regulatory, transcription/translation). This chapter serves as an introduction to metabolic and regulatory network reconstructions and models and gives a complete description of the core Escherichia coli metabolic model. This model can be analyzed in any computational format (such as MATLAB or Mathematica) based on the information given in this chapter. The core E. coli model is a small-scale model that can be used for educational purposes. It is meant to be used by senior undergraduate and first-year graduate students learning about constraint-based modeling and systems biology. This model has enough reactions and pathways to enable interesting and insightful calculations, but it is also simple enough that the results of such calculations can be understoodeasily.

Transcriptional Analysis of the Conjugal Transfer Genes of Rickettsia bellii RML 369-C

Rickettsia bellii is an obligate intracellular bacterium that is one of the few rickettsiae that encode a complete set of conjugative transfer (tra) genes involved in bacterial conjugation and has been shown to exhibit pili-like structures. The reductive genomes of rickettsiae beg the question whether the tra genes are nonfunctional or functioning to enhance the genetic plasticity and biology of rickettsiae. We characterized the transcriptional dynamics of R. bellii tra genes in comparison to genes transcribed stably and above the background level to understand when and at what levels the tra genes are active or whether the tra genes are degenerative. We determined that the best reference genes, out of 10 tested, were methionyl tRNA ligase (metG) or a combination of metG and ribonucleoside diphosphate reductase 2 subunit beta (nrdF), using statistical algorithms from two different programs: Normfinder and BestKeeper. To validate the use of metG with other rickettsial genes exhibiting variable transcriptional patterns we examined its use with sca2 and rickA, genes involved in actin based motility. Both were shown to be up-regulated at different times of replication in Vero cells, showing variable and stable transcription levels of rickA and sca2, respectively. traA_Ti was up-regulated at 72 hours post inoculation in the tick cell line ISE6, but showed no apparent changes in the monkey cell line Vero and mouse cell line L929. The transcription of tra genes was positively correlated with one another and up-regulated from 12 to 72 hours post inoculation (HPI) when compared to RBE_0422 (an inactivated transposase-derivative found within the tra cluster). Thus, the up-regulation of the tra genes indicated that the integrity and activity of each gene were intact and may facilitate the search for the optimal conditions necessary to demonstrate conjugation in rickettsiae.

Regulation of indole-3-acetic acid biosynthesis by branched-chain amino acids in Enterobacter cloacae UW5

The soil bacterium Enterobacter cloacae UW5 produces the rhizosphere signaling molecule indole-3-acetic acid (IAA) via the indolepyruvate pathway. Expression of indolepyruvate decarboxylase, a key pathway enzyme encoded by ipdC, is upregulated by the transcription factor TyrR in response to aromatic amino acids. Some members of the TyrR regulon may also be controlled by branched-chain amino acids and here we show that expression from the ipdC promoter and production of IAA are downregulated by valine, leucine and isoleucine. Regulation of the IAA synthesis pathway by both aromatic and branched-chain amino acids suggests a broader role for this pathway in bacterial physiology, beyond plant interactions.

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.

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

The regulators GadE, GadW and GadX (which we refer to as GadEWX) play a critical role in the transcriptional regulation of the glutamate-dependent acid resistance (GDAR) system in Escherichia coli K-12 MG1655. However, the genome-wide regulatory role of GadEWX is still unknown. Here we comprehensively reconstruct the genome-wide GadEWX transcriptional regulatory network and RpoS involvement in E. coli K-12 MG1655 under acidic stress. Integrative data analysis reveals that GadEWX regulons consist of 45 genes in 31 transcription units and 28 of these genes were associated with RpoS-binding sites. We demonstrate that GadEWX directly and coherently regulate several proton-generating/consuming enzymes with pairs of negative-feedback loops for pH homeostasis. In addition, GadEWX regulate genes with assorted functions, including molecular chaperones, acid resistance, stress response and other regulatory activities. These results show how GadEWX simultaneously coordinate many cellular processes to produce the overall response of E. coli to acid stress.

GadEWX regulons play a critical role in transcription regulation in response to acid stress. By reconstructing genome-wide GadEWX transcriptional network, here the authors show how GadEWX simultaneously coordinates many other cellular processes to produce the overall response of E. coli to acid stress.

Expanded roles of leucine-responsive regulatory protein in transcription regulation of the Escherichia coli genome: Genomic SELEX screening of the regulation targets

Leucine-responsive regulatory protein (Lrp) is a transcriptional regulator for the genes involved in transport, biosynthesis and catabolism of amino acids in Escherichia coli. In order to identify the whole set of genes under the direct control of Lrp, we performed Genomic SELEX screening and identified a total of 314 Lrp-binding sites on the E. coli genome. As a result, the regulation target of Lrp was predicted to expand from the hitherto identified genes for amino acid metabolism to a set of novel target genes for utilization of amino acids for protein synthesis, including tRNAs, aminoacyl-tRNA synthases and rRNAs. Northern blot analysis indicated alteration of mRNA levels for at least some novel targets, including the aminoacyl-tRNA synthetase genes. Phenotype MicroArray of the lrp mutant indicated significant alteration in utilization of amino acids and peptides, whilst metabolome analysis showed variations in the concentration of amino acids in the lrp mutant. From these two datasets we realized a reverse correlation between amino acid levels and cell growth rate: fast-growing cells contain low-level amino acids, whilst a high level of amino acids exists in slow-growing cells. Taken together, we propose that Lrp is a global regulator of transcription of a large number of the genes involved in not only amino acid transport and metabolism, but also amino acid utilization.

From ingestion to colonization: the influence of the host environment on regulation of the LEE encoded type III secretion system in enterohaemorrhagic Escherichia coli

Enterohaemorrhagic Escherichia coli (EHEC) binds to host tissue and intimately attaches to intestinal cells using a dedicated type III secretion system (T3SS). This complex multi-protein organelle is encoded within a large pathogenicity island called the locus of enterocyte effacement (LEE), which is subject to extensive regulatory control. Over the past 15 years we have gained a wealth of knowledge concerning how the LEE is regulated transcriptionally by specific, global and phage encoded regulators. More recently, significant advances have been made in our understanding of how specific signals, including host or microbiota derived metabolic products and various nutrient sources, can affect how the LEE-encoded T3SS is regulated. In this review we discuss regulation of the LEE, focusing on how these physiologically relevant signals are sensed and how they affect the expression of this major virulence factor. The implications for understanding the disease process by specific regulatory mechanisms are also discussed.

Cross talk between ABC transporter mRNAs via a target mRNA-derived sponge of the GcvB small RNA

There is an expanding list of examples by which one mRNA can posttranscriptionally influence the expression of others. This can involve RNA sponges that sequester regulatory RNAs of mRNAs in the same regulon, but the underlying molecular mechanism of such mRNA cross talk remains little understood. Here, we report sponge-mediated mRNA cross talk in the posttranscriptional network of GcvB, a conserved Hfq-dependent small RNA with one of the largest regulons known in bacteria. We show that mRNA decay from the gltIJKL locus encoding an amino acid ABC transporter generates a stable fragment (SroC) that base-pairs with GcvB. This interaction triggers the degradation of GcvB by RNase E, alleviating the GcvB-mediated mRNA repression of other amino acid-related transport and metabolic genes. Intriguingly, since the gltIJKL mRNA itself is a target of GcvB, the SroC sponge seems to enable both an internal feed-forward loop to activate its parental mRNA in cis and activation of many trans-encoded mRNAs in the same pathway. Disabling this mRNA cross talk affects bacterial growth when peptides are the sole carbon and nitrogen sources.

Modes of Escherichia coli Dps Interaction with DNA as Revealed by Atomic Force Microscopy

Multifunctional protein Dps plays an important role in iron assimilation and a crucial role in bacterial genome packaging. Its monomers form dodecameric spherical particles accumulating ~400 molecules of oxidized iron ions within the protein cavity and applying a flexible N-terminal ends of each subunit for interaction with DNA. Deposition of iron is a well-studied process by which cells remove toxic Fe²⁺ ions from the genetic material and store them in an easily accessible form. However, the mode of interaction with linear DNA remained mysterious and binary complexes with Dps have not been characterized so far. It is widely believed that Dps binds DNA without any sequence or structural preferences but several lines of evidence have demonstrated its ability to differentiate gene expression, which assumes certain specificity. Here we show that Dps has a different affinity for the two DNA fragments taken from the dps gene regulatory region. We found by atomic force microscopy that Dps predominantly occupies thermodynamically unstable ends of linear double-stranded DNA fragments and has high affinity to the central part of the branched DNA molecule self-assembled from three single-stranded oligonucleotides. It was proposed that Dps prefers binding to those regions in DNA that provide more contact pads for the triad of its DNA-binding bundle associated with one vertex of the protein globule. To our knowledge, this is the first study revealed the nucleoid protein with an affinity to branched DNA typical for genomic regions with direct and inverted repeats. As a ubiquitous feature of bacterial and eukaryotic genomes, such structural elements should be of particular care, but the protein system evolutionarily adapted for this function is not yet known, and we suggest Dps as a putative component of this system.

Nucleoid-associated proteins encoded on plasmids: Occurrence and mode of function

Nucleoid-associated proteins (NAPs) play a role in changing the shape of microbial DNA, making it more compact and affecting the regulation of transcriptional networks in host cells. Genes that encode NAPs include H-NS family proteins (H-NS, Ler, MvaT, BpH3, Bv3F, HvrA, and Lsr2), FIS, HU, IHF, Lrp, and NdpA, and are found in both microbial chromosomes and plasmid DNA. In the present study, NAP genes were distributed among 442 plasmids out of 4602 plasmid sequences, and many H-NS family proteins, and HU, IHF, Lrp, and NdpA were found in plasmids of Alpha-, Beta-, and Gammaproteobacteria, while HvrA, Lsr2, HU, and Lrp were found in other classes including Actinobacteria and Bacilli. Larger plasmids frequently carried multiple NAP genes. In addition, NAP genes were more frequently found in conjugative plasmids than non-transmissible plasmids. Several host cells carried the same types of H-NS family proteins on both their plasmids and chromosome(s), while this was not observed for other NAPs. Recent studies have shown that NAP genes on plasmids and chromosomes play important roles in the physical and regulatory integration of plasmids into the host cell.

Stochastic analysis of bistability in coherent mixed feedback loops combining transcriptional and posttranscriptional regulations

Mixed feedback loops combining transcriptional and posttranscriptional regulations are common in cellular regulatory networks. They consist of two genes, encoding a transcription factor and a small noncoding RNA (sRNA), which mutually regulate each other's expression. We present a theoretical and numerical study of coherent mixed feedback loops of this type, in which both regulations are negative. Under suitable conditions, these feedback loops are expected to exhibit bistability, namely, two stable states, one dominated by the transcriptional repressor and the other dominated by the sRNA. We use deterministic methods based on rate equation models, in order to identify the range of parameters in which bistability takes place. However, the deterministic models do not account for the finite lifetimes of the bistable states and the spontaneous, fluctuation-driven transitions between them. Therefore, we use stochastic methods to calculate the average lifetimes of the two states. It is found that these lifetimes strongly depend on rate coefficients such as the transcription rates of the transcriptional repressor and the sRNA. In particular, we show that the fraction of time the system spends in the sRNA-dominated state follows a monotonically decreasing sigmoid function of the transcriptional repressor transcription rate. The biological relevance of these results is discussed in the context of such mixed feedback loops in Escherichia coli. It is shown that the fluctuation-driven transitions and the dependence of some rate coefficients on the biological conditions enable the cells to switch to the state which is better suited for the existing conditions and to remain in that state as long as these conditions persist.

The architecture of ArgR-DNA complexes at the genome-scale in Escherichia coli

DNA-binding motifs that are recognized by transcription factors (TFs) have been well studied; however, challenges remain in determining the in vivo architecture of TF-DNA complexes on a genome-scale. Here, we determined the in vivo architecture of Escherichia coli arginine repressor (ArgR)-DNA complexes using high-throughput sequencing of exonuclease-treated chromatin-immunoprecipitated DNA (ChIP-exo). The ChIP-exo has a unique peak-pair pattern indicating 5′ and 3′ ends of ArgR-binding region. We identified 62 ArgR-binding loci, which were classified into three groups, comprising single, double and triple peak-pairs. Each peak-pair has a unique 93 base pair (bp)-long (±2 bp) ArgR-binding sequence containing two ARG boxes (39 bp) and residual sequences. Moreover, the three ArgR-binding modes defined by the position of the two ARG boxes indicate that DNA bends centered between the pair of ARG boxes facilitate the non-specific contacts between ArgR subunits and the residual sequences. Additionally, our approach may also reveal other fundamental structural features of TF-DNA interactions that have implications for studying genome-scale transcriptional regulatory networks.

Many means to a common end: the intricacies of (p)ppGpp metabolism and its control of bacterial homeostasis

In nearly all bacterial species examined so far, amino acid starvation triggers the rapid accumulation of the nucleotide second messenger (p)ppGpp, the effector of the stringent response. While for years the enzymes involved in (p)ppGpp metabolism and the significance of (p)ppGpp accumulation to stress survival were considered well defined, a recent surge of interest in the field has uncovered an unanticipated level of diversity in how bacteria metabolize and utilize (p)ppGpp to rapidly synchronize a variety of biological processes important for growth and stress survival. In addition to the classic activation of the stringent response, it has become evident that (p)ppGpp exerts differential effects on cell physiology in an incremental manner rather than simply acting as a biphasic switch that controls growth or stasis. Of particular interest is the intimate relationship of (p)ppGpp with persister cell formation and virulence, which has spurred the pursuit of (p)ppGpp inhibitors as a means to control recalcitrant infections. Here, we present an overview of the enzymes responsible for (p)ppGpp metabolism, elaborate on the intricacies that link basal production of (p)ppGpp to bacterial homeostasis, and discuss the implications of targeting (p)ppGpp synthesis as a means to disrupt long-term bacterial survival strategies.

Effective trapping of fruit flies with cultures of metabolically modified acetic acid bacteria

Acetoin in vinegar is an attractant to fruit flies when combined with acetic acid. To make vinegar more effective in attracting fruit flies with increased acetoin production, Komagataeibacter europaeus KGMA0119 was modified by specific gene disruption of the acetohydroxyacid isomeroreductase gene (ilvC). A previously constructed mutant lacking the putative ligand-sensing region in the leucine-responsive regulatory protein (KeLrp, encoded by Kelrp) was also used. The ilvC and Kelrp disruptants (KGMA5511 and KGMA7203, respectively) produced greater amounts of acetoin (KGMA5511, 0.11%; KGMA7203, 0.13%) than the wild-type strain KGMA0119 (0.069%). KGMA7203 produced a trace amount of isobutyric acid (0.007%), but the other strains did not. These strains produced approximately equal amounts of acetic acid (0.7%). The efficiency of fruit fly attraction was investigated with cultured Drosophila melanogaster. D. melanogaster flies (approximately 1,500) were released inside a cage (2.5 m by 2.5 m by 1.5 m) and were trapped with a device containing vinegar and a sticky sheet. The flies trapped on the sticky sheet were counted. The cell-free supernatant from KGMA7203 culture captured significantly more flies (19.36 to 36.96% of released flies) than did KGMA0119 (3.25 to 11.40%) and KGMA5511 (6.87 to 21.50%) cultures. Contrastingly, a 0.7% acetic acid solution containing acetoin (0.13%) and isobutyric acid (0.007%), which mimicked the KGMA7203 supernatant, captured significantly fewer flies (0.88 to 4.57%). Furthermore, the KGMA0119 supernatant with additional acetoin (0.13%) and isobutyric acid (0.007%) captured slightly more flies than the original KGMA0119 supernatant but fewer than the KGMA7203 supernatant, suggesting that the synergistic effects of acetic acid, acetoin, isobutyric acid, and unidentified metabolites achieved the efficient fly trapping of the KGMA7203 supernatant.

The metabolic regulator CodY links Listeria monocytogenes metabolism to virulence by directly activating the virulence regulatory gene prfA

Metabolic adaptations are critical to the ability of bacterial pathogens to grow within host cells and are normally preceded by sensing of host-specific metabolic signals, which in turn can influence the pathogen's virulence state. Previously, we reported that the intracellular bacterial pathogen Listeria monocytogenes responds to low availability of branched-chain amino acids (BCAAs) within mammalian cells by up-regulating both BCAA biosynthesis and virulence genes. The induction of virulence genes required the BCAA-responsive transcription regulator, CodY, but the molecular mechanism governing this mode of regulation was unclear. In this report, we demonstrate that CodY directly binds the coding sequence of the L. monocytogenes master virulence activator gene, prfA, 15 nt downstream of its start codon, and that this binding results in up-regulation of prfA transcription specifically under low concentrations of BCAA. Mutating this site abolished CodY binding and reduced prfA transcription in macrophages, and attenuated bacterial virulence in mice. Notably, the mutated binding site did not alter prfA transcription or PrfA activity under other conditions that are known to activate PrfA, such as during growth in the presence of glucose-1-phosphate. This study highlights the tight crosstalk between L. monocytogenes metabolism and virulence, while revealing novel features of CodY-mediated regulation.

Inference of expanded Lrp-like feast/famine transcription factor targets in a non-model organism using protein structure-based prediction

Widespread microbial genome sequencing presents an opportunity to understand the gene regulatory networks of non-model organisms. This requires knowledge of the binding sites for transcription factors whose DNA-binding properties are unknown or difficult to infer. We adapted a protein structure-based method to predict the specificities and putative regulons of homologous transcription factors across diverse species. As a proof-of-concept we predicted the specificities and transcriptional target genes of divergent archaeal feast/famine regulatory proteins, several of which are encoded in the genome of Halobacterium salinarum. This was validated by comparison to experimentally determined specificities for transcription factors in distantly related extremophiles, chromatin immunoprecipitation experiments, and cis-regulatory sequence conservation across eighteen related species of halobacteria. Through this analysis we were able to infer that Halobacterium salinarum employs a divergent local trans-regulatory strategy to regulate genes (carA and carB) involved in arginine and pyrimidine metabolism, whereas Escherichia coli employs an operon. The prediction of gene regulatory binding sites using structure-based methods is useful for the inference of gene regulatory relationships in new species that are otherwise difficult to infer.

Deciphering Fur transcriptional regulatory network highlights its complex role beyond iron metabolism in Escherichia coli

The ferric uptake regulator (Fur) plays a critical role in the transcriptional regulation of iron metabolism. However, the full regulatory potential of Fur remains undefined. Here we comprehensively reconstruct the Fur transcriptional regulatory network in Escherichia coli K-12 MG1655 in response to iron availability using genome-wide measurements (ChIP-exo and RNA-seq). Integrative data analysis reveals that a total of 81 genes in 42 transcription units are directly regulated by three different modes of Fur regulation, including apo- and holo-Fur activation and holo-Fur repression. We show that Fur connects iron transport and utilization enzymes with negative-feedback loop pairs for iron homeostasis. In addition, direct involvement of Fur in the regulation of DNA synthesis, energy metabolism, and biofilm development is found. These results show how Fur exhibits a comprehensive regulatory role affecting many fundamental cellular processes linked to iron metabolism in order to coordinate the overall response of E. coli to iron availability.

LeuO enhances butyrate-induced virulence expression through a positive regulatory loop in enterohaemorrhagic Escherichia coli

Enterohaemorrhagic Escherichia coli (EHEC) causes bloody diarrhoea and other severe symptoms such as haemorrhagic uraemic syndrome. The expression of virulence genes on the locus for enterocyte effacement (LEE) and associated genes is regulated by a variety of factors, including transcriptional regulators and environmental signals. Butyrate, one of the major short-chain fatty acids present in the intestine, enhances expression of LEE genes and flagella biosynthesis genes in EHEC O157:H7, resulting in increased bacterial adherence and motility. Here, we show that expression of the leuO gene, which encodes a LysR-type transcriptional regulator, is enhanced by butyrate via Lrp, which is also necessary for butyrate-induced responses of LEE genes. LeuO expression induces prolonged activation of the promoter of LEE1 operon, including the ler gene, as well as virulence mechanisms such as microcolony formation. Activation of the LEE1 promoter by LeuO depends on another regulator, called Pch. The response of the leuO promoter to butyrate requires two virulence regulators, Pch and Ler, in addition to Lrp. Pch, Ler and Lrp bind the upstream region of the leuO promoter. Thus, leuO is involved in butyrate-enhanced expression of LEE genes through a positive feedback mechanism, but its expression and action on the LEE1 promoter are dependent on the virulence regulators Pch and Ler.

Leucine-responsive regulatory protein Lrp and PapI homologues influence phase variation of CS31A fimbriae

CS31A, a K88-related surface antigen specified by the clp operon, is a member of the type P family of adhesive factors and plays a key role in the establishment of disease caused by septicemic and enterotoxigenic Escherichia coli strains. Its expression is under the control of methylation-dependent transcriptional regulation, for which the leucine-responsive regulatory protein (Lrp) is essential. CS31A is preferentially in the OFF state and exhibits distinct regulatory features compared to the regulation of other P family members. In the present study, surface plasmon resonance and DNase I protection assays showed that Lrp binds to the distal moiety of the clp regulatory region with low micromolar affinity compared to its binding to the proximal moiety, which exhibits stronger, nanomolar affinity. The complex formation was also influenced by the addition of PapI or FooI, which increased the affinity of Lrp for the clp distal and proximal regions and was required to induce phase variation. The influence of PapI or FooI, however, was predominantly associated with a more complete shutdown of clp expression, in contrast to what has previously been observed with AfaF (a PapI ortholog). Taken together, these results suggest that the preferential OFF state observed in CS31A cells is mainly due to the weak interaction of the leucine-responsive regulatory protein with the clp distal region and that the PapI homolog favors the OFF phase. Within the large repertoire of fimbrial variants in the P family, our study illustrates that having a fimbrial operon that lacks its own PapI ortholog allows it to be more flexibly regulated by other orthologs in the cell.

Determining the control circuitry of redox metabolism at the genome-scale

Determining how facultative anaerobic organisms sense and direct cellular responses to electron acceptor availability has been a subject of intense study. However, even in the model organism Escherichia coli, established mechanisms only explain a small fraction of the hundreds of genes that are regulated during electron acceptor shifts. Here we propose a qualitative model that accounts for the full breadth of regulated genes by detailing how two global transcription factors (TFs), ArcA and Fnr of E. coli, sense key metabolic redox ratios and act on a genome-wide basis to regulate anabolic, catabolic, and energy generation pathways. We first fill gaps in our knowledge of this transcriptional regulatory network by carrying out ChIP-chip and gene expression experiments to identify 463 regulatory events. We then interfaced this reconstructed regulatory network with a highly curated genome-scale metabolic model to show that ArcA and Fnr regulate >80% of total metabolic flux and 96% of differential gene expression across fermentative and nitrate respiratory conditions. Based on the data, we propose a feedforward with feedback trim regulatory scheme, given the extensive repression of catabolic genes by ArcA and extensive activation of chemiosmotic genes by Fnr. We further corroborated this regulatory scheme by showing a 0.71 r(2) (p<1e-6) correlation between changes in metabolic flux and changes in regulatory activity across fermentative and nitrate respiratory conditions. Finally, we are able to relate the proposed model to a wealth of previously generated data by contextualizing the existing transcriptional regulatory network.

Coordination of microbial metabolism

Beyond fuelling cellular activities with building blocks and energy, metabolism also integrates environmental conditions into intracellular signals. The underlying regulatory network is complex and multifaceted: it ranges from slow interactions, such as changing gene expression, to rapid ones, such as the modulation of protein activity via post-translational modification or the allosteric binding of small molecules. In this Review, we outline the coordination of common metabolic tasks, including nutrient uptake, central metabolism, the generation of energy, the supply of amino acids and protein synthesis. Increasingly, a set of key metabolites is recognized to control individual regulatory circuits, which carry out specific functions of information input and regulatory output. Such a modular view of microbial metabolism facilitates an intuitive understanding of the molecular mechanisms that underlie cellular decision making.

Genome-scale reconstruction of the sigma factor network in Escherichia coli: topology and functional states

Background

At the beginning of the transcription process, the RNA polymerase (RNAP) core enzyme requires a σ-factor to recognize the genomic location at which the process initiates. Although the crucial role of σ-factors has long been appreciated and characterized for many individual promoters, we do not yet have a genome-scale assessment of their function.

Results

Using multiple genome-scale measurements, we elucidated the network of σ-factor and promoter interactions in Escherichia coli. The reconstructed network includes 4,724 σ-factor-specific promoters corresponding to transcription units (TUs), representing an increase of more than 300% over what has been previously reported. The reconstructed network was used to investigate competition between alternative σ-factors (the σ⁷⁰ and σ³⁸ regulons), confirming the competition model of σ substitution and negative regulation by alternative σ-factors. Comparison with σ-factor binding in Klebsiella pneumoniae showed that transcriptional regulation of conserved genes in closely related species is unexpectedly divergent.

Conclusions

The reconstructed network reveals the regulatory complexity of the promoter architecture in prokaryotic genomes, and opens a path to the direct determination of the systems biology of their transcriptional regulatory networks.

PortEco: a resource for exploring bacterial biology through high-throughput data and analysis tools

PortEco (http://porteco.org) aims to collect, curate and provide data and analysis tools to support basic biological research in Escherichia coli (and eventually other bacterial systems). PortEco is implemented as a ‘virtual’ model organism database that provides a single unified interface to the user, while integrating information from a variety of sources. The main focus of PortEco is to enable broad use of the growing number of high-throughput experiments available for E. coli, and to leverage community annotation through the EcoliWiki and GONUTS systems. Currently, PortEco includes curated data from hundreds of genome-wide RNA expression studies, from high-throughput phenotyping of single-gene knockouts under hundreds of annotated conditions, from chromatin immunoprecipitation experiments for tens of different DNA-binding factors and from ribosome profiling experiments that yield insights into protein expression. Conditions have been annotated with a consistent vocabulary, and data have been consistently normalized to enable users to find, compare and interpret relevant experiments. PortEco includes tools for data analysis, including clustering, enrichment analysis and exploration via genome browsers. PortEco search and data analysis tools are extensively linked to the curated gene, metabolic pathway and regulation content at its sister site, EcoCyc.

The expression of stlA in Photorhabdus luminescens is controlled by nutrient limitation

Photorhabdus is a genus of Gram-negative entomopathogenic bacteria that also maintain a mutualistic association with nematodes from the family Heterorhabditis. Photorhabdus has an extensive secondary metabolism that is required for the interaction between the bacteria and the nematode. A major component of this secondary metabolism is a stilbene molecule, called ST. The first step in ST biosynthesis is the non-oxidative deamination of phenylalanine resulting in the production of cinnamic acid. This reaction is catalyzed by phenylalanine-ammonium lyase, an enzyme encoded by the stlA gene. In this study we show, using a stlA-gfp transcriptional fusion, that the expression of stlA is regulated by nutrient limitation through a regulatory network that involves at least 3 regulators. We show that TyrR, a LysR-type transcriptional regulator that regulates gene expression in response to aromatic amino acids in E. coli, is absolutely required for stlA expression. We also show that stlA expression is modulated by σ^S and Lrp, regulators that are implicated in the regulation of the response to nutrient limitation in other bacteria. This work is the first that describes pathway-specific regulation of secondary metabolism in Photorhabdus and, therefore, our study provides an initial insight into the complex regulatory network that controls secondary metabolism, and therefore mutualism, in this model organism.

The Escherichia coli Cpx envelope stress response regulates genes of diverse function that impact antibiotic resistance and membrane integrity

The Cpx envelope stress response mediates adaptation to stresses that cause envelope protein misfolding. Adaptation is partly conferred through increased expression of protein folding and degradation factors. The Cpx response also plays a conserved role in the regulation of virulence determinant expression and impacts antibiotic resistance. We sought to identify adaptive mechanisms that may be involved in these important functions by characterizing changes in the transcriptome of two different Escherichia coli strains when the Cpx response is induced. We show that, while there is considerable strain- and condition-specific variability in the Cpx response, the regulon is enriched for proteins and functions that are inner membrane associated under all conditions. Genes that were changed by Cpx pathway induction under all conditions were involved in a number of cellular functions and included several intergenic regions, suggesting that posttranscriptional regulation is important during Cpx-mediated adaptation. Some Cpx-regulated genes are centrally involved in energetics and play a role in antibiotic resistance. We show that a number of small, uncharacterized envelope proteins are Cpx regulated and at least two of these affect phenotypes associated with membrane integrity. Altogether, our work suggests new mechanisms of Cpx-mediated envelope stress adaptation and antibiotic resistance.

Phylogenetic group-associated differences in regulation of the common colonization factor Mat fimbria in Escherichia coli

Heterogeneity of cell population is a key component behind the evolutionary success of Escherichia coli. The heterogeneity supports species adaptation and mainly results from lateral gene transfer. Adaptation may also involve genomic alterations that affect regulation of conserved genes. Here we analysed regulation of the mat (or ecp) genes that encode a conserved fimbrial adhesin of E. coli. We found that the differential and temperature-sensitive expression control of the mat operon is dependent on mat promoter polymorphism and closely linked to phylogenetic grouping of E. coli. In the mat promoter lineage favouring fimbriae expression, the mat operon-encoded regulator MatA forms a positive feedback loop that overcomes the repression by H-NS and stabilizes the fimbrillin mRNA under low growth temperature, acidic pH or elevated levels of acetate. The study exemplifies phylogenetic group-associated expression of a highly common surface organelle in E. coli.

Expanded target and cofactor repertoire for the transcriptional activator LysM from Sulfolobus

Previously, Lrp-like transcriptional regulator LysM from the hyperthermoacidophilic crenarchaeon Sulfolobus solfataricus was proposed to have a single target, the lysWXJK operon of lysine biosynthesis, and a single effector molecule, l-lysine. Here we identify ∼70 novel binding sites for LysM in the S. solfataricus genome with a LysM-specific nanobody-based chromatin immunoprecipitation assay coupled to microarray hybridization (ChIP-chip) and in silico target site prediction using an energy-based position weight matrix, and validate these findings with in vitro binding. LysM binds to intergenic and coding regions, including promoters of various amino acid biosynthesis and transport genes. We confirm that l-lysine is the most potent effector molecule that reduces, but does not completely abolish, LysM binding, and show that several other amino acids and derivatives, including d-lysine, l-arginine, l-homoarginine, l-glutamine and l-methionine and branched-chain amino acids l-leucine, l-isoleucine and l-valine, significantly affect DNA-binding properties of LysM. Therefore, it appears from this study that LysM is a much more versatile regulator than previously thought, and that it uses a variety of amino acids to sense nutritional quality of the environment and to modulate expression of the metabolic machinery of Sulfolobus accordingly.

Sa-Lrp from Sulfolobus acidocaldarius is a versatile, glutamine-responsive, and architectural transcriptional regulator

Sa-Lrp is a member of the leucine-responsive regulatory protein (Lrp)-like family of transcriptional regulators in Sulfolobus acidocaldarius. Previously, we demonstrated the binding of Sa-Lrp to the control region of its own gene in vitro. However, the function and cofactor of Sa-Lrp remained an enigma. In this work, we demonstrate that glutamine is the cofactor of Sa-Lrp by inducing the formation of octamers and increasing the DNA-binding affinity and sequence specificity. In vitro protein-DNA interaction assays indicate that Sa-Lrp binds to promoter regions of genes with a variety of functions including ammonia assimilation, transcriptional control, and UV-induced pili synthesis. DNA binding occurs with a specific affinity for AT-rich binding sites, and the protein induces DNA bending and wrapping upon binding, indicating an architectural role of the regulator. Furthermore, by analyzing an Sa-lrp deletion mutant, we demonstrate that the protein affects transcription of some of the genes of which the promoter region is targeted and that it is an important determinant of the cellular aggregation phenotype. Taking all these results into account, we conclude that Sa-Lrp is a glutamine-responsive global transcriptional regulator with an additional architectural role.

PreCisIon: PREdiction of CIS-regulatory elements improved by gene's positION

Conventional approaches to predict transcriptional regulatory interactions usually rely on the definition of a shared motif sequence on the target genes of a transcription factor (TF). These efforts have been frustrated by the limited availability and accuracy of TF binding site motifs, usually represented as position-specific scoring matrices, which may match large numbers of sites and produce an unreliable list of target genes. To improve the prediction of binding sites, we propose to additionally use the unrelated knowledge of the genome layout. Indeed, it has been shown that co-regulated genes tend to be either neighbors or periodically spaced along the whole chromosome. This study demonstrates that respective gene positioning carries significant information. This novel type of information is combined with traditional sequence information by a machine learning algorithm called PreCisIon. To optimize this combination, PreCisIon builds a strong gene target classifier by adaptively combining weak classifiers based on either local binding sequence or global gene position. This strategy generically paves the way to the optimized incorporation of any future advances in gene target prediction based on local sequence, genome layout or on novel criteria. With the current state of the art, PreCisIon consistently improves methods based on sequence information only. This is shown by implementing a cross-validation analysis of the 20 major TFs from two phylogenetically remote model organisms. For Bacillus subtilis and Escherichia coli, respectively, PreCisIon achieves on average an area under the receiver operating characteristic curve of 70 and 60%, a sensitivity of 80 and 70% and a specificity of 60 and 56%. The newly predicted gene targets are demonstrated to be functionally consistent with previously known targets, as assessed by analysis of Gene Ontology enrichment or of the relevant literature and databases.

Improved predictions of transcription factor binding sites using physicochemical features of DNA

Typical approaches for predicting transcription factor binding sites (TFBSs) involve use of a position-specific weight matrix (PWM) to statistically characterize the sequences of the known sites. Recently, an alternative physicochemical approach, called SiteSleuth, was proposed. In this approach, a linear support vector machine (SVM) classifier is trained to distinguish TFBSs from background sequences based on local chemical and structural features of DNA. SiteSleuth appears to generally perform better than PWM-based methods. Here, we improve the SiteSleuth approach by considering both new physicochemical features and algorithmic modifications. New features are derived from Gibbs energies of amino acid-DNA interactions and hydroxyl radical cleavage profiles of DNA. Algorithmic modifications consist of inclusion of a feature selection step, use of a nonlinear kernel in the SVM classifier, and use of a consensus-based post-processing step for predictions. We also considered SVM classification based on letter features alone to distinguish performance gains from use of SVM-based models versus use of physicochemical features. The accuracy of each of the variant methods considered was assessed by cross validation using data available in the RegulonDB database for 54 Escherichia coli TFs, as well as by experimental validation using published ChIP-chip data available for Fis and Lrp.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNAsynthetases (aaRSs) are modular enzymesglobally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation.Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g.,in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show hugestructural plasticity related to function andlimited idiosyncrasies that are kingdom or even speciesspecific (e.g.,the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS).Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably betweendistant groups such as Gram-positive and Gram-negative Bacteria.Thereview focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation,and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulatedin last two decades is reviewed,showing how thefield moved from essentially reductionist biologytowards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRSparalogs (e.g., during cellwall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointedthroughout the reviewand distinctive characteristics of bacterium-like synthetases from organelles are outlined.

The extracytoplasmic function sigma factor EcfO protects Bacteroides fragilis against oxidative stress

The anaerobe Bacteroides fragilis is a highly aerotolerant, opportunistic pathogen that is an important component of the human intestinal microbiota. Aerotolerance has been linked to a robust oxidative stress response, which in turn is necessary for maximal virulence in a mouse intra-abdominal abscess model. During oxidative stress, there is a dynamic change in gene expression that encompasses a third of the genome, but there is a paucity of information on factors that control this response. A large number of transcription regulators, including about 14 extracytoplasmic function (ECF) sigma factors, are affected by oxidative stress, and one of these, EcfO, was used as a model of ECF sigma factor activity during stress. Genetic and biochemical experiments showed that EcfO was located in an operon with a structurally unique anti-sigma factor, Reo. Cells deleted for EcfO were impaired during exposure to oxygen or other forms of oxidative stress, whereas reo mutants were more resistant to stress. Protein-protein interaction experiments demonstrated that Reo directly interacts with and regulates the activity of EcfO. Expression microarray and chromatin affinity precipitation assays were used to identify target genes regulated by EcfO, and an EcfO recognition sequence was identified. The results revealed that EcfO controls a regulon of novel lipoproteins whose distribution in nature is restricted to members of the Bacteroidetes phylum.