Quantifying dynamic pro-inflammatory gene expression and heterogeneity in single macrophage cells.

Macrophages must respond appropriately to pathogens and other pro-inflammatory stimuli in order to perform their roles in fighting infection. One way in which inflammatory stimuli can vary is in their dynamics-that is, the amplitude and duration of stimulus experienced by the cell. In this study, we performed long-term live cell imaging in a microfluidic device to investigate how the pro-inflammatory genes IRF1, CXCL10, and CXCL9 respond to dynamic interferon-gamma (IFNγ) stimulation. We found that IRF1 responds to low concentration or short duration IFNγ stimulation, whereas CXCL10 and CXCL9 require longer or higherconcentration stimulation to be expressed. We also investigated the heterogeneity in the expression of each gene and found that CXCL10 and CXCL9 have substantial cell-to-cell variability. In particular, the expression of CXCL10 appears to be largely stochastic with a subpopulation of nonresponding cells across all the stimulation conditions tested. We developed both deterministic and stochastic models for the expression of each gene. Our modeling analysis revealed that the heterogeneity in CXCL10 can be attributed to a slow chromatin-opening step that is on a similar timescale to that of adaptation of the upstream signal. In this way, CXCL10 expression in individual cells can remain stochastic in response to each pulse of repeated stimulation, which we also validated by experiments. Together, we conclude that pro-inflammatory genes in the same signaling pathway can respond to dynamic IFNγ stimulus with very different response features and that upstream signal adaptation can contribute to shaping heterogeneous gene expression.

Enhanced cellular longevity arising from environmental fluctuations.

Cellular longevity is regulated by both genetic and environmental factors. However, the interactions of these factors in the context of aging remain largely unclear. Here, we formulate a mathematical model for dynamic glucose modulation of a core gene circuit in yeast aging, which not only guided the design of pro-longevity interventions, but also revealed the theoretical principles underlying these interventions. We introduce the dynamical systems theory to capture two general means for promoting longevity - the creation of a stable fixed point in the "healthy" state of the cell and the dynamic stabilization of the system around this healthy state through environmental oscillations. Guided by the model, we investigate how both of these can be experimentally realized by dynamically modulating environmental glucose levels. The results establish a paradigm for theoretically analyzing the trajectories and perturbations of aging that can be generalized to aging processes in diverse cell types and organisms.

Engineering longevity-design of a synthetic gene oscillator to slow cellular aging.

Synthetic biology enables the design of gene networks to confer specific biological functions, yet it remains a challenge to rationally engineer a biological trait as complex as longevity. A naturally occurring toggle switch underlies fate decisions toward either nucleolar or mitochondrial decline during the aging of yeast cells. We rewired this endogenous toggle to engineer an autonomous genetic clock that generates sustained oscillations between the nucleolar and mitochondrial aging processes in individual cells. These oscillations increased cellular life span through the delay of the commitment to aging that resulted from either the loss of chromatin silencing or the depletion of heme. Our results establish a connection between gene network architecture and cellular longevity that could lead to rationally designed gene circuits that slow aging.

Statistical theory of asymmetric damage segregation in clonal cell populations.

Asymmetric damage segregation (ADS) is ubiquitous among unicellular organisms: After a mother cell divides, its two daughter cells receive sometimes slightly, sometimes strongly different fractions of damaged proteins accumulated in the mother cell. Previous studies demonstrated that ADS provides a selective advantage over symmetrically dividing cells by rejuvenating and perpetuating the population as a whole. In this work we focus on the statistical properties of damage in individual lineages and the overall damage distributions in growing populations for a variety of ADS models with different rules governing damage accumulation, segregation, and the lifetime dependence on damage. We show that for a large class of deterministic ADS rules the trajectories of damage along the lineages are chaotic, and the distributions of damage in cells born at a given time asymptotically becomes fractal. By exploiting the analogy of linear ADS models with the Iterated Function Systems known in chaos theory, we derive the Frobenius-Perron equation for the stationary damage density distribution and analytically compute the damage distribution moments and fractal dimensions. We also investigate nonlinear and stochastic variants of ADS models and show the robustness of the salient features of the damage distributions.

