Field theory of enzyme-substrate systems with restricted long-range interactions.

Enzyme-substrate kinetics form the basis of many biomolecular processes. The interplay between substrate binding and substrate geometry can give rise to long-range interactions between enzyme binding events. Here we study a general model of enzyme-substrate kinetics with restricted long-range interactions described by an exponent -γ. We employ a coherent-state path integral and renormalization group approach to calculate the first moment and two-point correlation function of the enzyme-binding profile. We show that starting from an empty substrate the average occupancy follows a power law with an exponent 1/(1-γ) over time. The correlation function decays algebraically with two distinct spatial regimes characterized by exponents -γ on short distances and -(2/3)(2-γ) on long distances. The crossover between both regimes scales inversely with the average substrate occupancy. Our work allows associating experimental measurements of bound enzyme locations with their binding kinetics and the spatial conformation of the substrate.

scEpiAge: an age predictor highlighting single-cell ageing heterogeneity in mouse blood.

Ageing is the accumulation of changes and decline of function of organisms over time. The concept and biomarkers of biological age have been established, notably DNA methylation-based clocks. The emergence of single-cell DNA methylation profiling methods opens the possibility of studying the biological age of individual cells. Here, we generate a large single-cell DNA methylation and transcriptome dataset from mouse peripheral blood samples, spanning a broad range of ages. The number of genes expressed increases with age, but gene-specific changes are small. We next develop scEpiAge, a single-cell DNA methylation age predictor, which can accurately predict age in (very sparse) publicly available datasets, and also in single cells. DNA methylation age distribution is wider than technically expected, indicating epigenetic age heterogeneity and functional differences. Our work provides a foundation for single-cell and sparse data epigenetic age predictors, validates their functionality and highlights epigenetic heterogeneity during ageing.

Clonal dynamics of surface-driven growing tissues.

The self-organization of cells into complex tissues relies on a tight coordination of cell behavior. Identifying the cellular processes driving tissue growth is key to understanding the emergence of tissue forms and devising targeted therapies for aberrant growth, such as in cancer. Inferring the mode of tissue growth, whether it is driven by cells on the surface or by cells in the bulk, is possible in cell culture experiments but difficult in most tissues in living organisms (in vivo). Genetic tracing experiments, where a subset of cells is labeled with inheritable markers, have become important experimental tools to study cell fate in vivo. Here we show that the mode of tissue growth is reflected in the size distribution of the progeny of marked cells. To this end, we derive the clone size distributions using analytical calculations in the limit of negligible cell migration and cell death, and we test our predictions with an agent-based stochastic sampling technique. We show that for surface-driven growth the clone size distribution takes a characteristic power-law form with an exponent determined by fluctuations of the tissue surface. Our results propose a possible way of determining the mode of tissue growth from genetic tracing experiments.

Compartmentalization and synergy of osteoblasts drive bone formation in the regenerating fin.

Zebrafish regenerate their fins which involves a component of cell plasticity. It is currently unclear how regenerate cells divide labor to allow for appropriate growth and patterning. Here, we studied lineage relationships of fluorescence-activated cell sorting-enriched epidermal, bone-forming (osteoblast), and (non-osteoblast) blastemal fin regenerate cells by single-cell RNA sequencing, lineage tracing, targeted osteoblast ablation, and electron microscopy. Most osteoblasts in the outgrowing regenerate derive from osterix+ osteoblasts, while mmp9+ cells reside at segment joints. Distal blastema cells contribute to distal osteoblast progenitors, suggesting compartmentalization of the regenerating appendage. Ablation of osterix+ osteoblasts impairs segment joint and bone matrix formation and decreases regenerate length which is partially compensated for by distal regenerate cells. Our study characterizes expression patterns and lineage relationships of rare fin regenerate cell populations, indicates inherent detection and compensation of impaired regeneration, suggests variable dependence on growth factor signaling, and demonstrates zonation of the elongating fin regenerate.

