Single-cell variation leads to population invariance in NF-κB signaling dynamics

Co-imaging of RelA and c-Rel reveals features of NF-κB signaling for ligand discrimination

Individual cell sensing of external cues has evolved through the temporal patterns in signaling. Since nuclear factor κB (NF-κB) signaling dynamics have been examined using a single subunit, RelA, it remains unclear whether more information might be transmitted via other subunits. Using NF-κB double-knockin reporter mice, we monitored both canonical NF-κB subunits, RelA and c-Rel, simultaneously in single macrophages by quantitative live-cell imaging. We show that signaling features of RelA and c-Rel convey more information about the stimuli than those of either subunit alone. Machine learning is used to predict the ligand identity accurately based on RelA and c-Rel signaling features without considering the co-activated factors. Ligand discrimination is achieved through selective non-redundancy of RelA and c-Rel signaling dynamics, as well as their temporal coordination. These results suggest a potential role of c-Rel in fine-tuning immune responses and highlight the need for approaches that will elucidate the mechanisms regulating NF-κB subunit specificity.

Evidence of collective influence in innate sensing using fluidic force microscopy

The innate immune system initiates early response to infection by sensing molecular patterns of infection through pattern-recognition receptors (PRRs). Previous work on PRR stimulation of macrophages revealed significant heterogeneity in single cell responses, suggesting the importance of individual macrophage stimulation. Current methods either isolate individual macrophages or stimulate a whole culture and measure individual readouts. We probed single cell NF-κB responses to localized stimuli within a naïve culture with Fluidic Force Microscopy (FluidFM). Individual cells stimulated in naïve culture were more sensitive compared to individual cells in uniformly stimulated cultures. In cluster stimulation, NF-κB activation decreased with increased cell density or decreased stimulation time. Our results support the growing body of evidence for cell-to-cell communication in macrophage activation, and limit potential mechanisms. Such a mechanism might be manipulated to tune macrophage sensitivity, and the density-dependent modulation of sensitivity to PRR signals could have relevance to biological situations where macrophage density increases.

Myddosome clustering in IL-1 receptor signaling regulates the formation of an NF-kB activating signalosome

IL-1 receptor (IL-1R) signaling can activate thresholded invariant outputs and proportional outputs that scale with the amount of stimulation. Both responses require the Myddosome, a multiprotein complex. The Myddosome is required for polyubiquitin chain formation and NF-kB signaling. However, how these signals are spatially and temporally regulated to drive switch-like and proportional outcomes is not understood. During IL-1R signaling, Myddosomes dynamically reorganize into multi-Myddosome clusters at the cell membrane. Blockade of clustering using nanoscale extracellular barriers reduces NF-kB activation. Myddosomes function as scaffolds that assemble an NF-kB signalosome consisting of E3-ubiquitin ligases TRAF6 and LUBAC, K63/M1-linked polyubiquitin chains, phospho-IKK, and phospho-p65. This signalosome preferentially assembles at regions of high Myddosome density, which enhances the recruitment of TRAF6 and LUBAC. Extracellular barriers that restrict Myddosome clustering perturbed the recruitment of both ligases. We find that LUBAC was especially sensitive to clustering with 10-fold lower recruitment to single Myddosomes than clustered Myddosomes. These data reveal that the clustering behavior of Myddosomes provides a basis for digital and analog IL-1R signaling.

Droplet microarray platforms for high-throughput drug screening

High-throughput screening platforms are fundamental for the rapid and efficient processing of large amounts of experimental data. Parallelization and miniaturization of experiments are important for improving their cost-effectiveness. The development of miniaturized high-throughput screening platforms is essential in the fields of biotechnology, medicine, and pharmacology. Currently, most laboratories use 96- or 384-well microtiter plates for screening; however, they have disadvantages, such as high reagent and cell consumption, low throughput, and inability to avoid cross-contamination, which need to be further optimized. Droplet microarrays, as novel screening platforms, can effectively avoid these shortcomings. Here, the preparation method of the droplet microarray, method of adding compounds in parallel, and means to read the results are briefly described. Next, the latest research on droplet microarray platforms in biomedicine is presented, including their application in high-throughput culture, cell screening, high-throughput nucleic acid screening, drug development, and individualized medicine. Finally, the challenges and future trends in droplet microarray technology are summarized.

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.

Insights on the NF-κB System Using Live Cell Imaging: Recent Developments and Future Perspectives

The transcription factor family of nuclear factor kappa B (NF-κB) proteins is widely recognized as a key player in inflammation and the immune responses, where it plays a fundamental role in translating external inflammatory cues into precise transcriptional programs, including the timely expression of a wide variety of cytokines/chemokines. Live cell imaging in single cells showed approximately 15 years ago that the canonical activation of NF-κB upon stimulus is very dynamic, including oscillations of its nuclear localization with a period close to 1.5 hours. This observation has triggered a fruitful interdisciplinary research line that has provided novel insights on the NF-κB system: how its heterogeneous response differs between cell types but also within homogeneous populations; how NF-κB dynamics translate external cues into intracellular signals and how NF-κB dynamics affects gene expression. Here we review the main features of this live cell imaging approach to the study of NF-κB, highlighting the key findings, the existing gaps of knowledge and hinting towards some of the potential future steps of this thriving research field.

