Dynamic spherical harmonics approach for shape classification of migrating cells

Epithelial uptake leads to fungal killing in vivo and is aberrant in COPD-derived epithelial cells

Hundreds of spores of Aspergillus fumigatus (Af) are inhaled daily by human beings, representing a constant, possibly fatal, threat to respiratory health. The small size of Af spores suggests that interactions with alveolar epithelial cells (AECs) are frequent; thus, we hypothesized that spore uptake by AECs is important for driving fungal killing and susceptibility to Aspergillus-related disease. Using single-cell approaches to measure spore uptake and its outcomes in vivo, we demonstrate that Af spores are internalized and killed by AECs during whole-animal infection. Moreover, comparative analysis of primary human AECs from healthy and chronic obstructive pulmonary disease (COPD) donors revealed significant alterations in the uptake and killing of spores in COPD-derived AECs. We conclude that AECs contribute to the killing of Af spores and that dysregulation of curative AEC responses in COPD may represent a driver of Aspergillus-related diseases.

CAJAL enables analysis and integration of single-cell morphological data using metric geometry

High-resolution imaging has revolutionized the study of single cells in their spatial context. However, summarizing the great diversity of complex cell shapes found in tissues and inferring associations with other single-cell data remains a challenge. Here, we present CAJAL, a general computational framework for the analysis and integration of single-cell morphological data. By building upon metric geometry, CAJAL infers cell morphology latent spaces where distances between points indicate the amount of physical deformation required to change the morphology of one cell into that of another. We show that cell morphology spaces facilitate the integration of single-cell morphological data across technologies and the inference of relations with other data, such as single-cell transcriptomic data. We demonstrate the utility of CAJAL with several morphological datasets of neurons and glia and identify genes associated with neuronal plasticity in C. elegans. Our approach provides an effective strategy for integrating cell morphology data into single-cell omics analyses.

4D reconstruction of murine developmental trajectories using spherical harmonics

Normal organogenesis cannot be recapitulated in vitro for mammalian organs, unlike in species including Drosophila and zebrafish. Available 3D data in the form of ex vivo images only provide discrete snapshots of the development of an organ morphology. Here, we propose a computer-based approach to recreate its continuous evolution in time and space from a set of 3D volumetric images. Our method is based on the remapping of shape data into the space of the coefficients of a spherical harmonics expansion where a smooth interpolation over time is simpler. We tested our approach on mouse limb buds and embryonic hearts. A key advantage of this method is that the resulting 4D trajectory can take advantage of all the available data while also being able to interpolate well through time intervals for which there are little or no data. This allows for a quantitative, data-driven 4D description of mouse limb morphogenesis.

•

Computer-based method recreating a 3D plus time evolution of a set of volumetric images

•

Technique based on the interpolation of the coefficients of spherical harmonics

•

Data-driven quantitative 4D description of limb and heart morphogenesis

•

Quantitatively reliable baseline description of organ development

T cell morphodynamics reveal periodic shape oscillations in three-dimensional migration

T cells use sophisticated shape dynamics (morphodynamics) to migrate towards and neutralize infected and cancerous cells. However, there is limited quantitative understanding of the migration process in three-dimensional extracellular matrices (ECMs) and across timescales. Here, we leveraged recent advances in lattice light-sheet microscopy to quantitatively explore the three-dimensional morphodynamics of migrating T cells at high spatio-temporal resolution. We first developed a new shape descriptor based on spherical harmonics, incorporating key polarization information of the uropod. We found that the shape space of T cells is low-dimensional. At the behavioural level, run-and-stop migration modes emerge at approximately 150 s, and we mapped the morphodynamic composition of each mode using multiscale wavelet analysis, finding 'stereotyped' motifs. Focusing on the run mode, we found morphodynamics oscillating periodically (every approx. 100 s) that can be broken down into a biphasic process: front-widening with retraction of the uropod, followed by a rearward surface motion and forward extension, where intercalation with the ECM in both of these steps likely facilitates forward motion. Further application of these methods may enable the comparison of T cell migration across different conditions (e.g. differentiation, activation, tissues and drug treatments) and improve the precision of immunotherapeutic development.

High-Throughput 3D Phenotyping of Plant Shoot Apical Meristems From Tissue-Resolution Data

