Aswin Sai Narain Seshasayee

Direct and indirect effects of H-NS and Fis on global gene expression control in Escherichia coli

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.

Complete Genome Sequence of Uropathogenic Proteus mirabilis, a Master of both Adherence and Motility

The gram-negative enteric bacterium Proteus mirabilis is a frequent cause of urinary tract infections in individuals with long-term indwelling catheters or with complicated urinary tracts (e.g., due to spinal cord injury or anatomic abnormality). P. mirabilis bacteriuria may lead to acute pyelonephritis, fever, and bacteremia. Most notoriously, this pathogen uses urease to catalyze the formation of kidney and bladder stones or to encrust or obstruct indwelling urinary catheters. Here we report the complete genome sequence of P. mirabilis HI4320, a representative strain cultured in our laboratory from the urine of a nursing home patient with a long-term (> or =30 days) indwelling urinary catheter. The genome is 4.063 Mb long and has a G+C content of 38.88%. There is a single plasmid consisting of 36,289 nucleotides. Annotation of the genome identified 3,685 coding sequences and seven rRNA loci. Analysis of the sequence confirmed the presence of previously identified virulence determinants, as well as a contiguous 54-kb flagellar regulon and 17 types of fimbriae. Genes encoding a potential type III secretion system were identified on a low-G+C-content genomic island containing 24 intact genes that appear to encode all components necessary to assemble a type III secretion system needle complex. In addition, the P. mirabilis HI4320 genome possesses four tandem copies of the zapE metalloprotease gene, genes encoding six putative autotransporters, an extension of the atf fimbrial operon to six genes, including an mrpJ homolog, and genes encoding at least five iron uptake mechanisms, two potential type IV secretion systems, and 16 two-component regulators.

Evidence of non-random mutation rates suggests an evolutionary risk management strategy

A central tenet in evolutionary theory is that mutations occur randomly with respect to their value to an organism; selection then governs whether they are fixed in a population. This principle has been challenged by long-standing theoretical models predicting that selection could modulate the rate of mutation itself. However, our understanding of how the mutation rate varies between different sites within a genome has been hindered by technical difficulties in measuring it. Here we present a study that overcomes previous limitations by combining phylogenetic and population genetic techniques. Upon comparing 34 Escherichia coli genomes, we observe that the neutral mutation rate varies by more than an order of magnitude across 2,659 genes, with mutational hot and cold spots spanning several kilobases. Importantly, the variation is not random: we detect a lower rate in highly expressed genes and in those undergoing stronger purifying selection. Our observations suggest that the mutation rate has been evolutionarily optimized to reduce the risk of deleterious mutations. Current knowledge of factors influencing the mutation rate—including transcription-coupled repair and context-dependent mutagenesis—do not explain these observations, indicating that additional mechanisms must be involved. The findings have important implications for our understanding of evolution and the control of mutations.

Comparative genomics of cyclic-di-GMP signalling in bacteria: post-translational regulation and catalytic activity

Cyclic-di-GMP is a bacterial second messenger that controls the switch between motile and sessile states. It is synthesized by proteins containing the enzymatic GGDEF domain and degraded by the EAL domain. Many bacterial genomes encode several copies of proteins containing these domains, raising questions on how the activities of parallel c-di-GMP signalling systems are segregated to avoid potentially deleterious cross-talk. Moreover, many ‘hybrid’ proteins contain both GGDEF and EAL domains; the relationship between the two apparently opposing enzymatic activities has been termed a ‘biochemical conundrum’. Here, we present a computational analysis of 11 248 GGDEF- and EAL-containing proteins in 867 prokaryotic genomes to address these two outstanding questions. Over half of these proteins contain a signal for cell-surface localization, and a majority accommodate a signal-sensing partner domain; these indicate widespread prevalence of post-translational regulation that may segregate the activities of proteins that are co-expressed. By examining the conservation of amino acid residues in the GGDEF and EAL catalytic sites, we show that there are predominantly two types of hybrid proteins. In the first, both sites are intact; an additional regulatory partner domain, present in most of these proteins, might determine the balance between the two enzymatic activities. In the second type, only the EAL catalytic site is intact; these—unlike EAL-only proteins—generally contain a signal-sensing partner domain, suggesting distinct modes of regulation for EAL activity under different sequence contexts. Finally, we discuss the role of proteins that have lost GGDEF and EAL catalytic sites as potential c-di-GMP-binding effectors. Our findings will serve as a genomic framework for interpreting ongoing molecular investigations of these proteins.

