This paper describes the quantification of hemocytometer and migration/invasion micrographs through two new open-source ImageJ plugins Cell Concentration Calculator and migration assay Counter. Furthermore, it describes image acquisition and calibration protocols as well as discusses in detail the input requirements of the plugins.
The National Institute of Health's ImageJ is a powerful, freely available image processing software suite. ImageJ has comprehensive particle analysis algorithms which can be used effectively to count various biological particles. When counting large numbers of cell samples, the hemocytometer presents a bottleneck with regards to time. Likewise, counting membranes from migration/invasion assays with the ImageJ plugin Cell Counter, although accurate, is exceptionally labor intensive, subjective, and infamous for causing wrist pain. To address this need, we developed two plugins within ImageJ for the sole task of automated hemocytometer (or known volume) and migration/invasion cell counting. Both plugins rely on the ability to acquire high quality micrographs with minimal background. They are easy to use and optimized for quick counting and analysis of large sample sizes with built-in analysis tools to help calibration of counts. By combining the core principles of Cell Counter with an automated counting algorithm and post-counting analysis, this greatly increases the ease with which migration assays can be processed without any loss of accuracy.
In vitro cell counting is an important basic technique in a wide range of tissue culture experiments. Accurately determining the number of cells in a culture is essential for experimental reproducibility and standardization1,2. Cell counting can be performed manually using a hemocytometer as well as using a variety of automated methods, each with their own advantages and disadvantages3,4,5. Most of the automated methods for cell counting belong to one of two classes, those that use the Coulter principle or flow cytometry. Coulter counters take advantage of cells electrical resistance to determine cell number and size. They are fast, accurate and cheaper than flow cytometers. However, they are rarely used for only cell counting due to their considerable cost compared to manual counting3. Flow cytometers, on the other hand, are expensive but they have many applications such as cell counting, analysis of the cells shape, structure and measuring internal cell markers4,5. Machines that use either of these two principles are available from many manufacturers. Manual counting is affordable but time-consuming and subject to bias while the automated methods come with a fraction of the time required for the manual counting but using expensive machines6.
Other common cell culture procedures are in vitro cell motility assays, namely, cell migration and invasion7. Migration and invasion assays are commonly used to investigate cell motility and invasiveness in response to a chemotactic response. In addition, they are widely used to study embryonic development, differentiation, inflammatory response, and metastasis of multiple cell types7-11. Cells that have migrated or invaded through the porous membrane of a migration assay can be quantified in two different ways. Firstly, by staining the cells with a fluorescent dye, dissociation from the membrane, and quantification using a fluorescent reader12. A limitation of this method of quantification is that no record can be retained of the membranes and there is no possibility for further analysis13. The second quantification method is for migrated/invaded cells to be fixed and stained with fluorescent dye or more commonly, with cytological dyes such as crystal violet, toluidine blue dye or hematoxylin; then cells are quantified manually using inverted microscopic images of these membranes which is a very time-consuming task12,13.
To overcome the drawbacks of manual cell counting, two reliable and accurate automated cell counters for cell concentration and for the migration assay were developed. These automated cell counter algorithms were developed for ImageJ as a plugin using Oracle's Java computer language. ImageJ is a public and widely-used image processing tool developed by the National Institute of Health (NIH)14,15; thus, writing these plugins for ImageJ facilitates easy integration into the biological community.
Automation of cell counting ensures high throughput and reproducibility compared to manual counting. Although other available software and plugins can be used to calculate cell concentration through image analysis5,16,17, Cell Concentration Calculator plugin is fast and can also handle dilutions of cells and treatments. Moreover, all results and calculations from these two counters can be saved and exported. The two plugins described in this paper are optimized for the use of a phase contrast microscope for live cell imaging and large field of view (entire membrane capture) imaging for migration assay membranes through the use of a dissecting scope. The plugins are freely available for download with installation instructions from: http://peng.lab.yorku.ca/imagej-plugins.
1. Compound Microscope and Camera Setup (Cell Concentration Calculator)
2. Image Volume Calibration
3. Camera Exposure Calibration
4. Image Acquisition
5. Image Counting and Dilutions
6. Migration and Invasion (Counter)
7. Dissecting Scope and Camera Setup
8. Image Acquisition and Flatfield
9. Configuration Settings
10. Counting Images and Calibration
11. Saving/Opening Results and Exporting to CSV
Cell Concentration Calculator
Figure 1 presents the overall process of CCC calibration and countable image acquisition. Figure 1A and 1B depict the P-square calibration image and calculation of P-square length in pixels. CCC determines cell concentration in a given volume using the formula:
A hemocytometer's P-square has a volume of 100 nL (1 mm x 1 mm x 0.1 mm) and given this constant, the total image volume can be calculated after converting pixels to mm. Figure 1C is an ideal countable image with in-focus cells displaying the characteristic phase illumination and no background hemocytometer graduations. Finally, the scatter plot of Figure 1D shows a highly correlated, manual versus automated count of cells from 57 images taken over various experiments and cell types including HTR8, ES-2, and Swan71. Of note, the upper range of concentration was ~5.6 x 106 cells/mL with a low end of ~2.3 x 103 cells/mL (≈ 5 cells per image). It is suggested that if counts are below 10-15 cells per image, the sample should be suspended in a smaller volume to increase statistical power.
Acquisition and Processing of Images for Migration Assay Counter
Hundreds of migration assay membranes were imaged over many experiments, all of which were seeded with the Human trophoblast cell line HTR8, Swan71, or ovarian cancer cell line ES-2. Of these images, multiple were chosen to represent a range of categories from very poor to excellent quality based on brightness, clarity and color of stained cells, and the degree of background staining and unwanted particles (noise). Using these images, the default RGB Threshold color settings (≈ 150, 120, 0) were determined (Figure 2C) and used as a baseline in all subsequent developments of the algorithm. The goal was to maximize the nuclear color to total color ratio, i.e., the majority of colored pixels should be within cell nuclei. The upper panel in Figure 2A depicts the ideal image brightness, cell nuclei clarity and color, and negligible background noise. The brightness of the image is important to ensure there is a great enough contrast between cell and background to produce a black and white binary image.
If this difference is not met, large or unpredictable areas may be counted by the ImageJ Analyze Particles function; the best case scenario is a completely white background. In opposition, the lower panel of Figure 2A has extreme background staining and cells that are nearly indistinguishable. Membranes with this degree of staining will most likely produce unreliable results from the migration assay counter.
In some instances, increasing brightness to the ideal level may affect image fidelity by overexposing the image. Similarly, high exposure or brightness may produce systemic chromatic effects with progressive or irregular light and dark regions. With flatfield correction, these effects can be minimized or removed entirely as shown in Figure 2B (upper vs. lower panel). Furthermore, flatfield correction is a good way to equalize the brightness of multiple membrane images.
Calibration and Validation of Migration Assay Counter
To assist the user in better determining if migration assay membrane images meet the required criteria for accurate counting, two qualifiers were designed, called image quality (Q), and calibration recommendation (CR). Importantly, both qualifiers, as the name CR suggests, are recommendations only and act as guides rather than absolute judges of the countability of each image. Both Q and CR are based on the metrics of the frequency scatter plot of particle area (10.3.2). For simplification, a metric can be regarded as the various shapes of the curve that make up the frequency plot. The desired metric of an adequate Q (≥0.5) was determined by the observed approximate normal distribution of cell size. Commonly, cells that are unresolvable from each other, over-stained, or just particularly large, shift the skewedness to a right-tailed normal distribution (Figure 3B). As such, this is the ideal metric of a typical calibrated image. With the addition of background noise, this normally leads to a large number of particles in the 1-5 pixel area range (Figure 3A). In order to calculate the plot metrics, the data are fitted with ten Savitzky-Golay smoothed curves of differing polynomial degrees produced by Dr. Michael Thomas Flanagan's Java Scientific Library (http://www.ee.ucl.ac.uk/~mflanaga/java/). Essentially, this creates multiple points in the approximate region of each extremum. Through a series of density cluster maps of local minima and maxima, a general picture of the plot metrics can be computed as an ordered list, i.e., minimum, maximum, minimum/maximum overlap, …, nth extremum. In brief, the degree to which the smoothed curves fit the frequency data determines how tightly packed the extrema points are. The greater the clustering of non-overlapping extrema, the greater Q becomes. In theory, Q represents overall image clarity based on the distribution of distinct particle sizes.
Whether the Boolean CR is true or not depends on the sequence of extrema. From Figure 3A, the ordered list would be minimum, maximum, and a sequence of overlapping extrema based on the properties of the Savitzky-Golay curve. Since Figure 3A represents an uncalibrated image, this general sequence of extrema flags CR as true, suggesting to the user that the image may need calibration. Moving forward, it can be seen that from Figure 3B, this list would be identical but with the exclusion of the first minimum. Thus the metrics of ideal uncalibrated and calibrated images were determined. Deviations from these metrics will generally flag an image for calibration and reduce Q ≤ 0, such as the case with the high background noise membrane in Figure 3C. Applying these qualifiers to a selection of modest to excellent images, with regards to brightness and background noise, it is clear that calibrated image counts are significantly closer to manual counts than uncalibrated (1.9% ± 0.3 vs. 21.7% ± 2.9, respectively; Figure 3D, E). Given the overall high quality of images chosen for analysis (uncalibrated = 0.71 ± 0.04; mean ± SE), there was an expected modest increase only in Q following calibration ( = 0.90 ± 0.04). Taken together, Q suggests the overall countability potential of an image whereas CR is a strong indicator of whether calibration is successful or not.
Timing Comparisons of Both Cell Concentration Counter and Migration Assay Counter
Two researchers (User 1 and User 2) were used to test the speed comparisons between manual and automated methods. Both researchers had substantial experience with cell counting and the migration assay but User 1 was experienced in the usage of both CCC and TC while User 2 had to follow written instructions. Figure 4A compares CCC calibration time between User 1 and 2. As expected, User 1 was substantially faster than User 2, together, taking on average approximately five minutes for CCC calibration. In order to compare manual hemocytometer counting and automated rates, the manual rate of counting using a tally counter was compared to the time it took to take nine images of the hemocytometer chamber and counted in CCC (Figure 4B). A typical cell concentration of 1.15 x 106 cells/mL was used to compare timings, leading to an average of 1x increase in throughput. This rate will vary depending on the number of cells loaded into the hemocytometer as the total time taken to capture images and process them is independent of cell number.
Lastly, timings encompassing image acquisition of membranes, quantification within TC, and manual adjustments of these images was measured in Figure 4C. Notably, 12 migration assay membranes with low cell number (Total = 10,571 cells) and substantial background staining and cellular debris were chosen to facilitate a worst case scenario that would require manual adjustments to cell number. This is reflected in Figure 4C adjustment (Adj) column where it took an average of 13 min to remove unwanted counts and add missed cells. For comparison to manual counting, optimal cell counting rates were determined with Cell Counter; high quality membrane images with high cell density were used (results not shown). This yielded an average maximal rate between User 1 and 2 of 9.1 x 103 cells/h (~2.5 cells/s). Using these numbers, the membranes from Figure 4C were counted 4.4x faster with TC than would be expected at maximal manual rate. The time savings are directly dependent on cell number and image quality, by counting migration membranes that required little or no adjustment and high cell density (~7,000 cells/membrane), TC generated cell counts 1,395x faster than the maximal manual rate.
Figure 1: Cell Concentration Calculator. (A) Phase contrast micrograph of the central primary square of the hemocytometer taken at 40X total magnification. Scale bar represents 200 µm. (B) A cropped version of the image from (A) depicting what length to measure on the hemocytometer. The Results window Length column was highlighted for easy identification. Yellow bar = 1 mm. (C) An ideal micrograph of a hemocytometer containing HTR8 cells after microscope and software calibration at 40X total magnification with a resolution of 1,600 x 1,200 px. Scale bar = 200 µm. (D) A scatter plot comparing manual counts using the ImageJ plugin Cell Counter compared to automated counts of the same images using Cell Concentration Calculator. n = 57 images. Please click here to view a larger version of this figure.
Figure 2: Migration assay counter representative images. (A) The upper panel is a one eighth portion of a total membrane depicting an ideal image with minimal background; Scale bar = 200 µm. The lower panel contains an image with significant background staining that severely affects the accuracy of the automatic counting; Scale bar = 585 µm. (B) The upper panel represents a dark image of good quality which produced inaccurate counts. The bottom panel is a countable version of the same image after flatfield correction. Scale bars = 593 µm. (C) The upper panel represents a zoomed-in view of a typical membrane. The lower panel is the same image followed by the desired ImageJ color thresholding (RGB Threshold = 150, 120, 0) to remove intercellular background and the majority of the visibly stained cytoplasm. All images were taken at 1.35X total magnification with a calibrated dissecting scope at a resolution of 2,592 x 1,944 px from multiple migration assay invasions of HTR8 stained with hematoxylin. Scale bars = 100 µm. Please click here to view a larger version of this figure.
Figure 3: Calibration and validation of migration assay counter. (A) Typical example of an uncalibrated image of high quality. (B) The calibrated version of (A) using migration assay counter's Recount > 'Suggested size'. (C) A typical frequency plot of a low image quality membrane with significant background stain or overall darkness. (D) A scatter plot comparing manual counts using the ImageJ plugin Cell Counter and migration assay counter automated counts. (E) A bar graph showing the average percent difference of calibrated versus uncalibrated automated counts. Data are Mean ± SE. P < 0.0001, n = 30, unpaired t-test with Welch's correction. The same images were used for D and E. Calibration of both was done with the Recount > 'Suggested size' and Recount > 'Manual settings' to further refine the minimum size counted and color threshold. No manual count adjustments were made using the 'Open image with counts' function. Please click here to view a larger version of this figure.
Figure 4: Timing of Cell Concentration Counter and migration assay counter usage. (A) Comparison of CCC calibration times (steps 1 and 2) between User 1 (CCC/TC expert) and User 2 (CCC/TC novice). (B) Manual cell counting times were averaged over five trials and expressed in cells counted per min. Automatic counts were calculated by averaging times to capture nine images of a hemocytometer chamber and counted with CCC. Cell concentration used was ~1.15 x 106 cells/mL. Data are Mean ± SE. n = 5. (C) Transwell Counter timings were compared over three categories: acquisition (Acq; steps 7 and 8), quantification (Qnt; steps 10.1-10.3.3 using default configuration settings), and manual adjustment (Adj; steps 10.4-1.4.2). The membranes quantified had a high degree of background staining and debris in order to necessitate substantial manual adjustment for accurate counts. Please click here to view a larger version of this figure.
Critical Steps, Troubleshooting, and Limitations
The very nature of automated computational methods, specifically those of particle analysis, necessitates the mathematical ability to define these particles. Consequently, the accuracy of both Cell Concentration Calculator and migration assay counter is majorly dependent on image fidelity, that is, how closely the captured image resembles the cell sample or migration assay membrane. It is therefore of the upmost importance to follow microscope and associated software calibration protocols as best as possible. This includes limiting background noise, reducing unwanted particles, capturing bright and uniformly in-focus images, and saving in non-lossy file formats such as tiff. Thus, both plugins are likely to produce erroneous results if these requirements are not met. It is therefore always good practice to double check cell counts to determine if they match visual expectations when in doubt. If indeed the results do not match, comparing the original image with the binary image may help elucidate the issue (10.4 note). In some cases, darker migration assay images may contain large areas that have been counted due to pixel values falling within the color threshold limits. In general, using a flatfield image will remove darkness and blotchiness and prevent most unwanted counts due to chromatic irregularities.
Given these possible and relatively common issues, it is important to follow standard image integrity guidelines like those proposed by the Nature Publishing Group and others13. To keep counts consistent, any adjustment to one image should be applied equally to all others. This includes but is not limited to contrast, saturation, brightness, gain, filters such as averaging and sharpness, and color filters. Adjustment of gamma should be avoided as it applies a non-linear change of pixel color and so may affect each image differently14. When applying a common exposure time to images with few cells (<1,000), this can lead to overexposure and loss of image fidelity. Conversely, a membrane with thousands of cells can require higher exposure times to prevent underexposure. Accordingly, it is best practice to adjust images as little as possible, and only when necessary, to maintain the highest image fidelity achievable.
Even with a high fidelity image, true particle counts can be severely skewed by unwanted particles. When using Cell Concentration Calculator, if a sample contains significant amounts of cell debris, specifically of areas similar to those of suspended cells, this may skew the result. In some cases, this could prevent automated analysis if the debris cannot be minimized. Likewise, migration assay membranes that contain high levels of background staining of similar color to stained cells or unremoved cells from the backside of the membrane, will likely produce inaccurate results. These unwanted particles can normally be removed in a few ways: increasing the light source brightness or exposure time, modifying the color thresholding to be more stringent, or manually removed via 'Open image with counts' (10.4). It is important to note that this protocol produces the optimal quality of image required for CCC and TC to maintain accuracy. However, TC is capable of counting images of significantly less quality, varying magnifications, or different color stains, generally only requiring more time spent on manual adjustments (10.4).
During calibration of any migration assay membrane image it is important to take into consideration the resolution of the image and the relative size of the cells. As previously described, Image Quality (Q) and Calibration recommendation (CR) are dependent on frequency plot metrics of particle area. Of the images analyzed, each cell takes up on average 1/100,000 of the total image area. When using images of just millions of pixels with a small cell size ratio, the variation in cell size is also small. That is, the variance may only range from 10-60 px. But as the cell area to image area ratio gets larger, the distribution of cell sizes increases, reducing the Kurtosis of the cell size normal distribution by decreasing the frequency of any given area. This in turn can make automated calculation of cell area more difficult or impossible because a definite minimum area cannot be determined. Similarly, this also applies to images with few numbers of cells where the frequency of different cell sizes may be very low (<50). As a result, when analyzing images with different resolutions than those used in this study or different distributions of cell size, manual identification of cell size range may be needed (10.3.3).
Significance and Future Directions
The ImageJ plugin Cell Counter was initially released in 2001 by Dr. Kurt De Vos and has served as a staple for manual cell counting to this day. To continue the trend of freely available tools for the biological community, Cell Concentration Calculator and migration assay counter offer the next step in free tools to help increase throughput and inter-experimental reproducibility of migration assays. The end result of migration assays is highly dependent on the number of cells seeded. Ergo a reliable cell count is a necessity for reproducibility. In this way, CCC offers both increased accuracy due to higher statistical certainty of cell concentration and shorter count times preserving cell health.
Moreover, the end result of both plugins is a count of cell number and not a relative index such as fluorescence. Various other protocols exist that measure fluorescence after staining but these methods suffer from lack of sensitivity13. Adobe's Photoshop offers its own particle analysis tool but the program must be purchased for use and does not offer the post-counting analysis available from migration assay counter20. Both plugins rely on the particle counting feature from ImageJ and consequently can be modified by the end user by editing the macros used by ImageJ. This offers greater flexibility for the user to expand the scope of the plugins to other particles by creating new macros. Further development to increase the breadth of countable particles and incorporate end user contributions is the next logical step. By combining the core principles of Cell Counter with an automated counting algorithm and post-counting analysis, this greatly increases the ease with which migration assays can be processed without any loss of accuracy. The plugins are freely available for download with installation instructions from: http://peng.lab.yorku.ca/imagej-plugins.
The authors have nothing to disclose.
This work was supported by the Canadian Institute of Health Research to CP (OR 142730 and OR 89931). We would like to thank Jelena Brkic for her initial idea of binary particle analysis in ImageJ.
HyClone Classical Liquid Media: RPMI 1640 – With L-Glutamine |
Fisher Scientific | SH3002702 | Cell culturing media |
Fetal bovian serum (FBS) | GIBCO BRL | P00015 | Media suppliment |
HTR8/SVneo trophoblast cell line | Cells were obtained from Dr. Charles Graham (Queen’s University, Kingston, Canada) | Software is designed to work with any cell line. | |
Trypsin | GIBCO BRL | 27250-018 | Prepared as 0.20% (w/v) in 10uM EDTA 1X PBS |
Accutase | Innovative Cell Technologies | AT104 | |
10 cm cell culture plates | SARSTEDT | 833902 | Any tissue culture treated plates will be suitable |
Transwell Polyester Membrane Inserts – 8.0µm Pore size | Costar 3422 ordered from Fisher Scientific | 7200150 | For 24-well plates; Pore size: 8.0µm; 6.5mm diameter; 0.33cm2 growth area |
HARLECO Hematology Stains and Reagents, EMD Millipore – Soluntions 1, 2 & 3 | EMD Millipore and ordered from VWR | 65044A, B, & C | Hemacolor stain set consists of three 500mL (16.9oz.) poly bottles & includes a methanol fixative (Solution 1), an eosin or acid stain (Solution 2), and a methylene blue or basic stain (Solution 3) |
Cotton Tipped Applicator | Puritan Medical | 806-WC | |
Single-edge industrial razor blades | VWR | 55411 – 055 | Thickness: 0.30 mm (0.012") |
Microscope Slides – Precleaned/Plain | Fisher Scientific | 12550A3 | Dimentions: 25 X 75 X 1.0 mm |
Fisherbrand Cover Glasses – Rectangles no. 1 | Fisher Scientific | 12-545E | Thickness: 0.13 to 0.17mm; Size: 50 x 22mm |
Fisher Chemical Permount Mounting Medium | Fisher Scientific | SP15-500 | |
Leica Stereo dissecting microscope | Leica Microsystems | The microsope is equipped with Leica microscope camera Model MC170 HD & camera software is Leica App. Suite (LAS E2) Version 3.1.1 [Build: 490]. Microscope parts: LED3000 Spot Light Illumination Model: MEB126, Leica M80 Optic Carrier Model M80, Objective achromat 1.0X, WD=90mm Model: MOB306 & Objective achromat 0.32X, WD=303mm Model: MOB315, Video Objective 0.5X Model: MTU-293, | |
Hemacytometer | Assistant Germany | 0.100 mm Depth – 0.0025 mm2 | |
Olympus inverted light microscope | Olympus Corporation | CKX41SF | The microsope is equipped with Lumenera Infinity 1-2 2.0 Megapixel CMOS Color Camera & camera software is Infinity analyze Version 6.5.2 |
Laminar flow cabinet 1300 Series A2 | Thermo Scientific | Model: 1375 | Any laminar flow cabinet for cell culture work will be suitable |
Cell culture incubator | Thermo Scientific | Model: 370 | Any cell culture incubator will be suitable – Cells were cultured under humidefied environment, 5% CO2, 37 °C |
ImageJ | NIH | Version 1.50e | Minimum version required |
Java Runtime Environment | Oracle | Version 1.8.0_66 | Minimum version required |