IHF and HU are two heterodimeric nucleoid-associated proteins (NAP) that belong to the same protein family but interact differently with the DNA. IHF is a sequence-specific DNA-binding protein that bends the DNA by over 160°. HU is the most conserved NAP, which binds non-specifically to duplex DNA with a particular preference for targeting nicked and bent DNA. Despite their importance, the in vivo interactions of the two proteins to the DNA remain to be described at a high resolution and on a genome-wide scale. Further, the effects of these proteins on gene expression on a global scale remain contentious. Finally, the contrast between the functions of the homo- and heterodimeric forms of proteins deserves the attention of further study. Here we present a genome-scale study of HU- and IHF binding to the Escherichia coli K12 chromosome using ChIP-seq. We also perform microarray analysis of gene expression in single- and double-deletion mutants of each protein to identify their regulons. The sequence-specific binding profile of IHF encompasses ∼30% of all operons, though the expression of <10% of these is affected by its deletion suggesting combinatorial control or a molecular backup. The binding profile for HU is reflective of relatively non-specific binding to the chromosome, however, with a preference for A/T-rich DNA. The HU regulon comprises highly conserved genes including those that are essential and possibly supercoiling sensitive. Finally, by performing ChIP-seq experiments, where possible, of each subunit of IHF and HU in the absence of the other subunit, we define genome-wide maps of DNA binding of the proteins in their hetero- and homodimeric forms.

Conflicts between replication and transcription challenge chromosome duplication. Escherichia coli replisome movement along transcribed DNA is promoted by Rep and UvrD accessory helicases with Δrep ΔuvrD cells being inviable under rapid growth conditions. We have discovered that mutations in a tRNA gene, aspT, in an aminoacyl tRNA synthetase, AspRS, and in a translation factor needed for efficient proline-proline bond formation, EF-P, suppress Δrep ΔuvrD lethality. Thus replication-transcription conflicts can be alleviated by the partial sacrifice of a mechanism that reduces replicative barriers, namely translating ribosomes that reduce RNA polymerase backtracking. Suppression depends on RelA-directed synthesis of (p)ppGpp, a signalling molecule that reduces replication-transcription conflicts, with RelA activation requiring ribosomal pausing. Levels of (p)ppGpp in these suppressors also correlate inversely with the need for Rho activity, an RNA translocase that can bind to emerging transcripts and displace transcription complexes. These data illustrate the fine balance between different mechanisms in facilitating gene expression and genome duplication and demonstrate that accessory helicases are a major determinant of this balance. This balance is also critical for other aspects of bacterial survival: the mutations identified here increase persistence indicating that similar mutations could arise in naturally occurring bacterial populations facing antibiotic challenge.

The goal of this group project has been to coordinate and bring up-to-date information on all genes of Escherichia coli K-12. Annotation of the genome of an organism entails identification of genes, the boundaries of genes in terms of precise start and end sites, and description of the gene products. Known and predicted functions were assigned to each gene product on the basis of experimental evidence or sequence analysis. Since both kinds of evidence are constantly expanding, no annotation is complete at any moment in time. This is a snapshot analysis based on the most recent genome sequences of two E.coli K-12 bacteria. An accurate and up-to-date description of E.coli K-12 genes is of particular importance to the scientific community because experimentally determined properties of its gene products provide fundamental information for annotation of innumerable genes of other organisms. Availability of the complete genome sequence of two K-12 strains allows comparison of their genotypes and mutant status of alleles.

