Interactions between the 2.4 and 4.2 regions of sigmaS, the stress-specific sigma factor of Escherichia coli, and the -10 and -35 promoter elements

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.

Phylogenies from unaligned proteomes using sequence environments of amino acid residues

Alignment-free methods for sequence comparison and phylogeny inference have attracted a great deal of attention in recent years. Several algorithms have been implemented in diverse software packages. Despite the great number of existing methods, most of them are based on word statistics. Although they propose different filtering and weighting strategies and explore different metrics, their performance may be limited by the phylogenetic signal preserved in these words. Herein, we present a different approach based on the species-specific amino acid neighborhood preferences. These differential preferences can be assessed in the context of vector spaces. In this way, a distance-based method to build phylogenies has been developed and implemented into an easy-to-use R package. Tests run on real-world datasets show that this method can reconstruct phylogenetic relationships with high accuracy, and often outperforms other alignment-free approaches. Furthermore, we present evidence that the new method can perform reliably on datasets formed by non-orthologous protein sequences, that is, the method not only does not require the identification of orthologous proteins, but also does not require their presence in the analyzed dataset. These results suggest that the neighborhood preference of amino acids conveys a phylogenetic signal that may be of great utility in phylogenomics.

Universal functions of the σ finger in alternative σ factors during transcription initiation by bacterial RNA polymerase

The bacterial σ factor plays the central role in promoter recognition by RNA polymerase (RNAP). The primary σ factor, involved in transcription of housekeeping genes, was also shown to participate in the initiation of RNA synthesis and promoter escape by RNAP. In the open promoter complex, the σ finger formed by σ region 3.2 directly interacts with the template DNA strand upstream of the transcription start site. Here, we analysed the role of the σ finger in transcription initiation by four alternative σ factors in Escherichia coli, σ³⁸, σ³², σ²⁸ and σ²⁴. We found that deletions of the σ finger to various extent compromise the activity of RNAP holoenzymes containing alternative σ factors, especially at low NTP concentrations. All four σs are able to utilize NADH as a noncanonical priming substrate but it has only mild effects on the efficiency of transcription initiation. The mediators of the stringent response, transcription factor DksA and the alarmone ppGpp decrease RNAP activity and promoter complex stability for all four σ factors on tested promoters. For all σs except σ³⁸, deletions of the σ finger conversely increase the stability of promoter complexes and decrease their sensitivity to DksA and ppGpp. The result suggests that the σ finger plays a universal role in transcription initiation by alternative σ factors and sensitizes promoter complexes to the action of global transcription regulators DksA and ppGpp by modulating promoter complex stability.

Dissecting the molecular evolution of fluoroquinolone-resistant Shigella sonnei

Shigella sonnei increasingly dominates the international epidemiological landscape of shigellosis. Treatment options for S. sonnei are dwindling due to resistance to several key antimicrobials, including the fluoroquinolones. Here we analyse nearly 400 S. sonnei whole genome sequences from both endemic and non-endemic regions to delineate the evolutionary history of the recently emergent fluoroquinolone-resistant S. sonnei. We reaffirm that extant resistant organisms belong to a single clonal expansion event. Our results indicate that sequential accumulation of defining mutations (gyrA-S83L, parC-S80I, and gyrA-D87G) led to the emergence of the fluoroquinolone-resistant S. sonnei population around 2007 in South Asia. This clone was then transmitted globally, resulting in establishments in Southeast Asia and Europe. Mutation analysis suggests that the clone became dominant through enhanced adaptation to oxidative stress. Experimental evolution reveals that under fluoroquinolone exposure in vitro, resistant S. sonnei develops further intolerance to the antimicrobial while the susceptible counterpart fails to attain complete resistance.

Shigella sonnei is one of the main species causing shigellosis worldwide. Here the authors analyse nearly 400 S. sonnei genome sequences and carry out experimental evolution experiments to shed light into the evolutionary processes underlying the recent emergence of fluoroquinolone resistance in this pathogen.

