Stochastic timing in gene expression for simple regulatory strategies

Timing is essential for many cellular processes, from cellular responses to external stimuli to the cell cycle and circadian clocks. Many of these processes are based on gene expression. For example, an activated gene may be required to reach in a precise time a threshold level of expression that triggers a specific downstream process. However, gene expression is subject to stochastic fluctuations, naturally inducing an uncertainty in this threshold-crossing time with potential consequences on biological functions and phenotypes. Here, we consider such "timing fluctuations", and we ask how they can be controlled. Our analytical estimates and simulations show that, for an induced gene, timing variability is minimal if the threshold level of expression is approximately half of the steady-state level. Timing fuctuations can be reduced by increasing the transcription rate, while they are insensitive to the translation rate. In presence of self-regulatory strategies, we show that self-repression reduces timing noise for threshold levels that have to be reached quickly, while selfactivation is optimal at long times. These results lay a framework for understanding stochasticity of endogenous systems such as the cell cycle, as well as for the design of synthetic trigger circuits.Comment: 10 pages, 5 figure

It takes two to tango: NAD+ and sirtuins in aging/longevity control

AbstractThe coupling of nicotinamide adenine dinucleotide (NAD+) breakdown and protein deacylation is a unique feature of the family of proteins called ‘sirtuins.’ This intimate connection between NAD+ and sirtuins has an ancient origin and provides a mechanistic foundation that translates the regulation of energy metabolism into aging and longevity control in diverse organisms. Although the field of sirtuin research went through intensive controversies, an increasing number of recent studies have put those controversies to rest and fully established the significance of sirtuins as an evolutionarily conserved aging/longevity regulator. The tight connection between NAD+ and sirtuins is regulated at several different levels, adding further complexity to their coordination in metabolic and aging/longevity control. Interestingly, it has been demonstrated that NAD+ availability decreases over age, reducing sirtuin activities and affecting the communication between the nucleus and mitochondria at a cellular level and also between the hypothalamus and adipose tissue at a systemic level. These dynamic cellular and systemic processes likely contribute to the development of age-associated functional decline and the pathogenesis of diseases of aging. To mitigate these age-associated problems, supplementation of key NAD+ intermediates is currently drawing significant attention. In this review article, we will summarize these important aspects of the intimate connection between NAD+ and sirtuins in aging/longevity control.</jats:p

Cell-type specific analysis of translating RNAs in developing flowers reveals new levels of control

Determining both the expression levels of mRNA and the regulation of its translation is important in understanding specialized cell functions. In this study, we describe both the expression profiles of cells within spatiotemporal domains of the Arabidopsis thaliana flower and the post-transcriptional regulation of these mRNAs, at nucleotide resolution. We express a tagged ribosomal protein under the promoters of three master regulators of flower development. By precipitating tagged polysomes, we isolated cell type specific mRNAs that are probably translating, and quantified those mRNAs through deep sequencing. Cell type comparisons identified known cell-specific transcripts and uncovered many new ones, from which we inferred cell type-specific hormone responses, promoter motifs and coexpressed cognate binding factor candidates, and splicing isoforms. By comparing translating mRNAs with steady-state overall transcripts, we found evidence for widespread post-transcriptional regulation at both the intron splicing and translational stages. Sequence analyses identified structural features associated with each step. Finally, we identified a new class of noncoding RNAs associated with polysomes. Findings from our profiling lead to new hypotheses in the understanding of flower development

Functional roles for noise in genetic circuits

The genetic circuits that regulate cellular functions are subject to stochastic fluctuations, or ‘noise’, in the levels of their components. Noise, far from just a nuisance, has begun to be appreciated for its essential role in key cellular activities. Noise functions in both microbial and eukaryotic cells, in multicellular development, and in evolution. It enables coordination of gene expression across large regulons, as well as probabilistic differentiation strategies that function across cell populations. At the longest timescales, noise may facilitate evolutionary transitions. Here we review examples and emerging principles that connect noise, the architecture of the gene circuits in which it is present, and the biological functions it enables. We further indicate some of the important challenges and opportunities going forward

Discovery of time-delayed gene regulatory networks based on temporal gene expression profiling

