Skip to main content
  • Methodology article
  • Open access
  • Published:

Local cell metrics: a novel method for analysis of cell-cell interactions



The regulation of many cell functions is inherently linked to cell-cell contact interactions. However, effects of contact interactions among adherent cells can be difficult to detect with global summary statistics due to the localized nature and noise inherent to cell-cell interactions. The lack of informatics approaches specific for detecting cell-cell interactions is a limitation in the analysis of large sets of cell image data, including traditional and combinatorial or high-throughput studies. Here we introduce a novel histogram-based data analysis strategy, termed local cell metrics (LCMs), which addresses this shortcoming.


The new LCM method is demonstrated via a study of contact inhibition of proliferation of MC3T3-E1 osteoblasts. We describe how LCMs can be used to quantify the local environment of cells and how LCMs are decomposed mathematically into metrics specific to each cell type in a culture, e.g., differently-labelled cells in fluorescence imaging. Using this approach, a quantitative, probabilistic description of the contact inhibition effects in MC3T3-E1 cultures has been achieved. We also show how LCMs are related to the naïve Bayes model. Namely, LCMs are Bayes class-conditional probability functions, suggesting their use for data mining and classification.


LCMs are successful in robust detection of cell contact inhibition in situations where conventional global statistics fail to do so. The noise due to the random features of cell behavior was suppressed significantly as a result of the focus on local distances, providing sensitive detection of cell-cell contact effects. The methodology can be extended to any quantifiable feature that can be obtained from imaging of cell cultures or tissue samples, including optical, fluorescent, and confocal microscopy. This approach may prove useful in interpreting culture and histological data in fields where cell-cell interactions play a critical role in determining cell fate, e.g., cancer, developmental biology, and tissue regeneration.


Cell-cell recognition is critical to a wide range of problems in biology and medicine [116]. The development of experimental approaches associated with cell-cell recognition has promoted advances in understanding these effects, e.g., biochemical assays for protein binding and transcription,. However, less attention has been focused on developing algorithms for the detection of cell-cell recognition from the structure and spatial distribution of cells. Such methods would offer complimentary benefits to biochemical assays, due to the relative ease of collecting microscopy data from cell cultures and tissues. This would be useful also in combinatorial and high-throughput screening of cell-cell and cell-material interactions [1723]. In adhesion dependent cells, cell-cell recognition is known to be a crucial step in initiating contact inhibition (CI) of cell migration[24, 25] and proliferation[14]. CI plays an important role in the proliferation, invasion, and metastasis of cancer cells, [2628] cardiovascular tissue homeostasis and development, [5, 29] embryonic development, [1] and wound healing, [16] among many other biological phenomena. Conversely, it has also been shown that under certain conditions cell-cell contact can promote cell proliferation, known as "density-dependent" contact stimulation of cell proliferation[12, 3032]. The investigation of mechanisms relating cell contact, cell proliferation, migration, and differentiation, in which cell adhesion molecules play a major role, is a rich area of research. Cadherins, e.g., VE-cadherin, [3] E-cadherin[13] and N-cadherin, [8, 9, 15, 33] and notch proteins[1] have been shown to mediate contact-dependent phenomena in a wide variety of cell types.

Here, CI of proliferation, a known cell-cell recognition phenomenon, is used as a model system for developing algorithms for the analysis of cell-cell recognition from microscopy data. Usually, the effects of cell density on proliferation are studied as relationships between global descriptions, such as average cell density, proliferation rate, and protein expression level[2, 11, 12, 14]. However, as we demonstrate below, these summary-statistic descriptions are only sensitive to the effects of very large changes in cell density. As a result, global metrics do not illuminate all of the information available from image data for cell contact phenomena. This is because cell-cell contacts are local interactions and are very sensitive to short-range cell-to-cell distance. When global metrics are used, such as cell density in a set of images, all cell-to-cell distances are treated equally. Critical information pertaining to nearest-neighbor cell-to-cell distances is "diluted" by the many other cell-to-cell distances in the data set, which is observed as noise in the data. Furthermore, the stochastic characteristics of cell behavior add to the noise so that only major trends of the responses of cells to neighbors is distinguishable over very large changes in cell density[15]. To address this "dilution" effect, Nelson and Chen studied contact stimulation effects on the growth of a single pair of cells by using specially-designed surfaces to decouple the effects of cell-cell contact from others[12].

In this paper, we introduce a complementary approach that allows focused analysis on nearest-neighbor cells, but permits sampling from cultures with high cell densities and use of any type of surface. We had previously applied this technique to screen large image databases from cell cultures on combinatorial libraries of biomaterials[22, 23]. Here, we outline the details of this method and discuss it's generalization to Bayesian modelling. The method is based upon defining local cell metrics (LCMs), which are histograms of cell properties. The use of these local variables expands the sample space considerably and allows separation of arbitrarily-defined short- and long-range effects. We show how the local cell metrics are then incorporated into a Bayesian model. The new method and model are examined quantitatively and compared with traditional summary approaches in a study of contact inhibition of osteoblast proliferation.



Surface Preparation

Poly (DL-lactic-glycolic acid) (PLGA, block copolymer, 50:50 ratio of PGA and PLA, 40,000~75,000 Da) and poly (ε-caprolactone) (PCL, 114,000 Da, Mw/Mn = 1.43) were obtained from Sigma Aldrich, St Louis, MO. PLGA and PCL, respectively, were dissolved in chloroform to 8% and 5% mass and spin coated on silicon chips (22 × 22 mm). To provide adhesion of these polymers to the silicon during cell culture, the silicon was pretreated with a Piranha etch (70% H2SO4/21% H2O/9% H2O2 at 80°C for 1 h) followed by 1 min in a HF acid bath and a final rinse in DI water (filtered at 0.2 μm).

