Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Highly Accessed Database

Benchmark for multi-cellular segmentation of bright field microscopy images

Assaf Zaritsky1*, Nathan Manor1, Lior Wolf1, Eshel Ben-Jacob234 and Ilan Tsarfaty5*

Author Affiliations

1 Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv, 69978, Israel

2 School of Physics and Astronomy, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv University, Tel-Aviv, 69978, Israel

3 Center for Theoretical Biological Physics, Rice University, Houston, TX, 77005-1827, USA

4 Research & Development Unit Assaf Harofeh Medical Center, Zerifin, 70300, Israel

5 Department of Clinical Microbiology and Immunology, Sackler School of Medicine, Tel Aviv University, Tel Aviv, 69978, Israel

For all author emails, please log on.

BMC Bioinformatics 2013, 14:319  doi:10.1186/1471-2105-14-319


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/14/319


Received:20 June 2013
Accepted:29 October 2013
Published:7 November 2013

© 2013 Zaritsky et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Multi-cellular segmentation of bright field microscopy images is an essential computational step when quantifying collective migration of cells in vitro. Despite the availability of various tools and algorithms, no publicly available benchmark has been proposed for evaluation and comparison between the different alternatives.

Description

A uniform framework is presented to benchmark algorithms for multi-cellular segmentation in bright field microscopy images. A freely available set of 171 manually segmented images from diverse origins was partitioned into 8 datasets and evaluated on three leading designated tools.

Conclusions

The presented benchmark resource for evaluating segmentation algorithms of bright field images is the first public annotated dataset for this purpose. This annotated dataset of diverse examples allows fair evaluations and comparisons of future segmentation methods. Scientists are encouraged to assess new algorithms on this benchmark, and to contribute additional annotated datasets.

Keywords:
Collective cell migration; Wound healing assay; Segmentation; Benchmarking

Background

Characterizing and quantifying collective migration phenotypes of a monolayer of cells in vitro is an important step in understanding physiological processes such as development, wound repair and cancer motility. The prevalent approach is to acquire still or time-lapse images using bright field microscopy, followed by manual or automated extraction of quantitative measures of cellular morphology or dynamics (e.g., [1-3]).

The vast numbers of microscopic images acquired in high throughput studies preclude manual annotation and hence automatic computational tools become indispensable. Indeed, several tools to tackle these tasks were recently reported; some exploit local motion-estimation to quantify dynamic intercellular phenomena [4,5], whereas others are designed to quantify only global motion of complete colonies or confluent monolayers [6-15]. The basic common computational step in all approaches is segmentation of an image into cellular and non-cellular regions, the accuracy of which is crucial for further analysis. It is inherently a foreground-background segmentation task: no explicit cell segmentation is performed; each pixel is rather assigned a binary label as being part of either a cellular or a non-cellular region.

The high variability in imaging conditions and cells’ appearance requires robust algorithms that can deal with this imaging diversity automatically, accurately and preferably without the need for parameter-tuning. It is difficult to systematically select the most appropriate segmentation tool from the available options [16,17]. Proposed methods are usually evaluated on in-house benchmarks that are not freely available to the public. These evaluations often compare accuracy to human-annotations and rarely to alternative computational methods, hence are not subjected to a thorough comparative assessment of extant methods [18].

We therefore propose a uniform framework to benchmark algorithms for multi-cellular segmentation in bright field microscopy images.

Construction and content

A set of 171 manually segmented images of 5 different cell lines at diverse confluence levels, acquired in several laboratories under different imaging conditions, were partitioned into 8 datasets as follows (example images are presented in Figure 1, detailed description of the cells and imaging conditions can be found on the benchmark website):

• TScratch: 24 bright field images of confluent cells available at the TScratch site, http://www.cse-lab.ethz.ch/index.php?&option=com_content&view=article&id=363 webcite[6];

• Melanoma: 20 bright field images of confluent populations of brain metastatic melanoma cells acquired during a wound healing experiment [19];

• Init: 28 differential interference contrast (DIC) images of confluent DA3 cells, derived from the mouse mammary adenocarcinoma line D1-DMBA-3, acquired during wound healing experiments;

• SN15: 54 DIC images of confluent DA3 cells acquired during a multi-well wound healing experiment;

• Scatter: 6 DIC images of Madin-Darby Canine Kidney (MDCK) epithelial cells acquired during a multi-well scatter experiment (unpublished data);

• Microfluidics: 13 DIC images of MDCK cells grown in a microfluidic plate acquired during a scatter assay experiment with a Hepatocyte growth factor/scatter factor gradient (unpublished data);

• HEK293: 12 DIC images of confluent HEK293T cells acquired during a multi-well wound healing experiment (unpublished data);

