Identifying global regulators in transcriptional regulatory networks in bacteria

Biological computation through recurrence

One of the defining features of living systems is their adaptability to changing environmental conditions. This requires organisms to extract temporal and spatial features of their environment, and use that information to compute the appropriate response. In the last two decades, a growing body of work, mainly coming from the machine learning and computational neuroscience fields, has shown that such complex information processing can be performed by recurrent networks. Temporal computations arise in these networks through the interplay between the external stimuli and the network's internal state. In this article we review our current understanding of how recurrent networks can be used by biological systems, from cells to brains, for complex information processing. Rather than focusing on sophisticated, artificial recurrent architectures such as long short-term memory (LSTM) networks, here we concentrate on simpler network structures and learning algorithms that can be expected to have been found by evolution. We also review studies showing evidence of naturally occurring recurrent networks in living organisms. Lastly, we discuss some relevant evolutionary aspects concerning the emergence of this natural computation paradigm.

A library of reporters of the global regulators of gene expression in Escherichia coli

The topology of the transcription factor network (TFN) of Escherichia coli is far from uniform, with 22 global regulator (GR) proteins controlling one-third of all genes. So far, their production rates cannot be tracked by comparable fluorescent proteins. We developed a library of fluorescent reporters for 16 GRs for this purpose. Each consists of a single-copy plasmid coding for green fluorescent protein (GFP) fused to the full-length copy of the native promoter. We tracked their activity in exponential and stationary growth, as well as under weak and strong stresses. We show that the reporters have high sensitivity and specificity to all stresses tested and detect single-cell variability in transcription rates. Given the influence of GRs on the TFN, we expect that the new library will contribute to dissecting global transcriptional stress-response programs of E. coli. Moreover, the library can be invaluable in bioindustrial applications that tune those programs to, instead of cell growth, favor productivity while reducing energy consumption.IMPORTANCECells contain thousands of genes. Many genes are involved in the control of cellular activities. Some activities require a few hundred genes to run largely synchronous transcriptional programs. To achieve this, cells have evolved global regulator (GR) proteins that can influence hundreds of genes simultaneously. We have engineered a library of Escherichia coli strains to track the levels over time of these, phenotypically critical, GRs. Each strain has a single-copy plasmid coding for a fast-maturing green fluorescent protein whose transcription is controlled by a copy of the natural GR promoter. By allowing the tracking of GR levels, with sensitivity and specificity, this library should become of wide use in scientific research on bacterial gene expression (from molecular to synthetic biology) and, later, be used in applications in therapeutics and bioindustries.

Mutations upstream from sdaC and malT in Escherichia coli uncover a complex interplay between the cAMP receptor protein and different sigma factors

In Escherichia coli, one of the best understood microorganisms, much can still be learned about the basic interactions between transcription factors and promoters. When a cAMP-deficient cya mutant is supplied with maltose as the main carbon source, mutations develop upstream from the two genes malT and sdaC. Here, we explore the regulation of the two promoters, using fluorescence-based genetic reporters in combination with both spontaneously evolved and systematically engineered cis-acting mutations. We show that in the cya mutant, regulation of malT and sdaC evolves toward cAMP-independence and increased expression in the stationary phase. Furthermore, we show that the location of the cAMP receptor protein (Crp) binding site upstream of malT is important for alternative sigma factor usage. This provides new insights into the architecture of bacterial promoters and the global interplay between Crp and sigma factors in different growth phases.IMPORTANCEThis work provides new general insights into (1) the architecture of bacterial promoters, (2) the importance of the location of Class I Crp-dependent promoters, and (3) the global interplay between Crp and sigma factors in different growth phases.

The Brucella abortus two-component system response regulator BvrR binds to three DNA regulatory boxes in the upstream region of omp25

Brucella abortus is a facultative extracellular-intracellular bacterial zoonotic pathogen worldwide. It is also a major cause of abortion in bovines, generating economic losses. The two-component regulatory system BvrR/BvrS modulates the expression of genes required to transition from extracellular to intracellular lifestyles. However, few regulatory regions of BvrR direct target genes have been studied. In this study, we characterized the regulatory region of omp25, a gene encoding an outer membrane protein that is positively regulated by TCS BvrR/BvrS. By omp25-lacZ reporter fusions and β-galactosidase activity assays, we found that the region between-262 and + 127 is necessary for transcriptional activity, particularly a 111-bp long fragment located from-262 to -152. In addition, we demonstrated the binding of P-BvrR to three sites within the -140 to +1 region. Two of these sites were delimited between -18 to +1 and - 99 to -76 by DNase I footprinting and called DNA regulatory boxes 1 and 2, respectively. The third binding site (box 3) was delimited from -140 to -122 by combining EMSA and fluorescence anisotropy results. A molecular docking analysis with HDOCK predicted BvrR-DNA interactions between 11, 13, and 12 amino acid residue-nucleotide pairs in boxes 1, 2, and 3, respectively. A manual sequence alignment of the three regulatory boxes revealed the presence of inverted and non-inverted repeats of five to eight nucleotides, partially matching DNA binding motifs previously described for BvrR. We propose that P-BvrR binds directly to up to three regulatory boxes and probably interacts with other transcription factors to regulate omp25 expression. This gene regulation model could apply to other BvrR target genes and to orthologs of the TCS BvrR/BvrS and Omp25 in phylogenetically closed Rhizobiales.

Retron-mediated multiplex genome editing and continuous evolution in Escherichia coli

While there are several genome editing techniques available, few are suitable for dynamic and simultaneous mutagenesis of arbitrary targeted sequences in prokaryotes. Here, to address these limitations, we present a versatile and multiplex retron-mediated genome editing system (REGES). First, through systematic optimization of REGES, we achieve efficiency of ∼100%, 85 ± 3%, 69 ± 14% and 25 ± 14% for single-, double-, triple- and quadruple-locus genome editing, respectively. In addition, we employ REGES to generate pooled and barcoded variant libraries with degenerate RBS sequences to fine-tune the expression level of endogenous and exogenous genes, such as transcriptional factors to improve ethanol tolerance and biotin biosynthesis. Finally, we demonstrate REGES-mediated continuous in vivo protein evolution, by combining retron, polymerase-mediated base editing and error-prone transcription. By these case studies, we demonstrate REGES as a powerful multiplex genome editing and continuous evolution tool with broad applications in synthetic biology and metabolic engineering.

Transcription factor expression levels and environmental signals constrain transcription factor innovation

