Limits on information transduction through amplitude and frequency regulation of transcription factor activity

Automated adjustment of metabolic niches enables the control of natural and engineered microbial co-cultures

Much attention has focused on understanding microbial interactions leading to stable co-cultures. In this work, substrate pulsing was performed to promote better control of the metabolic niches (MNs) corresponding to each species, leading to the continuous co-cultivation of diverse microbial organisms. We used a cell-machine interface, which allows adjustment of the temporal profile of two MNs according to a rhythm, ensuring the successive growth of two species, in our case, a yeast and a bacterium. The resulting approach, called 'automated adjustment of metabolic niches' (AAMN), was effective for stabilizing both cooperative and competitive co-cultures. AAMN can be considered an enabling technology for the deployment of co-cultures in bioprocesses, demonstrated here based on the continuous bioproduction of p-coumaric acid. The data accumulated suggest that AAMN could be used not only for a wider range of biological systems, but also to gain fundamental insights into microbial interaction mechanisms.

Information Processing in Biochemical Networks

Living systems are characterized by controlled flows of matter, energy, and information. While the biophysics community has productively engaged with the first two, addressing information flows has been more challenging, with some scattered success in evolutionary theory and a more coherent track record in neuroscience. Nevertheless, interdisciplinary work of the past two decades at the interface of biophysics, quantitative biology, and engineering has led to an emerging mathematical language for describing information flows at the molecular scale. This is where the central processes of life unfold: from detection and transduction of environmental signals to the readout or copying of genetic information and the triggering of adaptive cellular responses. Such processes are coordinated by complex biochemical reaction networks that operate at room temperature, are out of equilibrium, and use low copy numbers of diverse molecular species with limited interaction specificity. Here we review how flows of information through biochemical networks can be formalized using information-theoretic quantities, quantified from data, and computed within various modeling frameworks. Optimization of information flows is presented as a candidate design principle that navigates the relevant time, energy, crosstalk, and metabolic constraints to predict reliable cellular signaling and gene regulation architectures built of individually noisy components.

Amplitude and frequency encoding result in qualitatively distinct informational landscapes in cell signaling

Cells continuously sense their surroundings to detect modifications and generate responses. Very often changes in extracellular concentrations initiate signaling cascades that eventually result in changes in gene expression. Increasing stimulus strengths can be encoded in increasing concentration amplitudes or increasing activation frequencies of intermediaries of the pathway. In this paper we show that the different way in which amplitude and frequency encoding map environmental changes endow cells with qualitatively different information transmission capabilities. While amplitude encoding is optimal for a limited range of stimuli strengths, frequency encoding can transmit information with equal reliability over much broader ranges. The qualitative difference between the two strategies stems from the scale invariant discriminating power of the first transducing step in frequency codification. The apparently redundant combination of both strategies in some cell types may then serve the purpose of expanding the span over which stimulus strengths can be reliably discriminated. In this paper we discuss a possible example of this mechanism in yeast.

CellSAM: A Foundation Model for Cell Segmentation

Cells are a fundamental unit of biological organization, and identifying them in imaging data - cell segmentation - is a critical task for various cellular imaging experiments. While deep learning methods have led to substantial progress on this problem, most models are specialist models that work well for specific domains but cannot be applied across domains or scale well with large amounts of data. In this work, we present CellSAM, a universal model for cell segmentation that generalizes across diverse cellular imaging data. CellSAM builds on top of the Segment Anything Model (SAM) by developing a prompt engineering approach for mask generation. We train an object detector, CellFinder, to automatically detect cells and prompt SAM to generate segmentations. We show that this approach allows a single model to achieve human-level performance for segmenting images of mammalian cells, yeast, and bacteria collected across various imaging modalities. We show that CellSAM has strong zero-shot performance and can be improved with a few examples via few-shot learning. Additionally, we demonstrate how CellSAM can be applied across diverse bioimage analysis workflows. A deployed version of CellSAM is available at https://cellsam.deepcell.org/.

Regulated resource reallocation is transcriptionally hard wired into the yeast stress response

Many organisms maintain generalized stress responses activated by adverse conditions. Although details vary, a common theme is the redirection of transcriptional and translational capacity away from growth-promoting genes and toward defense genes. Yet the precise roles of these coupled programs are difficult to dissect. Here we investigated Saccharomyces cerevisiae responding to salt as a model stressor. We used molecular, genomic, and single-cell microfluidic methods to examine the interplay between transcription factors Msn2 and Msn4 that induce stress-defense genes and Dot6 and Tod6 that transiently repress growth-promoting genes during stress. Surprisingly, loss of Dot6/Tod6 led to slower acclimation to salt, whereas loss of Msn2/4 produced faster growth during stress. This supports a model where transient repression of growth-promoting genes accelerates the Msn2/4 response, which is essential for acquisition of subsequent peroxide tolerance. Remarkably, we find that Msn2/4 regulate DOT6 mRNA production, influence Dot6 activation dynamics, and are required for full repression of growth-promoting genes. Thus, Msn2/4 directly regulate resource reallocation needed to mount their own response. We discuss broader implications for common stress responses across organisms.

Dynamics of positional information in the vertebrate neural tube

In developing embryos, cells acquire distinct identities depending on their position in a tissue. Secreted signalling molecules, known as morphogens, act as long-range cues to provide the spatial information that controls these cell fate decisions. In several tissues, both the level and the duration of morphogen signalling appear to be important for determining cell fates. This is the case in the forming vertebrate nervous system where antiparallel morphogen gradients pattern the dorsal-ventral axis by partitioning the tissue into sharply delineated domains of molecularly distinct neural progenitors. How information in the gradients is decoded to generate precisely positioned boundaries of gene expression remains an open question. Here, we adopt tools from information theory to quantify the positional information in the neural tube and investigate how temporal changes in signalling could influence positional precision. The results reveal that the use of signalling dynamics, as well as the signalling level, substantially increases the precision possible for the estimation of position from morphogen gradients. This analysis links the dynamics of opposing morphogen gradients with precise pattern formation and provides an explanation for why time is used to impart positional information.

