Towards a digital model of zebrafish embryogenesis. Integration of Cell Tracking and Gene Expression Quantification

We elaborate on a general framework composed of a set of computational tools to accurately quantificate cellular position and gene expression levels throughout early zebrafish embryogenesis captured over a time-lapse series of in vivo 3D images. Our modeling strategy involves nuclei detection, cell geometries extraction, automatic gene levels quantification and cell tracking to reconstruct cell trajectories and lineage tree which describe the animal development. Each cell in the embryo is then precisely described at each given time t by a vector composed of the cell 3D spatial coordinates (x; y; z) along with its gene expression level g. This comprehensive description of the embryo development is used to assess the general connection between genetic expression and cell movement. We also investigate genetic expression propagation between a cell and its progeny in the lineage tree. More to the point, this paper focuses on the evolution of the expression pattern of transcriptional factor goosecoid (gsc) through the gastrulation process between 6 and 9 hours post fertilization (hpf

Efficient Algorithms for Moral Lineage Tracing

Lineage tracing, the joint segmentation and tracking of living cells as they move and divide in a sequence of light microscopy images, is a challenging task. Jug et al. have proposed a mathematical abstraction of this task, the moral lineage tracing problem (MLTP), whose feasible solutions define both a segmentation of every image and a lineage forest of cells. Their branch-and-cut algorithm, however, is prone to many cuts and slow convergence for large instances. To address this problem, we make three contributions: (i) we devise the first efficient primal feasible local search algorithms for the MLTP, (ii) we improve the branch-and-cut algorithm by separating tighter cutting planes and by incorporating our primal algorithms, (iii) we show in experiments that our algorithms find accurate solutions on the problem instances of Jug et al. and scale to larger instances, leveraging moral lineage tracing to practical significance.Comment: Accepted at ICCV 201

A Novel Validation Algorithm Allows for Automated Cell Tracking and the Extraction of Biologically Meaningful Parameters

Automated microscopy is currently the only method to non-invasively and label-free observe complex multi-cellular processes, such as cell migration, cell cycle, and cell differentiation. Extracting biological information from a time-series of micrographs requires each cell to be recognized and followed through sequential microscopic snapshots. Although recent attempts to automatize this process resulted in ever improving cell detection rates, manual identification of identical cells is still the most reliable technique. However, its tedious and subjective nature prevented tracking from becoming a standardized tool for the investigation of cell cultures. Here, we present a novel method to accomplish automated cell tracking with a reliability comparable to manual tracking. Previously, automated cell tracking could not rival the reliability of manual tracking because, in contrast to the human way of solving this task, none of the algorithms had an independent quality control mechanism; they missed validation. Thus, instead of trying to improve the cell detection or tracking rates, we proceeded from the idea to automatically inspect the tracking results and accept only those of high trustworthiness, while rejecting all other results. This validation algorithm works independently of the quality of cell detection and tracking through a systematic search for tracking errors. It is based only on very general assumptions about the spatiotemporal contiguity of cell paths. While traditional tracking often aims to yield genealogic information about single cells, the natural outcome of a validated cell tracking algorithm turns out to be a set of complete, but often unconnected cell paths, i.e. records of cells from mitosis to mitosis. This is a consequence of the fact that the validation algorithm takes complete paths as the unit of rejection/acceptance. The resulting set of complete paths can be used to automatically extract important biological parameters with high reliability and statistical significance. These include the distribution of life/cycle times and cell areas, as well as of the symmetry of cell divisions and motion analyses. The new algorithm thus allows for the quantification and parameterization of cell culture with unprecedented accuracy. To evaluate our validation algorithm, two large reference data sets were manually created. These data sets comprise more than 320,000 unstained adult pancreatic stem cells from rat, including 2592 mitotic events. The reference data sets specify every cell position and shape, and assign each cell to the correct branch of its genealogic tree. We provide these reference data sets for free use by others as a benchmark for the future improvement of automated tracking methods

Spatiotemporal Identification of Cell Divisions Using Symmetry Properties in Time-Lapse Phase Contrast Microscopy

A variety of biological and pharmaceutical studies, such as for anti-cancer drugs, require the quantification of cell responses over long periods of time. This is performed with time-lapse video microscopy that gives a long sequence of frames. For this purpose, phase contrast imaging is commonly used since it is minimally invasive. The cell responses of interest in this study are the mitotic cell divisions. Their manual measurements are tedious, subjective, and restrictive. This study introduces an automated method for these measurements. The method starts with preprocessing for restoration and reconstruction of the phase contrast time-lapse sequences. The data are first restored from intensity non-uniformities. Subsequently, the circular symmetry of the contour of the mitotic cells in phase contrast images is used by applying a Circle Hough Transform (CHT) to reconstruct the entire cells. The CHT is also enhanced with the ability to “vote” exclusively towards the center of curvature. The CHT image sequence is then registered for misplacements between successive frames. The sequence is subsequently processed to detect cell centroids in individual frames and use them as starting points to form spatiotemporal trajectories of cells along the positive as well as along the negative time directions, that is, anti-causally. The connectivities of different trajectories enhanced by the symmetry of the trajectories of the daughter cells provide as topological by-products the events of cell divisions together with the corresponding entries into mitoses as well as exits from cytokineses. The experiments use several experimental video sequences from three different cell lines with many cells undergoing mitoses and divisions. The quantitative validations of the results of the processing demonstrate the high performance and efficiency of the method

Cellular tracking in time-lapse phase contrast images

The quantitative analysis of live cells is a key issue in evaluating biological processes. The current clinical practice involves the application of a tedious and time consuming manual tracking procedure on large amount of data. As a result, automatic tracking systems are currently developed and evaluated. However, problems caused by cellular division, agglomeration, Brownian motion and topology changes are difficult issues that have to be accommodated by automatic tracking techniques. In this paper, we detail the development of a fully automated multi-target tracking system that is able to deal with Brownian motion and cellular division. During the tracking process our approach includes the neighbourhood relationship and motion history to enforce the cellular tracking continuity in the spatial and temporal domain. The experimental results reported in this paper indicate that our method is able to accurately track cellular structures in time-lapse data

Dynamics of living cells in a cytomorphological state space.

Cells are nonequilibrium systems that exchange matter and energy with the environment to sustain their metabolic needs. The nonequilibrium nature of this system presents considerable challenges to developing a general theory describing its behavior; however, when studied at appropriate spatiotemporal scales, the behavior of ensembles of nonequilibrium systems can resemble that of a system at equilibrium. Here we apply this principle to a population of cells within a cytomorphological state space and demonstrate that cellular transition dynamics within this space can be described using equilibrium formalisms. We use this framework to map the effective energy landscape underlying the cytomorphological state space of a population of mouse embryonic fibroblasts (MEFs) and identify topographical nonuniformity in this space, indicating nonuniform occupation of cytomorphological states within an isogenic population. The introduction of exogenous apoptotic agents fundamentally altered this energy landscape, inducing formation of additional energy minima that correlated directly with changes in sensitivity to apoptosis induction. An equilibrium framework allows us to describe the behavior of an ensemble of single cells, suggesting that although cells are complex nonequilibrium systems, the application of formalisms derived from equilibrium thermodynamics can provide insight into the basis of nongenetic heterogeneities within cell populations, as well as the relationship between cytomorphological and functional heterogeneity