Automated Cell Detection and Morphometry on Growth Plate Images of Mouse Bone.

Microscopy imaging of mouse growth plates is extensively used in biology to understand the effect of specific molecules on various stages of normal bone development and on bone disease. Until now, such image analysis has been conducted by manual detection. In fact, when existing automated detection techniques were applied, morphological variations across the growth plate and heterogeneity of image background color, including the faint presence of cells (chondrocytes) located deeper in tissue away from the image's plane of focus, and lack of cell-specific features, interfered with identification of cell. We propose the first method of automated detection and morphometry applicable to images of cells in the growth plate of long bone. Through ad hoc sequential application of the Retinex method, anisotropic diffusion and thresholding, our new cell detection algorithm (CDA) addresses these challenges on bright-field microscopy images of mouse growth plates. Five parameters, chosen by the user in respect of image characteristics, regulate our CDA. Our results demonstrate effectiveness of the proposed numerical method relative to manual methods. Our CDA confirms previously established results regarding chondrocytes' number, area, orientation, height and shape of normal growth plates. Our CDA also confirms differences previously found between the genetic mutated mouse Smad1/5CKO and its control mouse on fluorescence images. The CDA aims to aid biomedical research by increasing efficiency and consistency of data collection regarding arrangement and characteristics of chondrocytes. Our results suggest that automated extraction of data from microscopy imaging of growth plates can assist in unlocking information on normal and pathological development, key to the underlying biological mechanisms of bone growth.

[1]  Zia-ur Rahman,et al.  Retinex processing for automatic image enhancement , 2004, J. Electronic Imaging.

[2]  Anne E Carpenter,et al.  CellProfiler: image analysis software for identifying and quantifying cell phenotypes , 2006, Genome Biology.

[3]  Dennis Gabor,et al.  Theory of communication , 1946 .

[4]  M. Krein,et al.  On Regularly Convex Sets in the Space Conjugate to a Banach Space , 1940 .

[5]  Ryan Wilson,et al.  Effect of localization, length and orientation of chondrocytic primary cilium on murine growth plate organization. , 2011, Journal of theoretical biology.

[6]  Iman M Helmy,et al.  Efficacy of ImageJ in the assessment of apoptosis , 2012, Diagnostic Pathology.

[7]  Jean-Michel Morel,et al.  A PDE Formalization of Retinex Theory , 2010, IEEE Transactions on Image Processing.

[8]  Gjlles Aubert,et al.  Mathematical problems in image processing , 2001 .

[9]  Robert M. Haralick,et al.  Textural Features for Image Classification , 1973, IEEE Trans. Syst. Man Cybern..

[10]  Kevin W Eliceiri,et al.  NIH Image to ImageJ: 25 years of image analysis , 2012, Nature Methods.

[11]  E. Land,et al.  Lightness and retinex theory. , 1971, Journal of the Optical Society of America.

[12]  G. Karsenty,et al.  Genetic control of bone formation. , 2009, Annual review of cell and developmental biology.

[13]  L. R. Dice Measures of the Amount of Ecologic Association Between Species , 1945 .

[14]  B. Yoon,et al.  Connective tissue growth factor coordinates chondrogenesis and angiogenesis during skeletal development , 2003, Development.

[15]  M. Ascenzi,et al.  Analysis of the orientation of primary cilia in growth plate cartilage: a mathematical method based on multiphoton microscopical images. , 2007, Journal of structural biology.

[16]  H. J. Arnold Introduction to the Practice of Statistics , 1990 .

[17]  Michael Elad,et al.  A Variational Framework for Retinex , 2002, IS&T/SPIE Electronic Imaging.

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

[19]  Wotao Yin,et al.  Image Cartoon-Texture Decomposition and Feature Selection Using the Total Variation Regularized L1 Functional , 2005, VLSM.

[20]  Luiz Velho,et al.  Image moments-based structuring and tracking of objects , 2002, Proceedings. XV Brazilian Symposium on Computer Graphics and Image Processing.

[21]  Yves Meyer,et al.  Oscillating Patterns in Image Processing and Nonlinear Evolution Equations: The Fifteenth Dean Jacqueline B. Lewis Memorial Lectures , 2001 .

[22]  D. Ornitz,et al.  Development of the endochondral skeleton. , 2013, Cold Spring Harbor perspectives in biology.

[23]  C. Farnum,et al.  Visualization of living terminal hypertrophic chondrocytes of growth plate cartilage in situ by differential interference contrast microscopy and time‐lapse cinematography , 1990, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[24]  Nikos Paragios,et al.  Variational, Geometric, and Level Set Methods in Computer Vision, Third International Workshop, VLSM 2005, Beijing, China, October 16, 2005, Proceedings , 2005, VLSM.

[25]  Fabian J. Theis,et al.  An automatic method for robust and fast cell detection in bright field images from high-throughput microscopy , 2013, BMC Bioinformatics.

[26]  Joachim Weickert,et al.  Anisotropic diffusion in image processing , 1996 .

[27]  Yongjin Sung,et al.  Multiple Phases of Chondrocyte Enlargement Underlie Differences in Skeletal Proportions , 2013, Nature.

[28]  Sankar K. Pal,et al.  A review on image segmentation techniques , 1993, Pattern Recognit..

[29]  Stanley Osher,et al.  Modeling Textures with Total Variation Minimization and Oscillating Patterns in Image Processing , 2003, J. Sci. Comput..

[30]  Tony J Collins,et al.  ImageJ for microscopy. , 2007, BioTechniques.

[31]  K. Retting,et al.  Smad6 is essential to limit BMP signaling during cartilage development , 2011, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[32]  Stanley Osher,et al.  A Low Patch-Rank Interpretation of Texture , 2013, SIAM J. Imaging Sci..

[33]  K. Retting,et al.  BMP canonical Smad signaling through Smad1 and Smad5 is required for endochondral bone formation , 2009, Development.

[34]  H. Kronenberg,et al.  Developmental regulation of the growth plate , 2003, Nature.

[35]  Pierre Kornprobst,et al.  Mathematical problems in image processing - partial differential equations and the calculus of variations , 2010, Applied mathematical sciences.

[36]  S. P. Lloyd,et al.  Least squares quantization in PCM , 1982, IEEE Trans. Inf. Theory.

[37]  I. Shapiro,et al.  Fate of the hypertrophic chondrocyte: microenvironmental perspectives on apoptosis and survival in the epiphyseal growth plate. , 2005, Birth defects research. Part C, Embryo today : reviews.

[38]  D Fenistein,et al.  A fast, fully automated cell segmentation algorithm for high‐throughput and high‐content screening , 2008, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[39]  B. Song,et al.  Development of the post-natal growth plate requires intraflagellar transport proteins. , 2007, Developmental biology.