NfκB signaling dynamics and their target genes differ between mouse blood cell types and induce distinct cell behavior

Cells can use signaling pathway activity over time (i.e., dynamics) to control cell fates. However, little is known about the potential existence and function of signaling dynamics in primary hematopoietic stem and progenitor cells (HSPCs). Here, we use time-lapse imaging and tracking of single murine HSPCs from GFP-p65/H2BmCherry reporter mice to quantify their nuclear factor κB (NfκB) activity dynamics in response to TNFα and IL1β. We find response dynamics to be heterogeneous between individual cells, with cell type specific dynamics distributions. Transcriptome sequencing of single cells physically isolated after live dynamics quantification shows activation of different target gene programs in cells with different dynamics. Finally, artificial induction of oscillatory NfκB activity causes changes in GMP behavior. Thus, HSPC behavior can be influenced by signaling dynamics, which are tightly regulated during hematopoietic differentiation and enable cell type specific responses to the same signaling inputs.

Stimulus-specific responses in innate immunity: Multilayered regulatory circuits

Immune sentinel cells initiate immune responses to pathogens and tissue injury and are capable of producing highly stimulus-specific responses. Insight into the mechanisms underlying such specificity has come from the identification of regulatory factors and biochemical pathways, as well as the definition of signaling circuits that enable combinatorial and temporal coding of information. Here, we review the multi-layered molecular mechanisms that underlie stimulus-specific gene expression in macrophages. We categorize components of inflammatory and anti-pathogenic signaling pathways into five layers of regulatory control and discuss unifying mechanisms determining signaling characteristics at each layer. In this context, we review mechanisms that enable combinatorial and temporal encoding of information, identify recurring regulatory motifs and principles, and present strategies for integrating experimental and computational approaches toward the understanding of signaling specificity in innate immunity.

A forecast for large-scale, predictive biology: Lessons from meteorology

Quantitative systems biology, in which predictive mathematical models are constructed to guide the design of experiments and predict experimental outcomes, is at an exciting transition point, where the foundational scientific principles are becoming established, but the impact is not yet global. The next steps necessary for mathematical modeling to transform biological research and applications, in the same way it has already transformed other fields, is not completely clear. The purpose of this perspective is to forecast possible answers to this question-what needs to happen next-by drawing on the experience gained in another field, specifically meteorology. We review here a number of lessons learned in weather prediction that are directly relevant to biological systems modeling, and that we believe can enable the same kinds of global impact in our field as atmospheric modeling makes today.

Six distinct NFκB signaling codons convey discrete information to distinguish stimuli and enable appropriate macrophage responses

Macrophages initiate inflammatory responses via the transcription factor NFκB. The temporal pattern of NFκB activity determines which genes are expressed and thus, the type of response that ensues. Here, we examined how information about the stimulus is encoded in the dynamics of NFκB activity. We generated an mVenus-RelA reporter mouse line to enable high-throughput live-cell analysis of primary macrophages responding to host- and pathogen-derived stimuli. An information-theoretic workflow identified six dynamical features-termed signaling codons-that convey stimulus information to the nucleus. In particular, oscillatory trajectories were a hallmark of responses to cytokine but not pathogen-derived stimuli. Single-cell imaging and RNA sequencing of macrophages from a mouse model of Sjögren's syndrome revealed inappropriate responses to stimuli, suggestive of confusion of two NFκB signaling codons. Thus, the dynamics of NFκB signaling classify immune threats through six signaling codons, and signal confusion based on defective codon deployment may underlie the etiology of some inflammatory diseases.

Epigenetically regulated digital signaling defines epithelial innate immunity at the tissue level

To prevent damage to the host or its commensal microbiota, epithelial tissues must match the intensity of the immune response to the severity of a biological threat. Toll-like receptors allow epithelial cells to identify microbe associated molecular patterns. However, the mechanisms that mitigate biological noise in single cells to ensure quantitatively appropriate responses remain unclear. Here we address this question using single cell and single molecule approaches in mammary epithelial cells and primary organoids. We find that epithelial tissues respond to bacterial microbe associated molecular patterns by activating a subset of cells in an all-or-nothing (i.e. digital) manner. The maximum fraction of responsive cells is regulated by a bimodal epigenetic switch that licenses the TLR2 promoter for transcription across multiple generations. This mechanism confers a flexible memory of inflammatory events as well as unique spatio-temporal control of epithelial tissue-level immune responses. We propose that epigenetic licensing in individual cells allows for long-term, quantitative fine-tuning of population-level responses.

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.

