Regulation of plasmid-encoded isoprene metabolism in Rhodococcus, a representative of an important link in the global isoprene cycle.

Emissions of biogenic volatile organic compounds (VOCs) form an important part of the global carbon cycle, comprising a significant proportion of net ecosystem productivity. They impact atmospheric chemistry and contribute directly and indirectly to greenhouse gases. Isoprene, emitted largely from plants, comprises one third of total VOCs, yet in contrast to methane, which is released in similar quantities, we know little of its biodegradation. Here, we report the genome of an isoprene degrading isolate, Rhodococcus sp. AD45, and, using mutagenesis shows that a plasmid-encoded soluble di-iron centre isoprene monooxygenase (IsoMO) is essential for isoprene metabolism. Using RNA sequencing (RNAseq) to analyse cells exposed to isoprene or epoxyisoprene in a substrate-switch time-course experiment, we show that transcripts from 22 contiguous genes, including those encoding IsoMO, were highly upregulated, becoming among the most abundant in the cell and comprising over 25% of the entire transcriptome. Analysis of gene transcription in the wild type and an IsoMO-disrupted mutant strain showed that epoxyisoprene, or a subsequent product of isoprene metabolism, rather than isoprene itself, was the inducing molecule. We provide a foundation of molecular data for future research on the environmental biological consumption of this important, climate-active compound.

Design of orthogonal genetic switches based on a crosstalk map of σs, anti-σs, and promoters.

Cells react to their environment through gene regulatory networks. Network integrity requires minimization of undesired crosstalk between their biomolecules. Similar constraints also limit the use of regulators when building synthetic circuits for engineering applications. Here, we mapped the promoter specificities of extracytoplasmic function (ECF) σs as well as the specificity of their interaction with anti-σs. DNA synthesis was used to build 86 ECF σs (two from every subgroup), their promoters, and 62 anti-σs identified from the genomes of diverse bacteria. A subset of 20 σs and promoters were found to be highly orthogonal to each other. This set can be increased by combining the -35 and -10 binding domains from different subgroups to build chimeras that target sequences unrepresented in any subgroup. The orthogonal σs, anti-σs, and promoters were used to build synthetic genetic switches in Escherichia coli. This represents a genome-scale resource of the properties of ECF σs and a resource for synthetic biology, where this set of well-characterized regulatory parts will enable the construction of sophisticated gene expression programs.

Predicting the strength of UP-elements and full-length E. coli σE promoters.

Predicting the location and strength of promoters from genomic sequence requires accurate sequenced-based promoter models. We present the first model of a full-length bacterial promoter, encompassing both upstream sequences (UP-elements) and core promoter modules, based on a set of 60 promoters dependent on σ(E), an alternative ECF-type σ factor. UP-element contribution, best described by the length and frequency of A- and T-tracts, in combination with a PWM-based core promoter model, accurately predicted promoter strength both in vivo and in vitro. This model also distinguished active from weak/inactive promoters. Systematic examination of promoter strength as a function of RNA polymerase (RNAP) concentration revealed that UP-element contribution varied with RNAP availability and that the σ(E) regulon is comprised of two promoter types, one of which is active only at high concentrations of RNAP. Distinct promoter types may be a general mechanism for increasing the regulatory capacity of the ECF group of alternative σ's. Our findings provide important insights into the sequence requirements for the strength and function of full-length promoters and establish guidelines for promoter prediction and for forward engineering promoters of specific strengths.

Unexpected stress-reducing effect of PhaP, a poly(3-hydroxybutyrate) granule-associated protein, in Escherichia coli.

