The whole set of constitutive promoters recognized by RNA polymerase RpoD holoenzyme of Escherichia coli

AI and Knowledge-Based Method for Rational Design of Escherichia coli Sigma70 Promoters

Expanding sigma70 promoter libraries can support the engineering of metabolic pathways and enhance recombinant protein expression. Herein, we developed an artificial intelligence (AI) and knowledge-based method for the rational design of sigma70 promoters. Strong sigma70 promoters were identified by using high-throughput screening (HTS) with enhanced green fluorescent protein (eGFP) as a reporter gene. The features of these strong promoters were adopted to guide promoter design based on our previous reported deep learning model. In the following case study, the obtained strong promoters were used to express collagen and microbial transglutaminase (mTG), resulting in increased expression levels by 81.4% and 33.4%, respectively. Moreover, these constitutive promoters achieved soluble expression of mTG-activating protease and contributed to active mTG expression in Escherichia coli. The results suggested that the combined method may be effective for promoter engineering.

Alteration of the -35 and -10 sequences and deletion the upstream sequence of the -35 region of the promoter A1 of the phage T7 in dsDNA confirm the contribution of non-specific interactions with E. coli RNA polymerase to the transcription initiation process

Transcription initiation is a multi-step process, in which the RNA polymerase holoenzyme binds to the specific promoter sequences to form a closed complex, which, through intermediate stages, isomerizes into an open complex capable of initiating the productive phase of transcription. The aim of this work was to determine the contribution of the -10 and -35 regions of the promoter, as well as the role of non-specific interactions, in the binding of RNA polymerase and the formation of an active initiation complex capable of transcription. Therefore, fragments of promoter DNA, derived from the strong promoter A1 of the phage T7, containing completely and partially altered elements -35 and -10, and devoid of an upstream region, were constructed using genetic engineering methods. Functional analyses of modified promoter fragments were carried out, checking their ability to form binary complexes with Escherichia coli RNA polymerase (RNAP) and the efficiency of converting binary complexes into triple complexes characteristic of the productive phase of transcription. The obtained results suggest that, in relation to the A1 promoter of the T7 phage, the most important role of the -35 region is carrying the open complex through the next phases of transcription initiation. The weakening of specific impacts within the region -35 is a reason for the defect associated with the transformation of the open complex, formed by a DNA fragment containing the completely altered -35 region, into elongation and the impairment of RNA synthesis. This leads to breaking contacts with the RNA polymerase holoenzyme, and destabilization and disintegration of the complex in the initial phase of productive transcription. This confirms the hypothesis of the so-called stressed intermediate state associated with the stage of transition from the open complex to the elongation complex. The experiments carried out in this work confirm also that the process of promoter localization and recognition, as well as the formation of binary complexes, is sequential in nature, and that the region located upstream of the -35 hexamer, and the hexamer itself, plays here an additive role.

Genomic SELEX Screening of Regulatory Targets of Transcription Factors

The genome of Escherichia coli K-12 is transcribed by a single species of RNA polymerase. The selectivity of transcriptional targets is determined via interaction with one of seven species of the sigma subunit and a total of approximately 300 species of transcription factor (TFs). For comprehensive identification of the regulatory targets of these two groups of regulatory proteins on the genome, we developed an in vitro approach, "Genomic SELEX" (gSELEX) screening. Here we describe a detailed protocol of the gSELEX screening system, which uses purified regulatory proteins and fragments of genomic DNA from E. coli. Moreover, we describe methods and examples of results using cell-free synthetic proteins.

Deep Learning-Assisted Design of Novel Promoters in Escherichia coli

Deep learning (DL) approaches have the ability to accurately recognize promoter regions and predict their strength. Here, the potential for controllably designing active Escherichia coli promoter is explored by combining multiple deep learning models. First, "DRSAdesign," which relies on a diffusion model to generate different types of novel promoters is created, followed by predicting whether they are real or fake and strength. Experimental validation showed that 45 out of 50 generated promoters are active with high diversity, but most promoters have relatively low activity. Next, "Ndesign," which relies on generating random sequences carrying functional -35 and -10 motifs of the sigma70 promoter is introduced, and their strength is predicted using the designed DL model. The DL model is trained and validated using 200 and 50 generated promoters, and displays Pearson correlation coefficients of 0.49 and 0.43, respectively. Taking advantage of the DL models developed in this work, possible 6-mers are predicted as key functional motifs of the sigma70 promoter, suggesting that promoter recognition and strength prediction mainly rely on the accommodation of functional motifs. This work provides DL tools to design promoters and assess their functions, paving the way for DL-assisted metabolic engineering.

Sense and antisense RNA products of the uxuR gene can affect motility and chemotaxis acting independent of the UxuR protein

Small non-coding and antisense RNAs are widespread in all kingdoms of life, however, the diversity of their functions in bacteria is largely unknown. Here, we study RNAs synthesised from divergent promoters located in the 3'-end of the uxuR gene, encoding transcription factor regulating hexuronate metabolism in Escherichia coli. These overlapping promoters were predicted in silico with rather high scores, effectively bound RNA polymerase in vitro and in vivo and were capable of initiating transcription in sense and antisense directions. The genome-wide correlation between in silico promoter scores and RNA polymerase binding in vitro and in vivo was higher for promoters located on the antisense strands of the genes, however, sense promoters within the uxuR gene were more active. Both regulatory RNAs synthesised from the divergent promoters inhibited expression of genes associated with the E. coli motility and chemotaxis independent of a carbon source on which bacteria had been grown. Direct effects of these RNAs were confirmed for the fliA gene encoding σ²⁸ subunit of RNA polymerase. In addition to intracellular sRNAs, promoters located within the uxuR gene could initiate synthesis of transcripts found in the fraction of RNAs secreted in the extracellular medium. Their profile was also carbon-independent suggesting that intragenic uxuR transcripts have a specific regulatory role not directly related to the function of the protein in which gene they are encoded.

Sigma70Pred: A highly accurate method for predicting sigma70 promoter in Escherichia coli K-12 strains

Sigma70 factor plays a crucial role in prokaryotes and regulates the transcription of most of the housekeeping genes. One of the major challenges is to predict the sigma70 promoter or sigma70 factor binding site with high precision. In this study, we trained and evaluate our models on a dataset consists of 741 sigma70 promoters and 1,400 non-promoters. We have generated a wide range of features around 8,000, which includes Dinucleotide Auto-Correlation, Dinucleotide Cross-Correlation, Dinucleotide Auto Cross-Correlation, Moran Auto-Correlation, Normalized Moreau-Broto Auto-Correlation, Parallel Correlation Pseudo Tri-Nucleotide Composition, etc. Our SVM based model achieved maximum accuracy 97.38% with AUROC 0.99 on training dataset, using 200 most relevant features. In order to check the robustness of the model, we have tested our model on the independent dataset made by using RegulonDB10.8, which included 1,134 sigma70 and 638 non-promoters, and able to achieve accuracy of 90.41% with AUROC of 0.95. Our model successfully predicted constitutive promoters with accuracy of 81.46% on an independent dataset. We have developed a method, Sigma70Pred, which is available as webserver and standalone packages at https://webs.iiitd.edu.in/raghava/sigma70pred/. The services are freely accessible.

Identification of small RNAs associated with RNA chaperone Hfq reveals a new stress response regulator in Actinobacillus pleuropneumoniae

The RNA chaperone Hfq promotes the association of small RNAs (sRNAs) with cognate mRNAs, controlling the expression of bacterial phenotype. Actinobacillus pleuropneumoniae hfq mutants strains are attenuated for virulence in pigs, impaired in the ability to form biofilms, and more susceptible to stress, but knowledge of the extent of sRNA involvement is limited. Here, using A. pleuropneumoniae strain MIDG2331 (serovar 8), 14 sRNAs were identified by co-immunoprecipitation with Hfq and the expression of eight, identified as trans-acting sRNAs, were confirmed by Northern blotting. We focused on one of these sRNAs, named Rna01, containing a putative promoter for RpoE (stress regulon) recognition. Knockout mutants of rna01 and a double knockout mutant of rna01 and hfq, both had decreased biofilm formation and hemolytic activity, attenuation for virulence in Galleria mellonella, altered stress susceptibility, and an altered outer membrane protein profile. Rna01 affected extracellular vesicle production, size and toxicity in G. mellonella. qRT-PCR analysis of rna01 and putative cognate mRNA targets indicated that Rna01 is associated with the extracytoplasmic stress response. This work increases our understanding of the multilayered and complex nature of the influence of Hfq-dependent sRNAs on the physiology and virulence of A. pleuropneumoniae.