Dynamic phosphorylation of Hcm1 promotes fitness in chronic stress

Cell survival depends upon the ability to adapt to changing environments. Environmental stressors trigger an acute stress response program that rewires cell physiology, downregulates proliferation genes and pauses the cell cycle until the cell adapts. Here, we show that dynamic phosphorylation of the yeast cell cycle-regulatory transcription factor Hcm1 is required to maintain fitness in chronic stress. Hcm1 is activated by cyclin dependent kinase (CDK) and inactivated by the phosphatase calcineurin (CN) in response to stressors that signal through increases in cytosolic Ca2+. Expression of a constitutively active, phosphomimetic Hcm1 mutant reduces fitness in stress, suggesting Hcm1 inactivation is required. However, a comprehensive analysis of Hcm1 phosphomutants revealed that Hcm1 activity is also important to survive stress, demonstrating that Hcm1 activity must be toggled on and off to promote gene expression and fitness. These results suggest that dynamic control of cell cycle regulators is critical for survival in stressful environments.

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information-theoretic framework to quantify the distribution of sensing abilities from single-cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an 'average cell'. We verify these predictions using live-cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information-theoretic framework will significantly improve our understanding of how cells sense in their environment.

Information Transmission via Molecular Communication in Astrobiological Environments

The ubiquity of information transmission via molecular communication between cells is comprehensively documented on Earth; this phenomenon might even have played a vital role in the origin(s) and early evolution of life. Motivated by these considerations, a simple model for molecular communication entailing the diffusion of signaling molecules from transmitter to receiver is elucidated. The channel capacity C (maximal rate of information transmission) and an optimistic heuristic estimate of the actual information transmission rate ℐ are derived for this communication system; the two quantities, especially the latter, are demonstrated to be broadly consistent with laboratory experiments and more sophisticated theoretical models. The channel capacity exhibits a potentially weak dependence on environmental parameters, whereas the actual information transmission rate may scale with the intercellular distance d as ℐ ∝ d-4 and could vary substantially across settings. These two variables are roughly calculated for diverse astrobiological environments, ranging from Earth's upper oceans (C ∼ 3.1 × 103 bits/s; ℐ ∼ 4.7 × 10-2 bits/s) and deep sea hydrothermal vents (C ∼ 4.2 × 103 bits/s; ℐ ∼ 1.2 × 10-1 bits/s) to the hydrocarbon lakes and seas of Titan (C ∼ 3.8 × 103 bits/s; ℐ ∼ 2.6 × 10-1 bits/s).

A guide to ERK dynamics, part 2: downstream decoding

Signaling by the extracellular signal-regulated kinase (ERK) pathway controls many cellular processes, including cell division, death, and differentiation. In this second installment of a two-part review, we address the question of how the ERK pathway exerts distinct and context-specific effects on multiple processes. We discuss how the dynamics of ERK activity induce selective changes in gene expression programs, with insights from both experiments and computational models. With a focus on single-cell biosensor-based studies, we summarize four major functional modes for ERK signaling in tissues: adjusting the size of cell populations, gradient-based patterning, wave propagation of morphological changes, and diversification of cellular gene expression states. These modes of operation are disrupted in cancer and other related diseases and represent potential targets for therapeutic intervention. By understanding the dynamic mechanisms involved in ERK signaling, there is potential for pharmacological strategies that not only simply inhibit ERK, but also restore functional activity patterns and improve disease outcomes.

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information theoretic framework to quantify the distribution of sensing abilities from single cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an " average cell ". We verify these predictions using live cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information theoretic framework will significantly improve our understanding of how cells sense in their environment.

Advances in automated and reactive flow cytometry for synthetic biotechnology

Automated flow cytometry (FC) has been initially considered for bioprocess monitoring and optimization. More recently, new physical and software interfaces have been made available, facilitating the access to this technology for labs and industries. It also comes with new capabilities, such as being able to act on the cultivation conditions based on population data. This approach, known as reactive FC, extended the range of applications of automated FC to bioprocess control and the stabilization of cocultures, but also to the broad field of synthetic and systems biology for the characterization of gene circuits. However, several issues must be addressed before automated and reactive FC can be considered standard and modular technologies.

Modulating Biological Rhythms: A Noncomputational Strategy Harnessing Nonlinearity and Decoupling Frequency and Amplitude

Understanding and achieving concurrent modulation of amplitude and frequency, particularly adjusting one quantity and simultaneously sustaining the other at an invariant level, are of paramount importance for complex biophysical systems, including the signal pathway where different frequency indicates different upstream signal yielding a certain downstream physiological function while different amplitude further determines different efficacy of a downstream output. However, such modulators with clearly described and universally valid mechanisms are still lacking. Here, we rigorously propose an easy-to-use control strategy containing only one or two steps, leveraging the nonlinearity in the modulated systems to decouple frequency and amplitude in a noncomputational manner. The strategy's efficacy is demonstrated using representative biochemical systems and, thus, it could be potentially applicable to modulating rhythms in experiments of biochemistry and synthetic biology.

NFAT and NF-κB dynamically co-regulate TCR and CAR signaling responses in human T cells