FOXO1-mediated lipid metabolism maintains mammalian embryos in dormancy.

Mammalian developmental timing is adjustable in vivo by preserving pre-implantation embryos in a dormant state called diapause. Inhibition of the growth regulator mTOR (mTORi) pauses mouse development in vitro, yet how embryonic dormancy is maintained is not known. Here we show that mouse embryos in diapause are sustained by using lipids as primary energy source. In vitro, supplementation of embryos with the metabolite L-carnitine balances lipid consumption, puts the embryos in deeper dormancy and boosts embryo longevity. We identify FOXO1 as an essential regulator of the energy balance in dormant embryos and propose, through meta-analyses of dormant cell signatures, that it may be a common regulator of dormancy across adult tissues. Our results lift a constraint on in vitro embryo survival and suggest that lipid metabolism may be a critical metabolic transition relevant for longevity and stem cell function across tissues.

Self-organization of plasticity and specialization in a primitively social insect.

Biological systems have the capacity to not only build and robustly maintain complex structures but also to rapidly break up and rebuild such structures. Here, using primitive societies of Polistes wasps, we show that both robust specialization and rapid plasticity are emergent properties of multi-scale dynamics. We combine theory with experiments that, after perturbing the social structure by removing the queen, correlate time-resolved multi-omics with video recordings. We show that the queen-worker dimorphism relies on the balance between the development of a molecular queen phenotype in all insects and colony-scale inhibition of this phenotype via asymmetric interactions. This allows Polistes to be stable against intrinsic perturbations of molecular states while reacting plastically to extrinsic cues affecting the whole society. Long-term stability of the social structure is reinforced by dynamic DNA methylation. Our study provides a general principle of how both specialization and plasticity can be achieved in biological systems. A record of this paper's transparent peer review process is included in the supplemental information.

Making sense of fragmentation and merging in lineage tracing experiments.

Lineage tracing experiments give dynamic information on the functional behaviour of dividing cells. These experiments therefore have become an important tool for studying stem and progenitor cell fate behavior in vivo. When cell proliferation is high or the frequency of induced clones cannot be precisely controlled, the merging and fragmentation of clones renders the retrospective interpretation of clonal fate data highly ambiguous, potentially leading to unguarded interpretations about lineage relationships and fate behaviour. Here, we discuss and generalize statistical strategies to detect, resolve and make use of clonal fragmentation and merging. We first explain how to detect the rates of clonal fragmentation and merging using simple statistical estimates. We then discuss ways to restore the clonal provenance of labelled cells algorithmically and statistically and elaborate on how the process of clonal fragmentation can indirectly inform about cell fate. We generalize and extend results from the context of their original publication.

Single cell biology-a Keystone Symposia report.

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Four-pronged negative feedback of DSB machinery in meiotic DNA-break control in mice.

In most taxa, halving of chromosome numbers during meiosis requires that homologous chromosomes (homologues) pair and form crossovers. Crossovers emerge from the recombination-mediated repair of programmed DNA double-strand breaks (DSBs). DSBs are generated by SPO11, whose activity requires auxiliary protein complexes, called pre-DSB recombinosomes. To elucidate the spatiotemporal control of the DSB machinery, we focused on an essential SPO11 auxiliary protein, IHO1, which serves as the main anchor for pre-DSB recombinosomes on chromosome cores, called axes. We discovered that DSBs restrict the DSB machinery by at least four distinct pathways in mice. Firstly, by activating the DNA damage response (DDR) kinase ATM, DSBs restrict pre-DSB recombinosome numbers without affecting IHO1. Secondly, in their vicinity, DSBs trigger IHO1 depletion mainly by another DDR kinase, ATR. Thirdly, DSBs enable homologue synapsis, which promotes the depletion of IHO1 and pre-DSB recombinosomes from synapsed axes. Finally, DSBs and three DDR kinases, ATM, ATR and PRKDC, enable stage-specific depletion of IHO1 from all axes. We hypothesize that these four negative feedback pathways protect genome integrity by ensuring that DSBs form without excess, are well-distributed, and are restricted to genomic locations and prophase stages where DSBs are functional for promoting homologue pairing and crossover formation.