Riboswitch-mediated regulation of riboflavin biosynthesis genes in prokaryotes

Prokaryotic organisms frequently use riboswitches to quantify intracellular metabolite concentration via high-affinity metabolite receptors. Riboswitches possess a metabolite-sensing system that controls gene regulation in a cis-acting fashion at the initiation of transcriptional/translational level by binding with a specific metabolite and controlling various biochemical pathways. Riboswitch binds with flavin mononucleotide (FMN), a phosphorylated form of riboflavin and controls gene expression involved in riboflavin biosynthesis and transport pathway. The first step of the riboflavin biosynthesis pathway is initiated by the conversion of guanine nucleotide triphosphate (GTP), which is an intermediate of the purine biosynthesis pathway. An alternative pentose phosphate pathway of riboflavin biosynthesis includes the enzymatic conversion of ribulose-5-phosphate into 3, 4 dihydroxy-2-butanone-4-phosphates by DHBP synthase. The product of ribAB interferes with both GTP cyclohydrolase II as well as DHBP synthase activities, which catalyze the cleavage of GTP and converts DHBP Ribu5P in the initial steps of both riboflavin biosynthesis branches. Riboswitches are located in the 5' untranslated region (5' UTR) of messenger RNAs and contain an aptamer domain (highly conserved in sequence) where metabolite binding leads to a conformational change in an aptamer domain, which modulate the regulation of gene expression located on bacterial mRNA. In this review, we focus on how riboswitch regulates the riboflavin biosynthesis pathway in Bacillus subtilis and Lactobacillus plantarum.

The archetypal gene transfer agent RcGTA is regulated via direct interaction with the enigmatic RNA polymerase omega subunit

Gene transfer agents (GTAs) are small virus-like particles that indiscriminately package and transfer any DNA present in their host cell, with clear implications for bacterial evolution. The first transcriptional regulator that directly controls GTA expression, GafA, was recently discovered, but its mechanism of action has remained elusive. Here, we demonstrate that GafA controls GTA gene expression via direct interaction with the RNA polymerase omega subunit (Rpo-ω) and also positively autoregulates its own expression by an Rpo-ω-independent mechanism. We show that GafA is a modular protein with distinct DNA and protein binding domains. The functional domains we observe in Rhodobacter GafA also correspond to two-gene operons in Hyphomicrobiales pathogens. These data allow us to produce the most complete regulatory model for a GTA and point toward an atypical mechanism for RNA polymerase recruitment and specific transcriptional activation in the Alphaproteobacteria.

How Do Transposable Elements Activate Expression of Transcriptionally Silent Antibiotic Resistance Genes?

The rapidly emerging phenomenon of antibiotic resistance threatens to substantially reduce the efficacy of available antibacterial therapies. Dissemination of resistance, even between phylogenetically distant bacterial species, is mediated mainly by mobile genetic elements, considered to be natural vectors of horizontal gene transfer. Transposable elements (TEs) play a major role in this process-due to their highly recombinogenic nature they can mobilize adjacent genes and can introduce them into the pool of mobile DNA. Studies investigating this phenomenon usually focus on the genetic load of transposons and the molecular basis of their mobility. However, genes introduced into evolutionarily distant hosts are not necessarily expressed. As a result, bacterial genomes contain a reservoir of transcriptionally silent genetic information that can be activated by various transposon-related recombination events. The TEs themselves along with processes associated with their transposition can introduce promoters into random genomic locations. Thus, similarly to integrons, they have the potential to convert dormant genes into fully functional antibiotic resistance determinants. In this review, we describe the genetic basis of such events and by extension the mechanisms promoting the emergence of new drug-resistant bacterial strains.

Dual role of CsrA in regulating the hemolytic activity of Escherichia coli O157:H7

Post-transcriptional global carbon storage regulator A (CsrA) is a sequence-specific RNA-binding protein involved in the regulation of multiple bacterial processes. Hemolysin is an important virulence factor in the enterohemorrhagic Escherichia coli O157:H7 (EHEC). Here, we show that CsrA plays a dual role in the regulation of hemolysis in EHEC. CsrA significantly represses plasmid-borne enterohemolysin (EhxA)-mediated hemolysis and activates chromosome-borne hemolysin E (HlyE)-mediated hemolysis through different mechanisms. RNA structure prediction revealed a well-matched stem-loop structure with two potential CsrA binding sites located on the 5' untranslated region (UTR) of ehxB, which encodes a translocator required for EhxA secretion. CsrA inhibits EhxA secretion by directly binding to the RNA leader sequence of ehxB to repress its expression in two different ways: CsrA either binds to the Shine–Dalgarno sequence of ehxB to block ribosome access or to ehxB transcript to promote its mRNA decay. The predicted CsrA-binding site 1 of ehxB is essential for its regulation. There is a single potential CsrA-binding site at the 5'-end of the hlyE transcript, and its mutation completely abolishes CsrA-dependent activation. CsrA can also stabilize hlyE mRNA by directly binding to its 5' UTR. Overall, our results indicate that CsrA acts as a hemolysis modulator to regulate pathogenicity under certain conditions.

CRISPR/dCas9-RpoD-Mediated Simultaneous Transcriptional Activation and Repression in Shewanella oneidensis MR-1

Extracellular electron transfer (EET) of electroactive microorganisms (EAMs) is the dominating factor for versatile applications of bio-electrochemical systems. Shewanella oneidensis MR-1 is one of the model EAMs for the study of EET, which is associated with a variety of cellular activities. However, due to the lack of a transcriptional activation tool, regulation of multiple genes is labor-intensive and time-consuming, which hampers the advancement of improving the EET efficiency in S. oneidensis. In this study, we developed an easily operated and multifunctional regulatory tool, that is, a simultaneous clustered regularly interspaced short palindromic repeats (CRISPR)-mediated transcriptional activation (CRISPRa) and interference (CRISPRi) system, for application in S. oneidensis. First, a large number of activators were screened, and RpoD (σ⁷⁰) was determined as the optimal activator. Second, the effective activation range was identified to be 190-216 base upstream of the transcriptional start site. Third, up- and downregulation was achieved in concert by two orthogonal single guide RNAs targeting different positions. The activation of the cell division gene (minCDE) and repression of the cytotoxic gene (SO_3166) were concurrently implemented, increasing the power density by 2.5-fold and enhancing the degradation rate of azo dyes by 2.9-fold. The simultaneous CRISPRa and CRISPRi system enables simultaneous multiplex genetic regulation, offering the potential to further advance studies of the EET mechanism and application in S. oneidensis.

Importance of RpoD- and Non-RpoD-Dependent Expression of Horizontally Acquired Genes in Cupriavidus metallidurans

The genome of the metal-resistant, hydrogen-oxidizing bacterium Cupriavidus metallidurans contains a large number of horizontally acquired plasmids and genomic islands that were integrated into its chromosome or chromid. For the C. metallidurans CH34 wild-type strain growing under nonchallenging conditions, 5,763 transcriptional starting sequences (TSSs) were determined. Using a custom-built motif discovery software based on hidden Markov models, patterns upstream of the TSSs were identified. The pattern TTGACA, −35.6 ± 1.6 bp upstream of the TSSs, in combination with a TATAAT sequence 15.8 ± 1.4 bp upstream occurred frequently, especially upstream of the TSSs for 48 housekeeping genes, and these were assigned to promoters used by RNA polymerase containing the main housekeeping sigma factor RpoD. From patterns upstream of the housekeeping genes, a score for RpoD-dependent promoters in C. metallidurans was derived and applied to all 5,763 TSSs. Among these, 2,572 TSSs could be associated with RpoD with high probability, 373 with low probability, and 2,818 with no probability. In a detailed analysis of horizontally acquired genes involved in metal resistance and not involved in this process, the TSSs responsible for the expression of these genes under nonchallenging conditions were assigned to RpoD- or non-RpoD-dependent promoters. RpoD-dependent promoters occurred frequently in horizontally acquired metal resistance and other determinants, which should allow their initial expression in a new host. However, other sigma factors and sense/antisense effects also contribute—maybe to mold in subsequent adaptation steps the assimilated gene into the regulatory network of the cell.