While it has been established that the responses of T cells to antigens are combinatorially regulated by multiple signaling pathways, it remains elusive what mechanisms cells utilize to quantitatively modulate T cell responses during pathway integration. Here, we show that two key pathways in T cell signaling, calcium/nuclear factor of activated T cells (NFAT) and protein kinase C (PKC)/nuclear factor κB (NF-κB), integrate through a dynamic and combinatorial strategy to fine-tune T cell response genes. At the cis-regulatory level, the two pathways integrate through co-binding of NFAT and NF-κB to immune response genes. Pathway integration is further regulated temporally, where T cell receptor (TCR) and chimeric antigen receptor (CAR) activation signals modulate the temporal relationships between the nuclear localization dynamics of NFAT and NF-κB. Such physical and temporal integrations together contribute to distinct modes of expression modulation for genes. Thus, the temporal relationships between regulators can be modulated to affect their co-targets during immune responses, underscoring the importance of dynamic combinatorial regulation in cellular signaling.

Nearly maximal information gain due to time integration in central dogma reactions

Living cells process information about their environment through the central dogma processes of transcription and translation, which drive the cellular response to stimuli. Here, we study the transfer of information from environmental input to the transcript and protein expression levels. Evaluation of both experimental and analogous simulation data reveals that transcription and translation are not two simple information channels connected in series. Instead, we demonstrate that the central dogma reactions often create a time-integrating information channel, where the translation channel receives and integrates multiple outputs from the transcription channel. This information channel model of the central dogma provides new information-theoretic selection criteria for the central dogma rate constants. Using the data for four well-studied species we show that their central dogma rate constants achieve information gain because of time integration while also keeping the loss because of stochasticity in translation relatively low (<0.5 bits).

Transcription factor localization dynamics and DNA binding drive distinct promoter interpretations

Environmental information may be encoded in the temporal dynamics of transcription factor (TF) activation and subsequently decoded by gene promoters to enact stimulus-specific gene expression programs. Previous studies of this behavior focused on the encoding and decoding of information in TF nuclear localization dynamics, yet cells control the activity of TFs in myriad ways, including by regulating their ability to bind DNA. Here, we use light-controlled mutants of the yeast TF Msn2 as a model system to investigate how promoter decoding of TF localization dynamics is affected by changes in the ability of the TF to bind DNA. We find that yeast promoters directly decode the light-controlled localization dynamics of Msn2 and that the effects of changing Msn2 affinity on that decoding behavior are highly promoter dependent, illustrating how cells could regulate TF localization dynamics and DNA binding in concert for improved control of gene expression.

Modulation of transcription factor dynamics allows versatile information transmission

Cells detect changes in their environment and generate responses, often involving changes in gene expression. In this paper we use information theory and a simple transcription model to analyze whether the resulting gene expression serves to identify extracellular stimuli and assess their intensity when they are encoded in the amplitude, duration or frequency of pulses of a transcription factor's nuclear concentration (or activation state). We find, for all cases, that about three ranges of input strengths can be distinguished and that maximum information transmission occurs for fast and high activation threshold promoters. The three input modulation modes differ in the sensitivity to changes in the promoters parameters. Frequency modulation is the most sensitive and duration modulation, the least. This is key for signal identification: there are promoter parameters that yield a relatively high information transmission for duration or amplitude modulation and a much smaller value for frequency modulation. The reverse situation cannot be found with a single promoter transcription model. Thus, pulses of transcription factors can selectively activate the "frequency-tuned" promoter while prolonged nuclear accumulation would activate promoters of all three modes simultaneously. Frequency modulation is therefore highly selective and better suited than the other encoding modes for signal identification without requiring other mediators of the transduction process.

Epigenetics as a mediator of plasticity in cancer

The concept of an epigenetic landscape describing potential cellular fates arising from pluripotent cells, first advanced by Conrad Waddington, has evolved in light of experiments showing nondeterministic outcomes of regulatory processes and mathematical methods for quantifying stochasticity. In this Review, we discuss modern approaches to epigenetic and gene regulation landscapes and the associated ideas of entropy and attractor states, illustrating how their definitions are both more precise and relevant to understanding cancer etiology and the plasticity of cancerous states. We address the interplay between different types of regulatory landscapes and how their changes underlie cancer progression. We also consider the roles of cellular aging and intrinsic and extrinsic stimuli in modulating cellular states and how landscape alterations can be quantitatively mapped onto phenotypic outcomes and thereby used in therapy development.

The regulatory mechanism of the yeast osmoresponse under different glucose concentrations

Cells constantly respond to environmental changes by modulating gene expression programs. These responses may demand substantial costs and, thus, affect cell growth. Understanding the regulation of these processes represents a key question in biology and biotechnology. Here, we studied the responses to osmotic stress in glucose-limited environments. By analyzing seventeen osmotic stress-induced genes and stress-activated protein kinase Hog1, we found that cells exhibited stronger osmotic gene expression response and larger integral of Hog1 nuclear localization during adaptation to osmotic stress under glucose-limited conditions than under glucose-rich conditions. We proposed and verified that in glucose-limited environment, glycolysis intermediates (representing "reserve flux") were limited, which required cells to express more glycerol-production enzymes for stress adaptation. Consequently, the regulatory mechanism of osmoresponse was derived in the presence and absence of such reserve flux. Further experiments suggested that this reserve flux-dependent stress-defense strategy may be a general principle under nutrient-limited environments.

Quantifying information of intracellular signaling: progress with machine learning

Cells convey information about their extracellular environment to their core functional machineries. Studying the capacity of intracellular signaling pathways to transmit information addresses fundamental questions about living systems. Here, we review how information-theoretic approaches have been used to quantify information transmission by signaling pathways that are functionally pleiotropic and subject to molecular stochasticity. We describe how recent advances in machine learning have been leveraged to address the challenges of complex temporal trajectory datasets and how these have contributed to our understanding of how cells employ temporal coding to appropriately adapt to environmental perturbations.

Exploiting Information and Control Theory for Directing Gene Expression in Cell Populations

