CellProfiler Analyst: data exploration and analysis software for complex image-based screens
BMC Bioinformatics volume 9, Article number: 482 (2008)
Image-based screens can produce hundreds of measured features for each of hundreds of millions of individual cells in a single experiment.
Here, we describe CellProfiler Analyst, open-source software for the interactive exploration and analysis of multidimensional data, particularly data from high-throughput, image-based experiments.
The system enables interactive data exploration for image-based screens and automated scoring of complex phenotypes that require combinations of multiple measured features per cell.
Visual analysis of cell samples has played a dominant role in the history of biology. The scientific community has only begun to scratch the surface of computationally extracting the rich information visible in fluorescence microscopy images of cell samples . This capability is increasingly important given the ease now to systematically perturb cells with libraries of chemicals or gene-perturbing reagents like RNA interference or gene overexpression and collect hundreds of thousands of images of these cell samples [2, 3]. We recently developed open-source image analysis software, CellProfiler, which measures a rich set of cellular features in images, such as size, shape, and staining patterns including intensity, texture, and colocalization [4, 5]http://www.cellprofiler.org. This tool has been useful for extracting image-based measurements to score sophisticated screens [6–8], with many more in progress.
The volume and richness of individual-cell data from large image-based screens is unprecedented and existing software is inadequate for the challenge of data analysis. For analysis of small or very simple experiments, spreadsheet programs like Microsoft Excel are sufficient, and useful open-source tools exist for analysis and exploration of data from high throughput screens in general [9–12]. Existing software packages targeted for image-based screening, however, have one or more limitations which prevent sophisticated visualization and extraction of information from image-based screens: (a) they are not designed for the hierarchical data structure inherent in image-based data (each treatment condition is replicated in several samples, each sample is usually represented by several images, each image contains a population of cells, and each cell has hundreds of associated measures), (b) they ignore the inherent biological variability of cell populations such that assays requiring subpopulation analysis cannot be scored, (c) they cannot handle the volumes of data typical in image-based experiments (e.g., ~500 measurements for each of ~100 million individual cells), (d) they provide limited linking to raw or processed image data or chemical structure data, (e) they allow only limited statistical analyses of the data, (f) they are proprietary and new methods cannot be easily added, (g) they are limited to data from a particular image analysis package, (h) they require expertise in statistics or programming, and/or (i) they require intense hands-on data management.
Given that no existing tools meet the specific needs of image-based screens, researchers have needed computational expertise to directly query databases of image-based information using command-line tools. Often, the researchers best able to explore and interpret the data lack these computational skills. These researchers are therefore less likely to make serendipitous discoveries (or identify quality-control issues) in their image-based screens, which inherently contain enormous amounts of information beyond that which is pertinent to the original, intended biological question. It is critical to provide exploration tools to screening researchers, tools that employ their understanding of the experiment in question and their creativity and ability to recognize and interpret patterns and relationships within data. These capabilities flourish when united with a computer's unique ability to store, retrieve, display, and quantitatively analyze billions of data points.
We therefore sought to develop a software system that would make high-dimensional image-based data exploration feasible for researchers who lack computational skills, and flexible for computer scientists who want to develop and add advanced new methods for image-based screening, such as machine learning-based phenotype scoring. We describe here the result of our work, an open-source software package called "CellProfiler Analyst".
Results and discussion
Four types of plots are the starting points for exploration of large, multi-dimensional image-based screens in CellProfiler Analyst (Figure 1). Importantly, these tools are compatible with the scale of data typically acquired in image-based screens, which can be hundreds of features for each of hundreds of millions of cells. Histograms display the distribution of values for one measured feature by grouping image or object data into evenly spaced bins, on a linear or logarithmic scale (Figure 1a). Such plots can be helpful, for example, to examine the cell cycle status of samples (by plotting per-cell DNA content) or to examine outliers for quality control purposes (e.g., by plotting per-image cell counts). Two measured features per image or object can be displayed on the same chart via a scatterplot (Figure 1b), which is also useful for identifying hits and for quality control purposes. For example, the researcher can readily exclude out-of-focus images from analysis based on measurements made by CellProfiler's "Measure Image Quality" module. Because data points in scatterplots can occlude each other, they are typically unsuitable for individual cell data where hundreds of millions of data points are examined to identify interesting subpopulations. For these cases, a density plot is more appropriate (Figure 1c). Every pixel in the plot represents a histogram "bin" and the color of the pixel represents the number of data points in the bin. These plots are useful, for example, for establishing thresholds at which to classify individual cells as "positive" or "negative" based on two features (e.g., based on two intensity measures as in flow cytometry). To explore more than two measured features of each image or data point, a parallel coordinate plot is used. Parallel coordinate plots  allow analysis of multiple dimensions of data, whereby each measured feature's scaled (0–1) values are given a separate y-axis and individual data points are connected across these multiple axes (Figure 1d).
Each data point in a plot can represent an individual cell or, by contrast, the mean value of the population of cells within an image. Data can also be grouped by characteristics the samples have in common (e.g., chemical name or dose). Multiple experiments that investigate the same set of treatment conditions (e.g., chemical compounds or RNA interference reagents) can be grouped together, which eases analysis of replicates. For all types of plots, the data to be displayed can be filtered, for example to plot data only from a single image, from a sample of data points at specified equal intervals, or data that satisfies certain criteria (specified in SQL "where" clauses like "CellCount > 100").
Exploring relationships among data
Data points selected and highlighted in one plot are immediately highlighted in all other open plots (a technique often called "brushing" ) such that a sample or set of samples can be examined in the context of other sets of samples (Figure 2). This allows, for example, the comparison of measurements from samples of interest vs. all samples in the experiment. Brushing helps the user to more easily examine relationships in the data, especially when the data has a large number of attributes or items, when the data spans multiple experiments (including, for example, replicates), or when it is natural to examine different parts of the data using different views. The brushing concept is extended in CellProfiler Analyst for situations where multiple experiments are being simultaneously explored: when a point corresponding to a particular image is highlighted, all points corresponding to that experimental treatment condition can be highlighted, even if the data comes from multiple experiments that are being examined together. In the scatterplot in Figure 1b, for example, four data points are blue because one was originally selected and the user requested that replicates for that sample be highlighted.
Interesting data points or sets of data points can be investigated by drilling down into the data in several ways (Figure 3). For plots showing data points representing image measurements, a data point or set of data points can be selected and the original images that produced the data point can be displayed (Figure 3d). This can reveal artifacts in sample preparation or imaging, such as fluorescent test compounds, aggregates or overabundance of staining reagents, fibers, or debris (Figure 3g). These artifacts not only occlude actual cells in images but can also disrupt the proper identification and measurement of remaining cells in the image. For these and other reasons, images showing identification outlines resulting from image analysis (if available) can also be shown for selected data points (Figure 3e), to identify whether the identification of cells occurred properly. This is an important consideration given that no segmentation algorithms are flawless.
Additionally, a data point or set of data points can be selected and a plot of the measurements of individual cells that were present in those images can be displayed as a separate subplot. This allows, for example, a DNA content histogram indicating cell cycle distribution of the cell population to be displayed for a particular image or set of images of interest (Figure 2c and Figure 3b). To investigate the identity of interesting samples, a simple list of the treatment conditions that produced a set of data points can be displayed to get an overview (Figure 2d). For further investigation, web-based information about each image's treatment condition can be launched in an external web browser (Figure 3f), if web addresses associated with each sample are stored in the database. All available measurements and other information for a particular sample can be displayed in a simple table and saved as a comma-delimited text file for analysis in another software package (Figure 3c).
Gating individual cell data to score complex phenotypes
Image-based data is tremendously valuable in that multiple single-cell measurements are available. Responses of individual cells to a treatment are usually inhomogeneous because of cell cycle variations or differences in protein levels due to memory or stochastic noise [15, 16]. In many cases, a single measured feature (e.g., the total intensity of red stain within the nucleus) can be used to score individual cells and the only challenge is to identify a suitable threshold for scoring positive cells. This can be accomplished in CellProfiler Analyst using histograms of individual cell data. For complex phenotypes, several features of each cell may be required for effective scoring. In these cases, a density plot showing individual cells (Figure 4a) can be useful for identifying interesting cell subpopulations, by delineating a section of the plot (often called "gating"). Whether the gate contains the cells of interest can be tested using two features: the "Show Object Montage" feature to see what individual cells within the gate look like (Figure 4b), and the "Show Image" feature to see whether cells within a particular sample are appropriately marked as inside or outside the gate (Figure 4c). Once the final, desired subpopulation of cells is gated, the number of cells that fall within that subpopulation is calculated for each image, for further statistical analysis (Figure 4d). As an example, when DNA and phosphorylated Serine 10 of histone H3 are both stained, a simple two-feature gate in CellProfiler Analyst enables scoring mitotic subphases in human HT29 cells (Figure 4e). Many software systems perform image analysis on the fly during image acquisition; in such cases, a threshold value for a feature of interest must be chosen in advance to score the screen. By contrast, these tools in CellProfiler Analyst allow testing the efficacy of scoring based on different features and different measurement thresholds.
If more than two features are needed to score a phenotype, sequential gates can be used upon the cell data. This approach is applied as follows: (1) display the entire population of cells from an experiment in a density plot, (2) draw a gate around the data points representing potential cells of interest, (3) adjust the gate to include nearly all positive cells and exclude as many negative cells as possible, (4) plot the resulting gated subpopulation in a new density plot with two new measurement features as axes, (5) gate the subpopulation again based on these new features, and (6) calculate the percentage of each image's cells that fall within the final gate.
Case study: mitotic subphase screen
We wanted to test CellProfiler Analyst's ability to plot, explore, and filter individual cell data to identify subpopulations defined by several morphological features. We chose to identify Drosophila melanogaster Kc167 cells in telophase and metaphase of the cell cycle, using only a DNA stain. Identification of samples with perturbed cell cycle regulation is of clear importance to normal cell biology as well as cancer studies. Regulators of the cell cycle have been sought intensively for decades via traditional and high-throughput screens for changes in overall cell cycle distribution or for increased phospho-histone H3 staining, a marker of cells in late G2 and M phase (e.g.,  and references therein). We reasoned that additional genes might exist which, when perturbed, yield increased numbers of metaphase- or telophase-stage nuclei without substantially affecting the overall mitotic index (phospho-histone H3 staining) or cell cycle distribution. While we were not aware of any positive controls with such a phenotype, we suspected such genes might have been previously overlooked because we noticed that not all metaphase nuclei stain brightly for phospho-histone H3 (Figure 5a), for unknown reasons. Identifying genes whose RNAi produces cells appearing to be in particular subphases of mitosis, regardless of concomitant phospho-histone H3 staining, would be a first step towards understanding these phenomena.
Several groups have tested automated methods for scoring mitotic subphases [18–20]; these studies were accomplished by computational tools tailored to the specific assay and often relied on multiple cellular stains. Machine learning methods have been explored by our own group and others [21–26] (and see Conclusions), but we also wanted to explore allowing the user to manually select a small number of features of known biological relevance, followed by sequential gating on those features. This would give the researcher full control over the features used in the scoring, and the scoring would be more readily transferable from one experiment to the next because a small number of features are selected. We therefore wanted to score mitotic subphases using a DNA stain only, using supervised selection of measurements followed by sequential gating on those measurements, in the context of a software package usable by a non-computer scientist.
Image scoring by sequential gating of individual cell data
We screened genes using Drosophila RNA interference living cell microarrays [27–29] to identify gene "knockdowns" that yield a disproportionate number of cells in two sub-phases of mitosis: metaphase and anaphase/telophase (referred to as telophase for simplicity). We created and analyzed 5 replicates of a Drosophila array, with 1120 spots of dsRNA on a single microscope slide (Figure 5b), including three replicate spots for each of 288 genes (mostly kinases and phosphatases), plus 256 negative control spots lacking dsRNA. Some phenotypes produced in these Drosophila Kc167 cells (e.g. cell death) are visible at low resolution (5× lens; Figure 5c), but to identify telophase and metaphase nuclei we collected individual high resolution images within each spot on each slide (40× lens; small portion of one image shown in Figure 5d).
We began with the telophase phenotype. To determine which measured cellular features would be most effective for scoring, we handpicked representative telophase nuclei and normal G2-phase nuclei from random screening images and created image montages for these two classes (Figure 6a) using Adobe Photoshop. We used CellProfiler to measure nuclear features in these montage images, then exported the results to Excel and selected five features to use for sequential gating, based on a combination of biological intuition plus the quantitative ability of each feature to discriminate telophase from normal nuclei, using simple statistical tests in Excel. The selected features included DNA content, intensity, shape, and texture features (additional data file 1).
We then interactively developed sequential gates using density plots of these features in CellProfiler Analyst (see "Gating individual cell data to score complex phenotypes" section). To accomplish this task, CellProfiler was used to process the full set of screening images and load the resulting data into a database (2.8 million cells × 396 features/cell = 1.1 billion measurements total). This allowed us to display all individual cells in the experiment in an initial density plot with two of our selected features as axes, i.e., DNA content and size (area) of the nucleus. We drew an initial gate around the 2N DNA content peak and small nuclear area, and empirically refined the gate for telophase cells by examining images of the gated nuclei and adjusting the gate's boundaries accordingly. While automated approaches could certainly identify a boundary based on a researcher-provided training set, this manual approach allows the biologist to specifically assess many cells near the relevant boundaries. Once the appropriate gate was selected for the initial density plot, the subpopulation was transferred to a new density plot with two new features used as axes and the next gate was created, again finding the optimal parameters to distinguish telophase nuclei from all other nuclei. This procedure was repeated for the fifth, and final, selected feature. Once the final gate was refined, we applied the sequential gates to a new set of images and confirmed that their scoring was effective (Table 1 and Figure 6b), successfully differentiating telophase from other nuclei. In creating the gates, we tried to minimize the false positive rate while accepting a higher false negative rate (Table 1). We reasoned that true hits would have enough positives to overcome this intentionally stringent selection procedure. At this point, we applied the final sequential gates to all the cells in order to score the entire screen for the telophase phenotype. We found that the gates must typically be adjusted slightly between different replicate slides due to inter-experiment variability (e.g., staining intensity), although experiment-to-experiment normalization methods could be explored to reduce this effect.
We separately performed the same procedure for the metaphase phenotype (using four features to distinguish metaphase nuclei from all other nuclei); a complete list of the 288 genes tested and their scores for telophase and metaphase is shown in additional data file 2.
Rank-ordering samples by the percentage of telophase nuclei revealed 4 gene knockdowns with a significant increase in telophase nuclei (Figure 6c, first 4 rows). Validating the approach, two of the genes are PP2A complex subunits that have been previously associated with mitosis: the PP2A-C catalytic subunit mts (CG7109/microtubule star) and a PP2A-A family regulatory subunit (CG17291/CG33297/CG13383, Note: dicistronic with CSN8). RNAi against both genes increased the percentage of cells that were phospho-histone H3-positive (Figure 6c, fifth column). A third hit, Ck1α (Casein kinase 1α/CG2028), has also previously been linked to mitosis (Figure 6c, last column). We noticed that its knockdown by RNAi produced nuclei whose chromatin appeared to be slightly less condensed than typical telophase nuclei (Figure 7), while still more condensed than interphase nuclei. The percentage of cells that were phospho-histone H3-positive was normal (Figure 6c, fifth column). Together, these observations suggest that this defect occurs in late-stage telophase/anaphase. The fourth hit was a predicted kinase with no functional annotation (CG8878). Visual inspection revealed that nearly all nuclei in these samples appeared brighter and more compact than controls, a subtle but reproducible effect (Figure 7). This understandably resulted in more of the 2N nuclei being counted as having telophase-like morphology. We found that these cells were not enriched for phospho-histone H3-positivity (Figure 6c, fifth column); without further experimentation, it is unclear whether this is a true late-stage mitotic phenotype or rather a condensed nuclei phenotype.
Interestingly, the only metaphase hit in this screen (Figure 6c, last row) is the B'/B56 subfamily regulatory subunit of PP2A (CG5643/widerborst), which at the time of our screen had not been linked to cell cycle regulation. The percentage of cells that were phospho-histone H3-positive was not much higher than normal (Figure 6c, fifth column). We confirmed by eye the metaphase-inducing phenotype of widerborst knockdown in the original images and in separate experiments with two other dsRNAs, including one that was non-overlapping with the original (Figure 8a). Widerborst is an essential gene involved in planar cell polarization  and apoptosis [31, 32]. Notably, in other contexts (circadian clock protein cycling  and sensory organ development ) widerborst is indirectly linked to the B/PR55 subfamily member twins/aar, which is itself known to be required for metaphase to anaphase transition . Our work therefore confirms, with non-overlapping dsRNAs, a recently reported cell cycle regulation role for widerborst  and together indicates that it is unlikely this phenotype is due to off-target effects [37, 38].
The closest human homolog of widerborst is PPP2R5E, the epsilon isoform of a subfamily of PP2R5 (a.k.a. B'/PR61/B56) regulatory subunits of the PP2A complex. As yet, no particular function has been associated with PPP2R5E. We wondered whether PPP2R5E might be a B' regulatory subunit that modulates the known role for PP2A in mitosis, given our finding of its homolog widerborst's role in Drosophila. PPP2R5E knockdown did not increase the mitotic index significantly in recent RNA interference screens for increased phospho-histone H3 [7, 39, 40]. However, when we scored these same PPP2R5E-knockdown images for metaphase morphology, rather than phospho-histone H3 levels, we discovered a metaphase-arrest phenotype for PPP2R5E knockdown, confirmed by two different shRNAs (Figure 8b), consistent with the phenotype seen for widerborst in Drosophila. Whether widerborst/PPP2R5E are themselves required for metaphase-to-anaphase transition or whether their depletion causes the phenotype by specifically disrupting the stoichiometry of the relevant PP2A complex remains to be determined. Recent findings that PPP2R5E localizes to centromeres and that the B' subfamily of regulatory subunits are required for proper meiotic sister chromatid separation in fission and budding yeast  support the idea that this family of subunits is indeed important for proper chromatin dynamics during cell division.
We have described here a software system for exploration and analysis of large, hierarchical, multi-dimensional data sets. While it is compatible with any type of data (e.g., players on teams, trees within forests), it is particularly capable of high-end exploration and analysis of measured features from high-throughput image-based screens for both quality control and identifying hits in a screen. Researchers are welcome to download the Java source code and add new types of plots and analysis tools (e.g., for normalizing screen data [9, 10]) to the system.
We have demonstrated the utility of this software for interactive data exploration and analysis – especially for intentionally selecting cells with particular measurement values in order to score complex visual phenotypes. Of course, often the features that successfully specify a particular phenotype are either unknown or so numerous as to make the sequential plotting shown here impractical, and choosing decision boundaries empirically may not be optimal to score the phenotype. For these reasons, we recently added machine-learning methodology to CellProfiler Analyst (TRJ, AEC, DMS, PG, unpublished data). Nonetheless, the complete control over features and thresholds offered by sequential gating is quite useful in some cases. Often a researcher needs to ignore certain features of positive control cells (for example, when a positive control treatment has pleiotropic effects on cells) and emphasize other, better-understood cellular features. Interactive observation of the original cellular images while making gating decisions to define a phenotype also leverages the biologist's intuition about a phenotype. Within the same open-source software infrastructure, both approaches (sequential gating and machine learning) can now be applied to large-scale imaging screens.
Availability and requirements
* Project name: CellProfiler Analyst
* Project home page: http://www.cellprofiler.org
* Operating systems: Platform independent (Mac, Windows, and Unix)
* Programming language: Java
* Other requirements: Java 1.4.2 or greater. For full functionality, CellProfiler Analyst requires Java 1.5.0_6 or greater, Python version 2.5 or greater http://www.python.org/ and the NumPy Python package http://scipy.org.
* License: GNU GENERAL PUBLIC LICENSE, Version 2
* No additional restrictions to use by non-academics
CellProfiler Analyst can be downloaded for Mac, Windows, and Unix operating systems from the CellProfiler Project website http://www.cellprofiler.org, where it is distributed under an open-source license (GNU General Public License, version 2). An archived version is also available as additional data file 3, submitted with this article. The Examples page of the website provides demonstration movies showing the software in use, an example database and images, and links to an online forum where questions about the software are answered.
CellProfiler Analyst is designed to explore and analyze any MySQL database of image-based screening data that follows a simple format: at least one image table with rows corresponding to images and columns of image data (examples of columns are: the name of the treatment condition, total intensity of the entire image, mean cell area averaged over all cells in the image, path to the original image), and at least one object table with rows corresponding to objects (e.g., cells) and columns of object data (examples of columns are: area of the cell, intensity of DNA stain in the nucleus, location of the cell in the original image – the latter being important for viewing individual cells during exploration). This data format is automatically produced if images are analyzed with CellProfiler open-source cell image analysis software http://www.cellprofiler.org, using its ExportToDatabase module. The data should be normalized for plate-to-plate or spatial-layout variations prior to exploration in CellProfiler Analyst. While the software is designed to access remote databases because typical data sets are far too large to be stored in physical memory, the "Make Local Object Table" option allows particularly relevant measurements to be stored locally in memory to speed analysis while still allowing access to the full dataset in the remote database.
We prepared Drosophila Kc167 cells as previously described . In brief, cells were grown on living cell microarrays with spots of double stranded RNA for 3 days. For confirmation of phenotypes in Drosophila, we grew cells on plain slides for 3 days, after being pre-treated with dsRNA for 2 days. We used images of human HT29 cells as previously described .
For the screen of the metaphase and telophase phenotypes, each gene was tested in three replicate spots on five independently prepared cell array slides, and the results for all genes are shown in Additional file 2. Because the three replicate spots were near each other, cell counts for the groups of three were accumulated and not treated as independent samples. A p-value for each gene on each of the five slides was calculated based on the number of metaphase nuclei found and the number of cells total, relative to the average percentage of metaphase nuclei on the entire slide (i.e., as a Bernoulli random variable). To add stringency, we report results for second- and third-strongest scoring replicates only (shown on two separate sheets of Additional file 2). We required that two or three of the five scores were above a threshold that results in a combined p-value below 0.01. For Bonferroni-adjusted p-values from single experiments, these thresholds are 0.6 for two experiments (out of 5), and 5.2 for three (out of 5). De-enriched samples are listed with a p-value of 1 and samples with a p-value of 1 are ordered by enrichment. "Enrichment" is the fold-enrichment of the sample relative to all the samples.
In the bar charts in the fourth column of Figure 6c, we statistically analyzed the cell cycle distribution and cell count for the screens' hits. To do this, we first gathered DNA content data (i.e., integrated nuclear DNA intensity) from the database for all cells on the slides where the hits occurred. Then, to normalize for illumination and staining variation between slides and between images, the DNA content measurements were log2-transformed and shifted so that the mode of the DNA content for each image (calculated by binning the log2-transformed DNA into 50 bins) was equal to 1. Based on this normalized log2(DNA intensity), cells were then counted as 2N, 4N, and 8N as follows:
[-0.5, 0.5) was categorized as "2N"
[0.5, 1.5) was categorized as " 4N"
[1.5, 2.5) was categorized as " 8N", although this includes ~6N to ~11N
Using the resulting cell counts for each subpopulation (2N, 4N, 8N), we calculated p-values as follows: first, subpopulation counts were converted to fractions for each image by dividing the subpopulation counts by the total number of cells in the image (taken as the sum of the 2N, 4N, and 8N subpopulations). Each fraction was then normalized by the median fraction for that subpopulation on that slide, to account for any per-slide biases in cell-cycle distribution. These normalized fractions were averaged across replicate samples for each gene. Lastly, these averaged normalized fractions were used to calculate p-values for each subpopulation in a 10,000-trial permutation test (where labels were permuted within slides, but not between slides, to ensure that the same number of images was taken from the slide as in the experiment). Cell-count p-values were calculated similarly: the total number of cells in each image was normalized by the median per-image cell count on that slide, to prevent biases for more densely populated slides, and a permutation test was performed on the average normalized cell-count. For cell cycle and cell count, p-values were Bonferroni-corrected for 20 experiments (5 genes examined for 4 populations: 2N, 4N, 8N, and count).
Wollman R, Stuurman N: High throughput microscopy: from raw images to discoveries. J Cell Sci 2007, 120(Pt 21):3715–3722. 10.1242/jcs.013623
Carpenter AE: Image-based chemical screening. Nat Chem Biol 2007, 3(8):461–465. 10.1038/nchembio.2007.15
Carpenter AE, Sabatini DM: Systematic genome-wide screens of gene function. Nat Rev Genet 2004, 5(1):11–22. 10.1038/nrg1248
Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang IH, Friman O, Guertin DA, Chang JH, Lindquist RA, Moffat J, et al.: CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol 2006, 7(10):R100. 10.1186/gb-2006-7-10-r100
Lamprecht MR, Sabatini DM, Carpenter AE: CellProfiler: free, versatile software for automated biological image analysis. Biotechniques 2007, 42(1):71–75. 10.2144/000112257
Kuijl C, Savage ND, Marsman M, Tuin AW, Janssen L, Egan DA, Ketema M, Nieuwendijk R, Eeden SJ, Geluk A, et al.: Intracellular bacterial growth is controlled by a kinase network around PKB/AKT1. Nature 2007, 450(7170):725–730. 10.1038/nature06345
Moffat J, Grueneberg DA, Yang X, Kim SY, Kloepfer AM, Hinkle G, Piqani B, Eisenhaure TM, Luo B, Grenier JK, et al.: A Lentiviral RNAi Library for Human and Mouse Genes Applied to an Arrayed Viral High-Content Screen. Cell 2006, 124(6):1283–1298. 10.1016/j.cell.2006.01.040
Dobbelaere J, Josue F, Suijkerbuijk S, Baum B, Tapon N, Raff J: A Genome-Wide RNAi Screen to Dissect Centriole Duplication and Centrosome Maturation in Drosophila. PLoS Biol 2008, 6(9):e224. 10.1371/journal.pbio.0060224
Malo N, Hanley JA, Cerquozzi S, Pelletier J, Nadon R: Statistical practice in high-throughput screening data analysis. Nat Biotechnol 2006, 24(2):167–175. 10.1038/nbt1186
Boutros M, Bras LP, Huber W: Analysis of cell-based RNAi screens. Genome Biol 2006, 7(7):R66. 10.1186/gb-2006-7-7-r66
Cook D, Swayne DF: Interactive and Dynamic Graphics for Data Analysis: With Examples Using R and GGobi: Springer. 2007.
Gee AG, Li H, Yu M, Smrtic MB, Cvek U, Goodell H, Gupta V, Lawrence C, Zhou J, Chiang C, et al.: The Universal Visualization Platform. SPIE Visualization and Data Analysis 2005, 5669: 274–283.
Inselberg A: Multidimensional Detective. In Readings in Information Visualization. Edited by: Card SK, Mackinlay JD, Shneiderman B. Morgan Kaufman; 1999.
Becker RA, Cleveland WS: Brushing scatterplots. Technometrics 1987, 29(2):127–142. 10.2307/1269768
Levsky JM, Singer RH: Gene expression and the myth of the average cell. Trends Cell Biol 2003, 13(1):4–6. 10.1016/S0962-8924(02)00002-8
Sigal A, Milo R, Cohen A, Geva-Zatorsky N, Klein Y, Alaluf I, Swerdlin N, Perzov N, Danon T, Liron Y, et al.: Dynamic proteomics in individual human cells uncovers widespread cell-cycle dependence of nuclear proteins. Nat Methods 2006, 3(7):525–531. 10.1038/nmeth892
Bjorklund M, Taipale M, Varjosalo M, Saharinen J, Lahdenpera J, Taipale J: Identification of pathways regulating cell size and cell-cycle progression by RNAi. Nature 2006, 439(7079):1009–1013. 10.1038/nature04469
Gasparri F, Cappella P, Galvani A: Multiparametric Cell Cycle Analysis by Automated Microscopy. J Biomol Screen 2006.
Wilson CJ, Si Y, Thompson CM, Smellie A, Ashwell MA, Liu JF, Ye P, Yohannes D, Ng SC: Identification of a Small Molecule That Induces Mitotic Arrest Using a Simplified High-Content Screening Assay and Data Analysis Method. J Biomol Screen 2006, 11(1):21–28. 10.1177/1087057105280726
Goshima G, Wollman R, Goodwin SS, Zhang N, Scholey JM, Vale RD, Stuurman N: Genes required for mitotic spindle assembly in Drosophila S2 cells. Science 2007, 316(5823):417–421. 10.1126/science.1141314
Neumann B, Held M, Liebel U, Erfle H, Rogers P, Pepperkok R, Ellenberg J: High-throughput RNAi screening by time-lapse imaging of live human cells. Nat Methods 2006, 3(5):385–390. 10.1038/nmeth876
Harder N, Neumann B, Held M, Liebel U, Erfle H, Ellenberg J, Eils R, Rohr K: Automated Recognition of Mitotic Patterns in Fluorescence Microscopy Images of Human Cells. IEEE ISBI 2006.
Wang M, Zhou X, Li F, Huckins J, King RW, Wong ST: Novel Cell Segmentation and Online SVM for Cell Cycle Phase Identification in Automated Microscopy. Bioinformatics 2008, 24(1):94–101. 10.1093/bioinformatics/btm530
Tao CY, Hoyt J, Yan F: A support vector machine classifier for recognizing mitotic subphases using high-content screening data. J Biomol Screen 2007, 12(4):490–496. 10.1177/1087057107300707
Gambe AE, Ono RM, Matsunaga S, Kutsuna N, Higaki T, Higashi T, Hasezawa S, Uchiyama S, Fukui K: Development of a multistage classifier for a monitoring system of cell activity based on imaging of chromosomal dynamics. Cytometry A 2007, 71(5):286–296.
Mukherji M, Bell R, Supekova L, Wang Y, Orth AP, Batalov S, Miraglia L, Huesken D, Lange J, Martin C, et al.: Genome-wide functional analysis of human cell-cycle regulators. Proc Natl Acad Sci USA 2006, 103(40):14819–14824. 10.1073/pnas.0604320103
Ziauddin J, Sabatini DM: Microarrays of cells expressing defined cDNAs. Nature 2001, 411(6833):107–110. 10.1038/35075114
Wheeler DB, Carpenter AE, Sabatini DM: Cell microarrays and RNA interference chip away at gene function. Nat Genet 2005, 37(Suppl):S25–30. 10.1038/ng1560
Wheeler DB, Bailey SN, Guertin DA, Carpenter AE, Higgins CO, Sabatini DM: RNAi living-cell microarrays for loss-of-function screens in Drosophila melanogaster cells. Nat Methods 2004, 1(2):127–132. 10.1038/nmeth711
Hannus M, Feiguin F, Heisenberg CP, Eaton S: Planar cell polarization requires Widerborst, a B' regulatory subunit of protein phosphatase 2A. Development 2002, 129(14):3493–3503.
Li X, Scuderi A, Letsou A, Virshup DM: B56-associated protein phosphatase 2A is required for survival and protects from apoptosis in Drosophila melanogaster. Mol Cell Biol 2002, 22(11):3674–3684. 10.1128/MCB.22.11.3674-3684.2002
Liu W, Silverstein AM, Shu H, Martinez B, Mumby MC: A functional genomics analysis of the B56 isoforms of Drosophila protein phosphatase 2A. Mol Cell Proteomics 2007, 6(2):319–332.
Sathyanarayanan S, Zheng X, Xiao R, Sehgal A: Posttranslational regulation of Drosophila PERIOD protein by protein phosphatase 2A. Cell 2004, 116(4):603–615. 10.1016/S0092-8674(04)00128-X
Abdelilah-Seyfried S, Chan YM, Zeng C, Justice NJ, Younger-Shepherd S, Sharp LE, Barbel S, Meadows SA, Jan LY, Jan YN: A gain-of-function screen for genes that affect the development of the Drosophila adult external sensory organ. Genetics 2000, 155(2):733–752.
Mayer-Jaekel RE, Ohkura H, Gomes R, Sunkel CE, Baumgartner S, Hemmings BA, Glover DM: The 55 kd regulatory subunit of Drosophila protein phosphatase 2A is required for anaphase. Cell 1993, 72(4):621–633. 10.1016/0092-8674(93)90080-A
Chen F, Archambault V, Kar A, Lio P, D'Avino PP, Sinka R, Lilley K, Laue ED, Deak P, Capalbo L, et al.: Multiple protein phosphatases are required for mitosis in Drosophila. Curr Biol 2007, 17(4):293–303. 10.1016/j.cub.2007.01.068
Ma Y, Creanga A, Lum L, Beachy PA: Prevalence of off-target effects in Drosophila RNA interference screens. Nature 2006, 443(7109):359–363. 10.1038/nature05179
Kulkarni MM, Booker M, Silver SJ, Friedman A, Hong P, Perrimon N, Mathey-Prevot B: Evidence of off-target effects associated with long dsRNAs in Drosophila melanogaster cell-based assays. Nat Methods 2006, 3(10):833–838.
Rines DR, Gomez-Ferreria MA, Zhou Y, DeJesus P, Grob S, Batalov S, Labow M, Huesken D, Mickanin C, Hall J, et al.: Whole genome functional analysis identifies novel components required for mitotic spindle integrity in human cells. Genome Biol 2008, 9: R44. 10.1186/gb-2008-9-2-r44
Kittler R, Pelletier L, Heninger AK, Slabicki M, Theis M, Miroslaw L, Poser I, Lawo S, Grabner H, Kozak K, et al.: Genome-scale RNAi profiling of cell division in human tissue culture cells. Nat Cell Biol 2007, 9(12):1401–1412. 10.1038/ncb1659
Riedel CG, Katis VL, Katou Y, Mori S, Itoh T, Helmhart W, Galova M, Petronczki M, Gregan J, Cetin B, et al.: Protein phosphatase 2A protects centromeric sister chromatid cohesion during meiosis I. Nature 2006.
Bettencourt-Dias M, Giet R, Sinka R, Mazumdar A, Lock WG, Balloux F, Zafiropoulos PJ, Yamaguchi S, Winter S, Carthew RW, et al.: Genome-wide survey of protein kinases required for cell cycle progression. Nature 2004, 432(7020):980–987. 10.1038/nature03160
Eggert US, Kiger AA, Richter C, Perlman ZE, Perrimon N, Mitchison TJ, Field CM: Parallel chemical genetic and genome-wide RNAi screens identify cytokinesis inhibitors and targets. PLoS Biol 2004, 2(12):e379. 10.1371/journal.pbio.0020379
Snaith HA, Armstrong CG, Guo Y, Kaiser K, Cohen PT: Deficiency of protein phosphatase 2A uncouples the nuclear and centrosome cycles and prevents attachment of microtubules to the kinetochore in Drosophila microtubule star (mts) embryos. J Cell Sci 1996, 109(Pt 13):3001–3012.
Kitajima TS, Sakuno T, Ishiguro KI, Iemura SI, Natsume T, Kawashima SA, Watanabe Y: Shugoshin collaborates with protein phosphatase 2A to protect cohesin. Nature 2006.
The authors are grateful to Michael R. Lamprecht, David A. Guertin, Vebjørn Ljoså, Peggy Anthony, and Jason Moffat for technology development and technical support that made possible the software and biological experimentation presented in this paper.
This work was supported by the Broad Institute, the MIT EECS/Whitehead/Broad Training Program in Computational Biology (NIH grant DK070069-01) supporting TRJ, DOD TSC research program grant W81XWH-05-1-0318-DS (DMS), NIH NIGMS R01 GM0725555 (DMS), NIH NIAD RO1 AI047389 (DMS), NSF CAREER award 0642971 (PG), a Merck/CSBi postdoctoral fellowship (AEC), a L'Oreal for Women in Science fellowship (AEC), a Novartis fellowship from the Life Sciences Research Foundation (AEC), a Society for Biomolecular Screening Academic grant (AEC).
The software project was conceived, designed, and developed by TRJ, AEC, PG, DMS, and IHK. IHK wrote most of the software's code, RAL and DBW prepared the cell microarrays and Drosophila follow-up experiments, and TRJ and AP made adjustments to the software. AEC and TRJ used the software to perform the screen, and together analyzed and interpreted the data. The manuscript was prepared by AEC, then revised and approved by all authors.
Electronic supplementary material
Additional file 3: Source code for CellProfiler Analyst (most recent version is available at http://www.cellprofiler.org). Additional File 3 must be uncompressed. Doubleclicking the file often activates the computer's native decompression software; otherwise download free software online (e.g., StuffIt Expander) for this purpose. (ZIP 319 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Jones, T.R., Kang, I.H., Wheeler, D.B. et al. CellProfiler Analyst: data exploration and analysis software for complex image-based screens. BMC Bioinformatics 9, 482 (2008). https://doi.org/10.1186/1471-2105-9-482
- Additional Data File
- Cell Cycle Distribution
- Sequential Gating
- Parallel Coordinate Plot
- Telophase Nucleus