Droplet-based methodology for investigating bacterial population dynamics in response to phage exposure.

An alarming rise in antimicrobial resistance worldwide has spurred efforts into the search for alternatives to antibiotic treatments. The use of bacteriophages, bacterial viruses harmless to humans, represents a promising approach with potential to treat bacterial infections (phage therapy). Recent advances in microscopy-based single-cell techniques have allowed researchers to develop new quantitative methodologies for assessing the interactions between bacteria and phages, especially the ability of phages to eradicate bacterial pathogen populations and to modulate growth of both commensal and pathogen populations. Here we combine droplet microfluidics with fluorescence time-lapse microscopy to characterize the growth and lysis dynamics of the bacterium Escherichia coli confined in droplets when challenged with phage. We investigated phages that promote lysis of infected E. coli cells, specifically, a phage species with DNA genome, T7 (Escherichia virus T7) and two phage species with RNA genomes, MS2 (Emesvirus zinderi) and Qß (Qubevirus durum). Our microfluidic trapping device generated and immobilized picoliter-sized droplets, enabling stable imaging of bacterial growth and lysis in a temperature-controlled setup. Temporal information on bacterial population size was recorded for up to 25 h, allowing us to determine growth rates of bacterial populations and helping us uncover the extent and speed of phage infection. In the long-term, the development of novel microfluidic single-cell and population-level approaches will expedite research towards fundamental understanding of the genetic and molecular basis of rapid phage-induced lysis and eco-evolutionary aspects of bacteria-phage dynamics, and ultimately help identify key factors influencing the success of phage therapy.

Quantifying heterologous gene expression during ectopic MazF production in Escherichia coli.

OBJECTIVE: MazF is a sequence-specific endoribonuclease-toxin of the MazEF toxin-antitoxin system. MazF cleaves single-stranded ribonucleic acid (RNA) regions at adenine-cytosine-adenine (ACA) sequences in the bacterium Escherichia coli. The MazEF system has been used in various biotechnology and synthetic biology applications. In this study, we infer how ectopic mazF overexpression affects production of heterologous proteins. To this end, we quantified the levels of fluorescent proteins expressed in E. coli from reporters translated from the ACA-containing or ACA-less messenger RNAs (mRNAs). Additionally, we addressed the impact of the 5'-untranslated region of these reporter mRNAs under the same conditions by comparing expression from mRNAs that comprise (canonical mRNA) or lack this region (leaderless mRNA). RESULTS: Flow cytometry analysis indicates that during mazF overexpression, fluorescent proteins are translated from the canonical as well as leaderless mRNAs. Our analysis further indicates that longer mazF overexpression generally increases the concentration of fluorescent proteins translated from ACA-less mRNAs, however it also substantially increases bacterial population heterogeneity. Finally, our results suggest that the strength and duration of mazF overexpression should be optimized for each experimental setup, to maximize the heterologous protein production and minimize the amount of phenotypic heterogeneity in bacterial populations, which is unfavorable in biotechnological processes.

Nutrient and salt depletion synergistically boosts glucose metabolism in individual Escherichia coli cells.

The interaction between a cell and its environment shapes fundamental intracellular processes such as cellular metabolism. In most cases growth rate is treated as a proximal metric for understanding the cellular metabolic status. However, changes in growth rate might not reflect metabolic variations in individuals responding to environmental fluctuations. Here we use single-cell microfluidics-microscopy combined with transcriptomics, proteomics and mathematical modelling to quantify the accumulation of glucose within Escherichia coli cells. In contrast to the current consensus, we reveal that environmental conditions which are comparatively unfavourable for growth, where both nutrients and salinity are depleted, increase glucose accumulation rates in individual bacteria and population subsets. We find that these changes in metabolic function are underpinned by variations at the translational and posttranslational level but not at the transcriptional level and are not dictated by changes in cell size. The metabolic response-characteristics identified greatly advance our fundamental understanding of the interactions between bacteria and their environment and have important ramifications when investigating cellular processes where salinity plays an important role.

Autoregulation of bacterial gene expression: lessons from the MazEF toxin-antitoxin system.

