Synergistic and antagonistic effects of deterministic and stochastic cell-cell variations.

By diversifying, cells in a clonal population can together overcome the limits of individuals. Diversity in single-cell growth rates allows the population to survive environmental stresses, such as antibiotics, and grow faster than the undiversified population. These functional cell-cell variations can arise stochastically, from noise in biochemical reactions, or deterministically, by asymmetrically distributing damaged components. While each of the mechanisms is well understood, the effect of the combined mechanisms is unclear. To evaluate the contribution of the deterministic component we developed a mathematical model by mapping the growing population to the Ising model. To analyze the combined effects of stochastic and deterministic contributions we introduced the analytical results of the Ising-mapping into an Euler-Lotka framework. Model results, confirmed by simulations and experimental data, show that deterministic cell-cell variations increase near-linearly with stress. As a consequence, we predict that the gain in population doubling time from cell-cell variations is primarily stochastic at low stress but may cross over to deterministic at higher stresses. Furthermore, we find that while the deterministic component minimizes population damage, stochastic variations antagonize this effect. Together our results may help identifying stress-tolerant pathogenic cells and thus inspire novel antibiotic strategies.

Lineage tracing identifies heterogeneous hepatoblast contribution to cell lineages and postembryonic organ growth dynamics.

To meet the physiological demands of the body, organs need to establish a functional tissue architecture and adequate size as the embryo develops to adulthood. In the liver, uni- and bipotent progenitor differentiation into hepatocytes and biliary epithelial cells (BECs), and their relative proportions, comprise the functional architecture. Yet, the contribution of individual liver progenitors at the organ level to both fates, and their specific proportion, is unresolved. Combining mathematical modelling with organ-wide, multispectral FRaeppli-NLS lineage tracing in zebrafish, we demonstrate that a precise BEC-to-hepatocyte ratio is established (i) fast, (ii) solely by heterogeneous lineage decisions from uni- and bipotent progenitors, and (iii) independent of subsequent cell type-specific proliferation. Extending lineage tracing to adulthood determined that embryonic cells undergo spatially heterogeneous three-dimensional growth associated with distinct environments. Strikingly, giant clusters comprising almost half a ventral lobe suggest lobe-specific dominant-like growth behaviours. We show substantial hepatocyte polyploidy in juveniles representing another hallmark of postembryonic liver growth. Our findings uncover heterogeneous progenitor contributions to tissue architecture-defining cell type proportions and postembryonic organ growth as key mechanisms forming the adult liver.

A bipartite function of ESRRB can integrate signaling over time to balance self-renewal and differentiation.

Cooperative DNA binding of transcription factors (TFs) integrates the cellular context to support cell specification during development. Naive mouse embryonic stem cells are derived from early development and can sustain their pluripotent identity indefinitely. Here, we ask whether TFs associated with pluripotency evolved to directly support this state or if the state emerges from their combinatorial action. NANOG and ESRRB are key pluripotency factors that co-bind DNA. We find that when both factors are expressed, ESRRB supports pluripotency. However, when NANOG is absent, ESRRB supports a bistable culture of cells with an embryo-like primitive endoderm identity ancillary to pluripotency. The stoichiometry between NANOG and ESRRB allows quantitative titration of this differentiation, and in silico modeling of bipartite ESRRB activity suggests it safeguards plasticity in differentiation. Thus, the concerted activity of cooperative TFs can transform their effect to sustain intermediate cell identities and allow ex vivo expansion of immortal stem cells. A record of this paper's transparent peer review process is included in the supplemental information.

Jag1-Notch cis-interaction determines cell fate segregation in pancreatic development.

The Notch ligands Jag1 and Dll1 guide differentiation of multipotent pancreatic progenitor cells (MPCs) into unipotent pro-acinar cells (PACs) and bipotent duct/endocrine progenitors (BPs). Ligand-mediated trans-activation of Notch receptors induces oscillating expression of the transcription factor Hes1, while ligand-receptor cis-interaction indirectly represses Hes1 activation. Despite Dll1 and Jag1 both displaying cis- and trans-interactions, the two mutants have different phenotypes for reasons not fully understood. Here, we present a mathematical model that recapitulates the spatiotemporal differentiation of MPCs into PACs and BPs. The model correctly captures cell fate changes in Notch pathway knockout mice and small molecule inhibitor studies, and a requirement for oscillatory Hes1 expression to maintain the multipotent state. Crucially, the model entails cell-autonomous attenuation of Notch signaling by Jag1-mediated cis-inhibition in MPC differentiation. The model sheds light on the underlying mechanisms, suggesting that cis-interaction is crucial for exiting the multipotent state, while trans-interaction is required for adopting the bipotent fate.

