A novel view on stem cell development: analysing the shape of cellular genealogies

CellTrackVis: analyzing the performance of cell tracking algorithms

Live-cell imaging is a common data acquisition technique used by biologists to analyze cell behavior. Since manually tracking cells in a video sequence is extremely time-consuming, many automatic algorithms have been developed in the last twenty years to accomplish the task. However, none of these algorithms can yet claim robust tracking performance at the varying of acquisition conditions (e.g., cell type, acquisition device, cell treatments). While many visualization tools exist to help with cell behavior analysis, there are no tools to help with the algorithm's validation. This paper proposes CellTrackVis, a new visualization tool for evaluating cell tracking algorithms. CellTrackVis allows comparing automatically generated cell tracks with ground truth data to help biologists select the best-suited algorithm for their experimented pipeline. Moreover, CellTackVis can be used as a debugging tool while developing a new cell tracking algorithm to investigate where, when, and why each tracking error occurred.

Ontology patterns for the representation of quality changes of cells in time

Background

Cell tracking experiments, based on time-lapse microscopy, have become an important tool in biomedical research. The goal is the reconstruction of cell migration patterns, shape and state changes, and, comprehensive genealogical information from these data. This information can be used to develop process models of cellular dynamics. However, so far there has been no structured, standardized way of annotating and storing the tracking results, which is critical for comparative analysis and data integration. The key requirement to be satisfied by an ontology is the representation of a cell’s change over time. Unfortunately, popular ontology languages, such as Web Ontology Language (OWL), have limitations for the representation of temporal information. The current paper addresses the fundamental problem of modeling changes of qualities over time in biomedical ontologies specified in OWL.

Results

The presented analysis is a result of the lessons learned during the development of an ontology, intended for the annotation of cell tracking experiments. We present, discuss and evaluate various representation patterns for specifying cell changes in time. In particular, we discuss two patterns of temporally changing information: n-ary relation reification and 4d fluents. These representation schemes are formalized within the ontology language OWL and are aimed at the support for annotation of cell tracking experiments. We analyze the performance of each pattern with respect to standard criteria used in software engineering and data modeling, i.e. simplicity, scalability, extensibility and adequacy. We further discuss benefits, drawbacks, and the underlying design choices of each approach.

Conclusions

We demonstrate that patterns perform differently depending on the temporal distribution of modeled information. The optimal model can be constructed by combining two competitive approaches. Thus, we demonstrate that both reification and 4d fluents patterns can work hand in hand in a single ontology. Additionally, we have found that 4d fluents can be reconstructed by two patterns well known in the computer science community, i.e. state modeling and actor-role pattern.

Statistical and Mathematical Modeling of Spatiotemporal Dynamics of Stem Cells

Statistical and mathematical modeling are crucial to describe, interpret, compare, and predict the behavior of complex biological systems including the organization of hematopoietic stem and progenitor cells in the bone marrow environment. The current prominence of high-resolution and live-cell imaging data provides an unprecedented opportunity to study the spatiotemporal dynamics of these cells within their stem cell niche and learn more about aberrant, but also unperturbed, normal hematopoiesis. However, this requires careful quantitative statistical analysis of the spatial and temporal behavior of cells and the interaction with their microenvironment. Moreover, such quantification is a prerequisite for the construction of hypothesis-driven mathematical models that can provide mechanistic explanations by generating spatiotemporal dynamics that can be directly compared to experimental observations. Here, we provide a brief overview of statistical methods in analyzing spatial distribution of cells, cell motility, cell shapes, and cellular genealogies. We also describe cell-based modeling formalisms that allow researchers to simulate emergent behavior in a multicellular system based on a set of hypothesized mechanisms. Together, these methods provide a quantitative workflow for the analytic and synthetic study of the spatiotemporal behavior of hematopoietic stem and progenitor cells.

Estimation methods for heterogeneous cell population models in systems biology

Heterogeneity among individual cells is a characteristic and relevant feature of living systems. A range of experimental techniques to investigate this heterogeneity is available, and multiple modelling frameworks have been developed to describe and simulate the dynamics of heterogeneous populations. Measurement data are used to adjust computational models, which results in parameter and state estimation problems. Methods to solve these estimation problems need to take the specific properties of data and models into account. The aim of this review is to give an overview on the state of the art in estimation methods for heterogeneous cell population data and models. The focus is on models based on the population balance equation, but stochastic and individual-based models are also discussed. It starts with a brief discussion of common experimental approaches and types of measurement data that can be obtained in this context. The second part describes computational modelling frameworks for heterogeneous populations and the types of estimation problems occurring for these models. The third part starts with a discussion of observability and identifiability properties, after which the computational methods to solve the various estimation problems are described.

