CRP (cAMP receptor protein), the global regulator of genes for carbon source utilization in the absence of glucose, is the best-studied prokaryotic transcription factor. A total of 195 target promoters on the Escherichia coli genome have been proposed to be under the control of cAMP-bound CRP. Using the newly developed Genomic SELEX screening system of transcription factor-binding sequences, however, we have identified a total of at least 254 CRP-binding sites. Based on their location on the E. coli genome, we predict a total of at least 183 novel regulation target operons, altogether with the 195 hitherto known targets, reaching to the minimum of 378 promoters as the regulation targets of cAMP-CRP. All the promoters selected from the newly identified targets and examined by using the lacZ reporter assay were found to be under the control of CRP, indicating that the Genomic SELEX screening allowed to identify the CRP targets with high accuracy. Based on the functions of novel target genes, we conclude that CRP plays a key regulatory role in the whole processes from the selective transport of carbon sources, the glycolysis-gluconeogenesis switching to the metabolisms downstream of glycolysis, including tricarboxylic acid (TCA) cycle, pyruvate dehydrogenase (PDH) pathway and aerobic respiration. One unique regulation mode is that a single and the same CRP molecule bound within intergenic regions often regulates both of divergently transcribed operons.

We analyzed the transcriptome of Escherichia coli K-12 by strand-specific RNA sequencing at single-nucleotide resolution during steady-state (logarithmic-phase) growth and upon entry into stationary phase in glucose minimal medium. To generate high-resolution transcriptome maps, we developed an organizational schema which showed that in practice only three features are required to define operon architecture: the promoter, terminator, and deep RNA sequence read coverage. We precisely annotated 2,122 promoters and 1,774 terminators, defining 1,510 operons with an average of 1.98 genes per operon. Our analyses revealed an unprecedented view of E. coli operon architecture. A large proportion (36%) of operons are complex with internal promoters or terminators that generate multiple transcription units. For 43% of operons, we observed differential expression of polycistronic genes, despite being in the same operons, indicating that E. coli operon architecture allows fine-tuning of gene expression. We found that 276 of 370 convergent operons terminate inefficiently, generating complementary 3′ transcript ends which overlap on average by 286 nucleotides, and 136 of 388 divergent operons have promoters arranged such that their 5′ ends overlap on average by 168 nucleotides. We found 89 antisense transcripts of 397-nucleotide average length, 7 unannotated transcripts within intergenic regions, and 18 sense transcripts that completely overlap operons on the opposite strand. Of 519 overlapping transcripts, 75% correspond to sequences that are highly conserved in E. coli (>50 genomes). Our data extend recent studies showing unexpected transcriptome complexity in several bacteria and suggest that antisense RNA regulation is widespread. IMPORTANCE We precisely mapped the 5′ and 3′ ends of RNA transcripts across the E. coli K-12 genome by using a single-nucleotide analytical approach. Our resulting high-resolution transcriptome maps show that ca. one-third of E. coli operons are complex, with internal promoters and terminators generating multiple transcription units and allowing differential gene expression within these operons. We discovered extensive antisense transcription that results from more than 500 operons, which fully overlap or extensively overlap adjacent divergent or convergent operons. The genomic regions corresponding to these antisense transcripts are highly conserved in E. coli (including Shigella species), although it remains to be proven whether or not they are functional. Our observations of features unearthed by single-nucleotide transcriptome mapping suggest that deeper layers of transcriptional regulation in bacteria are likely to be revealed in the future.

Bacterial genomes are transcribed by DNA-dependent RNA polymerase (RNAP), which achieves gene selectivity through interaction with sigma factors that recognize promoters, and transcription factors (TFs) that control the activity and specificity of RNAP holoenzyme. To understand the molecular mechanisms of transcriptional regulation, the identification of regulatory targets is needed for all these factors. We then performed genomic SELEX screenings of targets under the control of each sigma factor and each TF. Here we describe the assembly of 156 SELEX patterns of a total of 116 TFs performed in the presence and absence of effector ligands. The results reveal several novel concepts: (i) each TF regulates more targets than hitherto recognized; (ii) each promoter is regulated by more TFs than hitherto recognized; and (iii) the binding sites of some TFs are located within operons and even inside open reading frames. The binding sites of a set of global regulators, including cAMP receptor protein, LeuO and Lrp, overlap with those of the silencer H-NS, suggesting that certain global regulators play an anti-silencing role. To facilitate sharing of these accumulated SELEX datasets with the research community, we compiled a database, ‘Transcription Profile of Escherichia coli’ (www.shigen.nig.ac.jp/ecoli/tec/).