Small non-coding RNAs (sRNAs) are an emerging class of regulators of bacterial gene expression. Most of the regulatory Escherichia coli sRNAs known to date modulate translation of trans-encoded target mRNAs. We studied the specificity of sRNA target interactions using gene fusions to green fluorescent protein (GFP) as a novel reporter of translational control by bacterial sRNAs in vivo. Target sequences were selected from both monocistronic and polycistronic mRNAs. Upon expression of the cognate sRNA (DsrA, GcvB, MicA, MicC, MicF, RprA, RyhB, SgrS and Spot42), we observed highly specific translation repression/activation of target fusions under various growth conditions. Target regulation was also tested in mutants that lacked Hfq or RNase III, or which expressed a truncated RNase E (rne701). We found that translational regulation by these sRNAs was largely independent of full-length RNase E, e.g. despite the fact that ompA fusion mRNA decay could no longer be promoted by MicA. This is the first study in which multiple well-defined E.coli sRNA target pairs have been studied in a uniform manner in vivo. We expect our GFP fusion approach to be applicable to sRNA targets of other bacteria, and also demonstrate that Vibrio RyhB sRNA represses a Vibrio sodB fusion when co-expressed in E.coli.

In this study, we examined the intracellular whereabouts of Mrr, a cryptic type IV restriction endonuclease of Escherichia coli K12, in response to different conditions. In absence of stimuli triggering its activity, Mrr was found to be strongly associated with the nucleoid as a number of discrete foci, suggesting the presence of Mrr hotspots on the chromosome. Previously established elicitors of Mrr activity, such as exposure to high (hydrostatic) pressure (HP) or expression of the HhaII methyltransferase, both caused nucleoid condensation and an unexpected coalescence of Mrr foci. However, although the resulting Mrr/nucleoid complex was stable when triggered with HhaII, it tended to be only short-lived when elicited with HP. Moreover, HP-mediated activation of Mrr typically led to cellular blebbing, suggesting a link between chromosome and cellular integrity. Interestingly, Mrr variants could be isolated that were specifically compromised in either HhaII- or HP-dependent activation, underscoring a mechanistic difference in the way both triggers activate Mrr. In general, our results reveal that Mrr can take part in complex spatial distributions on the nucleoid and can be engaged in distinct modes of activity.

The gene expression state of exponentially growing Escherichia coli cells is manifested by high expression of essential and growth-associated genes and low levels of stress-related and horizontally acquired genes. An important player in maintaining this homeostasis is the H-NS-StpA gene silencing system. A Δhns-stpA deletion mutant results in high expression of otherwise-silent horizontally acquired genes, many located in the terminus-half of the chromosome, and an indirect downregulation of many highly expressed genes. The Δhns-stpA double mutant displays slow growth. Using laboratory evolution we address the evolutionary strategies that E. coli would adopt to redress this gene expression imbalance. We show that two global gene regulatory mutations-(i) point mutations inactivating the stress-responsive sigma factor RpoS or σ38 and (ii) an amplification of ∼40% of the chromosome centred around the origin of replication-converge in partially reversing the global gene expression imbalance caused by Δhns-stpA. Transcriptome data of these mutants further show a three-way link amongst the global gene regulatory networks of H-NS and σ38, as well as chromosome architecture. Increasing gene expression around the terminus of replication results in a decrease in the expression of genes around the origin and vice versa; this appears to be a persistent phenomenon observed as an association across ∼300 publicly-available gene expression data sets for E. coli. These global suppressor effects are transient and rapidly give way to more specific mutations, whose roles in reversing the growth defect of H-NS mutations remain to be understood.

