Spherical harmonics texture extraction for versatile analysis of biological objects

Abstract

The characterization of phenotypes in cells or organisms from microscopy data largely depends on differences in the spatial distribution of image intensity. Multiple methods exist for quantifying the intensity distribution - or image texture - across objects in natural images. However, many of these texture extraction methods do not directly adapt to 3D microscopy data. Here, we present Spherical Texture extraction, which measures the variance in intensity per angular wavelength by calculating the Spherical Harmonics or Fourier power spectrum of a spherical or circular projection of the angular mean intensity of the object. This method provides a 20-value characterization that quantifies the scale of features in the spherical projection of the intensity distribution, giving a different signal if the intensity is, for example, clustered in parts of the volume or spread across the entire volume. We apply this method to different systems and demonstrate its ability to describe various biological problems through feature extraction. The Spherical Texture extraction characterizes biologically defined gene expression patterns in Drosophila melanogaster embryos, giving a quantitative read-out for pattern formation. Our method can also quantify morphological differences in Caenorhabditis elegans germline nuclei, which lack a predefined pattern. We show that the classification of germline nuclei using their Spherical Texture outperforms a convolutional neural net when training data is limited. Additionally, we use a similar pipeline on 2D cell migration data to extract polarization direction, quantifying the alignment of fluorescent markers to the migration direction. We implemented the Spherical Texture method as a plugin in ilastik, making it easy to install and apply to any segmented 3D or 2D dataset. Additionally, this technique can also easily be applied through a Python package to provide extra feature extraction for any object classification pipeline or downstream analysis. Author summaryWe introduce a novel method to extract quantitative data from microscopy images by precisely measuring the distribution of intensities within objects in both 3D or 2D. This method is easily accessible through the object classification workflow of ilastik, provided the original image is segmented into separate objects. The method is specifically designed to analyze mostly convex objects, focusing on the variation in fluorescence intensity caused by differences in their shapes or patterns. We demonstrate the versatility of our method by applying it to very different biological samples. Specifically, we showcase its effectiveness in quantifying the patterning in D. melanogaster embryos, in classifying the nuclei in C. elegans germlines, and in extracting polarization information from individual migratory cells. Through these examples, we illustrate that our technique can be employed across different biological scales. Furthermore, we highlight the multiple ways in which the data generated by our method can be used, including quantifying the strength of a specific pattern, employing machine learning to classify diverse morphologies, or extracting directionality or polarization from fluorescence intensity.