Hindawi Publishing Corporation BioMed Research International Volume 2014, Article ID 863564, 6 pages http://dx.doi.org/10.1155/2014/863564

Research Article CELLCOUNTER: Novel Open-Source Software for Counting Cell Migration and Invasion In Vitro
Xiaoni Li,1 Hongshun Yang,2 Hailiang Huang,3,4 and Tao Zhu1
1

Hefei National Laboratory for Physical Sciences at Microscale and School of Life Sciences, USTC, Hefei, Anhui, China School of Information Science and Technology, USTC, Hefei, Anhui, China 3 Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, MA, USA 4 Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
2

Correspondence should be addressed to Hailiang Huang; hhuang@atgu.mgh.harvard.edu and Tao Zhu; zhut@ustc.edu.cn Received 16 April 2014; Revised 7 June 2014; Accepted 7 June 2014; Published 26 June 2014 Academic Editor: Calvin Yu-Chian Chen Copyright © 2014 Xiaoni Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Transwell Boyden chamber based migration/invasion assay is a simple and extensively used approach for the characterization of cell motility in vitro. Cell motility is quantified by counting the number of cells that pass through the filter membrane. The counting is usually performed manually, which is laborious and error prone. We have therefore developed CELLCOUNTER, an application that is capable of recognizing and counting the total number of cells through an intuitive graphical user interface. The counting can be performed in batch, and the counting results can be visualized and further curated manually. CELLCOUNTER will be helpful in streamlining the experimental process and improving the reliability of the data acquisition.

1. Introduction
Cell migration is the movement of cells from one location to another generally in response to and toward specific external chemical signals. Cell invasion is similar to cell migration, except that it requires the cell to migrate through an extracellular matrix or basement membrane barrier by enzymatically degrading the barrier. Cell migration/invasion is central to many physiological and pathological processes such as embryonic development, wound repair, and tumor metastasis [1–3]. Transwell Boyden chamber [4] based cell migration/invasion assay is a simple and extensively used approach for the quantitation of cell motility in vitro [3, 5]. The number of cells that pass through the filter membrane from Boyden chamber is usually counted manually from the inverted microscopic images. Such images may contain hundreds of cells and manually counting the number of them is not a trivial work. Existing image analysis programs, for example, CELLPROFILER and IMAGEJ [6–8], require pipeline/macro/plugin files that are specific to the cell/assay types, which is not

available yet for the Transwell assays. Although it is possible to create a CELLPROFILER/IMAGEJ pipeline/macro/plugin for the Transwell assays, an independent program taking into account the specific characteristics of the assay can most likely perform better [9]. Another related application [10], although available for some migration assays, requires fluorescently stained cells and does not work for the Transwell assays. We have therefore developed CELLCOUNTER, an application specialized in automatically counting the number of cells in Transwell assays. This application supports all major image formats (JPEG, PNG, and GIF) and offers an intuitive graphical user interface (GUI). The default parameters in CELLCOUNTER work well for most users, while the more advanced users have the option to set the parameters of their choice to further improve the counting accuracy (Figure 1). CELLCOUNTER is capable of processing hundreds of images in a batch and exporting the results to a plain text file for statistical analysis. The recognized cells will be marked on the assay images for curation purpose. Users may manually add or remove recognized cells through the GUI.

2

BioMed Research International

1 6

2

3

4

5

7

Figure 1: The graphical user interface of CELLCOUNTER. Green circles indicate automatically recognized cells and red circles indicate manually labeled cells. Functions for buttons/inputs are (1) manual counting/curation; (2) zooming in; (3) zooming out; (4) restoring to the original image size; (5) saving the counting results (including images); (6) setting the radius of the label (for visualization only); (7) setting the radius of small wells (in px, see the main text for more details, default recommended).

The counting process, including the automatic recognition and manual curations, can be saved to a disk file and loaded back to resume the work. CELLCOUNTER will prove to be useful in streamlining the analysis of Transwell assays, reducing human errors, and improving the reliability of the assay quantitation.

the non-stimulated cells, we added ethyl alcohol as a control compound. After the stimulation we incubated the cells for 4 hours. 2.4. Existing Software. CELLPROFILER and IMAGEJ are versatile and flexible image processing software for scientific purposes. For comparison purpose, we used a na¨ve setup and ı choice of parameters. CELLPROFILER uses a customized pipeline, including input and analysis modules, to process images. The input modules in our pipeline are “Images,” “Metadata,” “NamesAndType,” and “Groups.” The analysis modules are “ColorToGray,” “ApplyThreshold,” and “IdentifyPrimaryObjects.” There are many tunable parameters in each module. We chose the parameters that make the most sense to us. The project configuration file used in this study (CellProfiler.cpproj) with all parameters is available to download at https://bitbucket.org/linora/cellcounter/downloads. For IMAGEJ, after loading the images, we first adjusted the color threshold using the default thresholding method (Image → Adjust → Color Threshold → select). We then analyzed the images using Analyze → Analyze particle → summarize. 2.5. Implementation. CELLCOUNTER was developed using the C++ programming language. Some image processing