Phasins (PhaP) are proteins normally associated with granules of poly(3-hydroxybutyrate) (PHB), a biodegradable polymer accumulated by many bacteria as a reserve molecule. These proteins enhance growth and polymer production in natural and recombinant PHB producers. It has been shown that the production of PHB causes stress in recombinant Escherichia coli, revealed by an increase in the concentrations of several heat stress proteins. In this work, quantitative reverse transcription (qRT)-PCR analysis was used to study the effect of PHB accumulation, and that of PhaP from Azotobacter sp. strain FA8, on the expression of stress-related genes in PHB-producing E. coli. While PHB accumulation was found to increase the transcription of dnaK and ibpA, the expression of these genes and of groES, groEL, rpoH, dps, and yfiD was reduced, when PhaP was coexpressed, to levels even lower than those detected in the non-PHB-accumulating control. These results demonstrated the protective role of PhaP in PHB-synthesizing E. coli and linked the effects of the protein to the expression of stress-related genes, especially ibpA. The effect of PhaP was also analyzed in non-PHB-synthesizing strains, showing that expression of this heterologous protein has an unexpected protective effect in E. coli, under both normal and stress conditions, resulting in increased growth and higher resistance to both heat shock and superoxide stress by paraquat. In addition, PhaP expression was shown to reduce RpoH protein levels during heat shock, probably by reducing or titrating the levels of misfolded proteins.

Small RNAs endow a transcriptional activator with essential repressor functions for single-tier control of a global stress regulon.

The Escherichia coli σ(E) envelope stress response monitors and repairs the outer membrane, a function central to the life of Gram-negative bacteria. The σ(E) stress response was characterized as a single-tier activation network comprised of ~100 genes, including the MicA and RybB noncoding sRNAs. These highly expressed sRNAs were thought to carry out the specialized function of halting de novo synthesis of several abundant porins when envelope homeostasis was perturbed. Using a systematic target profiling and validation approach we discovered that MicA and RybB are each global mRNA repressors of both distinct and shared targets, and that the two sRNAs constitute a posttranscriptional repression arm whose regulatory scope rivals that of the protein-based σ(E) activation arm. Intriguingly, porin mRNAs constitute only ~1/3 of all targets and new nonporin targets predict roles for MicA and RybB in crosstalk with other regulatory responses. This work also provides an example of evolutionarily unrelated sRNAs that are coinduced and bind the same targets, but at different sites. Our finding that expression of either MicA or RybB sRNA protects the cell from the loss of viability experienced when σ(E) activity is inadequate illustrates the importance of the posttranscriptional repression arm of the response. σ(E) is a paradigm of a single-tier stress response with a clear division of labor in which highly expressed noncoding RNAs (MicA, RybB) endow a transcriptional factor intrinsically restricted to gene activation (σ(E)) with the opposite repressor function.

Using DNA microarrays to assay part function.

In recent years, the capability of synthetic biology to design large genetic circuits has dramatically increased due to rapid advances in DNA synthesis technology and development of tools for large-scale assembly of DNA fragments. Large genetic circuits require more components (parts), especially regulators such as transcription factors, sigma factors, and viral RNA polymerases to provide increased regulatory capability, and also devices such as sensors, receivers, and signaling molecules. All these parts may have a potential impact upon the host that needs to be considered when designing and fabricating circuits. DNA microarrays are a well-established technique for global monitoring of gene expression and therefore are an ideal tool for systematically assessing the impact of expressing parts of genetic circuits in host cells. Knowledge of part impact on the host enables the user to design circuits from libraries of parts taking into account their potential impact and also to possibly modify the host to better tolerate stresses induced by the engineered circuit. In this chapter, we present the complete methodology of performing microarrays from choice of array platform, experimental design, preparing samples for array hybridization, and associated data analysis including preprocessing, normalization, clustering, identifying significantly differentially expressed genes, and interpreting the data based on known biology. With these methodologies, we also include lists of bioinformatic resources and tools for performing data analysis. The aim of this chapter is to provide the reader with the information necessary to be able to systematically catalog the impact of genetic parts on the host and also to optimize the operation of fully engineered genetic circuits.

Predicting strength and function for promoters of the Escherichia coli alternative sigma factor, sigmaE.

Sequenced bacterial genomes provide a wealth of information but little understanding of transcriptional regulatory circuits largely because accurate prediction of promoters is difficult. We examined two important issues for accurate promoter prediction: (1) the ability to predict promoter strength and (2) the sequence properties that distinguish between active and weak/inactive promoters. We addressed promoter prediction using natural core promoters recognized by the well-studied alternative sigma factor, Escherichia coli sigma(E), as a representative of group 4 sigmas, the largest sigma group. To evaluate the contribution of sequence to promoter strength and function, we used modular position weight matrix models comprised of each promoter motif and a penalty score for suboptimal motif location. We find that a combination of select modules is moderately predictive of promoter strength and that imposing minimal motif scores distinguished active from weak/inactive promoters. The combined -35/-10 score is the most important predictor of activity. Our models also identified key sequence features associated with active promoters. A conserved "AAC" motif in the -35 region is likely to be a general predictor of function for promoters recognized by group 4 sigmas. These results provide valuable insights into sequences that govern promoter strength, distinguish active and inactive promoters for the first time, and are applicable to both in vivo and in vitro measures of promoter strength.