Lineage marker synchrony in hematopoietic genealogies refutes the PU.1/GATA1 toggle switch paradigm

Molecular regulation of cell fate decisions underlies health and disease. To identify molecules that are active or regulated during a decision, and not before or after, the decision time point is crucial. However, cell fate markers are usually delayed and the time of decision therefore unknown. Fortunately, dividing cells induce temporal correlations in their progeny, which allow for retrospective inference of the decision time point. We present a computational method to infer decision time points from correlated marker signals in genealogies and apply it to differentiating hematopoietic stem cells. We find that myeloid lineage decisions happen generations before lineage marker onsets. Inferred decision time points are in agreement with data from colony assay experiments. The levels of the myeloid transcription factor PU.1 do not change during, but long after the predicted lineage decision event, indicating  that the PU.1/GATA1 toggle switch paradigm cannot explain the initiation of early myeloid lineage choice.

Dissecting mechanisms of mouse embryonic stem cells heterogeneity through a model-based analysis of transcription factor dynamics

Pluripotent mouse embryonic stem cells (mESCs) show heterogeneous expression levels of transcription factors (TFs) involved in pluripotency regulation, among them Nanog and Rex1. The expression of both TFs can change dynamically between states of high and low activity, correlating with the cells' capacity for self-renewal. Stochastic fluctuations as well as sustained oscillations in gene expression are possible mechanisms to explain this behaviour, but the lack of suitable data hampered their clear distinction. Here, we present a systems biology approach in which novel experimental data on TF heterogeneity is complemented by an agent-based model of mESC self-renewal. Because the model accounts for intracellular interactions, cell divisions and heredity structures, it allows for evaluating the consistency of the proposed mechanisms with data on population growth and on TF dynamics after cell sorting. Our model-based analysis revealed that a bistable, noise-driven network model fulfils the minimal requirements to consistently explain Nanog and Rex1 expression dynamics in heterogeneous and sorted mESC populations. Moreover, we studied the impact of TF-related proliferation capacities on the frequency of state transitions and demonstrate that cellular genealogies can provide insights into the heredity structures of mESCs.

Elucidating functional heterogeneity in hematopoietic progenitor cells: a combined experimental and modeling approach

A detailed understanding of the mechanisms maintaining the hierarchical balance of cell types in hematopoiesis will be important for the therapeutic manipulation of normal and leukemic cells. Mathematical modeling is expected to make an important contribution to this area, but the iterative development of increasingly accurate models will rely on repeated validation using experimental data of sufficient resolution to distinguish between alternative model scenarios. The multipotent hematopoietic progenitor FDCP-Mix cells maintain a hierarchy from self-renewal to post-mitotic differentiation in vitro and are accessible to detailed analysis. Here, we report the development of a combined mathematical modeling and experimental approach to study the principles underlying heterogeneity in FDCP-Mix cultures. We adapt a single-cell based model of hematopoiesis to the conditions of cell culture and describe an association between proliferative history and phenotype of FDCP-Mix cells. While data derived from population studies are incapable of distinguishing between three mechanistically different model scenarios, statistical analysis of single cell tracking data provides a resolution sufficient to select one of them. This scenario favors differences between granulocytic and monocytic lineage with respect to their proliferative behavior and death rates as a mechanistic explanation for the observed heterogeneity. Our results demonstrate the power of a combined experimental/modeling approach in which single cell fate analysis is the key to revealing regulatory principles at the cellular level.

Tumor growth dynamics: insights into evolutionary processes

Identifying the types of event that drive tumor evolution and progression is crucial for understanding cancer. We suggest that the analysis of tumor growth dynamics can provide a window into tumor biology and evolution by connecting them with the types of genetic change that have occurred. Although fundamentally important, the documentation of tumor growth kinetics is more sparse in the literature than is the molecular analysis of cells. Here, we provide a historical summary of tumor growth patterns and argue that they can be classified into five basic categories. We then illustrate how those categories can provide insights into events that drive tumor progression, by discussing a particular evolutionary model as an example and encouraging such analysis in a more general setting.

Quantitative semi-automated analysis of morphogenesis with single-cell resolution in complex embryos

A quantitative understanding of tissue morphogenesis requires description of the movements of individual cells in space and over time. In transparent embryos, such as C. elegans, fluorescently labeled nuclei can be imaged in three-dimensional time-lapse (4D) movies and automatically tracked through early cleavage divisions up to ~350 nuclei. A similar analysis of later stages of C. elegans development has been challenging owing to the increased error rates of automated tracking of large numbers of densely packed nuclei. We present Nucleitracker4D, a freely available software solution for tracking nuclei in complex embryos that integrates automated tracking of nuclei in local searches with manual curation. Using these methods, we have been able to track >99% of all nuclei generated in the C. elegans embryo. Our analysis reveals that ventral enclosure of the epidermis is accompanied by complex coordinated migration of the neuronal substrate. We can efficiently track large numbers of migrating nuclei in 4D movies of zebrafish cardiac morphogenesis, suggesting that this approach is generally useful in situations in which the number, packing or dynamics of nuclei present challenges for automated tracking.