Transcriptional heterogeneity and cell cycle regulation as central determinants of Primitive Endoderm priming.

During embryonic development cells acquire identity as they proliferate, implying that an intrinsic facet of cell fate choice requires coupling lineage decisions to cell division. How is the cell cycle regulated to promote or suppress heterogeneity and differentiation? We explore this question combining time lapse imaging with single-cell RNA-seq in the contexts of self-renewal, priming, and differentiation of mouse embryonic stem cells (ESCs) towards the Primitive Endoderm (PrE) lineage. Since ESCs are derived from the inner cell mass (ICM) of the mammalian blastocyst, ESCs in standard culture conditions are transcriptionally heterogeneous containing dynamically interconverting subfractions primed for either of the two ICM lineages, Epiblast and PrE. Here, we find that differential regulation of cell cycle can tip the balance between these primed populations, such that naïve ESC culture promotes Epiblast-like expansion and PrE differentiation stimulates the selective survival and proliferation of PrE-primed cells. In endoderm differentiation, this change is accompanied by a counter-intuitive increase in G1 length, also observed in vivo. While fibroblast growth factor/extracellular signal-regulated kinase (FGF/ERK) signalling is a key regulator of ESC differentiation and PrE specification, we find it is not just responsible for ESCs heterogeneity, but also the inheritance of similar cell cycles between sisters and cousins. Taken together, our results indicate a tight relationship between transcriptional heterogeneity and cell cycle regulation in lineage specification, with primed cell populations providing a pool of flexible cell types that can be expanded in a lineage-specific fashion while allowing plasticity during early determination.

Spatiotemporal NF-κB dynamics encodes the position, amplitude, and duration of local immune inputs.

Infected cells communicate through secreted signaling molecules like cytokines, which carry information about pathogens. How differences in cytokine secretion affect inflammatory signaling over space and how responding cells decode information from propagating cytokines are not understood. By computationally and experimentally studying NF-κB dynamics in cocultures of signal-sending cells (macrophages) and signal-receiving cells (fibroblasts), we find that cytokine signals are transmitted by wave-like propagation of NF-κB activity and create well-defined activation zones in responding cells. NF-κB dynamics in responding cells can simultaneously encode information about cytokine dose, duration, and distance to the cytokine source. Spatially resolved transcriptional analysis reveals that responding cells transmit local cytokine information to distance-specific proinflammatory gene expression patterns, creating "gene expression zones." Despite single-cell variability, the size and duration of the signaling zone are tightly controlled by the macrophage secretion profile. Our results highlight how macrophages tune cytokine secretion to control signal transmission distance and how inflammatory signaling interprets these signals in space and time.

Identification of the central intermediate in the extra-embryonic to embryonic endoderm transition through single-cell transcriptomics.

High-resolution maps of embryonic development suggest that acquisition of cell identity is not limited to canonical germ layers but proceeds via alternative routes. Despite evidence that visceral organs are formed via embryonic and extra-embryonic trajectories, the production of organ-specific cell types in vitro focuses on the embryonic one. Here we resolve these differentiation routes using massively parallel single-cell RNA sequencing to generate datasets from FOXA2Venus reporter mouse embryos and embryonic stem cell differentiation towards endoderm. To relate cell types in these datasets, we develop a single-parameter computational approach and identify an intermediate en route from extra-embryonic identity to embryonic endoderm, which we localize spatially in embryos at embryonic day 7.5. While there is little evidence for this cell type in embryonic stem cell differentiation, by following the extra-embryonic trajectory starting with naïve extra-embryonic endoderm stem cells we can generate embryonic gut spheroids. Exploiting developmental plasticity therefore offers alternatives to pluripotent cells and opens alternative avenues for in vitro differentiation.

Establishment of heterochromatin in domain-size-dependent bursts.