• MDCK: 14 DIC images of confluent MDCK cells acquired in a multi-well wound healing experiment (unpublished data).

thumbnailFigure 1. Examples of images from the presented benchmark and their corresponding manual segmentations.

Four out of eight datasets (“TScratch”, “Melanoma”, “Init”, “SN15”) were previously reported before in our study comparing TScratch to MultiCellSeg [15]. The other four are published here for the first time. “TScratch” is the only dataset that was already publically available.

Two freely available tools (TScratch, [6]; MultiCellSeg, [15]) and one implementation of a designated algorithm [10] were evaluated on these datasets; a brief description of these algorithms is found in the Additional file 1: Supporting Text. Each of the algorithms was evaluated using the same set of parameters on all datasets to assess robustness. All raw images, manual segmentations, algorithm segmentation results, performance measurements and an evaluation Matlab script are freely available at the Broad Bioimage Benchmark Collection (BBBC, http://www.broadinstitute.org/bbbc/ webcite) [20].

Additional file 1. Supporting Text. This file contains a brief description of the evaluated algorithms, notes on parameter tuning, details on evaluation of Topman’s thresholding method, and details on assessing the baseline variance in the annotated data.

Format: PDF Size: 496KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The algorithms were ranked by mean F-measure (i.e., the harmonic mean of precision and recall) of the pixel labeled across all images in each dataset <a onClick="popup('http://www.biomedcentral.com/1471-2105/14/319/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/319/mathml/M1">View MathML</a>, where n is the number of images in the given dataset). F-measure is the evaluation measure used for foreground-background datasets in the BBBC. Evaluation results are found in Table 1, which also present the median F-measure as a more robust (less sensitive) evaluation. Additional file 2: Table S1 displays the average precision and recall measures, Additional file 3: Figure S1 plots the performance on each of the images in each dataset.

Table 1. Evaluation of the three designated tools on the eight available datasets

Additional file 2: Table S1. Precision/recall. Precision/recall of all algorithms on all datasets.

Format: PPTX Size: 43KB Download fileOpen Data

Additional file 3: Figure S1. Direct comparison of algorithms on all images. Image-by-image evaluation. Scatter plots displaying for each image the F-measure produced by the 3 algorithms. Each x-axis entry represents an image (ordered by the filename), y-axis is the F-measure. Red – Tscratch, Green – MultiCellSeg, Cyan – Topman’s algorithm. a, Init. b, NN15. c, Melanoma. d, TScratch. e, Scatter. f, Microfluidics. g, HEK293. h, MDCK.

Format: TIFF Size: 2.3MB Download fileOpen Data

To assess the baseline variance that should be expected when scoring the results of an algorithm, a second expert annotated an arbitrary partial set of the images (64 images from all datasets, excluding the “Scatter” dataset). The two annotations were evaluated using mean and median F-measures compared with the primary annotated ground truth. The results are reported in Additional file 4: Table S2, and in the Additional file 1: Supporting Text. The annotators were instructed to pay special attention to small isolated cells and voids between groups of cells; this is important for some applications and can be later excluded via post-processing if not relevant. Most inconsistencies found were in defining the exact border contours of the cells, as described in the Additional file 1: Supporting Text and in Additional file 5: Figure S2.

Additional file 4: Table S2. Baseline variance. An arbitrary partial set of the images (62 images from all datasets, excluding the “Scatter” dataset) was selected to be annotated by another expert. This annotation was compared with the primary annotated ground truth by calculating the mean F-measure to assess the baseline variance of each dataset.

Format: PPTX Size: 41KB Download fileOpen Data

Additional file 5: Figure S2. Baseline variance examples. Visualization of inconsistencies between manual annotations by different experts. Annotations shown were selected from the dataset with higher baseline variance (“Melanoma”, “Miscrofluidics”). The green channel is the raw image, the blue channel is the official annotation of cells, and the red channel is the second annotation. Thus, light-magenta represents agreement in annotation of cells, green represents agreement in annotation of non-cellular regions, light-red represents regions annotated as non-cellular in the ground truth but as cellular by the second expert, light blue represents regions that were annotated as cellular according to the ground truth but non-cellular according to the second expert. It is clear from this visualization that most inconsistencies appear at cell borders.

Format: TIFF Size: 5.1MB Download fileOpen Data

Twenty arbitrary images were selected as a training set for algorithms that apply supervised learning [15], or for adjusting parameters’ values. Use of different arbitrarily selected training images did not significantly change the algorithms’ performance. Comprehensive assessment of [10] is presented in Additional file 6: Table S3 and discussed in the Additional file 1: Supporting Text; Table 1 contains the optimal results achieved considering a single set of parameters for all datasets.