Escherichia coli's cAMP receptor protein (CRP), the archetypal bacterial transcription factor, regulates over a hundred promoters by binding 22 bp symmetrical sites with the consensus core half-site TGTGA. However, Haemophilus influenzae has two types of CRP sites, one like E.coli's and one with the core sequence TGCGA that regulates genes required for DNA uptake (natural competence). Only the latter 'CRP-S' sites require both CRP and the coregulator Sxy for activation. To our knowledge, the TGTGA and TGCGA motifs are the first example of one transcription factor having two distinct binding-site motifs. Here we show that CRP-S promoters are widespread in the gamma-proteobacteria and demonstrate their Sxy-dependence in E.coli. Orthologs of most H.influenzae CRP-S-regulated genes are ubiquitous in the five best-studied gamma-proteobacteria families, Enterobacteriaceae, Pasteurellaceae, Pseudomonadaceae, Vibrionaceae and Xanthomonadaceae. Phylogenetic footprinting identified CRP-S sites in the promoter regions of the Enterobacteriaceae, Pasteurellaceae and Vibrionaceae orthologs, and canonical CRP sites in orthologs of genes known to be Sxy-independent in H.influenzae. Bandshift experiments confirmed that E.coli CRP-S sequences are low affinity binding sites for CRP, and mRNA analysis showed that they require CRP, cAMP (CRP's allosteric effector) and Sxy for gene induction. This work suggests not only that the gamma-proteobacteria share a common DNA uptake mechanism, but also that, in the three best studied families, their competence regulons share both CRP-S specificity and Sxy dependence.

Nucleoid-associated proteins (NAPs) are global regulators of gene expression in Escherichia coli, which affect DNA conformation by bending, wrapping and bridging the DNA. Two of these--H-NS and Fis--bind to specific DNA sequences and structures. Because of their importance to global gene expression, the binding of these NAPs to the DNA was previously investigated on a genome-wide scale using ChIP-chip. However, variation in their binding profiles across the growth phase and the genome-scale nature of their impact on gene expression remain poorly understood. Here, we present a genome-scale investigation of H-NS and Fis binding to the E. coli chromosome using chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-seq). By performing our experiments under multiple time-points during growth in rich media, we show that the binding regions of the two proteins are mutually exclusive under our experimental conditions. H-NS binds to significantly longer tracts of DNA than Fis, consistent with the linear spread of H-NS binding from high- to surrounding lower-affinity sites; the length of binding regions is associated with the degree of transcriptional repression imposed by H-NS. For Fis, a majority of binding events do not lead to differential expression of the proximal gene; however, it has a significant indirect effect on gene expression partly through its effects on the expression of other transcription factors. We propose that direct transcriptional regulation by Fis is associated with the interaction of tandem arrays of Fis molecules to the DNA and possible DNA bending, particularly at operon-upstream regions. Our study serves as a proof-of-principle for the use of ChIP-seq for global DNA-binding proteins in bacteria, which should become significantly more economical and feasible with the development of multiplexing techniques.

We devised software tools to systematically investigate the contents and contexts of bacterial tRNA and tmRNA genes, which are known insertion hotspots for genomic islands (GIs). The strategy, based on MAUVE-facilitated multigenome comparisons, was used to examine 87 Escherichia coli MG1655 tRNA and tmRNA genes and their orthologues in E.coli EDL933, E.coli CFT073 and Shigella flexneri Sf301. Our approach identified 49 GIs occupying approximately 1.7 Mb that mapped to 18 tRNA genes, missing 2 but identifying a further 30 GIs as compared with Islander [Y. Mantri and K. P. Williams (2004), Nucleic Acids Res., 32, D55-D58]. All these GIs had many strain-specific CDS, anomalous GC contents and/or significant dinucleotide biases, consistent with foreign origins. Our analysis demonstrated marked conservation of sequences flanking both empty tRNA sites and tRNA-associated GIs across all four genomes. Remarkably, there were only 2 upstream and 5 downstream deletions adjacent to the 328 loci investigated. In silico PCR analysis based on conserved flanking regions was also used to interrogate hotspots in another eight completely or partially sequenced E.coli and Shigella genomes. The tools developed are ideal for the analysis of other bacterial species and will lead to in silico and experimental discovery of new genomic islands.