Methylation of histone H3K9 is a hallmark of epigenetic silencing in eukaryotes. Nucleosome modifications often rely on positive feedback where enzymes are recruited by modified nucleosomes. A combination of local and global feedbacks has been proposed to account for some dynamic properties of heterochromatin, but the range at which the global feedbacks operate and the exact mode of heterochromatin propagation are not known. We investigated these questions in fission yeast. Guided by mathematical modeling, we incrementally increased the size of the mating-type region and profiled heterochromatin establishment over time. We observed exponential decays in the proportion of cells with active reporters, with rates that decreased with domain size. Establishment periods varied from a few generations in wild type to >200 generations in the longest region examined, and highly correlated silencing of two reporters located outside the nucleation center was observed. On a chromatin level, this indicates that individual regions are silenced in sudden bursts. Mathematical modeling accounts for these bursts if heterochromatic nucleosomes facilitate a deacetylation or methylation reaction at long range, in a distance-independent manner. A likely effector of three-dimensional interactions is the evolutionarily conserved Swi6HP1 H3K9me reader, indicating the bursting behavior might be a general mode of heterochromatin propagation.

Model to Link Cell Shape and Polarity with Organogenesis.

How do flat sheets of cells form gut and neural tubes? Across systems, several mechanisms are at play: cells wedge, form actomyosin cables, or intercalate. As a result, the cell sheet bends, and the tube elongates. It is unclear to what extent each mechanism can drive tube formation on its own. To address this question, we computationally probe if one mechanism, either cell wedging or intercalation, may suffice for the entire sheet-to-tube transition. Using a physical model with epithelial cells represented by polarized point particles, we show that either cell intercalation or wedging alone can be sufficient and that each can both bend the sheet and extend the tube. When working in parallel, the two mechanisms increase the robustness of the tube formation. The successful simulations of the key features in Drosophila salivary gland budding, sea urchin gastrulation, and mammalian neurulation support the generality of our results.

Self-assembly, buckling and density-invariant growth of three-dimensional vascular networks.

The experimental actualization of organoids modelling organs from brains to pancreases has revealed that much of the diverse morphologies of organs are emergent properties of simple intercellular 'rules' and not the result of top-down orchestration. In contrast to other organs, the initial plexus of the vascular system is formed by aggregation of cells in the process known as vasculogenesis. Here we study this self-assembling process of blood vessels in three dimensions through a set of simple rules that align intercellular apical-basal and planar cell polarity. We demonstrate that a fully connected network of tubes emerges above a critical initial density of cells. Through planar cell polarity, our model demonstrates convergent extension, and this polarity furthermore allows for both morphology-maintaining growth and growth-induced buckling. We compare this buckling with the special vasculature of the islets of Langerhans in the pancreas and suggest that the mechanism behind the vascular density-maintaining growth of these islets could be the result of growth-induced buckling.

Chaperone-mediated reflux of secretory proteins to the cytosol during endoplasmic reticulum stress.

Diverse perturbations to endoplasmic reticulum (ER) functions compromise the proper folding and structural maturation of secretory proteins. To study secretory pathway physiology during such "ER stress," we employed an ER-targeted, redox-responsive, green fluorescent protein-eroGFP-that reports on ambient changes in oxidizing potential. Here we find that diverse ER stress regimes cause properly folded, ER-resident eroGFP (and other ER luminal proteins) to "reflux" back to the reducing environment of the cytosol as intact, folded proteins. By utilizing eroGFP in a comprehensive genetic screen in Saccharomyces cerevisiae, we show that ER protein reflux during ER stress requires specific chaperones and cochaperones residing in both the ER and the cytosol. Chaperone-mediated ER protein reflux does not require E3 ligase activity, and proceeds even more vigorously when these ER-associated degradation (ERAD) factors are crippled, suggesting that reflux may work in parallel with ERAD. In summary, chaperone-mediated ER protein reflux may be a conserved protein quality control process that evolved to maintain secretory pathway homeostasis during ER protein-folding stress.

The fitness cost and benefit of phase-separated protein deposits.

Phase separation of soluble proteins into insoluble deposits is associated with numerous diseases. However, protein deposits can also function as membrane-less compartments for many cellular processes. What are the fitness costs and benefits of forming such deposits in different conditions? Using a model protein that phase-separates into deposits, we distinguish and quantify the fitness contribution due to the loss or gain of protein function and deposit formation in yeast. The environmental condition and the cellular demand for the protein function emerge as key determinants of fitness. Protein deposit formation can influence cell-to-cell variation in free protein abundance between individuals of a cell population (i.e., gene expression noise). This results in variable manifestation of protein function and a continuous range of phenotypes in a cell population, favoring survival of some individuals in certain environments. Thus, protein deposit formation by phase separation might be a mechanism to sense protein concentration in cells and to generate phenotypic variability. The selectable phenotypic variability, previously described for prions, could be a general property of proteins that can form phase-separated assemblies and may influence cell fitness.

