Embryonic development is a crucial period in the life of a multicellular organism, during which limited sets of embryonic progenitors produce all cells in the adult body. Determining which fate these progenitors acquire in adult tissues requires the simultaneous measurement of clonal history and cell identity at single-cell resolution, which has been a major challenge. Clonal history has traditionally been investigated by microscopically tracking cells during development, monitoring the heritable expression of genetically encoded fluorescent proteins and, more recently, using next-generation sequencing technologies that exploit somatic mutations, microsatellite instability, transposon tagging, viral barcoding, CRISPR-Cas9 genome editing and Cre-loxP recombination. Single-cell transcriptomics provides a powerful platform for unbiased cell-type classification. Here we present ScarTrace, a single-cell sequencing strategy that enables the simultaneous quantification of clonal history and cell type for thousands of cells obtained from different organs of the adult zebrafish. Using ScarTrace, we show that a small set of multipotent embryonic progenitors generate all haematopoietic cells in the kidney marrow, and that many progenitors produce specific cell types in the eyes and brain. In addition, we study when embryonic progenitors commit to the left or right eye. ScarTrace reveals that epidermal and mesenchymal cells in the caudal fin arise from the same progenitors, and that osteoblast-restricted precursors can produce mesenchymal cells during regeneration. Furthermore, we identify resident immune cells in the fin with a distinct clonal origin from other blood cell types. We envision that similar approaches will have major applications in other experimental systems, in which the matching of embryonic clonal origin to adult cell type will ultimately allow reconstruction of how the adult body is built from a single cell.

The body plan of the mammalian embryo is shaped through the process of gastrulation, an early developmental event that transforms an isotropic group of cells into an ensemble of tissues ordered with reference to three orthogonal axes 1 . While model organisms have provided much insight into this process, we know very little about gastrulation in humans due to the difficulty of obtaining embryos at such early stages of development, as well as to the ethical and technical restrictions that limit the feasibility of observing gastrulation ex vivo 2 . Here we show that human embryonic stem cells can be used to generate gastruloids: three dimensional multicellular aggregates that differentiate to derivatives of the three germ layers organised spatiotemporally, without additional extra-embryonic tissues. Human gastruloids undergo elongation along an anteroposterior axis and, using spatial transcriptomics, we show that they exhibit patterned gene expression. This includes a somitogenesis signature that suggests that 72 hour human gastruloids exhibit features of Carnegie Stage 9 embryos 3 . Our study represents a new, experimentally tractable model system to reveal and probe human-specific regulatory processes occurring during axial organisation in early development.The body plan of mammalian embryos emerges through interactions of sequential cell fate decisions and morphogenetic events, which have hitherto been difficult to observe in humans.Human Embryonic Stem Cells (hESCs) 4 have opened up opportunities for studying early fate decisions, and have hinted at the existence of regulatory mechanisms specific to humans 5,6 . But, in contrast to the embryo, where proportionate populations interact with one another to generate tissues and organs, differentiation in adherent culture is heterogeneous and favours a limited number of cell types 7 . Seeding hESCs on micropatterned surfaces yields coordinated patterns of gene expression, but without the axial organization characteristic of embryos 8 . However, when mouse ESCs are aggregated in suspension under defined conditions, they generate 'gastruloids': a three-dimensional, in vitro model of mammalian development, which exhibits an embryo-like spatiotemporal organization of gene expression 9,10 . We hypothesised that similar human gastruloids could be derived from hESCs. Generation of human gastruloidsWhen hESCs in 2D culture were treated with Chiron, a Wnt agonist, for one day before seeding defined numbers in low-adherence plates in the presence of Chiron, they formed compact, spherical aggregates within a few hours (Fig. 1a, Extended Data Fig. 1a-c). These aggregates progressively broke symmetry and formed elongated structures, with maximal elongation at 72-96h (Fig. 1a-d). On average, ~66% of aggregates from each experiment displayed an elongated morphology at 72h. Although some of the structures remained TOP award (NWO-CW 714.016.001), and the Foundation for Fundamental Research on Matter, financially supported by NWO (FOM-14NOISE01) to S.C.v.d.B., A.A. and A.v.O.. This w...

Gastruloids are three-dimensional aggregates of embryonic stem cells (ESCs) that display key features of mammalian post-implantation development, including germ layer specification and axial organization 1-3. So far, the expression pattern of only a small number of genes in gastruloids has been explored with microscopy, but it is still unclear to what extent genome-wide expression patterns mimic those in embryos. Here, we compared mouse gastruloids with mouse embryos using single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (tomo-seq). We

Highlights d scRNA-seq of the healthy and injured liver epithelium reveals dynamic heterogeneity d Transcriptional heterogeneity is driven by fluctuating activation of YAP signaling d YAP is essential for adult biliary cell survival under homeostatic conditions d YAP is required in hepatocytes for the ductular response during regeneration