2. Materials and Methods
2.1. Cell Lines. 20 images from 6 human cancer cell lines were used in this study (Table 1): 17 images are from breast cancer cell lines (1 x T47D, 13 x MCF-7, 2 x MDA-MB-231, and 1 x MDA-MB-435s), 2 images are from hepatic carcinoma (Hep3B), and 1 image is from prostatic cancer (pc-3). All cell lines were obtained from the American Type Culture Collection (Rockville, MD, USA). 2.2. Transwell Assay. Assays were performed in BioCoat Matrigel invasion chambers (Corning Costar, Acton, MA) as described previously [11–14]. The assay images were captured using a Nikon camera with 5 M pixels. 2.3. Biology Stimulation. Five images (out of the 20 total images) are from stimulated cell (4 x MCF-7 and 1 x MDAMB-231). We stimulated the cells with 20 ng/mL epidermal growth factor (EGF) when the confluency reached 95%. For

BioMed Research International

3

Table 1: Images used in this study. Cells were stained with 0.1% crystal violet for 20 minutes for the deep staining group and 8 minutes for the light staining group. High density images refer to those having over 400 cells and low density images refer to those having less than 200 cells. Biology stimulation of cells was discussed in Section 2. ID 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Cell line T47D MCF-7 MDA-MB-231 pc-3 MCF-7 MCF-7 MDA-MB-231 Hep3B MDA-MB-435s Hep3B MCF-7 MCF-7 MCF-7 MCF-7 MCF-7 MCF-7 MCF-7 MCF-7 MCF-7 MCF-7 Stimulation No No No No — Yes Yes — — — — Yes No Yes — Yes — — — — Stain — Deep — — Deep Deep — Light — — Light — Deep Light Light — — Light Deep — Density Expert High High High — — — High — — — Low — — — Low Low Low — High Low 931 792 475 523 314 359 566 284 464 186 165 403 417 189 94 120 139 8 554 191 Number of cells CELLCOUNTER 957 748 479 507 322 364 572 279 461 178 159 403 395 187 98 118 140 8 553 192 Beginner 955 538 413 390 263 300 366 249 430 144 — — — — — — — — — —

algorithms were implemented with the help of Open Source Computer Vision Library [15], and the graphical user interface was designed using Nokia Qt (a cross platform application and UI framework). CELLCOUNTER is compatible with all major platforms including Windows (XP or newer), Mac OSX (10.6+), and Linux (Kernel 3.x). Assays images for CELLCOUNTER have to have an aspect ratio of 4 : 3 (default in most digital cameras) and a minimal resolution of 1280 px by 960 px. Images that have the correct aspect ratio but in higher resolutions will be resized to 1280 px by 960 px. Resizing to a lower resolution reduces the use of memory and CPU time in analyzing images. We have verified that the resolution of 1280 px by 960 px is high enough to retain all the necessary details for cell counting. In CELLCOUNTER, the original assay image is first converted to a grayscale image, and further partitioned into cell areas and background areas using an adaptive threshold. Areas that have grayscale below the threshold (close to white)

are background and areas that have grayscale above the threshold (close to black) are cells. The threshold is established by using the maximum entropy method [16]. A small value can be added to the threshold (threshold adjustment) to fine-tune the counting accuracy. Users may adjust this value in GUI and use the visual feedback to choose the optimal value. Alternately, the default value works very well in most circumstances. Standard image processing procedures, such as contrast enhancement, smoothing, eroding, and dilating [17], are also performed to remove noises in the image. Due to the design of the assay, small wells in the equipment can be captured in the images and counted falsely as cells. We solve this problem by further partitioning the cell areas into small wells and true cells. Another adaptive threshold is established by applying the OSTU method [18] to only the cell areas identified in the previous step. Because the small wells are generally darker than the true

4
1200 1000 Cell counting results 800 600 400 200 0 1 2 3 4 5 6 7 8 9 10