BACKGROUND: It is one of the ultimate goals for modern biological research to fully elucidate the intricate interplays and the regulations of the molecular determinants that propel and characterize the progression of versatile life phenomena, to name a few, cell cycling, developmental biology, aging, and the progressive and recurrent pathogenesis of complex diseases. The vast amount of large-scale and genome-wide time-resolved data is becoming increasing available, which provides the golden opportunity to unravel the challenging reverse-engineering problem of time-delayed gene regulatory networks. RESULTS: In particular, this methodological paper aims to reconstruct regulatory networks from temporal gene expression data by using delayed correlations between genes, i.e., pairwise overlaps of expression levels shifted in time relative each other. We have thus developed a novel model-free computational toolbox termed TdGRN (Time-delayed Gene Regulatory Network) to address the underlying regulations of genes that can span any unit(s) of time intervals. This bioinformatics toolbox has provided a unified approach to uncovering time trends of gene regulations through decision analysis of the newly designed time-delayed gene expression matrix. We have applied the proposed method to yeast cell cycling and human HeLa cell cycling and have discovered most of the underlying time-delayed regulations that are supported by multiple lines of experimental evidence and that are remarkably consistent with the current knowledge on phase characteristics for the cell cyclings. CONCLUSION: We established a usable and powerful model-free approach to dissecting high-order dynamic trends of gene-gene interactions. We have carefully validated the proposed algorithm by applying it to two publicly available cell cycling datasets. In addition to uncovering the time trends of gene regulations for cell cycling, this unified approach can also be used to study the complex gene regulations related to the development, aging and progressive pathogenesis of a complex disease where potential dependences between different experiment units might occurs

Mapping Dynamic Histone Acetylation Patterns to Gene Expression in Nanog-depleted Murine Embryonic Stem Cells

Embryonic stem cells (ESC) have the potential to self-renew indefinitely and to differentiate into any of the three germ layers. The molecular mechanisms for self-renewal, maintenance of pluripotency and lineage specification are poorly understood, but recent results point to a key role for epigenetic mechanisms. In this study, we focus on quantifying the impact of histone 3 acetylation (H3K9,14ac) on gene expression in murine embryonic stem cells. We analyze genome-wide histone acetylation patterns and gene expression profiles measured over the first five days of cell differentiation triggered by silencing Nanog, a key transcription factor in ESC regulation. We explore the temporal and spatial dynamics of histone acetylation data and its correlation with gene expression using supervised and unsupervised statistical models. On a genome-wide scale, changes in acetylation are significantly correlated to changes in mRNA expression and, surprisingly, this coherence increases over time. We quantify the predictive power of histone acetylation for gene expression changes in a balanced cross-validation procedure. In an in-depth study we focus on genes central to the regulatory network of Mouse ESC, including those identified in a recent genome-wide RNAi screen and in the PluriNet, a computationally derived stem cell signature. We find that compared to the rest of the genome, ESC-specific genes show significantly more acetylation signal and a much stronger decrease in acetylation over time, which is often not reflected in an concordant expression change. These results shed light on the complexity of the relationship between histone acetylation and gene expression and are a step forward to dissect the multilayer regulatory mechanisms that determine stem cell fate.Comment: accepted at PLoS Computational Biolog

Rate-limiting Steps in Transcription Initiation are Key Regulatory Mechanisms of Escherichia coli Gene Expression Dynamics

