# A boundary delimitation algorithm to approximate cell soma volumes of bipolar cells from topographical data obtained by scanning probe microscopy

- Patrick Happel
^{1, 2}Email author, - Kerstin Möller
^{2, 3}, - Ralf Kunz
^{1}and - Irmgard D Dietzel
^{2}

**11**:323

**DOI: **10.1186/1471-2105-11-323

© Happel et al; licensee BioMed Central Ltd. 2010

**Received: **10 March 2010

**Accepted: **15 June 2010

**Published: **15 June 2010

## Abstract

### Background

Cell volume determination plays a pivotal role in the investigation of the biophysical mechanisms underlying various cellular processes. Whereas light microscopy in principle enables one to obtain three dimensional data, the reconstruction of cell volume from *z*-stacks is a time consuming procedure. Thus, three dimensional topographic representations of cells are easier to obtain by scanning probe microscopical measurements.

### Results

We present a method of separating the cell soma volume of bipolar cells in adherent cell cultures from the contributions of the cell processes from data obtained by scanning ion conductance microscopy. Soma volume changes between successive scans obtained from the same cell can then be computed even if the cell is changing its position within the observed area. We demonstrate that the estimation of the cell volume on the basis of the width and the length of a cell may lead to erroneous determination of cell volume changes.

### Conclusions

We provide a new algorithm to repeatedly determine single cell soma volume and thus to quantify cell volume changes during cell movements occuring over a time range of hours.

## Background

Cell volume regulation occurs in a wide variety of tissues from kidney to brain [1–4]. Although much is known about ion and water fluxes involved in many regulatory processes, no method has so far been designed to investigate potential volume changes in moving cells. Light microscopy enables one to estimate the cellular volume via different techniques, ranging from extrapolation on the basis of the width and length of the cell [5], changes in light intensity and light scattering [6, 7], various staining techniques [8, 9] to quantitative phase microscopy [10]. All these techniques fail when it is required to investigate the volume of a cell undergoing notable changes in shape such as occur during cell migration [11, 12] since they require constant parameters such as height or refractive index and some have additional disadvantages such as bleaching of the dye [13].

A promising approach to circumvent these problems is to measure volume directly with a scanning probe microscope. Direct measurements of cellular volume have been performed by scanning ion conductance microscopy (SICM) on cellular layers [14] and single cells [15, 16] and by atomic force microscopy (AFM) on living and fixed cells [17, 18]. The volume determined by SICM of cells forming a confluent layer has been validated by scanning confocal laser microscopy [14].

We have previously proposed a method for distinguishing between cellular somata and processes during investigations of the surface of oligodendrocyte cell bodies at different developmental stages [19]. Here, every pixel exceeding a certain height was assigned to the cell soma. In cells undergoing marked changes in shape this method fails since it would result in different soma volumes if a cell flattens but performs a compensatory widening thus maintaining its volume. Hence, to estimate volume changes of single bipolar cell somata changing their shape and position we have now developed a novel procedure that allows us to separate the cell soma volume from the extended peripheral membrane processes of bipolar cells.

## Results and Discussion

### Approximation of the position of the nucleus

Atomic force microscopy measurements on hippocampal neurons revealed that the higher parts of the cell body form a harder structure and correspond most likely to the nucleus [20]. In order to determine a single point that represents the location of the nucleus the following procedure was employed: We stained the nucleus using Hoechst 33342 dye and recorded an epifluorescence as well as a phase contrast image.

*P*

_{ i }= (

*x*

_{ i },

*y*

_{ i },

*z*

_{ i }) (with

*i*denoting the number of the pixel) exceeding a certain height

*T z*

_{max}, where

*T*denotes a predefined threshold and

*z*

_{max}denotes the maximum cell height. To calculate the position of the centroid

*C*

_{ T }we reduced the

*z*-coordinates of

*P*

_{ i }to boolean values

*z*

_{ T,i }= [

*z*

_{ i }>

*T z*

_{max}]. The square brackets indicate a Heaviside-like function that yields 1 if the enclosed condition is true and 0 otherwise [21, 22]. Furthermore, we assumed constant step sizes between the pixels and thus calculated the

*x*-coordinate of

*C*

_{ T }, , as:

We next investigated the distance between *C*_{
T
}and the centroid of the Hoechst 33342 staining of the nucleus (see Methods section) for various thresholds *T*.

*C*

_{ T }for

*T*= 0.1, 0.15,..., 0.9 and the centroid of the nucleus obtained from the epifluorescence staining (

*C*

_{fluo}) are depicted in Figure 3D.