Statistical Theory of Asymmetric Damage Segregation in Clonal Cell Populations.

Asymmetric damage segregation (ADS) is ubiquitous among unicellular organisms: After a mother cell divides, its two daughter cells receive sometimes slightly, sometimes strongly different fractions of damaged proteins accumulated in the mother cell. Previous studies demonstrated that ADS provides a selective advantage over symmetrically dividing cells by rejuvenating and perpetuating the population as a whole. In this work we focus on the statistical properties of damage in individual lineages and the overall damage distributions in growing populations for a variety of ADS models with different rules governing damage accumulation, segregation, and the lifetime dependence on damage. We show that for a large class of deterministic ADS rules the trajectories of damage along the lineages are chaotic, and the distributions of damage in cells born at a given time asymptotically becomes fractal. By exploiting the analogy of linear ADS models with the Iterated Function Systems known in chaos theory, we derive the Frobenius-Perron equation for the stationary damage density distribution and analytically compute the damage distribution moments and fractal dimensions. We also investigate nonlinear and stochastic variants of ADS models and show the robustness of the salient features of the damage distributions.

Age-dependent aggregation of ribosomal RNA-binding proteins links deterioration in chromatin stability with challenges to proteostasis.

Chromatin instability and protein homeostasis (proteostasis) stress are two well-established hallmarks of aging, which have been considered largely independent of each other. Using microfluidics and single-cell imaging approaches, we observed that, during the replicative aging of Saccharomyces cerevisiae, a challenge to proteostasis occurs specifically in the fraction of cells with decreased stability within the ribosomal DNA (rDNA). A screen of 170 yeast RNA-binding proteins identified ribosomal RNA (rRNA)-binding proteins as the most enriched group that aggregate upon a decrease in rDNA stability induced by inhibition of a conserved lysine deacetylase Sir2. Further, loss of rDNA stability induces age-dependent aggregation of rRNA-binding proteins through aberrant overproduction of rRNAs. These aggregates contribute to age-induced proteostasis decline and limit cellular lifespan. Our findings reveal a mechanism underlying the interconnection between chromatin instability and proteostasis stress and highlight the importance of cell-to-cell variability in aging processes.

Optimal transcriptional regulation of dynamic bacterial responses to sudden drug exposures.

Cellular responses to the presence of toxic compounds in their environment require prompt expression of the correct levels of the appropriate enzymes, which are typically regulated by transcription factors that control gene expression for the duration of the response. The characteristics of each response dictate the choice of regulatory parameters such as the affinity of the transcription factor to its binding sites and the strength of the promoters it regulates. Although much is known about the dynamics of cellular responses, we still lack a framework to understand how different regulatory strategies evolved in natural systems relate to the selective pressures acting in each particular case. Here, we analyze a dynamical model of a typical antibiotic response in bacteria, where a transcriptionally repressed enzyme is induced by a sudden exposure to the drug that it processes. We identify strategies of gene regulation that optimize this response for different types of selective pressures, which we define as a set of costs associated with the drug, enzyme, and repressor concentrations during the response. We find that regulation happens in a limited region of the regulatory parameter space. While responses to more costly (toxic) drugs favor the usage of strongly self-regulated repressors, responses where expression of enzyme is more costly favor the usage of constitutively expressed repressors. Only a very narrow range of selective pressures favor weakly self-regulated repressors. We use this framework to determine which costs and benefits are most critical for the evolution of a variety of natural cellular responses that satisfy the approximations in our model and to analyze how regulation is optimized in new environments with different demands.

Publisher's Note: Suppression of Beneficial Mutations in Dynamic Microbial Populations [Phys. Rev. Lett. 118, 028102 (2017)].

This corrects the article DOI: 10.1103/PhysRevLett.118.028102.