Time-lapse microscopy and classification of 2D human mesenchymal stem cells based on cell shape picks up myogenic from osteogenic and adipogenic differentiation

Current methods to characterize mesenchymal stem cells (MSCs) are limited to CD marker expression, plastic adherence and their ability to differentiate into adipogenic, osteogenic and chondrogenic precursors. It seems evident that stem cells undergoing differentiation should differ in many aspects, such as morphology and possibly also behaviour; however, such a correlation has not yet been exploited for fate prediction of MSCs. Primary human MSCs from bone marrow were expanded and pelleted to form high-density cultures and were then randomly divided into four groups to differentiate into adipogenic, osteogenic chondrogenic and myogenic progenitor cells. The cells were expanded as heterogeneous and tracked with time-lapse microscopy to record cell shape, using phase-contrast microscopy. The cells were segmented using a custom-made image-processing pipeline. Seven morphological features were extracted for each of the segmented cells. Statistical analysis was performed on the seven-dimensional feature vectors, using a tree-like classification method. Differentiation of cells was monitored with key marker genes and histology. Cells in differentiation media were expressing the key genes for each of the three pathways after 21 days, i.e. adipogenic, osteogenic and chondrogenic, which was also confirmed by histological staining. Time-lapse microscopy data were obtained and contained new evidence that two cell shape features, eccentricity and filopodia (= 'fingers') are highly informative to classify myogenic differentiation from all others. However, no robust classifiers could be identified for the other cell differentiation paths. The results suggest that non-invasive automated time-lapse microscopy could potentially be used to predict the stem cell fate of hMSCs for clinical application, based on morphology for earlier time-points. The classification is challenged by cell density, proliferation and possible unknown donor-specific factors, which affect the performance of morphology-based approaches.

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.

Novel regulators of stem cell fates identified by a multivariate phenotype screen of small compounds on human embryonic stem cell colonies

Understanding the complex mechanisms that govern the fate decisions of human embryonic stem cells (hESCs) is fundamental to their use in cell replacement therapies. The progress of dissecting these mechanisms will be facilitated by the availability of robust high-throughput screening assays on hESCs. In this study, we report an image-based high-content assay for detecting compounds that affect hESC survival or pluripotency. Our assay was designed to detect changes in the phenotype of hESC colonies by quantifying multiple parameters, including the number of cells in a colony, colony area and shape, intensity of nuclear staining, and the percentage of cells in the colony that express a marker of pluripotency (TRA-1-60), as well as the number of colonies per well. We used this assay to screen 1040 compounds from two commercial compound libraries, and identified 17 that promoted differentiation, as well as 5 that promoted survival of hESCs. Among the novel small compounds we identified with activity on hESC are several steroids that promote hESC differentiation and the antihypertensive drug, pinacidil, which affects hESC survival. The analysis of overlapping targets of pinacidil and the other survival compounds revealed that activity of PRK2, ROCK, MNK1, RSK1, and MSK1 kinases may contribute to the survival of hESCs.

Expertomica Cells: analysis of cell monolayer development

Expertomica Cells is a program for the creation and analysis of pedigree plots from time-lapse micrographs of cell monolayers. It enables recording the basic events in a cell cycle, cell neighbourhoods and spatial migration. The output is both numeric and graphical. The software helps to lower main hurdles in the manual analysis of cell monolayer development to practical limits; it reduces the operator processing time of typical experiment containing 5000 consecutive images from the usual 3 months to 3-10 h.

AVAILABILITY AND IMPLEMENTATION

Freely available on the web at http://www.expertomicacells.tk or http://www.expertomicacells.wu.cz. The source code is implemented in JAVA 6 and supported by Linux, Mac and MS Windows.

SUPPLEMENTARY INFORMATION

Supplementary data available at Bioinformatics online.

Designing materials to direct stem-cell fate

Proper tissue function and regeneration rely on robust spatial and temporal control of biophysical and biochemical microenvironmental cues through mechanisms that remain poorly understood. Biomaterials are rapidly being developed to display and deliver stem-cell-regulatory signals in a precise and near-physiological fashion, and serve as powerful artificial microenvironments in which to study and instruct stem-cell fate both in culture and in vivo. Further synergism of cell biological and biomaterials technologies promises to have a profound impact on stem-cell biology and provide insights that will advance stem-cell-based clinical approaches to tissue regeneration.