Confocal imaging is a well-established method for investigating plant phenotypes on the tissue and organ level. However, many differences are difficult to assess by visual inspection and researchers rely extensively on ad hoc manual quantification techniques and qualitative assessment. Here we present a method for quantitatively phenotyping large samples of plant tissue morphologies using triangulated isosurfaces. We successfully demonstrate the applicability of the approach using confocal imaging of aerial organs in Arabidopsis thaliana. Automatic identification of flower primordia using the surface curvature as an indication of outgrowth allows for high-throughput quantification of divergence angles and further analysis of individual flowers. We demonstrate the throughput of our method by quantifying geometric features of 1065 flower primordia from 172 plants, comparing auxin transport mutants to wild type. Additionally, we find that a paraboloid provides a simple geometric parameterisation of the shoot inflorescence domain with few parameters. We utilise parameterisation methods to provide a computational comparison of the shoot apex defined by a fluorescent reporter of the central zone marker gene CLAVATA3 with the apex defined by the paraboloid. Finally, we analyse the impact of mutations which alter mechanical properties on inflorescence dome curvature and compare the results with auxin transport mutants. Our results suggest that region-specific expression domains of genes regulating cell wall biosynthesis and local auxin transport can be important in maintaining the wildtype tissue shape. Altogether, our results indicate a general approach to parameterise and quantify plant development in 3D, which is applicable also in cases where data resolution is limited, and cell segmentation not possible. This enables researchers to address fundamental questions of plant development by quantitative phenotyping with high throughput, consistency and reproducibility.

Enhanced segmentation of label-free cells for automated migration and interaction tracking

In biomedical research, the migration behavior of cells and interactions between various cell types are frequently studied subjects. An automated and quantitative analysis of time-lapse microscopy data is an essential component of these studies, especially when characteristic migration patterns need to be identified. Plenty of software tools have been developed to serve this need. However, the majority of algorithms is designed for fluorescently labeled cells, even though it is well-known that fluorescent labels can substantially interfere with the physiological behavior of interacting cells. We here present a fully revised version of our algorithm for migration and interaction tracking (AMIT), which includes a novel segmentation approach. This approach allows segmenting label-free cells with high accuracy and also enables detecting almost all cells within the field of view. With regard to cell tracking, we designed and implemented a new method for cluster detection and splitting. This method does not rely on any geometrical characteristics of individual objects inside a cluster but relies on monitoring the events of cell-cell fusion from and cluster fission into single cells forward and backward in time. We demonstrate that focusing on these events provides accurate splitting of transient clusters. Furthermore, the substantially improved quantitative analysis of cell migration by the revised version of AMIT is more than two orders of magnitude faster than the previous implementation, which makes it feasible to process video data at higher spatial and temporal resolutions.

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.

Parametric-based feature selection via spherical harmonic coefficients for the left ventricle myocardial infarction screening

Computer-aided diagnosis (CAD) of heart diseases using machine learning techniques has recently received much attention. In this study, we present a novel parametric-based feature selection method using the three-dimensional spherical harmonic (SHs) shape descriptors of the left ventricle (LV) for intelligent myocardial infarction (MI) classification. The main hypothesis is that the SH coefficients of the parameterized endocardial shapes in MI patients are recognizable and distinguishable from healthy subjects. The SH parameterization, expansion, and registration of the LV endocardial shapes were performed, then parametric-based features were extracted. The proposed method performance was investigated by varying considered phases (i.e., the end-systole (ES) or the end-diastole (ED) frames), the spatial alignment procedures based on three modes (i.e., the center of the apical (CoA), the center of mass (CoM), and the center of the basal (CoB)), and considered orders of SH coefficients. After applying principal component analysis (PCA) on the feature vectors, support vector machine (SVM), K-nearest neighbors (K-NN), and random forest (RF) were trained and tested using the leave-one-out cross-validation (LOOCV). The proposed method validation was performed via a dataset containing healthy and MI subjects selected from the automated cardiac diagnosis challenge (ACDC) database. The promising results show the effectiveness of the proposed classification model. SVM reached the best performance with accuracy, sensitivity, specificity, and F-score of 97.50%, 95.00%, 100.00%, and 97.56%, respectively, using the introduced optimum feature set. This study demonstrates the robustness of combining the SH coefficients and machine learning techniques. We also quantify and notably highlight the contribution of different parameters in the classification and finally introduce an optimal feature set with maximum discriminant strength for the MI classification task. Moreover, the obtained results confirm that the proposed method performs more accurately than conventional point-based methods and also the current start-of-the-art, i.e., clinical measures. We showed our method's generalizability using employing it in dilated cardiomyopathy (DCM) detection and achieving promising results too. Parametric-based feature selection via spherical harmonics coefficients for the left ventricle myocardial infarction screening.