Active turnover of DNA methylation during cell fate decisions.

DNA methylation is a key layer of epigenetic regulation. The deposition of methylation marks relies on the catalytic activity of DNA methyltransferases (DNMTs), and their active removal relies on the activity of ten-eleven translocation (TET) enzymes. Paradoxically, in important biological contexts these antagonistic factors are co-expressed and target overlapping genomic regions. The ensuing cyclic biochemistry of cytosine modifications gives rise to a continuous, out-of-thermal equilibrium transition through different methylation states. But what is the purpose of this intriguing turnover of DNA methylation? Recent evidence demonstrates that methylation turnover is enriched at gene distal regulatory elements, including enhancers, and can give rise to large-scale oscillatory dynamics. We discuss this phenomenon and propose that DNA methylation turnover might facilitate key lineage decisions.

Lineage hierarchies and stochasticity ensure the long-term maintenance of adult neural stem cells.

The cellular basis and extent of neural stem cell (NSC) self-renewal in adult vertebrates, and their heterogeneity, remain controversial. To explore the functional behavior and dynamics of individual NSCs, we combined genetic lineage tracing, quantitative clonal analysis, intravital imaging, and global population assessments in the adult zebrafish telencephalon. Our results are compatible with a model where adult neurogenesis is organized in a hierarchy in which a subpopulation of deeply quiescent reservoir NSCs with long-term self-renewal potential generate, through asymmetric divisions, a pool of operational NSCs activating more frequently and taking stochastic fates biased toward neuronal differentiation. Our data further suggest the existence of an additional, upstream, progenitor population that supports the continuous generation of new reservoir NSCs, thus contributing to their overall expansion. Hence, we propose that the dynamics of vertebrate neurogenesis relies on a hierarchical organization where growth, self-renewal, and neurogenic functions are segregated between different NSC types.

Single cell sequencing of radial glia progeny reveals the diversity of newborn neurons in the adult zebrafish brain.

Zebrafish display widespread and pronounced adult neurogenesis, which is fundamental for their regeneration capability after central nervous system injury. However, the cellular identity and the biological properties of adult newborn neurons are elusive for most brain areas. Here, we have used short-term lineage tracing of radial glia progeny to prospectively isolate newborn neurons from the her4.1+ radial glia lineage in the homeostatic adult forebrain. Transcriptome analysis of radial glia, newborn neurons and mature neurons using single cell sequencing identified distinct transcriptional profiles, including novel markers for each population. Specifically, we detected two separate newborn neuron types, which showed diversity of cell fate commitment and location. Further analyses showed that these cell types are homologous to neurogenic cells in the mammalian brain, identified neurogenic commitment in proliferating radial glia and indicated that glutamatergic projection neurons are generated in the adult zebrafish telencephalon. Thus, we prospectively isolated adult newborn neurons from the adult zebrafish forebrain, identified markers for newborn and mature neurons in the adult brain, and revealed intrinsic heterogeneity among adult newborn neurons and their homology with mammalian adult neurogenic cell types.

Lgr5+ stem and progenitor cells reside at the apex of a heterogeneous embryonic hepatoblast pool.