Autoregulation is the direct modulation of gene expression by the product of the corresponding gene. Autoregulation of bacterial gene expression has been mostly studied at the transcriptional level, when a protein acts as the cognate transcriptional repressor. A recent study investigating dynamics of the bacterial toxin-antitoxin MazEF system has shown how autoregulation at both the transcriptional and post-transcriptional levels affects the heterogeneity of Escherichia coli populations. Toxin-antitoxin systems hold a crucial but still elusive part in bacterial response to stress. This perspective highlights how these modules can also serve as a great model system for investigating basic concepts in gene regulation. However, as the genomic background and environmental conditions substantially influence toxin activation, it is important to study (auto)regulation of toxin-antitoxin systems in well-defined setups as well as in conditions that resemble the environmental niche.

Autoregulation of mazEF expression underlies growth heterogeneity in bacterial populations.

The MazF toxin sequence-specifically cleaves single-stranded RNA upon various stressful conditions, and it is activated as a part of the mazEF toxin-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 mazF expression 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 mazF mRNA 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 ectopic mazF overexpression, as explained by a stochastic model indicating that mazEF mRNA cleavage underlies temporal fluctuations in MazF levels during stress.

Cell-to-cell variation and specialization in sugar metabolism in clonal bacterial populations.

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.

MazF activation promotes translational heterogeneity of the grcA mRNA in Escherichia coli populations.

Bacteria adapt to adverse environmental conditions by altering gene expression patterns. Recently, a novel stress adaptation mechanism has been described that allows Escherichia coli to alter gene expression at the post-transcriptional level. The key player in this regulatory pathway is the endoribonuclease MazF, the toxin component of the toxin-antitoxin module mazEF that is triggered by various stressful conditions. In general, MazF degrades the majority of transcripts by cleaving at ACA sites, which results in the retardation of bacterial growth. Furthermore, MazF can process a small subset of mRNAs and render them leaderless by removing their ribosome binding site. MazF concomitantly modifies ribosomes, making them selective for the translation of leaderless mRNAs. In this study, we employed fluorescent reporter-systems to investigate mazEF expression during stressful conditions, and to infer consequences of the mRNA processing mediated by MazF on gene expression at the single-cell level. Our results suggest that mazEF transcription is maintained at low levels in single cells encountering adverse conditions, such as antibiotic stress or amino acid starvation. Moreover, using the grcA mRNA as a model for MazF-mediated mRNA processing, we found that MazF activation promotes heterogeneity in the grcA reporter expression, resulting in a subpopulation of cells with increased levels of GrcA reporter protein.

Analysis of fluorescent reporters indicates heterogeneity in glucose uptake and utilization in clonal bacterial populations.

BACKGROUND: In this study, we aimed at investigating heterogeneity in the expression of metabolic genes in clonal populations of Escherichia coli growing on glucose as the sole carbon source. Different metabolic phenotypes can arise in these clonal populations through variation in the expression of glucose transporters and metabolic enzymes. First, we focused on the glucose transporters PtsG and MglBAC to analyze the diversity of glucose uptake strategies. Second, we analyzed phenotypic variation in the expression of genes involved in gluconeogenesis and acetate scavenging (as acetate is formed and excreted during bacterial growth on glucose), which can reveal, for instance, phenotypic subpopulations that cross-feed through the exchange of acetate. In these experiments, E. coli MG1655 strains containing different transcriptional GFP reporters were grown in chemostats and reporter expression was measured with flow cytometry. RESULTS: Our results suggest heterogeneous expression of metabolic genes in bacterial clonal populations grown in glucose environments. The two glucose transport systems exhibited different level of heterogeneity. The majority of the bacterial cells expressed the reporters for both glucose transporters MglBAC and PtsG and a small fraction of cells only expressed the reporter for Mgl. At a low dilution rate, signals from transcriptional reporters for acetyl-CoA synthetase Acs and phosphoenolpyruvate carboxykinase Pck indicated that almost all cells expressed the genes that are part of acetate utilization and the gluconeogenesis pathway, respectively. Possible co-existence of two phenotypic subpopulations differing in acs expression occurred at the threshold of the switch to overflow metabolism. The overflow metabolism results in the production of acetate and has been previously reported to occur at intermediate dilution rates in chemostats with high concentration of glucose in the feed. CONCLUSIONS: Analysis of the heterogeneous expression of reporters for genes involved in glucose and acetate metabolism raises new question whether different metabolic phenotypes are expressed in clonal populations growing in continuous cultures fed on glucose as the initially sole carbon source.

