Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria.

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

Early fate of exogenous promoters in E. coli.

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Using stochastic cell division and death to probe minimal units of cellular replication.

The invariant cell initiation mass measured in bacterial growth experiments has been interpreted as a minimal unit of cellular replication. Here we argue that the existence of such minimal units induces a coupling between the rates of stochastic cell division and death. To probe this coupling we tracked live and dead cells in Escherichia coli populations treated with a ribosome-targeting antibiotic. We find that the growth exponent from macroscopic cell growth or decay measurements can be represented as the difference of microscopic first-order cell division and death rates. The boundary between cell growth and decay, at which the number of live cells remains constant over time, occurs at the minimal inhibitory concentration (MIC) of the antibiotic. This state appears macroscopically static but is microscopically dynamic: division and death rates exactly cancel at MIC but each is remarkably high, reaching 60% of the antibiotic-free division rate. A stochastic model of cells as collections of minimal replicating units we term 'widgets' reproduces both steady-state and transient features of our experiments. Sub-cellular fluctuations of widget numbers stochastically drive each new daughter cell to one of two alternate fates, division or death. First-order division or death rates emerge as eigenvalues of a stationary Markov process, and can be expressed in terms of the widget's molecular properties. High division and death rates at MIC arise due to low mean and high relative fluctuations of widget number. Isolating cells at the threshold of irreversible death might allow molecular characterization of this minimal replication unit.

The genome-scale interplay amongst xenogene silencing, stress response and chromosome architecture in Escherichia coli.

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.

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.

Cysteine desulfurase (IscS)-mediated fine-tuning of bioenergetics and SUF expression prevents Mycobacterium tuberculosis hypervirulence.

Iron-sulfur (Fe-S) biogenesis requires multiprotein assembly systems, SUF and ISC, in most prokaryotes. M. tuberculosis (Mtb) encodes a complete SUF system, the depletion of which was bactericidal. The ISC operon is truncated to a single gene iscS (cysteine desulfurase), whose function remains uncertain. Here, we show that MtbΔiscS is bioenergetically deficient and hypersensitive to oxidative stress, antibiotics, and hypoxia. MtbΔiscS resisted killing by nitric oxide (NO). RNA sequencing indicates that IscS is important for expressing regulons of DosR and Fe-S-containing transcription factors, WhiB3 and SufR. Unlike wild-type Mtb, MtbΔiscS could not enter a stable persistent state, continued replicating in mice, and showed hypervirulence. The suf operon was overexpressed in MtbΔiscS during infection in a NO-dependent manner. Suppressing suf expression in MtbΔiscS either by CRISPR interference or upon infection in inducible NO-deficient mice arrests hypervirulence. Together, Mtb redesigned the ISC system to "fine-tune" the expression of SUF machinery for establishing persistence without causing detrimental disease in the host.

Draft genome sequence of Staphylococcus aureus ST672, an emerging disease clone from India.

We report the draft genome sequence of methicillin-resistant Staphylococcus aureus (MRSA) strain ST672, an emerging disease clone in India, from a septicemia patient. The genome size is about 2.82 Mb with 2,485 open reading frames (ORFs). The staphylococcal cassette chromosome mec (SCCmec) element (type V) and immune evasion cluster appear to be different from those of strain ST772 on preliminary examination.

Draft genome sequence of Staphylococcus aureus 118 (ST772), a major disease clone from India.

We report the draft genome sequence of an ST772 Staphylococcus aureus disease isolate carrying staphylococcal cassette chromosome mec (SCCmec) type V from a pyomyositis patient. Our de novo short read assembly is ∼2.8 Mb and encodes a unique Panton-Valentine leukocidin (PVL) phage with structural genes similar to those of ϕ7247PVL and novel lysogenic genes at the N termini.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes.

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes.

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models.

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, Emax, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with Emax but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and Emax values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Emergence of networks of shared restriction-modification systems in phage-bacteria ecosystems.

Restriction-modification (RM) systems are the most ubiquitous bacterial defence systems against bacteriophages. Using genome sequence data, we showed that RM systems are often shared among bacterial strains in a structured way. Examining the network of interconnections between bacterial strains within genera, we found that many strains share more RM systems than expected compared with a suitable null model. We also found that many genera have a larger than expected number of bacterial strains with unique RM systems. We used population dynamics models of closed and open phage-bacteria ecosystems to qualitatively understand the selection pressures that could lead to such network structures with enhanced overlap or uniqueness. In our models, we found that the phages impose a selection pressure that favours bacteria with greater number of RM systems, and higher overlap of RM systems with other strains, but in bacteria-dominated states, this is opposed by the increased cost-to-growth rate of these bacteria. Similar to what we observed in the genome data, we found that two distinct bacterial strategies emerge - strains either have a greater overlap than expected, or, at the other extreme, have unique RM systems. The former strategy appears to dominate when the repertoire of available RM systems is smaller but the average number of RM systems per strain is larger.

