Improved method for the quantification of motility in glia and other morphologically complex cells

Cells such as astrocytes and radial glia with many densely ramified, fine processes pose particular challenges for the quantification of structural motility. Here we report the development of a method to calculate a motility index for individual cells with complex, dynamic morphologies. This motility index relies on boxcar averaging of the difference images generated by subtraction of images collected at consecutive time points. An image preprocessing step involving 2D projection, edge detection, and dilation of the raw images is first applied in order to binarize the images. The boxcar averaging of difference images diminishes the impact of artifactual pixel fluctuations while accentuating the group-wise changes in pixel values which are more likely to represent real biological movement. Importantly, this provides a value that correlates with mean process elongation and retraction rates without requiring detailed reconstructions of very complex cells. We also demonstrate that additional increases in the sensitivity of the method can be obtained by denoising images using the temporal frequency power spectra, based on the fact that rapid intensity fluctuations over time are mainly due to imaging artifact. The MATLAB programs implementing these motility analysis methods, complete with user-friendly graphical interfaces, have been made publicly available for download.