IMPORTANCE In their natural environment, bacteria are constantly acquiring genes by horizontal gene transfer. To be of any benefit, these genes should be expressed. We show here that the main housekeeping sigma factor RpoD plays an important role in the expression of horizontally acquired genes in the metal-resistant hydrogen-oxidizing bacterium C. metallidurans. By conservation of the RpoD recognition consensus sequence, a newly arriving gene has a high probability to be expressed in the new host cell. In addition to integrons and genes travelling together with that for their sigma factor, conservation of the RpoD consensus sequence may be an important contributor to the overall evolutionary success of horizontal gene transfer in bacteria. Using C. metallidurans as an example, this publication sheds some light on the fate and function of horizontally acquired genes in bacteria.

Whole set of constitutive promoters for RpoN sigma factor and the regulatory role of its enhancer protein NtrC in Escherichia coli K-12

The promoter selectivity of Escherichia coli RNA polymerase (RNAP) is determined by its promoter-recognition sigma subunit. The model prokaryote E. coli K-12 contains seven species of the sigma subunit, each recognizing a specific set of promoters. Using genomic SELEX (gSELEX) screening in vitro, we identified the whole set of ‘constitutive’ promoters recognized by the reconstituted RNAP holoenzyme alone, containing RpoD (σ⁷⁰), RpoS (σ³⁸), RpoH (σ³²), RpoF (σ²⁸) or RpoE (σ²⁴), in the absence of other supporting regulatory factors. In contrast, RpoN sigma (σ⁵⁴), involved in expression of nitrogen-related genes and also other cellular functions, requires an enhancer (or activator) protein, such as NtrC, for transcription initiation. In this study, a series of gSELEX screenings were performed to search for promoters recognized by the RpoN RNAP holoenzyme in the presence and absence of the major nitrogen response enhancer NtrC, the best-characterized enhancer. Based on the RpoN holoenzyme-binding sites, a total of 44 to 61 putative promoters were identified, which were recognized by the RpoN holoenzyme alone. In the presence of the enhancer NtrC, the recognition target increased to 61–81 promoters. Consensus sequences of promoters recognized by RpoN holoenzyme in the absence and presence of NtrC were determined. The promoter activity of a set of NtrC-dependent and -independent RpoN promoters was verified in vivo under nitrogen starvation, in the presence and absence of RpoN and/or NtrC. The promoter activity of some RpoN-recognized promoters increased in the absence of RpoN or NtrC, supporting the concept that the promoter-bound NtrC-enhanced RpoN holoenzyme functions as a repressor against RpoD holoenzyme. Based on our findings, we propose a model in which the RpoN holoenzyme fulfils the dual role of repressor and transcriptase for the same set of genes. We also propose that the promoter recognized by RpoN holoenzyme in the absence of enhancers is the ‘repressive’ promoter. The presence of high-level RpoN sigma in growing E. coli K-12 in rich medium may be related to the repression role of a set of genes needed for the utilization of ammonia as a nitrogen source in poor media. The list of newly identified regulatory targets of RpoN provides insight into E. coli survival under nitrogen-depleted conditions in nature.

Temporal evolution of master regulator Crp identifies pyrimidines as catabolite modulator factors

The evolution of microorganisms often involves changes of unclear relevance, such as transient phenotypes and sequential development of multiple adaptive mutations in hotspot genes. Previously, we showed that ageing colonies of an E. coli mutant unable to produce cAMP when grown on maltose, accumulated mutations in the crp gene (encoding a global transcription factor) and in genes involved in pyrimidine metabolism such as cmk; combined mutations in both crp and cmk enabled fermentation of maltose (which usually requires cAMP-mediated Crp activation for catabolic pathway expression). Here, we study the sequential generation of hotspot mutations in those genes, and uncover a regulatory role of pyrimidine nucleosides in carbon catabolism. Cytidine binds to the cytidine regulator CytR, modifies the expression of sigma factor 32 (RpoH), and thereby impacts global gene expression. In addition, cytidine binds and activates a Crp mutant directly, thus modulating catabolic pathway expression, and could be the catabolite modulating factor whose existence was suggested by Jacques Monod and colleagues in 1976. Therefore, transcription factor Crp appears to work in concert with CytR and RpoH, serving a dual role in sensing both carbon availability and metabolic flux towards DNA and RNA. Our findings show how certain alterations in metabolite concentrations (associated with colony ageing and/or due to mutations in metabolic or regulatory genes) can drive the evolution in non-growing cells.

Microbial evolution often involves transient phenotypes and sequential development of multiple mutations of unclear relevance. Here, the authors show that the evolution of non-growing E. coli cells can be driven by alterations in pyrimidine nucleoside levels associated with colony ageing and/or due to mutations in metabolic or regulatory genes.

Insertion of prpoD_rpoS fragment enhances expression of recombinant protein by dps auto-inducible promoter in Escherichia coli

BACKGROUND

Nowadays, recombinant therapeutic proteins have been widely produced and consumed. For the safety and effectiveness of the protein production, an auto-inducible expression vector is required to replace inducer interference, which is uneconomic and could be harmful. In this research, an auto-inducible expression plasmid, pCAD₂_sod (a pBR322 derivate plasmid), which was under dps (RpoS-dependent gene) promoter control, was modified to provide RpoS at earlier phase. Hence, accumulates more target protein and resulting a new plasmid, pCAD₂⁺_sod. pCAD₂_sod had been constructed to automatically induces the expression of recombinant superoxide dismutase (SOD) from Staphylococcus equorum (rMnSODSeq) in the stationary growth phase of Escherichia coli. This work aimed to obtain pCAD₂⁺_sod and determine the expression level of rMnSODSeq on mRNA and protein level.

METHOD AND RESULTS

A synthetic rpoS coding region under rpoD promoter control (prpoD_rpoS) was inserted to pCAD₂_sod and generated pCAD₂⁺_sod. The rMnSODSeq (24.3 kDa) produced from pCAD₂⁺_sod was ~ 1.5 fold higher at 37 °C and more intense at 43 °C compared to that from pCAD₂_sod, likewise shifted to earlier phase (after 1 h of incubation), as shown in the SDS-PAGE. The dismutase activity was also retained after zymography assay. The mRNA level from pCAD₂⁺_sod was determined by qPCR and gave quantification cycle (Cq) values of cDNA lowest among others. It made the relative quantification (RQ) of the mRNA expression towards rho reference gene were high.

CONCLUSIONS

The prpoD_rpoS insertion shifts and increases the rMnSODSeq production from stationary to exponential phase. The pCAD₂⁺_sod plasmid is potential for further recombinant protein productions.

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.

The Context-Dependent Influence of Promoter Sequence Motifs on Transcription Initiation Kinetics and Regulation

The fitness of an individual bacterial cell is highly dependent upon the temporal tuning of gene expression levels when subjected to different environmental cues. Kinetic regulation of transcription initiation is a key step in modulating the levels of transcribed genes to promote bacterial survival. The initiation phase encompasses the binding of RNA polymerase (RNAP) to promoter DNA and a series of coupled protein-DNA conformational changes prior to entry into processive elongation. The time required to complete the initiation phase can vary by orders of magnitude and is ultimately dictated by the DNA sequence of the promoter. In this review, we aim to provide the required background to understand how promoter sequence motifs may affect initiation kinetics during promoter recognition and binding, subsequent conformational changes which lead to DNA opening around the transcription start site, and promoter escape. By calculating the steady-state flux of RNA production as a function of these effects, we illustrate that the presence/absence of a consensus promoter motif cannot be used in isolation to make conclusions regarding promoter strength. Instead, the entire series of linked, sequence-dependent structural transitions must be considered holistically. Finally, we describe how individual transcription factors take advantage of the broad distribution of sequence-dependent basal kinetics to either increase or decrease RNA flux.

Designing next generation recombinant protein expression platforms by modulating the cellular stress response in Escherichia coli

BACKGROUND

A cellular stress response (CSR) is triggered upon recombinant protein synthesis which acts as a global feedback regulator of protein expression. To remove this key regulatory bottleneck, we had previously proposed that genes that are up-regulated post induction could be part of the signaling pathways which activate the CSR. Knocking out some of these genes which were non-essential and belonged to the bottom of the E. coli regulatory network had provided higher expression of GFP and L-asparaginase.