sigmaS, a major player in the response to environmental stresses in Escherichia coli: role, regulation and mechanisms of promoter recognition

Bacterial cells often face hostile environmental conditions, to which they adapt by activation of stress responses. In Escherichia coli, environmental stresses resulting in significant reduction in growth rate stimulate the expression of the rpoS gene, encoding the alternative σ factor σ(S). The σ(S) protein associates with RNA polymerase, and through transcription of genes belonging to the rpoS regulon allows the activation of a 'general stress response', which protects the bacterial cell from harmful environmental conditions. Each step of this process is finely tuned in order to cater to the needs of the bacterial cell: in particular, selective promoter recognition by σ(S) is achieved through small deviations from a common consensus DNA sequence for both σ(S) and the housekeeping σ(70). Recognition of specific DNA elements by σ(S) is integrated with the effects of environmental signals and the interaction with regulatory proteins, in what represents a fascinating example of multifactorial regulation of gene expression. In this report, we discuss the function of the rpoS gene in the general stress response, and review the current knowledge on regulation of rpoS expression and on promoter recognition by σ(S).

The E. coli anti-sigma factor Rsd: studies on the specificity and regulation of its expression

Background

Among the seven different sigma factors in E. coli σ⁷⁰ has the highest concentration and affinity for the core RNA polymerase. The E. coli protein Rsd is regarded as an anti-sigma factor, inhibiting σ⁷⁰-dependent transcription at the onset of stationary growth. Although binding of Rsd to σ⁷⁰ has been shown and numerous structural studies on Rsd have been performed the detailed mechanism of action is still unknown.

Methodology/Principal Findings

We have performed studies to unravel the function and regulation of Rsd expression in vitro and in vivo. Cross-linking and affinity binding revealed that Rsd is able to interact with σ⁷⁰, with the core enzyme of RNA polymerase and is able to form dimers in solution. Unexpectedly, we find that Rsd does also interact with σ³⁸, the stationary phase-specific sigma factor. This interaction was further corroborated by gel retardation and footprinting studies with different promoter fragments and σ³⁸- or σ⁷⁰-containing RNA polymerase in presence of Rsd. Under competitive in vitro transcription conditions, in presence of both sigma factors, a selective inhibition of σ⁷⁰-dependent transcription was prevailing, however. Analysis of rsd expression revealed that the nucleoid-associated proteins H-NS and FIS, StpA and LRP bind to the regulatory region of the rsd promoters. Furthermore, the major promoter P2 was shown to be down-regulated in vivo by RpoS, the stationary phase-specific sigma factor and the transcription factor DksA, while induction of the stringent control enhanced rsd promoter activity. Most notably, the dam-dependent methylation of a cluster of GATC sites turned out to be important for efficient rsd transcription.

Conclusions/Significance

The results contribute to a better understanding of the intricate mechanism of Rsd-mediated sigma factor specificity changes during stationary phase.

In vitro transcription profiling of the σS subunit of bacterial RNA polymerase: re-definition of the σS regulon and identification of σS-specific promoter sequence elements

Specific promoter recognition by bacterial RNA polymerase is mediated by σ subunits, which assemble with RNA polymerase core enzyme (E) during transcription initiation. However, σ⁷⁰ (the housekeeping σ subunit) and σ^S (an alternative σ subunit mostly active during slow growth) recognize almost identical promoter sequences, thus raising the question of how promoter selectivity is achieved in the bacterial cell. To identify novel sequence determinants for selective promoter recognition, we performed run-off/microarray (ROMA) experiments with RNA polymerase saturated either with σ⁷⁰ (Eσ⁷⁰) or with σ^S (Eσ^S) using the whole Escherichia coli genome as DNA template. We found that Eσ⁷⁰, in the absence of any additional transcription factor, preferentially transcribes genes associated with fast growth (e.g. ribosomal operons). In contrast, Eσ^S efficiently transcribes genes involved in stress responses, secondary metabolism as well as RNAs from intergenic regions with yet-unknown function. Promoter sequence comparison suggests that, in addition to different conservation of the −35 sequence and of the UP element, selective promoter recognition by either form of RNA polymerase can be affected by the A/T content in the −10/+1 region. Indeed, site-directed mutagenesis experiments confirmed that an A/T bias in the −10/+1 region could improve promoter recognition by Eσ^S.