*C*

_{90}(note that we use

*T*in percent when indexing or labeling, thus

*C*

_{T = 0.9}≡

*C*

_{90}) exhibited the minimal distance to

*C*

_{fluo}. Figure 3E shows the average distances between

*C*

_{ T }and

*C*

_{fluo}obtained from three different recordings. This confirms that

*C*

_{90}is located closest to

*C*

_{fluo}. Note that representations determined by using a larger threshold such as

*C*

_{95}often base on disjunct areas and were not investigated in detail. Thus we used

*C*

_{90}to approximate the position of the nucleus in the following.

### Estimation of the heading direction of the cell

*θ*

_{h}of an OPC by rotating a straight line

*C*

_{90}as the approximation of a straight line through the nucleus. In order to determine the heading direction of the cell we considered the arcs from each pixel representing the cell to

*y*(

*x*,

*θ*) Let

*ϕ*

_{ i }(

*θ*) denote the smallest angle between

*P*

_{ i }and

*y*(

*x*,

*θ*) and

*r*

_{ i }denote the distance from

*C*

_{90}to

*P*

_{ i }. Then the length

*s*

_{ i }(

*θ*) of the corresponding arc is calculated as

*s*

_{ i }(

*θ*) =

*ϕ*

_{ i }(

*θ*)

*r*

_{ i }. Figure 4A illustrates the relations between the introduced angles, lines and points for two different pixels

*P*

_{ i }located at opposite sides of

*C*

_{90}. We now defined the angle

*θ*

_{h}, that minimized the sum of

*s*

_{ i }(

*θ*) and thus matched the condition

as the heading direction of the cell. Here we assumed that pixels that exhibited a height of ≤1 μm represented the cell culture dish rather than the cell. Equation (3) was solved numerically by testing all angles 0 ≤ *θ* ≤ *π* in steps of ∆*θ* = 2π/360.

### Rotating and interpolating the data

*C*

_{90}was shifted into the origin of the new coordinate system. We denote the axes of the new coordinate system as

*x'*-,

*y'-*and

*z'*- axes and a rotated and translated pixel as , with

*j*indicating the number of the pixel in the rotated scan. To determine the lateral extent of the rotated scan we considered the distances of the vertices of the original scan and

*y*(

*x*,

*θ*

_{h}) or a straight line through

*C*

_{90}perpendicular to

*y*(

*x*,

*θ*

_{h}) as illustrated in Figure 4C. Since the approximation of the single line boundaries of the cell soma required lines of data points parallel to the heading direction of the cell, we defined the grid consisting of the projections of

*Q*

_{ j }to the

*x'*,

*y'*plane of the rotated and translated scan such that

Here
is the negative representation of the length
as a coordinate, Δ*x* and Δ*y* denote the step sizes of the original scan in the *x*- and *y*-directions, respectively, and the truncated square brackets represent the *ceil* and the *floor* functions [22, 23].

*z'*- coordinate of a pixel

*Q*

_{ j }we rotated its projection into the original scan dataset by applying the inverse rotation matrix

*P*

_{ i }to the

*x*-,

*y*-plane) that surrounded as depicted in Figure 5B. The

*z*-coordinates of the corresponding pixels were known from the original data. Each set of three out of these four projections defines a triangle as indicated by the dotted lines in Figure 5B. In the following we refer to the four triangles as

**M**

_{ k }(with

*k*= 1, 2, 3, 4) and to the vertices of one triangle as with

*l*= 1, 2, 3. We selected

*l*such that the right angle was located at and furthermore such that and . An example is shown in Figure 5C. If and only if was located inside

**M**

_{ k }the sum ζ

_{ k }of the angles at to the vertices of

**M**

_{ k }amounted to 2

*π*[24].

*M*

_{ k,l }that corresponded to the projections . The

*z*-value

*z*

_{ k }(

*x, y*) of this plane at a position (

*x, y*) is given by

**M**

_{k}:

### Approximation of the contour of a single data row

*y'*) separately. The corresponding

*y'*-values were defined by equation (4). Figure 6 shows sketches of the contours of two characteristic cell shapes; an almost circular cell body that is easy to distinguish from the cell processes (Figure 6A) and a cell soma that protruded into the direction of one of the extensions (Figure 6B). Thus, as indicated in Figure 6B, we assumed that a polynomial of third degree was convenient to approximate the cell soma contour but still suitable to crop the cell process.

*y'*-level we subdivided the data into positive and negative, or frontal and rear, parts with respect to the corresponding