RESULTS

We chose the best performing double knockout E. coli BW25113ΔelaAΔcysW and demonstrated its ability to enhance the expression of the toxic Rubella E1 glycoprotein by 2.5-fold by tagging it with sfGFP at the C-terminal end to better quantify expression levels. Transcriptomic analysis of this hyper-expressing mutant showed that a significantly lower proportion of genes got down-regulated post induction, which included genes for transcription, translation, protein folding and sorting, ribosome biogenesis, carbon metabolism, amino acid and ATP synthesis. This down-regulation which is a typical feature of the CSR was clearly blocked in the double knockout strain leading to its enhanced expression capability. Finally, we supplemented the expression of substrate uptake genes glpK and glpD whose down-regulation was not prevented in the double knockout, thus ameliorating almost all the negative effects of the CSR and obtained a further doubling in recombinant protein yields.

CONCLUSION

The study validated the hypothesis that these up-regulated genes act as signaling messengers which activate the CSR and thus, despite having no casual connection with recombinant protein synthesis, can improve cellular health and protein expression capabilities. Combining gene knockouts with supplementing the expression of key down-regulated genes can counter the harmful effects of CSR and help in the design of a truly superior host platform for recombinant protein expression.

Cyclic-di-GMP Regulates Autoaggregation Through the Putative Peptidoglycan Hydrolase, EagA, and Regulates Transcription of the znuABC Zinc Uptake Gene Cluster in Erwinia amylovora

Erwinia amylovora is the causal agent of fire blight, an economically impactful disease that affects apple and pear production worldwide. E. amylovora pathogenesis is comprised of distinct type III secretion-dependent and biofilm-dependent stages. Alterations in the intracellular levels of cyclic-di-GMP (c-di-GMP) regulate the transition between the different stages of infection in E. amylovora. We previously reported that hyper-elevation of c-di-GMP levels in E. amylovora Ea1189, resulting from the deletion of all three c-di-GMP specific phosphodiesterase genes (Ea1189ΔpdeABC), resulted in an autoaggregation phenotype. The two major exopolysaccharides, amylovoran and cellulose, were also shown to partially contribute to autoaggregation. In this study, we aimed to identify the c-di-GMP dependent factor(s) that contributes to autoaggregation. We conducted a transposon mutant screen in Ea1189ΔpdeABC and selected for loss of autoaggregation. Our search identified a peptidoglycan hydrolase, specifically, a D, D-endopeptidase of the metallopeptidase class, EagA (Erwiniaaggregation factor A), that was found to physiologically contribute to autoaggregation in a c-di-GMP dependent manner. The production of amylovoran was also positively affected by EagA levels. An eagA deletion mutant (Ea1189ΔeagA) was significantly reduced in virulence compared to the wild type E. amylovora Ea1189. eagA is part of the znuABC zinc uptake gene cluster and is located within an operon downstream of znuA. The znuAeagA/znuCB gene cluster was transcriptionally regulated by elevated levels of c-di-GMP as well as by the zinc-dependent transcriptional repressor Zur. We also observed that with an influx of Zn²⁺ in the environment, the transcription of the znuAeagA/znuBC gene cluster is regulated by both Zur and a yet to be characterized c-di-GMP dependent pathway.

Control of Transcription Initiation by Biased Thermal Fluctuations on Repetitive Genomic Sequences

In the process of transcription initiation by RNA polymerase, promoter DNA sequences affect multiple reaction pathways determining the productivity of transcription. However, the question of how the molecular mechanism of transcription initiation depends on the sequence properties of promoter DNA remains poorly understood. Here, combining the statistical mechanical approach with high-throughput sequencing results, we characterize abortive transcription and pausing during transcription initiation by Escherichia coli RNA polymerase at a genome-wide level. Our results suggest that initially transcribed sequences, when enriched with thymine bases, contain the signal for inducing abortive transcription, whereas certain repetitive sequence elements embedded in promoter regions constitute the signal for inducing pausing. Both signals decrease the productivity of transcription initiation. Based on solution NMR and in vitro transcription measurements, we suggest that repetitive sequence elements within the promoter DNA modulate the nonlocal base pair stability of its double-stranded form. This stability profoundly influences the reaction coordinates of the productive initiation via pausing.

A Decrease in Transcription Capacity Limits Growth Rate upon Translation Inhibition

Exposure of bacteria to sublethal concentrations of antibiotics can lead to bacterial adaptation and survival at higher doses of inhibitors, which in turn can lead to the emergence of antibiotic resistance. The presence of sublethal concentrations of antibiotics targeting translation results in an increase in the amount of ribosomes per cell but nonetheless a decrease in the cells’ growth rate. In this work, we have found that inhibition of ribosome activity can result in a decrease in the amount of free RNA polymerase available for transcription, thus limiting the protein expression rate via a different pathway than what was expected. This result can be explained by our observation that long genes, such as those coding for RNA polymerase subunits, have a higher probability of premature translation termination in the presence of ribosome inhibitors, while expression of short ribosomal genes is affected less, consistent with their increased concentration.

In bacterial cells, inhibition of ribosomes by sublethal concentrations of antibiotics leads to a decrease in the growth rate despite an increase in ribosome content. The limitation of ribosomal activity results in an increase in the level of expression from ribosomal promoters; this can deplete the pool of RNA polymerase (RNAP) that is available for the expression of nonribosomal genes. However, the magnitude of this effect remains to be quantified. Here, we use the change in the activity of constitutive promoters with different affinities for RNAP to quantify the change in the concentration of free RNAP. The data are consistent with a significant decrease in the amount of RNAP available for transcription of both ribosomal and nonribosomal genes. Results obtained with different reporter genes reveal an mRNA length dependence on the amount of full-length translated protein, consistent with the decrease in ribosome processivity affecting more strongly the translation of longer genes. The genes coding for the β and β' subunits of RNAP are among the longest genes in the Escherichia coli genome, while the genes coding for ribosomal proteins are among the shortest genes. This can explain the observed decrease in transcription capacity that favors the expression of genes whose promoters have a high affinity for RNAP, such as ribosomal promoters.

IMPORTANCE Exposure of bacteria to sublethal concentrations of antibiotics can lead to bacterial adaptation and survival at higher doses of inhibitors, which in turn can lead to the emergence of antibiotic resistance. The presence of sublethal concentrations of antibiotics targeting translation results in an increase in the amount of ribosomes per cell but nonetheless a decrease in the cells’ growth rate. In this work, we have found that inhibition of ribosome activity can result in a decrease in the amount of free RNA polymerase available for transcription, thus limiting the protein expression rate via a different pathway than what was expected. This result can be explained by our observation that long genes, such as those coding for RNA polymerase subunits, have a higher probability of premature translation termination in the presence of ribosome inhibitors, while expression of short ribosomal genes is affected less, consistent with their increased concentration.

Novel regulators of the csgD gene encoding the master regulator of biofilm formation in Escherichia coli K-12

Under stressful conditions, Escherichia coli forms biofilm for survival by sensing a variety of environmental conditions. CsgD, the master regulator of biofilm formation, controls cell aggregation by directly regulating the synthesis of Curli fimbriae. In agreement of its regulatory role, as many as 14 transcription factors (TFs) have so far been identified to participate in regulation of the csgD promoter, each monitoring a specific environmental condition or factor. In order to identify the whole set of TFs involved in this typical multi-factor promoter, we performed in this study ‘promoter-specific transcription-factor’ (PS-TF) screening in vitro using a set of 198 purified TFs (145 TFs with known functions and 53 hitherto uncharacterized TFs). A total of 48 TFs with strong binding to the csgD promoter probe were identified, including 35 known TFs and 13 uncharacterized TFs, referred to as Y-TFs. As an attempt to search for novel regulators, in this study we first analysed a total of seven Y-TFs, including YbiH, YdcI, YhjC, YiaJ, YiaU, YjgJ and YjiR. After analysis of curli fimbriae formation, LacZ-reporter assay, Northern-blot analysis and biofilm formation assay, we identified at least two novel regulators, repressor YiaJ (renamed PlaR) and activator YhjC (renamed RcdB), of the csgD promoter.

Antimicrobial resistance gene shuffling and a three-element mobilisation system in the monophasic Salmonella typhimurium strain ST1030