Microbial populations can adapt to adverse environmental conditions either by appropriately sensing and responding to the changes in their surroundings or by stochastically switching to an alternative phenotypic state. Recent data point out that these two strategies can be exhibited by the same cellular system, depending on the amplitude/frequency of the environmental perturbations and on the architecture of the genetic circuits involved in the adaptation process. Accordingly, several mitigation strategies have been designed for the effective control of microbial populations in different contexts, ranging from biomedicine to bioprocess engineering. Technically, such control strategies have been made possible by the advances made at the level of computational and synthetic biology combined with control theory. However, these control strategies have been applied mostly to synthetic gene circuits, impairing the applicability of the approach to natural circuits. In this review, we argue that it is possible to expand these control strategies to any cellular system and gene circuits based on a metric derived from this information theory, i.e., mutual information (MI). Indeed, based on this metric, it should be possible to characterize the natural frequency of any gene circuits and use it for controlling gene circuits within a population of cells.

Synthetic gene networks recapitulate dynamic signal decoding and differential gene expression

Cells live in constantly changing environments and employ dynamic signaling pathways to transduce information about the signals they encounter. However, the mechanisms by which dynamic signals are decoded into appropriate gene expression patterns remain poorly understood. Here, we devise networked optogenetic pathways that achieve dynamic signal processing functions that recapitulate cellular information processing. Exploiting light-responsive transcriptional regulators with differing response kinetics, we build a falling edge pulse detector and show that this circuit can be employed to demultiplex dynamically encoded signals. We combine this demultiplexer with dCas9-based gene networks to construct pulsatile signal filters and decoders. Applying information theory, we show that dynamic multiplexing significantly increases the information transmission capacity from signal to gene expression state. Finally, we use dynamic multiplexing for precise multidimensional regulation of a heterologous metabolic pathway. Our results elucidate design principles of dynamic information processing and provide original synthetic systems capable of decoding complex signals for biotechnological applications.

MYC amplifies gene expression through global changes in transcription factor dynamics

The MYC oncogene has been studied for decades, yet there is still intense debate over how this transcription factor controls gene expression. Here, we seek to answer these questions with an in vivo readout of discrete events of gene expression in single cells. We engineered an optogenetic variant of MYC (Pi-MYC) and combined this tool with single-molecule RNA and protein imaging techniques to investigate the role of MYC in modulating transcriptional bursting and transcription factor binding dynamics in human cells. We find that the immediate consequence of MYC overexpression is an increase in the duration rather than in the frequency of bursts, a functional role that is different from the majority of human transcription factors. We further propose that the mechanism by which MYC exerts global effects on the active period of genes is by altering the binding dynamics of transcription factors involved in RNA polymerase II complex assembly and productive elongation.

Mapping the dynamic transfer functions of eukaryotic gene regulation

Biological information can be encoded within the dynamics of signaling components, which has been implicated in a broad range of physiological processes including stress response, oncogenesis, and stem cell differentiation. To study the complexity of information transfer across the eukaryotic promoter, we screened 119 dynamic conditions-modulating the pulse frequency, amplitude, and pulse width of light-regulating the binding of an epigenome editor to a fluorescent reporter. This system revealed tunable gene expression and filtering behaviors and provided a quantification of the limit to the amount of information that can be reliably transferred across a single promoter as ∼1.7 bits. Using a library of over 100 orthogonal chromatin regulators, we further determined that chromatin state could be used to tune mutual information and expression levels, as well as completely alter the input-output transfer function of the promoter. This system unlocks the information-rich content of eukaryotic gene regulation.

Evaluation of Lee et al.: Clarity and interpretation of mutual information in promoter transfer functions

One snapshot of the peer review process for "Mapping the dynamic transfer functions of eukaryotic gene regulation" (Lee et al., 2021) appears below.

A frequency-amplitude coordinator and its optimal energy consumption for biological oscillators

Biorhythm including neuron firing and protein-mRNA interaction are fundamental activities with diffusive effect. Their well-balanced spatiotemporal dynamics are beneficial for healthy sustainability. Therefore, calibrating both anomalous frequency and amplitude of biorhythm prevents physiological dysfunctions or diseases. However, many works were devoted to modulate frequency exclusively whereas amplitude is usually ignored, although both quantities are equally significant for coordinating biological functions and outputs. Especially, a feasible method coordinating the two quantities concurrently and precisely is still lacking. Here, for the first time, we propose a universal approach to design a frequency-amplitude coordinator rigorously via dynamical systems tools. We consider both spatial and temporal information. With a single well-designed coordinator, they can be calibrated to desired levels simultaneously and precisely. The practical usefulness and efficacy of our method are demonstrated in representative neuronal and gene regulatory models. We further reveal its fundamental mechanism and optimal energy consumption providing inspiration for biorhythm regulation in future.

Information transmission in a two-step cascade: interplay of activation and repression

We present an information-theoretic formalism to study signal transduction in four architectural variants of a model two-step cascade with increasing input population. Our results categorize these four types into two classes depending upon the effect of activation and repression on mutual information, net synergy, and signal-to-noise ratio. Using the Gaussian framework and linear noise approximation, we derive the analytic expressions for these metrics to establish their underlying relationships in terms of the biochemical parameters. We also verify our approximations through stochastic simulations.

A functionally divergent intrinsically disordered region underlying the conservation of stochastic signaling