Proteomic properties reveal phyloecological clusters of Archaea.

In this study, we propose a novel way to describe the variety of environmental adaptations of Archaea. We have clustered 57 Archaea by using a non-redundant set of proteomic features, and verified that the clusters correspond to environmental adaptations to the archaeal habitats. The first cluster consists dominantly of hyperthermophiles and hyperthermoacidophilic aerobes. The second cluster joins together halophilic and extremely halophilic Archaea, while the third cluster contains mesophilic (mostly methanogenic) Archaea together with thermoacidophiles. The non-redundant subset of proteomic features was found to consist of five features: the ratio of charged residues to uncharged, average protein size, normalized frequency of beta-sheet, normalized frequency of extended structure and number of hydrogen bond donors. We propose this clustering to be termed phyloecological clustering. This approach could give additional insights into relationships among archaeal species that may be hidden by sole phylogenetic analysis.

A genome-wide analysis of promoter-mediated phenotypic noise in Escherichia coli.

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as "phenotypic noise." In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alone.

Proteome sequence features carry signatures of the environmental niche of prokaryotes.

BACKGROUND: Prokaryotic environmental adaptations occur at different levels within cells to ensure the preservation of genome integrity, proper protein folding and function as well as membrane fluidity. Although specific composition and structure of cellular components suitable for the variety of extreme conditions has already been postulated, a systematic study describing such adaptations has not yet been performed. We therefore explored whether the environmental niche of a prokaryote could be deduced from the sequence of its proteome. Finally, we aimed at finding the precise differences between proteome sequences of prokaryotes from different environments. RESULTS: We analyzed the proteomes of 192 prokaryotes from different habitats. We collected detailed information about the optimal growth conditions of each microorganism. Furthermore, we selected 42 physico-chemical properties of amino acids and computed their values for each proteome. Further, on the same set of features we applied two fundamentally different machine learning methods, Support Vector Machines and Random Forests, to successfully classify between bacteria and archaea, halophiles and non-halophiles, as well as mesophiles, thermophiles and mesothermophiles. Finally, we performed feature selection by using Random Forests. CONCLUSIONS: To our knowledge, this is the first time that three different classification cases (domain of life, halophilicity and thermophilicity) of proteome adaptation are successfully performed with the same set of 42 features. The characteristic features of a specific adaptation constitute a signature that may help understanding the mechanisms of adaptation to extreme environments.

Unstructured hydrophilic sequences in prokaryotic proteomes correlate with dehydration tolerance and host association.

Water loss or desiccation is among the most life-threatening stresses. It leads to DNA double-strand breakage, protein aggregation, cell shrinkage, and low water activity precluding all biological functions. Yet, in all kingdoms of life, rare organisms are resistant to desiccation through prevention or reversibility of such damage. Here, we explore possible hallmarks of prokaryotic desiccation tolerance in their proteomes. The content of unstructured, low complexity (LC) regions was analyzed in a total of 460 bacterial and archaeal proteomes. It appears that species endowed with proteomes abundant in unstructured hydrophilic LC regions are desiccation-tolerant or sporulating bacteria, halophilic archaea and bacteria, or host-associated species. In the desiccation- and radiation-resistant bacterium Deinococcus radiodurans, most proteins that contain large hydrophilic LC regions have unassigned function, but those with known function are mostly involved in diverse cellular recovery processes. Such LC regions are typically absent in orthologous proteins in desiccation-sensitive species. D. radiodurans encodes also special LC proteins, akin to those associated with desiccation resistance of plant seeds and some plants and animals. Therefore, we postulate that large unstructured hydrophilic LC regions and proteins provide for cellular resistance to dehydration and we discuss mechanisms of their protective activity.