A programmable fate decision landscape underlies single-cell aging in yeast.

Chromatin instability and mitochondrial decline are conserved processes that contribute to cellular aging. Although both processes have been explored individually in the context of their distinct signaling pathways, the mechanism that determines which process dominates during aging of individual cells is unknown. We show that interactions between the chromatin silencing and mitochondrial pathways lead to an epigenetic landscape of yeast replicative aging with multiple equilibrium states that represent different types of terminal states of aging. The structure of the landscape drives single-cell differentiation toward one of these states during aging, whereby the fate is determined quite early and is insensitive to intracellular noise. Guided by a quantitative model of the aging landscape, we genetically engineered a long-lived equilibrium state characterized by an extended life span.

Genetically engineered control of phenotypic structure in microbial colonies.

Rapid advances in cellular engineering1,2 have positioned synthetic biology to address therapeutic3,4 and industrial5 problems, but a substantial obstacle is the myriad of unanticipated cellular responses in heterogeneous real-world environments such as the gut6,7, solid tumours8,9, bioreactors10 or soil11. Complex interactions between the environment and cells often arise through non-uniform nutrient availability, which generates bidirectional coupling as cells both adjust to and modify their local environment through phenotypic differentiation12,13. Although synthetic spatial gene expression patterns14-17 have been explored under homogeneous conditions, the mutual interaction of gene circuits, growth phenotype and the environment remains a challenge. Here, we design gene circuits that sense and control phenotypic structure in microcolonies containing both growing and dormant bacteria. We implement structure modulation by coupling different downstream modules to a tunable sensor that leverages Escherichia coli's stress response and is activated on growth arrest. One is an actuator module that slows growth and thereby alters nutrient gradients. Environmental feedback in this circuit generates robust cycling between growth and dormancy in the interior of the colony, as predicted by a spatiotemporal computational model. We also use the sensor to drive an inducible gating module for selective gene expression in non-dividing cells, which allows us to radically alter population structure by eliminating the dormant phenotype with a 'stress-gated lysis circuit'. Our results establish a strategy to leverage and control microbial colony structure for synthetic biology applications in complex environments.

Advances in quantitative biology methods for studying replicative aging in Saccharomyces cerevisiae.

Aging is a complex, yet pervasive phenomenon in biology. As human cells steadily succumb to the deteriorating effects of aging, so too comes a host of age-related ailments such as neurodegenerative disorders, cardiovascular disease and cancer. Therefore, elucidation of the molecular networks that drive aging is of paramount importance to human health. Progress toward this goal has been aided by studies from simple model organisms such as Saccharomyces cerevisiae. While work in budding yeast has already revealed much about the basic biology of aging as well as a number of evolutionarily conserved pathways involved in this process, recent technological advances are poised to greatly expand our knowledge of aging in this simple eukaryote. Here, we review the latest developments in microfluidics, single-cell analysis and high-throughput technologies for studying single-cell replicative aging in S. cerevisiae. We detail the challenges each of these methods addresses as well as the unique insights into aging that each has provided. We conclude with a discussion of potential future applications of these techniques as well as the importance of single-cell dynamics and quantitative biology approaches for understanding cell aging.

Divergent Aging of Isogenic Yeast Cells Revealed through Single-Cell Phenotypic Dynamics.

Although genetic mutations that alter organisms' average lifespans have been identified in aging research, our understanding of the dynamic changes during aging remains limited. Here, we integrate single-cell imaging, microfluidics, and computational modeling to investigate phenotypic divergence and cellular heterogeneity during replicative aging of single S. cerevisiae cells. Specifically, we find that isogenic cells diverge early in life toward one of two aging paths, which are characterized by distinct age-associated phenotypes. We captured the dynamics of single cells along the paths with a stochastic discrete-state model, which accurately predicts both the measured heterogeneity and the lifespan of cells on each path within a cell population. Our analysis suggests that genetic and environmental factors influence both a cell's choice of paths and the kinetics of paths themselves. Given that these factors are highly conserved throughout eukaryotes, divergent aging might represent a general scheme in cellular aging of other organisms.

