Interoperability of time series cytometric data: a cross platform approach for modeling tumor heterogeneity

Loon: Using Exemplars to Visualize Large-Scale Microscopy Data

Which drug is most promising for a cancer patient? A new microscopy-based approach for measuring the mass of individual cancer cells treated with different drugs promises to answer this question in only a few hours. However, the analysis pipeline for extracting data from these images is still far from complete automation: human intervention is necessary for quality control for preprocessing steps such as segmentation, adjusting filters, removing noise, and analyzing the result. To address this workflow, we developed Loon, a visualization tool for analyzing drug screening data based on quantitative phase microscopy imaging. Loon visualizes both derived data such as growth rates and imaging data. Since the images are collected automatically at a large scale, manual inspection of images and segmentations is infeasible. However, reviewing representative samples of cells is essential, both for quality control and for data analysis. We introduce a new approach for choosing and visualizing representative exemplar cells that retain a close connection to the low-level data. By tightly integrating the derived data visualization capabilities with the novel exemplar visualization and providing selection and filtering capabilities, Loon is well suited for making decisions about which drugs are suitable for a specific patient.

Biologically Relevant Heterogeneity: Metrics and Practical Insights

Heterogeneity is a fundamental property of biological systems at all scales that must be addressed in a wide range of biomedical applications, including basic biomedical research, drug discovery, diagnostics, and the implementation of precision medicine. There are a number of published approaches to characterizing heterogeneity in cells in vitro and in tissue sections. However, there are no generally accepted approaches for the detection and quantitation of heterogeneity that can be applied in a relatively high-throughput workflow. This review and perspective emphasizes the experimental methods that capture multiplexed cell-level data, as well as the need for standard metrics of the spatial, temporal, and population components of heterogeneity. A recommendation is made for the adoption of a set of three heterogeneity indices that can be implemented in any high-throughput workflow to optimize the decision-making process. In addition, a pairwise mutual information method is suggested as an approach to characterizing the spatial features of heterogeneity, especially in tissue-based imaging. Furthermore, metrics for temporal heterogeneity are in the early stages of development. Example studies indicate that the analysis of functional phenotypic heterogeneity can be exploited to guide decisions in the interpretation of biomedical experiments, drug discovery, diagnostics, and the design of optimal therapeutic strategies for individual patients.

FogBank: a single cell segmentation across multiple cell lines and image modalities

Background

Many cell lines currently used in medical research, such as cancer cells or stem cells, grow in confluent sheets or colonies. The biology of individual cells provide valuable information, thus the separation of touching cells in these microscopy images is critical for counting, identification and measurement of individual cells. Over-segmentation of single cells continues to be a major problem for methods based on morphological watershed due to the high level of noise in microscopy cell images. There is a need for a new segmentation method that is robust over a wide variety of biological images and can accurately separate individual cells even in challenging datasets such as confluent sheets or colonies.

Results

We present a new automated segmentation method called FogBank that accurately separates cells when confluent and touching each other. This technique is successfully applied to phase contrast, bright field, fluorescence microscopy and binary images. The method is based on morphological watershed principles with two new features to improve accuracy and minimize over-segmentation.

First, FogBank uses histogram binning to quantize pixel intensities which minimizes the image noise that causes over-segmentation. Second, FogBank uses a geodesic distance mask derived from raw images to detect the shapes of individual cells, in contrast to the more linear cell edges that other watershed-like algorithms produce.

We evaluated the segmentation accuracy against manually segmented datasets using two metrics. FogBank achieved segmentation accuracy on the order of 0.75 (1 being a perfect match). We compared our method with other available segmentation techniques in term of achieved performance over the reference data sets. FogBank outperformed all related algorithms. The accuracy has also been visually verified on data sets with 14 cell lines across 3 imaging modalities leading to 876 segmentation evaluation images.

Conclusions

FogBank produces single cell segmentation from confluent cell sheets with high accuracy. It can be applied to microscopy images of multiple cell lines and a variety of imaging modalities. The code for the segmentation method is available as open-source and includes a Graphical User Interface for user friendly execution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0431-x) contains supplementary material, which is available to authorized users.

Dynamic rendering of the heterogeneous cell response to anticancer treatments

The antiproliferative response to anticancer treatment is the result of concurrent responses in all cell cycle phases, extending over several cell generations, whose complexity is not captured by current methods. In the proposed experimental/computational approach, the contemporary use of time-lapse live cell microscopy and flow cytometric data supported the computer rendering of the proliferative process through the cell cycle and subsequent generations during/after treatment. The effects of treatments were modelled with modules describing the functional activity of the main pathways causing arrest, repair and cell death in each phase. A framework modelling environment was created, enabling us to apply different types of modules in each phase and test models at the complexity level justified by the available data. We challenged the method with time-course measures taken in parallel with flow cytometry and time-lapse live cell microscopy in X-ray-treated human ovarian cancer cells, spanning a wide range of doses. The most suitable model of the treatment, including the dose-response of each effect, was progressively built, combining modules with a rational strategy and fitting simultaneously all data of different doses and platforms. The final model gave for the first time the complete rendering in silico of the cycling process following X-ray exposure, providing separate and quantitative measures of the dose-dependence of G₁, S and G₂M checkpoint activities in subsequent generations, reconciling known effects of ionizing radiations and new insights in a unique scenario.

The antiproliferative response to anticancer treatment is the result of concurrent effects in all cell cycle phases, where molecular control pathways (checkpoints) are activated and cells may be arrested to repair DNA damage or killed if not able to succeed in the repair process. The complexity and inter-cell variability of these phenomena are not captured by the available methods, and the origin of the dose-dependence of the response remains elusive. In this work, we present an experimental-computational method that discloses and measures the individual responses of cell cycle controls in each phase and generation. We demonstrate that the method, exploiting jointly data sets obtained by flow cytometry and time-lapse in vivo imaging with a suitable experimental design, is able to achieve a full reconstruction in silico of the actual movement of cell cohorts following X-ray exposure, providing separate and quantitative measures of the dose-dependence of G₁, S and G₂M checkpoint activities in subsequent generations. Best fit parameters values are actual measures of the probability of activation of the specific pathways of arrest, repair or death within the cell population, linking the molecular scale to the “macroscopic” response, with full appreciation of its dynamics and inter-cell heterogeneity.