In all living organisms, the “blueprints of life” are documented in the genetic material. This material is composed of genes, which are regions of DNA coding for proteins. To produce proteins, cells read the information on the DNA with the help of molecular machines, such as RNAp holoenzymes and a factors. Proteins carry out the cellular functions required for survival and, as such, cells deal with challenging environments by adjusting their gene expression pattern. For this, cells constantly perform decision- making processes of whether or not to actively express a protein, based on intracellular and environmental cues. In Escherichia coli, gene expression is mostly regulated at the stage of transcription initiation. Although most of its regulatory molecules have been identified, the dynamics and regulation of this step remain elusive. Due to a limited number of specific regulatory molecules in the cells, the stochastic fluctuations of these molecular numbers can result in a sizeable temporal change in the numbers of transcription outputs (RNA and proteins) and have consequences on the phenotype of the cells. To understand the dynamics of this process, one should study the activity of the gene by tracking mRNA and protein production events at a detailed level. Recent advancements in single-molecule detection techniques have been used to image and track individually labeled fluorescent macromolecules of living cells. This allows investigating the intermolecular dynamics under any given condition. In this thesis, by using in vivo, single-RNA time-lapse microscopy techniques along with stochastic modelling techniques, we studied the kinetics of multi-rate limiting steps in the transcription process of multiple promoters, in various conditions. Specifically, first, we established a novel method of dissecting transcription in Escherichia coli that combines state-of-the-art microscopy measurements and model fitting techniques to construct detailed models of the rate-limiting steps governing the in vivo transcription initiation of a synthetic Lac-ara-1 promoter. After that, we estimated the duration of the closed and open complex formation, accounting for the rate of reversibility of the first step. From this, we also estimated the duration of periods of promoter inactivity, from which we were able to determine the contribution from each step to the distribution of intervals between consecutive RNA productions in individual cells. Second, using the above method, we studied the a factor selective mechanisms for indirect regulation of promoters whose transcription is primarily initiated by RNAp holoenzymes carrying a70. From the analysis, we concluded that, in E. coli, a promoter’s responsiveness to indirect regulation by a factor competition is determined by its sequence-dependent, dynamically regulated multi-step initiation kinetics. Third, we investigated the effects of extrinsic noise, arising from cell-to-cell variability in cellular components, on the single-cell distribution of RNA numbers, in the context of cell lineages. For this, first, we used stochastic models to predict the variability in the numbers of molecules involved in upstream processes. The models account for the intake of inducers from the environment, which acts as a transient source of variability in RNA production numbers, as well as for the variability in the numbers of molecular species controlling transcription of an active promoter, which acts as a constant source of variability in RNA numbers. From measurement analysis, we demonstrated the existence of lineage-to-lineage variability in gene activation times and mean transcription rates. Finally, we provided evidence that this can be explained by differences in the kinetics of the rate-limiting steps in transcription and of the induction scheme, from which it is possible to conclude that these variabilities differ between promoters and inducers used. Finally, we studied how the multi-rate limiting steps in the transcription initiation are capable of tuning the asymmetry and tailedness of the distribution of time intervals between consecutive RNA production events in individual cells. For this, first, we considered a stochastic model of transcription initiation and predicted that the asymmetry and tailedness in the distribution of intervals between consecutive RNA production events can differ by tuning the rate-limiting steps in transcription. Second, we validated the model with measurements from single-molecule RNA microscopy of transcription kinetics of multiple promoters in multiple conditions. Finally, from our results, we concluded that the skewness and kurtosis in RNA and protein production kinetics are subject to regulation by the kinetics of the steps in transcription initiation and affect the single-cell distributions of RNAs and, thus, proteins. We further showed that this regulation can significantly affect the probability of RNA and protein numbers to cross specific thresholds. Overall, the studies conducted in this thesis are expected to contribute to a better understanding of the dynamic process of bacterial gene expression. The advanced data and image analysis techniques and novel stochastic modeling approaches that we developed during the course of these studies, will allow studying in detail the in vivo regulation of multi-rate limiting steps of transcription initiation of any given promoter. In addition, by tuning the kinetics of the rate-limiting steps in the transcription initiation as executed here should allow engineering new promoters, with predefined RNA and, thus, protein production dynamics in Escherichia coli

Bridging the gap between omics and earth system science to better understand how environmental change impacts marine microbes

The advent of genomic-, transcriptomic- and proteomic-based approaches has revolutionized our ability to describe marine microbial communities, including biogeography, metabolic potential and diversity, mechanisms of adaptation, and phylogeny and evolutionary history. New interdisciplinary approaches are needed to move from this descriptive level to improved quantitative, process-level understanding of the roles of marine microbes in biogeochemical cycles and of the impact of environmental change on the marine microbial ecosystem. Linking studies at levels from the genome to the organism, to ecological strategies and organism and ecosystem response, requires new modelling approaches. Key to this will be a fundamental shift in modelling scale that represents micro-organisms from the level of their macromolecular components. This will enable contact with omics data sets and allow acclimation and adaptive response at the phenotype level (i.e. traits) to be simulated as a combination of fitness maximization and evolutionary constraints. This way forward will build on ecological approaches that identify key organism traits and systems biology approaches that integrate traditional physiological measurements with new insights from omics. It will rely on developing an improved understanding of ecophysiology to understand quantitatively environmental controls on microbial growth strategies. It will also incorporate results from experimental evolution studies in the representation of adaptation. The resulting ecosystem-level models can then evaluate our level of understanding of controls on ecosystem structure and function, highlight major gaps in understanding and help prioritize areas for future research programs. Ultimately, this grand synthesis should improve predictive capability of the ecosystem response to multiple environmental drivers