In this study we describe the genetic elements and the antimicrobial resistance units (RUs) harboured by the Salmonella Typhimurium monophasic variant 1,4,[5],12:i:- strain ST1030. Of the three identified RUs two were chromosomal, RU1 (IS26-bla_TEM-1-IS26-strAB-sul2- IS26) and RU2 (IS26-tetR(B)-tetA(B)-ΔIS26), and one, RU3 (a sul3-associated class 1 integron with cassette array dfrA12-orfF-aadA2-cmlA1-aadA1), was embedded in a Tn21-derived element harboured by the conjugative I1 plasmid pST1030-1A. IS26 elements mediated the antimicrobial resistance gene (ARG) shuffling and this gave rise to pST1030-1A derivatives with different sets of ARGs. ST1030 also harboured two ColE1-like plasmids of which one, pST1030-2A, was mobilisable and the target of an intracellular translocation of the Tn21-derived element; the second (pST1030-3) was an orphan mob-associated oriT plasmid co-transferred with pST1030-1A and pST1030-2A. pST1030-2A and pST1030-3 also carried a parA gene and a type III restriction modification system, respectively. Overall analysis of our data reinforces the role played by IS26, Tn21-derived elements and non-conjugative plasmids in the spread of ARGs and supplies the first evidence, at least in Salmonella, for the identification of a natural isolate harbouring a three-element mobilisation system in the same cell.

Benchmarking Bacterial Promoter Prediction Tools: Potentialities and Limitations

The correct mapping of promoter elements is a crucial step in microbial genomics. Also, when combining new DNA elements into synthetic sequences, predicting the potential generation of new promoter sequences is critical. Over the last years, many bioinformatics tools have been created to allow users to predict promoter elements in a sequence or genome of interest. Here, we assess the predictive power of some of the main prediction tools available using well-defined promoter data sets. Using Escherichia coli as a model organism, we demonstrated that while some tools are biased toward AT-rich sequences, others are very efficient in identifying real promoters with low false-negative rates. We hope the potentials and limitations presented here will help the microbiology community to choose promoter prediction tools among many available alternatives.

The promoter region is a key element required for the production of RNA in bacteria. While new high-throughput technology allows massively parallel mapping of promoter elements, we still mainly rely on bioinformatics tools to predict such elements in bacterial genomes. Additionally, despite many different prediction tools having become popular to identify bacterial promoters, no systematic comparison of such tools has been performed. Here, we performed a systematic comparison between several widely used promoter prediction tools (BPROM, bTSSfinder, BacPP, CNNProm, IBBP, Virtual Footprint, iPro70-FMWin, 70ProPred, iPromoter-2L, and MULTiPly) using well-defined sequence data sets and standardized metrics to determine how well those tools performed related to each other. For this, we used data sets of experimentally validated promoters from Escherichia coli and a control data set composed of randomly generated sequences with similar nucleotide distributions. We compared the performance of the tools using metrics such as specificity, sensitivity, accuracy, and Matthews correlation coefficient (MCC). We show that the widely used BPROM presented the worse performance among the compared tools, while four tools (CNNProm, iPro70-FMWin, 70ProPred, and iPromoter-2L) offered high predictive power. Of these tools, iPro70-FMWin exhibited the best results for most of the metrics used. We present here some potentials and limitations of available tools, and we hope that future work can build upon our effort to systematically characterize this useful class of bioinformatics tools.

IMPORTANCE The correct mapping of promoter elements is a crucial step in microbial genomics. Also, when combining new DNA elements into synthetic sequences, predicting the potential generation of new promoter sequences is critical. Over the last years, many bioinformatics tools have been created to allow users to predict promoter elements in a sequence or genome of interest. Here, we assess the predictive power of some of the main prediction tools available using well-defined promoter data sets. Using Escherichia coli as a model organism, we demonstrated that while some tools are biased toward AT-rich sequences, others are very efficient in identifying real promoters with low false-negative rates. We hope the potentials and limitations presented here will help the microbiology community to choose promoter prediction tools among many available alternatives.

FEDS: a Novel Fluorescence-Based High-Throughput Method for Measuring DNA Supercoiling In Vivo

DNA supercoiling (DS) is essential for life because it controls critical processes, including transcription, replication, and recombination. Current methods to measure DNA supercoiling in vivo are laborious and unable to examine single cells. Here, we report a method for high-throughput measurement of bacterial DNA supercoiling in vivoFluorescent evaluation of DNA supercoiling (FEDS) utilizes a plasmid harboring the gene for a green fluorescent protein transcribed by a discovered promoter that responds exclusively to DNA supercoiling and the gene for a red fluorescent protein transcribed by a constitutive promoter as the internal standard. Using FEDS, we uncovered single-cell heterogeneity in DNA supercoiling and established that, surprisingly, population-level decreases in DNA supercoiling result from a low-mean/high-variance DNA supercoiling subpopulation rather than from a homogeneous shift in supercoiling of the whole population. In addition, we identified a regulatory loop in which a gene that decreases DNA supercoiling is transcriptionally repressed when DNA supercoiling increases.IMPORTANCE DNA represents the chemical support of genetic information in all forms of life. In addition to its linear sequence of nucleotides, it bears critical information in its structure. This information, called DNA supercoiling, is central to all fundamental DNA processes, such as transcription and replication, and defines cellular physiology. Unlike reading of a nucleotide sequence, DNA supercoiling determinations have been laborious. We have now developed a method for rapid measurement of DNA supercoiling and established its utility by identifying a novel regulator of DNA supercoiling in the bacterium Salmonella enterica as well as behaviors that could not have been discovered with current methods.

Activity evaluation of glycolytic promoters from Escherichia coli and application for mevalonate biosynthesis

Promoters are the most important tools to control and regulate gene expression in synthetic biology and metabolic engineering. The expression of target genes in Escherichia coli is usually controlled by inducible promoters which may cause excessive metabolic load on the host and may be uneconomical due to inducer cost. Therefore, it is important to identify more constitutive promoters that are capable of avoiding these limitations. In this study, using the monomeric red fluorescent protein (mRFP) as the reporter gene, ten promoters from the glycolytic pathway of E. coli were cloned and characterized. We found that glycolytic promoters exhibited the advantages of constitutive promoters and higher strength compared with the commonly used inducible promoter Plac. We further introduced glycolytic promoters into the mevalonate biosynthesis system. The maximum mevalonate titer produced by engineered E. coli under the control of glycolytic promoters was obviously higher than that under the control of Plac, indicating the superiority of glycolytic promoter for the metabolic engineering of E. coli. This set of glycolytic promoters significantly expands the range of engineering tools available for E. coli and can be applied in future metabolic engineering studies.

Tools to map target genes of bacterial two-component system response regulators

Studies on bacterial physiology are incomplete without knowledge of the signalling and regulatory systems that a bacterium uses to sense and respond to its environment. Two-component systems (TCSs) are among the most prevalent bacterial signalling systems, and they control essential and secondary physiological processes; however, even in model organisms, we lack a complete understanding of the signals sensed, the phosphotransfer partners and the functions regulated by these systems. In this review, we discuss several tools to map the genes targeted by transcriptionally acting TCSs. Many of these tools have been used for studying individual TCSs across diverse species, but systematic approaches to delineate entire signalling networks have been very few. Since genome sequences and high-throughput technologies are now readily available, the methods presented here can be applied to characterize the entire DNA-binding TCS signalling network in any bacterial species and are especially useful for non-model environmental bacteria.

Excessive Promoters as Silencers of Genes Horizontally Acquired by Escherichia coli

Horizontally acquired genes are usually transcriptionally inactive, although most of them are associated with genomic loci enriched with promoter-like sequences forming “promoter islands.” We hypothesized that lateral DNA transfer induces local mutagenesis, accumulating AT base pairs and creating promoter-like sequences, whose occupancy with RNA polymerase and a specific silencer H-NS suppresses the transcription of foreign genes. Error-prone mutagenesis was implemented for the “promoter island” of a foreign gene appY and the promoter region of an inherent gene dps. Derivatives with changed transcriptional activity were selected using a reporter plasmid pET28_eGFP. Only one cycle of mutagenesis with negative selection suppressed the activity of the main dps promoter to the background level due to a single substitution in its -10 element, while positive selection gave a sequence with improved -35 element, thus testifying feasibility of the approach. The same suppression for appY was achieved by three cycles, while eightfold transcription activation required nine iterations of mutagenesis. In both cases, the number of potential start points decreased resulting in an ordinary regulatory region with only one dominant promoter in the case of positive selection. Efficiency of H-NS binding remained virtually unchanged in all mutant constructs. Based on these findings we conclude that excessive promoters can adversely affect transcription by providing a platform for interference between several RNA polymerase molecules, which can act as a silencer at promoter-dense regions.