Cell Culture

Established from newborn mouse calvaria, [34] the MC3T3-E1 cell line has been shown capable of differentiating into osteoblast and osteocytes in vitro[35]. MC3T3-E1 cells have been shown to exhibit specific bone related protein expression patterns, under different developmental stages, similar to primary mouse calvaria cells[35, 36]. This cell line is thus a suitable in vitro model for investigating cell behaviors, regulations of such behaviors, and underlying mechanisms in different osteoblast maturation stages[37]. Since the original MC3T3-E1 cell line has been found phenotypically heterogeneous with regard to cell differentiation, more homogeneous subclones have been established[38]. In this work, MC3T3-E1 subclone 4 (from ATCC, VA), which shows homogenous capabilities of osteogenesis both in vitro and in vivo, [38] was chosen in order to minimize variations due to phenotypical heterogeneities.

Cell proliferation was assayed by BrdU immunohistochemistry. Briefly, PLGA- and PCL-coated wafers were mounted into Costar® 6-Well TC-Treated Microplates (Corning, NY). The tissue culture treated polystyrene (TCPS) surfaces of the microplate wells were used as controls. After sterilization (70% ethanol solution, 30 min), MC3T3-E1 cells (passage 6) were seeded onto the coated wafers at 4 × 104 cells/cm2. This relatively high seeding density was selected to highlight the effects of contact inhibition of cell growth and other space-sensitive cell-to-cell interactions. After seeding, microplates were shaken for 10 min on a shaker (Instrument model, operation frequency) to obtain uniform seeding. Cells were cultured in DMEM (Cellgro® Dulbecco's Modification of Eagle's Medium, Mediatech, Inc., VA) with 10% fetal bovine serum (ATCC® SCRC-1002™, ATCC, VA), L-glutamine and streptomycin at 37°C in a humidified 5% CO2 atmosphere. At 5 h post seeding, surfaces were washed with Dulbecco's Phosphate-Buffered Saline (DPBS, with Ca++ and Mg++) to remove non-attached cells, and fresh culture medium was then added. At 18 h post seeding, 2 mM BrdU (5-bromo-2-φ-deoxyuridine, Sigma, MO) in PBS was added to the culture medium to reach a final concentration of 20 μM. After 6 h of BrdU incorporation, cells were fixed with 3.6% paraformaldehyde and BrdU incorporation was assayed by immunohistochemistry (primary antibody: mouse anti-BrdU, BD Biosciences, CA; secondary antibody: goat anti-mouse, Rhodamine conjugated, Rockland Immunochemicals, Inc., PA; counter staining: Hoechst 33342, Molecular Probes, Invitrogen Corporation, CA).

Low calcium concentration suppresses contact inhibition of cell growth by deactivating calcium-dependent cadherins[7, 39]. This phenomenon was used in this study to validate the local cell metrics, and at the same time the dependency of contact inhibition on calcium was quantitatively studied. In order to investigate the role of Ca++ on cell spreading and proliferation, BrdU incorporation experiments in low Ca++ medium were performed on TCPS surfaces. Fifteen minutes before the introduction of BrdU, cells were rinsed twice with DPBS (without Ca++ and Mg++) and afterword cultured in the low Ca++ medium (0.5% FBS in Ca++ and Mg++ free DPBS)[39]. The rest of the protocol was the same as previously described.

Image Acquisition

Cell locations and proliferation were quantified using fluorescent microscopy (Olympus BX51 Clinical Microscope). A robotic translation stage was used to image predetermined locations on each culture surface using a MicroFire™ monochromic digital camera (SKU S99826, Optronics, CA). The image locations were fixed on a 16 × 20 grid with horizontal and vertical spacing of 1280 μm and vertical spacing of 960 μm. For each location a 1189 × 892 μm2 BrdU staining image and Hoechst counter staining image were acquired at a resolution of 1600 × 1200 pixels2. All images and contextual information were organized and stored in an Oracle® 10 g (Oracle, CA) database for further image processing and data analysis.

Image Processing

The Image Processing Toolbox of Matlab™ R14 (MathWorks, MA) was employed for image processing. Due to the volume of image data dynamic, self-adapting algorithms were developed for automated image processing. Binary images of both surface lateral patterns of cell nuclei counter staining were obtained from raw grayscale microscopic images by a variation-adjusted iterative selection method (VAIS), which was modified from the original iterative selection method [4045] (details in the BrdU thresholding part below).

Binary images of cell nuclei were segmented by the marker-controlled watershed method[42] to separate images of closely-spaced cell nuclei. This process was critical because the nuclei of a pair of recently-proliferated cells were often too close to be distinguished with thresholding alone. Resultant black-and-white cell nuclei images were used as masks by overlaying them with corresponding BrdU staining images to determine the fluorescence intensity of incorporated BrdU. The histogram of BrdU staining intensity per nucleus [see Additional File 1: Figure S.1] is composed of two major peaks: the low intensity peak (background) represents cells at rest, while the high intensity peak (foreground) indicates proliferating cells. The optimal threshold between these two peaks was determined automatically by VAIS. Briefly, starting at an initial threshold T i = 0.5, the histogram was divided into resting (background) and proliferating (foreground) parts. Means and standard deviations of the foreground and background, respectively denoted as M bi , M fi , σ bi , and σ fi , were determined by fitting each peak to a Gaussian curve. A new threshold was calculated as Ti+1= (σ fi M fi + σ bi M bi )/(σ fi + σ bi ) and was repeated until convergence on a stable threshold. Compared with more common iterative selection methods, which use a simple mean intensity, the modified VAIS procedure is more robust when background and foreground intensities have different variances. Indeed, the variance of the BrdU signal intensity from non-proliferating cells was significantly greater than that of the proliferating cells [see Additional File 1: Figure S.1]. During image processing, data washing in the form of median filtering was performed to remove noise below a certain threshold. Image processing was supervised in order to assure the performance of self-adaptive algorithms and images of poor quality not permitting quantification were occasionally discarded. Proliferation behaviors were determined for every cell and stored in the database along with the cell location on the surfaces.