Microfluidic platform enables live-cell imaging of signaling and transcription combined with multiplexed secretion measurements in the same single cells

Innate immune cells, including macrophages and dendritic cells, protect the host from pathogenic assaults in part through secretion of a program of cytokines and chemokines (C/Cs). Cell-to-cell variability in C/C secretion appears to contribute to the regulation of the immune response, but the sources of secretion variability are largely unknown. To begin to track the biological sources that control secretion variability, we developed and validated a microfluidic device to integrate live-cell imaging of fluorescent reporter proteins with a single-cell assay of protein secretion. We used this device to image NF-κB RelA nuclear translocation dynamics and Tnf transcription dynamics in macrophages in response to stimulation with the bacterial component lipopolysaccharide (LPS), followed by quantification of secretion of TNF, CCL2, CCL3, and CCL5. We found that the timing of the initial peak of RelA signaling in part determined the relative level of TNF and CCL3 secretion, but not CCL2 and CCL5 secretion. Our results support evidence that differences in timing across cell processes partly account for cell-to-cell variability in downstream responses, but that other factors introduce variability at each biological step.

Macrophages employ quorum licensing to regulate collective activation

Macrophage-initiated inflammation is tightly regulated to eliminate threats such as infections while suppressing harmful immune activation. However, individual cells’ signaling responses to pro-inflammatory cues are heterogeneous, with subpopulations emerging with high or low activation states. Here, we use single-cell tracking and dynamical modeling to develop and validate a revised model for lipopolysaccharide (LPS)-induced macrophage activation that invokes a mechanism we term quorum licensing. The results show that bimodal phenotypic partitioning of macrophages is primed during the resting state, dependent on cumulative history of cell density, predicted by extrinsic noise in transcription factor expression, and independent of canonical LPS-induced intercellular feedback in the tumor necrosis factor (TNF) response. Our analysis shows how this density-dependent coupling produces a nonlinear effect on collective TNF production. We speculate that by linking macrophage density to activation, this mechanism could amplify local responses to threats and prevent false alarms.

Macrophage activation is tightly regulated to maintain immune homeostasis, yet activation is also heterogeneous. Here, the authors show that macrophages coordinate activation by partitioning into two phenotypes that can nonlinearly amplify collective inflammatory cytokine production as a function of cell density.

Dual Feedforward Loops Modulate Type I Interferon Responses and Induce Selective Gene Expression during TLR4 Activation

Although the dynamic features of type I coherent feedforward loop (C1-FFL) has been well studied, how C1-FFL shapes cell-to-cell heterogeneity remains unclear. Here, we found that C1-FFL with OR logic serves as "noise reducer," whereas C1-FFL with AND logic functions as "noise propagator" to fine-tune the heterogeneity of signaling molecule's activation. Within Toll-like receptor 4 (TLR4) signaling pathway, we demonstrated that MyD88 together with TRIF generates a C1-FFL to control TBK1 phosphorylation and reduce its cell-to-cell heterogeneity, whereas noisy TRIF activation induced high heterogeneity of IRF3 activation through another C1-FFL. We further developed a mathematical model with dual C1-FFLs to uncover how MyD88 and TRIF encoded differential dynamics for TBK1 and IRF3 activation. Integration of dual FFLs drives MyD88-TBK1 axis to determine the specificity of IFN-stimulated genes transcription. Collectively, our work elucidates a paradigm that tunable TLR4-mediated type I IFN responses are subtly controlled by dual FFLs.

What Will B Will B: Identifying Molecular Determinants of Diverse B-Cell Fate Decisions Through Systems Biology

B-cells are the poster child for cellular diversity and heterogeneity. The diverse repertoire of B lymphocytes, each expressing unique antigen receptors, provides broad protection against pathogens. However, B-cell diversity goes beyond unique antigen receptors. Side-stepping B-cell receptor (BCR) diversity through BCR-independent stimuli or engineered organisms with monoclonal BCRs still results in seemingly identical B-cells reaching a wide variety of fates in response to the same challenge. Identifying to what extent the molecular state of a B-cell determines its fate is key to gaining a predictive understanding of B-cells and consequently the ability to control them with targeted therapies. Signals received by B-cells through transmembrane receptors converge on intracellular molecular signaling networks, which control whether each B-cell divides, dies, or differentiates into a number of antibody-secreting distinct B-cell subtypes. The signaling networks that interpret these signals are well known to be susceptible to molecular variability and noise, providing a potential source of diversity in cell fate decisions. Iterative mathematical modeling and experimental studies have provided quantitative insight into how B-cells achieve distinct fates in response to pathogenic stimuli. Here, we review how systems biology modeling of B-cells, and the molecular signaling networks controlling their fates, is revealing the key determinants of cell-to-cell variability in B-cell destiny.