Evolutionary innovation of transcription factors frequently drives phenotypic diversification and adaptation to environmental change. Transcription factors can gain or lose connections to target genes, resulting in novel regulatory responses and phenotypes. However the frequency of functional adaptation varies between different regulators, even when they are closely related. To identify factors influencing propensity for innovation, we utilise a Pseudomonas fluorescens SBW25 strain rendered incapable of flagellar mediated motility in soft-agar plates via deletion of the flagellar master regulator (fleQ ). This bacterium can evolve to rescue flagellar motility via gene regulatory network rewiring of an alternative transcription factor to rescue activity of FleQ. Previously, we have identified two members (out of 22) of the RpoN-dependent enhancer binding protein (RpoN-EBP) family of transcription factors (NtrC and PFLU1132) that are capable of innovating in this way. These two transcription factors rescue motility repeatably and reliably in a strict hierarchy – with NtrC the only route in a ∆fleQ background, and PFLU1132 the only route in a ∆fleQ ∆ntrC background. However, why other members in the same transcription factor family have not been observed to rescue flagellar activity is unclear. Previous work shows that protein homology cannot explain this pattern within the protein family (RpoN-EBPs), and mutations in strains that rescued motility suggested high levels of transcription factor expression and activation drive innovation. We predict that mutations that increase expression of the transcription factor are vital to unlock evolutionary potential for innovation. Here, we construct titratable expression mutant lines for 11 of the RpoN-EBPs in P. fluorescens . We show that in five additional RpoN-EBPs (FleR, HbcR, GcsR, DctD, AauR and PFLU2209), high expression levels result in different mutations conferring motility rescue, suggesting alternative rewiring pathways. Our results indicate that expression levels (and not protein homology) of RpoN-EBPs are a key constraining factor in determining evolutionary potential for innovation. This suggests that transcription factors that can achieve high expression through few mutational changes, or transcription factors that are active in the selective environment, are more likely to innovate and contribute to adaptive gene regulatory network evolution.

What remains from living cells in bacterial lysate-based cell-free systems

Because they mimic cells while offering an accessible and controllable environment, lysate-based cell-free systems (CFS) have emerged as valuable biotechnology tools for synthetic biology. Historically used to uncover fundamental mechanisms of life, CFS are nowadays used for a multitude of purposes, including protein production and prototyping of synthetic circuits. Despite the conservation of fundamental functions in CFS like transcription and translation, RNAs and certain membrane-embedded or membrane-bound proteins of the host cell are lost when preparing the lysate. As a result, CFS largely lack some essential properties of living cells, such as the ability to adapt to changing conditions, to maintain homeostasis and spatial organization. Regardless of the application, shedding light on the black-box of the bacterial lysate is necessary to fully exploit the potential of CFS. Most measurements of the activity of synthetic circuits in CFS and in vivo show significant correlations because these only require processes that are preserved in CFS, like transcription and translation. However, prototyping circuits of higher complexity that require functions that are lost in CFS (cell adaptation, homeostasis, spatial organization) will not show such a good correlation with in vivo conditions. Both for prototyping circuits of higher complexity and for building artificial cells, the cell-free community has developed devices to reconstruct cellular functions. This mini-review compares bacterial CFS to living cells, focusing on functional and cellular process differences and the latest developments in restoring lost functions through complementation of the lysate or device engineering.

Quantitative model for genome-wide cyclic AMP receptor protein binding site identification and characteristic analysis

Cyclic AMP receptor proteins (CRPs) are important transcription regulators in many species. The prediction of CRP-binding sites was mainly based on position-weighted matrixes (PWMs). Traditional prediction methods only considered known binding motifs, and their ability to discover inflexible binding patterns was limited. Thus, a novel CRP-binding site prediction model called CRPBSFinder was developed in this research, which combined the hidden Markov model, knowledge-based PWMs and structure-based binding affinity matrixes. We trained this model using validated CRP-binding data from Escherichia coli and evaluated it with computational and experimental methods. The result shows that the model not only can provide higher prediction performance than a classic method but also quantitatively indicates the binding affinity of transcription factor binding sites by prediction scores. The prediction result included not only the most knowns regulated genes but also 1089 novel CRP-regulated genes. The major regulatory roles of CRPs were divided into four classes: carbohydrate metabolism, organic acid metabolism, nitrogen compound metabolism and cellular transport. Several novel functions were also discovered, including heterocycle metabolic and response to stimulus. Based on the functional similarity of homologous CRPs, we applied the model to 35 other species. The prediction tool and the prediction results are online and are available at: https://awi.cuhk.edu.cn/∼CRPBSFinder.

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.

ReViTA: A novel in vitro transcription system to study gene regulation

ReViTA (Reverse in VitroTranscription Assay) is a novel in vitro transcription-based method to study gene expression under the regulation of specific transcription factors. The ReViTA system uses a plasmid with a control sequence, the promoter region of the studied gene, the transcription factor of interest, and an RNA polymerase saturated with σ⁷⁰. The main objective of this study was to evaluate the method; thus, as a proof of concept, two different transcription factors were used, a transcriptional inducer, AlgR, and a repressor, LexA, from Pseudomonas aeruginosa. After the promoters were incubated with the transcription factors, the plasmid was transcribed into RNA and reverse transcribed to cDNA. Gene expression was measured using qRTPCR. Using the ReViTA plasmid, transcription induction of 55% was observed when AlgR protein was added and a 27% transcription reduction with the repressor LexA, compared with the samples without transcription factors. The results demonstrated the correct functioning of ReViTA as a novel method to study transcription factors and gene expression. Thus, ReViTA could be a rapid and accessible in vitro method to evaluate genes and regulators of various species.

The prophage-encoded transcriptional regulator AppY has pleiotropic effects on E. coli physiology

Bacterial genome diversity is influenced by prophages, which are viral genomes integrated into the bacterial chromosome. Most prophage genes are silent but those that are expressed can provide unexpected properties to their host. Using as a model E. coli K-12 that carries 9 defective prophages in its genome, we aimed at highlighting the impact of genes encoded by prophages on host physiology. We focused our work on AppY, a transcriptional regulator encoded on the DLP12 prophage. By performing RNA-Seq experiments, we showed that AppY production modulates the expression of more than 200 genes. Among them, 11 were identified by ChIP-Seq as direct AppY targets. AppY directly and positively regulates several genes involved in the acid stress response including the master regulator gene gadE but also nhaR and gadY, two genes important for biofilm formation. Moreover, AppY indirectly and negatively impacts bacterial motility by favoring the degradation of FlhDC, the master regulator of the flagella biosynthesis. As a consequence of these regulatory effects, AppY increases acid stress resistance and biofilm formation while also causing a strong defect in motility. Our research shed light on the importance to consider the genetic interactions occurring between prophages and bacteria to fully understand bacterial physiology. It also highlights how a prophage-encoded transcriptional regulator integrates in a complex manner into the host regulatory network and how it benefits its host, allowing it to cope with changing environmental conditions.

Structural and Functional Characterization of Rv0792c from Mycobacterium tuberculosis: Identifying Small Molecule Inhibitor against HutC Protein

