No Arabic abstract
Numerous biological approaches are available to characterise the mechanisms which govern the formation of human embryonic stem cell (hESC) colonies. To understand how the kinematics of single and pairs of hESCs impact colony formation, we study their mobility characteristics using time-lapse imaging. We perform a detailed statistical analysis of their speed, survival, directionality, distance travelled and diffusivity. We confirm that single and pairs of cells migrate as a diffusive random walk. Moreover, we show that the presence of Cell Tracer significantly reduces hESC mobility. Our results open the path to employ the theoretical framework of the diffusive random walk for the prognostic modelling and optimisation of the growth of hESC colonies. Indeed, we employ this random walk model to estimate the seeding density required to minimise the occurrence of hESC colonies arising from more than one founder cell and the minimal cell number needed for successful colony formation. We believe that our prognostic model can be extended to investigate the kinematic behaviour of somatic cells emerging from hESC differentiation and to enable its wide application in phenotyping of pluripotent stem cells for large scale stem cell culture expansion and differentiation platforms.
The deluge of single-cell data obtained by sequencing, imaging and epigenetic markers has led to an increasingly detailed description of cell state. However, it remains challenging to identify how cells transition between different states, in part because data are typically limited to snapshots in time. A prerequisite for inferring cell state transitions from such snapshots is to distinguish whether transitions are coupled to cell divisions. To address this, we present two minimal branching process models of cell division and differentiation in a well-mixed population. These models describe dynamics where differentiation and division are coupled or uncoupled. For each model, we derive analytic expressions for each subpopulations mean and variance and for the likelihood, allowing exact Bayesian parameter inference and model selection in the idealised case of fully observed trajectories of differentiation and division events. In the case of snapshots, we present a sample path algorithm and use this to predict optimal temporal spacing of measurements for experimental design. We then apply this methodology to an textit{in vitro} dataset assaying the clonal growth of epiblast stem cells in culture conditions promoting self-renewal or differentiation. Here, the larger number of cell states necessitates approximate Bayesian computation. For both culture conditions, our inference supports the model where cell state transitions are coupled to division. For culture conditions promoting differentiation, our analysis indicates a possible shift in dynamics, with these processes becoming more coupled over time.
A quantitative description of the flagellar dynamics in the procyclic T. brucei is presented in terms of stationary oscillations and traveling waves. By using digital video microscopy to quantify the kinematics of trypanosome flagellar waveforms. A theoretical model is build starting from a Bernoulli-Euler flexural-torsional model of an elastic string with internal distribution of force and torque. The dynamics is internally driven by the action of the molecular motors along the string, which is proportional to the local shift and consequently to the local curvature. The model equation is a nonlinear partial differential wave equation of order four, containing nonlinear terms specific to the Korteweg-de Vries (KdV) equation and the modified-KdV equation. For different ranges of parameters we obtained kink-like solitons, breather solitons, and a new class of solutions constructed by smoothly piece-wise connected conic functions arcs (e.g. ellipse). The predicted amplitude and wavelengths are in good match with experiments. We also present the hypotheses for a step-wise kinematical model of swimming of procyclic African trypanosome.
Bacteria such as Escherichia coli move about in a series of runs and tumbles: while a run state (straight motion) entails all the flagellar motors spinning in counterclockwise mode, a tumble is caused by a shift in the state of one or more motors to clockwise spinning mode. In the presence of an attractant gradient in the environment, runs in the favourable direction are extended, and this results in a net drift of the organism in the direction of the gradient. The underlying signal transduction mechanism produces directed motion through a bi-lobed response function which relates the clockwise bias of the flagellar motor to temporal changes in the attractant concentration. The two lobes (positive and negative) of the response function are separated by a time interval of $sim 1$s, such that the bacterium effectively compares the concentration at two different positions in space and responds accordingly. We present here a novel path-integral method which allows us to address this problem in the most general way possible, including multi-step CW-CCW transitions, directional persistence and power-law waiting time distributions. The method allows us to calculate quantities such as the effective diffusion coefficient and drift velocity, in a power series expansion in the attractant gradient. Explicit results in the lowest order in the expansion are presented for specific models, which, wherever applicable, agree with the known results. New results for gamma-distributed run interval distributions are also presented.
Antimicrobial resistance is an emerging global health crisis that is undermining advances in modern medicine and, if unmitigated, threatens to kill 10 million people per year worldwide by 2050. Research over the last decade has demonstrated that the differences between genetically identical cells in the same environment can lead to drug resistance. Fluctuations in gene expression, modulated by gene regulatory networks, can lead to non-genetic heterogeneity that results in the fractional killing of microbial populations causing drug therapies to fail; this non-genetic drug resistance can enhance the probability of acquiring genetic drug resistance mutations. Mathematical models of gene networks can elucidate general principles underlying drug resistance, predict the evolution of resistance, and guide drug resistance experiments in the laboratory. Cells genetically engineered to carry synthetic gene networks regulating drug resistance genes allow for controlled, quantitative experiments on the role of non-genetic heterogeneity in the development of drug resistance. In this perspective article, we emphasize the contributions that mathematical, computational, and synthetic gene network models play in advancing our understanding of antimicrobial resistance to discover effective therapies against drug-resistant infections.
The growth of several biological tissues is known to be controlled in part by local geometrical features, such as the curvature of the tissue interface. This control leads to changes in tissue shape that in turn can affect the tissues evolution. Understanding the cellular basis of this control is highly significant for bioscaffold tissue engineering, the evolution of bone microarchitecture, wound healing, and tumour growth. While previous models have proposed geometrical relationships between tissue growth and curvature, the role of cell density and cell vigor remains poorly understood. We propose a cell-based mathematical model of tissue growth to investigate the systematic influence of curvature on the collective crowding or spreading of tissue-synthesising cells induced by changes in local tissue surface area during the motion of the interface. Depending on the strength of diffusive damping, the model exhibits complex growth patterns such as undulating motion, efficient smoothing of irregularities, and the generation of cusps. We compare this model with in-vitro experiments of tissue deposition in bioscaffolds of different geometries. By accounting for the depletion of active cells, the model is able to capture both smoothing of initial substrate geometry and tissue deposition slowdown as observed experimentally.