Methodologies of Data Analysis

Global Metrics

Cell density and proliferation were described with summary statistics such as number of resting and proliferated cells computed for each image. This provides a set of global metrics for features in each image. As indicated in Figure 1a, global metrics are most naturally understood in terms of conventional summary-statistics, exploratory data analysis, and well-known methods for estimating confidence and significance levels based on an assumed probability distribution. The ability to detect contact inhibition of cell proliferation, a known phenomenon, was used as an indicator of the effectiveness of the global metrics cell density and proliferation averages.

Figure 1
figure 1

Schematic comparing global versus local metrics. Schematic comparing global versus local metrics. (a) global point of view leading to summary statistics, illustrated by bar graph, (b) local or individual-cell point of view leading to histograms, illustrated by frequency plot.

Definition of Local Cell Metrics

Source codes that implement the algorithms presented in this section have been made available by the authors. [see Additional File 2] Consider that the collection of all cells (A) is composed of either proliferated (P) or resting cells (R), such that A = P + R. The symbol A represents any cell chosen at random, regardless of proliferative status. The proliferating-resting cell distance, PR, is used here to illustrate the definition and properties of local cell metrics, as indicated in Figure 1b. The definitions below are generalizable to any type of cell-cell distance, or any other spatial or temporal metric of cells. Assume that in the kth image the number of P-class and R-class cells is n Pk and n Rk , the distance PR ijk between the centroids of the nuclei of the ith P-cell and the jth R-cell can be calculated readily from the results of image analysis. In the kth image, the set of all such distances, PR k is defined as

And for all images an overall set PR can be defined as

A set of N+1 distance bins is defined as

where d0, d1, ..., d N is a user-defined distance scale over which analysis is to be performed. The centroid of each interval in bin dist is defined as

and the resultant centroid set for bin dist is

The bin dist is used to sort set PR into an N-bin histogram

where is the number of PR distances that fall in the interval [di-1, d i ), which is centered at .

The total number of elements in set PR is

After normalizing by n PR , the frequency function LCM is

and represents the R-type cell environment of the P-type cells observed over distances . Frequency functions, denoted as , , , and , may also calculated for cell-to-cell distances PP, AA, RR, and PA in similar manner.

LCM Normalization

Normalization is necessary to interpret LCMs in a meaningful manner and to compare the probability of cell responses under different cell environments. One method of normalization is to relate observed occurrences to random occurrences. Given the finite image size and generally non-overlapping nature of cultured cells, the distribution of random cell occurrences is not Gaussian. The random distribution for cell-cell distance, f std , was calculated as the any cell-any cell distribution (f AA ) of 1× 1010 randomly-chosen nuclei positions on a simulated image 1600 pixels by 1200 pixels. The normalized LCM is

Other LCMs (f AA , f PA , f RR ) are normalized similarly, which allows direct comparisons of different types of cell distances on different surfaces.

In addition to normalizing by the standard distribution, f std , direct ratios between LCMs are used also in our analysis, in which case f std cancels, as indicated in the next equation.

The ratio rPR|PAhighlights the specific effects of non-proliferated cells on the central proliferating cell relative to the effects of any given cell. Thus, the probability of cell responses under different cell environments can be compared meaningfully. Furthermore, each set of cell-to-cell distances can be decomposed into subsets, which allows investigation of the contribution of each subset to the overall effect. Therefore, ratios of cell backgrounds may be constructed and used as classifiers for screening and identifying significant cell environment patterns. These ratios also define posterior odds (PO) of observing certain proliferation behaviors. For example, consider , and using the subscript i to signify the evaluation at a certain distance , the ratio is calculated as


Applying equations (12) and (13) to equation (11),

Thus, is a posterior odds that quantifies how the probability of cell proliferation is changed by the presence of a second cell located at distance , relative to the average proliferation for all cell-cell distances. Computationally, to promote the efficiency of the codes, we defined A k = {P k , R k }={P1k, P2k, ..., P NPk , R1k, R2k, ..., R NRk }, and removed self-to-self cell pairs (PA ijk where i = j) and identical cell pairs (PA ijk where i >j) from PA.

LCM Decomposition

Furthermore, each set of cell-to-cell distances can be decomposed into subsets, which allows isolation of each subset's contribution. For example, consider rPR|PAdefined above. As described graphically in Figure 2, since the denominator PA = PP PR (the union of distance sets PP and PR) and PPPR = ( = the empty set), one may remove the PR component from PA, and the denominator becomes PP. The result is that the ratio rPR|PAis transformed into rPR|PP. By removing the shared, or overlapping, component PR from the denominator, rPR|PPhas higher "contrast" for observing effects of R-cells on P-cells than rPR|PA.

Figure 2
figure 2

Schematic indicating the decomposition of local metrics. Schematic indicating the decomposition of local metrics into groups of cell-cell distances based on cell proliferative status. A = any cell, P = proliferated cell, R = resting cell (nonproliferated).