Valproic acid-induced changes of 4D nuclear morphology in astrocyte cells

Histone deacetylase inhibitors, such as valproic acid (VPA), have important clinical therapeutic and cellular reprogramming applications. They induce chromatin reorganization that is associated with altered cellular morphology. However, there is a lack of comprehensive characterization of VPA-induced changes of nuclear size and shape. Here, we quantify 3D nuclear morphology of primary human astrocyte cells treated with VPA over time (hence, 4D). We compared volumetric and surface-based representations and identified seven features that jointly discriminate between normal and treated cells with 85% accuracy on day 7. From day 3, treated nuclei were more elongated and flattened and then continued to morphologically diverge from controls over time, becoming larger and more irregular. On day 7, most of the size and shape descriptors demonstrated significant differences between treated and untreated cells, including a 24% increase in volume and 6% reduction in extent (shape regularity) for treated nuclei. Overall, we show that 4D morphometry can capture how chromatin reorganization modulates the size and shape of the nucleus over time. These nuclear structural alterations may serve as a biomarker for histone (de-)acetylation events and provide insights into mechanisms of astrocytes-to-neurons reprogramming.

Integrated analysis of cell shape and movement in moving frame

The cell's movement and morphological change are two interrelated cellular processes. An integrated analysis is needed to explore the relationship between them. However, it has been challenging to investigate them as a whole. The cell's trajectory can be described by its speed, curvature, and torsion. On the other hand, the three-dimensional (3D) cell shape can be studied by using a shape descriptor such as spherical harmonic (SH) descriptor, which is an extension of a Fourier transform in 3D space. We propose a novel method using parallel-transport (PT) to integrate these shape-movement data by using moving frames as the 3D-shape coordinate system. This moving frame is purely determined by the velocity vector. On this moving frame, the movement change will influence the coordinate system for shape analysis. By analyzing the change of the SH coefficients over time in the moving frame, we can observe the relationship between shape and movement. We illustrate the application of our approach using simulated and real datasets in this paper.

In Situ Classification of Cell Types in Human Kidney Tissue Using 3D Nuclear Staining

To understand the physiology and pathology of disease, capturing the heterogeneity of cell types within their tissue environment is fundamental. In such an endeavor, the human kidney presents a formidable challenge because its complex organizational structure is tightly linked to key physiological functions. Advances in imaging-based cell classification may be limited by the need to incorporate specific markers that can link classification to function. Multiplex imaging can mitigate these limitations, but requires cumulative incorporation of markers, which may lead to tissue exhaustion. Furthermore, the application of such strategies in large scale 3-dimensional (3D) imaging is challenging. Here, we propose that 3D nuclear signatures from a DNA stain, DAPI, which could be incorporated in most experimental imaging, can be used for classifying cells in intact human kidney tissue. We developed an unsupervised approach that uses 3D tissue cytometry to generate a large training dataset of nuclei images (NephNuc), where each nucleus is associated with a cell type label. We then devised various supervised machine learning approaches for kidney cell classification and demonstrated that a deep learning approach outperforms classical machine learning or shape-based classifiers. Specifically, a custom 3D convolutional neural network (NephNet3D) trained on nuclei image volumes achieved a balanced accuracy of 80.26%. Importantly, integrating NephNet3D classification with tissue cytometry allowed in situ visualization of cell type classifications in kidney tissue. In conclusion, we present a tissue cytometry and deep learning approach for in situ classification of cell types in human kidney tissue using only a DNA stain. This methodology is generalizable to other tissues and has potential advantages on tissue economy and non-exhaustive classification of different cell types.