Stochastic signaling dynamics expand living cells' information processing capabilities. An increasing number of studies report that regulators encode information in their pulsatile dynamics. The evolutionary mechanisms that lead to complex signaling dynamics remain uncharacterized, perhaps because key interactions of signaling proteins are encoded in intrinsically disordered regions (IDRs), whose evolution is difficult to analyze. Here we focused on the IDR that controls the stochastic pulsing dynamics of Crz1, a transcription factor in fungi downstream of the widely conserved calcium signaling pathway. We find that Crz1 IDRs from anciently diverged fungi can all respond transiently to calcium stress; however, only Crz1 IDRs from the Saccharomyces clade support pulsatility, encode extra information, and rescue fitness in competition assays, while the Crz1 IDRs from distantly related fungi do none of the three. On the other hand, we find that Crz1 pulsing is conserved in the distantly related fungi, consistent with the evolutionary model of stabilizing selection on the signaling phenotype. Further, we show that a calcineurin docking site in a specific part of the IDRs appears to be sufficient for pulsing and show evidence for a beneficial increase in the relative calcineurin affinity of this docking site. We propose that evolutionary flexibility of functionally divergent IDRs underlies the conservation of stochastic signaling by stabilizing selection.

Gene expression at a single-molecule level: implications for myelodysplastic syndromes and acute myeloid leukemia

Nongenetic heterogeneity, or gene expression stochasticity, is an important source of variability in biological systems. With the advent and improvement of single molecule resolution technologies, it has been shown that transcription dynamics and resultant transcript number fluctuations generate significant cell-to-cell variability that has important biological effects and may contribute substantially to both tissue homeostasis and disease. In this respect, the pathophysiology of stem cell-derived malignancies such as acute myeloid leukemia and myelodysplastic syndromes, which has historically been studied at the ensemble level, may require reevaluation. To that end, it is our aim in this review to highlight the results of recent single-molecule, biophysical, and systems studies of gene expression dynamics, with the explicit purpose of demonstrating how the insights from these basic science studies may help inform and progress the field of leukemia biology and, ultimately, research into novel therapies.

Approximating gene transcription dynamics using steady-state formulas

Understanding how genes in a single cell respond to dynamically changing signals has been a central question in stochastic gene transcription research. Recent studies have generated massive steady-state or snapshot mRNA distribution data of individual cells, and inferred a large spectrum of kinetic transcription parameters under varying conditions. However, there have been few algorithms to convert these static data into the temporal variation of kinetic rates. Real-time imaging has been developed to monitor stochastic transcription processes at the single-cell level, but the immense technicality has prevented its application to most endogenous loci in mammalian cells. In this article, we introduced a stochastic gene transcription model with variable kinetic rates induced by unstable cellular conditions. We approximated the transcription dynamics using easily obtained steady-state formulas in the model. We tested the approximation against experimental data in both prokaryotic and eukaryotic cells and further solidified the conditions that guarantee the robustness of the method. The method can be easily implemented to provide convenient tools for quantifying dynamic kinetics and mechanisms underlying the widespread static transcription data, and may shed a light on circumventing the limitation of current bursting data on transcriptional real-time imaging.

Towards a 'Spot On' Understanding of Transcription in the Nucleus

Regulation of transcription by RNA Polymerase II (RNAPII) is a rapidly evolving area of research. Technological developments in microscopy have revealed insight into the dynamics, structure, and localization of transcription components within single cells. A frequent observation in many studies is the appearance of 'spots' in cell nuclei associated with the transcription process. In this review we highlight studies that characterize the temporal and spatial characteristics of these spots, examine possible pitfalls in interpreting these kind of imaging data, and outline directions where single-cell imaging may advance in ways to further our understanding of transcription regulation.

Molecular signatures of silencing suppression degeneracy from a complex RNA virus

As genomic architectures become more complex, they begin to accumulate degenerate and redundant elements. However, analyses of the molecular mechanisms underlying these genetic architecture features remain scarce, especially in compact but sufficiently complex genomes. In the present study, we followed a proteomic approach together with a computational network analysis to reveal molecular signatures of protein function degeneracy from a plant virus (as virus-host protein-protein interactions). We employed affinity purification coupled to mass spectrometry to detect several host factors interacting with two proteins of Citrus tristeza virus (p20 and p25) that are known to function as RNA silencing suppressors, using an experimental system of transient expression in a model plant. The study was expanded by considering two different isolates of the virus, and some key interactions were confirmed by bimolecular fluorescence complementation assays. We found that p20 and p25 target a common set of plant proteins including chloroplastic proteins and translation factors. Moreover, we noted that even specific targets of each viral protein overlap in function. Notably, we identified argonaute proteins (key players in RNA silencing) as reliable targets of p20. Furthermore, we found that these viral proteins preferentially do not target hubs in the host protein interactome, but elements that can transfer information by bridging different parts of the interactome. Overall, our results demonstrate that two distinct proteins encoded in the same viral genome that overlap in function also overlap in their interactions with the cell proteome, thereby highlighting an overlooked connection from a degenerate viral system.

Paralogs in the PKA Regulon Traveled Different Evolutionary Routes to Divergent Expression in Budding Yeast

Functional divergence of duplicate genes, or paralogs, is an important driver of novelty in evolution. In the model yeast Saccharomyces cerevisiae, there are 547 paralog gene pairs that survive from an interspecies Whole Genome Hybridization (WGH) that occurred ~100MYA. In this work, we report that ~1/6th (110) of these WGH paralogs pairs (or ohnologs) are differentially expressed with a striking pattern upon Protein Kinase A (PKA) inhibition. One member of each pair in this group has low basal expression that increases upon PKA inhibition, while the other has moderate and unchanging expression. For these genes, expression of orthologs upon PKA inhibition in the non-WGH species Kluyveromyces lactis and for PKA-related stresses in other budding yeasts shows unchanging expression, suggesting that lack of responsiveness to PKA was likely the typical ancestral phenotype prior to duplication. Promoter sequence analysis across related budding yeast species further revealed that the subsequent emergence of PKA-dependence took different evolutionary routes. In some examples, regulation by PKA and differential expression appears to have arisen following the WGH, while in others, regulation by PKA appears to have arisen in one of the two parental lineages prior to the WGH. More broadly, our results illustrate the unique opportunities presented by a WGH event for generating functional divergence by bringing together two parental lineages with separately evolved regulation into one species. We propose that functional divergence of two ohnologs can be facilitated through such regulatory divergence.