Flavin-based metabolic cycles are integral features of growth and division in single yeast cells.

The yeast metabolic cycle (YMC) is a fascinating example of biological organization, in which cells constrain the function of specific genetic, protein and metabolic networks to precise temporal windows as they grow and divide. However, understanding the intracellular origins of the YMC remains a challenging goal, as measuring the oxygen oscillations traditionally associated with it requires the use of synchronized cultures growing in nutrient-limited chemostat environments. To address these limitations, we used custom-built microfluidic devices and time-lapse fluorescence microscopy to search for metabolic cycling in the form of endogenous flavin fluorescence in unsynchronized single yeast cells. We uncovered robust and pervasive metabolic cycles that were synchronized with the cell division cycle (CDC) and oscillated across four different nutrient conditions. We then studied the response of these metabolic cycles to chemical and genetic perturbations, showing that their phase synchronization with the CDC can be altered through treatment with rapamycin, and that metabolic cycles continue even in respiratory deficient strains. These results provide a foundation for future studies of the physiological importance of metabolic cycles in processes such as CDC control, metabolic regulation and cell aging.

Coexistence and Pattern Formation in Bacterial Mixtures with Contact-Dependent Killing.

Multistrain microbial communities often exhibit complex spatial organization that emerges because of the interplay of various cooperative and competitive interaction mechanisms. One strong competitive mechanism is contact-dependent neighbor killing enabled by the type VI secretion system. It has been previously shown that contact-dependent killing can result in bistability of bacterial mixtures so that only one strain survives and displaces the other. However, it remains unclear whether stable coexistence is possible in such mixtures. Using a population dynamics model for two interacting bacterial strains, we found that coexistence can be made possible by the interplay of contact-dependent killing and long-range growth inhibition, leading to the formation of various cellular patterns. These patterns emerge in a much broader parameter range than that required for the linear Turing-like instability, suggesting this may be a robust mechanism for pattern formation.

Rational engineering of synthetic microbial systems: from single cells to consortia.

One promise of synthetic biology is to provide solutions for biomedical and industrial problems by rational design of added functionality in living systems. Microbes are at the forefront of this biological engineering endeavor due to their general ease of handling and their relevance in many potential applications from fermentation to therapeutics. In recent years, the field has witnessed an explosion of novel regulatory tools, from synthetic orthogonal transcription factors to posttranslational mechanisms for increased control over the behavior of synthetic circuits. Tool development has been paralleled by the discovery of principles that enable increased modularity and the management of host-circuit interactions. Engineered cell-to-cell communication bridges the scales from intracellular to population-level coordination. These developments facilitate the translation of more than a decade of circuit design into applications.

Multigenerational silencing dynamics control cell aging.

Cellular aging plays an important role in many diseases, such as cancers, metabolic syndromes, and neurodegenerative disorders. There has been steady progress in identifying aging-related factors such as reactive oxygen species and genomic instability, yet an emerging challenge is to reconcile the contributions of these factors with the fact that genetically identical cells can age at significantly different rates. Such complexity requires single-cell analyses designed to unravel the interplay of aging dynamics and cell-to-cell variability. Here we use microfluidic technologies to track the replicative aging of single yeast cells and reveal that the temporal patterns of heterochromatin silencing loss regulate cellular life span. We found that cells show sporadic waves of silencing loss in the heterochromatic ribosomal DNA during the early phases of aging, followed by sustained loss of silencing preceding cell death. Isogenic cells have different lengths of the early intermittent silencing phase that largely determine their final life spans. Combining computational modeling and experimental approaches, we found that the intermittent silencing dynamics is important for longevity and is dependent on the conserved Sir2 deacetylase, whereas either sustained silencing or sustained loss of silencing shortens life span. These findings reveal that the temporal patterns of a key molecular process can directly influence cellular aging, and thus could provide guidance for the design of temporally controlled strategies to extend life span.

