Population-scale three-dimensional reconstruction and quantitative profiling of microglia arbors.

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.

Quantitative 3-D analysis of GFAP labeled astrocytes from fluorescence confocal images.

BACKGROUND: There is a need for effective computational methods for quantifying the three-dimensional (3-D) spatial distribution, cellular arbor morphologies, and the morphological diversity of brain astrocytes to support quantitative studies of astrocytes in health, injury, and disease. NEW METHOD: Confocal fluorescence microscopy of multiplex-labeled (GFAP, DAPI) brain tissue is used to perform imaging of astrocytes in their tissue context. The proposed computational method identifies the astrocyte cell nuclei, and reconstructs their arbors using a local priority based parallel (LPP) tracing algorithm. Quantitative arbor measurements are extracted using Scorcioni's L-measure, and profiled by unsupervised harmonic co-clustering to reveal the morphological diversity. RESULTS: The proposed method identifies astrocyte nuclei, generates 3-D reconstructions of their arbors, and extracts quantitative arbor measurements, enabling a morphological grouping of the cell population. COMPARISON WITH EXISTING METHODS: Our method enables comprehensive spatial and morphological profiling of astrocyte populations in brain tissue for the first time, and overcomes limitations of prior methods. Visual proofreading of the results indicate a >95% accuracy in identifying astrocyte nuclei. The arbor reconstructions exhibited 3.2% fewer erroneous jumps in tracing, and 17.7% fewer false segments compared to the widely used fast-marching method that resulted in 9% jumps and 20.8% false segments. CONCLUSIONS: The proposed method can be used for large-scale quantitative studies of brain astrocyte distribution and morphology.

Large-scale automated image analysis for computational profiling of brain tissue surrounding implanted neuroprosthetic devices using Python.

In this article, we describe the use of Python for large-scale automated server-based bio-image analysis in FARSIGHT, a free and open-source toolkit of image analysis methods for quantitative studies of complex and dynamic tissue microenvironments imaged by modern optical microscopes, including confocal, multi-spectral, multi-photon, and time-lapse systems. The core FARSIGHT modules for image segmentation, feature extraction, tracking, and machine learning are written in C++, leveraging widely used libraries including ITK, VTK, Boost, and Qt. For solving complex image analysis tasks, these modules must be combined into scripts using Python. As a concrete example, we consider the problem of analyzing 3-D multi-spectral images of brain tissue surrounding implanted neuroprosthetic devices, acquired using high-throughput multi-spectral spinning disk step-and-repeat confocal microscopy. The resulting images typically contain 5 fluorescent channels. Each channel consists of 6000 × 10,000 × 500 voxels with 16 bits/voxel, implying image sizes exceeding 250 GB. These images must be mosaicked, pre-processed to overcome imaging artifacts, and segmented to enable cellular-scale feature extraction. The features are used to identify cell types, and perform large-scale analysis for identifying spatial distributions of specific cell types relative to the device. Python was used to build a server-based script (Dell 910 PowerEdge servers with 4 sockets/server with 10 cores each, 2 threads per core and 1TB of RAM running on Red Hat Enterprise Linux linked to a RAID 5 SAN) capable of routinely handling image datasets at this scale and performing all these processing steps in a collaborative multi-user multi-platform environment. Our Python script enables efficient data storage and movement between computers and storage servers, logs all the processing steps, and performs full multi-threaded execution of all codes, including open and closed-source third party libraries.