Connection of LCMs to Bayesian Statistics

Local cell metrics are naturally connected to Bayesian analysis, which is a powerful statistical method used for classification[46, 47]. Specifically for the PR distance, the Bayesian approach allows one to quantify the local environment of P cells, as the conditional probability of finding an R cell a certain distance PR from a P cell. Based on the definition of f PR in equation (8), a naïve Bayes model can be established as follows. Consider a "test" cell chosen at random. It is desired to predict the possibility this cell will be in proliferating status, based upon the local environment of non-proliferating cells, which is given by the following conditional probability function

where represents the probability of finding a non-proliferating cell at a distance of from the central, randomly chosen, test cell. Using Bayes's theorem,

In the above function, the components p() and p(prolif) are constants that can be determined from their frequency in the data. The only non-constant component, the class-conditional probability, is given by

Assuming the occurrence probabilities around the non-proliferating cell distances are conditionally independent (uncorrelated), then where ij (the naïve Bayes assumption). Hence, equation (17) reduces to

A key development is to notice that , which means that under the naïve Bayes assumption the LCM is in fact a class-conditional probability function. The term p() represents the probability of locating the R cells at distances () from any cell, which is . Hence, the Bayes conditional probability function from equation (16) becomes

The naïve Bayes model allows prediction of the probability of proliferation as a function of the LCMs, and , which are easily computed from a training data set, as is p(prolif). The evaluation of this modeling approach will be the subject of forthcoming work.

Results and Discussion

Traditional metrics

To provide a benchmark for establishing the effectiveness of local metrics, contact inhibition of cell proliferation was studied using global analysis first. For each image in the database, the overall cell proliferation is plotted versus cell density, shown in Figure 3 and Table 1. Although a trend towards lower proliferation at higher cell density is somewhat apparent in Figure 3, global analysis, does not allow for quantitative detection of CI effects on cell proliferation at a statistically-significant level. The linear regression (Table 1) yields in an adjusted R2 of 0.128 (on PLGA) or 0.109 (on PCL), indicating that the contact inhibition effect masked by "noise" in the data. Furthermore, it is obvious that in Figure 3, no regression function can be fit satisfactorily to the global statistics, since the noise is too high relative to the CI effects. The use of a larger range may allow the global analysis to distinguish contact effects from natural variance in cell properties. This has been demonstrated in repeated experiments on PLGA, PCL, and TCPS surfaces. [see Additional File 1: Figures S.2, S.3, S.4]. However, there are drawbacks to the use of larger ranges, such as the introduction of seeding-density effects that mask or alter the cell-cell interactions.

Figure 3
figure 3

Effects of global cell density on global cell proliferation. Effects of global cell density on global cell proliferation ratio for MC3T3-E1 cultured on (a) PLGA and (b) PCL. Number of images used was (a) 353 and (b) 288.

Table 1 Linear Regression from Global Analysis Results from Figure 3

Local Cell Based Metrics

The noise level inherent to proliferation measurements, which are normally carried out over a small seeding density range, make contact inhibition a robust test-case for comparing local vs. global metrics. A contact phenomenon is detected when a relevant metric changes significantly relative to the data sampling noise. For global statistics, the assumed distribution (usually Normal) provides the random noise reference. For local metrics, the random cell-cell distance frequency distribution was calculated using a Monte-Carlo approach, termed the standard frequency distribution, f std . The reference f std is shown in Figure 4 together with the experimental f AA for MC3T3-E1 osteoblasts on PLGA. The profile of f std is similar to a beta- or chi-distribution with asymmetry due to the non-overlapping nature of the nuclei centers at close distances. The computed f std distribution is nearly identical to the experimental f AA distribution at large distances (> 100 μm). This is expected since f AA indicates the likelihood of finding any two cells (whether proliferating or not) separated by a given distance, which should in principle be random. Figure 4 also shows the distance distribution f PA , which is the likelihood of finding a proliferated cell a certain distance from any cell. If cell-cell distance has any relation to proliferative status then f PA and f AA should differ from one another and from f std , but only at close distances where cell-cell contact is likely to occur.

Figure 4
figure 4

Comparison of experimental and calculated distributions. Comparison of experimentally determined f AA and f PA for MC3T3-E1 on a PLGA surface and the computed standard curve, f std . These data represent the frequency at which any cell is located a certain distance from any other cell (f AA ) or from a proliferated cell (f PA ) in the experiments. The standard curve was computed from a Monte Carlo simulation and represents a uniformly random probability of locating any cell a given distance from any other cell on the same sized area as our microscope images. Number of images used was 353.

Figures 5, 6, and 7, which show the normalized distributions f PA /f std , f AA /f std and f PP /f std at close distances, indicate the non-random effects of contact inhibition when the values become less than one. Specifically, CI occurs when the distance between cell nuclei becomes less than about 50 μm. The typical mean cell area was around 2500 μm2, resulting in a mean diameter of 56 μm, which corresponds closely to the onset of CI. Representative images of cultured MC3T3-E1 cells on these surfaces have been presented in previous work[22, 23]. In addition, local fine structure in the contact inhibition region is observed as a local maximum peak between 10 and 20 μm. This peak indicates enhanced local proliferation at very close distance, even when overall proliferation is being inhibited. Interestingly the local peak magnitudes at 10 to 20 μm always follow the order f PP >f AA >f PA on each of the three surfaces examined, TCPS, PLGA and PCL. We hypothesize that the local enhancement peak is due to two daughter cells (from the same parent cell) that are very close, which have not had enough time to migrate away during the BrdU staining time period. If so, then this cell division peak should appear on the f PP curve but not the f PR curve, which was observed comparing Figures 5, 6, and 7 (f PP ) to Figures 8, 9, 10, and 11 (f PR ). In addition, in the Monte Carlo simulation of random cell positions (f std ), with no proliferation, this local peak is absent.)

