ABSTRACTBacteria have evolved a wide range of defense strategies to protect themselves against bacterial viruses (phages). However, the known mechanisms almost exclusively target phages with DNA genomes. While several bacterial toxin-antitoxin systems have been considered to cleave single-stranded bacterial RNA in response to stressful conditions, their role in protecting bacteria against phages with single-stranded RNA genomes has not been studied. Here we investigate the role of a representative toxin-antitoxin system, MazEF, in protectingEscherichia coliagainst two RNA phages – MS2 and Qβ. Our population-level experiments revealed that amazEFdeletion strain is more susceptible to RNA phage infection than the wild-type. At the single-cell level, deletion of themazEFlocus significantly shortened the time to lysis of individual bacteria challenged with RNA phage. At the genomic level, we found that the adenine-cytosine-adenine sequence, directly recognized and cleaved by the MazF toxin, is systematically underrepresented in the genomes of RNA phages that are known to infectE. coli, indicating selection for decreased probability of cleavage. These results suggest that in addition to other physiological roles, RNA-degrading toxin-antitoxin modules can function as a primitive immune system against RNA phages.

Fisher suggested advantageous genes would spread through populations as a wave so we sought genetic waves in evolving populations, as follows. By fusing a fluorescent marker to a drug efflux protein (AcrB) whose expression providesEscherichia coliwith resistance to some antibiotics, we quantified the evolution and spread of drug-resistantE. colithrough spacetime using image analysis and quantitative PCR. As is done in hospitals routinely, we exposed the bacterium to a gradient of antibiotic in a ‘disk diffusion’ drug susceptibility test that we videoed. The videos show complex spatio-genomic patterns redolent of, yet more complex than, Fisher’s predictions whereby a decelerating wave front of advantageous genes colonises towards the antibiotic source, forming bullseye patterns en route and leaving a wave back of bacterial sub-populations expressing AcrB at decreasing levels away from the drug source. qPCR data show thatE. colisited at rapidly-adapting spatial hotspots gain 2 additional copies ofacr, the operon that encodes AcrB, within 24h and imaging data show resistant sub-populations thrive most near the antibiotic source due to non-monotone relationships between inhibition due to antibiotic and distance from the source. In the spirit of Fisher, we provide an explicitly spatial nonlinear diffusion equation that exhibits these properties too. Finally, linear diffusion theory quantifies how the spatial extent of bacterial killing scales with increases in antibiotic dosage, predicting that microbes can survive chemotherapies that have been escalated to 250× the clinical dosage if the antibiotic is diffusion-limited.

AbstractSince bacteria lack a nucleus, the location of mRNA molecules is determined by the different characteristics of the encoded proteins, and the transcriptome is spatially arranged into cytosolic and membrane-associated mRNA. While translation of membrane protein-encoding mRNA has been studied in great mechanistic detail using biochemical methods, the spatiotemporal dynamics of this process remains poorly understood at the subcellular level. Here, we investigate the dynamics of individual fluorescently labelled mRNA molecules encoding the transmembrane serine chemoreceptor Tsr, to probe the mechanism of membrane protein translation. Analysis oftsrmRNA diffusion in the proximity of the plasma membrane revealed distinct diffusive modes that reflect the state of the mRNA molecule and its involvement in the process of active translation into the Sec secretion system. We find that the composition, and hence the fluidity of the membrane affects diffusion of membrane targeted mRNAs. Moreover, Tsr translation occurs in localized membrane regions, similar to eukaryotic hotspots. The hotspot localization coincides with the physical location of the transcribed gene, which itself is displaced towards the inner membrane. These findings suggest that inner membrane protein translation is a spatially defined process that occurs in hotspots, indicative of long-lived transertion sites. Our results show an additional layer of spatio-temporal structuring within the bacterial cell, thus revealing a qualitatively different understanding of the basic process of transcription and translation in bacteria.Significance StatementA large fraction of the bacterial proteome is directly synthesized into the inner membrane, and this process shapes the overall distribution of mRNA transcripts within the cell. Although highly dynamic in their nature, bacterial transcriptomes have mostly been studied in fixed cells. Here, we track individual mRNA molecules encoding the serine chemoreceptor in living bacterial cells and find that translation occurs in membrane hotspots that were previously exclusive to eukaryotes. Our results indicate an additional layer of spatio-temporal structuring within the bacterial cell that impacts our understanding of transcription and translation in bacteria.