Theoretical tool bridging cell polarities with development of robust morphologies.

Despite continual renewal and damages, a multicellular organism is able to maintain its complex morphology. How is this stability compatible with the complexity and diversity of living forms? Looking for answers at protein level may be limiting as diverging protein sequences can result in similar morphologies. Inspired by the progressive role of apical-basal and planar cell polarity in development, we propose that stability, complexity, and diversity are emergent properties in populations of proliferating polarized cells. We support our hypothesis by a theoretical approach, developed to effectively capture both types of polar cell adhesions. When applied to specific cases of development - gastrulation and the origins of folds and tubes - our theoretical tool suggests experimentally testable predictions pointing to the strength of polar adhesion, restricted directions of cell polarities, and the rate of cell proliferation to be major determinants of morphological diversity and stability.

Stochastic priming and spatial cues orchestrate heterogeneous clonal contribution to mouse pancreas organogenesis.

Spatiotemporal balancing of cellular proliferation and differentiation is crucial for postnatal tissue homoeostasis and organogenesis. During embryonic development, pancreatic progenitors simultaneously proliferate and differentiate into the endocrine, ductal and acinar lineages. Using in vivo clonal analysis in the founder population of the pancreas here we reveal highly heterogeneous contribution of single progenitors to organ formation. While some progenitors are bona fide multipotent and contribute progeny to all major pancreatic cell lineages, we also identify numerous unipotent endocrine and ducto-endocrine bipotent clones. Single-cell transcriptional profiling at E9.5 reveals that endocrine-committed cells are molecularly distinct, whereas multipotent and bipotent progenitors do not exhibit different expression profiles. Clone size and composition support a probabilistic model of cell fate allocation and in silico simulations predict a transient wave of acinar differentiation around E11.5, while endocrine differentiation is proportionally decreased. Increased proliferative capacity of outer progenitors is further proposed to impact clonal expansion.

Four simple rules that are sufficient to generate the mammalian blastocyst.

Early mammalian development is both highly regulative and self-organizing. It involves the interplay of cell position, predetermined gene regulatory networks, and environmental interactions to generate the physical arrangement of the blastocyst with precise timing. However, this process occurs in the absence of maternal information and in the presence of transcriptional stochasticity. How does the preimplantation embryo ensure robust, reproducible development in this context? It utilizes a versatile toolbox that includes complex intracellular networks coupled to cell-cell communication, segregation by differential adhesion, and apoptosis. Here, we ask whether a minimal set of developmental rules based on this toolbox is sufficient for successful blastocyst development, and to what extent these rules can explain mutant and experimental phenotypes. We implemented experimentally reported mechanisms for polarity, cell-cell signaling, adhesion, and apoptosis as a set of developmental rules in an agent-based in silico model of physically interacting cells. We find that this model quantitatively reproduces specific mutant phenotypes and provides an explanation for the emergence of heterogeneity without requiring any initial transcriptional variation. It also suggests that a fixed time point for the cells' competence of fibroblast growth factor (FGF)/extracellular signal-regulated kinase (ERK) sets an embryonic clock that enables certain scaling phenomena, a concept that we evaluate quantitatively by manipulating embryos in vitro. Based on these observations, we conclude that the minimal set of rules enables the embryo to experiment with stochastic gene expression and could provide the robustness necessary for the evolutionary diversification of the preimplantation gene regulatory network.

Impact of Zygosity on Bimodal Phenotype Distributions.

Allele number, or zygosity, is a clear determinant of gene expression in diploid cells. However, the relationship between the number of copies of a gene and its expression can be hard to anticipate, especially when the gene in question is embedded in a regulatory circuit that contains feedback. Here, we study this question making use of the natural genetic variability of human populations, which allows us to compare the expression profiles of a receptor protein in natural killer cells among donors infected with human cytomegalovirus with one or two copies of the allele. Crucially, the distribution of gene expression in many of the donors is bimodal, which indicates the presence of a positive feedback loop somewhere in the regulatory environment of the gene. Three separate gene-circuit models differing in the location of the positive feedback loop with respect to the gene can all reproduce the homozygous data. However, when the resulting fitted models are applied to the hemizygous donors, one model (the one with the positive feedback located at the level of gene transcription) is superior in describing the experimentally observed gene-expression profile. In that way, our work shows that zygosity can help us relate the structure and function of gene regulatory networks.