Mobility of plasmids

Plasmids are key vectors of horizontal gene transfer and essential genetic engineering tools. They code for genes involved in many aspects of microbial biology, including detoxication, virulence, ecological interactions, and antibiotic resistance. While many studies have decorticated the mechanisms of mobility in model plasmids, the identification and characterization of plasmid mobility from genome data are unexplored. By reviewing the available data and literature, we established a computational protocol to identify and classify conjugation and mobilization genetic modules in 1,730 plasmids. This allowed the accurate classification of proteobacterial conjugative or mobilizable systems in a combination of four mating pair formation and six relaxase families. The available evidence suggests that half of the plasmids are nonmobilizable and that half of the remaining plasmids are conjugative. Some conjugative systems are much more abundant than others and preferably associated with some clades or plasmid sizes. Most very large plasmids are nonmobilizable, with evidence of ongoing domestication into secondary chromosomes. The evolution of conjugation elements shows ancient divergence between mobility systems, with relaxases and type IV coupling proteins (T4CPs) often following separate paths from type IV secretion systems. Phylogenetic patterns of mobility proteins are consistent with the phylogeny of the host prokaryotes, suggesting that plasmid mobility is in general circumscribed within large clades. Our survey suggests the existence of unsuspected new relaxases in archaea and new conjugation systems in cyanobacteria and actinobacteria. Few genes, e.g., T4CPs, relaxases, and VirB4, are at the core of plasmid conjugation, and together with accessory genes, they have evolved into specific systems adapted to specific physiological and ecological contexts.

Metagenomic analysis of the viral communities in fermented foods

Viruses are recognized as the most abundant biological components on Earth, and they regulate the structure of microbial communities in many environments. In soil and marine environments, microorganism-infecting phages are the most common type of virus. Although several types of bacteriophage have been isolated from fermented foods, little is known about the overall viral assemblages (viromes) of these environments. In this study, metagenomic analyses were performed on the uncultivated viral communities from three fermented foods, fermented shrimp, kimchi, and sauerkraut. Using a high-throughput pyrosequencing technique, a total of 81,831, 70,591 and 69,464 viral sequences were obtained from fermented shrimp, kimchi and sauerkraut, respectively. Moreover, 37 to 50% of these sequences showed no significant hit against sequences in public databases. There were some discrepancies between the prediction of bacteriophages hosts via homology comparison and bacterial distribution, as determined from 16S rRNA gene sequencing. These discrepancies likely reflect the fact that the viral genomes of fermented foods are poorly represented in public databases. Double-stranded DNA viral communities were amplified from fermented foods by using a linker-amplified shotgun library. These communities were dominated by bacteriophages belonging to the viral order Caudovirales (i.e., Myoviridae, Podoviridae, and Siphoviridae). This study indicates that fermented foods contain less complex viral communities than many other environmental habitats, such as seawater, human feces, marine sediment, and soil.

Genome organization and characterization of the virulent lactococcal phage 1358 and its similarities to Listeria phages

Virulent phage 1358 is the reference member of a rare group of phages infecting Lactococcus lactis. Electron microscopy revealed a typical icosahedral capsid connected to one of the smallest noncontractile tails found in a lactococcal phage of the Siphoviridae family. Microbiological characterization identified a burst size of 72 virions released per infected host cell and a latent period of 90 min. The host range of phage 1358 was limited to 3 out of the 60 lactococcal strains tested. Moreover, this phage was insensitive to four Abi systems (AbiK, AbiQ, AbiT, and AbiV). The genome of phage 1358 consisted of a linear, double-stranded, 36,892-bp DNA molecule containing 43 open reading frames (ORFs). At least 14 ORFs coded for structural proteins, as identified by SDS-PAGE coupled to liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses. The genomic organization was similar to those of other siphophages. All genes were on the same coding strand and in the same orientation. This lactococcal phage was unique, however, in its 51.4% GC content, much higher than those of other phages infecting this low-GC Gram-positive host. A bias for GC-rich codons was also observed. Comparative analyses showed that several phage 1358 structural proteins shared similarity with two Listeria monocytogenes phages, P35 and P40. The possible origin and evolution of lactococcal phage 1358 is discussed.