During mouse embryogenesis, progenitors within the liver known as hepatoblasts give rise to adult hepatocytes and cholangiocytes. Hepatoblasts, which are specified at E8.5-E9.0, have been regarded as a homogeneous progenitor population that initiate differentiation from E13.5. Recently, scRNA-seq analysis has identified sub-populations of transcriptionally distinct hepatoblasts at E11.5. Here, we show that hepatoblasts are not only transcriptionally but also functionally heterogeneous, and that a subpopulation of E9.5-E10.0 hepatoblasts exhibit a previously unidentified early commitment to cholangiocyte fate. Importantly, we also identify a subpopulation constituting 2% of E9.5-E10.0 hepatoblasts that express the adult stem cell marker Lgr5, and generate both hepatocyte and cholangiocyte progeny that persist for the lifespan of the mouse. Combining lineage tracing and scRNA-seq, we show that Lgr5 marks E9.5-E10.0 bipotent liver progenitors residing at the apex of a hepatoblast hierarchy. Furthermore, isolated Lgr5+ hepatoblasts can be clonally expanded in vitro into embryonic liver organoids, which can commit to either hepatocyte or cholangiocyte fates. Our study demonstrates functional heterogeneity within E9.5 hepatoblasts and identifies Lgr5 as a marker for a subpopulation of bipotent liver progenitors.

Deterministic fate assignment of Müller glia cells in the zebrafish retina suggests a clonal backbone during development.

The optic cup houses multipotent retinal progenitor cells that proliferate and differentiate to form the mature retina, containing five main types of neurons and a single glial cell type, the Müller cell. Progenitors of the zebrafish optic cup generate clones that vary regarding the number and types of neurons, a process we previously showed could be described by stochastic models. Here, we present data indicating that each retinal progenitor cell, in the 24 hrs post-fertilization optic cup, is predestined to form a single Müller cell. This striking fate assignment of Müller cells reveals a dual nature of retinal lineages where stochastic mechanisms produce variable numbers of neurons while there is a strong deterministic component governing the formation of glia cells. A possible mechanism for this stereotypic fate assignment could be the maintenance of a clonal backbone during retina development, which would be similar to invertebrate and rodent cortical neurogenesis.

Neurogenin3 phosphorylation controls reprogramming efficiency of pancreatic ductal organoids into endocrine cells.

ß-cell replacement has been proposed as an effective treatment for some forms of diabetes, and in vitro methods for ß-cell generation are being extensively explored. A potential source of ß-cells comes from fate conversion of exocrine pancreatic cells into the endocrine lineage, by overexpression of three regulators of pancreatic endocrine formation and ß-cell identity, Ngn3, Pdx1 and MafA. Pancreatic ductal organoid cultures have recently been developed that can be expanded indefinitely, while maintaining the potential to differentiate into the endocrine lineage. Here, using mouse pancreatic ductal organoids, we see that co-expression of Ngn3, Pdx1 and MafA are required and sufficient to generate cells that express insulin and resemble ß-cells transcriptome-wide. Efficiency of ß-like cell generation can be significantly enhanced by preventing phosphorylation of Ngn3 protein and further augmented by conditions promoting differentiation. Taken together, our new findings underline the potential of ductal organoid cultures as a source material for generation of ß-like cells and demonstrate that post-translational regulation of reprogramming factors can be exploited to enhance ß-cell generation.

Genome-Scale Oscillations in DNA Methylation during Exit from Pluripotency.

Pluripotency is accompanied by the erasure of parental epigenetic memory, with naïve pluripotent cells exhibiting global DNA hypomethylation both in vitro and in vivo. Exit from pluripotency and priming for differentiation into somatic lineages is associated with genome-wide de novo DNA methylation. We show that during this phase, co-expression of enzymes required for DNA methylation turnover, DNMT3s and TETs, promotes cell-to-cell variability in this epigenetic mark. Using a combination of single-cell sequencing and quantitative biophysical modeling, we show that this variability is associated with coherent, genome-scale oscillations in DNA methylation with an amplitude dependent on CpG density. Analysis of parallel single-cell transcriptional and epigenetic profiling provides evidence for oscillatory dynamics both in vitro and in vivo. These observations provide insights into the emergence of epigenetic heterogeneity during early embryo development, indicating that dynamic changes in DNA methylation might influence early cell fate decisions.

Defining Lineage Potential and Fate Behavior of Precursors during Pancreas Development.