Protein-based biorefining driven by nitrogen-responsive transcriptional machinery

BACKGROUND

Protein-based bioconversion has been demonstrated as a sustainable approach to produce higher alcohols and ammonia fertilizers. However, owing to the switchover from transcription mediated by the bacterial RNA polymerase σ⁷⁰ to that mediated by alternative σ factors, the biofuel production driven by σ⁷⁰-dependent promoters declines rapidly once cells enter the stationary phase or encounter stresses. To enhance biofuel production, in this study the growth phase-independent and nitrogen-responsive transcriptional machinery mediated by the σ⁵⁴ is exploited to drive robust protein-to-fuel conversion.

RESULTS

We demonstrated that disrupting the Escherichia coli ammonia assimilation pathways driven by glutamate dehydrogenase and glutamine synthetase could sustain the activity of σ⁵⁴-mediated transcription under ammonia-accumulating conditions. In addition, two σ⁵⁴-dependent promoters, argTp and glnAp2, were identified as suitable candidates for driving pathway expression. Using these promoters, biofuel production from proteins was shown to persist to the stationary phase, with the net production in the stationary phase being 1.7-fold higher than that derived from the optimal reported σ⁷⁰-dependent promoter P _LlacO₁. Biofuel production reaching levels 1.3- to 3.4-fold higher than those of the σ⁷⁰-dependent promoters was also achieved by argTp and glnAp2 under stressed conditions. Moreover, the σ⁵⁴-dependent promoters realized more rapid and stable production than that of σ⁷⁰-dependent promoters during fed-batch fermentation, producing up to 4.78 g L ^- 1 of total biofuels.

CONCLUSIONS

These results suggested that the nitrogen-responsive transcriptional machinery offers the potential to decouple production from growth, highlighting this system as a novel candidate to realize growth phase-independent and stress-resistant biofuel production.

Density of σ70 promoter-like sites in the intergenic regions dictates the redistribution of RNA polymerase during osmotic stress in Escherichia coli

RNA polymerase (RNAP), the transcription machinery, shows dynamic binding across the genomic DNA under different growth conditions. The genomic features that selectively redistribute the limited RNAP molecules to dictate genome-wide transcription in response to environmental cues remain largely unknown. We chose the bacterial osmotic stress response model to determine genomic features that direct genome-wide redistribution of RNAP during the stress. Genomic mapping of RNAP and transcriptome profiles corresponding to the different temporal states after salt shock were determined. We found rapid redistribution of RNAP across the genome, primarily at σ70 promoters. Three subsets of genes exhibiting differential salt sensitivities were identified. Sequence analysis using an information-theory based σ70 model indicates that the intergenic regions of salt-responsive genes are enriched with a higher density of σ70 promoter-like sites than those of salt-sensitive genes. In addition, the density of promoter-like sites has a positive linear correlation with RNAP binding at different salt concentrations. The RNAP binding contributed by the non-initiating promoter-like sites is important for gene transcription at high salt concentration. Our study demonstrates that hyperdensity of σ70 promoter-like sites in the intergenic regions of salt-responsive genes drives the RNAP redistribution for reprograming the transcriptome to counter osmotic stress.

Molecular dissection of the replication system of plasmid pIGRK encoding two in-frame Rep proteins with antagonistic functions

BACKGROUND

Gene overlapping is a frequent phenomenon in microbial genomes. Excluding so-called "trivial overlapping", there are significant implications of such genetic arrangements, including regulation of gene expression and modification of protein activity. It is also postulated that, besides gene duplication, the appearance of overlapping genes (OGs) is one of the most important factors promoting a genome's novelty and evolution. OGs coding for in-frame proteins with different functions are a particularly interesting case. In this study we identified and characterized two in-frame proteins encoded by OGs on plasmid pIGRK from Klebsiella pneumoniae, a representative of the newly distinguished pHW126 plasmid family.

RESULTS

A single repR locus located within the replication system of plasmid pIGRK encodes, in the same frame, two functional polypeptides: a full-length RepR protein and a RepR' protein (with N-terminal truncation) translated from an internal START codon. Both proteins form homodimers, and interact with diverse DNA regions within the plasmid replication origin and repR promoter operator. Interestingly, RepR and RepR' have opposing functions - RepR is crucial for initiation of pIGRK replication, while RepR' is a negative regulator of this process. Nevertheless, both proteins act cooperatively as negative transcriptional regulators of their own expression.

CONCLUSIONS

Regulation of the initiation of pIGRK replication is a complex process in which a major role is played by two in-frame proteins with antagonistic functions. In-frame encoded Rep proteins are uncommon, having been described in only a few plasmids. This is the first description of such proteins in a plasmid of the pHW126 family.

Contribution of σ70 and σN Factors to Expression of Class II pilE in Neisseria meningitidis

Neisseria meningitidis expresses multicomponent organelles called type four pili (Tfp), which are key virulence factors required for attachment to human cells during carriage and disease. Pilin (PilE) is the main component of Tfp, and N. meningitidis isolates either have a class I pilE locus and express pilins that undergo antigenic variation or have a class II pilE locus and express invariant pilins. The transcriptional regulation of class I pilE has been studied in both N. meningitidis and Neisseria gonorrhoeae, while the control of expression of class II pilE has been elucidated in the nonpathogenic species Neisseria elongata However, the factors that govern the regulation of the class II pilE gene in N. meningitidis are not known. In this work, we have bioinformatically and experimentally identified the class II pilE promoter. We confirmed the presence of conserved σ⁷⁰ and σ^N-dependent promoters upstream of pilE in a collection of meningococcal genomes and demonstrated that class II pilE expression initiates from the σ⁷⁰ family-dependent promoter. By deletion or overexpression of sigma factors, we showed that σ^N, σ^H, and σ^E do not affect class II pilin expression. These findings are consistent with a role of the housekeeping σ^D in expression of this important component of Tfp. Taken together, our data indicate that the σ-dependent network responsible for the expression of class II pilE has been selected to maintain pilE expression, consistent with the essential roles of Tfp in colonization and pathogenesis.IMPORTANCE The type four pilus (Tfp) of Neisseria meningitidis contributes to fundamental processes such as adhesion, transformation, and disease pathology. Meningococci express one of two distinct classes of Tfp (class I or class II), which can be distinguished antigenically or by the major subunit (pilE) locus and its genetic context. The factors that govern transcription of the class II pilE gene are not known, even though it is present in isolates that cause epidemic disease. Here we show that the transcription of class II pilE is maintained throughout growth and under different stress conditions and is driven by a σ⁷⁰-dependent promoter. This is distinct from Tfp regulation in nonpathogenic Neisseria spp. and may confer an advantage during host-cell interaction and infection.

Development of a Vanillate Biosensor for the Vanillin Biosynthesis Pathway in E. coli

The engineered de novo vanillin biosynthesis pathway constructed in Escherichia coli is industrially relevant but limited by the reaction catalyzed by catechol O-methyltransferase, which is intended to catalyze the conversion of protocatechuate to vanillate. To identify alternative O-methyltransferases, we constructed a vanillate sensor based on the Caulobacter crescentus VanR-VanO system. Using an E. coli promoter library, we achieved greater than 14-fold dynamic range in our best rationally constructed sensor. We found that this construct and an evolved variant demonstrate remarkable substrate selectivity, exhibiting no detectable response to the regioisomer byproduct isovanillate and minimal response to structurally similar pathway intermediates. We then harnessed the evolved biosensor to conduct rapid bioprospecting of natural catechol O-methyltransferases and identified three previously uncharacterized but active O-methyltransferases. Collectively, these efforts enrich our knowledge of how biosensing can aid metabolic engineering and constitute the foundation for future improvements in vanillin pathway productivity.

Sequence Determinants Spanning -10 Motif and Spacer Region Implicated in Unique Ehrlichia chaffeensis Sigma 32-Dependent Promoter Activity of dnaK Gene