In order to adapt in host tissues, microbial pathogens regulate their gene expression through a variety of transcription factors. Here, we have functionally characterized Rv0792c, a HutC homolog from Mycobacterium tuberculosis. In comparison to the parental strain, a strain of M. tuberculosis with a Rv0792c mutant was compromised for survival upon exposure to oxidative stress and infection in guinea pigs. RNA sequencing analysis revealed that Rv0792c regulates the expression of genes involved in stress adaptation and virulence of M. tuberculosis. Solution small-angle X-ray scattering (SAXS) data-steered model building confirmed that the C-terminal region plays a pivotal role in dimer formation. Systematic evolution of ligands by exponential enrichment (SELEX) resulted in the identification of single-strand DNA (ssDNA) aptamers that can be used as a tool to identify small-molecule inhibitors targeting Rv0792c. Using SELEX and SAXS data-based modeling, we identified residues essential for Rv0792c's aptamer binding activity. In this study, we also identified I-OMe-Tyrphostin as an inhibitor of Rv0792c's aptamer and DNA binding activity. The identified small molecule reduced the growth of intracellular M. tuberculosis in macrophages. The present study thus provides a detailed shape-function characterization of a HutC family of transcription factor from M. tuberculosis. IMPORTANCE Prokaryotes encode a large number of GntR family transcription factors that are involved in various fundamental biological processes, including stress adaptation and pathogenesis. Here, we investigated the structural and functional role of Rv0792c, a HutC homolog from M. tuberculosis. We demonstrated that Rv0792c is essential for M. tuberculosis to adapt to oxidative stress and establish disease in guinea pigs. Using a systematic evolution of ligands by exponential enrichment (SELEX) approach, we identified ssDNA aptamers from a random ssDNA library that bound to Rv0792c protein. These aptamers were thoroughly characterized using biochemical and biophysical assays. Using SAXS, we determined the structural model of Rv0792c in both the presence and absence of the aptamers. Further, using a combination of SELEX and SAXS methodologies, we identified I-OMe-Tyrphostin as a potential inhibitor of Rv0792c. Here we provide a detailed functional characterization of a transcription factor belonging to the HutC family from M. tuberculosis.

A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli

Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.

Novel organisation and regulation of the pic promoter from enteroaggregative and uropathogenic Escherichia coli

The serine protease autotransporters of the Enterobacteriaceae (SPATEs) are a large family of virulence factors commonly found in enteric bacteria. These secreted virulence factors have diverse functions during bacterial infection, including adhesion, aggregation and cell toxicity. One such SPATE, the Pic mucinase (protein involved in colonisation) cleaves mucin, allowing enteric bacterial cells to utilise mucin as a carbon source and to penetrate the gut mucus lining, thereby increasing mucosal colonisation. The pic gene is widely distributed within the Enterobacteriaceae, being found in human pathogens, such as enteroaggregative Escherichia coli (EAEC), uropathogenic E. coli (UPEC) and Shigella flexneri 2a. As the pic promoter regions from EAEC strain 042 and UPEC strain CFT073 differ, we have investigated the regulation of each promoter. Here, using in vivo and in vitro techniques, we show that both promoters are activated by the global transcription factor, CRP (cyclic AMP receptor protein), but the architectures of the EAEC and the UPEC pic promoter differ. Expression from both pic promoters is repressed by the nucleoid-associated factor, Fis, and maximal promoter activity occurs when cells are grown in minimal medium. As CRP activates transcription in conditions of nutrient depletion, whilst Fis levels are maximal in nutrient-rich environments, the regulation of the EAEC and UPEC pic promoters is consistent with Pic’s nutritional role in scavenging mucin as a suitable carbon source during colonisation and infection.

Engineering of global transcription factors in Bacillus, a genetic tool for increasing product yields: a bioprocess overview

Transcriptional factors are well studied in bacteria for their global interactions and the effects they produce at the phenotypic level. Particularly, Bacillus subtilis has been widely employed as a model Gram-positive microorganism used to characterize these network interactions. Bacillus species are currently used as efficient commercial microbial platforms to produce diverse metabolites such as extracellular enzymes, antibiotics, surfactants, industrial chemicals, heterologous proteins, among others. However, the pleiotropic effects caused by the genetic modification of specific genes that codify for global regulators (transcription factors) have not been implicated commonly from a bioprocess point of view. Recently, these strategies have attracted the attention in Bacillus species because they can have an application to increase production efficiency of certain commercial interest metabolites. In this review, we update the recent advances that involve this trend in the use of genetic engineering (mutations, deletion, or overexpression) performed to global regulators such as Spo0A, CcpA, CodY and AbrB, which can provide an advantage for the development or improvement of bioprocesses that involve Bacillus species as production platforms. Genetic networks, regulation pathways and their relationship to the development of growth stages are also discussed to correlate the interactions that occur between these regulators, which are important to consider for application in the improvement of commercial-interest metabolites. Reported yields from these products currently produced mostly under laboratory conditions and, in a lesser extent at bioreactor level, are also discussed to give valuable perspectives about their potential use and developmental level directed to process optimization at large-scale.

Basics on network theory to analyze biological systems: a hands-on outlook

Biological processes result from interactions among molecules and cell-to-cell communications. In the last 50 years, network theory has empowered advances in understanding molecular networks' structure and dynamics that regulate biological systems. Adopting a network data analysis point of view at more laboratories might enrich their research capacity to generate forward working hypotheses. This work briefly describes network theory origins and provides basic graph analysis principles in biological systems, specific centrality measurements, and the main models for network structures. Also, we describe a workflow employing user-friendly free platforms to process, construct, and analyze transcriptome data from a network perspective. With this assay, we expect to encourage the implementation of network theory analysis on biological data in everyday laboratory research.

The integration host factor regulates multiple virulence pathways in bacterial pathogen Dickeya zeae MS2

Dickeya zeae is an aggressive bacterial phytopathogen that infects a wide range of host plants. It has been reported that integration host factor (IHF), a nucleoid-associated protein consisting of IHFα and IHFβ subunits, regulates gene expression by influencing nucleoid structure and DNA bending. To define the role of IHF in the pathogenesis of D. zeae MS2, we deleted either and both of the IHF subunit encoding genes ihfA and ihfB, which significantly reduced the production of cell wall-degrading enzymes (CWDEs), an unknown novel phytotoxin and the virulence factor-modulating (VFM) quorum-sensing (QS) signal, cell motility, biofilm formation, and thereafter the infection ability towards both potato slices and banana seedlings. To characterize the regulatory pathways of IHF protein associated with virulence, IHF binding sites (consensus sequence 5'-WATCAANNNNTTR-3') were predicted and 272 binding sites were found throughout the genome. The expression of 110 tested genes was affected by IHF. Electrophoretic mobility shift assay (EMSA) showed direct interaction of IhfA protein with the promoters of vfmE, speA, pipR, fis, slyA, prtD, hrpL, hecB, hcp, indA, hdaA, flhD, pilT, gcpJ, arcA, arcB, and lysR. This study clarified the contribution of IHF in the pathogenic process of D. zeae by controlling the production of VFM and putrescine QS signals, phytotoxin, and indigoidine, the luxR-solo system, Fis, SlyA, and FlhD transcriptional regulators, and secretion systems from type I to type VI. Characterization of the regulatory networks of IHF in D. zeae provides a target for prevention and control of plant soft rot disease.

The regulon of Brucella abortus two-component system BvrR/BvrS reveals the coordination of metabolic pathways required for intracellular life