Reduced capacity of alternative sigmas to melt promoters ensures stringent promoter recognition.

In bacteria, multiple sigmas direct RNA polymerase to distinct sets of promoters. Housekeeping sigmas direct transcription from thousands of promoters, whereas most alternative sigmas are more selective, recognizing more highly conserved promoter motifs. For sigma(32) and sigma(28), two Escherichia coli Group 3 sigmas, altering a few residues in Region 2.3, the portion of sigma implicated in promoter melting, to those universally conserved in housekeeping sigmas relaxed their stringent promoter requirements and significantly enhanced melting of suboptimal promoters. All Group 3 sigmas and the more divergent Group 4 sigmas have nonconserved amino acids at these positions and rarely transcribe >100 promoters. We suggest that the balance of "melting" and "recognition" functions of sigmas is critical to setting the stringency of promoter recognition. Divergent sigmas may generally use a nonoptimal Region 2.3 to increase promoter stringency, enabling them to mount a focused response to altered conditions.

Promoter strength properties of the complete sigma E regulon of Escherichia coli and Salmonella enterica.

The sigma(E)-directed envelope stress response maintains outer membrane homeostasis and is an important virulence determinant upon host infection in Escherichia coli and related bacteria. sigma(E) is activated by at least two distinct mechanisms: accumulation of outer membrane porin precursors and an increase in the alarmone ppGpp upon transition to stationary phase. Expression of the sigma(E) regulon is driven from a suite of approximately 60 sigma(E)-dependent promoters. Using green fluorescent protein fusions to each of these promoters, we dissected promoter contributions to the output of the regulon under a variety of in vivo conditions. We found that the sigma(E) promoters exhibit a large dynamic range, with a few strong and many weak promoters. Interestingly, the strongest promoters control either transcriptional regulators or functions related to porin homeostasis, the very functions conserved among E. coli and its close relatives. We found that (i) the strength of most promoters is significantly affected by the presence of the upstream (-35 to -65) region of the promoter, which encompasses the UP element, a binding site for the C-terminal domain of the alpha-subunit of RNA polymerase; (ii) ppGpp generally activates sigma(E) promoters, and (iii) sigma(E) promoters are responsive to changing sigma(E) holoenzyme levels under physiological conditions, reinforcing the idea that the sigma(E) regulon is extremely dynamic, enabling cellular adaptation to a constantly changing environment.

Mutational analysis of Escherichia coli sigma28 and its target promoters reveals recognition of a composite -10 region, comprised of an 'extended -10' motif and a core -10 element.

Sigma28 controls the expression of flagella-related genes and is the most widely distributed alternative sigma factor, present in motile Gram-positive and Gram-negative bacteria. The distinguishing feature of sigma28 promoters is a long -10 region (GCCGATAA). Despite the fact that the upstream GC is highly conserved, previous studies have not indicated a functional role for this motif. Here we examine the functional relevance of the GCCG motif and determine which residues in sigma28 participate in its recognition. We find that the GCCG motif is a functionally important composite element. The upstream GC constitutes an extended -10 motif and is recognized by R91, a residue in Domain 3 of sigma28. The downstream CG is the upstream edge of -10 region of the promoter; two residues in Region 2.4, D81 and R84, participate in its recognition. Consistent with their role in base-specific recognition of the promoter, R91, D81 and D84 are universally conserved in sigma28 orthologues. Sigma28 is the second Group 3 sigma shown to use an extended -10 region in promoter recognition, raising the possibility that other Group 3 sigmas will do so as well.

Dissection of recognition determinants of Escherichia coli sigma32 suggests a composite -10 region with an 'extended -10' motif and a core -10 element.