Poising of Escherichia coli RNA polymerase and its release from the sigma 38 C-terminal tail for osmY transcription

Bacteria must adapt their transcription to overcome the osmotic stress associated with the gastrointestinal tract of their host. This requires the sigma 38 (rpoS) form of RNA polymerase. Here, chromatin immunoprecipitation experiments show that activation is associated with a poise-and-release mechanism in vivo. A C-terminal tail unique among sigma factors is shown to be required for in vivo recruitment of RNA polymerase to the promoter region prior to osmotic shock. C-terminal domain tail-dependent transcription in vivo can be mimicked by using the intracellular signaling molecule potassium glutamate in vitro. Following signaling, the barrier to elongation into the gene body is overcome and RNA polymerase is released to produce osmY mRNA.

Computational prediction of RpoS and RpoD regulatory sites in Geobacter sulfurreducens using sequence and gene expression information

RpoS, the sigma S subunit of RNA polymerase, is vital during the growth and survival of Geobacter sulfurreducens under conditions typically encountered in its native subsurface environments. We investigated the conservation of sites that may be important for RpoS function in G. sulfurreducens. We also employed sequence information and expression microarray data to predict G. sulfurreducens genome sites that may be related to RpoS regulation. Hierarchical clustering identified three clusters of significantly downregulated genes in the rpoS deletion mutant. The search for conserved overrepresented motifs in co-regulated operons identified likely -35 and -10 promoter elements upstream of a number of functionally important G. sulfurreducens operons that were downregulated in the rpoS deletion mutant. Putative -35/-10 promoter elements were also identified in the G. sulfurreducens genome using sequence similarity searches to matrices of -35/-10 promoter elements found in G. sulfurreducens and in Escherichia coli. Due to a sufficient degree of sequence similarity between -35/-10 promoter elements for RpoS, RpoD, and other sigma factors, both the sequence similarity searches and the search for conserved overrepresented motifs using microarray data may identify promoter elements for both RpoS and other sigma factors.

Osmolyte-induced transcription: -35 region elements and recognition by sigma38 (rpoS)

In order to meet osmotic challenges in the gastrointestinal tract, enteric bacteria rapidly accumulate salts of glutamate and other weak organic acids. The ensuing transcriptional activation is mediated by unknown elements at sigma38 (rpoS)-dependent promoters. Here we identify DNA elements needed for high levels of transcription in the presence of salt and acetate and show that they are associated with the -35 regions of target promoters. Unrelated -35 region sequences are shown to specify maximal salt-challenged transcription at the otsB promoter and maximal acetate-challenged transcription at the cfa promoter. Mutants in sigma38 are isolated that contribute to bypassing the salt response and most of these cluster in a small segment corresponding to the presumptive -35 DNA recognition determinant of the protein. Overall, the data suggest that an ensemble of -35 region elements exists at sigma38 promoters and these can help mediate responsiveness to physiological challenges through interactions involving region 4 of the sigma38 protein.

Role of the spacer between the -35 and -10 regions in sigmas promoter selectivity in Escherichia coli