The MazF toxin sequence-specifically cleaves single-stranded RNA upon various stressful conditions, and it is activated as a part of the mazEFtoxin-antitoxin module in Escherichia coli. Although autoregulation of mazEF expression through the MazE antitoxin-dependent transcriptional repression has been biochemically characterized, less is known about post-transcriptional autoregulation, as well as how both of these autoregulatory features affect growth of single cells during conditions that promote MazF production. Here, we demonstrate post-transcriptional autoregulation of mazFexpression dynamics by MazF cleaving its own transcript. Single-cell analyses of bacterial populations during ectopic MazF production indicated that two-level autoregulation of mazEF expression influences cell-to-cell growth rate heterogeneity. The increase in growth rate heterogeneity is governed by the MazE antitoxin, and tuned by the MazF-dependent mazFmRNA cleavage. Also, both autoregulatory features grant rapid exit from the stress caused by mazF overexpression. Time-lapse microscopy revealed that MazF-mediated cleavage of mazF mRNA leads to increased temporal variability in length of individual cells during ectopicmazFover-expression, as explained by a stochastic model indicating that mazEFmRNA cleavage underlies temporal fluctuations in MazF levels during stress.

Buffers are essential for diluting bacterial cultures for flow cytometry analysis in order to study bacterial physiology and gene expression parameters based on fluorescence signals. Using a variety of constitutively expressed fluorescent proteins in Escherichia coli K-12 strain MG1655, we found strong artifactual changes in fluorescence levels after dilution into the commonly used flow cytometry buffer phosphate-buffered saline (PBS) and two other buffer solutions, Tris-HCl and M9 salts. These changes appeared very rapidly after dilution, and were linked to increased membrane permeability and loss in cell viability. We observed buffer-related effects in several different E. coli strains, K-12, C and W, but not E. coli B, which can be partially explained by differences in lipopolysaccharide (LPS) and outer membrane composition. Supplementing the buffers with divalent cations responsible for outer membrane stability, Mg2+ and Ca2+, preserved fluorescence signals, membrane integrity and viability of E. coli. Thus, stabilizing the bacterial outer membrane is essential for precise and unbiased measurements of fluorescence parameters using flow cytometry.

Which properties of metabolic networks can be derived solely from stoichiometry? Predictive results have been obtained by flux balance analysis (FBA), by postulating that cells set metabolic fluxes to maximize growth rate. Here we consider a generalization of FBA to single-cell level using maximum entropy modeling, which we extend and test experimentally. Specifically, we define for Escherichia coli metabolism a flux distribution that yields the experimental growth rate: the model, containing FBA as a limit, provides a better match to measured fluxes and it makes a wide range of predictions: on flux variability, regulation, and correlations; on the relative importance of stoichiometry vs. optimization; on scaling relations for growth rate distributions. We validate the latter here with single-cell data at different sub-inhibitory antibiotic concentrations. The model quantifies growth optimization as emerging from the interplay of competitive dynamics in the population and regulation of metabolism at the level of single cells.

The bacterial flagellum is a self-assembling nanomachine. The external flagellar filament, several times longer than a bacterial cell body, is made of a few tens of thousands subunits of a single protein: flagellin. A fundamental problem concerns the molecular mechanism of how the flagellum grows outside the cell, where no discernible energy source is available. Here, we monitored the dynamic assembly of individual flagella using in situ labelling and real-time immunostaining of elongating flagellar filaments. We report that the rate of flagellum growth, initially ~1,700 amino acids per second, decreases with length and that the previously proposed chain mechanism does not contribute to the filament elongation dynamics. Inhibition of the proton motive force-dependent export apparatus revealed a major contribution of substrate injection in driving filament elongation. The combination of experimental and mathematical evidence demonstrates that a simple, injection-diffusion mechanism controls bacterial flagella growth outside the cell.

