Bacterial transcriptional networks consist of hundreds of transcription factors and thousands of promoters. However, the true complexity of transcription in a bacterial pathogen and the effect of the environments encountered during infection remain to be established. We present a simplified approach for global promoter identification in bacteria using RNA-seq-based transcriptomic analyses of 22 distinct infection-relevant environmental conditions. Individual RNA samples were combined to identify most of the 3,838 Salmonella enterica serovar Typhimurium promoters in just two RNA-seq runs. Individual in vitro conditions stimulated characteristic transcriptional signatures, and the suite of 22 conditions induced transcription of 86% of all S. Typhimurium genes. We highlight the environmental conditions that induce the Salmonella pathogenicity islands and present a small RNA expression landscape of 280 sRNAs. This publicly available compendium of environmentally controlled expression of every transcriptional feature of S. Typhimurium constitutes a useful resource for the bacterial research community.

SummaryThe stringent response to amino acid starvation, whereby stable RNA synthesis is curtailed in favour of transcription of amino acid biosynthetic genes, is controlled by the alarmone ppGpp. To elucidate the extent of gene expression effected by ppGpp, we designed an experimental system based on starvation for isoleucine, which could be applied to both wild-type Escherichia coli and the multiauxotrophic relA spoT mutant (ppGpp 0 ). We used microarrays to profile the response to amino acid starvation in both strains. The wild-type response included induction of the general stress response, downregulation of genes involved in production of macromolecular structures and comprehensive restructuring of metabolic gene expression, but not induction of amino acid biosynthesis genes en masse. This restructuring of metabolism was confirmed using kinetic Biolog assays. These responses were profoundly altered in the ppGpp 0 strain. Furthermore, upon isoleucine starvation, the ppGpp 0 strain exhibited a larger cell size and continued growth, ultimately producing 50% more biomass than the wildtype, despite producing a similar amount of protein. This mutant phenotype correlated with aberrant gene expression in diverse processes, including DNA replication, cell division, and fatty acid and membrane biosynthesis. We present a model that expands and functionally integrates the ppGpp-mediated stringent response to include control of virtually all macromolecular synthesis and intermediary metabolism.

Whole-genome expression profiling revealed Escherichia coli MG1655 genes induced by growth on mucus, conditions designed to mimic nutrient availability in the mammalian intestine. Most were nutritional genes corresponding to catabolic pathways for nutrients found in mucus. We knocked out several pathways and tested the relative fitness of the mutants for colonization of the mouse intestine in competition with their wild-type parent. We found that only mutations in sugar pathways affected colonization, not phospholipid and amino acid catabolism, not gluconeogenesis, not the tricarboxylic acid cycle, and not the pentose phosphate pathway. Gluconate appeared to be a major carbon source used by E. coli MG1655 to colonize, having an impact on both the initiation and maintenance stages. N-acetylglucosamine and N-acetylneuraminic acid appeared to be involved in initiation, but not maintenance. Glucuronate, mannose, fucose, and ribose appeared to be involved in maintenance, but not initiation. The in vitro order of preference for these seven sugars paralleled the relative impact of the corresponding metabolic lesions on colonization: gluconate > N-acetylglucosamine > N-acetylneuraminic acid ‫؍‬ glucuronate > mannose > fucose > ribose. The results of this systematic analysis of nutrients used by E. coli MG1655 to colonize the mouse intestine are intriguing in light of the nutrientniche hypothesis, which states that the ecological niches within the intestine are defined by nutrient availability. Because humans are presumably colonized with different commensal strains, differences in nutrient availability may provide an open niche for infecting E. coli pathogens in some individuals and a barrier to infection in others.