Brucella abortus is a facultative intracellular pathogen causing a severe zoonotic disease worldwide. The two-component regulatory system (TCS) BvrR/BvrS of B. abortus is conserved in members of the Alphaproteobacteria class. It is related to the expression of genes required for host interaction and intracellular survival. Here we report that bvrR and bvrS are part of an operon composed of 16 genes encoding functions related to nitrogen metabolism, DNA repair and recombination, cell cycle arrest, and stress response. Synteny of this genomic region within close Alphaproteobacteria members suggests a conserved role in coordinating the expression of carbon and nitrogen metabolic pathways. In addition, we performed a ChIP-Seq analysis after exposure of bacteria to conditions that mimic the intracellular environment. Genes encoding enzymes at metabolic crossroads of the pentose phosphate shunt, gluconeogenesis, cell envelope homeostasis, nucleotide synthesis, cell division, and virulence are BvrR/BvrS direct targets. A 14 bp DNA BvrR binding motif was found and investigated in selected gene targets such as virB1, bvrR, pckA, omp25, and tamA. Understanding gene expression regulation is essential to elucidate how Brucella orchestrates a physiological response leading to a furtive pathogenic strategy.

Roles of Two-Component Signal Transduction Systems in Shigella Virulence

Two-component signal transduction systems (TCSs) are widespread types of protein machinery, typically consisting of a histidine kinase membrane sensor and a cytoplasmic transcriptional regulator that can sense and respond to environmental signals. TCSs are responsible for modulating genes involved in a multitude of bacterial functions, including cell division, motility, differentiation, biofilm formation, antibiotic resistance, and virulence. Pathogenic bacteria exploit the capabilities of TCSs to reprogram gene expression according to the different niches they encounter during host infection. This review focuses on the role of TCSs in regulating the virulence phenotype of Shigella, an intracellular pathogen responsible for severe human enteric syndrome. The pathogenicity of Shigella is the result of the complex action of a wide number of virulence determinants located on the chromosome and on a large virulence plasmid. In particular, we will discuss how five TCSs, EnvZ/OmpR, CpxA/CpxR, ArcB/ArcA, PhoQ/PhoP, and EvgS/EvgA, contribute to linking environmental stimuli to the expression of genes related to virulence and fitness within the host. Considering the relevance of TCSs in the expression of virulence in pathogenic bacteria, the identification of drugs that inhibit TCS function may represent a promising approach to combat bacterial infections.

Single mutation makes Escherichia coli an insect mutualist

Microorganisms often live in symbiosis with their hosts, and some are considered mutualists, where all species involved benefit from the interaction. How free-living microorganisms have evolved to become mutualists is unclear. Here we report an experimental system in which non-symbiotic Escherichia coli evolves into an insect mutualist. The stinkbug Plautia stali is typically associated with its essential gut symbiont, Pantoea sp., which colonizes a specialized symbiotic organ. When sterilized newborn nymphs were infected with E. coli rather than Pantoea sp., only a few insects survived, in which E. coli exhibited specific localization to the symbiotic organ and vertical transmission to the offspring. Through transgenerational maintenance with P. stali, several hypermutating E. coli lines independently evolved to support the host's high adult emergence and improved body colour; these were called 'mutualistic' E. coli. These mutants exhibited slower bacterial growth, smaller size, loss of flagellar motility and lack of an extracellular matrix. Transcriptomic and genomic analyses of 'mutualistic' E. coli lines revealed independent mutations that disrupted the carbon catabolite repression global transcriptional regulator system. Each mutation reproduced the mutualistic phenotypes when introduced into wild-type E. coli, confirming that single carbon catabolite repression mutations can make E. coli an insect mutualist. These findings provide an experimental system for future work on host-microbe symbioses and may explain why microbial mutualisms are omnipresent in nature.

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 transcription factor network of E. coli steers global responses to shifts in RNAP concentration

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

New xylose transporters support the simultaneous consumption of glucose and xylose in Escherichia coli

Glucose and xylose are two major components of lignocellulose. Simultaneous consumption of glucose and xylose is critical for engineered microorganisms to produce fuels and chemicals from lignocellulosic biomass. Although many production limitations have been resolved, glucose-induced inhibition of xylose transport remains a challenge. In this study, a cell growth-based screening strategy was designed to identify xylose transporters uninhibited by glucose. The glucose pathway was genetically blocked in Escherichia coli so that glucose functions only as an inhibitor and cells need xylose as the carbon source for survival. Through adaptive evolution, omics analysis and reverse metabolic engineering, a new phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS) galactitol transporter (GalABC, encoded by EcolC_1640, EcolC_1641, and EcolC_1642 genes) that is not inhibited by glucose was identified. Inactivation of adenylate cyclase led to increased expression of the EcolC_1642 gene, and a point mutation in gene EcolC_1642 (N13S) further enhanced xylose transport. During the second round of gene mining, AraE and a new ABC transporter (AraFGH) of xylose were identified. A point mutation in the transcription regulator araC (L156I) caused increased expression of araE and araFGH genes without arabinose induction, and a point mutation in araE (D223Y) further enhanced xylose transport. These newly identified xylose transporters can support the simultaneous consumption of glucose and xylose and have potential use in producing chemicals from lignocellulose.

ERMer: a serverless platform for navigating, analyzing, and visualizing Escherichia coli regulatory landscape through graph database

Cellular regulation is inherently complex, and one particular cellular function is often controlled by a cascade of different types of regulatory interactions. For example, the activity of a transcription factor (TF), which regulates the expression level of downstream genes through transcriptional regulation, can be regulated by small molecules through compound–protein interactions. To identify such complex regulatory cascades, traditional relational databases require ineffective additional operations and are computationally expensive. In contrast, graph databases are purposefully developed to execute such deep searches efficiently. Here, we present ERMer (E. coli Regulation Miner), the first cloud platform for mining the regulatory landscape of Escherichia coli based on graph databases. Combining the AWS Neptune graph database, AWS lambda function, and G6 graph visualization engine enables quick search and visualization of complex regulatory cascades/patterns. Users can also interactively navigate the E. coli regulatory landscape through ERMer. Furthermore, a Q&A module is included to showcase the power of graph databases in answering complex biological questions through simple queries. The backend graph model can be easily extended as new data become available. In addition, the framework implemented in ERMer can be easily migrated to other applications or organisms. ERMer is available at https://ermer.biodesign.ac.cn/.

The ArcAB Two-Component System: Function in Metabolism, Redox Control, and Infection

ArcAB, also known as the Arc system, is a member of the two-component system family of bacterial transcriptional regulators and is composed of sensor kinase ArcB and response regulator ArcA. In this review, we describe the structure and function of these proteins and assess the state of the literature regarding ArcAB as a sensor of oxygen consumption. The bacterial quinone pool is the primary modulator of ArcAB activity, but questions remain for how this regulation occurs. This review highlights the role of quinones and their oxidation state in activating and deactivating ArcB and compares competing models of the regulatory mechanism. The cellular processes linked to ArcAB regulation of central metabolic pathways and potential interactions of the Arc system with other regulatory systems are also reviewed. Recent evidence for the function of ArcAB under aerobic conditions is challenging the long-standing characterization of this system as strictly an anaerobic global regulator, and the support for additional ArcAB functionality in this context is explored. Lastly, ArcAB-controlled cellular processes with relevance to infection are assessed.