Order by chance: origins and benefits of stochasticity in immune cell fate control

To protect against diverse challenges, the immune system must continuously generate an arsenal of specialized cell types, each of which can mount a myriad of effector responses upon detection of potential threats. To do so, it must generate multiple differentiated cell populations with defined sizes and proportions, often from rare starting precursor cells. Here, we discuss the emerging view that inherently probabilistic mechanisms, involving rare, rate-limiting regulatory events in single cells, control fate decisions and population sizes and fractions during immune development and function. We first review growing evidence that key fate control points are gated by stochastic signaling and gene regulatory events that occur infrequently over decision-making timescales, such that initially homogeneous cells can adopt variable outcomes in response to uniform signals. We next discuss how such stochastic control can provide functional capabilities that are harder to achieve with deterministic control strategies, and may be central to robust immune system function.

Stimulus-specificity in the Responses of Immune Sentinel Cells

Innate immune sentinel cells must initiate and orchestrate appropriate immune responses for myriad pathogens. These stimulus-specific gene expression responses are mediated by combinatorial and temporal coding within a handful of immune response signaling pathways. We outline the scope of our current understanding and indicate pressing outstanding questions. The innate immune response is a first-line defense against invading pathogens and coordinates the activation and recruitment of specialized immune cells, thereby initiating the adaptive immune response. While the adaptive immune system is capable of highly pathogen-specific immunity through the process of genetic recombination and clonal selection, innate immunity is frequently viewed as a catch-all system that initiates general immune activation. In this review, we are re-examining this view, as we are distinguishing between immune sentinel functions mediated by macrophages and dendritic cells and innate immune effector functions mediated by cells such as neutrophils, NK cells, etc. Given pathogen diversity, including modes of entry, replication cycles, and strategies of immune evasion and spread, all successive waves of the immune response ought to be tailored to the specific immune threat, leading us to postulate that immune sentinel functions by macrophages and dendritic cells ought to be highly stimulus-specific. Here we review the experimental evidence for stimulus-specific responses by immune sentinel cells which initiate and coordinate immune responses, as well as the mechanisms by which this specificity may be achieved.

An NFκB Activity Calculator to Delineate Signaling Crosstalk: Type I and II Interferons Enhance NFκB via Distinct Mechanisms

Nuclear factor kappa B (NFκB) is a transcription factor that controls inflammation and cell survival. In clinical histology, elevated NFκB activity is a hallmark of poor prognosis in inflammatory disease and cancer, and may be the result of a combination of diverse micro-environmental constituents. While previous quantitative studies of NFκB focused on its signaling dynamics in single cells, we address here how multiple stimuli may combine to control tissue level NFκB activity. We present a novel, simplified model of NFκB (SiMoN) that functions as an NFκB activity calculator. We demonstrate its utility by exploring how type I and type II interferons modulate NFκB activity in macrophages. Whereas, type I IFNs potentiate NFκB activity by inhibiting translation of IκBα and by elevating viral RNA sensor (RIG-I) expression, type II IFN amplifies NFκB activity by increasing the degradation of free IκB through transcriptional induction of proteasomal cap components (PA28). Both cross-regulatory mechanisms amplify NFκB activation in response to weaker (viral) inducers, while responses to stronger (bacterial or cytokine) inducers remain largely unaffected. Our work demonstrates how the NFκB calculator can reveal distinct mechanisms of crosstalk on NFκB activity in interferon-containing microenvironments.

Controlling Nuclear NF-κB Dynamics by β-TrCP-Insights from a Computational Model

The canonical nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling pathway regulates central processes in mammalian cells and plays a fundamental role in the regulation of inflammation and immunity. Aberrant regulation of the activation of the transcription factor NF-κB is associated with severe diseases such as inflammatory bowel disease and arthritis. In the canonical pathway, the inhibitor IκB suppresses NF-κB’s transcriptional activity. NF-κB becomes active upon the degradation of IκB, a process that is, in turn, regulated by the β-transducin repeat-containing protein (β-TrCP). β-TrCP has therefore been proposed as a promising pharmacological target in the development of novel therapeutic approaches to control NF-κB’s activity in diseases. This study explores the extent to which β-TrCP affects the dynamics of nuclear NF-κB using a computational model of canonical NF-κB signaling. The analysis predicts that β-TrCP influences the steady-state concentration of nuclear NF-κB, as well as changes characteristic dynamic properties of nuclear NF-κB, such as fold-change and the duration of its response to pathway stimulation. The results suggest that the modulation of β-TrCP has a high potential to regulate the transcriptional activity of NF-κB.

NF-κB signaling dynamics is controlled by a dose-sensing autoregulatory loop