*x'*-coordinates. In the following we describe the fitting procedure for the positive part, thus

*x'*> 0. was defined as the projection of

*Q*

_{ j }to the

*x'*,

*z'*-plane and furthermore with

*p*= 0, 1, 2,... as the set of projections at a constant

*y'*such that for all

*p*> 0

*z'*-coordinates (since the data points were filtered this is equivalent to

*z'*> 1 μm, see Methods section). In general

*n*+ 2 data points are needed to fit a polynomial of

*n*th degree (

*n*+ 1 data points define the polynomial). Furthermore, we assumed that the cell body is represented by the data points whose

*x'*-coordinates are located close to zero. Thus we additionally tested whether there was no gap within and it therefore matched the condition

Otherwise, data points with *x'*-coordinates close to zero existed with *z'* = 0. This most likely occured at the borders of the cell soma in ± *y'*-direction and was treated as a special case described later in this section.

*n*th degree to the data we used the function

*fit*from Matlab's Curve Fitting Toolbox that implements a least square algorithm [25, 26]. It provides, among others, the value that represents the goodness of the fit considering the number of data points that were approximated by the fit. We investigated the goodness of the fits to an increasing number

*r*of data points. We refer to the subset of

**S**

_{ y' }that contains the first

*r*elements as and we denote the goodness of the fit to

**S**

_{ r,y' }as Additionally, we defined

*X*

_{ y' }(

*r*) to be the smallest, positive, non-complex root of the polynomial that was determined by the function

*fit*. We approximated the polynomial boundary of the cell soma for each line segment towards the direction of fitting as the

*X*

_{ y' }(

*r*) that matched the condition

*r*=

*n*+ 1,

*n*+ 2,...,

*p*

_{max}. Here

*p*

_{max}denotes the largest index

*p*of the projections included in

**S**

_{ y' }Figure 7 shows examples of the fitting procedure for

*r*= 4, 8, 9 and 14, respectively, with

*n*= 3, hence fitting polynomials of third degree. For

*r*= 4 and

*r*= 14 (Figure 7A and 7D)

*F*

_{ y' }(

*r*) had no real root with a corresponding positive

*x'*-coordinate, thus these fits were not taken into consideration. Since (Figure 7B and 7C)

*X*

_{ y' }(

*r*= 8) (indicated by the red arrow-head in Figure 7C) was used to approximate the cell soma boundary at the corresponding

*y'*-level. Note that the goodness of the fit to

**S**

_{ 8,y' }was larger than those of all other fits that exhibited

*X*

_{ y' }(

*r*≠ 8) but are not shown in Figure 7 for clarity.

If the procedure failed to determine a cell soma boundary for the investigated set of data points **S**_{
y'
}no *r* with a corresponding *X* _{
y'
}(*r*) existed. We then defined the boundary to be *X*_{
y'
}(*r = n*), if it existed. Note that
(*r* = *n*) is not defined [26]. If *X*_{
y'
}(*r* = *n*) did also not exist we repeated the procedure with *n* : = *n* - 1 as long as *n* > 1, thus fitting polynomials of a reduced degree. In all cases investigated this procedure led to detection of bordering pixels.

*r*=

*n*) existed as well as the test of whether

*n*> 1, indicated by the dotted arrow in the lower right part of the chart. This procedure was named

*fitBest*.

#### Special cases of the fitting procedure

As indicated in Figure 8 an error was returned if the investigated set of data points did not match the conditions listed in equation (9). In this case data points with a corresponding
existed within the first *n* + 1 data points in the fit direction. This most likely occurred at the borders of the cell soma in ± *y'*-direction. This special situation might occur under two conditions. In the first case the cell body approximates to a circular shape causing the boundary perpendicular to the direction of fitting to consist of only a few pixels. Furthermore, the number of pixels available to the fitting procedure as depicted in Figure 7 is decreased by the division of the cell into its frontal and its rear part. Secondly, OPCs in a later stage of development might exhibit small additional extensions that grow perpendicularly to the heading direction.

It was important to consider these cases in order to provide an errorless and thus automatic processing of the data. There are different strategies to determine the boundary of the cell soma at these locations depending either on the chosen degree of the polynomial fitted to **S**_{
y'
}as well as whether potential extensions at these sides of the cell soma should be included or excluded from the soma approximation. The most restrictive and simple solution would be to omit and thus to crop these lines.

To obtain a more accurate fit and to include potential cell extensions at these sides we introduced three more functions: *fitOnePoint*, *fitTwoPoints* and *fitThreePoints* that were executed depending on the number of data points with *z'* > 0. We considered the set of pixels
that matched all conditions listed in equation (8) except one: The *z'*-coordinate was not tested, thus **T**_{
y'
}might also include projections with *z'* = 0. Let
be the number of projections with a *z'*-coordinate exceeding zero. If *N* = 4 we executed the function *fitBest*. If *N* = 1, *N* = 2 or *N* = 3 we executed the functions *fitOnePoint*, *fitTwoPoints* or *fitThreePoints*, respectively. Note that these functions might result in more than one boundary for the particular *y'* level, thus the resulting approximated cell soma might appear jagged.

The simplest case is *N* = 1 and the corresponding function *fitOnePoint*. We refer to the non-zero data point as
and used the roots of a parabola through
as the boundary if *u* < 4, otherwise the line was cropped.

Let the two non-zero projections be
and
with *u* <*v* in the case of *N* = 2 (function *fitTwoPoints*). We first considered the case *v - u* = 1, hence, the two points were neighbors. We fitted a polynomial of third degree to
and
and used its roots as the boundary in this case except if *v* = 4 and
. In the latter case we cropped the structure assuming that it did not belong to the cell soma.

If *v* - *u* > 1 we only assigned
to the cell soma and approximated the contour of the cell soma by the roots of the parabola through
as in the function *fitOnePoint*.

The most complicated case was *N* = 3. We refer to the single projection with zero *z'*-coordinate as
In this case the approximation was performed differently for varying values of *u*. If *u* = 1 we considered the *z'*-coordinate of the projection
. If
we assumed that the cell soma exhibited an asymmetric shape and applied the function *fitBest*. Otherwise, if
we approximated the cell soma boundary for the particular *y'* by the roots of a polynomial of third degree fitted to
.

If *u* ∈ {2, 3} we applied *fitOnePoint* to the single, non-zero projection and *fitTwoPoints* to the two neighboring non-zero projections, respectively. If *u* = 4 we considered the *z'*-coordinate of the first point opposite to the direction of fitting,
. If
we applied *fitBest*, otherwise we approximated the cell soma boundary at the current *y'*-level by the roots of a polynomial of third degree fitted to
.

### Approximation of the volume of the cell soma

To approximate the cell soma volume we summed the *z*-coordinates of every pixel located within the approximated boundaries of the cell soma. This required that the height of every pixel located within the approximated cell soma boundary was known. Hence, if a single delimitation of the cell soma was located outside the original scan we were not able to approximate the cell soma volume and the recording was discarded. This happened if the cell body was in part located outside of the SICM image or very close to its borders.

### Evaluation of the procedure

To evaluate the BDA we simulated objects of known volume and applied the morphometric fitting procedure to investigate any potential effect of geometry on the volume determinations. We have previously determined the restrictions of scan size and resolution for the successful investigation of migrating OPCs [27]. In brief, to image migrating OPCs with a suitable frame rate using our present SICM the dimensions of the recordings had to be restricted to 30 μm squares with a lateral step size of 1 μm, limiting the SICM images to 900 pixels.

We first applied the BDA to a hemisphere with a radius of *r*_{0} = 5 pixels (since the length of the cell body of an OPC is approximately 10 μm) in a data set consisting of 900 pixels as depicted in Figure 9A. The volume *V*_{comp} computed by the BDA (omitting the determination of a heading direction as well as rotation and translation) was the same as the volume *V*_{sum} calculated by summing the volume of the columns above each pixel.

We next compared the determination of the volume of an half-ellipsoid with the two methods. A possible effect of the direction of fitting was tested by applying the BDA to an ellipsoid defined by the three radii *r*_{
x
}, *r*_{
y
}and *r*_{
z
}with *r*_{
x
}> r_{y} and vice versa, as depicted in Figure 9B and 9C (the corresponding radii are *r*_{
x
}= 0.8*r*_{0}, *r*_{
y
}= 1.25*r*_{0}, *r*_{
z
}= *r*_{0} and *r*_{
x
}= 1.25*r*_{0}, *r*_{
y
}= 0.8*r*_{0}, *r*_{
z
}= *r*_{0}). Again, no difference was found between *V*_{comp} and *V*_{sum}.

To investigate whether the BDA in principle allows one to determine the volume of an object that flattens but maintains its volume by a compensatory widening we computed the volumes of an ellipsoid defined by the radii *r*_{
y
}= *r*_{0}, *r*_{
x
}= *t r*_{0} and *r*_{
z
}= *r*_{0}/*t* with 1 ≤ *t* ≤ 2 in step sizes of Δ*t* = 0.05. Figure 9G (blue crosses) shows the computed volume normalized to *V*_{sum} for every investigated value of *t*. There is no difference between *V*_{comp} and *V*_{sum}, thus *V*_{
n
}= 1. In contrast, the computed volume did not match *V*_{sum} when it was determined by using the method that every pixel exceeding a predefined threshold was assigned to the cell soma [16, 19]. The normalized volumes are displayed in Figure 9G (red dots and cross-hairs) for an absolute and a relative threshold. In the following we only consider the determination using a relative threshold since it is clearly visible that the use of an absolute threshold leads to increasing differences in the determination of the soma with increasing elongation of the ellipsoid. Additionally, we observed no difference in the volume determined by the BDA and *V*_{sum} when varying *r*_{
y
}instead of *r*_{
x
}or when varying both lateral radii by defining *r*_{
x
}= *r*_{
y
}= *t*^{
1/2
}*r*_{0}.

To simulate a bipolar cell we added extensions in ± *x'*-direction to a hemisphere of radius *r*_{0} as well as to the ellipsoids. Images of the resulting objects are depicted in Figure 9D-F. The height of the extension was chosen as *r*_{0}/2 and its width as 2 *r*_{0}/5. Every *z*-value at the corresponding positions was adjusted to *r*_{0}/2 if the *z*-value calculated by equation (11) (see Methods section) was below *r*_{0}/2. This avoids a gap between the half-ellipsoid and the extension but also increases the *z*-value of some pixels of the half-ellipsoid such that the volume differs from the volume *V*_{sum} computed by summing the *z*-values of the mere half-ellipsoid as depicted in Figure 9I. To our knowledge no exact definition exists describing where the cell soma ends and the cell process starts. At positions where the soma merges into the neurite a gradual decline of the soma and a corresponding increase of soma volume most likely occurs (Figure 9I).

Here we chose to use the calculated volume of the half-ellipsoid without extension as reference. Since we calculated the soma volume by summing all *z*-values corresponding to pixels within the approximated soma boundary an overestimation of the soma volume at positions merging into the neurites (Figure 9I) could be induced by the BDA.

The approximated volume, normalized to *V*_{sum} of the corresponding hemisphere or ellipsoid without extension, is shown in Figure 9H. As expected, the BDA (blue bars in Figure 9H) overestimates the volume with respect to *V*_{sum}. In contrast, the approximation via the threshold method [16, 19], in this case applied using a threshold of *r*_{0}/2, underestimates the volume with respect to *V*_{sum} since it omits all sections of the ellipsoid with a height below the selected threshold. Putative cell shape changes as depicted in Figure 9E-G would result in detections of relative soma volumes as indicated in the gray boxes in Figure 9H. Erroneous changes due to different shapes are indicated by the arrows. Both methods lead to almost similar errors (about 5%) in the determination of soma volume changes.

*h r*

_{0}, normalized to

*V*

_{sum}of the hemisphere without extensions. As expected from the result shown in Figure 9H, our method overestimates the volume with increasing height. An increase in the height of the extension from

*h*= 0.2 to

*h*= 0.6 results in an erroneous detection of a soma volume increase of about 6%. Although the relation seems to be linear in the depicted range, it is more complex: A threefold increase in the height of the extension from

*h*= 0.3 to

*h*= 0.9 leads to an erroneous detection of a volume increase of about 9% whereas a threefold increase from

*h*= 0.1 to

*h*= 0.3 leads to an erroneous detection of a volume increase of about 3%. In contrast, the thresholding method (red cross-hairs) shows an underestimation of the soma volume that increases stepwise but maintains a constant volume over a range of heights. However, the stepwise decrease of the calculated volume and thus the determination of a constant volume over a certain range of heights results from the imprecision that occurs due to the rasterization of the sphere as shown by the investigation of a simulated scan with a tenfold resolution (red dots in Figure 10A). An increase in the height of the extension from

*h*= 0.2 to

*h*= 0.6 results in an erroneous detection of a volume decrease of about 16% for the low resolution simulation and of 32% in the high resolution simulation. We observed similar results when performing the same investigation on the objects depicted in Figure 9F and 9G with only slight differences in the amount of errors determined by the two methods.

We used *h r*_{0} as threshold in these investigations. Note that the height of the processes of a live cell is much more difficult to determine due to the more complex and irregular shape and thus adds additional uncertainties to the determination of the soma volume.

Figure 10B shows the impact of various widths of the extension on the soma volume determination. As expected, the overestimation of the soma volume increases with increasing extension width *w* 2 *r*_{0}. A fourfold widening of the extension leads to an erroneous determination of a soma volume increase of 11%. Since the height of the extension defines the threshold for the thresholding method the increasing width of the extension is not detected by this method. Thus it computes a constant volume under these conditions.

Figure 10C and 10D show the impact of a combined variation of the radii and the height of the extension. We investigated the radii *r*_{
x
}(*t*) = *r*_{
y
}(*t*) = *t*^{1/2} *r*_{0} and *r*_{
z
}(*t*) = *r*_{0}/*t* for 1 ≤ *t* ≤ 2 and the fraction *h* of the height *h* *r*_{z} (*t*) of the extension for 0.2 ≤ *h* ≤ 0.6. Particularly when minor changes in shape were simulated, the BDA (Figure 10C) detects smaller erroneous volume changes compared with the thresholding method (Figure 10).

### Application to live cells

In the second recording the cell changed its shape to be more circular and thus the determined basal area of the cell soma only shows a slight expansion as clearly visible by the comparison of Figure 12Ab and 12Bb. Note that it is known that migrating OPCs show an average velocity of 10 μm/h [28] and that migrating OPCs show notable changes in shape [27]. The detailed analysis of the parameters of the cell shape and soma shows that the cell swelled by approximately 29% and that this swelling was dominated by an increase in cell height whereas the length of the soma decreased. Most notably, this challenges the method to approximate the volume of a cell from light microscopic images by interpolation via the width and the length of its soma. This approximation, in contrast, would detect a slight cell shrinkage since the basal area covered by the cell soma was reduced as depicted in Figure 12E. The separate consideration of the frontal and rear soma volume by dividing the cell soma at the level of *C*_{90} perpendicular to *y* (*x*, *θ*_{h}) yields that the volume increase is dominated by an increase in the frontal volume (Figure 12E). Figure 12F summarizes the changes in the lateral dimensions as well as the changes of area and volume between both scans.

## Conclusions

We provide a new algorithm to approximate the basal area of the cell soma for volume determinations of bipolar cells undergoing temporal changes in shape. Simulations show that this method provides smaller errors in the detection of soma volume changes in moving cells than the method to separate soma areas from neurites via height thresholds [16, 19]. We show that the algorithm can be applied successfully to detect soma volume changes of living oligodendrocyte precursor cells. Additionally we demonstrate that the approximation of the volume using two-dimensional data such as cell body length and width may lead to erroneous results.

## Methods

### Scanning ion conductance microscopy

Cells were obtained as described in [27] and images were acquired using a pulse mode SICM operating in floating backstep mode as detailed in [15, 29]. Scans were performed in Leibovitz-15 medium using scanning probes with an access resistance of about 4 MΩ. Scanning probes were filled with extracellular saline containing (in mM): NaCl 110, KCl 5.4, CaCl_{2} 1.8, MgCl_{2} 0.8, HEPES 10, glucose 10. Step sizes were 1 μm in lateral and 100 nm in vertical direction, frame acquisition time was about 10 minutes.

### Determination of the position of the nucleus

Staining of the nucleus was performed by applying Hoechst 33342 dye (20 μg/mL, in Leibovitz-15 medium) for 5 minutes before the SICM measurement was started. Fluorescence images were obtained with a Zeiss Axiovert microscope at an epifluorescence wavelength of 365 nm. The boundaries of the nucleus were obtained by applying the Sobel algorithm and thresholding. Centroids were calculated in the same manner as described above (equation (1)) but using a color value as threshold. The position of the scanning area within the light microscopic image was determined as described previously [19].

### Data processing

SICM data was processed using Matlab (R2008a) including the Curve Fitting Toolbox [25]. Data shown were interpolated by cubic splines unless otherwise noted. All *z*-data were plane corrected if necessary and filtered with a threshold filter setting each *z*-value less than 1 μm to 0 μm before the BDA was applied.

### Creation of half-ellipsoids

*z*-coordinates of the half-ellipsoids were computed as

Here the operator ℜ was defined to return the real part of a complex number *a* + *b* i, hence *a* = ℜ(*a* + *b* i) [22] and *r*_{
x
}, *r*_{
y
}and *r*_{
z
}denote the radii of the ellipsoid in the corresponding directions.

### Availability and requirements

**Project name:** BDA - A boundary delimitation algorithm

**Project homepage:**
https://sourceforge.net/projects/behindlight/files/BDA/

**Operating system(s):** Platform independent

**Programming language:** Matlab

**Other requirements:** Matlab R2008a with Curve Fitting Toolbox

**License:** GNU-LGPL

**Restrictions for commercial usage:** Notify the corresponding author of this report and include a citation to this report in the accompanying documentation.

## Declarations

### Acknowledgements

Some of the equipment used was initially funded by the DFG-grant Schuh929/5-1 and we acknowledge the support of a part of the work by a fellowship of the *Allgemeines Graduiertenkolleg der Ruhr-Universität Bochum* to PH. We thank Heidrun Breuker-Siraj for help with the cell culture, Jonathan Coles, Gerd Hoffmann, Stefan A. Mann and Bernd Kind for critical comments on the manuscript as well as Rolf Heumann, Jan Meijer and the RUBION-team for continuous support.

## Authors’ Affiliations

## References

- Lang F, Ritter M, Völkl H, Häussinger D: The biological significance of cell volume.
*Ren Physiol Biochem*1993, 16(12):48–65.PubMedGoogle Scholar - Wehner F, Olsen H, Tinel H, Kinne-Saffran E, Kinne RKH: Cell volume regulation: osmolytes, osmolyte transport, and signal transduction.
*Rev Physiol Biochem Pharmacol*2003, 148: 1–80. full_textPubMedGoogle Scholar - Dietzel I, Heinemann U, Hofmeier G, Lux HD: Transient changes in the size of the extracellular space in the sensorimotor cortex of cats in relation to stimulus-induced changes in potassium concentration.
*Exp Brain Res*1980, 40(4):432–9. 10.1007/BF00236151View ArticlePubMedGoogle Scholar - Dietzel I, Heinemann U, Lux HD: Relations between slow extracellular potential changes, glial potassium buffering, and electrolyte and cellular volume changes during neuronal hyperactivity in cat brain.
*Glia*1989, 2: 25–44. 10.1002/glia.440020104View ArticlePubMedGoogle Scholar - Hall SK, Zhang J, Lieberman M: An early transient current is associated with hyposmotic swelling and volume regulation in embryonic chick cardiac myocytes.
*Exp Physiol*1997, 82: 43–54.View ArticlePubMedGoogle Scholar - Farinas J, Kneen M, Moore M, Verkman AS: Plasma membrane water permeability of cultured cells and epithelia measured by light microscopy with spatial filtering.
*J Gen Physiol*1997, 110(3):283–96. 10.1085/jgp.110.3.283View ArticlePubMedPubMed CentralGoogle Scholar - Echevarria M, Verkman AS: Optical measurement of osmotic water transport in cultured cells. Role of glucose transporters.
*J Gen Physiol*1992, 99(4):573–89. 10.1085/jgp.99.4.573View ArticlePubMedGoogle Scholar - Kao HP, Verkman AS: Tracking of single fluorescent particles in three dimensions: use of cylindrical optics to encode particle position.
*Biophys J*1994, 67(3):1291–300. 10.1016/S0006-3495(94)80601-0View ArticlePubMedPubMed CentralGoogle Scholar - Droste MS, Biel SS, Terstegen L, Wittern KP, Wenck H, Wepf R: Noninvasive measurement of cell volume changes by negative staining.
*J Biomed Opt*2005, 10(6):064017. 10.1117/1.2138011View ArticlePubMedGoogle Scholar - Curl CL, Bellair CJ, Harris PJ, Allman BE, Roberts A, Nugent KA, Delbridge LM: Single cell volume measurement by quantitative phase microscopy (QPM): a case study of erythrocyte morphology.
*Cell Physiol Biochem*2006, 17(5–6):193–200. 10.1159/000094124View ArticlePubMedGoogle Scholar - Saadoun S, Papadopoulos MC, Hara-Chikuma M, Verkman AS: Impairment of angiogenesis and cell migration by targeted aquaporin-1 gene disruption.
*Nature*2005, 434(7034):786–92. 10.1038/nature03460View ArticlePubMedGoogle Scholar - Dieterich P, Klages R, Preuss R, Schwab A: Anomalous dynamics of cell migration.
*Proc Natl Acad Sci USA*2008, 105(2):459–63. 10.1073/pnas.0707603105View ArticlePubMedPubMed CentralGoogle Scholar - Satoh H, Delbridge LM, Blatter LA, Bers DM: Surface:volume relationship in cardiac myocytes studied with confocal microscopy and membrane capacitance measurements: species-dependence and developmental effects.
*Biophys J*1996, 70(3):1494–504. 10.1016/S0006-3495(96)79711-4View ArticlePubMedPubMed CentralGoogle Scholar - Korchev Y, Gorelik J, Lab M, Sviderskaya E, Johnston C, Coombes C, Vodyanoy I, Edwards C: Cell volume measurement using scanning ion conductance microscopy.
*Biophys J*2000, 78: 451–7. 10.1016/S0006-3495(00)76607-0View ArticlePubMedPubMed CentralGoogle Scholar - Happel P, Hoffmann G, Mann S, Dietzel ID: Monitoring cell movements and volume changes with pulse-mode scanning ion conductance microscopy.
*J Microsc*2003, 212(Pt 2):144–51. 10.1046/j.1365-2818.2003.01248.xView ArticlePubMedGoogle Scholar - Mann SA, Versmold B, Marx R, Stahlhofen S, Dietzel ID, Heumann R, Berger R: Corticosteroids reverse cytokine-induced block of survival and differentiation of oligodendrocyte progenitor cells from rats.
*J Neuroinflammation*2008, 5: 39. 10.1186/1742-2094-5-39View ArticlePubMedPubMed CentralGoogle Scholar - Schneider SW, Pagel P, Rotsch C, Danker T, Oberleithner H, Radmacher M, Schwab A: Volume dynamics in migrating epithelial cells measured with atomic force microscopy.
*Pflugers Arch*2000, 439(3):297–303. 10.1007/s004240050943View ArticlePubMedGoogle Scholar - Fabian A, Fortmann T, Dieterich P, Riethmüller C, Schön P, Mally S, Nilius B, Schwab A: TRPC1 channels regulate directionality of migrating cells.
*Pflugers Arch*2008, 457(2):475–84. 10.1007/s00424-008-0515-4View ArticlePubMedGoogle Scholar - Mann SA, Meyer JW, Dietzel ID: Integration of a scanning ion conductance microscope into phase contrast optics and its application to the quantification of morphological parameters of selected cells.
*J Microsc*2006, 224(Pt 2):152–7. 10.1111/j.1365-2818.2006.01693.xView ArticlePubMedGoogle Scholar - Yunxu S, Danying L, Yanfang R, Dong H, Wanyun M: Three-dimensional structural changes in living hippocampal neurons imaged using magnetic AC mode atomic force microscopy.
*J Electron Microsc (Tokyo)*2006, 55(3):165–72. 10.1093/jmicro/dfl013View ArticleGoogle Scholar - Knuth DE: Two notes on notation.
*AMER MATH MONTHLY 99*1992, 5: 403. 10.2307/2325085View ArticleGoogle Scholar - Graham RL, Knuth DE, Patashnik O:
*Concrete mathematics: a foundation for computer science*. Boston, MA, USA: Addison-Wesley Longman Publishing Co., Inc; 1994.Google Scholar - Iverson KE:
*A programming language*. New York, NY, USA: John Wiley & Sons, Inc; 1962.View ArticleGoogle Scholar - Sutherland IE, Sproull RF, Robert , Schumacker A: A characterization of ten hidden-surface algorithms.
*ACM Computing Surveys*1974, 6: 1–55. 10.1145/356625.356626View ArticleGoogle Scholar - Mathworks: Curve Fitting Toolbox - Documentation.[http://www.mathworks.com/access/helpdesk/help/toolbox/curvefit/]
- Mathworks: Fit model to data - MATLAB.[http://www.mathworks.com/access/helpdesk/help/toolbox/curvefit/fit.html]
- Happel P, Dietzel ID: Backstep scanning ion conductance microscopy as a tool for long term investigation of single living cells.
*J Nanobiotechnology*2009, 7: 7. 10.1186/1477-3155-7-7View ArticlePubMedPubMed CentralGoogle Scholar - Schmidt C, Ohlemeyer C, Labrakakis C, Walter T, Kettenmann H, Schnitzer J: Analysis of motile oligodendrocyte precursor cells in vitro and in brain slices.
*Glia*1997, 20(4):284–98. 10.1002/(SICI)1098-1136(199708)20:4<284::AID-GLIA2>3.0.CO;2-6View ArticlePubMedGoogle Scholar - Mann SA, Hoffmann G, Hengstenberg A, Schuhmann W, Dietzel ID: Pulse-mode scanning ion conductance microscopy-a method to investigate cultured hippocampal cells.
*J Neurosci Methods*2002, 116(2):113–7. 10.1016/S0165-0270(02)00023-7View ArticlePubMedGoogle Scholar

## Copyright

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