Pulses and delays, anticipation and memory: seeing bacterial stress responses from a single-cell perspective

Stress responses are crucial for bacteria to survive harmful conditions that they encounter in the environment. Although gene regulatory mechanisms underlying stress responses in bacteria have been thoroughly characterised for decades, recent advances in imaging technologies helped to uncover previously hidden dynamics and heterogeneity that become visible at the single-cell level. Despite the diversity of stress response mechanisms, certain dynamic regulatory features are frequently seen in single cells, such as pulses, delays, stress anticipation and memory effects. Often, these dynamics are highly variable across cells. While any individual cell may not achieve an optimal stress response, phenotypic diversity can provide a benefit at the population level. In this review, we highlight microscopy studies that offer novel insights into how bacteria sense stress, regulate protective mechanisms, cope with response delays and prepare for future environmental challenges. These studies showcase developments in the single-cell imaging toolbox including gene expression reporters, FRET, super-resolution microscopy and single-molecule tracking, as well as microfluidic techniques to manipulate cells and create defined stress conditions.

The Oscillation Amplitude, Not the Frequency of Cytosolic Calcium, Regulates Apoptosis Induction

Although a rising concentration of cytosolic Ca²⁺ has long been recognized as an essential signal for apoptosis, the dynamical mechanisms by which Ca²⁺ regulates apoptosis are not clear yet. To address this, we constructed a computational model that integrates known biochemical reactions and can reproduce the dynamical behaviors of Ca²⁺-induced apoptosis as observed in experiments. Model analysis shows that oscillating Ca²⁺ signals first convert into gradual signals and eventually transform into a switch-like apoptotic response. Via the two processes, the apoptotic signaling pathway filters the frequency of Ca²⁺ oscillations effectively but instead responds acutely to their amplitude. Collectively, our results suggest that Ca²⁺ regulates apoptosis mainly via oscillation amplitude, rather than frequency, modulation. This study not only provides a comprehensive understanding of how oscillatory Ca²⁺ dynamically regulates the complex apoptotic signaling network but also presents a typical example of how Ca²⁺ controls cellular responses through amplitude modulation.

Optogenetic Control Reveals Differential Promoter Interpretation of Transcription Factor Nuclear Translocation Dynamics

Gene expression is thought to be affected not only by the concentration of transcription factors (TFs) but also the dynamics of their nuclear translocation. Testing this hypothesis requires direct control of TF dynamics. Here, we engineer CLASP, an optogenetic tool for rapid and tunable translocation of a TF of interest. Using CLASP fused to Crz1, we observe that, for the same integrated concentration of nuclear TF over time, changing input dynamics changes target gene expression: pulsatile inputs yield higher expression than continuous inputs, or vice versa, depending on the target gene. Computational modeling reveals that a dose-response saturating at low TF input can yield higher gene expression for pulsatile versus continuous input, and that multi-state promoter activation can yield the opposite behavior. Our integrated tool development and modeling approach characterize promoter responses to Crz1 nuclear translocation dynamics, extracting quantitative features that may help explain the differential expression of target genes.

CLASP is a modular optogenetic strategy to control the nuclear localization of transcription factors (TFs) and elicit gene expression from their cognate promoters. CLASP control of Crz1 nuclear localization, coupled with computational modeling, revealed how promoters can differentially decode dynamic transcription factor signals. The integrated strategy of CLASP development and modeling presents a generalized approach to causally investigate the transcriptional consequences of dynamic TF nuclear shuttling.

Information transmission from NFkB signaling dynamics to gene expression

The dynamic signal encoding paradigm suggests that information flows from the extracellular environment into specific signaling patterns (encoding) that are then read by downstream effectors to control cellular behavior. Previous work empirically quantified the information content of dynamic signaling patterns. However, whether this information can be faithfully transmitted to the gene expression level is unclear. Here we used NFkB signaling as a model to understand the accuracy of information transmission from signaling dynamics into gene expression. Using a detailed mathematical model, we simulated realistic NFkB signaling patterns with different degrees of variability. The NFkB patterns were used as an input to a simple gene expression model. Analysis of information transmission between ligand and NFkB and ligand and gene expression allows us to determine information loss in transmission between receptors to dynamic signaling patterns and between signaling dynamics to gene expression. Information loss could occur due to biochemical noise or due to a lack of specificity. We found that noise-free gene expression has very little information loss suggesting that gene expression can preserve specificity in NFkB patterns. As expected, the addition of noise to the gene expression model results in information loss. Interestingly, this effect can be mitigated by a specific choice of parameters that can substantially reduce information loss due to biochemical noise during gene expression. Overall our results show that the cellular capacity for information transmission from dynamic signaling patterns to gene expression can be high enough to preserve ligand specificity and thereby the accuracy of cellular response to environmental cues.

Coherent feed-forward loop acts as an efficient information transmitting motif

We present a theoretical formalism to study steady-state information transmission in a coherent type-1 feed-forward loop motif with an additive signal integration mechanism. Our construct allows a two-step cascade to be slowly transformed into a bifurcation network via a feed-forward loop, which is a prominent network motif. Using a Gaussian framework, we show that among these three network patterns, the feed-forward loop motif harnesses the maximum amount of Shannon mutual information fractions constructed between the final gene-product and each of the master and coregulators of the target gene. We also show that this feed-forward loop motif provides a substantially lower amount of noise in target gene expression, compared with the other two network structures. Our theoretical predictions, which remain invariant for a couple of parametric transformations, point out that the coherent type-1 feed-forward loop motif may qualify as a better decoder of environmental signals when compared with the other two network patterns in perspective.

Insights about collective decision-making at the genetic level