Sigma32 controls expression of heat shock genes in Escherichia coli and is widely distributed in proteobacteria. The distinguishing feature of sigma32 promoters is a long -10 region (CCCCATNT) whose tetra-C motif is important for promoter activity. Using alanine-scanning mutagenesis of sigma32 and in vivo and in vitro assays, we identified promoter recognition determinants of this motif. The most downstream C (-13) is part of the -10 motif; our work confirms and extends recognition determinants of -13C. Most importantly, our work suggests that the two upstream Cs (-16, -15) constitute an 'extended -10' recognition motif that is recognized by K130, a residue universally conserved in beta- and gamma-proteobacteria. This residue is located in the alpha-helix of sigmaDomain 3 that mediates recognition of the extended -10 promoter motif in other sigmas. K130 is not conserved in alpha- and delta-/epsilon-proteobacteria and we found that sigma32 from the alpha-proteobacterium Caulobacter crescentus does not need the extended -10 motif for high promoter activity. This result supports the idea that K130 mediates extended -10 recognition. Sigma32 is the first Group 3 sigma shown to use the 'extended -10' recognition motif.

Technical considerations in using DNA microarrays to define regulons.

Transcription is the major regulatory target of gene expression in bacteria, and is controlled by many regulatory proteins and RNAs. Microarrays are a powerful tool to study the regulation of transcription on a genomic scale. Here we describe the use of transcription profiling and ChIP-chip to study transcriptional regulation in bacteria. Transcription profiling determines the outcome of regulatory events whereas ChIP-chip identifies the protein-DNA interactions that determine these events. Together they can provide detailed information on transcriptional regulatory systems.

Convergence of molecular, modeling, and systems approaches for an understanding of the Escherichia coli heat shock response.

The heat shock response (HSR) is a homeostatic response that maintains the proper protein-folding environment in the cell. This response is universal, and many of its components are well conserved from bacteria to humans. In this review, we focus on the regulation of one of the most well-characterized HSRs, that of Escherichia coli. We show that even for this simple model organism, we still do not fully understand the central component of heat shock regulation, a chaperone-mediated negative feedback loop. In addition, we review other components that contribute to the regulation of the HSR in E. coli and discuss how these additional components contribute to regulation. Finally, we discuss recent genomic experiments that reveal additional functional aspects of the HSR.

SigmaE regulates and is regulated by a small RNA in Escherichia coli.

RybB is a small, Hfq-binding noncoding RNA originally identified in a screen of conserved intergenic regions in Escherichia coli. Fusions of the rybB promoter to lacZ were used to screen plasmid genomic libraries and genomic transposon mutants for regulators of rybB expression. A number of plasmids, including some carrying rybB, negatively regulated the fusion. An insertion in the rep helicase and one upstream of dnaK decreased expression of the fusion. Multicopy suppressors of these insertions led to identification of two plasmids that stimulated the fusion. One contained the gene for the response regulator OmpR; the second contained mipA, encoding a murein hydrolase. The involvement of MipA and OmpR in cell surface synthesis suggested that the rybB promoter might be dependent on sigma(E). The sequence upstream of the +1 of rybB contains a consensus sigma(E) promoter. The activity of rybB-lacZ was increased in cells lacking the RseA anti-sigma factor and when sigma(E) was overproduced from a heterologous promoter. The activity of rybB-lacZ and the detection of RybB were totally abolished in an rpoE-null strain. In vitro, sigma(E) efficiently transcribes from this promoter. Both a rybB mutation and an hfq mutation significantly increased expression of both rybB-lacZ and rpoE-lacZ fusions, consistent with negative regulation of the sigma(E) response by RybB and other small RNAs. Based on the plasmid screens, NsrR, a repressor sensitive to nitric oxide, was also found to negatively regulate sigma(E)-dependent promoters in an RseA-independent fashion.

Hfq modulates the sigmaE-mediated envelope stress response and the sigma32-mediated cytoplasmic stress response in Escherichia coli.