Figure 5
figure 5

Normalized distributions for cells on PLGA surfaces. Normalized distributions for any cell-any cell (AA), proliferating-any cell (PA) and proliferating-proliferating cell (PP) on a PLGA surface: , , and . Normalization was performed by dividing the experimental frequency distribution by f std , the random cell-cell distribution determined from Monte Carlo simulation. Number of images used was 353.

Figure 6
figure 6

Normalized distributions for cells on PCL surfaces. Normalized distributions as described in Figure 5 on a PCL surface: , , and . Number of images used was 288.

Figure 7
figure 7

Normalized distributions for cells on TCPS surfaces. Normalized distributions as described in Figure 5 on a TCPS surface: , , and . Number of images used was 291.

Figure 8
figure 8

Ratio of distributions PA/RR on PLGA surfaces. Ratio of distributions on a PLGA surface. Number of images used was 353.

Figure 9
figure 9

Ratio of distributions PR/RR on PLGA surfaces. Ratio of distributions on a PLGA Surface. Number of images used was 353.

Figure 10
figure 10

Ratio of distributions PR/RR on PCL surfaces. Ratio of distributions rPR|RRon a PCL surface. Number of images used was 288.

Figure 11
figure 11

Ratio of distributions PR/RR on TCPS surfaces. Ratio of distributions rPR|RRon a TCPS surface. Number of images used was 291.

Direct ratios between experimentally-determined distributions can be chosen specifically to illuminate the CI phenomena of interest. Specifically, common components in the numerator and denominator not related to CI phenomena may be removed, thus isolating the phenomena of interest. This process is examined in Figures 8, 9, 10, and 11. Figure 8 shows the rPA|AAprofile, which is classified into two regions: the proliferation suppressed region (0 to 40 μm) where cell proliferation was suppressed up to 6 fold (rPA|AAfalls to ~1/6) and the null region (beyond 40 μm) where cell proliferation was not noticeably affected by the contact of other cells. Based on the discussion for Figure 7, the ratio rPA|AAdoes not fully decouple the division and daughter-cell migration phenomena (indicated by PP) from the proliferation phenomena (indicated by PR), since PP PA. We illustrate how separation of the PP and PR components enhances the detection of contact inhibition of proliferation. By definition the various distances are related as follows


where is the empty set. The two shared components of AA and PA are PP and PR. The PP distance component represents distances between cells that have both proliferated, i.e., proliferation in those pairs was not contact inhibited. Removal of the common PP component from the numerator and denominator of rPA|AAleads to rPR|RR, shown in Figures 9, 10, and 11 for the PLGA, PCL, and TCPS surfaces, respectively. By definition rPR|RRshould be more sensitive to CI of proliferation, because non contact inhibited cell pairs have been removed. In Figure 9, the rPR|RRratio is classified into two regions: the NaN region (below 5 μm) where few pairs occur, and the contact inhibition region (5 to 40 μm). In the contact inhibition region, a clear trend of decreasing probability of finding a neighboring cell is seen as the distance between cells decreases. A minimum is observed at d min = 8 μm, where contact inhibition effects are maximized. To our knowledge, this is the first time both the magnitude and the range of contact inhibition of cell proliferation have been determined quantitatively in a single function.

The physical meaning of the LCM ratio rPR|RRcan be can be illustrated by recognizing that it is the posterior odds (PO) of proliferation as a function of cell-cell distance. Consider two cells that are well-separated at 40 μm, and another two cells that are at a close distance of 8 μm, where the extreme in contact inhibition behavior is found (minimum rPR|RRin Figure 9). The PO that one of the closely-spaced cells has proliferated is POPR/RR= rPR|RR(8)/rPR|RR(40) = 1/32. This means there is a 32 fold lower chance of proliferation at 8 μm than at a distance of 40 μm.

The profiles of rPR/RRfrom the other polymer surfaces are shown in Figure 10 (PCL surface) and Figure 11 (TCPS surface). The ratios are similar in shape but have different magnitudes for the minimum point as a function of the surface. Table 2 summarizes the variation of rPR|RR, minand d min on the different surfaces. The different location and strength of contact inhibition might be due to surface features such as roughness, crystallinity, hydrophobicity, surface charge, or protein adsorption, factors which are known to influence osteoblast proliferation[48]. For example, the surface roughness increases in the order TCPS < PLGA < PCL, and at the same time the POPP/PRis decreasing, and the d min is increasing. We illustrate this point, however, not to make a definite mechanistic argument about surface effects on proliferation, which is certainly more complicated than roughness alone. Rather, the point made is that the LCM method is capable of sensitive detection of differences in proliferation for cells cultured on different surfaces.

Table 2 Minima in rPR|RRCurve Indicating Maximum Contact Inhibition

The effect of calcium depletion on LCMs is presented in Figure 12. The distribution of the control (calcium +) was significantly lower than the low calcium (calcium -) case. In addition, the 'calcium -' curve stays close to unity except at very close distances, less than 20 μm, whereas the 'calcium +' curve falls below unity at 50 μm. Hence, the LCM detects the expected result: that low calcium should inhibit the cell-cell self-avoidance and contact inhibition[7, 39]. This effect is seen more clearly in examining the ratio rPR/RR, in Figure 13. Contact inhibition was very strong within a cell-to-cell distance of 30 μm when cells were cultured under physiological calcium concentration (calcium +). However, contact inhibition disappeared when calcium was depleted (calcium -).