Evolutionary and Comparative Analysis of Bacterial Nonhomologous End Joining Repair.

DNA double-strand breaks (DSBs) are a threat to genome stability. In all domains of life, DSBs are faithfully fixed via homologous recombination. Recombination requires the presence of an uncut copy of duplex DNA which is used as a template for repair. Alternatively, in the absence of a template, cells utilize error-prone nonhomologous end joining (NHEJ). Although ubiquitously found in eukaryotes, NHEJ is not universally present in bacteria. It is unclear as to why many prokaryotes lack this pathway. Toward understanding what could have led to the current distribution of bacterial NHEJ, we carried out comparative genomics and phylogenetic analysis across ∼6,000 genomes. Our results show that this pathway is sporadically distributed across the phylogeny. Ancestral reconstruction further suggests that NHEJ was absent in the eubacterial ancestor and can be acquired via specific routes. Integrating NHEJ occurrence data for archaea, we also find evidence for extensive horizontal exchange of NHEJ genes between the two kingdoms as well as across bacterial clades. The pattern of occurrence in bacteria is consistent with correlated evolution of NHEJ with key genome characteristics of genome size and growth rate; NHEJ presence is associated with large genome sizes and/or slow growth rates, with the former being the dominant correlate. Given the central role these traits play in determining the ability to carry out recombination, it is possible that the evolutionary history of bacterial NHEJ may have been shaped by requirement for efficient DSB repair.

Laboratory Evolution Experiments Help Identify a Predominant Region of Constitutive Stable DNA Replication Initiation.

The bacterium Escherichia coli can initiate replication in the absence of the replication initiator protein DnaA and/or the canonical origin of replication oriC in a ΔrnhA background. This phenomenon, which can be primed by R-loops, is called constitutive stable DNA replication (cSDR). Whether DNA replication during cSDR initiates in a stochastic manner through the length of the chromosome or at specific sites and how E. coli can find adaptations to loss of fitness caused by cSDR remain inadequately answered. We use laboratory evolution experiments of ΔrnhA-ΔdnaA strains followed by deep sequencing to show that DNA replication preferentially initiates within a broad region located ∼0.4 to 0.7 Mb clockwise of oriC. This region includes many bisulfite-sensitive sites, which have been previously defined as R-loop-forming regions, and includes a site containing sequence motifs that favor R-loop formation. Initiation from this region would result in head-on replication-transcription conflicts at rRNA loci. Inversions of these rRNA loci, which can partly resolve these conflicts, help the bacterium suppress the fitness defects of cSDR. These inversions partially restore the gene expression changes brought about by cSDR. The inversion, however, increases the possibility of conflicts at essential mRNA genes, which would utilize only a minuscule fraction of RNA polymerase molecules, most of which transcribe rRNA genes. Whether subsequent adaptive strategies would attempt to resolve these conflicts remains an open question.IMPORTANCE The bacterium E. coli can replicate its DNA even in the absence of the molecules that are required for canonical replication initiation. This often requires the formation of RNA-DNA hybrid structures and is referred to as constitutive stable DNA replication (cSDR). Where on the chromosome does cSDR initiate? We answer this question using laboratory evolution experiments and genomics and show that selection favors cSDR initiation predominantly at a region ∼0.6 Mb clockwise of oriC. Initiation from this site will result in more head-on collisions of DNA polymerase with RNA polymerase operating on rRNA loci. The bacterium adapts to this problem by inverting a region of the genome including several rRNA loci such that head-on collisions between the two polymerases are minimized. Understanding such evolutionary strategies in the context of cSDR can provide insights into the potential causes of resistance to antibiotics that target initiation of DNA replication.

Whole-Genome Sequencing of Sphingobium sp. Strain RSMS, a Highly Efficient Tributyl Phosphate-Degrading Bacterium.

Sphingobium sp. strain RSMS was described earlier as an efficient degrader of tributyl phosphate, an organic pollutant. This report describes the generation and annotation of the genome sequence of Sphingobium sp. strain RSMS, which will facilitate future studies to identify genetic elements responsible for the degradation of tributyl phosphate.

Targeting redox heterogeneity to counteract drug tolerance in replicating Mycobacterium tuberculosis.