Using a combination of mathematical modelling, statistical simulation and large-scale data analysis we study the properties of linear regulatory chains (LRCs) within gene regulatory networks (GRNs). Our modelling indicates that downstream genes embedded within LRCs are highly insulated from the variation in expression of upstream genes, and thus LRCs act as attenuators. This observation implies a progressively weaker functionality of LRCs as their length increases. When analyzing the preponderance of LRCs in the GRNs of Escherichia coli K12 and several other organisms, we find that very long LRCs are essentially absent. In both E. coli and M. tuberculosis we find that four-gene LRCs are intimately linked to identical feedback loops that are involved in potentially chaotic stress response, indicating that the dynamics of these potentially destabilising motifs are strongly restrained under homeostatic conditions. The same relationship is observed in a human cancer cell line (K562), and we postulate that four-gene LRCs act as 'universal attenuators'. These findings suggest a role for long LRCs in dampening variation in gene expression, thereby protecting cell identity, and in controlling dramatic shifts in cell-wide gene expression through inhibiting chaos-generating motifs.

The genetic code can be manipulated to reassign codons for the incorporation of non-standard amino acids (NSAA). Deletion of release factor 1 in Escherichia coli enhances translation of UAG (Stop) codons, yet may also extended protein synthesis at natural UAG terminated messenger RNAs. The fidelity of protein synthesis at reassigned UAG codons and the purity of the NSAA containing proteins produced require careful examination. Proteomics would be an ideal tool for these tasks, but conventional proteomic analyses cannot readily identify the extended proteins and accurately discover multiple amino acid (AA) insertions at a single UAG. To address these challenges, we created a new proteomic workflow that enabled the detection of UAG readthrough in native proteins in E. coli strains in which UAG was reassigned to encode phosphoserine. The method also enabled quantitation of NSAA and natural AA incorporation at UAG in a recombinant reporter protein. As a proof-of-principle, we measured the fidelity and purity of the phosphoserine orthogonal translation system (OTS) and used this information to improve its performance. Our results show a surprising diversity of natural AAs at reassigned stop codons. Our method can be used to improve OTSs and to quantify amino acid purity at reassigned codons in organisms with expanded genetic codes.

Sequencing small quantities of DNA is important for applications ranging from the assembly of uncultivable microbial genomes to the identification of cancer-associated mutations. To obtain sufficient quantities of DNA for sequencing, the small amount of starting material must be amplified significantly. However, existing methods often yield errors or non-uniform coverage, reducing sequencing data quality. Here, we describe digital droplet multiple displacement amplification, a method that enables massive amplification of low-input material while maintaining sequence accuracy and uniformity. The low-input material is compartmentalized as single molecules in millions of picoliter droplets. Because the molecules are isolated in compartments, they amplify to saturation without competing for resources; this yields uniform representation of all sequences in the final product and, in turn, enhances the quality of the sequence data. We demonstrate the ability to uniformly amplify the genomes of single Escherichia coli cells, comprising just 4.7 fg of starting DNA, and obtain sequencing coverage distributions that rival that of unamplified material. Digital droplet multiple displacement amplification provides a simple and effective method for amplifying minute amounts of DNA for accurate and uniform sequencing.

The rational design of synthetic regulatory circuits critically hinges on the availability of orthogonal and well-characterized building blocks. Here, we focus on extracytoplasmic function (ECF) σ factors, which are the largest group of alternative σ factors and hold extensive potential as synthetic orthogonal regulators. By assembling multiple ECF σ factors into regulatory cascades of varying length, we benchmark the scalability of the approach, showing that these 'autonomous timer circuits' feature a tuneable time delay between inducer addition and target gene activation. The implementation of similar timers in Escherichia coli and Bacillus subtilis shows strikingly convergent circuit behavior, which can be rationalized by a computational model. These findings not only reveal ECF σ factors as powerful building blocks for a rational, multi-layered circuit design, but also suggest that ECF σ factors are universally applicable as orthogonal regulators in a variety of bacterial species.

