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.

Defining network topologies that can achieve biochemical adaptation.

Many signaling systems show adaptation-the ability to reset themselves after responding to a stimulus. We computationally searched all possible three-node enzyme network topologies to identify those that could perform adaptation. Only two major core topologies emerge as robust solutions: a negative feedback loop with a buffering node and an incoherent feedforward loop with a proportioner node. Minimal circuits containing these topologies are, within proper regions of parameter space, sufficient to achieve adaptation. More complex circuits that robustly perform adaptation all contain at least one of these topologies at their core. This analysis yields a design table highlighting a finite set of adaptive circuits. Despite the diversity of possible biochemical networks, it may be common to find that only a finite set of core topologies can execute a particular function. These design rules provide a framework for functionally classifying complex natural networks and a manual for engineering networks. For a video summary of this article, see the PaperFlick file with the Supplemental Data available online.

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.

Real-time redox measurements during endoplasmic reticulum stress reveal interlinked protein folding functions.

Disruption of protein folding in the endoplasmic reticulum (ER) causes unfolded proteins to accumulate, triggering the unfolded protein response (UPR). UPR outputs in turn decrease ER unfolded proteins to close a negative feedback loop. However, because it is infeasible to directly measure the concentration of unfolded proteins in vivo, cells are generically described as experiencing "ER stress" whenever the UPR is active. Because ER redox potential is optimized for oxidative protein folding, we reasoned that measureable redox changes should accompany unfolded protein accumulation. To test this concept, we employed fluorescent protein reporters to dynamically measure ER redox status and UPR activity in single cells. Using these tools, we show that diverse stressors, both experimental and physiological, compromise ER protein oxidation when UPR-imposed homeostatic control is lost. Using genetic analysis we uncovered redox heterogeneities in isogenic cell populations, and revealed functional interlinks between ER protein folding, modification, and quality control systems.

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.

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.

Two portable recombination enhancers direct donor choice in fission yeast heterochromatin.

Mating-type switching in fission yeast results from gene conversions of the active mat1 locus by heterochromatic donors. mat1 is preferentially converted by mat2-P in M cells and by mat3-M in P cells. Here, we report that donor choice is governed by two portable recombination enhancers capable of promoting use of their adjacent cassette even when they are transposed to an ectopic location within the mat2-mat3 heterochromatic domain. Cells whose silent cassettes are swapped to mat2-M mat3-P switch mating-type poorly due to a defect in directionality but cells whose recombination enhancers were transposed together with the cassette contents switched like wild type. Trans-acting mutations that impair directionality affected the wild-type and swapped cassettes in identical ways when the recombination enhancers were transposed together with their cognate cassette, showing essential regulatory steps occur through the recombination enhancers. Our observations lead to a model where heterochromatin biases competitions between the two recombination enhancers to achieve directionality.

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.

Conditional cooperativity in toxin-antitoxin regulation prevents random toxin activation and promotes fast translational recovery.

Many toxin-antitoxin (TA) loci are known to strongly repress their own transcription. This auto-inhibition is often called 'conditional cooperativity' as it relies on cooperative binding of TA complexes to operator DNA that occurs only when toxins are in a proper stoichiometric relationship with antitoxins. There has recently been an explosion of interest in TA systems due to their role in bacterial persistence, however the role of conditional cooperativity is still unclear. We reveal the biological function of conditional cooperativity by constructing a mathematical model of the well studied TA system, relBE of Escherichia coli. We show that the model with the in vivo and in vitro established parameters reproduces experimentally observed response to nutritional stress. We further demonstrate that conditional cooperativity stabilizes the level of antitoxin in rapidly growing cells such that random induction of relBE is minimized. At the same time it enables quick removal of free toxin when the starvation is terminated.

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.

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.

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.

Targeted bacterial immunity buffers phage diversity.

Bacteria have evolved diverse defense mechanisms that allow them to fight viral attacks. One such mechanism, the clustered, regularly interspaced, short palindromic repeat (CRISPR) system, is an adaptive immune system consisting of genetic loci that can take up genetic material from invasive elements (viruses and plasmids) and later use them to reject the returning invaders. It remains an open question how, despite the ongoing evolution of attack and defense mechanisms, bacteria and viral phages manage to coexist. Using a simple mathematical model and a two-dimensional numerical simulation, we found that CRISPR adaptive immunity allows for robust phage-bacterium coexistence even when the number of virus species far exceeds the capacity of CRISPR-encoded genetic memory. Coexistence is predicted to be a consequence of the presence of many interdependent species that stress but do not overrun the bacterial defense system.

Circuit architecture explains functional similarity of bacterial heat shock responses.

Heat shock response is a stress response to temperature changes and a consecutive increase in amounts of unfolded proteins. To restore homeostasis, cells upregulate chaperones facilitating protein folding by means of transcription factors (TFs). We here investigate two heat shock systems: one characteristic to gram negative bacteria, mediated by transcriptional activator σ(32) in E. coli, and another characteristic to gram positive bacteria, mediated by transcriptional repressor HrcA in L. lactis. We construct simple mathematical models of the two systems focusing on the negative feedbacks, where free chaperones suppress σ(32) activation in the former, while they activate HrcA repression in the latter. We demonstrate that both systems, in spite of the difference at the TF regulation level, are capable of showing very similar heat shock dynamics. We find that differences in regulation impose distinct constraints on chaperone-TF binding affinities: the binding constant of free σ(32) to chaperone DnaK, known to be in 100 nM range, set the lower limit of amount of free chaperone that the system can sense the change at the heat shock, while the binding affinity of HrcA to chaperone GroE set the upper limit and have to be rather large extending into the micromolar range.

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.

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.

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.

Ecosystems with mutually exclusive interactions self-organize to a state of high diversity.

Ecological systems comprise an astonishing diversity of species that cooperate or compete with each other forming complex mutual dependencies. The minimum requirements to maintain a large species diversity on long time scales are in general unknown. Using lichen communities as an example, we propose a model for the evolution of mutually excluding organisms that compete for space. We suggest that chainlike or cyclic invasions open for creation of spatially separated subpopulations that subsequently can lead to increased diversity. In contrast to its nonspatial counterpart, our model predicts robust coexistence of a large number of species. It is demonstrated that large species diversity can be obtained on evolutionary time scales, provided that interactions between species have spatial constraints. In particular, a phase transition to a sustainable state of high diversity is identified.