[G-quadruplex structures in bacteria: functional properties and prospects for use as biotargets]

G-quadruplexes (G4), non-canonical secondary DNA structures, are intensively investigated for a long time. In eukaryotic organisms they play an important role in the regulation of gene expression and DNA repair. G4 have also been found in the genomes of numerous bacteria and archaea, but their functional role has not yet been fully explored. Nevertheless, their participation in the formation of antigenic variability, pathogenicity, antibiotic resistance and survival in extreme conditions has been established. Currently, many tools have been developed to detect potential G4 sequences and confirm their formation ability. Since the controlled formation and resolution of the quadruplex are significant means for the regulation of genes critical for survival, a promising direction is the search for ligands - compounds that can have a stabilizing effect on the quadruplex structure and thereby alter gene expression. Currently, a number of ligands are already known, their use stops the growth of pathogenic microorganisms. G4 ligands are of interest as potential antibiotics, which are extremely relevant due to the wide spread of drug resistant pathogens.

The Role of Integration Host Factor in Escherichia coli Persister Formation

Persisters represent a small subpopulation of cells that are tolerant of killing by antibiotics and are implicated in the recalcitrance of chronic infections to antibiotic therapy. One general theme has emerged regarding persisters formed by different bacterial species, namely, a state of relative dormancy characterized by diminished activity of antibiotic targets. Within this framework, a number of studies have linked persister formation to stochastic decreases in energy-generating components, leading to low ATP and target activity. In this study, we screen knockouts in the main global regulators of Escherichia coli for their effect on persisters. A knockout in integration host factor (IHF) had elevated ATP and a diminished level of persisters. This was accompanied by an overexpression of isocitrate dehydrogenase (Icd) and a downregulation of isocitrate lyase (AceA), two genes located at the bifurcation between the tricarboxylic acid (TCA) cycle and the glyoxylate bypass. Using a translational ihfA-mVenus fusion, we sort out rare bright cells, and this subpopulation is enriched in persisters. Our results suggest that noise in the expression of ihf produces rare cells with low Icd/high AceA, diverting substrates into the glyoxylate bypass, which decreases ATP, leading to antibiotic-tolerant persisters. We further examine noise in a simple model, the lac operon, and show that a knockout of the lacI repressor increases expression of the operon and decreases persister formation. Our results suggest that noise quenching by overexpression serves as a general approach to determine the nature of persister genes in a variety of bacterial species and conditions. IMPORTANCE Persisters are phenotypic variants that survive exposure to antibiotics through temporary dormancy. Mutants with increased levels of persisters have been identified in clinical isolates, and evidence suggests these cells contribute to chronic infections and antibiotic treatment failure. Understanding the underlying mechanism of persister formation and tolerance is important for developing therapeutic approaches to treat chronic infections. In this study, we examine a global regulator, IHF, that plays a role in persister formation. We find that noise in expression of IHF contributes to persister formation, likely by regulating the switch between the TCA cycle that efficiently produces energy and the glyoxylate bypass. We extend this study to a simple model lac operon and show that when grown on lactose as the sole carbon source, noise in its expression influences ATP levels and determines persister formation. This noise is quenched by overexpression of the lac operon, providing a simple approach to test the involvement of a gene in persister formation.

Adaptive strategies of uropathogenic Escherichia coli CFT073: from growth in lab media to virulence during host cell adhesion

Urinary tract infections (UTIs) are a major concern in public health. The prevalent uropathogenic bacterium in healthcare settings is Escherichia coli. The increasing rate of antibiotic-resistant strains demands studies to understand E. coli pathogenesis to drive the development of new therapeutic approaches. This study compared the gene expression profile of selected target genes in the prototype uropathogenic E. coli (UPEC) strain CFT073 grown in Luria Bertani (LB), artificial urine (AU), and during adhesion to host bladder cells by semi-quantitative real-time PCR (RT-PCR) assays. AU effectively supported the growth of strain CFT073 as well as other E. coli strains with different lifestyles, thereby confirming the appropriateness of this medium for in vitro models. Unexpectedly, gene expression of strain CFT073 in LB and AU was quite similar; conversely, during the adhesion assay, adhesins and porins were upregulated, while key global regulators were downregulated with respect to lab media. Interestingly, fimH and papGII genes were significantly expressed in all tested conditions. Taken together, these results provide for the first time insights of the metabolic and pathogenic profile of strain CFT073 during the essential phase of host cell adhesion.

Global pleiotropic effects in adaptively evolved Escherichia coli lacking CRP reveal molecular mechanisms that define the growth physiology

Evolution facilitates emergence of fitter phenotypes by efficient allocation of cellular resources in conjunction with beneficial mutations. However, system-wide pleiotropic effects that redress the perturbations to the apex node of the transcriptional regulatory networks remain unclear. Here, we elucidate that absence of global transcriptional regulator CRP in Escherichia coli results in alterations in key metabolic pathways under glucose respiratory conditions, favouring stress- or hedging-related functions over growth-enhancing functions. Further, we disentangle the growth-mediated effects from the CRP regulation-specific effects on these metabolic pathways. We quantitatively illustrate that the loss of CRP perturbs proteome efficiency, as evident from metabolic as well as ribosomal proteome fractions, that corroborated with intracellular metabolite profiles. To address how E. coli copes with such systemic defect, we evolved Δcrp mutant in the presence of glucose. Besides acquiring mutations in the promoter of glucose transporter ptsG, the evolved populations recovered the metabolic pathways to their pre-perturbed state coupled with metabolite re-adjustments, which altogether enabled increased growth. By contrast to Δcrp mutant, the evolved strains remodelled their proteome efficiency towards biomass synthesis, albeit at the expense of carbon efficiency. Overall, we comprehensively illustrate the genetic and metabolic basis of pleiotropic effects, fundamental for understanding the growth physiology.

Characterization of a thermophilic facultatively anaerobic bacterium Paenibacillus sp. strain DA-C8 that exhibits xylan degradation under anaerobic conditions

The screening, identification, and study of the functional properties of cellulolytic xylanolytic bacteria are crucial for the construction of applicable bioprocesses. The thermophilic facultatively anaerobic, xylanolytic bacterial strain DA-C8 (=JCM34211=DSM111723) exhibiting efficient xylan degradation was newly isolated from compost. Strain DA-C8 completely degraded 1% beechwood xylan within 4 days under anaerobic conditions. By 16S rRNA gene sequence homology and phylogenetic tree analysis, strain DA-C8 was closely related to Paenibacillus cisolokensis and Xylanibacillus composti; however, the average nucleotide identity and digital DNA-DNA hybridization values based on genome information and the carbon source utilization properties indicated that strain DA-C8 belongs to Paenibacillus rather than Xylanibacillus. The gene numbers of xylanase and endoglucanase of strain DA-C8 and X. composti were not different; however, strain DA-C8 had higher abundance of α-L-arabinofuranosidase, β-xylosidase, and β-glucosidase than X. composti. Strain DA-C8 showed decreased xylan and corn hull degradation abilities and growth on xylan medium under aerobic conditions. Quantitative PCR showed high expression of xylan and cellulose degradation genes under anaerobic conditions, but the genes were repressed under aerobic conditions, indicating that strain DA-C8 controls polysaccharide degradation depending on the aeration conditions. Strain DA-C8 is a new species of Paenibacillus with a unique polysaccharide degradation system.

