Molecular Origins of Transcriptional Heterogeneity in Diazotrophic Klebsiella oxytoca.

Phenotypic heterogeneity in clonal bacterial batch cultures has been shown for a range of bacterial systems; however, the molecular origins of such heterogeneity and its magnitude are not well understood. Under conditions of extreme low-nitrogen stress in the model diazotroph Klebsiella oxytoca, we found remarkably high heterogeneity of nifHDK gene expression, which codes for the structural genes of nitrogenase, one key enzyme of the global nitrogen cycle. This heterogeneity limited the bulk observed nitrogen-fixing capacity of the population. Using dual-probe, single-cell RNA fluorescent in situ hybridization, we correlated nifHDK expression with that of nifLA and glnK-amtB, which code for the main upstream regulatory components. Through stochastic transcription models and mutual information analysis, we revealed likely molecular origins for heterogeneity in nitrogenase expression. In the wild type and regulatory variants, we found that nifHDK transcription was inherently bursty, but we established that noise propagation through signaling was also significant. The regulatory gene glnK had the highest discernible effect on nifHDK variance, while noise from factors outside the regulatory pathway were negligible. Understanding the basis of inherent heterogeneity of nitrogenase expression and its origins can inform biotechnology strategies seeking to enhance biological nitrogen fixation. Finally, we speculate on potential benefits of diazotrophic heterogeneity in natural soil environments. IMPORTANCE Nitrogen is an essential micronutrient for both plant and animal life and naturally exists in both reactive and inert chemical forms. Modern agriculture is heavily reliant on nitrogen that has been "fixed" into a reactive form via the energetically expensive Haber-Bosch process, with significant environmental consequences. Nitrogen-fixing bacteria provide an alternative source of fixed nitrogen for use in both biotechnological and agricultural settings, but this relies on a firm understanding of how the fixation process is regulated within individual bacterial cells. We examined the cell-to-cell variability in the nitrogen-fixing behavior of Klebsiella oxytoca, a free-living bacterium. The significance of our research is in identifying not only the presence of marked variability but also the specific mechanisms that give rise to it. This understanding gives insight into both the evolutionary advantages of variable behavior as well as strategies for biotechnological applications.

The route to transcription initiation determines the mode of transcriptional bursting in E. coli.

Transcription is fundamentally noisy, leading to significant heterogeneity across bacterial populations. Noise is often attributed to burstiness, but the underlying mechanisms and their dependence on the mode of promotor regulation remain unclear. Here, we measure E. coli single cell mRNA levels for two stress responses that depend on bacterial sigma factors with different mode of transcription initiation (σ70 and σ54). By fitting a stochastic model to the observed mRNA distributions, we show that the transition from low to high expression of the σ70-controlled stress response is regulated via the burst size, while that of the σ54-controlled stress response is regulated via the burst frequency. Therefore, transcription initiation involving σ54 differs from other bacterial systems, and yields bursting kinetics characteristic of eukaryotic systems.

Exactly solvable models of stochastic gene expression.

Stochastic models are key to understanding the intricate dynamics of gene expression. However, the simplest models that only account for active and inactive states of a gene fail to capture common observations in both prokaryotic and eukaryotic organisms. Here, we consider multistate models of gene expression that generalize the canonical Telegraph process and are capable of capturing the joint effects of transcription factors, heterochromatin state, and DNA accessibility (or, in prokaryotes, sigma-factor activity) on transcript abundance. We propose two approaches for solving classes of these generalized systems. The first approach offers a fresh perspective on a general class of multistate models and allows us to "decompose" more complicated systems into simpler processes, each of which can be solved analytically. This enables us to obtain a solution of any model from this class. Next, we develop an approximation method based on a power series expansion of the stationary distribution for an even broader class of multistate models of gene transcription. We further show that models from both classes cannot have a heavy-tailed distribution in the absence of extrinsic noise. The combination of analytical and computational solutions for these realistic gene expression models also holds the potential to design synthetic systems and control the behavior of naturally evolved gene expression systems in guiding cell-fate decisions.

Extrinsic Noise and Heavy-Tailed Laws in Gene Expression.

Noise in gene expression is one of the hallmarks of life at the molecular scale. Here we derive analytical solutions to a set of models describing the molecular mechanisms underlying transcription of DNA into RNA. Our ansatz allows us to incorporate the effects of extrinsic noise-encompassing factors external to the transcription of the individual gene-and discuss the ramifications for heterogeneity in gene product abundance that has been widely observed in single cell data. Crucially, we are able to show that heavy-tailed distributions of RNA copy numbers cannot result from the intrinsic stochasticity in gene expression alone, but must instead reflect extrinsic sources of variability.

Transition state characteristics during cell differentiation.

Models describing the process of stem-cell differentiation are plentiful, and may offer insights into the underlying mechanisms and experimentally observed behaviour. Waddington's epigenetic landscape has been providing a conceptual framework for differentiation processes since its inception. It also allows, however, for detailed mathematical and quantitative analyses, as the landscape can, at least in principle, be related to mathematical models of dynamical systems. Here we focus on a set of dynamical systems features that are intimately linked to cell differentiation, by considering exemplar dynamical models that capture important aspects of stem cell differentiation dynamics. These models allow us to map the paths that cells take through gene expression space as they move from one fate to another, e.g. from a stem-cell to a more specialized cell type. Our analysis highlights the role of the transition state (TS) that separates distinct cell fates, and how the nature of the TS changes as the underlying landscape changes-change that can be induced by e.g. cellular signaling. We demonstrate that models for stem cell differentiation may be interpreted in terms of either a static or transitory landscape. For the static case the TS represents a particular transcriptional profile that all cells approach during differentiation. Alternatively, the TS may refer to the commonly observed period of heterogeneity as cells undergo stochastic transitions.