BioMed Research International the numbers counted by experts (2-way ANOVA ������ value = 0.91). In contrast, the beginners gave slightly smaller cell counts (2-way ANOVA ������ value = 0.04), and these cell counts have a significantly higher standard deviation (Wilcoxon test ������ value = 0.004) than the cell counts from the experts. Therefore, we conclude that CELLCOUNTER is able to perform accurate cell counting and improve the stability of the counting results. Additional tests using unpublished images have confirmed this conclusion (results not shown). Systematic biases can lead to false positive findings. For example, if software systematically reports a larger (comparing with the true) number of cells if they are stained for a longer time, and if the person who stains the case samples tends to stain them longer, we might falsely claim a significant case/control difference. Therefore, we determined CELLCOUNTER for systematic biases towards the staining, the density of cells, and the biology stimulation. We found no systematic bias in all 3 conditions (Table 1). The ������ values are 0.42 for deep versus light stains, 0.91 for high versus low cell densities, and 0.96 for stimulated versus nonstimulated cells. CELLPROFILER and IMAGEJ are powerful scientific image processing software. They are versatile and flexible but also have steep learning curves and require efforts to set up the right pipeline/parameters. Using a na¨ve setup ı as discussed in Section 2, we noticed that both methods perform worse comparing with CELLCOUNTER (Figure 3). The correlation coefficients (������2 ) for CELLPROFILER and IMAGEJ are 0.49 and 0.12, respectively, while the ������2 for CELLCOUNTER reached almost 1. CELLCOUNTER performs better because it was designed for images from these assays, while CELLPROFILER and IMAGEJ have more general purposes. We are aware that more sophisticated setup or choices of parameters for CELLPROFILER and IMAGEJ may improve their performance. However, we argue that our choices are representative of average users with basic training and knowledge in image processing.

Experts CELLCOUNTER Beginners

Figure 2: Cell counting results for 10 assay images from 3 experts, CELLCOUNTER, and 3 beginners. Error bars are standard deviations (only available for human counting results).

cells, applying this adaptive threshold, followed by image smoothing, eroding, and dilating, can successfully remove the small wells from the cell areas. One cell area may contain multiple overlapping cells. We count the number of cells in a particular cell area by using the radius of its maximum inscribed circle. If the radius is smaller than the empirical threshold of the cell radius r (default to 6 pixels) and greater than the empirical threshold for noise (default to 4 pixels), we count this cell area as a single cell. If the radius is greater than the empirical threshold (r), we count the number of cells in this cell area as ℎ×[ ������ ������ + ������] × [ + ������] , 2������ 2������ (1)

in which L and W are the length and width of the minimum bounding rectangle, h is the number of layers that the cells stack, and a is a parameter to account for cells at the boundary of the rectangle. h and a are unknown and have been estimated using the human counted cell numbers as the training dataset (ℎ ≈ 4.0 and ������ ≈ 0.5). The total number of cells in the image is then the summation of cell numbers in each cell area.

4. Conclusion
We have developed CELLCOUNTER, a program that features an intuitive graphical user interface to count the number of cells in Transwell Boyden chamber based migration/invasion assays. This program allows high-throughput analysis of a large number of assay images. The counted cells are visibly marked on the assay images and can be manually curated. The accuracy of the counting results has been validated using expert counted cell numbers as the gold standard. CELLCOUNTER significantly simplifies the data acquisition in the Transwell assays, reduces human errors, and improves the stability of counting results. It will prove to be a helpful tool in the study of cell invasion and metastasis in vitro. CELLCOUNTER is currently available for Windows, Mac OSX, and Linux platforms and can be downloaded from https://bitbucket.org/linora/cellcounter/downloads.

3. Results and Discussion
The user interface and a counting example using CELLCOUNTER is provided in Figure 1. We used CELLCOUNTER to analyze 10 assay images (randomly selected from the 20 images in Section 2). The cell numbers counted by the experts (3 researchers that are proficient in cell counting) are used as the gold standard and compared with the cell numbers counted by CELLCOUNTER (default parameters, without manual curation) and beginners (3 researchers with basic trainings in cell counting). As shown in Figure 2, we found that the cell numbers counted by CELLCOUNTER are statistically the same as

BioMed Research International
1200 1400 1200 1000 800 600 400 200 0 R2 = 0.490

5

CELLCOUNTER

CELLPROFILER

Cell counting results of CELLCOUNTER

R = 0.996 800

2

600

400

200

0

0

200

400 600 800 Cell counting results of experts
(a)

1000

Cell counting results of CELLPROFILER

1000

0

200

400 600 800 Cell counting results of experts
(b)

1000

2000 1800 1600 Cell counting results of IMAGEJ 1400 1200 1000 800 600 400 200 0 0 200

IMAGEJ

R2 = 0.117