The Dynamics of Mycoplasma gallisepticum Nucleoid Structure at the Exponential and Stationary Growth Phases

The structure and dynamics of bacterial nucleoids play important roles in regulating gene expression. Bacteria of class Mollicutes and, in particular, mycoplasmas feature extremely reduced genomes. They lack multiple structural proteins of the nucleoid, as well as regulators of gene expression. We studied the organization of Mycoplasma gallisepticum nucleoids in the stationary and exponential growth phases at the structural and protein levels. The growth phase transition results in the structural reorganization of M. gallisepticum nucleoid. In particular, it undergoes condensation and changes in the protein content. The observed changes corroborate with the previously identified global rearrangement of the transcriptional landscape in this bacterium during the growth phase transition. In addition, we identified that the glycolytic enzyme enolase functions as a nucleoid structural protein in this bacterium. It is capable of non-specific DNA binding and can form fibril-like complexes with DNA.

Construction of an alkaline protease overproducer strain based on Bacillus licheniformis 2709 using an integrative approach

Bacillus licheniformis 2709 is a potential cell factory for the production of alkaline protease AprE, which has important value in industrial application but still lacks sufficient production capacity. To address this problem, we investigated the effects of the secretory viscous materials on the synthesis of AprE, which might seriously affect the industrial fermentation. Furthermore, an iterative chromosomal integration strategy at various chromosomal loci was implemented to achieve stable high-level expression of AprE in B. licheniformis 2709. The host was genetically modified by disrupting the native pgs cluster controlling the biosynthesis of viscous poly-glutamic acid identified in the study by GC/MS, generating a mutant with significantly higher biomass and better bioreactor performance. We further enhanced the expression of alkaline protease by integrating two additional aprE expression cassettes into the genome, generating the integration mutant BL ∆UEP-3 with three aprE expression cassettes, whose AprE enzyme activity in shake flasks reached 25,736 ± 997 U/mL, which was 136% higher than that of the original strain, while the aprE transcription level increased 4.05 times. Thus, an AprE high-yielding strain with excellent fermentation traits was engineered, which was more suitable for bulk-production. Finally, the AprE titer was further increased in a 5-L fermenter, reaching 57,763 ± 1039 U/mL. In summary, genetic modification is an enabling technology for enhancing enzyme production by eliminating the unfavorable characteristics of the host and optimizing the expression of aprE through iterative chromosomal integration. We believe that the protocol developed in this study provides a valuable reference for chromosomal overexpression of proteins or bioactive molecules in other Bacillus species.

Post-translational Lysine Ac(et)ylation in Bacteria: A Biochemical, Structural, and Synthetic Biological Perspective

Ac(et)ylation is a post-translational modification present in all domains of life. First identified in mammals in histones to regulate RNA synthesis, today it is known that is regulates fundamental cellular processes also in bacteria: transcription, translation, metabolism, cell motility. Ac(et)ylation can occur at the ε-amino group of lysine side chains or at the α-amino group of a protein. Furthermore small molecules such as polyamines and antibiotics can be acetylated and deacetylated enzymatically at amino groups. While much research focused on N-(ε)-ac(et)ylation of lysine side chains, much less is known about the occurrence, the regulation and the physiological roles on N-(α)-ac(et)ylation of protein amino termini in bacteria. Lysine ac(et)ylation was shown to affect protein function by various mechanisms ranging from quenching of the positive charge, increasing the lysine side chains’ size affecting the protein surface complementarity, increasing the hydrophobicity and by interfering with other post-translational modifications. While N-(ε)-lysine ac(et)ylation was shown to be reversible, dynamically regulated by lysine acetyltransferases and lysine deacetylases, for N-(α)-ac(et)ylation only N-terminal acetyltransferases were identified and so far no deacetylases were discovered neither in bacteria nor in mammals. To this end, N-terminal ac(et)ylation is regarded as being irreversible. Besides enzymatic ac(et)ylation, recent data showed that ac(et)ylation of lysine side chains and of the proteins N-termini can also occur non-enzymatically by the high-energy molecules acetyl-coenzyme A and acetyl-phosphate. Acetyl-phosphate is supposed to be the key molecule that drives non-enzymatic ac(et)ylation in bacteria. Non-enzymatic ac(et)ylation can occur site-specifically with both, the protein primary sequence and the three dimensional structure affecting its efficiency. Ac(et)ylation is tightly controlled by the cellular metabolic state as acetyltransferases use ac(et)yl-CoA as donor molecule for the ac(et)ylation and sirtuin deacetylases use NAD⁺ as co-substrate for the deac(et)ylation. Moreover, the accumulation of ac(et)yl-CoA and acetyl-phosphate is dependent on the cellular metabolic state. This constitutes a feedback control mechanism as activities of many metabolic enzymes were shown to be regulated by lysine ac(et)ylation. Our knowledge on lysine ac(et)ylation significantly increased in the last decade predominantly due to the huge methodological advances that were made in fields such as mass-spectrometry, structural biology and synthetic biology. This also includes the identification of additional acylations occurring on lysine side chains with supposedly different regulatory potential. This review highlights recent advances in the research field. Our knowledge on enzymatic regulation of lysine ac(et)ylation will be summarized with a special focus on structural and mechanistic characterization of the enzymes, the mechanisms underlying non-enzymatic/chemical ac(et)ylation are explained, recent technological progress in the field are presented and selected examples highlighting the important physiological roles of lysine ac(et)ylation are summarized.

Restoring Global Gene Regulation through Experimental Evolution Uncovers a NAP (Nucleoid-Associated Protein)-Like Behavior of Crp/Cap

How do hierarchical gene regulation networks evolve in bacteria? Nucleoid-associated proteins (NAPs) influence the overall structure of bacterial genomes, sigma factors and global transcription factors (TFs) control thousands of genes, and many operons are regulated by highly specific TFs that in turn are controlled allosterically by cellular metabolites. These regulatory hierarchies have been shaped by millions of years of evolution to optimize fitness in response to changing environmental conditions, but it is unclear how NAPs and TFs relate and have evolved together. Cyclic AMP (cAMP) receptor protein (Crp) is the paradigmatic global TF in Escherichia coli, and here we report that mutations in the topA gene compensate for loss of cAMP, showing that the interplay between Crp and the supercoiling status of promoters is key to global stress response. Furthermore, we observed an effect of apoCrp on gene expression in the absence of its effector cAMP. This provides support for the proposed NAP-like role for Crp, suggesting that it represents an intermediate point in the evolution of a ligand-controlled TF from a NAP.

Blocking, Bending, and Binding: Regulation of Initiation of Chromosome Replication During the Escherichia coli Cell Cycle by Transcriptional Modulators That Interact With Origin DNA