The capacity of Mycobacterium tuberculosis (Mtb) to tolerate multiple antibiotics represents a major problem in tuberculosis (TB) management. Heterogeneity in Mtb populations is one of the factors that drives antibiotic tolerance during infection. However, the mechanisms underpinning this variation in bacterial population remain poorly understood. Here, we show that phagosomal acidification alters the redox physiology of Mtb to generate a population of replicating bacteria that display drug tolerance during infection. RNA sequencing of this redox-altered population revealed the involvement of iron-sulfur (Fe-S) cluster biogenesis, hydrogen sulfide (H2S) gas, and drug efflux pumps in antibiotic tolerance. The fraction of the pH- and redox-dependent tolerant population increased when Mtb infected macrophages with actively replicating HIV-1, suggesting that redox heterogeneity could contribute to high rates of TB therapy failure during HIV-TB coinfection. Pharmacological inhibition of phagosomal acidification by the antimalarial drug chloroquine (CQ) eradicated drug-tolerant Mtb, ameliorated lung pathology, and reduced postchemotherapeutic relapse in in vivo models. The pharmacological profile of CQ (C max and AUClast) exhibited no major drug-drug interaction when coadministered with first line anti-TB drugs in mice. Our data establish a link between phagosomal pH, redox metabolism, and drug tolerance in replicating Mtb and suggest repositioning of CQ to shorten TB therapy and achieve a relapse-free cure.

Regulation of Global Transcription in Escherichia coli by Rsd and 6S RNA.

In Escherichia coli, the sigma factor σ70 directs RNA polymerase to transcribe growth-related genes, while σ38 directs transcription of stress response genes during stationary phase. Two molecules hypothesized to regulate RNA polymerase are the protein Rsd, which binds to σ70, and the non-coding 6S RNA which binds to the RNA polymerase-σ70 holoenzyme. Despite multiple studies, the functions of Rsd and 6S RNA remain controversial. Here we use RNA-Seq in five phases of growth to elucidate their function on a genome-wide scale. We show that Rsd and 6S RNA facilitate σ38 activity throughout bacterial growth, while 6S RNA also regulates widely different genes depending upon growth phase. We discover novel interactions between 6S RNA and Rsd and show widespread expression changes in a strain lacking both regulators. Finally, we present a mathematical model of transcription which highlights the crosstalk between Rsd and 6S RNA as a crucial factor in controlling sigma factor competition and global gene expression.

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 A, 10.2 A, 13.1 A 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.

Nucleoid-Associated Proteins: Genome Level Occupancy and Expression Analysis.

The advent of Chromatin Immunoprecipitation sequencing (ChIP-Seq) has allowed the identification of genomic regions bound by a DNA binding protein in-vivo on a genome-wide scale. The impact of the DNA binding protein on gene expression can be addressed using transcriptome experiments in appropriate genetic settings. Overlaying the above two sources of data enables us to dissect the direct and indirect effects of a DNA binding protein on gene expression. Application of these techniques to Nucleoid Associated Proteins (NAPs) and Global Transcription Factors (GTFs) has underscored the complex relationship between DNA-protein interactions and gene expression change, highlighting the role of combinatorial control. Here, we demonstrate the usage of ChIP-Seq to infer binding properties and transcriptional effects of NAPs such as Fis and HNS, and the GTF CRP in the model organism Escherichia coli K12 MG1655 (E. coli).

Modulation of Global Transcriptional Regulatory Networks as a Strategy for Increasing Kanamycin Resistance of the Translational Elongation Factor-G Mutants in Escherichia coli.

Evolve and resequence experiments have provided us a tool to understand bacterial adaptation to antibiotics. In our previous work, we used short-term evolution to isolate mutants resistant to the ribosome targeting antibiotic kanamycin, and reported that Escherichia coli develops low cost resistance to kanamycin via different point mutations in the translation Elongation Factor-G (EF-G). Furthermore, we had shown that the resistance of EF-G mutants could be increased by second site mutations in the genes rpoD/cpxA/topA/cyaA Mutations in three of these genes had been discovered in earlier screens for aminoglycoside resistance. In this work, we expand our understanding of these second site mutations, the goal being to understand how these mutations affect the activities of the mutated gene products to confer resistance. We show that the mutation in cpxA most likely results in an active Cpx stress response. Further evolution of an EF-G mutant in a higher concentration of kanamycin than what was used in our previous experiments identified the cpxA locus as a primary target for a significant increase in resistance. The mutation in cyaA results in a loss of catalytic activity and probably results in resistance via altered CRP function. Despite a reduction in cAMP levels, the CyaAN600Y mutant has a transcriptome indicative of increased CRP activity, pointing to an unknown role for CyaA and / or cAMP in gene expression. From the transcriptomes of double and single mutants, we describe the epistasis between the mutation in EF-G and these second site mutations. We show that the large scale transcriptomic changes in the topoisomerase I (FusAA608E-TopAS180L) mutant likely result from increased negative supercoiling in the cell. Finally, genes with known roles in aminoglycoside resistance were present among the misregulated genes in the mutants.