Figure 12
figure 12

Effects of calcium depletion: loss of contact inhibition. Normalized distribution examining effects of calcium depletion on proliferation of MC3T3-E1 cultured on TCPS surfaces. Number of images used was 198.

Figure 13
figure 13

Effects of calcium depletion: loss of contact inhibition. Ratio of distributions rPR|RRexamining effects of calcium depletion on proliferation of MC3T3-E1 cultured on TCPS surfaces. Number of images used was 204.


We have shown that global summary statistics are not adequate metrics for detecting local cell interactions, due to noise and non-local effects inherent to cell-cell contact phenomena. A novel data analysis strategy, local cell metrics, has been introduced in this paper. LCMs, which are cell-cell distance histograms, describe cell environments from the "point of view" of individual cells. These metrics allow focusing of analysis onto arbitrarily-defined close distances. In addition, LCMs can be 'tuned' to be sensitive to specific contact phenomena by decomposing the distributions into specific cell-types (proliferating vs. nonproliferating) and removing unwanted components. Local metrics as defined herein are not limited to proliferation analysis, nor to cell-cell interactions alone. The metrics are generic and can be, in principle, applied to any type of quantifiable cell assay, and can be applied to cell-biomaterial and cell-tissue interactions as well. We have also shown how LCMs are related to the naïve Bayes model, which makes them useful for data mining and classification (the subject of forthcoming work.)