Over the last decade, multiple studies have shown that signaling proteins activated in different temporal patterns, such as oscillatory, transient, and sustained, can result in distinct gene expression patterns or cell fates. However, the molecular events that ensure appropriate stimulus- and dose-dependent dynamics are not often understood and are difficult to investigate. Here, we used single-cell analysis to dissect the mechanisms underlying the stimulus- and dose-encoding patterns in the innate immune signaling network. We found that Toll-like receptor (TLR) and interleukin-1 receptor (IL-1R) signaling dynamics relied on a dose-dependent, autoinhibitory loop that rendered cells refractory to further stimulation. Using inducible gene expression and optogenetics to perturb the network at different levels, we identified IL-1R-associated kinase 1 (IRAK1) as the dose-sensing node responsible for limiting signal flow during the innate immune response. Although the kinase activity of IRAK1 was not required for signal propagation, it played a critical role in inhibiting the nucleocytoplasmic oscillations of the transcription factor NF-κB. Thus, protein activities that may be "dispensable" from a topological perspective can nevertheless be essential in shaping the dynamic response to the external environment.

NF-κB Signaling in Macrophages: Dynamics, Crosstalk, and Signal Integration

The nuclear factor-κB (NF-κB) signaling pathway is one of the best understood immune-related pathways thanks to almost four decades of intense research. NF-κB signaling is activated by numerous discrete stimuli and is a master regulator of the inflammatory response to pathogens and cancerous cells, as well as a key regulator of autoimmune diseases. In this regard, the role of NF-κB signaling in immunity is not unlike that of the macrophage. The dynamics by which NF-κB proteins shuttle between the cytoplasm and the nucleus to initiate transcription have been studied rigorously in fibroblasts and other non-hematopoietic cells, but many questions remain as to how current models of NF-κB signaling and dynamics can be translated to innate immune cells such as macrophages. In this review, we will present recent research on the dynamics of NF-κB signaling and focus especially on how these dynamics vary in different cell types, while discussing why these characteristics may be important. We will end by looking ahead to how new techniques and technologies should allow us to analyze these signaling processes with greater clarity, bringing us closer to a more complete understanding of inflammatory transcription factor dynamics and how different cellular contexts might allow for appropriate control of innate immune responses.

A Two-Pulse Cellular Stimulation Test Elucidates Variability and Mechanisms in Signaling Pathways

Mammalian cells respond in a variable manner when provided with physiological pulses of ligand, such as low concentrations of acetylcholine present for just tens of seconds or TNFα for just tens of minutes. For a two-pulse stimulation, some cells respond to both pulses, some do not respond, and yet others respond to only one or the other pulse. Are these different response patterns the result of the small number of ligands being able to only stochastically activate the pathway at random times or an output pattern from a deterministic algorithm responding differently to different stimulation intervals? If the response is deterministic in nature, what parameters determine whether a response is generated or skipped? To answer these questions, we developed a two-pulse test that utilizes different rest periods between stimulation pulses. This "rest-period test" revealed that cells skip responses predictably as the rest period is shortened. By combining these experimental results with a mathematical model of the pathway, we further obtained mechanistic insight into potential sources of response variability. Our analysis indicates that in both intracellular calcium and NFκB signaling, response variability is consistent with extrinsic noise (cell-to-cell variability in protein levels), a short-term memory of stimulation, and high Hill coefficient processes. Furthermore, these results support recent works that have emphasized the role of deterministic processes for explaining apparently stochastic cellular response variability and indicate that even weak stimulations likely guide mammalian cells to appropriate fates rather than leaving outcomes to chance. We envision that the rest-period test can be applied to other signaling pathways to extract mechanistic insight.

Mechanisms for the epigenetic inheritance of stress response in single cells

Cells have evolved to dynamically respond to different types of environmental and physiological stress conditions. The information about a previous stress stimulus experience by a mother cell can be passed to its descendants, allowing them to better adapt to and survive in new environments. In recent years, live-cell imaging combined with cell-lineage tracking approaches has elucidated many important principles that guide stress inheritance at the single-cell and population level. In this review, we summarize different strategies that cells can employ to pass the 'memory' of previous stress responses to their descendants. Among these strategies, we focus on a recent discovery of how specific features of Msn2 nucleo-cytoplasmic shuttling dynamics could be inherited across cell lineages. We also discuss how stress response can be transmitted to progenies through changes in chromatin and through partitioning of anti-stress factors and/or damaged macromolecules between mother and daughter cells during cell division. Finally, we highlight how emergent technologies will help address open questions in the field.

Combinatorial processing of bacterial and host-derived innate immune stimuli at the single-cell level

During the course of a bacterial infection, cells are exposed simultaneously to a range of bacterial and host factors, which converge on the central transcription factor nuclear factor (NF)-κB. How do single cells integrate and process these converging stimuli? Here we tackle the question of how cells process combinatorial signals by making quantitative single-cell measurements of the NF-κB response to combinations of bacterial lipopolysaccharide and the stress cytokine tumor necrosis factor. We found that cells encode the presence of both stimuli via the dynamics of NF-κB nuclear translocation in individual cells, suggesting the integration of NF-κB activity for these stimuli occurs at the molecular and pathway level. However, the gene expression and cytokine secretion response to combinatorial stimuli were more complex, suggesting that other factors in addition to NF-κB contribute to signal integration at downstream layers of the response. Taken together, our results support the theory that during innate immune threat assessment, a pathogen recognized as both foreign and harmful will recruit an enhanced immune response. Our work highlights the remarkable capacity of individual cells to process multiple input signals and suggests that a deeper understanding of signal integration mechanisms will facilitate efforts to control dysregulated immune responses.

A Bayesian approach to modelling heterogeneous calcium responses in cell populations

Calcium responses have been observed as spikes of the whole-cell calcium concentration in numerous cell types and are essential for translating extracellular stimuli into cellular responses. While there are several suggestions for how this encoding is achieved, we still lack a comprehensive theory. To achieve this goal it is necessary to reliably predict the temporal evolution of calcium spike sequences for a given stimulus. Here, we propose a modelling framework that allows us to quantitatively describe the timing of calcium spikes. Using a Bayesian approach, we show that Gaussian processes model calcium spike rates with high fidelity and perform better than standard tools such as peri-stimulus time histograms and kernel smoothing. We employ our modelling concept to analyse calcium spike sequences from dynamically-stimulated HEK293T cells. Under these conditions, different cells often experience diverse stimulus time courses, which is a situation likely to occur in vivo. This single cell variability and the concomitant small number of calcium spikes per cell pose a significant modelling challenge, but we demonstrate that Gaussian processes can successfully describe calcium spike rates in these circumstances. Our results therefore pave the way towards a statistical description of heterogeneous calcium oscillations in a dynamic environment.

Upon stimulation a large number of cell types respond with transient increases of the intracellular calcium concentration, which often take the form of repetitive spikes. It is therefore believed that calcium spikes play a central role in cellular signal transduction. A critical feature of these calcium spikes is that they occur randomly, which raises the question of how we can predict the timing of calcium spikes. We here show that by using Bayesian ideas and concepts from stochastic processes, we can quantitatively compute the calcium spike rate for a given stimulus. Our analysis also demonstrates that traditional methods for spike rate estimation perform less favourably compared to a Bayesian approach when small numbers of cells are investigated. To test our methodology under conditions that closely mimic those experienced in vivo we challenged cells with agonist concentrations that vary both in space and time. We find that cells that experience similar stimulus profiles are described by similar calcium spike rates. This suggests that calcium spike rates may constitute a quantitative description of whole-cell calcium spiking that reflects both the randomness and the spatiotemporal organisation of the calcium signalling machinery.

Measuring Signaling and RNA-Seq in the Same Cell Links Gene Expression to Dynamic Patterns of NF-κB Activation

Signaling proteins display remarkable cell-to-cell heterogeneity in their dynamic responses to stimuli, but the consequences of this heterogeneity remain largely unknown. For instance, the contribution of the dynamics of the innate immune transcription factor nuclear factor κB (NF-κB) to gene expression output is disputed. Here we explore these questions by integrating live-cell imaging approaches with single-cell sequencing technologies. We used this approach to measure both the dynamics of lipopolysaccharide-induced NF-κB activation and the global transcriptional response in the same individual cell. Our results identify multiple, distinct cytokine expression patterns that are correlated with NF-κB activation dynamics, establishing a functional role for NF-κB dynamics in determining cellular phenotypes. Applications of this approach to other model systems and single-cell sequencing technologies have significant potential for discovery, as it is now possible to trace cellular behavior from the initial stimulus, through the signaling pathways, down to genome-wide changes in gene expression, all inside of a single cell.

New insights into mammalian signaling pathways using microfluidic pulsatile inputs and mathematical modeling

Temporally modulated input mimics physiology. This chemical communication strategy filters the biochemical noise through entrainment and phase-locking. Under laboratory conditions, it also expands the observability space for downstream responses. A combined approach involving microfluidic pulsatile stimulation and mathematical modeling has led to deciphering of hidden/unknown temporal motifs in several mammalian signaling pathways and has provided mechanistic insights, including how these motifs combine to form distinct band-pass filters and govern fate regulation under dynamic microenvironment. This approach can be utilized to understand signaling circuit architectures and to gain mechanistic insights for several other signaling systems. Potential applications include synthetic biology and biotechnology, in developing pharmaceutical interventions, and in developing lab-on-chip models.

Droplet Microarray Based on Superhydrophobic-Superhydrophilic Patterns for Single Cell Analysis