Cra (or FruR), a global transcription factor with both repression and activation activities, controls a large number of the genes for glycolysis and gluconeogenesis. To get insights into the entire network of transcription regulation of the E. coli genome by Cra, we isolated a set of Cra‐binding sequences using an improved method of genomic SELEX. From the DNA sequences of 97 independently isolated DNA fragments by SELEX, the Cra‐binding sequences were identified in a total of ten regions on the E. coli genome, including promoters of six known genes and four hitherto‐unidentified genes. All six known promoters are repressed by Cra, but none of the activation‐type promoters were cloned after two cyles of SELEX, because the Cra‐binding affinity to the repression‐type promoters is higher than the activation‐type promoters, as determined by the quantitative gel shift assay. Of a total of four newly identified Cra‐binding sequences, two are associated with promoter regions of the gapA (glyceraldehyde 3‐phosphate dehydrogenase) and eno (enolase) genes, both involved in sugar metabolism. The regulation of newly identified genes by Cra was confirmed by the in vivo promoter strength assay using a newly developed TFP (two‐fluorescent protein) vector for promoter assay or by in vitro transcription assay in the presence of Cra protein.

The Escherichia coli RutR protein is the master regulator of genes involved in pyrimidine catabolism. Here we have used chromatin immunoprecipitation in combination with DNA microarrays to measure the binding of RutR across the chromosome of exponentially growing E. coli cells. Twenty RutR-binding targets were identified and analysis of these targets generated a DNA consensus logo for RutR binding. Complementary in vitro binding assays showed high-affinity RutR binding to 16 of the 20 targets, with the four low-affinity RutR targets lacking predicted key binding determinants. Surprisingly, most of the DNA targets for RutR are located within coding segments of the genome and appear to have little or no effect on transcript levels in the conditions tested. This contrasts sharply with other E. coli transcription factors whose binding sites are primarily located in intergenic regions. We suggest that either RutR has yet undiscovered function or that evolution has been slow to eliminate non-functional DNA sites for RutR because they do not have an adverse effect on cell fitness.

The promoter selectivity of Escherichia coli RNA polymerase is determined by the sigma subunit with promoter recognition activity. The model prokaryote Escherichia coli contains seven species of the sigma subunit, each recognizing a specific set of promoters. The major sigma subunit, sigma-70 encoded by rpoD, plays a major role in transcription of growth-related genes. Concomitant with the increase in detection of promoters functioning in vivo under various stressful conditions, the variation is expanding in the consensus sequence of RpoD promoters. In order to identify the canonical sequence of “constitutive promoters” that are recognized by the RNA polymerase holoenzyme containing RpoD sigma in the absence of supporting transcription factors, an in vitro mixed transcription assay was carried out using a whole set of variant promoters, each harboring one base replacement, within the model promoter with the conserved -35 and -10 sequences of RpoD promoters. The consensus sequences, TTGACA(-35) and TATAAT(-10), were identified to be ideal for the maximum level of open complex formation and the highest rate of promoter opening, respectively. For identification of the full range of constitutive promoters on the E. coli genome, a total of 2,701 RpoD holoenzyme-binding sites were identified by Genomic SELEX screening, and using the reconfirmed consensus promoter sequence, a total of maximum 669 constitutive promoters were identified, implying that the majority of hitherto identified promoters represents the TF-dependent “inducible promoters”. One unique feature of the constitutive promoters is the high level of promoter sequence conservation, about 85% carrying five-out-of-six agreements with -35 or -10 consensus sequence. The list of constitutive promoters provides the community resource toward estimation of the inducible promoters that operate under various stressful conditions in nature.