Genome duplication is a critical event in the reproduction cycle of every cell. Because all daughter cells must inherit a complete genome, chromosome replication is tightly regulated, with multiple mechanisms focused on controlling when chromosome replication begins during the cell cycle. In bacteria, chromosome duplication starts when nucleoprotein complexes, termed orisomes, unwind replication origin (oriC) DNA and recruit proteins needed to build new replication forks. Functional orisomes comprise the conserved initiator protein, DnaA, bound to a set of high and low affinity recognition sites in oriC. Orisomes must be assembled each cell cycle. In Escherichia coli, the organism in which orisome assembly has been most thoroughly examined, the process starts with DnaA binding to high affinity sites after chromosome duplication is initiated, and orisome assembly is completed immediately before the next initiation event, when DnaA interacts with oriC’s lower affinity sites, coincident with origin unwinding. A host of regulators, including several transcriptional modulators, targets low affinity DnaA-oriC interactions, exerting their effects by DNA bending, blocking access to recognition sites, and/or facilitating binding of DnaA to both DNA and itself. In this review, we focus on orisome assembly in E. coli. We identify three known transcriptional modulators, SeqA, Fis (factor for inversion stimulation), and IHF (integration host factor), that are not essential for initiation, but which interact directly with E. coli oriC to regulate orisome assembly and replication initiation timing. These regulators function by blocking sites (SeqA) and bending oriC DNA (Fis and IHF) to inhibit or facilitate cooperative low affinity DnaA binding. We also examine how the growth rate regulation of Fis levels might modulate IHF and DnaA binding to oriC under a variety of nutritional conditions. Combined, the regulatory mechanisms mediated by transcriptional modulators help ensure that at all growth rates, bacterial chromosome replication begins once, and only once, per cell cycle.

Hierarchy of transcription factor network in Escherichia coli K-12: H-NS-mediated silencing and Anti-silencing by global regulators

Transcriptional regulation for genome expression determines growth and adaptation of single-cell bacteria that are directly exposed to environment. The transcriptional apparatus in Escherichia coli K-12 is composed of RNA polymerase core enzyme and two groups of its regulatory proteins, seven species of promoter-recognition subunit sigma and about 300 species of transcription factors. The identification of regulatory targets for all these regulatory proteins is critical toward understanding the genome regulation as a whole. For this purpose, we performed a systematic search in vitro of the whole set of binding sites for each factor by gSELEX system. This review summarizes the accumulated knowledge of regulatory targets for more than 150 TFs from E. coli K-12. Overall TFs could be classified into four families: nucleoid-associated bifunctional TFs; global regulators; local regulators; and single-target regulators, in which the regulatory functions remain uncharacterized for the nucleoid-associated TFs. Here we overview the regulatory targets of two nucleoid-associated TFs, H-NS and its paralog StpA, both together playing the silencing role of a set of non-essential genes. Participation of LeuO and other global regulators have been indicated for the anti-silencing. Finally, we propose the hierarchy of TF network as a key framework of the bacterial genome regulation.

Single-Target Regulators Constitute the Minority Group of Transcription Factors in Escherichia coli K-12

The identification of regulatory targets of all transcription factors (TFs) is critical for understanding the entire network of genome regulation. A total of approximately 300 TFs exist in the model prokaryote Escherichia coli K-12, but the identification of whole sets of their direct targets is impossible with use of in vivo approaches. For this end, the most direct and quick approach is to identify the TF-binding sites in vitro on the genome. We then developed and utilized the gSELEX screening system in vitro for identification of more than 150 E. coli TF-binding sites along the E. coli genome. Based on the number of predicted regulatory targets, we classified E. coli K-12 TFs into four groups, altogether forming a hierarchy ranging from a single-target TF (ST-TF) to local TFs, global TFs, and nucleoid-associated TFs controlling as many as 1,000 targets. Using the collection of purified TFs and a library of genome DNA segments from a single and the same E. coli K-12, we identified here a total of 11 novel ST-TFs, CsqR, CusR, HprR, NorR, PepA, PutA, QseA, RspR, UvrY, ZraR, and YqhC. The regulation of single-target promoters was analyzed in details for the hitherto uncharacterized QseA and RspR. In most cases, the ST-TF gene and its regulatory target genes are adjacently located on the E. coli K-12 genome, implying their simultaneous transfer in the course of genome evolution. The newly identified 11 ST-TFs and the total of 13 hitherto identified altogether constitute the minority group of TFs in E. coli K-12.

Enhanced production of nonanedioic acid from nonanoic acid by engineered Escherichia coli

In this study, whole-cell biotransformation was conducted to produce nonanedioic acid from nonanoic acid by expressing the alkane hydroxylating system (AlkBGT) from Pseudomonas putida GPo1 in Escherichia coli. Following adaptive laboratory evolution, an efficient E. coli mutant strain, designated as MRE, was successfully obtained, demonstrating the fastest growth (27-fold higher) on nonanoic acid as the sole carbon source compared to the wild-type strain. Additionally, the MRE strain was engineered to block nonanoic acid degradation by deleting fadE. The resulting strain exhibited a 12.8-fold increase in nonanedioic acid production compared to the wild-type strain. Six mutations in acrR, P_crp , dppA, P_fadD , e14, and yeaR were identified in the mutant MRE strain, which was characterized using genomic modifications and RNA-sequencing. The acquired mutations were found to be beneficial for rapid growth and nonanedioic acid production.

Construction of recombinant Escherichia coli for production of L-phenylalanine-derived compounds

L-phenylalanine is an important amino acid that is widely used in the fields of food flavors and pharmaceuticals. Apart from L-phenylalanine itself, various commercially valuable chemical compounds can also be generated via the L-phenylalanine biosynthesis pathway. Compared with direct extraction from plants or synthesis by chemical reaction, microbial production of L-phenylalanine -derived compounds can overcome the drawbacks of environmental pollution, low yield, and mixtures of stereoisomeric products. Accordingly, increasing intracellular levels of precursors, deregulating feedback inhibition and transcription repression, engineering global regulators and other effective strategies have been implemented to produce different L-phenylalanine -derived compounds in the excellent chassis host Escherichia coli. Finally, this review highlights principal strategies for improving the production of L-phenylalanine and/or its derivatives in E. coli, and discusses the future outlook for further enhancing the titer and yields of these compounds.

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.

Piggybacking on Niche Adaptation Improves the Maintenance of Multidrug-Resistance Plasmids

The persistence of plasmids in bacterial populations represents a puzzling evolutionary problem with serious clinical implications due to their role in the ongoing antibiotic resistance crisis. Recently, major advancements have been made toward resolving this “plasmid paradox” but mainly in a nonclinical context. Here, we propose an additional explanation for the maintenance of multidrug‐resistance plasmids in clinical Escherichia coli strains. After coevolving two multidrug‐resistance plasmids encoding resistance to last resort carbapenems with an extraintestinal pathogenic E. coli strain, we observed that chromosomal media adaptive mutations in the global regulatory systems CCR (carbon catabolite repression) and ArcAB (aerobic respiration control) pleiotropically improved the maintenance of both plasmids. Mechanistically, a net downregulation of plasmid gene expression reduced the fitness cost. Our results suggest that global chromosomal transcriptional rewiring during bacterial niche adaptation may facilitate plasmid maintenance.