The molecular mechanisms underlying phenotypic variation in isogenic bacterial populations remain poorly understood.We report that AcrAB-TolC, the main multidrug efflux pump of Escherichia coli, exhibits a strong partitioning bias for old cell poles by a segregation mechanism that is mediated by ternary AcrAB-TolC complex formation. Mother cells inheriting old poles are phenotypically distinct and display increased drug efflux activity relative to daughters. Consequently, we find systematic and long-lived growth differences between mother and daughter cells in the presence of subinhibitory drug concentrations. A simple model for biased partitioning predicts a population structure of long-lived and highly heterogeneous phenotypes. This straightforward mechanism of generating sustained growth rate differences at subinhibitory antibiotic concentrations has implications for understanding the emergence of multidrug resistance in bacteria.

While we have good understanding of bacterial metabolism at the population level, we know little about the metabolic behavior of individual cells: do single cells in clonal populations sometimes specialize on different metabolic pathways? Such metabolic specialization could be driven by stochastic gene expression and could provide individual cells with growth benefits of specialization. We measured the degree of phenotypic specialization in two parallel metabolic pathways, the assimilation of glucose and arabinose. We grew Escherichia coli in chemostats, and used isotope-labeled sugars in combination with nanometer-scale secondary ion mass spectrometry and mathematical modeling to quantify sugar assimilation at the single-cell level. We found large variation in metabolic activities between single cells, both in absolute assimilation and in the degree to which individual cells specialize in the assimilation of different sugars. Analysis of transcriptional reporters indicated that this variation was at least partially based on cell-to-cell variation in gene expression. Metabolic differences between cells in clonal populations could potentially reduce metabolic incompatibilities between different pathways, and increase the rate at which parallel reactions can be performed.

Bacteria in groups vary individually, and interact with other bacteria and the environment to produce population-level patterns of gene expression. Investigating such behavior in detail requires measuring and controlling populations at the single-cell level alongside precisely specified interactions and environmental characteristics. Here we present an automated, programmable platform that combines image-based gene expression and growth measurements with on-line optogenetic expression control for hundreds of individual Escherichia coli cells over days, in a dynamically adjustable environment. This integrated platform broadly enables experiments that bridge individual and population behaviors. We demonstrate: (i) population structuring by independent closed-loop control of gene expression in many individual cells, (ii) cell-cell variation control during antibiotic perturbation, (iii) hybrid bio-digital circuits in single cells, and freely specifiable digital communication between individual bacteria. These examples showcase the potential for real-time integration of theoretical models with measurement and control of many individual cells to investigate and engineer microbial population behavior.

Parasitism creates selection for resistance mechanisms in host populations and is hypothesized to promote increased host evolvability. However, the influence of these traits on host evolution when parasites are no longer present is unclear. We used experimental evolution and whole-genome sequencing of Escherichia coli to determine the effects of past and present exposure to parasitic viruses (phages) on the spread of mutator alleles, resistance, and bacterial competitive fitness. We found that mutator alleles spread rapidly during adaptation to any of four different phage species, and this pattern was even more pronounced with multiple phages present simultaneously. However, hypermutability did not detectably accelerate adaptation in the absence of phages and recovery of fitness costs associated with resistance. Several lineages evolved phage resistance through elevated mucoidy, and during subsequent evolution in phage-free conditions they rapidly reverted to nonmucoid, phage-susceptible phenotypes. Genome sequencing revealed that this phenotypic reversion was achieved by additional genetic changes rather than by genotypic reversion of the initial resistance mutations. Insertion sequence (IS) elements played a key role in both the acquisition of resistance and adaptation in the absence of parasites; unlike single nucleotide polymorphisms, IS insertions were not more frequent in mutator lineages. Our results provide a genetic explanation for rapid reversion of mucoidy, a phenotype observed in other bacterial species including human pathogens. Moreover, this demonstrates that the types of genetic change underlying adaptation to fitness costs, and consequently the impact of evolvability mechanisms such as increased point-mutation rates, depend critically on the mechanism of resistance.

