Search

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 FOXA2

SummaryCooperative DNA binding of transcription factors (TFs) integrates external stimuli and context across tissues and time. Naïve mouse embryonic stem cells are derived from early development and can sustain the 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 not present, ESRRB supports a bistable culture of cells with an embryo-like primitive endoderm identity ancillary to pluripotency. The stoichiometry between NANOG and ESRRB quantitatively influences differentiation, and in silico modeling of bipartite TF activity suggests ESRRB 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 highly stable stem cell models.

Source code release

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., Science Advances, 8 (35), ISSN:2375-2548

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 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 mechanism is well understood, the effect of the combined mechanisms is unclear. To evaluate the contribution of the deterministic component we mapped the growing population to the Ising model. Model results, confirmed by simulations and experimental data, show that cell-cell variations increase near-linearly with stress. As a consequence, we predict that the entropic gain — the gain in population doubling time compared to an “average” cell — is primarily stochastic at low stress but crosses 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.

Significance How repressive heterochromatic states propagate along chromosomes and are subsequently maintained for many cell divisions remain important unanswered questions in biology. Combining mathematical modeling and single-cell measurements, we find that heterochromatin does not propagate in a purely linear manner as often assumed. Rather, sudden transitions can affect large domains globally at once. We suggest this is due to long-range interactions that bring distant nucleosomes close to each other to facilitate their modification. This supports the notion that the compartmentalization of chromatin components within the nucleus and the folding of chromatin into loops dynamically control heterochromatin propagation in three dimensions by bringing nucleosomes and their modifiers into close proximity., 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.

The differentiation of ICM cells into epiblast (EPI) and primitive endoderm (PE) is central in embryonic development. It is known that FGF4 signaling is important in this process, but it remains unclear how cells can be correctly partitioned. Here we model the NANOG-GATA6-FGF4 network, and test all 64 logical regulatory combinations for their ability to partition a group of cells. We found that nearly all the logic combinations allowed for correct partitioning, including a minimal network where self-activation of NANOG and GATA6 was inactivated. However such self-activation increased the robustness of the system. Furthermore, the model also captured the reported changes in cell proportions in response to FGF perturbations. This constrains the possible regulatory logic and predicts the presence of an “OR” gate in cell-cell communication. We repeatedly found that FGF4 coordinated the decision in two phases: A convergence and a bifurcation phase. First FGF4 negative feedback drives the cells to a balanced “battle” state where most cells have intermediate levels of both regulators, thus being double positive. Subsequent bifurcation happens at constant FGF4 level. Together our results suggest that the frequently observed state of multipotency during differentiation may be an emergent phenomenon resulting from inter-cellular negative feedbacks.

