Search

One of the best known ways bacteria cells understand and respond to the environment are through Two-Component Systems (TCS). These signalling systems are highly diverse in function and can detect a range of physical stimuli including molecular concentrations and temperature, with a range of responses including chemotaxis and anaerobic energy production.TCS exhibit a range of different molecular structures and energy costs, and multiple types co-exist in the same cell. TCSs that incur relatively high energy cost are abundant in biology, despite strong evolutionary pressure to efficiently spend energy.We are motivated to discern what benefits, if any, the more energetically expensive variants had for a cell.We seek to answer this question by modelling energy flow through two variants of TCS. This was accomplished using bond graphs, a physics-based modelling framework that accurately models energy transfer through different physical domains. Our analysis demonstrates that energy availability can affect a cell’s signal sensitivity, noise filtering effectiveness, and the stimulus level where cell response is maximal. We also found that these properties are determined not by the molecular parameters themselves, but the reaction rate parameters that govern the reaction systems as a whole.This suggests possible connections between the molecular structure and evolutionary purpose of any two-component system. This opens the door to new synthetic circuit design in systems biology, and we propose new hypotheses about this link between structure and purpose that could be experimentally verified.Author summaryTwo-component systems are the main way many bacteria sense and respond to their environment. They exist in such well-studied bacteria asE. coliwhere they have been shown to detect a range of stimuli including nutrients, temperature, acidity, and pressure.Two-component systems are ubiquitous in bacteria yet have a deceptively simple structure. Knowing how they operate and the purpose of variations in signalling structure is helpful to our understanding of cellular biology and the design of synthetic biological circuits. Critical unanswered questions remain about the energy usage and functional benefits of these systems.We sought to improve our understanding of two-component systems by applying a physics-based modelling framework. We found that tracking energy flow through the cell reveals new energy-dependent behaviour in signalling sensitivity, noise filtering, and maximal cell response. We also found that these properties are not strictly dependent on the molecular properties themselves, but from the configuration of the reaction system as a whole.

Background In many cancers, microRNAs (miRs) contribute to metastatic progression by modulating phenotypic reprogramming processes such as epithelial-mesenchymal plasticity. This can be driven by miRs targeting multiple mRNA transcripts, inducing regulated changes across large sets of genes. The miR-target databases TargetScan and DIANA-microT predict putative relationships by examining sequence complementarity between miRs and mRNAs. However, it remains a challenge to identify which miR-mRNA interactions are active at endogenous expression levels, and of biological consequence. Methods We developed a workflow to integrate TargetScan and DIANA-microT predictions into the analysis of data-driven associations calculated from transcript abundance (RNASeq) data, specifically the mutual information and Pearson’s correlation metrics. We use this workflow to identify putative relationships of miR-mediated mRNA repression with strong support from both lines of evidence. Applying this approach systematically to a large, published collection of unique melanoma cell lines – the Ludwig Melbourne melanoma (LM-MEL) cell line panel – we identified putative miR-mRNA interactions that may contribute to invasiveness. This guided the selection of interactions of interest for further in vitro validation studies. Results Several miR-mRNA regulatory relationships supported by TargetScan and DIANA-microT demonstrated differential activity across cell lines of varying matrigel invasiveness. Strong negative statistical associations for these putative regulatory relationships were consistent with target mRNA inhibition by the miR, and suggest that differential activity of such miR-mRNA relationships contribute to differences in melanoma invasiveness. Many of these relationships were reflected across the skin cutaneous melanoma TCGA dataset, indicating that these observations also show graded activity across clinical samples. Several of these miRs are implicated in cancer progression (miR-211, -340, -125b, −221, and -29b). The specific role for miR-29b-3p in melanoma has not been well studied. We experimentally validated the predicted miR-29b-3p regulation of LAMC1 and PPIC and LASP1, and show that dysregulation of miR-29b-3p or these mRNA targets can influence cellular invasiveness in vitro. Conclusions This analytic strategy provides a comprehensive, systems-level approach to identify miR-mRNA regulation in high-throughput cancer data, identifies novel putative interactions with functional phenotypic relevance, and can be used to direct experimental resources for subsequent experimental validation. Computational scripts are available: http://github.com/uomsystemsbiology/LMMEL-miR-miner Electronic supplementary material The online version of this article (doi:10.1186/s12943-016-0554-y) contains supplementary material, which is available to authorized users.