Pancreas development involves a coordinated process in which an early phase of cell segregation is followed by a longer phase of lineage restriction, expansion, and tissue remodeling. By combining clonal tracing and whole-mount reconstruction with proliferation kinetics and single-cell transcriptional profiling, we define the functional basis of pancreas morphogenesis. We show that the large-scale organization of mouse pancreas can be traced to the activity of self-renewing precursors positioned at the termini of growing ducts, which act collectively to drive serial rounds of stochastic ductal bifurcation balanced by termination. During this phase of branching morphogenesis, multipotent precursors become progressively fate-restricted, giving rise to self-renewing acinar-committed precursors that are conveyed with growing ducts, as well as ductal progenitors that expand the trailing ducts and give rise to delaminating endocrine cells. These findings define quantitatively how the functional behavior and lineage progression of precursor pools determine the large-scale patterning of pancreatic sub-compartments.

Clonal analysis of Notch1-expressing cells reveals the existence of unipotent stem cells that retain long-term plasticity in the embryonic mammary gland.

Recent lineage tracing studies have revealed that mammary gland homeostasis relies on unipotent stem cells. However, whether and when lineage restriction occurs during embryonic mammary development, and which signals orchestrate cell fate specification, remain unknown. Using a combination of in vivo clonal analysis with whole mount immunofluorescence and mathematical modelling of clonal dynamics, we found that embryonic multipotent mammary cells become lineage-restricted surprisingly early in development, with evidence for unipotency as early as E12.5 and no statistically discernable bipotency after E15.5. To gain insights into the mechanisms governing the switch from multipotency to unipotency, we used gain-of-function Notch1 mice and demonstrated that Notch activation cell autonomously dictates luminal cell fate specification to both embryonic and basally committed mammary cells. These functional studies have important implications for understanding the signals underlying cell plasticity and serve to clarify how reactivation of embryonic programs in adult cells can lead to cancer.

Universality of clone dynamics during tissue development.

The emergence of complex organs is driven by the coordinated proliferation, migration and differentiation of precursor cells. The fate behaviour of these cells is reflected in the time evolution their progeny, termed clones, which serve as a key experimental observable. In adult tissues, where cell dynamics is constrained by the condition of homeostasis, clonal tracing studies based on transgenic animal models have advanced our understanding of cell fate behaviour and its dysregulation in disease (1, 2). But what can be learned from clonal dynamics in development, where the spatial cohesiveness of clones is impaired by tissue deformations during tissue growth? Drawing on the results of clonal tracing studies, we show that, despite the complexity of organ development, clonal dynamics may converge to a critical state characterized by universal scaling behaviour of clone sizes. By mapping clonal dynamics onto a generalization of the classical theory of aerosols, we elucidate the origin and range of scaling behaviours and show how the identification of universal scaling dependences may allow lineage-specific information to be distilled from experiments. Our study shows the emergence of core concepts of statistical physics in an unexpected context, identifying cellular systems as a laboratory to study non-equilibrium statistical physics.

Muscle Stem Cells Exhibit Distinct Clonal Dynamics in Response to Tissue Repair and Homeostatic Aging.

The clonal complexity of adult stem cell pools is progressively lost during homeostatic turnover in several tissues, suggesting a decrease in the number of stem cells with distinct clonal origins. The functional impact of reduced complexity on stem cell pools, and how different tissue microenvironments may contribute to such a reduction, are poorly understood. Here, we performed clonal multicolor lineage tracing of skeletal muscle stem cells (MuSCs) to address these questions. We found that MuSC clonal complexity is maintained during aging despite heterogenous reductions in proliferative capacity, allowing aged muscle to mount a clonally diverse, albeit diminished, response to injury. In contrast, repeated bouts of tissue repair cause a progressive reduction in MuSC clonal complexity indicative of neutral drift. Consistently, biostatistical modeling suggests that MuSCs undergo symmetric expansions with stochastic fate acquisition during tissue repair. These findings establish distinct principles that underlie stem cell dynamics during homeostatic aging and muscle regeneration.