We have demonstrated the new local metrics by considering the contact inhibition of proliferation of the osteoblast cell line MC3T3-E1. A quantitative and probabilistic description of the contact inhibition effect as a function of cell-cell distance has been achieved. In fact, the probability of proliferation is shown to be strongly dependent on the distance to, and proliferative state of, neighboring cells. The LCMs were also sensitive to effects of the culture surface, and of calcium composition in the culture media, on proliferation.


  1. Artavanis-Tsakonas S, Matsuno K, Fortini M: Notch signaling. Science 1995, 5208: 225–232. 10.1126/science.7716513

    Article  Google Scholar 

  2. Castor LN: Contact Inhibitions of Cell Division and Cell Movement. J Investig Dermatol 1972, 59(1):27–32. 10.1111/1523-1747.ep12625732

    Article  CAS  PubMed  Google Scholar 

  3. Caveda L, Martin-Padura I, Navarro P, Breviario F, Corada M, Gulino D, Lampugnani MG, Dejana E: Inhibition of Cultured Cell Growth by Vascular Endothelial Cadherin (Cadherin-5/VE-Cadherin). Journal of Clinical Investigation 1996, 98(4):886–893. 10.1172/JCI118870

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Charrasse S, Comunale F, Gilbert E, Delattre O, Gauthier-Rouviere C: Variation in cadherins and catenins expression is linked to both proliferation and transformation of Rhabdomyosarcoma. Oncogene 2004, 23(13):2420–2430. 10.1038/sj.onc.1207382

    Article  CAS  PubMed  Google Scholar 

  5. Grego-Bessa L, Luna-Zurita G, del Monte V, Bolós P, Melgar A, Arandilla A, Garratt H, Zang Y, Mukouyama HC: Notch Signaling Is Essential for Ventricular Chamber Development. Developmental Cell 2007, 12: 415–429. 10.1016/j.devcel.2006.12.011

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Hirano S, Nose A, Hatta K, Kawakami A, Takeichi M: Calcium-dependent cell-cell adhesion molecules (cadherins): subclass specificities and possible involvement of actin bundles. J Cell Biol 1987, 105(6):2501–2510. 10.1083/jcb.105.6.2501

    Article  CAS  PubMed  Google Scholar 

  7. Kandikonda S, Oda D, Niederman R, Sorkin BC: Cadherin-Mediated Adhesion Is Required for Normal Growth Regulation of Human Gingival Epithelial Cells. Cell Communication & Adhesion 1996, 4(1):13–24.

    Article  CAS  Google Scholar 

  8. Levenberg S, Katz BZ, Yamada KM, Geiger B: Long-range and selective autoregulation of cell-cell or cell-matrix adhesions by cadherin or integrin ligands. Journal of Cell Science 1998, 111(3):347–357.

    CAS  PubMed  Google Scholar 

  9. Levenberg S, Sadot E, Goichberg P, Geiger B: Cadherin-mediated transmembrane interactions. Cell Adhes Commun 1998, 6(2–3):161–170. 10.3109/15419069809004473

    Article  CAS  PubMed  Google Scholar 

  10. Mary S, Charrasse S, Meriane M, Comunale F, Travo P, Blangy A, Gauthier-Rouviere C: Biogenesis of N-Cadherin-dependent Cell-Cell Contacts in Living Fibroblasts Is a Microtubule-dependent Kinesin-driven Mechanism. Molecular Biology of the Cell 2002, 13(1):285–301. 10.1091/mbc.01-07-0337

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Nakatsuji Y, Miller RH: Density dependent modulation of cell cycle protein expression in astrocytes. Journal of Neuroscience Research 2001, 66(3):487–496. 10.1002/jnr.1240

    Article  CAS  PubMed  Google Scholar 

  12. Nelson CM, Chen CS: Cell-cell signaling by direct contact increases cell proliferation via a PI3K-dependent signal. FEBS Letters 2002, 514: 238–242. 10.1016/S0014-5793(02)02370-0

    Article  CAS  PubMed  Google Scholar 

  13. St Croix B, Sheehan C, Rak JW, Florenes VA, Slingerland JM, Kerbel RS: E-Cadherin-dependent Growth Suppression is Mediated by the Cyclin-dependent Kinase Inhibitor p27KIP1. J Cell Biol 1998, 142(2):557–571. 10.1083/jcb.142.2.557

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Stoker MGP, Rubin H: Density Dependent Inhibition of Cell Growth in Culture. Nature 1967, 215(5097):171–172. 10.1038/215171a0

    Article  CAS  PubMed  Google Scholar 

  15. Warchol ME: Cell Density and N-Cadherin Interactions Regulate Cell Proliferation in the Sensory Epithelia of the Inner Ear. Journal of Neuroscience 2002, 22(7):2607–2616.

    CAS  PubMed  Google Scholar 

  16. Zegers MMP, Forget M-A, Chernoff J, Mostov KE, Beest MBAt, Hansen SH: Pak1 and PIX regulate contact inhibition during epithelial wound healing. The EMBO Journal 2003, 22: 4155–4165. 10.1093/emboj/cdg398

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Broderick S, Suh C, Nowers J, Vogel B, Mallapragada S, Narasimhan B, Rajan K: Informatics for combinatorial materials science. Journal of the Minerals, Metals and Materials Society 2008, 60: 56–59.

    Article  CAS  Google Scholar 

  18. Kennedy S, Washburn N, Simon C, Amis E: Combinatorial screen of the effect of surface energy on fibronectin-mediated osteoblast adhesion, spreading and proliferation. Biomaterials 2006, 27: 3817–3824. 10.1016/j.biomaterials.2006.02.044

    Article  CAS  PubMed  Google Scholar 

  19. Kohn J: New approaches to biomaterials design. Nature Materials 2004, 3: 745–747. 10.1038/nmat1249

    Article  CAS  PubMed  Google Scholar 

  20. Meredith JC: Advances in Combinatorial and High-Throughput Screening of Biofunctional Polymers for Gene Delivery, Tissue Engineering and Anti-Fouling Coatings. Journal of Materials Chemistry 2009, 19: 34–45. 10.1039/b808649d

    Article  CAS  Google Scholar 

  21. Meredith JC, Sormana JL, Keselowsky BG, Garcia AJ, Tona A, Karim A, Amis EJ: Combinatorial characterization of cell interactions with polymer surfaces. J Biomed Mater Res A 2003, 66(3):483–490. 10.1002/jbm.a.10004

    Article  PubMed  Google Scholar 

  22. Su J, Meredith JC: Local Histogram Analysis: Detecting Cell-Microstructure Interactions on Combinatorial Biomaterial Libraries. Combinatorial Chemistry and High-Throughput Screening 2009, 12: 626–633. 10.2174/138620709788681880

    Article  CAS  PubMed  Google Scholar 

  23. Su J, Zapata P, Meredith JC: Knowledge Discovery Applications in High-Throughput Polymer Characterization. In Combinatorial Methods and Informatics in Materials Science. Edited by: Fasolka M. Warrendale, PA: Materials Research Society; 2006. 0894-LL0805–0804 0894-LL0805-0804

    Google Scholar 

  24. Abercrombie M, Ambrose EJ: Interference microscope studies of cell contacts in tissue culture. Experimental Cell Research 1958, 15(2):332–345. 10.1016/0014-4827(58)90034-X

    Article  CAS  PubMed  Google Scholar 

  25. Abercrombie M, Heaysman JEM: Observations on the social behaviour of cells in tissue culture: II. "Monolayering" of fibroblasts. Experimental Cell Research 1954, 6(2):293–306. 10.1016/0014-4827(54)90176-7

    Article  CAS  PubMed  Google Scholar 

  26. Abercrombie M, Heaysman JEM, Karthauser HM: Social behaviour of cells in tissue culture III. Mutual influence of sarcoma cells and fibroblasts. Experimental Cell Research 1957., 13(2): 10.1016/0014-4827(57)90007-1

  27. Bremnes RM, Veve R, Hirsch FR, Franklin WA: The E-cadherin cell-cell adhesion complex and lung cancer invasion, metastasis, and prognosis. Lung Cancer 2002, 36(2):115–124. 10.1016/S0169-5002(01)00471-8

    Article  PubMed  Google Scholar 

  28. Sharma VM, Draheim KM, Kelliher MA: The Notch1/c-Myc pathway in T cell leukemia. Cell Cycle 2007, 6: 927–930.

    Article  CAS  PubMed  Google Scholar 

  29. Liu Z-J, Shirakawa T, L Y, Soma A, Oka M, Dotto GP, Fairman RM, Velazquez OC, Herlyn M: Regulation of Notch1 and Dll4 by Vascular Endothelial Growth Factor in Arterial Endothelial Cells: Implications for Modulating Arteriogenesis and Angiogenesis. Mol Cell Biol 2003, 23: 14–25. 10.1128/MCB.23.1.14-25.2003

    Article  PubMed Central  PubMed  Google Scholar 

  30. Burke JM: Cell-cell contact promotes DNA synthesis in retinal glia but not in fibroblasts. Exp Cell Res 1983, 146(1):204–206. 10.1016/0014-4827(83)90339-7

    Article  CAS  PubMed  Google Scholar 

  31. Burke JM: Growth in retinal glial cells in vitro is affected differentially by two types of cell contact-mediated interactions. Exp Cell Res 1989, 180(1):13–19. 10.1016/0014-4827(89)90207-3

    Article  CAS  PubMed  Google Scholar 

  32. Gaunt SJ, Subak-Sharpe JH: Cell cycle variation associated with feeder effects in cultures of Chinese hamster fibroblasts. Exp Cell Res 1977, 109(2):341–348. 10.1016/0014-4827(77)90013-1

    Article  CAS  PubMed  Google Scholar 

  33. Polyak K, Kato JY, Solomon MJ, Sherr CJ, Massague J, Roberts JM, Koff A: p27Kip1, a cyclin-Cdk inhibitor, links transforming growth factor-beta and contact inhibition to cell cycle arrest. Genes Dev 1994, 8(1):9–22. 10.1101/gad.8.1.9

    Article  CAS  PubMed  Google Scholar 

  34. Kodama H, Amagai Y, Sudo H, Kasai S, Yamamoto S: Establishment of a clonal osteogenic cell line from newborn mouse calvaria. Japanese Journal of Oral Biology 1981, 23: 899–901.

    Article  Google Scholar 

  35. Sudo H, Kodama HA, Amagai Y, Yamamoto S, Kasai S: In vitro differentiation and calcification in a new clonal osteogenic cell line derived from newborn mouse calvaria. The Journal of Cell Biology 1983, 96(1):191–198. 10.1083/jcb.96.1.191

    Article  CAS  PubMed  Google Scholar 

  36. Choi JY, Lee BH, Song KB, Park RW, Kim IS, Sohn KY, Jo JS, Ryoo HM: Expression patterns of bone-related proteins during osteoblastic differentiation in MC3T3-E1 cells. Journal Of Cellular Biochemistry 1996, 61(4):609–618. 10.1002/(SICI)1097-4644(19960616)61:4<609::AID-JCB15>3.0.CO;2-A

    Article  CAS  PubMed  Google Scholar 

  37. Sodek J, Berkman FA: Bone cell cultures. Methods Enzymol 1987, 145: 303–324. full_text

    Article  CAS  PubMed  Google Scholar 

  38. Wang D, Christensen K, Chawla K, Xiao G, Krebsbach PH, Franceschi RT: Isolation and characterization of MC3T3-E1 preosteoblast subclones with distinct in vitro and in vivo differentiation/mineralization potential. J Bone Miner Res 1999, 14(6):893–903. 10.1359/jbmr.1999.14.6.893

    Article  CAS  PubMed  Google Scholar 

  39. Qian X, Karpova T, Sheppard AM, McNally J, Lowy DR: E-cadherin-mediated adhesion inhibits ligand-dependent activation of diverse receptor tyrosine kinases. EMBO J 2004, 23(8):1739–1748. 10.1038/sj.emboj.7600136

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Magid A, Rotman SR, Weiss AM: Comments on Picture thresholding using an iterative selectionmethod. ystems, Man and Cybernetics, IEEE Transactions on 1990, 20(5):1238–1239. 10.1109/21.59988

    Article  Google Scholar 

  41. Meyer F, Beucher S: Morphological segmentation. Journal of Visual Communication and Image Representation 1990, 1(1):21–46. 10.1016/1047-3203(90)90014-M

    Article  Google Scholar 

  42. Ridler TW, Calvard S: Picture thresholding using an iterative selection method. IEEE Transactions on Systems, Man and Cybernetics 1978, SMC-8(8):630–632.

    Google Scholar 

  43. Sezgin M, Sankur B: Survey over image thresholding techniques and quantitative performance evaluation. Journal of Electronic Imaging 2004, 13(1):146–168. 10.1117/1.1631315

    Article  Google Scholar 

  44. Trussell HJ: Comments on 'Picture thresholding using an iterative selection method'. IEEE Transactions on Systems, Man and Cybernetics 1979, SMC-9(5):311. 10.1109/TSMC.1979.4310204

    Article  Google Scholar 

  45. Yanni MK, Horne E: A new approach to dynamic thresholding. Eur Assoc Signal Process: Edinburgh, UK 1994, 34–41.

    Google Scholar 

  46. Carlin BP, Louis TA: Bayes and Empirical Bayes methods for data analysis. Volume 47. Taylor & Francis CRC Press; 2000.

    Chapter  Google Scholar 

  47. French S, Smith JQ: The practice of Bayesian analysis. New York: J. Wiley & Sons; 1997.

    Google Scholar 

  48. Martin JY, Schwartz Z, Hummert TW, Schraub DM, Simpson J, Lankford J, Dean DD, Cochran DL, Boyan BD: Effect of Titanium Surface-Roughness on Proliferation, Differentiation, and Protein-Synthesis of Human Osteoblast-Like Cells (Mg63). Journal of Biomedical Materials Research 1995, 29(3):389–401. 10.1002/jbm.820290314

    Article  CAS  PubMed  Google Scholar 

Download references


We gratefully acknowledge support from NIH Grant Numbers RR17425 and HK072039.

Author information

Authors and Affiliations


Corresponding author

Correspondence to J Carson Meredith.

Additional information

Authors' contributions

JCM directed the experimental design, data analysis method development, and writing of the manuscript. JS performed most cell culture experiments and implemented the LCM method, including programming, and wrote the manuscript. PJZ helped in data interpretation and statistical analysis. CCC performed the calcium-depletion experiments and analysis. All authors have read and approved this manuscript.

Electronic supplementary material


Additional file 1: Supplementary materials. Supplementary theoretical and calculation details, as well as additional data on image analysis and global analysis are provided. (DOC 476 KB)


Additional file 2: Source codes for local cell metrics. This .zip archive folder contains sources codes for implementing local cell metrics provided under terms of the GNU open source license. (ZIP 487 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Su, J., Zapata, P.J., Chen, CC. et al. Local cell metrics: a novel method for analysis of cell-cell interactions. BMC Bioinformatics 10, 350 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: