ShapeMetrics: A userfriendly pipeline for 3D cell segmentation and spatial tissue analysis

ShapeMetrics: A 3D Cell Segmentation Pipeline for Single-Cell Spatial Morphometric Analysis

There is a growing need for single-cell level data analysis in correlation with the advancements of microscopy techniques. Morphology-based statistics gathered from individual cells are essential for detection and quantification of even subtle changes within the complex tissues, yet the information available from high-resolution imaging is oftentimes sub-optimally utilized due to the lack of proper computational analysis software. Here we present ShapeMetrics, a 3D cell segmentation pipeline that we have developed to identify, analyze, and quantify single cells in an image. This MATLAB-based script enables users to extract morphological parameters, such as ellipticity, longest axis, cell elongation, or the ratio between cell volume and surface area. We have specifically invested in creating a user-friendly pipeline, aimed for biologists with a limited computational background. Our pipeline is presented with detailed stepwise instructions, starting from the establishment of machine learning-based prediction files of immuno-labeled cell membranes followed by the application of 3D cell segmentation and parameter extraction script, leading to the morphometric analysis and spatial visualization of cell clusters defined by their morphometric features.

Fusome topology and inheritance during insect gametogenesis

From insects to mammals, oocytes and sperm develop within germline cysts comprising cells connected by intercellular bridges (ICBs). In numerous insects, formation of the cyst is accompanied by growth of the fusome-a membranous organelle that permeates the cyst. Fusome composition and function are best understood in Drosophila melanogaster: during oogenesis, the fusome dictates cyst topology and size and facilitates oocyte selection, while during spermatogenesis, the fusome synchronizes the cyst's response to DNA damage. Despite its distinct and sex-specific roles during insect gametogenesis, elucidating fusome growth and inheritance in females and its structure and connectivity in males has remained challenging. Here, we take advantage of advances in three-dimensional (3D) confocal microscopy and computational image processing tools to reconstruct the topology, growth, and distribution of the fusome in both sexes. In females, our experimental findings inform a theoretical model for fusome assembly and inheritance and suggest that oocyte selection proceeds through an 'equivalency with a bias' mechanism. In males, we find that cell divisions can deviate from the maximally branched pattern observed in females, leading to greater topological variability. Our work consolidates existing disjointed experimental observations and contributes a readily generalizable computational approach for quantitative studies of gametogenesis within and across species.

Novel image analysis tool for rapid screening of cell morphology in preclinical animal models of disease

The field of cell biology has seen major advances in both cellular imaging modalities and the development of automated image analysis platforms that increase rigor, reproducibility, and throughput for large imaging data sets. However, there remains a need for tools that provide accurate morphometric analysis of single cells with complex, dynamic cytoarchitecture in a high-throughput and unbiased manner. We developed a fully automated image-analysis algorithm to rapidly detect and quantify changes in cellular morphology using microglia cells, an innate immune cell within the central nervous system, as representative of cells that exhibit dynamic and complex cytoarchitectural changes. We used two preclinical animal models that exhibit robust changes in microglia morphology: (1) a rat model of acute organophosphate intoxication, which was used to generate fluorescently labeled images for algorithm development; and (2) a rat model of traumatic brain injury, which was used to validate the algorithm using cells labeled using chromogenic detection methods. All ex vivo brain sections were immunolabeled for IBA-1 using fluorescence or diaminobenzidine (DAB) labeling, images were acquired using a high content imaging system and analyzed using a custom-built algorithm. The exploratory data set revealed eight statistically significant and quantitative morphometric parameters that distinguished between phenotypically distinct groups of microglia. Manual validation of single-cell morphology was strongly correlated with the automated analysis and was further supported by a comparison with traditional stereology methods. Existing image analysis pipelines rely on high-resolution images of individual cells, which limits sample size and is subject to selection bias. However, our fully automated method integrates quantification of morphology and fluorescent/chromogenic signals in images from multiple brain regions acquired using high-content imaging. In summary, our free, customizable image analysis tool provides a high-throughput, unbiased method for accurately detecting and quantifying morphological changes in cells with complex morphologies.

Usability of deep learning pipelines for 3D nuclei identification with Stardist and Cellpose

Segmentation of 3D images to identify cells and their molecular outputs can be difficult and tedious. Machine learning algorithms provide a promising alternative to manual analysis as emerging 3D image processing technology can save considerable time. For those unfamiliar with machine learning or 3D image analysis, the rapid advancement of the field can make navigating the newest software options confusing. In this paper, two open-source machine learning algorithms, Cellpose and Stardist, are compared in their application on a 3D light sheet dataset counting fluorescently stained proliferative cell nuclei. The effects of image tiling and background subtraction are shown through image analysis pipelines for both algorithms. Based on our analysis, the relative ease of use of Cellpose and the absence of need to train a model leaves it a strong option for 3D cell segmentation despite relatively longer processing times. When Cellpose's pretrained model yields results that are not of sufficient quality, or the analysis of a large dataset is required, Stardist may be more appropriate. Despite the time it takes to train the model, Stardist can create a model specialized to the users' dataset that can be iteratively improved until predictions are satisfactory with far lower processing time relative to other methods.