400 600 800 Cell counting results of experts
(c)

1000

Figure 3: Counting results from CELLCOUNTER (a), CELLPROFILER (b), and IMAGEJ (c) comparing with human counting results.

Conflict of Interests
The authors declare that there is no conflict of interests.

supported by the National Key Scientific Program of China (2012CB934002, 2010CB912804) and National Natural Science Foundation of China (81272925).

Acknowledgments
The authors thank Dr. Zhiming Pi for helpful advice to improve this software. They would also like to thank the members of their laboratories for contributing to the development of the software and this paper. This work was

References
[1] A. J. Ridley, M. A. Schwartz, K. Burridge et al., “Cell migration: integrating signals from front to back,” Science, vol. 302, no. 5651, pp. 1704–1709, 2003.

6
[2] S. A. Courtneidge, “Cell migration and invasion in human disease: the Tks adaptor proteins,” Biochemical Society Transactions, vol. 40, no. 1, pp. 129–132, 2012. [3] E. A. Goncharova, D. A. Goncharov, and V. P. Krymskaya, “Assays for in vitro monitoring of human airway smooth muscle (ASM) and human pulmonary arterial vascular smooth muscle (VSM) cell migration,” Nature Protocols, vol. 1, no. 6, pp. 2933– 2939, 2006. [4] H. Chen, “Boyden chamber assay,” Methods in Molecular Biology, vol. 294, pp. 15–22, 2005. [5] X. H. Liu, K. H. Lu, K. M. Wang et al., “MicroRNA-196a promotes non-small cell lung cancer cell proliferation and invasion through targeting HOXA5,” BMC Cancer, vol. 12, article 348, 2012. [6] A. E. Carpenter, T. R. Jones, M. R. Lamprecht et al., “CellProfiler: image analysis software for identifying and quantifying cell phenotypes,” Genome Biology, vol. 7, no. 10, article R100, 2006. [7] M. R. Lamprecht, D. M. Sabatini, and A. E. Carpenter, “CellProfiler: free, versatile software for automated biological image analysis,” BioTechniques, vol. 42, no. 1, pp. 71–75, 2007. [8] C. A. Schneider, W. S. Rasband, and K. W. Eliceiri, “NIH Image to ImageJ: 25 years of image analysis,” Nature Methods, vol. 9, no. 7, pp. 671–675, 2012. [9] N. N. Kachouie, L. Kang, and A. Khademhosseini, “Arraycount, an algorithm for automatic cell counting in microwell arrays,” BioTechniques, vol. 47, no. 3, pp. 1–100, 2009. [10] B. K. Al-Khazraji, P. J. Medeiros, N. M. Novielli, and D. N. Jackson, “An automated cell-counting algorithm for fluorescentlystained cells in migration assays,” Biological Procedures Online, vol. 13, no. 1, article 9, pp. 1–9, 2011. [11] S. Mukhina, H. C. Mertani, K. Guo, K. Lee, P. D. Gluckman, and P. E. Lobie, “Phenotypic conversion of human mammary carcinoma cells by autocrine human growth hormone,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 42, pp. 15166–15171, 2004. [12] X. Li, X. Liu, W. Xu et al., “C-MYC-regulated miR-23a/242/27a cluster promotes mammary carcinoma cell invasion and hepatic metastasis by targeting sprouty2,” Journal of Biological Chemistry, vol. 288, no. 25, pp. 18121–18133, 2013. [13] P. Qian, Z. Zuo, Z. Wu et al., “Pivotal role of reduced let-7g expression in breast cancer invasion and metastasis,” Cancer Research, vol. 71, no. 20, pp. 6463–6474, 2011. [14] J. Chen, Y. Yao, C. Gong et al., “CCL18 from tumor-associated macrophages promotes breast cancer metastasis via PITPNM3,” Cancer Cell, vol. 19, no. 4, pp. 541–555, 2011. [15] G. Bradski, “The OpenCV library,” Dr Dobbs Journal, vol. 25, no. 11, pp. 120–126, 2000. [16] S. F. Gull and J. Skilling, “Maximum-entropy method in imageprocessing,” IEE Proceedings F: Communications, Radar and Signal Processing, vol. 131, no. 6, pp. 646–659, 1984. [17] M. Sonka, V. Hlavac, and R. Boyle, Image Processing, Analysis, and Machine Vision, PWS, Pacific Grove, Calif, USA, 2nd edition, 1999. [18] H. Ng, “Automatic thresholding for defect detection,” Pattern Recognition Letters, vol. 27, no. 14, pp. 1644–1649, 2006.

BioMed Research International