By living in a collective, individuals can share and aggregate information to base their decisions on the many rather than on the one, thereby increasing accuracy. But a collective can also be defined at the molecular level. In the following, we reason that genes, by working collectively, share fundamental features with social organisms, which ends, without invoking cognition, in wiser responses. For that, we compile into a single picture the terms redundancy, stochastic resonance, intrinsic and extrinsic noise, and cross-regulation.

On Information Extraction and Decoding Mechanisms Improved by Noisy Amplification in Signaling Pathways

The cells need to process information about extracellular stimuli. They encode, transmit and decode the information to elicit an appropriate response. Studies aimed at understanding how such information is decoded in the signaling pathways to generate a specific cellular response have become essential. Eukaryotic cells decode information through two different mechanisms: the feed-forward loop and the promoter affinity. Here, we investigate how these two mechanisms improve information transmission. A detailed comparison is made between the stochastic model of the MAPK/ERK pathway and a stochastic minimal decoding model. The maximal amount of transmittable information was computed. The results suggest that the decoding mechanism of the MAPK/ERK pathway improve the channel capacity because it behaves as a noisy amplifier. We show a positive dependence between the noisy amplification and the amount of information extracted. Additionally, we show that the extrinsic noise can be tuned to improve information transmission. This investigation has revealed that the feed-forward loop and the promoter affinity motifs extract information thanks to processes of amplification and noise addition. Moreover, the channel capacity is enhanced when both decoding mechanisms are coupled. Altogether, these findings suggest novel characteristics in how decoding mechanisms improve information transmission.

Complex genetic and epigenetic regulation deviates gene expression from a unifying global transcriptional program

Environmental or genetic perturbations lead to gene expression changes. While most analyses of these changes emphasize the presence of qualitative differences on just a few genes, we now know that changes are widespread. This large-scale variation has been linked to the exclusive influence of a global transcriptional program determined by the new physiological state of the cell. However, given the sophistication of eukaryotic regulation, we expect to have a complex architecture of specific control affecting this program. Here, we examine this architecture. Using data of Saccharomyces cerevisiae expression in different nutrient conditions, we first propose a five-sector genome partition, which integrates earlier models of resource allocation, as a framework to examine the deviations from the global control. In this scheme, we recognize invariant genes, whose regulation is dominated by physiology, specific genes, which substantially depart from it, and two additional classes that contain the frequently assumed growth-dependent genes. Whereas the invariant class shows a considerable absence of specific regulation, the rest is enriched by regulation at the level of transcription factors (TFs) and epigenetic modulators. We nevertheless find markedly different strategies in how these classes deviate. On the one hand, there are TFs that act in a unique way between partition constituents, and on the other, the action of chromatin modifiers is significantly diverse. The balance between regulatory strategies ultimately modulates the action of the general transcription machinery and therefore limits the possibility of establishing a unifying program of expression change at a genomic scale.

How can we understand expression changes observed as a result of environmental or genetic perturbations? This issue has been conventionally answered by examining small groups of genes whose expression becomes qualitatively altered after these perturbations. But this approach is too simplistic, as we now know that extensive variation is typically observed. To explain this variation, recent works proposed a model in which genome-wide changes are explained by the action of a general program of transcription. Our manuscript reasons that given the complexity of eukaryotic transcriptional control, a unifying program of regulation cannot be achievable. Instead, we propose within an integrated framework of resource allocation that a rich structure of deviations from it exists and that by characterizing these deviations we can fully appreciate large-scale expression change.

Estimating information in time-varying signals

Across diverse biological systems-ranging from neural networks to intracellular signaling and genetic regulatory networks-the information about changes in the environment is frequently encoded in the full temporal dynamics of the network nodes. A pressing data-analysis challenge has thus been to efficiently estimate the amount of information that these dynamics convey from experimental data. Here we develop and evaluate decoding-based estimation methods to lower bound the mutual information about a finite set of inputs, encoded in single-cell high-dimensional time series data. For biological reaction networks governed by the chemical Master equation, we derive model-based information approximations and analytical upper bounds, against which we benchmark our proposed model-free decoding estimators. In contrast to the frequently-used k-nearest-neighbor estimator, decoding-based estimators robustly extract a large fraction of the available information from high-dimensional trajectories with a realistic number of data samples. We apply these estimators to previously published data on Erk and Ca2+ signaling in mammalian cells and to yeast stress-response, and find that substantial amount of information about environmental state can be encoded by non-trivial response statistics even in stationary signals. We argue that these single-cell, decoding-based information estimates, rather than the commonly-used tests for significant differences between selected population response statistics, provide a proper and unbiased measure for the performance of biological signaling networks.

Information-theoretic analysis of multivariate single-cell signaling responses

Mathematical methods of information theory appear to provide a useful language to describe how stimuli are encoded in activities of signaling effectors. Exploring the information-theoretic perspective, however, remains conceptually, experimentally and computationally challenging. Specifically, existing computational tools enable efficient analysis of relatively simple systems, usually with one input and output only. Moreover, their robust and readily applicable implementations are missing. Here, we propose a novel algorithm, SLEMI—statistical learning based estimation of mutual information, to analyze signaling systems with high-dimensional outputs and a large number of input values. Our approach is efficient in terms of computational time as well as sample size needed for accurate estimation. Analysis of the NF-κB single—cell signaling responses to TNF-α reveals that NF-κB signaling dynamics improves discrimination of high concentrations of TNF-α with a relatively modest impact on discrimination of low concentrations. Provided R-package allows the approach to be used by computational biologists with only elementary knowledge of information theory.