Restriction-modification (RM) systems represent a minimal and ubiquitous biological system of self/non-self discrimination in prokaryotes [1], which protects hosts from exogenous DNA [2]. The mechanism is based on the balance between methyltransferase (M) and cognate restriction endonuclease (R). M tags endogenous DNA as self by methylating short specific DNA sequences called restriction sites, whereas R recognizes unmethylated restriction sites as non-self and introduces a double-stranded DNA break [3]. Restriction sites are significantly underrepresented in prokaryotic genomes [4-7], suggesting that the discrimination mechanism is imperfect and occasionally leads to autoimmunity due to self-DNA cleavage (self-restriction) [8]. Furthermore, RM systems can promote DNA recombination [9] and contribute to genetic variation in microbial populations, thus facilitating adaptive evolution [10]. However, cleavage of self-DNA by RM systems as elements shaping prokaryotic genomes has not been directly detected, and its cause, frequency, and outcome are unknown. We quantify self-restriction caused by two RM systems of Escherichia coli and find that, in agreement with levels of restriction site avoidance, EcoRI, but not EcoRV, cleaves self-DNA at a measurable rate. Self-restriction is a stochastic process, which temporarily induces the SOS response, and is followed by DNA repair, maintaining cell viability. We find that RM systems with higher restriction efficiency against bacteriophage infections exhibit a higher rate of self-restriction, and that this rate can be further increased by stochastic imbalance between R and M. Our results identify molecular noise in RM systems as a factor shaping prokaryotic genomes.

In bacteria, replicative aging manifests as a difference in growth or survival between the two cells emerging from division. One cell can be regarded as an aging mother with a decreased potential for future survival and division, the other as a rejuvenated daughter. Here, we aimed at investigating some of the processes involved in aging in the bacterium Escherichia coli, where the two types of cells can be distinguished by the age of their cell poles. We found that certain changes in the regulation of the carbohydrate metabolism can affect aging. A mutation in the carbon storage regulator gene, csrA, leads to a dramatically shorter replicative lifespan; csrA mutants stop dividing once their pole exceeds an age of about five divisions. These old-pole cells accumulate glycogen at their old cell poles; after their last division, they do not contain a chromosome, presumably because of spatial exclusion by the glycogen aggregates. The new-pole daughters produced by these aging mothers are born young; they only express the deleterious phenotype once their pole is old. These results demonstrate how manipulations of nutrient allocation can lead to the exclusion of the chromosome and limit replicative lifespan in E. coli, and illustrate how mutations can have phenotypic effects that are specific for cells with old poles. This raises the question how bacteria can avoid the accumulation of such mutations in their genomes over evolutionary times, and how they can achieve the long replicative lifespans that have recently been reported.

Gene duplication is important in evolution, because it provides new raw material for evolutionary adaptations. Several existing hypotheses about the causes of duplicate retention and diversification differ in their emphasis on gene dosage, subfunctionalization, and neofunctionalization. Little experimental data exist on the relative importance of gene expression changes and changes in coding regions for the evolution of duplicate genes. Furthermore, we do not know how strongly the environment could affect this importance. To address these questions, we performed evolution experiments with the TEM-1 beta lactamase gene in Escherichia coli to study the initial stages of duplicate gene evolution in the laboratory. We mimicked tandem duplication by inserting two copies of the TEM-1 gene on the same plasmid. We then subjected these copies to repeated cycles of mutagenesis and selection in various environments that contained antibiotics in different combinations and concentrations. Our experiments showed that gene dosage is the most important factor in the initial stages of duplicate gene evolution, and overshadows the importance of point mutations in the coding region. © 2014 the Society for the Study of Evolution.

High relatedness among interacting individuals has generally been considered a precondition for the evolution of altruism. However, kin-selection theory also predicts the evolution of altruism when relatedness is low, as long as the cost of the altruistic act is minor compared with its benefit. Here, we demonstrate evidence for a low-cost altruistic act in bacteria. We investigated Escherichia coli responding to the attack of an obligately lytic phage by committing suicide in order to prevent parasite transmission to nearby relatives. We found that bacterial suicide provides large benefits to survivors at marginal costs to committers. The cost of suicide was low, because infected cells are moribund, rapidly dying upon phage infection, such that no more opportunity for reproduction remains. As a consequence of its marginal cost, host suicide was selectively favoured even when relatedness between committers and survivors approached zero. Altogether, our findings demonstrate that low-cost suicide can evolve with ease, represents an effective host-defence strategy, and seems to be widespread among microbes. Moreover, low-cost suicide might also occur in higher organisms as exemplified by infected social insect workers leaving the colony to die in isolation. © 2013 the Author(s) Published by the Royal Society. All rights reserved.