In vitro, the sigma(s) subunit of RNA polymerase (RNAP), RpoS, recognizes nearly identical -35 and -10 promoter consensus sequences as the vegetative sigma70. In vivo, promoter selectivity of RNAP holoenzyme containing either sigma(s) (Esigma(s)) or sigma70 (Esigma70) seems to be achieved by the differential ability of the two holoenzymes to tolerate deviations from the promoter consensus sequence. In this study, we suggest that many natural sigma(s)-dependent promoters possess a -35 element, a feature that has been considered as not conserved among sigma(s)-dependent promoters. These -35 hexamers are mostly non-optimally spaced from the -10 region, but nevertheless functional. A +/- 2 bp deviation from the optimal spacer length of 17 bp or the complete absence of a -35 consensus sequence decreases overall promoter activity, but at the same time favours Esigma(s) in its competition with Esigma70 for promoter recognition. On the other hand, the reduction of promoter activity due to shifting of the -35 element can be counterbalanced by an activity-stimulating feature such as A/T-richness of the spacer region without compromising Esigma(s) selectivity. Based on mutational analysis of sigma(s), we suggest a role of regions 2.5 and 4 of sigma(s) in sensing sub-optimally located -35 elements.

Substitutions in Region 2.4 of σ70 Allow Recognition of the σS-Dependent aidB Promoter

The strict dependence of transcription from the aidB promoter (PaidB) on the Esigma(S) form of RNA polymerase is because of the presence of a C nucleotide as the first residue of the -10 promoter sequence (-12C), which does not allow an open complex formation by Esigma(70). In this report, sigma(70) mutants carrying either the Q437H or the T440I single amino acid substitutions, which allow -12C recognition by sigma(70), were tested for their ability to carry out transcription from PaidB. The Gln-437 and Thr-440 residues are located in region 2.4 of sigma(70) and correspond to Gln-152 and Glu-155 in sigma(S). Interestingly, the Q437H mutant of sigma(70), but not T440I, was able to promote an open complex formation and to initiate transcription at PaidB. In contrast to T440I, a T440E mutant was proficient in carrying out transcription from PaidB. No sigma(70) mutant displayed significantly increased interaction with a PaidB mutant in which the -12C was substituted by a T (PaidB((C12T))), which is also efficiently recognized by wild type sigma(70). The effect of the T440E mutation suggests that the corresponding Glu-155 residue in sigma(S) might be involved in -12C recognition. However, substitution to alanine of the Glu-155 residue, as well as of Gln-152, in the sigma(S) protein did not significantly affect Esigma(S) interaction with PaidB. Our results reiterate the importance of the -12C residue for sigma(S)-specific promoter recognition and strongly suggest that interaction with the -10 sequence and open complex formation are carried out by different determinants in the two sigma factors.

SigmaS-dependent gene expression at the onset of stationary phase in Escherichia coli: function of sigmaS-dependent genes and identification of their promoter sequences

The sigma(S) subunit of RNA polymerase, the product of the rpoS gene, controls the expression of genes responding to starvation and cellular stresses. Using gene array technology, we investigated rpoS-dependent expression at the onset of stationary phase in Escherichia coli grown in rich medium. Forty-one genes were expressed at significantly lower levels in an rpoS mutant derived from the MG1655 strain; for 10 of these, we also confirmed rpoS and stationary-phase dependence by reverse transcription-PCR. Only seven genes (dps, osmE, osmY, sodC, rpsV, wrbA, and yahO) had previously been recognized as rpoS dependent. Several newly identified rpoS-dependent genes are involved in the uptake and metabolism of amino acids, sugars, and iron. Indeed, the rpoS mutant strain shows severely impaired growth on some sugars such as fructose and N-acetylglucosamine. The rpoS gene controls the production of indole, which acts as a signal molecule in stationary-phase cells, via regulation of the tnaA-encoded tryptophanase enzyme. Genes involved in protein biosynthesis, encoding the ribosome-associated protein RpsV (sra) and the initiation factor IF-1 (infA), were also induced in an rpoS-dependent fashion. Using primer extension, we determined the promoter sequences of a selection of rpoS-regulated genes representative of different functional classes. Significant fractions of these promoters carry sequence features specific for Esigma(S) recognition of the -10 region, such as cytosines at positions -13 (70%) and -12 (30%) as well as a TG motif located upstream of the -10 region (50%), thus supporting the TGN(0-2)C(C/T)ATA(C/A)T consensus sequence recently proposed for sigma(S).