Additional file 6: Table S3. Adjusting Tompan’s algorithm. The automatic threshold extraction method in Topman’s algorithm was evaluated compared to a constant threshold. Evaluation of different values demonstrated that a constant threshold surpasses the automatic adjustment for most datasets. The best value found was used to evaluate this algorithm’s performance in the main text.

Format: TIFF Size: 734KB Download fileOpen Data

Utility

The benchmark includes two directories:

• train: 20 images (images directory) and the corresponding manual annotations (manual directory). These images can be used as a training set for algorithms that apply supervised learning or for adjusting parameters’ values.

• datasets: 8 different datasets, each consists of images (images directory), ground truth manual annotations (manual directory), results masks of the 3 algorithms we compare (tscratch, multiCellSeg, topman directories), a measures file (measures.mat) containing evaluation summary, and a second annotation for a partial set of the images.

∘ A documented Matlab script (bbbcCalcMeasures.m) that compares the different algorithms (and can easily be adjusted for comparing new algorithms).

Benchmark structure is described in the README file in the main directory.

Discussion

The wound healing assay (aka scratch assay), the traditional method used to study collective cell motility and migration [21,22] in the life sciences, is performed by inducing a sudden injury created by removal of a sheet of cells from a confluent monolayer [19]. This assay can be performed using multi-well plates, with up to 384-wells [23], providing a large amount of data for high-quality quantitative analysis. The scratch is imaged and measured periodically during the healing process, and rate of change in the wound area is recorded and can be compared with other cell lines, environmental conditions or chemical treatments. Quantifying wound healing assays is a natural application of multi-cellular segmentation algorithms. The availability of a benchmark to evaluate algorithms on a variety of cell lines and imaging conditions will enable educated algorithm selection. The general segmentation of cell clusters in bright-field images has additional applications (e.g., quantifying scatter assays [15]), thus emphasizing the importance of evaluating the segmentation of non-confluent cells images. The dataset provided is diverse in terms of cell lines, image acquisition parameters, cellular confluence levels, and was collected from several laboratories, and can thus address the need for public access to image repositories [24] as well as the general concern regarding poor algorithmic comparisons [17,18,25,26].

Conclusions

A variety of software tools and imaging apparatuses exist to enable high throughput multi-cellular segmentation in bright field images. This is the first and currently only freely available public annotated dataset for evaluations. We encourage scientists to evaluate new algorithms and to contribute additional annotated datasets to this benchmark.

Availability and requirements

All raw images, manual segmentations and evaluations are freely available at the Broad Bioimage Benchmark Collection (BBBC), http://www.broadinstitute.org/bbbc/BBBC019/ webcite.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AZ collected the data, performed the manual segmentation, evaluated the algorithms, published the data sets and wrote the manuscript. NM performed the second annotation and the ground truth evaluation. NM, LW, EBJ and IT revised the manuscript. All authors read and approved the final manuscript.

Acknowledgments