Genomic analysis of DNA binding and gene regulation by homologous nucleoid-associated proteins IHF and HU in Escherichia coli K12

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.

Genomics of DNA cytosine methylation in Escherichia coli reveals its role in stationary phase transcription

DNA cytosine methylation regulates gene expression in mammals. In bacteria, its role in gene expression and genome architecture is less understood. Here we perform high-throughput sequencing of bisulfite-treated genomic DNA from Escherichia coli K12 to describe, for the first time, the extent of cytosine methylation of bacterial DNA at single-base resolution. Whereas most target sites (C(m)CWGG) are fully methylated in stationary phase cells, many sites with an extended CC(m)CWGG motif are only partially methylated in exponentially growing cells. We speculate that these partially methylated sites may be selected, as these are slightly correlated with the risk of spontaneous, non-synonymous conversion of methylated cytosines to thymines. Microarray analysis in a cytosine methylation-deficient mutant of E. coli shows increased expression of the stress response sigma factor RpoS and many of its targets in stationary phase. Thus, DNA cytosine methylation is a regulator of stationary phase gene expression in E. coli.

Principles of transcriptional regulation and evolution of the metabolic system in E. coli

Organisms must adapt to make optimal use of the metabolic system in response to environmental changes. In the long-term, this involves evolution of the genomic repertoire of enzymes; in the short-term, transcriptional control ensures that appropriate enzymes are expressed in response to transitory extracellular conditions. Unicellular organisms are particularly susceptible to environmental changes; however, genome-scale impact of these modulatory effects has not been explored so far in bacteria. Here, we integrate genome-scale data to investigate the evolutionary trends and transcriptional control of metabolism in Escherichia coli K12. Globally, the regulatory system is organized in a clear hierarchy of general and specific transcription factors (TFs) that control differing ranges of metabolic functions. Further, catabolic, anabolic, and central metabolic pathways are targeted by distinct combinations of these TFs. Locally, enzymes catalyzing sequential reactions in a metabolic pathway are co-regulated by the same TFs. Regulation is more complex at junctions: General TFs control the overall activity of all connecting reactions, whereas specific TFs control individual enzymes. Divergent junctions play a special role in delineating metabolic pathways and decouple the regulation of incoming and outgoing reactions. We find little evidence for differential usage of isozymes, which are generally co-expressed in similar conditions, and thus are likely to reinforce the metabolic system through redundancy. Finally, we show that enzymes controlled by the same TFs have a strong tendency to co-evolve, suggesting a significant constraint to maintain similar regulatory regimes during evolution. Catabolic, anabolic, and central energy pathways evolve differently, emphasizing the role of the environment in shaping the metabolic system. Many of the observations also occur in yeast, and our findings may apply across large evolutionary distances.

Activation of Candida rugosa lipase at alkane–aqueous interfaces: A molecular dynamics study

The effect of solvent hydrophobicity on activation of Candida rugosa lipase (CRL) was investigated by performing molecular dynamics simulations for four nano seconds (ns). The closed/inactive conformer of CRL (PDB code 1TRH) was solvated in three alkane–aqueous environments. The alkanes aggregated in a predominantly aqueous environment and by 1 ns a stable spherical alkane–aqueous interface had formed. This led to the interfacial activation of CRL. On analyzing the simulated conformers with the closed conformer of CRL, the flap was found to have opened from a closed state by 7.7 Å, 10.2 Å, 13.1 Å at hexane–aqueous, octane–aqueous, and decane–aqueous interfaces. Further, essential dynamics analysis revealed that major anharmonic fluctuations were confined to residues 64–81, the flap of CRL.