Evolution of a G protein-coupled receptor response by mutations in regulatory network interactions.

All cellular functions depend on the concerted action of multiple proteins organized in complex networks. To understand how selection acts on protein networks, we used the yeast mating receptor Ste2, a pheromone-activated G protein-coupled receptor, as a model system. In Saccharomyces cerevisiae, Ste2 is a hub in a network of interactions controlling both signal transduction and signal suppression. Through laboratory evolution, we obtained 21 mutant receptors sensitive to the pheromone of a related yeast species and investigated the molecular mechanisms behind this newfound sensitivity. While some mutants show enhanced binding affinity to the foreign pheromone, others only display weakened interactions with the network's negative regulators. Importantly, the latter changes have a limited impact on overall pathway regulation, despite their considerable effect on sensitivity. Our results demonstrate that a new receptor-ligand pair can evolve through network-altering mutations independently of receptor-ligand binding, and suggest a potential role for such mutations in disease.

Asymmetric Damage Segregation Constitutes an Emergent Population-Level Stress Response.

Asymmetric damage segregation (ADS) is a mechanism for increasing population fitness through non-random, asymmetric partitioning of damaged macromolecules at cell division. ADS has been reported across multiple organisms, though the measured effects on fitness of individuals are often small. Here, we introduce a cell-lineage-based framework that quantifies the population-wide effects of ADS and then verify our results experimentally in E. coli under heat and antibiotic stress. Using an experimentally validated mathematical model, we find that the beneficial effect of ADS increases with stress. In effect, low-damage subpopulations divide faster and amplify within the population acting like a positive feedback loop whose strength scales with stress. Analysis of protein aggregates shows that the degree of asymmetric inheritance is damage dependent in single cells. Together our results indicate that, despite small effects in single cell, ADS exerts a strong beneficial effect on the population level and arises from the redistribution of damage within a population, through both single-cell and population-level feedback.

Nucleation and spreading of a heterochromatic domain in fission yeast.

Outstanding questions in the chromatin field bear on how large heterochromatin domains are formed in space and time. Positive feedback, where histone-modifying enzymes are attracted to chromosomal regions displaying the modification they catalyse, is believed to drive the formation of these domains; however, few quantitative studies are available to assess this hypothesis. Here we quantified the de novo establishment of a naturally occurring ∼20-kb heterochromatin domain in fission yeast through single-cell analyses, measuring the kinetics of heterochromatin nucleation in a region targeted by RNAi and its subsequent expansion. We found that nucleation of heterochromatin is stochastic and can take from one to ten cell generations. Further silencing of the full region takes another one to ten generations. Quantitative modelling of the observed kinetics emphasizes the importance of local feedback, where a nucleosome-bound enzyme modifies adjacent nucleosomes, combined with a feedback where recruited enzymes can act at a distance.

Dynamics of the DNA repair proteins WRN and BLM in the nucleoplasm and nucleoli.

We have investigated the mobility of two EGFP-tagged DNA repair proteins, WRN and BLM. In particular, we focused on the dynamics in two locations, the nucleoli and the nucleoplasm. We found that both WRN and BLM use a "DNA-scanning" mechanism, with rapid binding-unbinding to DNA resulting in effective diffusion. In the nucleoplasm WRN and BLM have effective diffusion coefficients of 1.62 and 1.34 µm(2)/s, respectively. Likewise, the dynamics in the nucleoli are also best described by effective diffusion, but with diffusion coefficients a factor of ten lower than in the nucleoplasm. From this large reduction in diffusion coefficient we were able to classify WRN and BLM as DNA damage scanners. In addition to WRN and BLM we also classified other DNA damage proteins and found they all fall into one of two categories. Either they are scanners, similar to WRN and BLM, with very low diffusion coefficients, suggesting a scanning mechanism, or they are almost freely diffusing, suggesting that they interact with DNA only after initiation of a DNA damage response.