Single-cell analysis provides fundamental information on individual cell response to different environmental cues and is a growing interest in cancer and stem cell research. However, current existing methods are still facing challenges in performing such analysis in a high-throughput manner whilst being cost-effective. Here we established the Droplet Microarray (DMA) as a miniaturized screening platform for high-throughput single-cell analysis. Using the method of limited dilution and varying cell density and seeding time, we optimized the distribution of single cells on the DMA. We established culturing conditions for single cells in individual droplets on DMA obtaining the survival of nearly 100% of single cells and doubling time of single cells comparable with that of cells cultured in bulk cell population using conventional methods. Our results demonstrate that the DMA is a suitable platform for single-cell analysis, which carries a number of advantages compared with existing technologies allowing for treatment, staining and spot-to-spot analysis of single cells over time using conventional analysis methods such as microscopy.

Mapping paths: new approaches to dissect eukaryotic signaling circuitry

Eukaryotic cells are precisely “wired” to coordinate changes in external and intracellular signals with corresponding adjustments in the output of complex and often interconnected signaling pathways. These pathways are critical in understanding cellular growth and function, and several experimental trends are emerging with applicability toward more fully describing the composition and topology of eukaryotic signaling networks. In particular, recent studies have implemented CRISPR/Cas-based screens in mouse and human cell lines for genes involved in various cell growth and disease phenotypes. Proteomic methods using mass spectrometry have enabled quantitative and dynamic profiling of protein interactions, revealing previously undiscovered complexes and allele-specific protein interactions. Methods for the single-cell study of protein localization and gene expression have been integrated with computational analyses to provide insight into cell signaling in yeast and metazoans. In this review, we present an overview of exemplary studies using the above approaches, relevant for the analysis of cell signaling and indeed, more broadly, for many modern biological applications.

Signal transduction controls heterogeneous NF-κB dynamics and target gene expression through cytokine-specific refractory states

Cells respond dynamically to pulsatile cytokine stimulation. Here we report that single, or well-spaced pulses of TNFα (>100 min apart) give a high probability of NF-κB activation. However, fewer cells respond to shorter pulse intervals (<100 min) suggesting a heterogeneous refractory state. This refractory state is established in the signal transduction network downstream of TNFR and upstream of IKK, and depends on the level of the NF-κB system negative feedback protein A20. If a second pulse within the refractory phase is IL-1β instead of TNFα, all of the cells respond. This suggests a mechanism by which two cytokines can synergistically activate an inflammatory response. Gene expression analyses show strong correlation between the cellular dynamic response and NF-κB-dependent target gene activation. These data suggest that refractory states in the NF-κB system constitute an inherent design motif of the inflammatory response and we suggest that this may avoid harmful homogenous cellular activation.

Dynamic NF-κB and E2F interactions control the priority and timing of inflammatory signalling and cell proliferation

Dynamic cellular systems reprogram gene expression to ensure appropriate cellular fate responses to specific extracellular cues. Here we demonstrate that the dynamics of Nuclear Factor kappa B (NF-κB) signalling and the cell cycle are prioritised differently depending on the timing of an inflammatory signal. Using iterative experimental and computational analyses, we show physical and functional interactions between NF-κB and the E2 Factor 1 (E2F-1) and E2 Factor 4 (E2F-4) cell cycle regulators. These interactions modulate the NF-κB response. In S-phase, the NF-κB response was delayed or repressed, while cell cycle progression was unimpeded. By contrast, activation of NF-κB at the G1/S boundary resulted in a longer cell cycle and more synchronous initial NF-κB responses between cells. These data identify new mechanisms by which the cellular response to stress is differentially controlled at different stages of the cell cycle.

DOI:http://dx.doi.org/10.7554/eLife.10473.001

Investigating how cells adapt to the constantly changing environment inside the body is vitally important for understanding how the body responds to an injury or infection. One of the ways in which human cells adapt is by dividing to produce new cells. This takes place in a repeating pattern of events, known as the cell cycle, through which a cell copies its DNA (in a stage known as S-phase) and then divides to make two daughter cells. Each stage of the cell cycle is tightly controlled; for example, a family of proteins called E2 factors control the entry of the cell into S phase.

“Inflammatory” signals produced by a wound or during an infection can activate a protein called Nuclear Factor-kappaB (NF-κB), which controls the activity of genes that allow cells to adapt to the situation. Research shows that the activity of NF-κB is also regulated by the cell cycle, but it has not been clear how this works. Here, Ankers et al. investigated whether the stage of the cell cycle might affect how NF-κB responds to inflammatory signals.

The experiments show that the NF-κB response was stronger in cells that were just about to enter S-phase than in cells that were already copying their DNA. An E2 factor called E2F-1 –which accumulates in the run up to S-phase – interacts with NF-κB and can alter the activity of certain genes. However, during S-phase, another E2 factor family member called E2F-4 binds to NF-κB and represses its activation. Next, Ankers et al. used a mathematical model to understand how these protein interactions can affect the response of cells to inflammatory signals.