We are grateful to all colleagues who participated in data collection: Dr. Sivan Izraeli from Prof. Isaac P. Witz’s laboratory, who acquired the “Melanoma” dataset; Prof. Petros Komoutsakos, who approved publication of the images available in the TScratch website as part of this benchmark (“TScratch” dataset); Sari Natan, Dr. Doron Kaplan and Yaniv Goikhman from Prof. Ilan Tsarfay’s laboratory, who acquired the rest of the datasets. Doron Kaplan labeled regions of interest in the “MDCK” images. We thank Prof. Arieh Zaritsky for proofreading the manuscript. The work was supported in part by grants from the Breast Cancer Research Foundation; the Federico Foundation Grants; the US - Israel Binational Science Foundation and the Tauber Family Foundation at Tel Aviv University, the Center for Theoretical Biological Physics sponsored by the NSF (# PHY-0822283), and by the Cancer Prevention and Research Institute of Texas (CPRIT) at Rice University. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

References

  1. Simpson KJ, Selfors LM, Bui J, Reynolds A, Leake D, Khvorova A, Brugge JS: Identification of genes that regulate epithelial cell migration using an siRNA screening approach.

    Nat Cell Biol 2008, 10(9):1027-1038. OpenURL

  2. Petitjean L, Reffay M, Grasland-Mongrain E, Poujade M, Ladoux B, Buguin A, Silberzan P: Velocity fields in a collectively migrating epithelium.

    Biophys J 2010, 98(9):1790-1800. OpenURL

  3. Peruani F, Starruss J, Jakovljevic V, Sogaard-Andersen L, Deutsch A, Bar M: Collective motion and nonequilibrium cluster formation in colonies of gliding bacteria.

    Phys Rev Lett 2012, 108(9):098102. OpenURL

  4. Deforet M, Parrini MC, Petitjean L, Biondini M, Buguin A, Camonis J, Silberzan P: Automated velocity mapping of migrating cell populations (AVeMap).

    Nat Methods 2012, 9(11):1081-1083. OpenURL

  5. Milde F, Franco D, Ferrari A, Kurtcuoglu V, Poulikakos D, Koumoutsakos P: Cell Image Velocimetry (CIV): boosting the automated quantification of cell migration in wound healing assays.

    Integr Biol 2012, 4(11):1437-1447. OpenURL

  6. Geback T, Schulz MM, Koumoutsakos P, Detmar M: TScratch: a novel and simple software tool for automated analysis of monolayer wound healing assays.

    Biotechniques 2009, 46(4):265-274. OpenURL

  7. Korzynska A, Strojny W, Hoppe A, Wertheim D, Hoser P: Segmentation of microscope images of living cells.

    Pattern Anal Appl 2007, 10(4):301-319. OpenURL

  8. Selinummi J, Ruusuvuori P, Podolsky I, Ozinsky A, Gold E, Yli-Harja O, Aderem A, Shmulevich I: Bright field microscopy as an alternative to whole cell fluorescence in automated analysis of macrophage images.

    Plos One 2009, 4(10):e7497. OpenURL

  9. Ali R, Gooding M, Christlieb M, Brady M: Advanced phase-based segmentation of multiple cells from brightfield microscopy images. In 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro: 2008. Paris: ; 2008:181-184. OpenURL

  10. Topman G, Sharabani-Yosef O, Gefen A: A standardized objective method for continuously measuring the kinematics of cultures covering a mechanically damaged site.

    Med Eng Phys 2011, 34:225-232. OpenURL

  11. Zordan MD, Mill CP, Riese DJ, Leary JF: A high throughput, interactive imaging, bright-field wound healing assay.

    Cytometry A 2011, 79A(3):227-232. OpenURL

  12. Placidi G, Sabatini M, Gallucci M, Sotgiu A: A calculation method for semi automatic follow up of multiple sclerosis by magnetic resonance eco planar perfusion imaging.

    Stud Health Technol Inform 2003, 95:164-169. OpenURL

  13. Qin J, Chittenden TW, Gao L, Pearlman JD: Automated migration analysis based on cell texture: method & reliability.

    BMC Cell Biol 2005, 6(1):9. OpenURL

  14. Skopin MD, Molitor SC: Effects of near-infrared laser exposure in a cellular model of wound healing.

    Photodermatol Photoimmunol Photomed 2009, 25(2):75-80. OpenURL

  15. Zaritsky A, Natan S, Horev J, Hecht I, Wolf L, Ben-Jacob E, Tsarfaty I: Cell motility dynamics: a novel segmentation algorithm to quantify multi-cellular bright field microscopy images.

    Plos One 2011, 6(11):e27593. OpenURL

  16. Evanko D: All things being equal.

    Nat Methods 2012, 9(2):111. OpenURL

  17. Treloar KK, Simpson MJ: Sensitivity of edge detection methods for quantifying cell migration assays.

    PLoS One 2013, 8(6):e67389. OpenURL

  18. Smith R, Ventura D, Prince JT: Novel algorithms and the benefits of comparative validation.

    Bioinformatics 2013, 29(12):1583-1585. OpenURL

  19. Liang CC, Park AY, Guan JL: In vitro scratch assay: a convenient and inexpensive method for analysis of cell migration in vitro.

    Nat Protoc 2007, 2(2):329-333. OpenURL

  20. Ljosa V, Sokolnicki KL, Carpenter AE: Annotated high-throughput microscopy image sets for validation.

    Nat Methods 2012, 9(7):637. OpenURL

  21. Lampugnani MG: Cell migration into a wounded area in vitro.

    Methods Mol Biol 1999, 96:177-182. OpenURL

  22. Zahm JM, Kaplan H, Herard AL, Doriot F, Pierrot D, Somelette P, Puchelle E: Cell migration and proliferation during the in vitro wound repair of the respiratory epithelium.

    Cell Motil Cytoskeleton 1997, 37(1):33-43. OpenURL

  23. Yarrow JC, Perlman ZE, Westwood NJ, Mitchison TJ: A high-throughput cell migration assay using scratch wound healing, a comparison of image-based readout methods.

    BMC Biotechnol 2004, 4:21. OpenURL

  24. Swedlow JR: Finding an image in a haystack: the case for public image repositories.

    Nat Cell Biol 2011, 13(3):183. OpenURL

  25. Boulesteix AL: On representative and illustrative comparisons with real data in bioinformatics: response to the letter to the editor by Smith et al.

    Bioinformatics 2013, 29(20):2664-2666. OpenURL

  26. Boulesteix AL, Lauer S, Eugster MJ: A plea for neutral comparison studies in computational sciences.

    PLoS One 2013, 8(4):e61562. OpenURL