Ehrlichia chaffeensis is an obligate intracellular tick-borne bacterium that causes human monocytic ehrlichiosis. Studying Ehrlichia gene regulation is challenge, as this and related rickettsiales lack natural plasmids and mutagenesis experiments are of a limited scope. E. chaffeensis contains only two sigma factors, σ³² and σ⁷⁰. We previously developed Escherichia coli surrogate system to study transcriptional regulation from RNA polymerase (RNAP) containing Ehrlichia σ³² or σ⁷⁰. We reported that RNAP binding motifs of E. chaffeensis genes recognized by σ³² or σ⁷⁰ share extensive homology and that transcription may be initiated by either one of the sigma factors, although transcriptional efficiencies differ. In the current study, we investigated mapping the E. chaffeensis dnaK gene promoter using the pathogen σ³² expressed in E. coli lacking its native σ³². The E. coli surrogate system and our previously described in vitro transcription system aided in defining the unique −10 motif and spacer sequence of the dnaK promoter. We also mapped σ³² amino acids/domains engaged in its promoter regulation in E. chaffeensis. The data reported in this study demonstrate that the −10 and −35 motifs and spacer sequence located between the two motifs of dnaK promoter are critical for the RNAP function. Further, we mapped the importance of all six nucleotide positions of the −10 motif and identified critical determinants within it. In addition, we reported that the lack of C-rich sequence upstream to the −10 motif is unique in driving the pathogen-specific transcription by its σ³² from dnaK gene promoter. This is the first study in defining an E. chaffeensis σ³²-dependent promoter and it offers insights about how this and other related rickettsial pathogens regulate stress response genes.

Single-target regulators form a minor group of transcription factors in Escherichia coli K-12

The identification of regulatory targets of all TFs is critical for understanding the entire network of the genome regulation. The lac regulon of Escherichia coli K-12 W3110 is composed of the lacZYA operon and its repressor lacI gene, and has long been recognized as the seminal model of transcription regulation in bacteria with only one highly preferred target. After the Genomic SELEX screening in vitro of more than 200 transcription factors (TFs) from E. coli K-12, however, we found that most TFs regulate multiple target genes. With respect to the number of regulatory targets, a total of these 200 E. coli TFs form a hierarchy ranging from a single target to as many as 1000 targets. Here we focus a total of 13 single-target TFs, 9 known TFs (BetI, KdpE, LacI, MarR, NanR, RpiR, TorR, UlaR and UxuR) and 4 uncharacterized TFs (YagI, YbaO, YbiH and YeaM), altogether forming only a minor group of TFs in E. coli. These single-target TFs were classified into three groups based on their functional regulation.

Deep genome annotation of the opportunistic human pathogen Streptococcus pneumoniae D39