RNase E is an endoribonuclease found in many bacteria, including important human pathogens. Within Escherichia coli, it has been shown to have a major role in both the maturation of all classes of RNA involved in translation and the initiation of mRNA degradation. Thus, knowledge of the major determinants of RNase E cleavage is central to our understanding and manipulation of bacterial gene expression. We show here that the binding of RNase E to structured RNA elements is crucial for the processing of tRNA, can activate catalysis and may be important in mRNA degradation. The recognition of structured elements by RNase E is mediated by a recently discovered groove that is distant from the domains associated with catalysis. The functioning of this groove is shown here to be essential for E. coli cell viability and may represent a key point of evolutionary divergence from the paralogous RNase G family, which we show lack amino acid residues conserved within the RNA-binding groove of members of the RNase E family. Overall, this work provides new insights into the recognition and cleavage of RNA by RNase E and provides further understanding of the basis of RNase E essentiality in E. coli.

CRISPR/Cas9 is a promising tool in prokaryotic genome engineering, but its success is limited by the widely varying on-target activity of single guide RNAs (sgRNAs). Based on the association of CRISPR/Cas9-induced DNA cleavage with cellular lethality, we systematically profiled sgRNA activity by co-expressing a genome-scale library (∼70 000 sgRNAs) with Cas9 or its specificity-improved mutant in Escherichia coli. Based on this large-scale dataset, we constructed a comprehensive and high-density sgRNA activity map, which enables selecting highly active sgRNAs for any locus across the genome in this model organism. We also identified 'resistant' genomic loci with respect to CRISPR/Cas9 activity, notwithstanding the highly accessible DNA in bacterial cells. Moreover, we found that previous sgRNA activity prediction models that were trained on mammalian cell datasets were inadequate when coping with our results, highlighting the key limitations and biases of previous models. We hence developed an integrated algorithm to accurately predict highly effective sgRNAs, aiming to facilitate CRISPR/Cas9-based genome engineering, screenings and antimicrobials design in bacteria. We also isolated the important sgRNA features that contribute to DNA cleavage and characterized their key differences among wild type Cas9 and its mutant, shedding light on the biophysical mechanisms of the CRISPR/Cas9 system.

The Escherichia coli RecQ DNA helicase participates in a pathway of DNA repair that operates in parallel to the recombination pathway driven by the multisubunit helicase-nuclease machine RecBCD. The model mycobacterium Mycobacterium smegmatis executes homologous recombination in the absence of its helicase-nuclease machine AdnAB, though it lacks a homolog of E. coli RecQ. Here, we identify and characterize M. smegmatis RqlH, a RecQ-like helicase with a distinctive domain structure. The 691-amino acid RqlH polypeptide consists of a RecQ-like ATPase domain (amino acids 1-346) and tetracysteine zinc-binding domain (amino acids 435-499), separated by an RqlH-specific linker. RqlH lacks the C-terminal HRDC domain found in E. coli RecQ. Rather, the RqlH C-domain resembles bacterial ComF proteins and includes a phosphoribosyltransferase-like module. We show that RqlH is a DNA-dependent ATPase/dATPase that translocates 3'-5' on single-stranded DNA and has 3'-5' helicase activity. These functions inhere to RqlH-(1-505), a monomeric motor unit comprising the ATPase, linker and zinc-binding domains. RqlH homologs are distributed widely among bacterial taxa. The mycobacteria that encode RqlH lack a classical RecQ, though many other Actinobacteria have both RqlH and RecQ. Whereas E. coli K12 encodes RecQ but lacks a homolog of RqlH, other strains of E. coli have both RqlH and RecQ.

