Population-scale three-dimensional reconstruction and quantitative profiling of microglia arbors - PubMed (original) (raw)
Population-scale three-dimensional reconstruction and quantitative profiling of microglia arbors
Murad Megjhani et al. Bioinformatics. 2015.
Abstract
Motivation: The arbor morphologies of brain microglia are important indicators of cell activation. This article fills the need for accurate, robust, adaptive and scalable methods for reconstructing 3-D microglial arbors and quantitatively mapping microglia activation states over extended brain tissue regions.
Results: Thick rat brain sections (100-300 µm) were multiplex immunolabeled for IBA1 and Hoechst, and imaged by step-and-image confocal microscopy with automated 3-D image mosaicing, producing seamless images of extended brain regions (e.g. 5903 × 9874 × 229 voxels). An over-complete dictionary-based model was learned for the image-specific local structure of microglial processes. The microglial arbors were reconstructed seamlessly using an automated and scalable algorithm that exploits microglia-specific constraints. This method detected 80.1 and 92.8% more centered arbor points, and 53.5 and 55.5% fewer spurious points than existing vesselness and LoG-based methods, respectively, and the traces were 13.1 and 15.5% more accurate based on the DIADEM metric. The arbor morphologies were quantified using Scorcioni's L-measure. Coifman's harmonic co-clustering revealed four morphologically distinct classes that concord with known microglia activation patterns. This enabled us to map spatial distributions of microglial activation and cell abundances.
Availability and implementation: Experimental protocols, sample datasets, scalable open-source multi-threaded software implementation (C++, MATLAB) in the electronic supplement, and website (www.farsight-toolkit.org). http://www.farsight-toolkit.org/wiki/Population-scale\_Three-dimensional\_Reconstruction\_and\_Quanti-tative\_Profiling\_of\_Microglia\_Arbors
Contact: broysam@central.uh.edu
Supplementary information: Supplementary data are available at Bioinformatics online.
© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.
Figures
Fig. 1.
Illustrating the problem of large-scale microglia reconstruction. (A) Maximum-intensity projection of a 3-D two-channel mosaic image (size 5903 × 9874 × 229 voxels, 3117 microglia, yellow: IBA1, blue: Hoechst). (A1) close-up of the boxed region. (B) Automated reconstructions with microglial soma in red and arbor reconstructions in green. (B1) Close-up of the reconstructions in the boxed region. (C1–C5) Close-ups of individual microglial arbor reconstructions illustrating the diversity of cellular morphologies
Fig. 2.
A subset of the 3-D image patch dictionary entries displayed as maximum-intensity projections. It is clear from this illustration that some patches (highlighted by red boxes) clearly belong to parts of microglial arbors, whereas the other patches are either ambiguous or part of the background
Fig. 3.
Illustrating 3-D seed points used for arbor reconstruction. (A) Sample IBA1+ cell (grayscale) overlaid with initially detected seed points (red). (B) Close-up of boxed region in Panel A. The white arrows indicate some valid seed points that lie on the cellular arbor, and the yellow arrows some indicate invalid seed points. Only a subset of the seed points are marked by arrows here to avoid clutter
Fig. 4.
Performance of our algorithm compared to manual reconstruction as measured by the DIADEM metric. (A) The mean DIADEM metric for our algorithm (0.77) was higher than the value for the algorithm of (Galbreath, 2011)[1](0.65) and the GVF based Vesselness-(Wang et al., 2011)[2] (0.68) based algorithm. (B) Our algorithm was insensitive to the size of the dictionary parameter K. (C) Our algorithm was insensitive to variations in the sparsity constraint T. (D) Increasing the cost threshold τ led to reduced performance, as expected. A threshold of 200 was chosen empirically for all the subsequent experiments based on this analysis
Fig. 5.
Quantitative analysis of the microglial cell population by harmonic co-clustering of the L-measure data. (A) 3-D rendering of the reconstructed microglial field in Figure 1, with the arbor traces color coded to indicate the automatically identified groups (red: Group 1, Yellow: Group 2, Green: Group 3, Blue: Group 4). (B) Heatmap representation of the co-clustering results. Each row in the heatmap corresponds to one cell, and each column corresponds to one feature in the L-measure. The co-clustering identified four principal groups. Representative cells from these groups are shown in Panels G1–G4, and correspond to ramified cells (G1), elongated cells (G2), activated cells (G3) and amoeboid/round cells (G4), respectively
Similar articles
- Quantitative arbor analytics: unsupervised harmonic co-clustering of populations of brain cell arbors based on L-measure.
Lu Y, Carin L, Coifman R, Shain W, Roysam B. Lu Y, et al. Neuroinformatics. 2015 Jan;13(1):47-63. doi: 10.1007/s12021-014-9237-2. Neuroinformatics. 2015. PMID: 25086878 - Quantitative 3-D analysis of GFAP labeled astrocytes from fluorescence confocal images.
Kulkarni PM, Barton E, Savelonas M, Padmanabhan R, Lu Y, Trett K, Shain W, Leasure JL, Roysam B. Kulkarni PM, et al. J Neurosci Methods. 2015 May 15;246:38-51. doi: 10.1016/j.jneumeth.2015.02.014. Epub 2015 Mar 5. J Neurosci Methods. 2015. PMID: 25745860 - Morphologically constrained spectral unmixing by dictionary learning for multiplex fluorescence microscopy.
Megjhani M, Correa de Sampaio P, Leigh Carstens J, Kalluri R, Roysam B. Megjhani M, et al. Bioinformatics. 2017 Jul 15;33(14):2182-2190. doi: 10.1093/bioinformatics/btx108. Bioinformatics. 2017. PMID: 28334208 Free PMC article. - The DIADEM data sets: representative light microscopy images of neuronal morphology to advance automation of digital reconstructions.
Brown KM, Barrionuevo G, Canty AJ, De Paola V, Hirsch JA, Jefferis GS, Lu J, Snippe M, Sugihara I, Ascoli GA. Brown KM, et al. Neuroinformatics. 2011 Sep;9(2-3):143-57. doi: 10.1007/s12021-010-9095-5. Neuroinformatics. 2011. PMID: 21249531 Free PMC article. Review. - Stereology, morphometry, and mapping: the whole is greater than the sum of its parts.
Glaser JR, Glaser EM. Glaser JR, et al. J Chem Neuroanat. 2000 Oct;20(1):115-26. doi: 10.1016/s0891-0618(00)00073-9. J Chem Neuroanat. 2000. PMID: 11074348 Review.
Cited by
- Preserving Retinal Structure and Function with the Novel Nitroxide Antioxidant, DCTEIO.
Rayner CL, Bottle SE, Martyn AP, Barnett NL. Rayner CL, et al. Neurochem Res. 2023 Nov;48(11):3402-3419. doi: 10.1007/s11064-023-03978-w. Epub 2023 Jul 14. Neurochem Res. 2023. PMID: 37450210 Free PMC article. - More than meets the eye: The role of microglia in healthy and diseased retina.
Murenu E, Gerhardt MJ, Biel M, Michalakis S. Murenu E, et al. Front Immunol. 2022 Nov 29;13:1006897. doi: 10.3389/fimmu.2022.1006897. eCollection 2022. Front Immunol. 2022. PMID: 36524119 Free PMC article. Review. - Assaying Microglia Functions In Vitro.
Maguire E, Connor-Robson N, Shaw B, O'Donoghue R, Stöberl N, Hall-Roberts H. Maguire E, et al. Cells. 2022 Oct 28;11(21):3414. doi: 10.3390/cells11213414. Cells. 2022. PMID: 36359810 Free PMC article. Review. - Cell morphologies in the nervous system: Glia steal the limelight.
Ascoli GA. Ascoli GA. J Comp Neurol. 2023 Feb;531(3):338-343. doi: 10.1002/cne.25429. Epub 2022 Oct 31. J Comp Neurol. 2023. PMID: 36316800 Free PMC article. - Machine learning classification reveals robust morphometric biomarker of glial and neuronal arbors.
Akram MA, Wei Q, Ascoli GA. Akram MA, et al. J Neurosci Res. 2023 Jan;101(1):112-129. doi: 10.1002/jnr.25131. Epub 2022 Oct 5. J Neurosci Res. 2023. PMID: 36196621 Free PMC article.
References
- Aharon M., et al. . (2006) K-SVD: an algorithm for designing overcomplete dictionaries for sparse representation. IEEE Trans. Signal Process., 54, 4311–4322.
- Al-Kofahi Y., et al. . (2010) Improved automatic detection and segmentation of cell nuclei in histopathology images. IEEE Trans. Biomed. Eng., 57, 841–852. - PubMed
- Al-Kofahi Y., et al. . (2008) Improved detection of branching points in algorithms for automated neuron tracing from 3D confocal images. Cytom. Part A, 73, 36–43. - PubMed
- Bas E., Erdogmus D. (2010) Piecewise linear cylinder models for 3-dimensional axon segmentation in brainbow imagery. In: Proceedings of the IEEE International Symposium on Biomedical Imaging: From Nano to Macro., pp. 1297–1300.
Publication types
MeSH terms
LinkOut - more resources
Full Text Sources
Other Literature Sources
Miscellaneous