Elucidation of the Pathogenicity-Associated Regulatory Network in Xanthomonas oryzae pv. oryzae

Xanthomonas is a notorious plant pathogen causing serious diseases in hundreds of plant hosts. Xanthomonas species are equipped with an array of signal transduction systems that regulate gene expression to survive in various harsh environments and successfully infect hosts. Although certain pathogenicity-associated regulators have been functionally characterized, signal transduction systems always function as a regulatory network which remains to be elucidated in Xanthomonas. This study used a systematic approach to characterize all identified pathogenicity-associated regulators in Xanthomonas oryzae pv. oryzae (Xoo), including a transcriptional regulator with unknown function, and their interactive regulatory network. RNA sequencing was used in elucidating the patterns of the 10 pathogenicity-associated regulators identified. Results revealed that each pathogenicity-associated regulator has cross talk with others and all these regulators function as a regulatory network, with VemR and PXO_RS20790 being the master pathogenicity-associated regulators and HrpX being the final executant. Moreover, regulome analysis showed that numerous genes other than genes in pathogenicity islands are finely regulated within the regulatory network. Given that most of the pathogenicity-associated regulators are conserved in Xanthomonadales, our findings suggest a global network of gene regulation in this evolutionarily conserved pathogen. In conclusion, our study provides essential basic information about the regulatory network in Xoo, suggesting that this complicated regulatory network is one of the reasons for the robustness and fitness of Xanthomonas spp.

IMPORTANCE The host plant infection process of pathogenic bacteria is a coordinating cellular behavior, which requires dynamic regulation at several levels in response to variations in host plants or fluctuations in the external environment. As one of the most important genera of plant-pathogenic bacteria, Xanthomonas has been studied as a model. Although certain pathogenicity-associated regulators have been functionally characterized, interactions among them remain to be elucidated. This study systematically characterized pathogenicity-associated regulators in Xoo and revealed that cross talk exists among pathogenicity-associated regulators and function as a regulatory network in which a hierarchy exists among the regulators. Our study elucidated the landscape of the pathogenicity-associated regulatory network in Xanthomonas, promoting understanding of the infection process of pathogenic bacteria.

Devil in the details: Mechanistic variations impact information transfer across models of transcriptional cascades

The transcriptional network determines a cell’s internal state by regulating protein expression in response to changes in the local environment. Due to the interconnected nature of this network, information encoded in the abundance of various proteins will often propagate across chains of noisy intermediate signaling events. The data-processing inequality (DPI) leads us to expect that this intracellular game of “telephone” should degrade this type of signal, with longer chains losing successively more information to noise. However, a previous modeling effort predicted that because the steps of these signaling cascades do not truly represent independent stages of data processing, the limits of the DPI could seemingly be surpassed, and the amount of transmitted information could actually increase with chain length. What that work did not examine was whether this regime of growing information transmission was attainable by a signaling system constrained by the mechanistic details of more complex protein-binding kinetics. Here we address this knowledge gap through the lens of information theory by examining a model that explicitly accounts for the binding of each transcription factor to DNA. We analyze this model by comparing stochastic simulations of the fully nonlinear kinetics to simulations constrained by the linear response approximations that displayed a regime of growing information. Our simulations show that even when molecular binding is considered, there remains a regime wherein the transmitted information can grow with cascade length, but ends after a critical number of links determined by the kinetic parameter values. This inflection point marks where correlations decay in response to an oversaturation of binding sites, screening informative transcription factor fluctuations from further propagation down the chain where they eventually become indistinguishable from the surrounding levels of noise.

Rewiring the specificity of extracytoplasmic function sigma factors

Bacterial genomes are being sequenced at an exponentially increasing rate, but our inability to decipher their transcriptional wiring limits our ability to derive new biology from these sequences. De novo determination of regulatory interactions requires accurate prediction of regulators' DNA binding and precise determination of biologically significant binding sites. Here we address these challenges by solving the DNA-specificity code of extracytoplasmic function sigma factors (ECF σs), a major family of bacterial regulators, and determining their putative regulons. We generated an aligned collection of ECF σs and their promoters by leveraging the autoregulatory nature of ECF σs as a means of promoter discovery and analyzed it to identify and characterize the conserved amino acid-nucleotide interactions that determine promoter specificity. This enabled de novo prediction of ECF σ specificity, which we combined with a statistically rigorous phylogenetic footprinting pipeline based on precomputed orthologs to predict the direct targets of ∼67% of ECF σs. This global survey indicated that some ECF σs are conserved global regulators controlling many genes throughout the genome, which are important under many conditions, while others are local regulators, controlling a few closely linked genes in response to specific stimuli in select species. This analysis reveals important organizing principles of bacterial gene regulation and presents a conceptual and computational framework for deciphering gene regulatory networks.

RegulomePA: a database of transcriptional regulatory interactions in Pseudomonas aeruginosa PAO1

We present RegulomePA, a database that contains biological information on regulatory interactions between transcription factors (TFs), sigma factor (SFs) and target genes in Pseudomonas aeruginosa PAO1. RegulomePA consists of 4827 regulatory interactions between 2831 nodes, which represent the interactions of TFs and SFs with their target genes, from the total of predicted RegulomePA including 27.27% of the TFs, 54.16% of SFs and 50.8% of the total genes. Each entry in the database corresponds to one node in the network and provides comprehensive details about the gene and its regulatory interactions such as gene description, nucleotide sequence, genome-strand position and links to other databases as well as the type of regulation it exerts or to which it is being subject (repression or activation), the associated experimental evidence and references, and topological information. Additionally, RegulomePA provides a way to recover information on the regulatory circuits of the network to which a gene pertains and also makes available the source codes to analyze the topology of any other regulatory network. The database will be updated yearly, by our team, with the contributions from ourselves and users, since the users are provided with an interactive platform where they can add interactions to the regulatory network feeding it with their respective references. Database URL: www.regulome.pcyt.unam.mx.

Contextual centrality: going beyond network structure

Centrality is a fundamental network property that ranks nodes by their structural importance. However, the network structure alone may not predict successful diffusion in many applications, such as viral marketing and political campaigns. We propose contextual centrality, which integrates structural positions, the diffusion process, and, most importantly, relevant node characteristics. It nicely generalizes and relates to standard centrality measures. We test the effectiveness of contextual centrality in predicting the eventual outcomes in the adoption of microfinance and weather insurance. Our empirical analysis shows that the contextual centrality of first-informed individuals has higher predictive power than that of other standard centrality measures. Further simulations show that when the diffusion occurs locally, contextual centrality can identify nodes whose local neighborhoods contribute positively. When the diffusion occurs globally, contextual centrality signals whether diffusion may generate negative consequences. Contextual centrality captures more complicated dynamics on networks than traditional centrality measures and has significant implications for network-based interventions.