These findings suggest that direct interactions between E2 factor proteins and NF-κB enable cells to decide whether to divide or react in different ways to inflammatory signals. The research tools developed in this study, combined with other new experimental techniques, will allow researchers to accurately predict how cells will respond to inflammatory signals at different points in the cell cycle.

DOI:http://dx.doi.org/10.7554/eLife.10473.002

Redefining Signaling Pathways with an Expanding Single-Cell Toolbox

Genetically identical cells respond heterogeneously to uniform environmental stimuli. Consequently, investigating the signaling networks that control these cell responses using 'average' bulk cell measurements can obscure underlying mechanisms and misses information emerging from cell-to-cell variability. Here we review recent technological advances including live-cell fluorescence imaging-based approaches and microfluidic devices that enable measurements of signaling networks, dynamics, and responses in single cells. We discuss how these single-cell tools have uncovered novel mechanistic insights for canonical signaling pathways that control cell proliferation (ERK), DNA-damage responses (p53), and innate immune and stress responses (NF-κB). Future improvements in throughput and multiplexing, analytical pipelines, and in vivo applicability will all significantly expand the biological information gained from single-cell measurements of signaling pathways.

Digital signaling decouples activation probability and population heterogeneity

Digital signaling enhances robustness of cellular decisions in noisy environments, but it is unclear how digital systems transmit temporal information about a stimulus. To understand how temporal input information is encoded and decoded by the NF-κB system, we studied transcription factor dynamics and gene regulation under dose- and duration-modulated inflammatory inputs. Mathematical modeling predicted and microfluidic single-cell experiments confirmed that integral of the stimulus (or area, concentration × duration) controls the fraction of cells that activate NF-κB in the population. However, stimulus temporal profile determined NF-κB dynamics, cell-to-cell variability, and gene expression phenotype. A sustained, weak stimulation lead to heterogeneous activation and delayed timing that is transmitted to gene expression. In contrast, a transient, strong stimulus with the same area caused rapid and uniform dynamics. These results show that digital NF-κB signaling enables multidimensional control of cellular phenotype via input profile, allowing parallel and independent control of single-cell activation probability and population heterogeneity.

Paracrine communication maximizes cellular response fidelity in wound signaling

Population averaging due to paracrine communication can arbitrarily reduce cellular response variability. Yet, variability is ubiquitously observed, suggesting limits to paracrine averaging. It remains unclear whether and how biological systems may be affected by such limits of paracrine signaling. To address this question, we quantify the signal and noise of Ca(2+) and ERK spatial gradients in response to an in vitro wound within a novel microfluidics-based device. We find that while paracrine communication reduces gradient noise, it also reduces the gradient magnitude. Accordingly we predict the existence of a maximum gradient signal to noise ratio. Direct in vitro measurement of paracrine communication verifies these predictions and reveals that cells utilize optimal levels of paracrine signaling to maximize the accuracy of gradient-based positional information. Our results demonstrate the limits of population averaging and show the inherent tradeoff in utilizing paracrine communication to regulate cellular response fidelity.

Distinct single-cell signaling characteristics are conferred by the MyD88 and TRIF pathways during TLR4 activation

Toll-like receptors (TLRs) recognize specific pathogen-associated molecular patterns and initiate innate immune responses through signaling pathways that depend on the adaptor proteins MyD88 (myeloid differentiation marker 88) or TRIF (TIR domain-containing adaptor protein-inducing interferon-β). TLR4, in particular, uses both adaptor proteins to activate the transcription factor nuclear factor κB (NF-κB); however, the specificity and redundancy of these two pathways remain to be elucidated. We developed a mathematical model to show how each pathway encodes distinct dynamical features of NF-κB activity and makes distinct contributions to the high variability observed in single-cell measurements. The assembly of a macromolecular signaling platform around MyD88 associated with receptors at the cell surface determined the timing of initial responses to generate a reliable, digital NF-κB signal. In contrast, ligand-induced receptor internalization into endosomes produced noisy, delayed, yet sustained NF-κB signals through TRIF. With iterative mathematical model development, we predicted the molecular mechanisms by which the MyD88- and TRIF-mediated pathways provide ligand concentration-dependent signaling dynamics that transmit information about the pathogen threat.

Automated co-culture system for spatiotemporal analysis of cell-to-cell communication

We present a microfluidic co-culture system that generates localized and precisely formulated immune signals among a population of cells, enabling spatiotemporal analysis of paracrine signal transmission between different cell types. The automated system allows us to create temporally modulated chemical inputs that can be delivered to single signal-transmitting and receiving cells in a highly controlled way. Using this system we stimulated a single macrophage with brief pulses of bacterial LPS and observed the macrophage transmitted TNF signal propagating in a population of fibroblasts via NF-κB activation. The signal receiving fibroblasts transformed the TNF signal into a spatiotemporally distributed NF-κB output, recapitulating the initiation of immune response to bacterial infection.