3D cell morphology detection by association for embryo heart morphogenesis

Advances in tissue engineering for cardiac regenerative medicine require cellular-level understanding of the mechanism of cardiac muscle growth during embryonic developmental stage. Computational methods to automatize cell segmentation in 3D and deliver accurate, quantitative morphology of cardiomyocytes, are imperative to provide insight into cell behavior underlying cardiac tissue growth. Detecting individual cells from volumetric images of dense tissue, poised with low signal-to-noise ratio and severe intensity in homogeneity, is a challenging task. In this article, we develop a robust segmentation tool capable of extracting cellular morphological parameters from 3D multifluorescence images of murine heart, captured via light-sheet microscopy. The proposed pipeline incorporates a neural network for 2D detection of nuclei and cell membranes. A graph-based global association employs the 2D nuclei detections to reconstruct 3D nuclei. A novel optimization embedding the network flow algorithm in an alternating direction method of multipliers is proposed to solve the global object association problem. The associated 3D nuclei serve as the initialization of an active mesh model to obtain the 3D segmentation of individual myocardial cells. The efficiency of our method over the state-of-the-art methods is observed via various qualitative and quantitative evaluation.

A novel deep learning-based 3D cell segmentation framework for future image-based disease detection

Cell segmentation plays a crucial role in understanding, diagnosing, and treating diseases. Despite the recent success of deep learning-based cell segmentation methods, it remains challenging to accurately segment densely packed cells in 3D cell membrane images. Existing approaches also require fine-tuning multiple manually selected hyperparameters on the new datasets. We develop a deep learning-based 3D cell segmentation pipeline, 3DCellSeg, to address these challenges. Compared to the existing methods, our approach carries the following novelties: (1) a robust two-stage pipeline, requiring only one hyperparameter; (2) a light-weight deep convolutional neural network (3DCellSegNet) to efficiently output voxel-wise masks; (3) a custom loss function (3DCellSeg Loss) to tackle the clumped cell problem; and (4) an efficient touching area-based clustering algorithm (TASCAN) to separate 3D cells from the foreground masks. Cell segmentation experiments conducted on four different cell datasets show that 3DCellSeg outperforms the baseline models on the ATAS (plant), HMS (animal), and LRP (plant) datasets with an overall accuracy of 95.6%, 76.4%, and 74.7%, respectively, while achieving an accuracy comparable to the baselines on the Ovules (plant) dataset with an overall accuracy of 82.2%. Ablation studies show that the individual improvements in accuracy is attributable to 3DCellSegNet, 3DCellSeg Loss, and TASCAN, with the 3DCellSeg demonstrating robustness across different datasets and cell shapes. Our results suggest that 3DCellSeg can serve a powerful biomedical and clinical tool, such as histo-pathological image analysis, for cancer diagnosis and grading.

Origami: Single-cell 3D shape dynamics oriented along the apico-basal axis of folding epithelia from fluorescence microscopy data

A common feature of morphogenesis is the formation of three-dimensional structures from the folding of two-dimensional epithelial sheets, aided by cell shape changes at the cellular-level. Changes in cell shape must be studied in the context of cell-polarised biomechanical processes within the epithelial sheet. In epithelia with highly curved surfaces, finding single-cell alignment along a biological axis can be difficult to automate in silico. We present 'Origami', a MATLAB-based image analysis pipeline to compute direction-variant cell shape features along the epithelial apico-basal axis. Our automated method accurately computed direction vectors denoting the apico-basal axis in regions with opposing curvature in synthetic epithelia and fluorescence images of zebrafish embryos. As proof of concept, we identified different cell shape signatures in the developing zebrafish inner ear, where the epithelium deforms in opposite orientations to form different structures. Origami is designed to be user-friendly and is generally applicable to fluorescence images of curved epithelia.

Red blood cell phenotyping from 3D confocal images using artificial neural networks

The investigation of cell shapes mostly relies on the manual classification of 2D images, causing a subjective and time consuming evaluation based on a portion of the cell surface. We present a dual-stage neural network architecture for analyzing fine shape details from confocal microscopy recordings in 3D. The system, tested on red blood cells, uses training data from both healthy donors and patients with a congenital blood disease, namely hereditary spherocytosis. Characteristic shape features are revealed from the spherical harmonics spectrum of each cell and are automatically processed to create a reproducible and unbiased shape recognition and classification. The results show the relation between the particular genetic mutation causing the disease and the shape profile. With the obtained 3D phenotypes, we suggest our method for diagnostics and theragnostics of blood diseases. Besides the application employed in this study, our algorithms can be easily adapted for the 3D shape phenotyping of other cell types and extend their use to other applications, such as industrial automated 3D quality control.