Recent advances in non-equilibrium statistical mechanics and single molecule technologies make it possible to extract free energy differences from irreversible work measurements in pulling experiments. To date, free energy recovery has been focused on native or equilibrium molecular states, whereas free energy measurements of kinetic states (i.e. finite lifetime states that are generated dynamically and are metastable) have remained unexplored. Kinetic states can play an important role in various domains of physics, such as nanotechnology or condensed matter physics. In biophysics, there are many examples where they determine the fate of molecular reactions: protein and peptide-nucleic acid binding, specific cation binding, antigen-antibody interactions, transient states in enzymatic reactions or the formation of transient intermediates and non-native structures in molecular folders. Here we demonstrate that it is possible to obtain free energies of kinetic states by applying extended fluctuation relations. This is shown by using optical tweezers to mechanically unfold and refold DNA structures exhibiting intermediate and misfolded kinetic states.Kinetic states are observed under non-equilibrium conditions and have higher free energies than native states. Yet, they can be crucial, as shown by the role that misfolded proteins play in numerous severe diseases [1]. The measurement of the free energy of formation of kinetic states is therefore a central question in biophysics. Recent theoretical developments known as fluctuation relations [13, 12, 4, 5, 6] have been applied to extract free energy differences of equilibrium states from irreversible work measurements. Applications include the measurement of the free energy of formation of RNA and DNA hairpins [7]; the determination of the stability of native domains in proteins [8]; the measurement of mechanical torque in rotary motors [9]; the conversion of information into work in systems under feedback control [10]; or the recovery of free energy landscapes from unidirectional work measurements [11, 12].The characterization of kinetic states under non-equilibrium conditions remains a challenging problem. Here we use a recently introduced extended fluctuation relation (EFR) to extract free energies of kinetic states and thermodynamic branches using irreversible work measurements [13, 10]. In the EFR, a kinetic state is a partially equilibrated region of configurational space, meaning that during a finite timescale the system is confined and thermalized within that region [15]. This is mathematically described by a Boltzmann-Gibbs distribution restricted to configurations contained in that region ( Fig. 1a).Let A, B denote any two kinetic states and λ a control parameter. We consider a forward (F) non-equilibrium process, where the system starts in partial equilibrium in A at λ 0 , and its time-reversed (R), where the partial equilibrium condition is required over B at λ 1 . In the F process λ varies from λ 0 to λ 1 during a time τ according to a predetermined proto...

RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to both monovalent and multivalent ions. Despite recent advances, we do not have general methods to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, and the thermodynamic parameters for secondary structure prediction have only been derived at 1M [Na+]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed monovalent/Mg2+ salt conditions. We measure the unfolding and folding rupture forces and apply Kramers theory to extract accurate information about the hairpin free energy landscape under tension at a wide range of ionic conditions. We obtain non-specific corrections for the free energy of formation of the RNA hairpin and measure how the distance of the transition state to the folded state changes with force and ionic strength. We experimentally validate the Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test the approximate rule by which the non-specific binding affinity of divalent cations at a given concentration is equivalent to that of monovalent cations taken at 100-fold concentration for small molecular constructs.

Most methods for single-cell transcriptome sequencing amplify the termini of polyadenylated transcripts, capturing only a small fraction of the total cellular transcriptome. This precludes the detection of many long non-coding, short non-coding and non-polyadenylated protein-coding transcripts and hinders alternative splicing analysis. We, therefore, developed VASA-seq to detect the total transcriptome in single cells, which is enabled by fragmenting and tailing all RNA molecules subsequent to cell lysis. The method is compatible with both plate-based formats and droplet microfluidics. We applied VASA-seq to more than 30,000 single cells in the developing mouse embryo during gastrulation and early organogenesis. Analyzing the dynamics of the total single-cell transcriptome, we discovered cell type markers, many based on non-coding RNA, and performed in vivo cell cycle analysis via detection of non-polyadenylated histone genes. RNA velocity characterization was improved, accurately retracing blood maturation trajectories. Moreover, our VASA-seq data provide a comprehensive analysis of alternative splicing during mammalian development, which highlighted substantial rearrangements during blood development and heart morphogenesis.

The unfolding and folding of protein barnase has been extensively investigated in bulk conditions under the effect of denaturant and temperature. These experiments provided information about structural and kinetic features of both the native and the unfolded states of the protein, and debates about the possible existence of an intermediate state in the folding pathway have arisen. Here, we investigate the folding/unfolding reaction of protein barnase under the action of mechanical force at the single-molecule level using optical tweezers. We measure unfolding and folding force-dependent kinetic rates from pulling and passive experiments, respectively, and using Kramers-based theories (e.g., Bell-Evans and Dudko-Hummer-Szabo models), we extract the position of the transition state and the height of the kinetic barrier mediating unfolding and folding transitions, finding good agreement with previous bulk measurements. Measurements of the force-dependent kinetic barrier using the continuous effective barrier analysis show that protein barnase verifies the Leffler-Hammond postulate under applied force and allow us to extract its free energy of folding, ΔG0. The estimated value of ΔG0 is in agreement with our predictions obtained using fluctuation relations and previous bulk studies. To address the possible existence of an intermediate state on the folding pathway, we measure the power spectrum of force fluctuations at high temporal resolution (50 kHz) when the protein is either folded or unfolded and, additionally, we study the folding transition-path time at different forces. The finite bandwidth of our experimental setup sets the lifetime of potential intermediate states upon barnase folding/unfolding in the submillisecond timescale.