Hfq, a chaperone for small noncoding RNAs, regulates many processes in Escherichia coli, including the sigma(S)-mediated general stress response. Here we used microarray analysis to identify the changes in gene expression resulting from lack of Hfq. We identify several potential new targets for Hfq regulation, including genes encoding outer membrane proteins, enzymes, factors, and transporters. Many of these genes are involved in amino acid uptake and biosynthesis, sugar uptake and metabolism, and cell energetics. In addition, we find altered regulation of the sigma(E)- and sigma(32)-mediated stress responses, which we analyze further. We show that cells lacking Hfq induce the sigma(E)-mediated envelope stress response and are defective in sigma(E)-mediated repression of outer membrane proteins. We also show that the sigma(32)-mediated cytoplasmic stress response is repressed in cells lacking Hfq due to increased expression of DnaK. Furthermore, we show that cells lacking Hfq are defective in the "long-term adaptation" of sigma(32) to chronic chaperone overexpression. Together, our results indicate that Hfq may play a general role in stress response regulation in E. coli.

Regulon and promoter analysis of the E. coli heat-shock factor, sigma32, reveals a multifaceted cellular response to heat stress.

The heat-shock response (HSR), a universal cellular response to heat, is crucial for cellular adaptation. In Escherichia coli, the HSR is mediated by the alternative sigma factor, sigma32. To determine its role, we used genome-wide expression analysis and promoter validation to identify genes directly regulated by sigma32 and screened ORF overexpression libraries to identify sigma32 inducers. We triple the number of genes validated to be transcribed by sigma32 and provide new insights into the cellular role of this response. Our work indicates that the response is propagated as the regulon encodes numerous global transcriptional regulators, reveals that sigma70 holoenzyme initiates from 12% of sigma32 promoters, which has important implications for global transcriptional wiring, and identifies a new role for the response in protein homeostasis, that of protecting complex proteins. Finally, this study suggests that the response protects the cell membrane and responds to its status: Fully 25% of sigma32 regulon members reside in the membrane and alter its functionality; moreover, a disproportionate fraction of overexpressed proteins that induce the response are membrane localized. The intimate connection of the response to the membrane rationalizes why a major regulator of the response resides in that cellular compartment.

Conserved and variable functions of the sigmaE stress response in related genomes.

Bacteria often cope with environmental stress by inducing alternative sigma (sigma) factors, which direct RNA polymerase to specific promoters, thereby inducing a set of genes called a regulon to combat the stress. To understand the conserved and organism-specific functions of each sigma, it is necessary to be able to predict their promoters, so that their regulons can be followed across species. However, the variability of promoter sequences and motif spacing makes their prediction difficult. We developed and validated an accurate promoter prediction model for Escherichia coli sigmaE, which enabled us to predict a total of 89 unique sigmaE-controlled transcription units in E. coli K-12 and eight related genomes. SigmaE controls the envelope stress response in E. coli K-12. The portion of the regulon conserved across genomes is functionally coherent, ensuring the synthesis, assembly, and homeostasis of lipopolysaccharide and outer membrane porins, the key constituents of the outer membrane of Gram-negative bacteria. The larger variable portion is predicted to perform pathogenesis-associated functions, suggesting that sigmaE provides organism-specific functions necessary for optimal host interaction. The success of our promoter prediction model for sigmaE suggests that it will be applicable for the prediction of promoter elements for many alternative sigma factors.

Uses and pitfalls of microarrays for studying transcriptional regulation.

Microarrays provide a powerful new tool for understanding the regulation of gene expression in bacteria. Many recent publications have used microarrays for identifying regulon members and stimulons that describe the complex organismal responses to environmental perturbations. The use of bioinformatics to identify DNA binding sites of transcription factors greatly facilitates the interpretation of these experiments. Understanding the transcriptome of an organism includes identifying all transcripts and mapping their 5' and 3' ends. High-density oligonucleotide arrays have enabled the identification of many new transcripts, including small RNAs and antisense RNAs.

Architectural requirements for optimal activation by tandem CRP molecules at a class I CRP-dependent promoter.

The Escherichia coli cyclic AMP receptor protein (CRP) activates transcription at target promoters by interacting with the C-terminal domain of the RNA polymerase alpha subunit. We have constructed a set of promoters carrying tandem DNA sites for CRP with one site centred at position -61.5 and the other site located at different upstream positions. Optimal CRP-dependent activation of transcription is observed when the upstream DNA site for CRP is located at position -93.5 or at position -103.5. Evidence is presented to suggest that activation by the upstream-bound CRP molecule is due to interaction with the C-terminal domain of the RNA polymerase alpha subunit.