In light of single-cell, live-imaging experiments understanding of how cells transmit information about identity and quantity of stimuli is incomplete. When exposed to the same stimulus individual cells exhibit substantial cell-to-cell heterogeneity. Besides, stimuli have been shown to regulate temporal profiles of signaling effectors. Therefore, it is, for instance, not entirely clear whether single-cell responses are binary or contain more information about the quantity of stimuli. The above questions resulted in a considerable interest to study cellular signaling within the framework of information theory. Unfortunately, the utilization of the information-theoretic perspective is handicapped in part by the lack of suitable methods that account for multivariate signaling data. Here, we propose a novel algorithm that breaks a considerable computational barrier by allowing the effective information-theoretic analysis of highly-dimensional single-cell measurements. Our approach is computationally efficient, robust and straightforward to use. Moreover, we provide a simple R-package implementation.

Multivariate information processing characterizes fitness of a cascaded gene-transcription machinery

We report that a genetic two-step activation cascade processes diverse flavors of information, e.g., synergy, redundancy, and unique information. Our computations measuring reduction in Shannon entropies and reduction in variances produce differently behaving absolute magnitudes of these informational flavors. We find that similarity can be brought in if these terms are evaluated in fractions with respect to corresponding total information. Each of the input signal and final gene-product is found to generate common or redundant information fractions (mostly) to predict each other, whereas they also complement one another to harness synergistic information fraction, predicting the intermediate biochemical species. For an optimally growing signal to maintain fixed steady-state abundance of activated downstream gene-products, the interaction information fractions for this cascade module shift from net-redundancy to information-independence.

Techniques for Studying Decoding of Single Cell Dynamics

Cells must be able to interpret signals they encounter and reliably generate an appropriate response. It has long been known that the dynamics of transcription factor and kinase activation can play a crucial role in selecting an individual cell's response. The study of cellular dynamics has expanded dramatically in the last few years, with dynamics being discovered in novel pathways, new insights being revealed about the importance of dynamics, and technological improvements increasing the throughput and capabilities of single cell measurements. In this review, we highlight the important developments in this field, with a focus on the methods used to make new discoveries. We also include a discussion on improvements in methods for engineering and measuring single cell dynamics and responses. Finally, we will briefly highlight some of the many challenges and avenues of research that are still open.

Limits to the rate of information transmission through the MAPK pathway

Two important signalling pathways of NF-κB and ERK transmit merely 1 bit of information about the level of extracellular stimulation. It is thus unclear how such systems can coordinate complex cell responses to external cues. We analyse information transmission in the MAPK/ERK pathway that converts both constant and pulsatile EGF stimulation into pulses of ERK activity. Based on an experimentally verified computational model, we demonstrate that, when input consists of sequences of EGF pulses, transmitted information increases nearly linearly with time. Thus, pulse-interval transcoding allows more information to be relayed than the amplitude–amplitude transcoding considered previously for the ERK and NF-κB pathways. Moreover, the information channel capacity C, or simply bitrate, is not limited by the bandwidth B = 1/τ, where τ ≈ 1 h is the relaxation time. Specifically, when the input is provided in the form of sequences of short binary EGF pulses separated by intervals that are multiples of τ/n (but not shorter than τ), then for n = 2, C ≈ 1.39 bit h⁻¹; and for n = 4, C ≈ 1.86 bit h⁻¹. The capability to respond to random sequences of EGF pulses enables cells to propagate spontaneous ERK activity waves across tissue.

Computational Strategies for Niche-Dependent Cell Conversion to Assist Stem Cell Therapy

The field of regenerative medicine has blossomed in recent decades. However, the ultimate goal of tissue regeneration - replacing damaged or aged cells with healthy functioning cells - still faces a number of challenges. In particular, better understanding of the role of the cellular niche in shaping stem cell phenotype and conversion would aid in improving current protocols for stem cell therapies. In this regard, the implementation of novel computational approaches that consider the niche effect on stem cells would be valuable. Here we discuss current problems in stem cell transplantation and rejuvenation, and we propose computational strategies to control niche-dependent cell conversion to overcome them.

A Noisy Analog-to-Digital Converter Connects Cytosolic Calcium Bursts to Transcription Factor Nuclear Localization Pulses in Yeast

Several examples of transcription factors that show stochastic, unsynchronized pulses of nuclear localization have been described. Here we show that under constant calcium stress, nuclear localization pulses of the transcription factor Crz1 follow stochastic variations in cytosolic calcium concentration. We find that the size of the stochastic calcium bursts is positively correlated with the number of subsequent Crz1 pulses. Based on our observations, we propose a simple stochastic model of how the signaling pathway converts a constant external calcium concentration into a digital number of Crz1 pulses in the nucleus, due to the time delay from nuclear transport and the stochastic decoherence of individual Crz1 molecule dynamics. We find support for several additional predictions of the model and suggest that stochastic input to nuclear transport may produce noisy digital responses to analog signals in other signaling systems.

Robustness and Information Transfer within IL-6-induced JAK/STAT Signalling

Cellular communication via intracellular signalling pathways is crucial. Expression and activation of signalling proteins is heterogenous between isogenic cells of the same cell-type. However, mechanisms evolved to enable sufficient communication and to ensure cellular functions. We use information theory to clarify mechanisms facilitating IL-6-induced JAK/STAT signalling despite cell-to-cell variability. We show that different mechanisms enabling robustness against variability complement each other. Early STAT3 activation is robust as long as cytokine concentrations are low. Robustness at high cytokine concentrations is ensured by high STAT3 expression or serine phosphorylation. Later the feedback-inhibitor SOCS3 increases robustness. Channel Capacity of JAK/STAT signalling is limited by cell-to-cell variability in STAT3 expression and is affected by the same mechanisms governing robustness. Increasing STAT3 amount increases Channel Capacity and robustness, whereas increasing STAT3 tyrosine phosphorylation reduces robustness but increases Channel Capacity. In summary, we elucidate mechanisms preventing dysregulated signalling by enabling reliable JAK/STAT signalling despite cell-to-cell heterogeneity.