Summary 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., Graphical Abstract, Highlights • Cell wedging and intercalation are modeled using a polarized point-particle approach • Cell intercalation is sufficient for tube budding • Tube budding is more robust when intercalation is complemented by wedging • Wedging and differential proliferation are sufficient for mammalian neurulation, Developmental Biology; Experimental Models in Systems Biology

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 agents 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 inS. cerevisiae, we show that ER protein reflux during ER stress requires specific chaperones and co-chaperones 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.SIGNIFICANCEApproximately one third of eukaryotic proteins are synthesized on ribosomes attached to the endoplasmic reticulum (ER) membrane. Many of these polypeptides co- or post-translationally translocate into the ER, wherein they fold and mature. An ER quality-control system proofreads these proteins by facilitating their folding and modification, while eliminating misfolded proteins through ER-associated degradation (ERAD). Yet, the fate of many secretory proteins during ER stress is not completely understood. Here, we uncovered an ER-stress induced “protein reflux” system that delivers intact, folded ER luminal proteins back to the cytosol without degrading them. We found that ER protein reflux works in parallel to ERAD and requires distinct ER-resident and cytosolic chaperones and co-chaperones.

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 lead to cell-to-cell differences in free protein abundance between individuals. 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.Stand-first textUsing 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.Bullet pointsThe presented approach identifies and quantifies different fitness effects associated with protein deposit formation due to phase separationThe environmental condition and the cellular demand for the protein function emerge as key determinants of fitness upon protein deposit formationVariability in protein deposit formation can lead to cell-to-cell differences in free protein abundance between individualsProtein phase separation can generate a continuous range of phenotypes in a cell population

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., eLife digest Cells have the power to organise themselves to form complex and stable structures, whether it is to create a fully shaped baby from a single egg, or to allow adult salamanders to grow a new limb after losing a leg. This ability has been scrutinised at many different levels. For example, researchers have looked at the chemical messages exchanged by cells, or they have recorded the different shapes an embryo goes through during development. However, it is still difficult to reconcile the information from these approaches into a description that makes sense at multiple scales. When an embryo develops, sheets of cells fold and unfold to create complex 3D shapes, like the tubes that make our lungs. Moulding sheets into tubes relies on interactions between cells that are not the same in all directions. In fact, two types of asymmetry (or polarity) guide these interactions. Apical-basal polarity runs across a sheet of cells, which means that the top surface of the sheet differs from the bottom. Planar cell polarity runs along the sheet and distinguishes one end from the other. For instance, apical-basal polarity marks the inner and outer surfaces of our skin, while planar cell polarity controls the direction in which our hair grows. Nissen et al. set out to investigate how these polarities help cells in an embryo organise themselves to form complicated folds and tubes. To do this, simple mathematical representations of both apical-basal and planar cell polarities were designed. The representations were then combined to create computer simulations of groups of cells as these divide and interact with each other. Simulations of ‘cells’ with only apical-basal polarity were able to generate different shapes in the ‘tissues’ produced, including many found in living organisms. External conditions, such as how cells were arranged to start with, determined the resulting shape. With both apical-basal and planar cell polarities, the simulations reproduced an important change that occurs during early development. They also replicated how the tubes that transport nutrients and oxygen form. These results show that simple properties of individual cells, such as polarities, can produce different shapes in developing tissues and organs, without the need for a complicated overarching program. Abnormal changes in cell polarity are also associated with diseases such as cancer. The mathematical model developed by Nissen et al. could therefore be a useful tool to study these events.

Allele number, or zygosity, is a clear determinant of gene expression in diploid cells. But 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 feedbacks. 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 between donors infected with human cytomegalovirus (HCMV) with one or two copies of the allele. Crucially, the distribution of gene expression in many of the donors is bimodal, indicative of the presence of a positive feedback somewhere in the regulatory environment of the gene. Three separate gene-circuit models differing in the location of the positive feedback with respect to the gene can all reproduce well the homozygous data. However, when the resulting fitted models are applied to the hemizygous donors, only one model (the one with the positive feedback located at the level of gene transcription) reproduces the experimentally observed gene-expression profile. In that way, our work shows that zygosity can help us relate structure and function of gene regulatory networks.Author SummaryNearly all mammalian cells, including human cells, have two copies of each chromosome, and thus possess two potentially different copies of each gene (which might be in some cases non-functional or even absent). Naïively one might expect that two identical copies of the gene would lead to the protein being expressed at twice the rate, but many factors can alter this simple calculation. One of these factors is the existence of feedback mechanisms affecting in one way or another the regulatory circuit in which our gene of interest is embedded. Here we study the relationship between the number of gene copies and the expression of a receptor protein that plays a crucial role in the recognition of pathogens by natural killer cells, which are important elements of the innate immune system. Experimental data of virus-infected donors reveals a bimodal expression profile of this receptor, typical of a positive feedback, and a clear difference between donors with one or two copies of the gene. Mathematical modeling allows us to find the likely location of the feedback loop within the gene’s regulatory circuit, by requiring the correct model to reproduce the expression profiles of both types of donors.

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.