Genome and Transcriptome Adaptation Accompanying Emergence of the Definitive Type 2 Host-Restricted Salmonella enterica Serovar Typhimurium Pathovar

ABSTRACT Salmonella enterica serovar Typhimurium definitive type 2 (DT2) is host restricted to Columba livia (rock or feral pigeon) but is also closely related to S. Typhimurium isolates that circulate in livestock and cause a zoonosis characterized by gastroenteritis in humans. DT2 isolates formed a distinct phylogenetic cluster within S. Typhimurium based on whole-genome-sequence polymorphisms. Comparative genome analysis of DT2 94-213 and S. Typhimurium SL1344, DT104, and D23580 identified few differences in gene content with the exception of variations within prophages. However, DT2 94-213 harbored 22 pseudogenes that were intact in other closely related S. Typhimurium strains. We report a novel in silico approach to identify single amino acid substitutions in proteins that have a high probability of a functional impact. One polymorphism identified using this method, a single-residue deletion in the Tar protein, abrogated chemotaxis to aspartate in vitro. DT2 94-213 also exhibited an altered transcriptional profile in response to culture at 42°C compared to that of SL1344. Such differentially regulated genes included a number involved in flagellum biosynthesis and motility. IMPORTANCE Whereas Salmonella enterica serovar Typhimurium can infect a wide range of animal species, some variants within this serovar exhibit a more limited host range and altered disease potential. Phylogenetic analysis based on whole-genome sequences can identify lineages associated with specific virulence traits, including host adaptation. This study represents one of the first to link pathogen-specific genetic signatures, including coding capacity, genome degradation, and transcriptional responses to host adaptation within a Salmonella serovar. We performed comparative genome analysis of reference and pigeon-adapted definitive type 2 (DT2) S. Typhimurium isolates alongside phenotypic and transcriptome analyses, to identify genetic signatures linked to host adaptation within the DT2 lineage. Whereas Salmonella enterica serovar Typhimurium can infect a wide range of animal species, some variants within this serovar exhibit a more limited host range and altered disease potential. Phylogenetic analysis based on whole-genome sequences can identify lineages associated with specific virulence traits, including host adaptation. This study represents one of the first to link pathogen-specific genetic signatures, including coding capacity, genome degradation, and transcriptional responses to host adaptation within a Salmonella serovar. We performed comparative genome analysis of reference and pigeon-adapted definitive type 2 (DT2) S. Typhimurium isolates alongside phenotypic and transcriptome analyses, to identify genetic signatures linked to host adaptation within the DT2 lineage.

Context-dependent conservation of DNA methyltransferases in bacteria

DNA methytransferases (MTs) in bacteria are best understood in the context of restriction–modification (R–M) systems, which act as bacterial immune systems against incoming DNA including phages, but have also been described as selfish elements. But several orphan MTs, which are not associated with any restriction enzyme, have also been characterized and may protect against parasitism by R–M systems. The occurrence of MTs in these two contexts, namely as part of R–M systems or as orphans, is poorly understood. Here we report the results of a comparative genomic survey of DNA MTs across ∼1000 bacterial genomes. We show that orphan MTs overwhelm R–M systems in their occurrence. In general, R–M MTs are poorly conserved, whereas orphans are nearly as conserved within a genus as any average gene. However, oligonucleotide usage and conservation patterns across genera suggest that both forms of MTs might have been horizontally acquired. We suggest that many orphan MTs might be ‘degradation’ products of R–M systems, based on the properties of orphan MTs encoded adjacent to highly diverged REs. In addition, several fully degraded R–M systems exist in which both the MT and the RE are highly divergent from their corresponding reference R–M pair. Despite their sporadic occurrence, conserved R–M systems are present in strength in two highly transformable genera, in which they may contribute to selection against integration of foreign DNA.