A precise understanding of the genomic organization into transcriptional units and their regulation is essential for our comprehension of opportunistic human pathogens and how they cause disease. Using single-molecule real-time (PacBio) sequencing we unambiguously determined the genome sequence of Streptococcus pneumoniae strain D39 and revealed several inversions previously undetected by short-read sequencing. Significantly, a chromosomal inversion results in antigenic variation of PhtD, an important surface-exposed virulence factor. We generated a new genome annotation using automated tools, followed by manual curation, reflecting the current knowledge in the field. By combining sequence-driven terminator prediction, deep paired-end transcriptome sequencing and enrichment of primary transcripts by Cappable-Seq, we mapped 1015 transcriptional start sites and 748 termination sites. We show that the pneumococcal transcriptional landscape is complex and includes many secondary, antisense and internal promoters. Using this new genomic map, we identified several new small RNAs (sRNAs), RNA switches (including sixteen previously misidentified as sRNAs), and antisense RNAs. In total, we annotated 89 new protein-encoding genes, 34 sRNAs and 165 pseudogenes, bringing the S. pneumoniae D39 repertoire to 2146 genetic elements. We report operon structures and observed that 9% of operons are leaderless. The genome data are accessible in an online resource called PneumoBrowse (https://veeninglab.com/pneumobrowse) providing one of the most complete inventories of a bacterial genome to date. PneumoBrowse will accelerate pneumococcal research and the development of new prevention and treatment strategies.

Recent Advances in Understanding σ70-Dependent Transcription Initiation Mechanisms

Prokaryotic transcription is one of the most studied biological systems, with relevance to many fields including the development and use of antibiotics, the construction of synthetic gene networks, and the development of many cutting-edge methodologies. Here, we discuss recent structural, biochemical, and single-molecule biophysical studies targeting the mechanisms of transcription initiation in bacteria, including the formation of the open complex, the reaction of initial transcription, and the promoter escape step that leads to elongation. We specifically focus on the mechanisms employed by the RNA polymerase holoenzyme with the housekeeping sigma factor σ⁷⁰. The recent progress provides answers to long-held questions, identifies intriguing new behaviours, and opens up fresh questions for the field of transcription.

Role and Function of Class III LitR, a Photosensor Homolog from Burkholderia multivorans

The LitR/CarH protein family is an adenosyl B₁₂ (AdoB₁₂)-dependent photoreceptor family with DNA-binding activity, and its homologs are widely distributed in the genomes of diverse bacterial genera. In this investigation, we studied the role and functions of a LitR homolog from a Gram-negative soil bacterium, Burkholderia multivorans, which does not possess an AdoB₁₂-binding domain. Transcriptome analysis indicated the existence of 19 light-induced genes, including folE2, cfaB, litS, photolyase gene phrB2, and cryB, located in the region flanking litR Disruption of litR caused constitutive expression of all the light-inducible genes, while mutation in the light-induced sigma factor gene, litS, abolished the transcription of the phrB2 operon and the cfa operon, indicating that LitR and LitS play a central role in light-inducible transcription. A gel shift assay showed that recombinant protein LitR specifically binds to the promoter regions of litR and the folE2 operon, and its binding was weakened by UV-A illumination. LitR absorbs light at maximally near 340 nm and exhibited a photocyclic response and light-dependent dissociation of multimer into tetramer. The litR mutant produced a 20-fold-higher intracellular level of folate than that of the wild-type strain. Thus, the evidence suggests that LitR light-dependently regulates the transcription of litR itself and the folE2 operon, resulting in the production of folate, and then the expressed RNA polymerase complex containing σ^LitS directs the transcription of the phrB2 operon and the cfa operon. These light-dependent characteristics suggest that class III LitR, in complex with a UV-A-absorbing molecule, follows a novel light-sensing mechanism.IMPORTANCE Members of the LitR/CarH family are adenosyl B₁₂-based photosensory transcriptional regulator involved in light-inducible carotenoid production in nonphototrophic bacteria. Our study provides the first evidence of the involvement of a class III LitR, which lacks an adenosyl B₁₂-binding domain in the light response of Burkholderia multivorans belonging to betaproteobacteria. Our biochemical analysis suggests that class III LitR protein exhibits features as a photosensor including absorption of light at the UV-A region (λ_max = ca. 340 nm), photocyclic response, and light-dependent dissociation. This suggests that class III LitR associates with a UV-A-absorbing molecule, and it has a photosensing mechanism distinguishable from that of the B₁₂-based type.

Added value of autoregulation and multi-step kinetics of transcription initiation

Bacterial gene expression regulation occurs mostly during transcription, which has two main rate-limiting steps: the close complex formation, when the RNA polymerase binds to an active promoter, and the subsequent open complex formation, after which it follows elongation. Tuning these steps' kinetics by the action of e.g. transcription factors, allows for a wide diversity of dynamics. For example, adding autoregulation generates single-gene circuits able to perform more complex tasks. Using stochastic models of transcription kinetics with empirically validated parameter values, we investigate how autoregulation and the multi-step transcription initiation kinetics of single-gene autoregulated circuits can be combined to fine-tune steady state mean and cell-to-cell variability in protein expression levels, as well as response times. Next, we investigate how they can be jointly tuned to control complex behaviours, namely, time counting, switching dynamics and memory storage. Overall, our finding suggests that, in bacteria, jointly regulating a single-gene circuit's topology and the transcription initiation multi-step dynamics allows enhancing complex task performance.

Coordinated Hibernation of Transcriptional and Translational Apparatus during Growth Transition of Escherichia coli to Stationary Phase

During the growth transition of E. coli from exponential phase to stationary, the genome expression pattern is altered markedly. For this alteration, the transcription apparatus is altered by binding of anti-sigma factor Rsd to the RpoD sigma factor for sigma factor replacement, while the translation machinery is modulated by binding of RMF to 70S ribosome to form inactive ribosome dimer. Using the PS-TF screening system, a number of TFs were found to bind to both the rsd and rmf promoters, of which the regulatory roles of 5 representative TFs (one repressor ArcA and the four activators McbR, RcdA, SdiA, and SlyA) were analyzed in detail. The results altogether indicated the involvement of a common set of TFs, each sensing a specific environmental condition, in coordinated hibernation of the transcriptional and translational apparatus for adaptation and survival under stress conditions.

In the process of Escherichia coli K-12 growth from exponential phase to stationary, marked alteration takes place in the pattern of overall genome expression through modulation of both parts of the transcriptional and translational apparatus. In transcription, the sigma subunit with promoter recognition properties is replaced from the growth-related factor RpoD by the stationary-phase-specific factor RpoS. The unused RpoD is stored by binding with the anti-sigma factor Rsd. In translation, the functional 70S ribosome is converted to inactive 100S dimers through binding with the ribosome modulation factor (RMF). Up to the present time, the regulatory mechanisms of expression of these two critical proteins, Rsd and RMF, have remained totally unsolved. In this study, attempts were made to identify the whole set of transcription factors involved in transcription regulation of the rsd and rmf genes using the newly developed promoter-specific transcription factor (PS-TF) screening system. In the first screening, 74 candidate TFs with binding activity to both of the rsd and rmf promoters were selected from a total of 194 purified TFs. After 6 cycles of screening, we selected 5 stress response TFs, ArcA, McbR, RcdA, SdiA, and SlyA, for detailed analysis in vitro and in vivo of their regulatory roles. Results indicated that both rsd and rmf promoters are repressed by ArcA and activated by McbR, RcdA, SdiA, and SlyA. We propose the involvement of a number of TFs in simultaneous and coordinated regulation of the transcriptional and translational apparatus. By using genomic SELEX (gSELEX) screening, each of the five TFs was found to regulate not only the rsd and rmf genes but also a variety of genes for growth and survival.

IMPORTANCE During the growth transition of E. coli from exponential phase to stationary, the genome expression pattern is altered markedly. For this alteration, the transcription apparatus is altered by binding of anti-sigma factor Rsd to the RpoD sigma factor for sigma factor replacement, while the translation machinery is modulated by binding of RMF to 70S ribosome to form inactive ribosome dimer. Using the PS-TF screening system, a number of TFs were found to bind to both the rsd and rmf promoters, of which the regulatory roles of 5 representative TFs (one repressor ArcA and the four activators McbR, RcdA, SdiA, and SlyA) were analyzed in detail. The results altogether indicated the involvement of a common set of TFs, each sensing a specific environmental condition, in coordinated hibernation of the transcriptional and translational apparatus for adaptation and survival under stress conditions.

Mining Novel Constitutive Promoter Elements in Soil Metagenomic Libraries in Escherichia coli

Although functional metagenomics has been widely employed for the discovery of genes relevant to biotechnology and biomedicine, its potential for assessing the diversity of transcriptional regulatory elements of microbial communities has remained poorly explored. Here, we experimentally mined novel constitutive promoter sequences in metagenomic libraries by combining a bi-directional reporter vector, high-throughput fluorescence assays and predictive computational methods. Through the expression profiling of fluorescent clones from two independent soil sample libraries, we have analyzed the regulatory dynamics of 260 clones with candidate promoters as a set of active metagenomic promoters in the host Escherichia coli. Through an in-depth analysis of selected clones, we were able to further explore the architecture of metagenomic fragments and to report the presence of multiple promoters per fragment with a dominant promoter driving the expression profile. These approaches resulted in the identification of 33 novel active promoters from metagenomic DNA originated from very diverse phylogenetic groups. The in silico and in vivo analysis of these individual promoters allowed the generation of a constitutive promoter consensus for exogenous sequences recognizable by E. coli in metagenomic studies. The results presented here demonstrates the potential of functional metagenomics for exploring environmental bacterial communities as a source of novel regulatory genetic parts to expand the toolbox for microbial engineering.

Regulation of Global Transcription in Escherichia coli by Rsd and 6S RNA

In Escherichia coli, the sigma factor σ⁷⁰ directs RNA polymerase to transcribe growth-related genes, while σ³⁸ directs transcription of stress response genes during stationary phase. Two molecules hypothesized to regulate RNA polymerase are the protein Rsd, which binds to σ⁷⁰, and the non-coding 6S RNA which binds to the RNA polymerase-σ⁷⁰ holoenzyme. Despite multiple studies, the functions of Rsd and 6S RNA remain controversial. Here we use RNA-Seq in five phases of growth to elucidate their function on a genome-wide scale. We show that Rsd and 6S RNA facilitate σ³⁸ activity throughout bacterial growth, while 6S RNA also regulates widely different genes depending upon growth phase. We discover novel interactions between 6S RNA and Rsd and show widespread expression changes in a strain lacking both regulators. Finally, we present a mathematical model of transcription which highlights the crosstalk between Rsd and 6S RNA as a crucial factor in controlling sigma factor competition and global gene expression.

Genome-wide relationship between R-loop formation and antisense transcription in Escherichia coli

Transcription termination by Rho is essential for viability in various bacteria, including some major pathogens. Since Rho acts by targeting nascent RNAs that are not simultaneously translated, it also regulates antisense transcription. Here we show that RNase H-deficient mutants of Escherichia coli exhibit heightened sensitivity to the Rho inhibitor bicyclomycin, and that Rho deficiency provokes increased formation of RNA-DNA hybrids (R-loops) which is ameliorated by expression of the phage T4-derived R-loop helicase UvsW. We also provide evidence that in Rho-deficient cells, R-loop formation blocks subsequent rounds of antisense transcription at more than 500 chromosomal loci. Hence these antisense transcripts, which can extend beyond 10 kb in their length, are only detected when Rho function is absent or compromised and the UvsW helicase is concurrently expressed. Thus the potential for antisense transcription in bacteria is much greater than hitherto recognized; and the cells are able to retain viability even when nearly one-quarter of their total non-rRNA abundance is accounted for by antisense transcripts, provided that R-loop formation from them is curtailed.

Random sequences rapidly evolve into de novo promoters

How new functions arise de novo is a fundamental question in evolution. We studied de novo evolution of promoters in Escherichia coli by replacing the lac promoter with various random sequences of the same size (~100 bp) and evolving the cells in the presence of lactose. We find that ~60% of random sequences can evolve expression comparable to the wild-type with only one mutation, and that ~10% of random sequences can serve as active promoters even without evolution. Such a short mutational distance between random sequences and active promoters may improve the evolvability, yet may also lead to accidental promoters inside genes that interfere with normal expression. Indeed, our bioinformatic analyses indicate that E. coli was under selection to reduce accidental promoters inside genes by avoiding promoter-like sequences. We suggest that a low threshold for functionality balanced by selection against undesired targets can increase the evolvability by making new beneficial features more accessible.

Bacterial promoters initiate gene transcription and have distinct sequence features. Here, the authors show that random sequences that contain no information are just on the verge of functioning as promoters in Escherichia coli.

Piggy: a rapid, large-scale pan-genome analysis tool for intergenic regions in bacteria

Background

The concept of the "pan-genome," which refers to the total complement of genes within a given sample or species, is well established in bacterial genomics. Rapid and scalable pipelines are available for managing and interpreting pan-genomes from large batches of annotated assemblies. However, despite overwhelming evidence that variation in intergenic regions in bacteria can directly influence phenotypes, most current approaches for analyzing pan-genomes focus exclusively on protein-coding sequences.

Findings

To address this we present Piggy, a novel pipeline that emulates Roary except that it is based only on intergenic regions. A key utility provided by Piggy is the detection of highly divergent ("switched") intergenic regions (IGRs) upstream of genes. We demonstrate the use of Piggy on large datasets of clinically important lineages of Staphylococcus aureus and Escherichia coli.

Conclusions

For S. aureus, we show that highly divergent (switched) IGRs are associated with differences in gene expression and we establish a multilocus reference database of IGR alleles (igMLST; implemented in BIGSdb).