Computer-assisted quantification of axo-somatic boutons at the cell membrane of motoneurons

This paper presents a system for computer-assisted quantification of axe-somatic boutons at motoneuron cell-surface membranes. Different immunohistochemical stains can be used to prepare tissue of the spinal cord. Based on micrographs displaying single neurons, a finite element balloon model has been applied to determine the exact location of the cell membrane. A synaptic profile is extracted next to the cell membrane and normalized with reference to the intracellular brightness. Furthermore, a manually selected reference cell is used to normalize settings of the microscope as well as variations in histochemical processing for each stain. Thereafter, staining, homogeneity, and allocation of boutons are determined automatically from the synaptic profiles. The system is evaluated by applying the coefficient of variation (C/sub v/) to repeated measurements of a quantity. Based on 1856 motoneuronal images acquired from four animals with three stains, 93% of the images are analyzed correctly. The others were rejected, based on process protocols. Using only rabbit anti-synaptophysin as primary antibody, the correctness increases above 96%, C/sub v/ values are below 3%, 5%, and 6% for all measures with respect to stochastic optimization, cell positioning, and a large range of microscope settings, respectively. A sample size of about 100 is required to validate a significant reduction of staining in motoneurons below a hemi-section (Wilcoxon rank-sum test, /spl alpha/=0.05, /spl beta/=0.9). The authors' system yields statistically robust results from light micrographs. In future, it is hoped that this system will substitute for the expensive and time-consuming analysis of spinal cord injury at the ultra-structural level, such as by manual interpretation of nonoverlapping electron micrographs.

[1]  N. Otsu A threshold selection method from gray level histograms , 1979 .

[2]  J. Rafols,et al.  The ultrastructure and synaptic architecture of phrenic motor neurons in the spinal cord of the adult rat , 1984, Journal of neurocytology.

[3]  I. T. Young,et al.  Quantitative Microscopy , 1984, Definitions.

[4]  K Mottaghy,et al.  A novel method for quantifying shape deformation applied to biocompatibility testing. , 1999, ASAIO journal.

[5]  Daniel Rueckert,et al.  Automatic tracking of the aorta in cardiovascular MR images using deformable models , 1997, IEEE Transactions on Medical Imaging.

[6]  Kenong Wu,et al.  Live cell image segmentation , 1995, IEEE Transactions on Biomedical Engineering.

[7]  Laurent D. Cohen,et al.  Finite-Element Methods for Active Contour Models and Balloons for 2-D and 3-D Images , 1993, IEEE Trans. Pattern Anal. Mach. Intell..

[8]  J. Noth,et al.  Structural changes of anterior horn neurons and their synaptic input caudal to a low thoracic spinal cord hemisection in the adult rat: a light and electron microscopic study , 2004, Acta Neuropathologica.

[9]  Xin Zheng,et al.  A particle segmentation method based on nucleation and growth of the background , 1994 .

[10]  G. Horcholle-Bossavit,et al.  Distribution of glycinergic terminals on lumbar motoneurons of the adult cat: an ultrastructural study , 1992, Brain Research.

[11]  Q Tai,et al.  Ultrastructural quantitative analysis of glutamatergic and GABAergic synaptic terminals in the phrenic nucleus after spinal cord injury , 1996, The Journal of comparative neurology.

[12]  D A McCrea,et al.  Effects of chronic spinalization on ankle extensor motoneurons. I. Composite monosynaptic Ia EPSPs in four motoneuron pools. , 1994, Journal of neurophysiology.

[13]  J.-M. Chassery,et al.  A new approach to man-machine communication for computerized microscopy , 1994, IEEE Transactions on Biomedical Engineering.

[14]  R. Bloch,et al.  Statistical and analytical techniques , 1984 .

[15]  M. E. Goldberger,et al.  The Recovery of Postural Reflexes and Locomotion Following Low Thoracic Hemisection in Adult Cats Involves Compensation by Undamaged Primary Afferent Pathways , 1993, Experimental Neurology.

[16]  Geir Storvik,et al.  A Bayesian Approach to Dynamic Contours Through Stochastic Sampling and Simulated Annealing , 1994, IEEE Trans. Pattern Anal. Mach. Intell..

[17]  Petros Maragos,et al.  Pattern Spectrum and Multiscale Shape Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[18]  J. Noth,et al.  Trends in the pathophysiology and pharmacotherapy of spasticity , 1991, Journal of Neurology.

[19]  Thomas Martin Deserno,et al.  Automatic parameter setting for balloon models , 2000, Medical Imaging: Image Processing.

[20]  Demetri Terzopoulos,et al.  Deformable models in medical image analysis: a survey , 1996, Medical Image Anal..

[21]  Thomas Martin Deserno,et al.  Scale-independent shape analysis for quantitative cytology using mathematical morphology , 2000, Comput. Biol. Medicine.

[22]  Xinhua Zhuang,et al.  Image Analysis Using Mathematical Morphology , 1987, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[23]  Xiao-Jun Yu,et al.  Aging Enhances Synaptic Efficacy in a Latent Motor Pathway Following Spinal Cord Hemisection in Adult Rats , 1993, Experimental Neurology.

[24]  Gene H. Golub,et al.  Matrix computations , 1983 .

[25]  Thomas Martin Deserno,et al.  Survey: interpolation methods in medical image processing , 1999, IEEE Transactions on Medical Imaging.

[26]  Thomas Martin Deserno,et al.  General finite element model for segmentation in 2, 3, and 4 dimensions , 2000, Medical Imaging: Image Processing.

[27]  S. Bramble Image analysis for the biological sciences , 1996 .

[28]  B. Macq,et al.  Morphological feature extraction for the classification of digital images of cancerous tissues , 1996, IEEE Transactions on Biomedical Engineering.

[29]  J. Barba,et al.  A parametric fitting algorithm for segmentation of cell images , 1998, IEEE Transactions on Biomedical Engineering.

[30]  Thomas Martin Deserno,et al.  Deformable membrane for the segmentation of cytological samples , 1998, Medical Imaging.

[31]  Jayasooriah,et al.  Image analysis of tissue sections , 1996, Comput. Biol. Medicine.

[32]  Hai-Shan Wu,et al.  An efficient semi‐automatic algorithm for cell contour extraction , 1995 .

[33]  Vassili A. Kovalev,et al.  Segmentation technique of complex image scene for an automatic blood-cell-counting system , 1996, Medical Imaging.

[34]  A Elmoataz,et al.  Generalized region growing operator with optimal scanning: application to segmentation of breast cancer images , 1997, Journal of microscopy.