With the advent of high-throughput sequencing technologies, the rapid generation and accumulation of large amounts of sequencing data pose an insurmountable demand for efficient algorithms for constructing whole-genome phylogenies. The existing phylogenomic methods all use assembled sequences, which are often not available owing to the difficulty of assembling short-reads; this obstructs phylogenetic investigations on species without a reference genome. In this report, we present co-phylog, an assembly-free phylogenomic approach that creates a 'micro-alignment' at each 'object' in the sequence using the 'context' of the object and calculates pairwise distances before reconstructing the phylogenetic tree based on those distances. We explored the parameters' usages and the optimal working range of co-phylog, assessed co-phylog using the simulated next-generation sequencing (NGS) data and the real NGS raw data. We also compared co-phylog method with traditional alignment and alignment-free methods and illustrated the advantages and limitations of co-phylog method. In conclusion, we demonstrated that co-phylog is efficient algorithm and that it delivers high resolution and accurate phylogenies using whole-genome unassembled sequencing data, especially in the case of closely related organisms, thereby significantly alleviating the computational burden in the genomic era.

We present a computational method for the prediction of functional modules encoded in microbial genomes. In this work, we have also developed a formal measure to quantify the degree of consistency between the predicted and the known modules, and have carried out statistical significance analysis of consistency measures. We first evaluate the functional relationship between two genes from three different perspectives--phylogenetic profile analysis, gene neighborhood analysis and Gene Ontology assignments. We then combine the three different sources of information in the framework of Bayesian inference, and we use the combined information to measure the strength of gene functional relationship. Finally, we apply a threshold-based method to predict functional modules. By applying this method to Escherichia coli K12, we have predicted 185 functional modules. Our predictions are highly consistent with the previously known functional modules in E.coli. The application results have demonstrated that our approach is highly promising for the prediction of functional modules encoded in a microbial genome.

Targeted PCR amplification and high-throughput sequencing (amplicon sequencing) of 16S rRNA gene fragments is widely used to profile microbial communities. New long-read sequencing technologies can sequence the entire 16S rRNA gene, but higher error rates have limited their attractiveness when accuracy is important. Here we present a high-throughput amplicon sequencing methodology based on PacBio circular consensus sequencing and the DADA2 sample inference method that measures the full-length 16S rRNA gene with single-nucleotide resolution and a near-zero error rate. In two artificial communities of known composition, our method recovered the full complement of full-length 16S sequence variants from expected community members without residual errors. The measured abundances of intra-genomic sequence variants were in the integral ratios expected from the genuine allelic variants within a genome. The full-length 16S gene sequences recovered by our approach allowed Escherichia coli strains to be correctly classified to the O157:H7 and K12 sub-species clades. In human fecal samples, our method showed strong technical replication and was able to recover the full complement of 16S rRNA alleles in several E. coli strains. There are likely many applications beyond microbial profiling for which high-throughput amplicon sequencing of complete genes with single-nucleotide resolution will be of use.

In a wide range of genomes, it was observed that the usage of synonymous codons is biased toward specific codons and codon patterns. Factors that are implicated in the selection for codon usage include facilitation of fast and accurate translation. There are two types of translational errors: missense errors and processivity errors. There is considerable evidence in support of the hypothesis that codon usage is optimized to minimize missense errors. In contrast, little is known about the relationship between codon usage and frameshifting errors, an important form of processivity errors, which appear to occur at frequencies comparable to the frequencies of missense errors. Based on the recently proposed pause-and-slip model of frameshifting, we developed Frameshifting Robustness Score (FRS). We used this measure to test if the pattern of codon usage indicates optimization against frameshifting errors. We found that the FRS values of protein-coding sequences from four analyzed genomes (the bacteria Bacillus subtilis and Escherichia coli, and the yeasts Saccharomyces cerevisiae and Schizosaccharomyce pombe) were typically higher than expected by chance. Other properties of FRS patterns observed in B. subtilis, S. cerevisiae and S. pombe, such as the tendency of FRS to increase from the 5'- to 3'-end of protein-coding sequences, were also consistent with the hypothesis of optimization against frameshifting errors in translation. For E. coli, the results of different tests were less consistent, suggestive of a much weaker optimization, if any. Collectively, the results fit the concept of selection against mistranslation-induced protein misfolding being one of the factors shaping the evolution of both coding and non-coding sequences.