Gene Conversion Facilitates Adaptive Evolution on Rugged Fitness Landscapes.

Gene conversion is a ubiquitous phenomenon that leads to the exchange of genetic information between homologous DNA regions and maintains coevolving multi-gene families in most prokaryotic and eukaryotic organisms. In this paper, we study its implications for the evolution of a single functional gene with a silenced duplicate, using two different models of evolution on rugged fitness landscapes. Our analytical and numerical results show that, by helping to circumvent valleys of low fitness, gene conversion with a passive duplicate gene can cause a significant speedup of adaptation, which depends nontrivially on the frequency of gene conversion and the structure of the landscape. We find that stochastic effects due to finite population sizes further increase the likelihood of exploiting this evolutionary pathway. A universal feature appearing in both deterministic and stochastic analysis of our models is the existence of an optimal gene conversion rate, which maximizes the speed of adaptation. Our results reveal the potential for duplicate genes to act as a "scratch paper" that frees evolution from being limited to strictly beneficial mutations in strongly selective environments.

A stabilized microbial ecosystem of self-limiting bacteria using synthetic quorum-regulated lysis.

Microbial ecologists are increasingly turning to small, synthesized ecosystems1-5 as a reductionist tool to probe the complexity of native microbiomes6,7. Concurrently, synthetic biologists have gone from single-cell gene circuits8-11 to controlling whole populations using intercellular signalling12-16. The intersection of these fields is giving rise to new approaches in waste recycling17, industrial fermentation18, bioremediation19 and human health16,20. These applications share a common challenge7 well-known in classical ecology21,22-stability of an ecosystem cannot arise without mechanisms that prohibit the faster-growing species from eliminating the slower. Here, we combine orthogonal quorum-sensing systems and a population control circuit with diverse self-limiting growth dynamics to engineer two 'ortholysis' circuits capable of maintaining a stable co-culture of metabolically competitive Salmonella typhimurium strains in microfluidic devices. Although no successful co-cultures are observed in a two-strain ecology without synthetic population control, the 'ortholysis' design dramatically increases the co-culture rate from 0% to approximately 80%. Agent-based and deterministic modelling reveal that our system can be adjusted to yield different dynamics, including phase-shifted, antiphase or synchronized oscillations, as well as stable steady-state population densities. The 'ortholysis' approach establishes a paradigm for constructing synthetic ecologies by developing stable communities of competitive microorganisms without the need for engineered co-dependency.

Coupled feedback loops control the stimulus-dependent dynamics of the yeast transcription factor Msn2.

Information about environmental stimuli often can be encoded by the dynamics of signaling molecules or transcription factors. In the yeast Saccharomyces cerevisiae, different types of stresses induce distinct nuclear translocation dynamics of the general stress-responsive transcription factor Msn2, but the underlying mechanisms remain unclear. Using deterministic and stochastic modeling, we reproduced in silico the different dynamic responses of Msn2 to glucose limitation and osmotic stress observed in vivo and found that a positive feedback loop on protein kinase A mediated by the AMP-activated protein kinase Snf1 is coupled with a negative feedback loop to generate the characteristic pulsatile dynamics of Msn2. The model predicted that the stimulus-specific positive feedback loop could be responsible for the difference between Msn2 dynamics induced by glucose limitation and osmotic stress. This prediction was further verified experimentally by time-lapse microscopic examinations of the snf1Δ strain. In this mutant lacking the Snf1-mediated positive feedback loop, Msn2 responds similarly to glucose limitation and osmotic stress, and its pulsatile translocation is largely abrogated. Our combined computational and experimental analysis reveals a regulatory mechanism by which cells can encode information about environmental cues into distinct signaling dynamics through stimulus-specific network architectures.

Synthetic Gene Circuits Learn to Classify.

An efficient computational algorithm is developed to design microRNA-based synthetic cell classifiers and to optimize their performance.