Essential genes code for fundamental cellular functions required for the viability of an organism. For this reason, essential genes are often highly conserved across organisms. However, this is not always the case: orthologues of genes that are essential in one organism are sometimes not essential in other organisms or are absent from their genomes. This suggests that, in the course of evolution, essential genes can be rendered nonessential. How can a gene become non-essential? Here we used genetic manipulation to deplete the products of 26 different essential genes in Escherichia coli. This depletion results in a lethal phenotype, which could often be rescued by the overexpression of a non-homologous, non-essential gene, most likely through replacement of the essential function. We also show that, in a smaller number of cases, the essential genes can be fully deleted from the genome, suggesting that complete functional replacement is possible. Finally, we show that essential genes whose function can be replaced in the laboratory are more likely to be non-essential or not present in other taxa. These results are consistent with the notion that patterns of evolutionary conservation of essential genes are influenced by their compensability-that is, by how easily they can be functionally replaced, for example through increased expression of other genes. © 2012 Bergmiller et al.

A number of recent experiments at the single-cell level have shown that genetically identical bacteria that live in homogeneous environments often show a substantial degree of phenotypic variation between cells. Often, this variation is attributed to stochastic aspects of biology-the fact that many biological processes involve small numbers of molecules and are thus inherently variable. However, not all variation between cells needs to be stochastic in nature; one deterministic process that could be important for cell variability in some bacterial species is the age of the cell poles. Working with the alphaproteobacterium Methylobacterium extorquens, we monitored individuals in clonally growing populations over several divisions and determined the pole age, cell size, and interdivision intervals of individual cells. We observed the high levels of variation in cell size and the timing of cell division that have been reported before. A substantial fraction of this variation could be explained by each cell's pole age and the pole age of its mother: cell size increased with increasing pole age, and the interval between cell divisions decreased. A theoretical model predicted that populations governed by such processes will quickly reach a stable distribution of different age and size classes. These results show that the pole age distribution in bacterial populations can contribute substantially to cellular individuality. In addition, they raise questions about functional differences between cells of different ages and the coupling of cell division to cell size. © 2011, American Society for Microbiology.

BACKGROUND: the essential Escherichia coli gene ygjD belongs to a universally conserved group of genes whose function has been the focus of a number of recent studies. Here, we put ygjD under control of an inducible promoter, and used time-lapse microscopy and single cell analysis to investigate the phenotypic consequences of the depletion of YgjD protein from growing cells. RESULTS: We show that loss of YgjD leads to a marked decrease in cell size and termination of cell division. The transition towards smaller size occurs in a controlled manner: cell elongation and cell division remain coupled, but cell size at division decreases. We also find evidence that depletion of YgjD leads to the synthesis of the intracellular signaling molecule (p)ppGpp, inducing a cellular reaction resembling the stringent response. Concomitant deletion of the relA and spoT genes - leading to a strain that is uncapable of synthesizing (p)ppGpp - abrogates the decrease in cell size, but does not prevent termination of cell division upon YgjD depletion. CONCLUSIONS: Depletion of YgjD protein from growing cells leads to a decrease in cell size that is contingent on (p)ppGpp, and to a termination of cell division. The combination of single-cell timelapse microscopy and statistical analysis can give detailed insights into the phenotypic consequences of the loss of essential genes, and can thus serve as a new tool to study the function of essential genes.

Strains of Escherichia coli lacking MalQ (maltodextrin glucanotransferase or amylomaltase) are endogenously induced for the maltose regulon by maltotriose that is derived from the degradation of glycogen (glycogen-dependent endogenous induction). A high level of induction was dependent on the presence of MalP, maltodextrin phosphorylase, while expression was counteracted by MalZ, maltodextrin glucosidase. Glycogen-derived endogenous induction was sensitive to high osmolarity. This osmodependence was caused by MalZ. malZ, the gene encoding this enzyme, was found to be induced by high osmolarity even in the absence of MalT, the central regulator of all mal genes. The osmodependent expression of malZ was neither RpoS nor OmpR dependent. In contrast, the malPQ operon, whose expression was also increased at a high osmolarity, was partially dependent on RpoS. In the absence of glycogen, residual endogenous induction of the mal genes that is sensitive to increasing osmolarity can still be observed. This glycogen-independent endogenous induction is not understood, and it is not affected by altering the expression of MalP, MalQ, and MalZ. In particular, its independence from MalZ suggests that the responsible inducer is not maltotriose. Copyright © 2005, American Society for Microbiology. All Rights Reserved.