An early step of evaluating a nanomaterial’s potential for biological applications is investigating its interactions with cells and biological fluids. These experiments generate complex data, which...

Understanding the interactions between nano-engineered particles and cells is necessary for the rational design of particles for therapeutic, diagnostic and imaging purposes. In particular, the informed design of particles relies on the quantification of the relationship between the physicochemical properties of the particles and the rate at which cells interact with, and subsequently internalise, particles. Quantitative models, both mathematical and computational, provide a powerful tool for elucidating this relationship, as well as for understanding the mechanisms governing the intertwined processes of interaction and internalisation. Here we review the different types of mathematical and computational models that have been used to examine particle-cell interactions and particle internalisation. We detail the mathematical methodology for each type of model, the benefits and limitations associated with the different types of models, and highlight the advances in understanding gleaned from the application of these models to experimental observations of particle internalisation. We discuss the recent proposal and ongoing community adoption of standardised experimental reporting, and how this adoption is an important step toward unlocking the full potential of modelling approaches. Finally, we consider future directions in quantitative models of particle-cell interactions and highlight the need for hybrid experimental and theoretical investigations to address hitherto unanswered questions.

Calcium (Ca

The Systems Biology Markup Language (SBML) is a popular software-independent XML-based format for describing models of biological phenomena. The BioModels Database is the largest online repository of SBML models. Several tools and platforms are available to support the reuse and composition of SBML models. However, these tools do not explicitly assess whether models are physically plausibile or thermodynamically consistent. This often leads to ill-posed models that are physically impossible, impeding the development of realistic complex models in biology. Here, we present a framework that can automatically convert SBML models into bond graphs, which imposes energy conservation laws on these models. The new bond graph models are easily mergeable, resulting in physically plausible coupled models. We illustrate this by automatically converting and coupling a model of pyruvate distribution to a model of the pentose phosphate pathway.Graphical AbstractHighlightsA framework to convert suitable SBML models of biochemical networks into bond graphs is developed.The framework is applied here to two interconnecting models of metabolism pathways.We automatically integrate the generated bond graph modules.We qualitatively illustrate the functionality of the composed model.

The maintenance of tissue and organ structures during dynamic homeostasis is often not well understood. In order for a system to be stable, cell renewal, cell migration and cell death must be finely balanced. Moreover, a tissue's shape must remain relatively unchanged. Simple epithelial tissues occur in various structures throughout the body, such as the endothelium, mesothelium, linings of the lungs, saliva and thyroid glands, and gastrointestinal tract. Despite the prevalence of simple epithelial tissues, there are few models which accurately describe how these tissues maintain a stable structure. Here, we present a novel, 3D, deformable, multilayer, cell-centre model of a simple epithelium. Cell movement is governed by the minimisation of a bending potential across the epithelium, cell-cell adhesion, and viscous effects. We show that the model is capable of maintaining a consistent tissue structure while undergoing self renewal. We also demonstrate the model's robustness under tissue renewal, cell migration and cell removal. The model presented here is a valuable advancement towards the modelling of tissues and organs with complex and generalised structures.

Calcium (Ca2＋) plays a critical role in the excitation contraction coupling (ECC) process that mediates the contraction of cardiomyocytes during each heartbeat. While ryanodine receptors (RyRs) are the primary Ca2＋ channels responsible for generating the cell-wide Ca2＋ transients during ECC, Ca2＋ release, via inositol 1,4,5-trisphosphate (IP3) receptors (IP3Rs) are also reported in cardiomyocytes to elicit ECC-modulating effects. Recent studies suggest that the localization of IP3Rs at dyads grant their ability to modify the occurrence of Ca2＋ sparks (elementary Ca2＋ release events that constitute cell wide Ca2＋ releases associated with ECC) which may underlie their modulatory influence on ECC. Here, we aim to uncover the mechanism by which dyad-localized IP3Rs influence Ca2＋ spark dynamics. To this end, we developed a mathematical model of the dyad that incorporates the behaviour of IP3Rs, in addition to RyRs, to reveal the impact of their activity on local Ca2＋ handling and consequent Ca2＋ spark occurrence and its properties. Consistent with published experimental data, our model predicts that the propensity for Ca2＋ spark formation increases in the presence of IP3R activity. Our simulations support the hypothesis that IP3Rs elevate Ca2＋ in the dyad, sensitizing proximal RyRs towards activation and hence Ca2＋ spark formation. The stochasticity of IP3R gating is an important aspect of this mechanism. However, dyadic IP3R activity lowers the Ca2＋ available in the junctional sarcoplasmic reticulum (JSR) for release, thus resulting in Ca2＋ sparks with similar durations but lower amplitudes. ispartof: Mathematical Biosciences vol:355 ispartof: location:United States status: published

Like all physical systems, biological systems are constrained by the laws of physics. However, mathematical models of biochemistry frequently neglect the conservation of energy, leading to unrealistic behaviour. Energy-based models that are consistent with conservation of mass, charge and energy have the potential to aid the understanding of complex interactions between biological components, and are becoming easier to develop with recent advances in experimental measurements and databases. In this paper, we motivate the use of bond graphs (a modelling tool from engineering) for energy-based modelling and introduce, BondGraphTools, a Python library for constructing and analysing bond graph models. We use examples from biochemistry to illustrate how BondGraphTools can be used to automate model construction in systems biology while maintaining consistency with the laws of physics.

Hierarchical modelling is essential to achieving complex, large-scale models. However, not all modelling schemes support hierarchical composition, and correctly mapping points of connection between models requires comprehensive knowledge of each model’s components and assumptions. To address these challenges in integrating biosimulation models, we propose an approach to automatically and confidently compose biosimulation models. The approach uses bond graphs to combine aspects of physical and thermodynamics-based modelling with biological semantics. We improved on existing approaches by using semantic annotations to automate the recognition of common components. The approach is illustrated by coupling a model of the Ras-MAPK cascade to a model of the upstream activation of EGFR. Through this methodology, we aim to assist researchers and modellers in readily having access to more comprehensive biological systems models.Author summaryDetailed, multi-scale computational models bridging from biomolecular processes to entire organs and bodies have the potential to revolutionise medicine by enabling personalised treatments. One of the key challenges to achieving these models is connecting together the vast number of isolated biosimulation models into a coherent whole. Using recent advances in both modelling techniques and biological standards in the scientific community, we developed an approach to integrate and compose models in a physics-based environment. This provides significant advantages, including the automation of model composition and post-model-composition adjustments. We anticipate that our approach will enable the faster development of realistic and accurate models to understand complex biological systems.

Nanoparticles hold great preclinical promise in cancer therapy but continue to suffer attrition through clinical trials. Advanced, three dimensional (3D) cellular models such as tumor spheroids can recapitulate elements of the tumor environment and are considered the superior model to evaluate nanoparticle designs. However, there is an important need to better understand nanoparticle penetration kinetics and determine how different cell characteristics may influence this nanoparticle uptake. A key challenge with current approaches for measuring nanoparticle accumulation in spheroids is that they are often static, losing spatial and temporal information which may be necessary for effective nanoparticle evaluation in 3D cell models. To overcome this challenge, we developed an analysis platform, termed the Determination of Nanoparticle Uptake in Tumor Spheroids (DONUTS), which retains spatial and temporal information during quantification, enabling evaluation of nanoparticle uptake in 3D tumor spheroids. Outperforming linear profiling methods, DONUTS was able to measure silica nanoparticle uptake to 10 μm accuracy in both isotropic and irregularly shaped cancer cell spheroids. This was then extended to determine penetration kinetics, first by a forward-in-time, center-in-space model, and then by mathematical modelling, which enabled the direct evaluation of nanoparticle penetration kinetics in different spheroid models. Nanoparticle uptake was shown to inversely relate to particle size and varied depending on the cell type, cell stiffness and density of the spheroid model. The automated analysis method we have developed can be applied to live spheroids in situ, for the advanced evaluation of nanoparticles as delivery agents in cancer therapy.

BackgroundThe advent of next-generation sequencing revealed extensive transcription beyond protein-coding genes, identifying tens of thousands of long non-coding RNAs (lncRNAs). Selected functional examples raised the possibility that lncRNAs, as a class, may maintain broad regulatory roles. Compellingly, lncRNA expression is strongly linked with adjacent protein-coding gene expression, suggesting a potential cis-regulatory function. Evidence for these regulatory roles may be obtained through careful examination of the precise timing of lncRNA expression relative to adjacent protein-coding genes.ResultsWhere causal cis-regulatory relationships exist, lncRNA activation is expected to precede changes in adjacent target gene expression. Using an RNA-seq time course of uniquely high temporal resolution, we profiled the expression dynamics of several thousand lncRNAs and protein-coding genes in synchronized, transitioning human cells. Our findings reveal lncRNAs are expressed synchronously with adjacent protein-coding genes. Analysis of lipopolysaccharide-activated mouse dendritic cells revealed the same temporal relationship observed in transitioning human cells.ConclusionOur findings suggest broad-scale cis-regulatory roles for lncRNAs are not common. The strong association between lncRNAs and adjacent genes may instead indicate an origin as transcriptional by-products from active protein-coding gene promoters and enhancers.

Simulating complex biological and physiological systems and predicting their behaviours under different conditions remains challenging. Breaking systems into smaller and more manageable modules can address this challenge, assisting both model development and simulation. Nevertheless, existing computational models in biology and physiology are often not modular and therefore difficult to assemble into larger models. Even when this is possible, the resulting model may not be useful due to inconsistencies either with the laws of physics or the physiological behaviour of the system. Here, we propose a general methodology for composing models, combining the energy-based bond graph approach with semantics-based annotations. This approach improves model composition and ensures that a composite model is physically plausible. As an example, we demonstrate this approach to automated model composition using a model of human arterial circulation. The major benefit is that modellers can spend more time on understanding the behaviour of complex biological and physiological systems and less time wrangling with model composition., Author summary Biological and physiological systems usually involve multiple underlying processes, mechanisms, structures, and phenomena, referred to here as sub-systems. Modelling the whole system every time from scratch requires a huge amount of effort. An alternative is to model each sub-system in a modular fashion, i.e., containing meaningful interfaces for connecting to other modules. Such modules are readily combined to produce a whole-system model. For the combined model to be consistent, modules must be described using the same modelling scheme. One way to achieve this is to use energy-based models that are consistent with the conservation laws of physics. Here, we present an approach that achieves this using bond graphs, which allows modules to be combined faster and more efficiently. First, physically plausible modules are generated using a small number of template modules. Then a meaningful interface is added to each module to automate connection. This approach is illustrated by applying this method to an existing model of the circulatory system and verifying the results against the reference model.

Like all physical systems, biological systems are constrained by the laws of physics. However, mathematical models of biochemistry frequently neglect the conservation of energy, leading to unrealistic behaviour. Energy-based models that are consistent with conservation of mass, charge and energy have the potential to aid the understanding of complex interactions between biological components, and are becoming easier to develop with recent advances in experimental measurements and databases. In this paper, we motivate the use of bond graphs (a modelling tool from engineering) for energy-based modelling and introduce, BondGraphTools, a Python library for constructing and analysing bond graph models. We use examples from biochemistry to illustrate how BondGraphTools can be used to automate model construction in systems biology while maintaining consistency with the laws of physics.

Calcium (Ca2+) plays a central role in mediating both contractile function and hypertrophic signaling in ventricular cardiomyocytes. L-type Ca2+ channels trigger release of Ca2+ from ryanodine receptors for cellular contraction, whereas signaling downstream of G-protein-coupled receptors stimulates Ca2+ release via inositol 1,4,5-trisphosphate receptors (IP3Rs), engaging hypertrophic signaling pathways. Modulation of the amplitude, duration, and duty cycle of the cytosolic Ca2+ contraction signal and spatial localization have all been proposed to encode this hypertrophic signal. Given current knowledge of IP3Rs, we develop a model describing the effect of functional interaction (cross talk) between ryanodine receptor and IP3R channels on the Ca2+ transient and examine the sensitivity of the Ca2+ transient shape to properties of IP3R activation. A key result of our study is that IP3R activation increases Ca2+ transient duration for a broad range of IP3R properties, but the effect of IP3R activation on Ca2+ transient amplitude is dependent on IP3 concentration. Furthermore we demonstrate that IP3-mediated Ca2+ release in the cytosol increases the duty cycle of the Ca2+ transient, the fraction of the cycle for which [Ca2+] is elevated, across a broad range of parameter values and IP3 concentrations. When coupled to a model of downstream transcription factor (NFAT) activation, we demonstrate that there is a high correspondence between the Ca2+ transient duty cycle and the proportion of activated NFAT in the nucleus. These findings suggest increased cytosolic Ca2+ duty cycle as a plausible mechanism for IP3-dependent hypertrophic signaling via Ca2+-sensitive transcription factors such as NFAT in ventricular cardiomyocytes. ispartof: BIOPHYSICAL JOURNAL vol:119 issue:6 pages:1178-1192 ispartof: location:United States status: published

The Na$^+$/K$^+$ ATPase is an essential component of cardiac electrophysiology, maintaining physiological Na$^+$ and K$^+$ concentrations over successive heart beats. Terkildsen et al. (2007) developed a model of the ventricular myocyte Na$^+$/K$^+$ ATPase to study extracellular potassium accumulation during ischaemia, demonstrating the ability to recapitulate a wide range of experimental data, but unfortunately there was no archived code associated with the original manuscript. Here we detail an updated version of the model and provide CellML and MATLAB code to ensure reproducibility and reusability. We note some errors within the original formulation which have been corrected to ensure that the model is thermodynamically consistent, and although this required some reparameterisation, the resulting model still provides a good fit to experimental measurements that demonstrate the dependence of Na$^+$/K$^+$ ATPase pumping rate upon membrane voltage and metabolite concentrations. To demonstrate thermodynamic consistency we also developed a bond graph version of the model. We hope that these models will be useful for community efforts to assemble a whole-cell cardiomyocyte model which facilitates the investigation of cellular energetics.

Calcium plays critical roles in cardiac cells, coupling electrical excitation to mechanical contraction with each heartbeat, while simultaneously mediating biochemical signals that regulate cell growth. While ryanodine receptors (RyRs) are fundamental to generation of elementary calcium release events (sparks) and global calcium elevations that underlie excitation-contraction coupling (ECC), calcium release via inositol 1,4,5-trisphosphate receptors (IP 3 Rs) is also reported in cardiomyocytes. IP 3 R calcium release modifies ECC as well as contributing to downstream regulation of hypertrophic gene expression. Recent studies suggest that proximal localisation of IP 3 Rs with RyRs contributes to their ability to modify Ca 2+ handling during ECC. Here we aim to determine the mechanism by which IP 3 Rs modify Ca 2+ handling in cardiomyocytes. We develop a mathematical model incorporating the stochastic behaviour of receptor opening that allows for the parametric tuning of the system to reveal the impact of IP 3 Rs on spark activation. By testing multiple spark initiation mechanisms, we find that Ca 2+ release via IP 3 Rs result in increased propensity for spark initiation within the cardiac dyad. Our simulations suggest that opening of IP 3 Rs elevates Ca 2+ within the dyad, which increase the probability of spark initiation. Finally, we find that while increasing the number of IP 3 Rs increases the probability of spark formation, it has little effect on spark amplitude, duration, or overall shape. Our study therefore suggests that IP 3 R play a critical role in modulating Ca 2+ signaling for excitation contraction coupling Author summary While Ca 2+ release through ryanodine receptors (RyRs) initiates contraction in cardiomyocytes, the role of inositol 1,4,5-trisphosphate receptors (IP 3 Rs) in cardiomyocytes is less clear with Ca 2+ release through these channels being invoked in regulating ECC and hypertrophic signalling. RyRs generate cytosolic Ca 2+ signals through elemental Ca 2+ release events called sparks. The mechanisms by which IP 3 Rs influence cytosolic Ca 2+ are not well understood. We created a 1D model of calcium spark formation in a cardiomyocyte dyad—the primary site of elemental RyR-based calcium release. We investigated possible behaviours of IP 3 Rs and their interaction with RyRs in generating Ca 2+ sparks. We show that for high IP 3 concentration, a large number of IP 3 Rs and high IP 3 R affinity are required to noticeably affect spark shape. At lower IP 3 concentration IP 3 Rs can increase Ca 2+ spark activity, but do not significantly alter the spark shape. Finally our simulations suggest that spark frequency can be reliably increased when IP 3 Rs activity is such that a small continuous Ca 2+ flux is introduced to the dyad to elevate Ca 2+ , and not via brief but high Ca 2+ release from these receptors.

In experimental studies on cardiac tissue, the end-systolic force-length relation (ESFLR) has been shown to depend on the mode of contraction: isometric or isotonic. The isometric ESFLR is derived from isometric contractions spanning a range of muscle lengths while the isotonic ESFLR is derived from shortening contractions across a range of afterloads. The ESFLR of isotonic contractions consistently lies below its isometric counterpart. Despite the passing of over a hundred years since the first insight by Otto Frank, the mechanism(s) underlying this protocol-dependent difference in the ESFLR remain incompletely explained. Here, we investigate the role of mechano-calcium feedback in accounting for the difference between these two ESFLRs. Previous studies have compared the dynamics of isotonic contractions to those of a single isometric contraction at a length that produces maximum force, without considering isometric contractions at shorter muscle lengths. We used a mathematical model of cardiac excitation-contraction to simulate isometric and force-length work-loop contractions (the latter being the 1D equivalent of the whole-heart pressure-volume loop), and compared Ca2+ transients produced under equivalent force conditions. We found that the duration of the simulated Ca2+ transient increases with decreasing sarcomere length for isometric contractions, and increases with decreasing afterload for work-loop contractions. At any given force, the Ca2+ transient for an isometric contraction is wider than that during a work-loop contraction. By driving simulated work-loops with wider Ca2+ transients generated from isometric contractions, we show that the duration of muscle shortening was prolonged, thereby shifting the work-loop ESFLR toward the isometric ESFLR. These observations are explained by an increase in the rate of binding of Ca2+ to troponin-C with increasing force. However, the leftward shift of the work-loop ESFLR does not superimpose on the isometric ESFLR, leading us to conclude that while mechano-calcium feedback does indeed contribute to the difference between the two ESFLRs, it does not completely account for it.

Calcium (Ca(2+)) plays a central role in mediating both contractile function and hypertrophic signaling in ventricular cardiomyocytes. L-type Ca(2+) channels trigger release of Ca(2+) from ryanodine receptors for cellular contraction, whereas signaling downstream of G-protein-coupled receptors stimulates Ca(2+) release via inositol 1,4,5-trisphosphate receptors (IP(3)Rs), engaging hypertrophic signaling pathways. Modulation of the amplitude, duration, and duty cycle of the cytosolic Ca(2+) contraction signal and spatial localization have all been proposed to encode this hypertrophic signal. Given current knowledge of IP(3)Rs, we develop a model describing the effect of functional interaction (cross talk) between ryanodine receptor and IP(3)R channels on the Ca(2+) transient and examine the sensitivity of the Ca(2+) transient shape to properties of IP(3)R activation. A key result of our study is that IP(3)R activation increases Ca(2+) transient duration for a broad range of IP(3)R properties, but the effect of IP(3)R activation on Ca(2+) transient amplitude is dependent on IP(3) concentration. Furthermore we demonstrate that IP(3)-mediated Ca(2+) release in the cytosol increases the duty cycle of the Ca(2+) transient, the fraction of the cycle for which [Ca(2+)] is elevated, across a broad range of parameter values and IP(3) concentrations. When coupled to a model of downstream transcription factor (NFAT) activation, we demonstrate that there is a high correspondence between the Ca(2+) transient duty cycle and the proportion of activated NFAT in the nucleus. These findings suggest increased cytosolic Ca(2+) duty cycle as a plausible mechanism for IP(3)-dependent hypertrophic signaling via Ca(2+)-sensitive transcription factors such as NFAT in ventricular cardiomyocytes.

The question of whether a population will persist or go extinct is of key interest throughout ecology and biology. Various mathematical techniques allow us to generate knowledge regarding individual behaviour, which can be analysed to obtain predictions about the ultimate survival or extinction of the population. A common model employed to describe population dynamics is the lattice-based random walk model with crowding (exclusion). This model can incorporate behaviour such as birth, death and movement, while including natural phenomena such as finite size effects. Performing sufficiently many realizations of the random walk model to extract representative population behaviour is computationally intensive. Therefore, continuum approximations of random walk models are routinely employed. However, standard continuum approximations are notoriously incapable of making accurate predictions about population extinction. Here, we develop a new continuum approximation, the state-space diffusion approximation , which explicitly accounts for population extinction. Predictions from our approximation faithfully capture the behaviour in the random walk model, and provides additional information compared to standard approximations. We examine the influence of the number of lattice sites and initial number of individuals on the long-term population behaviour, and demonstrate the reduction in computation time between the random walk model and our approximation.

Summary: A thorough understanding of cellular development is incumbent on assessing the complexities of fate and kinetics of individual clones within a population. Here, we develop a system for robust periodical assessment of lineage outputs of thousands of transient clones and establishment of bona fide cellular trajectories. We appraise the development of dendritic cells (DCs) in fms-like tyrosine kinase 3 ligand culture from barcode-labeled hematopoietic stem and progenitor cells (HSPCs) by serially measuring barcode signatures and visualize these multidimensional data using developmental interpolated t-distributed stochastic neighborhood embedding (DiSNE) time-lapse movies. We identify multiple cellular trajectories of DC development that are characterized by distinct fate bias and expansion kinetics and determine that these are intrinsically programmed. We demonstrate that conventional DC and plasmacytoid DC trajectories are largely separated already at the HSPC stage. This framework allows systematic evaluation of clonal dynamics and can be applied to other steady-state or perturbed developmental systems. : Lin et al. develop a framework to longitudinally track and visualize clonal cellular trajectories during dendritic cell development. The authors demonstrate that properties including fate bias and the timing and size of clonal bursts are heterogeneous within populations yet largely imprinted in single progenitors at early developmental stages. Keywords: dendritic cells, hematopoiesis, single cell, cellular barcoding, t-SNE, DiSNE, lineage, fate, kinetics, cellular trajectory