The carbon sources that support the growth of pathogenic Escherichia coli O157:H7 in the mammalian intestine have not previously been investigated. In vivo, the pathogenic E. coli EDL933 grows primarily as single cells dispersed within the mucus layer that overlies the mouse cecal epithelium. We therefore compared the pathogenic strain and the commensal E. coli strain MG1655 modes of metabolism in vitro, using a mixture of the sugars known to be present in cecal mucus, and found that the two strains used the 13 sugars in a similar order and cometabolized as many as 9 sugars at a time. We conducted systematic mutation analyses of E. coli EDL933 and E. coli MG1655 by using lesions in the pathways used for catabolism of 13 mucus-derived sugars and five other compounds for which the corresponding bacterial gene system was induced in the transcriptome of cells grown on cecal mucus. Each of 18 catabolic mutants in both bacterial genetic backgrounds was fed to streptomycin-treated mice, together with the respective wild-type parent strain, and their colonization was monitored by fecal plate counts. None of the mutations corresponding to the five compounds not found in mucosal polysaccharides resulted in colonization defects. Based on the mutations that caused colonization defects, we determined that both E. coli EDL933 and E. coli MG1655 used arabinose, fucose, and N-acetylglucosamine in the intestine. In addition, E. coli EDL933 used galactose, hexuronates, mannose, and ribose, whereas E. coli MG1655 used gluconate and N-acetylneuraminic acid. The colonization defects of six catabolic lesions were found to be additive with E. coli EDL933 but not with E. coli MG1655. The data indicate that pathogenic E. coli EDL933 uses sugars that are not used by commensal E. coli MG1655 to colonize the mouse intestine. The results suggest a strategy whereby invading pathogens gain advantage by simultaneously consuming several sugars that may be available because they are not consumed by the commensal intestinal microbiota.

Fig. 2. Culture reproducibility of whole-genome expression profiles. The normalized, averaged data for the first time-point of the diauxic culture (culture 1) and the H 2 O 2 -treated culture (culture 2) are compared. The correlation coefficient is 0.992.

Escherichia coli MG1655 acid-inducible genes were identified by whole-genome expression profiling. Cultures were grown to the mid-logarithmic phase on acidified glucose minimal medium, conditions that induce glutamate-dependent acid resistance (AR), while the other AR systems are either repressed or not induced. A total of 28 genes were induced in at least two of three experiments in which the gene expression profiles of cells grown in acid (pH 5.5 or 4.5) were compared to those of cells grown at pH 7.4. As expected, the genes encoding glutamate decarboxylase, gadA and gadB, were significantly induced. Interestingly, two acid-inducible genes code for small basic proteins with pIs of >10.5, and six code for small acidic proteins with pIs ranging from 5.7 to 4.0; the roles of these small basic and acidic proteins in acid resistance are unknown. The acid-induced genes represented only five functional grouping categories, including eight genes involved in metabolism, nine associated with cell envelope structures or modifications, two encoding chaperones, six regulatory genes, and six unknown genes. It is unlikely that all of these genes are involved in the glutamate-dependent AR. However, nine acid-inducible genes are clustered in the gadA region, including hdeA, which encodes a putative periplasmic chaperone, and four putative regulatory genes. One of these putative regulators, yhiE, was shown to significantly increase acid resistance when overexpressed in cells that had not been preinduced by growth at pH 5.5, and mutation of yhiE decreased acid resistance; yhiE could therefore encode an activator of AR genes. Thus, the acid-inducible genes clustered in the gadA region appear to be involved in glutatmate-dependent acid resistance, although their specific roles remain to be elucidated.

DNA arrays of the entire set of Escherichia coli genes were used to measure the genomic expression patterns of cells growing in late logarithmic phase on minimal glucose medium and on Luria broth containing glucose. Ratios of the transcript levels for all 4,290E. coli protein-encoding genes (cds) were obtained, and analysis of the expression ratio data indicated that the physiological state of the cells under the two growth conditions could be ascertained. The cells in the rich medium grew faster, and expression of the majority of the translation apparatus genes was significantly elevated under this growth condition, consistent with known patterns of growth rate-dependent regulation and increased rate of protein synthesis in rapidly growing cells. The cells grown on minimal medium showed significantly elevated expression of many genes involved in biosynthesis of building blocks, most notably the amino acid biosynthetic pathways. Nearly half of the known RpoS-dependent genes were expressed at significantly higher levels in minimal medium than in rich medium, and rpoS expression was similarly elevated. The role of RpoS regulation in these logarithmic phase cells was suggested by the functions of the RpoS dependent genes that were induced. The hallmark features of E. coli cells growing on glucose minimal medium appeared to be the formation and excretion of acetate, metabolism of the acetate, and protection of the cells from acid stress. A hypothesis invoking RpoS and UspA (universal stress protein, also significantly elevated in minimal glucose medium) as playing a role in coordinating these various aspects and consequences of glucose and acetate metabolism was generated. This experiment demonstrates that genomic expression assays can be applied in a meaningful way to the study of whole-bacterial-cell physiology for the generation of hypotheses and as a guide for more detailed studies of particular genes of interest.

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.