Abstract
Quantitative analysis of morphological changes in a cell nucleus is important for the understanding of nuclear architecture and its relationship with pathological conditions such as cancer. However, dimensionality of imaging data, together with a great variability of nuclear shapes, presents challenges for 3D morphological analysis. Thus, there is a compelling need for robust 3D nuclear morphometric techniques to carry out population-wide analysis. We propose a new approach that combines modeling, analysis, and interpretation of morphometric characteristics of cell nuclei and nucleoli in 3D. We used robust surface reconstruction that allows accurate approximation of 3D object boundary. Then, we computed geometric morphological measures characterizing the form of cell nuclei and nucleoli. Using these features, we compared over 450 nuclei with about 1,000 nucleoli of epithelial and mesenchymal prostate cancer cells, as well as 1,000 nuclei with over 2,000 nucleoli from serum-starved and proliferating fibroblast cells. Classification of sets of 9 and 15 cells achieved accuracy of 95.4% and 98%, respectively, for prostate cancer cells, and 95% and 98% for fibroblast cells. To our knowledge, this is the first attempt to combine these methods for 3D nuclear shape modeling and morphometry into a highly parallel pipeline workflow for morphometric analysis of thousands of nuclei and nucleoli in 3D.
Similar content being viewed by others

Introduction
Motivation
Cell nuclear morphology is regulated by complex underlying biological mechanisms related to cell differentiation, development, proliferation, and disease1,2,3. Changes in the nuclear form are associated with reorganization of chromatin architecture related to altered functional properties such as gene regulation and expression1,3. Moreover, many studies in mechanobiology show that geometric constraints and mechanical forces applied to a cell deform it and, conversely, affect nuclear and chromatin dynamics, as well as gene and pathway activation4,5. Thus, nuclear morphological quantification becomes of major relevance as studies of the reorganization of the chromatin and DNA architecture in the spatial and temporal framework, known as the 4D nucleome, emerge6,7. Cellular structures of interest in the context of the 4D nucleome include not only the nucleus itself, but also the nucleolus and nucleolar-associating domains, chromosome territories, topologically associating domains, lamina-associating domains, and loop domains in transcription factories6,8. Furthermore, understanding these processes through quantitative analysis of morphological changes also has many medical implications, for example, in detection, understanding, and treatment of pathological conditions such as cancer7,8,9,10.
While efforts have been made to develop cell and nuclear shape characteristics in 2D or pseudo-3D11,12, several studies have demonstrated that 3D morphometric measures provide better results for nuclear shape description and discrimination13,14,15. However, 3D shape descriptors that permit robust morphological analysis and facilitate human interpretation are still under active investigation16. Additionally, the dimensionality and volume of acquired data, various image acquisition conditions, and great variability of cell shapes in a population present challenges for 3D shape analysis methods that should be scalable, robust to noise, and specific enough across cell populations at the same time. Thus, there is a compelling need for robust 3D nuclear morphometric techniques to carry out population-wide analysis17.
3D shape representation and morphometric measures
The way cell nuclear shapes can be measured depends on their representation extracted from image data11. Many 3D morphometric measures are applied âas isâ to 3D geometric objects represented by volumetric data18. However, voxels-based shape representations are noisy, and they may lose fine geometric details or even break the objectâs topological structure. Moreover, these representations are not intrinsic, and vary when changing pose or deforming the object. A recent review of approaches to 3D cell shape description16 separated them into three categories in increasing order of complexity: landmark-based, graph-based, and moment-based. This last category includes approaches that are widely used in cellular morphology and allow the user to obtain a global representation that combines low-order moments describing the coarse conformation with high-order moments retaining information at higher frequency. Typically, before applying these methods, a binary mask or outline of the shape (surface) is first extracted from image data, which is done by most segmentation methods. These masks are assumed to have a sphere-like topology and can be projected onto an appropriate basis. Two popular approaches of this type are spherical harmonics (SPHARM)19 and spherical wavelets20. Both methods first map the surface of interest onto the sphere using appropriate spherical parameterization techniques, and then project it onto a reference function basis living on the sphere. SPHARM is arguably one of the most widely applied cell morphology modeling approaches21,22,23,24. In SPHARM the spherical signal is projected onto a basis of Legendre polynomials, extending the classical Fourier analysis to signals on the two-sphere. SPHARM coefficients describe general conformation of the shape of interest at different spatial scales, are rotation invariant, and can be directly used as features for further analysis25. However, SPHARM methods are most appropriate when low order approximation is satisfactory and become less effective in preserving surface details, as artificial oscillations start to appear when higher order basis functions are incorporated26. More robust smooth surface reconstruction can be obtained from a 3D binary mask via Laplace-Beltrami (LB) eigen-projection, followed by topology-preserving boundary deformation to remove various artifacts26. On a unit sphere, the LB eigen-functions correspond to spherical harmonics, so overall they can be viewed as a generalization of the SPHARM to the complex geometry manifold with local adaptation of the basis to the dataset at hand27. The proposed method has been demonstrated to produce smoother and more detailed surfaces compared to the SPHARM and the topology preserving level sets28. Extracted surfaces are smooth, accurately represent the shape of an object, and can be further used for morphometric analysis.
In order to extract shape geometric characteristics, boundary surfaces of binary masks are typically reconstructed from voxel data and discretized as meshes. At the next step, various useful morphometric descriptors can be computed based on this representation. Useful extrinsic and intrinsic geometric descriptors aim to distinguish between global and local shape features. Intrinsic measures capture shape properties that are invariant under transformations (e.g., affine: rotation, translation and scaling). Various shape morphometry measures, like surface area and Gaussian curvature, represent invariant metrics of complexity, which are stable under special transformations of the surface (e.g., bending) that do not affect the inner geometry of the boundary of the 3D volume29. Alternatively, shape metrics, e.g., mean L2-norm and the extrinsic curvature index, are sensitive to affine transformation and other shape morphology in the ambient space. Shape index and curvedness are morphometric descriptors that can capture local shape features, independently or in relation to the size of an object30. Combination of the object surface reconstruction with the extraction of such shape measures demonstrated high performance in recent neuroimaging studies for discriminatory morphometric analysis of complex 3D shapes of cortical and subcortical brain areas31,32,33.
Technical capabilities and interoperability of tools
When it comes to a choice of tools for 3D cell nuclear morphometrics, reproducibility and implementation availability are among major concerns in the field of bioimage analysis16. To date, many of the widely available software tools for cell shape morphometry were either developed for the analysis of 2D11,34,35,36,37,38 or pseudo-3D images39. Other tools only implement slice-by-slice or voxel-based morphometry40,41,42,43, providing a coarse approximation of the global cell shape that is sensitive to increasing amounts of noise and usually fails to characterize morphological variations occurring at different spatial scales. Other common limitations of many 3D cell morphology solutions include a lack of high-throughput processing capabilities or restrictions to the specific programming language or platform that dictate principles of tool implementation44,45,46. Implementations of methods in a bioimage analysis landscape are highly diverse. They range across programming languages, software libraries, and file formats, which increases module interoperability issues and makes code reuse extremely difficult. Re-implementing underlying methods is often very challenging, time-consuming, and error prone47. Some of the existing bioimage analysis frameworks, including ImageJ48, rely on a plugin architecture, which allows their extension via third-party contributions40,41,43. High-throughput capabilities of some of these software packages are limited to processing of multiple objects simultaneously within its graphical user interface (GUI), for example, Tango43. More advanced packages, such as CellProfiler 2.037, BioimageXD42, and Icy41, provide a basic graphical interface to assemble elementary tasks into reusable pipelines that make it possible to execute in GUI and batch modes. However, these solutions are still limited to specific scripting languages and libraries supported by the main software package. They also donât provide a straightforward way to take advantage of the growing number of parallel hardware configurations, such as clusters, clouds, and high-performance computing, which limits the scalability of these solutions.
An alternative to plugin-based solutions, software platforms with modular design allow integration of already existing solutions into workflows without re-implementing them in a specific language, and provide methods for optimizing module interaction, re-usage, and extension49. An example of an extensive and feature rich solution for building and executing complex workflows is the LONI Pipeline31,50. This client-server platform enables users to efficiently describe atomic modules and end-to-end protocols in a graphical canvas using a large library of powerful computational tools. The Pipeline back-end server has extensive support for parallel execution on a grid cluster, including automated data converting, formatting and transfer, optimal job submission and management, pausing execution, and combining local and remote software and data sources. Most importantly, it makes it very easy to create new custom modules from any software that supports a command line interface (CLI). The Pipeline allows users to take advantage of a highly diverse set of tools and connect them together as steps of a computational protocol that is then executed in a high-throughput, parallel fashion. Validated individual modules and end-to-end workflows may be saved, reused in other workflows, easily modified and repurposed. Additionally, the LONI Pipeline saves information about executed steps (such as software origin, version, and architecture) providing provenance information50,51.
Study aims
This study has two complementary aims. The first aim is to assess and validate 3D morphometry metrics for nuclear and nucleolar shape description and classification. Improving the discriminative performance in terms of statistical metrics has been the driving force behind our methodological efforts and selection of specific tools in this work. First, surfaces of 3D masks extracted from the microscopy data are reconstructed using Laplace-Beltrami eigen-projection and topology-preserving boundary deformation26. Then, we computed intrinsic and extrinsic geometric metrics, which are used as derived signature vectors (shape biomarkers) to characterize the complexity of the 3D shapes and discriminate between observed clinical and phenotypic traits. These metrics include volume, surface area, mean curvature, curvedness, shape index, and fractal dimension30,52,53. Although these methods were previously used in recent neuroimaging studies31,32,33, this is the first attempt, to our knowledge, to apply robust smooth LB-based surface reconstruction with intrinsic and extrinsic morphometric measure extraction to 3D cell nuclear and nucleolar shape modeling and morphometry. Suggested modeling and analysis methods are not restricted to nuclear and nucleolar shapes and can be used for the shape quantification of other cellular compartments, depending on their topology.
The second aim is to develop a reproducible pipeline workflow implementing the entire process that can be customized and expanded for deep exploration of associations between 3D nuclear and nucleolar shape phenotypes in health and disease. High-throughput imaging (HTI) can include automatization of liquid handling, microscopy-based image acquisition, image processing, and statistical data analysis17. Our work focuses on the last two aspects of this definition. We implemented a streamlined multi-step protocol using a diverse set of tools to achieve optimal performance compared to alternatives at each step of analysis. These tools are represented as individual modules seamlessly connected in the LONI Pipeline workflow. This workflow meets modern standards for high-throughput imaging processing and analysis and is mostly automated with a focus on validity and reproducibility. Our implementation is massively parallel, customizable, and provides fully automated execution and data provenance out-of-the-box. At the final step of the workflow, we employed machine learning methods to investigate the associations between cell phenotypes and treatment conditions using cell shape morphometric measures as features. We show that using a combination of 3D nuclear and nucleolar morphometry improves the discrimination between in vitro cell conditions of human fibroblast and human prostate cancer (PC3) cell lines.
To promote the reproducibility of results, facilitate open-scientific development, and enable collaborative validation, we made the pipeline workflows, together with underlying source code, documentation, and derived data from this study, available online54. The workflow will be made available via the LONI Pipeline, along with publicly available computational resources to showcase an online demonstration.
Methods
Figure 1 shows a high-level view of the end-to-end protocol. We start with a dataset of 3D binary nuclear and nucleolar masks. We modeled 3D nuclear and nucleolar boundaries by their surface reconstruction and extracted the derived morphometry measures. Finally, we computed statistical differences, identified shape morphometry-phenotype associations, and evaluated the results.
High-level schematic flow of the 3D image processing protocol: (A) 3D binary mask data; (B) mathematical representation and modeling of shape and size; (C) calculation of derived intrinsic and extrinsic geometric measures; and (D) machine learning based classification, feature ranking, and analysis.
Dataset description
In this study we used the 3D Cell Nuclear Morphology Microscopy Imaging Dataset, one of the largest publicly available 3D cell imaging datasets to date18. This dataset consists of two collections of 3D volumetric microscopic cell images with corresponding nuclear and nucleolar binary masks. Each collection includes images of cells in two phenotypic states, and thus poses a binary classification problem with image-level labels that can be used for the assessment of cell nuclear and nucleolar morphometric analysis. Binary masks in each collection are obtained by segmentation of the original data. Nuclear masks are extracted from a DAPI (4â²,6-diamidino-2-phenylindole) channel, while fibrillarin antibody-stained (anti-fibrillarin) and ethidium bromide-stained (EtBr) channels are both used for nucleolar binary mask extraction (see18 for details). Segmented binary masks are represented by 1024âÃâ1024âÃâZ 3D TIFF sub-volumes. For every mask sub-volume, accompanying vendor metadata extracted from the original data are available for analysis, as well.
Robust smooth surface reconstruction
To model the 3D shape of cell nuclei and nucleoli, boundaries of their 3D masks extracted from the microscopy data are modeled as genus zero two-dimensional manifolds (homeomorphic to a 2-sphere S2)55 that are embedded as triangulated surfaces in \({{\mathbb{R}}}^{3}\), Fig. 1B. Our approach uses an iterative Laplace-Beltrami eigen-projection and a topology-preserving boundary deformation algorithm26. This algorithm performs robust reconstruction of the objectsâ surfaces from their segmented masks using an iterative mask filtering process. First, a mesh representation is constructed from the boundary of an objectâs binary mask. Then, the boundary is projected onto the subspace of its LaplaceâBeltrami eigen-functions27, which allows the algorithm to automatically locate the position of spurious features by computing the metric distortion in eigen-projection. LB eigen-functions are intrinsically defined and can be easily computed from the boundary surface with no need of any parameterizations. They are also isometry invariant, and thus robust to the jagged nature of the boundary surface, which is desirable for biomedical shape analysis56. In our prior experience26, the discretized LB spectrum captures intrinsic shape characteristics (e.g., global shape transformations will preserve the spectral signature). The magnitude of the eigenvalues of the LB operator intuitively corresponds to the frequency in Fourier analysis, thus it provides a convenient mechanism to control the smoothness of the reconstructed surface. Using this information, the second step is a mask deformation process that only removes the spurious features while keeping the rest of the mask intact, thus preventing unintended volume shrinkage. This deformation is topology-preserving and well-composed such that the boundary surface of the mask is a manifold. The last two steps iterate until convergence and the method generates the final surface as the eigen-projection of the mask boundary, which is a smooth surface with genus zero topology26. These properties allow application of this algorithm to any shape, including, for example, crescent-shaped, multi-lobed, and folded, as long as shape topology is homeomorphic to a sphere. The exemplar results of this step performed on nuclear and nucleolar masks are shown in Fig. 2.
Robust smooth surface reconstruction. 3D visualization of: (A) a binary mask representation of a nucleus segmented from a Fibroblast cell image; (B) a mesh representation of a reconstructed smooth surface of a nucleus; (C) three binary masks for nucleoli segmented within this nucleus; and (D) three mesh representations of nucleolar surfaces, color-coded along the Z axis. Visualizations are produced with the SOCR Dynamic Visualization Toolkit web application65.
Morphometric measures
In this study, we used six shape measures as features quantifying geometric characteristics of the 3D surfaces, Fig. 1C. To calculate these measures, first the principal (min and max) curvatures \(({\kappa }_{1}\le {\kappa }_{2})\) were computed using triangulated surface models representing the boundaries of genus zero solids57. Then, shape morphometry measures can be expressed in terms of principal curvatures: mean curvature as \(MC=\frac{{\kappa }_{1}+{\kappa }_{2}}{2}\), shape index as \(SI=\frac{2}{\pi }\arctan (\frac{{\kappa }_{1}+{\kappa }_{2}}{{\kappa }_{2}-{\kappa }_{1}})\), and curvedness as \(CV=\sqrt{\frac{{\kappa }_{1}^{2}+{\kappa }_{2}^{2}}{2}}\). The principal curvatures of a surface are the eigenvalues of the Hessian matrix (second fundamental form), which solve for \(k\,,|H-kI|=0\), where I is the identity matrix. If S is a surface with second fundamental form H(X, Y), \(p\in M\) is a fixed point, and we denote an orthonormal basis u, v of tangent vectors at p, then the principal curvatures are the eigenvalues of the symmetric Hessian matrix, \(H=\,[\begin{array}{cc}{H}_{u,u} & {H}_{u,v}\\ {H}_{v,u} & {H}_{v,v}\end{array}]={H}_{uu}\partial {u}^{2}+2{H}_{u,v}\partial u\partial v+{H}_{v,v}\partial {v}^{2}\), a.k.a. shape tensor. Let râ=âr(u, v) be a parameterization of the surface \(S\subseteq {R}^{3}\), representing a smooth vector valued function of two variables with partial derivatives with respect to u and v denoted by ru and rv, Fig. 3. Then, the Hessian coefficients Hi,j at a given point (p) in the parametric u, v-plane are given by the projections of the second partial derivatives of r at that point onto the normal to \(S,\,n=\frac{{r}_{u}\times {r}_{v}}{|{r}_{u}\times {r}_{v}|}\), and can be computed using the dot product operator: \({H}_{u,u}={r}_{u,u}\cdot n,\) \({H}_{u,v}={H}_{v,u}={r}_{u,v}\cdot n,\,{H}_{v,v}={r}_{v,v}\cdot n\), Fig. 3.
Volume is the amount of 3D space enclosed by a closed boundary surface and can be expressed as \(V={\iiint }_{{{\mathbb{R}}}^{3}}{I}_{D}(x,y,z)dxdydz\), where \({I}_{D}(x,y,z)\) represents the indicator function of the region of interest (D)58. If r(u, v) is a continuously differentiable function and the normal vector to the surface over the appropriate region D in the parametric u, v plane is denoted by \({\overrightarrow{r}}_{u}\times {\overrightarrow{r}}_{v}\), then \({S}_{\Omega }:r\,=\,r(u,v),\,(u,v)\in \Omega \), is the parametric surface representation of the region boundary59. Then surface area can be expressed as \(SA={\iint }_{{\rm{\Omega }}}|{\mathop{r}\limits^{\rightharpoonup }}_{u}\times {\mathop{r}\limits^{\rightharpoonup }}_{v}|dudv\). The fractal dimension calculations are based on the fractal scaling down ratio, Ï, and the number of replacement parts, N60. Accurate, discrete approximations of these metrics are used to compute them on mesh-represented surfaces53,61. These discrete metrics were first introduced as a part of the shape analysis protocol31 and were further applied in neuroimaging studies32,33.
The extracted 3D morphometric measures serve as features for training a number of machine learning algorithms in order to assess classification performance, Fig. 1D. The number of detected nucleoli per nucleus is included as an individual feature. We merged nucleoli-level features within each nucleus by computing sample statistics (e.g., average, minimum, maximum, and higher moments) for each morphometry measure18. These statistics are used to augment the signature feature vectors of the corresponding parent nuclei such that all feature vectors are of the same length. Correspondingly, nuclei that do not have any automatically detected internally positioned nucleoli were excluded from further analysis, such that for each nucleus there was at least one nucleolus.
Visual analytics and machine learning for morphometric analysis
We performed exploratory visual analysis of extracted features using SOCRAT49, a web platform for interactive visual analytics. The goal of visual analytics is to support analytical reasoning and decision making with a combination of highly interactive visualizations and data analysis techniques49,62. SOCRAT implements a visual analytics workflow that encompasses an iterative process, in which data analysts can interactively interrogate extracted morphometric measures in the form of interactive dialogue supported by visualizations and data analysis components. In order to assess the variability of extracted morphometry data, we include t-Distributed Stochastic Neighbor Embedding (t-SNE)63 visualizations of the feature space generated by SOCRAT49. We also used SOCRAT to demonstrate interactions between the top-3 important features according to the best-performing classification algorithm. All derived morphometric datasets are made available within SOCRAT Web Demo application64. Finally, we demonstrated the ability to visualize volumetric images and extracted meshes online via SOCR Dynamic Visualization Toolkit web application65.
In general, correct classification of every single cell (type, stage, treatment, etc.) is a challenging task due to significant population heterogeneity of the observed cell phenotypes. For example, the same sample may contain a close mixture of intertwined âcancerousâ and ânon-cancerousâ cell phenotypes; or, both classes may include apoptotic cells exhibiting similar shapes or sizes. Given the nature of cell samples, culturing, preparation, and collection, we have considered classification of cell sets rather than single cells. The idea of classifying sets of cells, rather than individual samples, is not new and has been used in recent biomedical image classification studies12,66. The rationale behind this is based upon the observation that even if an algorithm misclassifies a few cells in a sample, the final (cell set) label will still be assigned correctly, as long as majority of cells are classified correctly. Using this strategy, we performed classification on small groups of cells, ranging from 3 to 31 cells per set. During each fold of the internal cross-validation, these small cell sets are randomized by bootstrapping procedure with 1,000 repetitions. Random uniform sub-sampling is used to resolve the sample-size imbalance between the classes. Due to the possible presence of batch effects in data, we employed the Leave-2-Opposite-Groups-Out (L2OGO) cross-validation scheme18. L2OGO ensures that: (1) all masks derived from one image fall either in the training or testing set, and (2) the testing set always contains masks from 2 images of different classes. We used scikit-learn, a popular Python machine learning toolkit67, to evaluate a number of supervised classification algorithms.
High-throughput workflow protocol
While the LONI Pipeline is a popular tool in neuroimaging and bioinformatics, it has been overlooked by the bioimage analysis community. In this work, we utilized the LONI Pipeline for the implementation of a streamlined multi-step protocol that relies on a diverse set of tools and solutions seamlessly connected in the LONI Pipeline workflow, Fig. 4. From a high-level perspective, every step of data processing and analysis protocol is wrapped as an individual module in the workflow that provides input and output specifications that allow the Pipeline to automatically connect and manage atomic modules. The modular structure of our implementation makes it highly flexible and not limited to specific tools included in the workflow. It can be repurposed for a wide range of different experiments by adjusting parameters, adding, removing, or replacing individual modules, while preserving high-throughput capabilities, as presented in the Discussion section. Every module represents an independent component that can be used in a stand-alone fashion. As a result, a distributed, massively parallel implementation of our protocol makes it possible to easily process thousands of nuclei and nucleoli simultaneously. The workflow does not depend on the total number of 3D objects, biological conditions, or a number of running instances since its execution is completely automated once the workflow configuration is fixed, including job scheduling and resource allocation. During the execution, our workflow provides a researcher with real-time information about progress and allows the viewing of intermediate results at every individual step. In addition, failed modules may easily be restarted.
Screenshots of the exemplar graphical workflow in the LONI Pipeline client interface that include: (left) overview of the validated workflow protocol showing nested groups of modules; (A) expanded Volume to Shape group that includes modules that perform 3D shape modeling refinement; and (B) expanded Morphometry group that includes a module that performs morphological measure extraction.
The workflow is configured in such a way that it can consume data in the specific format we used, i.e. 1024âÃâ1024âÃâZ 3D volumes in different channels as 16-bit 3D TIFF files. Each volume is processed independently, in parallel fashion, such that workflow automatically defines how many processes are needed to analyze all of the input data. 3D shape modeling and morphometric feature extraction are performed on individual masks independently, which allows us to simultaneously run up to 1,200 jobs on the cluster during our experiments, effectively reducing the computing time. Finally, the workflow collects morphometry information from each individual mask and combines them in the results table that is further used as an input to classification algorithm. These capabilities allow the user to take advantage of modern computational resources, lift the burden of low-level configuration from researchers, make it easier to control the execution process, and improve reproducibility of the whole process.
Results
Validation on synthetic data
To validate the shape morphometry metrics, we first applied them to synthetically generated 3D masks. We used the scikit-image Python library68 to create 3D solids representing cubes, octahedra, spheres, ellipsoids, and 3 overlapping spheres with linearly aligned centers (Supplementary Fig. S1). We processed these objects and compare the resulting shape morphometry measures. Specifically, we aimed to confirm the expected close relation between the analytically derived measures of volume and surface area computed using the corresponding shape parameters (e.g., radius, size), and their computationally derived counterparts reported by the processing pipeline workflow. Our results illustrate that for nucleus-like shapes, e.g., spheres and ellipsoids, the computational error is within 2%. For faceted objects, e.g., cubes and octahedrons, the calculation error is within 6%. The increased error in the latter case can be explained by the mesh smoothing the surface reconstruction algorithm applies at the shape vertices to resolve points of singularity (e.g., smooth but non-differentiable surface boundaries).
To demonstrate the detection of shape differences between different types of 3D objects, we also compared overlapping spheres against circumscribed ellipsoids. As expected, the average mean curvature and curvedness measures are lower and shape index values are higher for spheres compared to ellipsoids. We observed a similar trend when comparing changes in these shape morphometry measures for spheres, ellipsoids, and overlapping spheres. For example, average mean curvature and curvedness were highest for overlapping spheres and lowest for spheres, which is expected based on definitions of these measures (see Supplementary Table S2). This simulation confirms our ability to accurately measure size and shape characteristics of 3D objects, which forms the basis for machine-learning based object classification based on boundary shapes. Exemplar results of synthetic data morphometry are available in Supplementary Table S1.
Comparison with SPHARM for fibroblast nuclei shape classification
Classification of single cell nuclei from the fibroblast collection of the 3D Cell Nuclear Morphology Microscopy Imaging Dataset18 may be assessed using shape morphometry metrics as salient discriminatory features, which we did and compared against their corresponding SPHARM coefficients16,25. We used images of primary human fibroblast cells that were subjected to a G0/G1 Serum Starvation Protocol for cell cycle synchronization69. This protocol has previously been shown to alter nuclear organization, which was reflected in the observed morphology changes such as nuclear size and shape70. As a result, this collection contains 962 3D nuclear binary masks in the following phenotypic classes: (1) proliferating fibroblasts (PROLIF), and (2) cell cycle synchronized by the serum-starvation protocol cells (SS). We used these binary nuclear masks to calculate both SPHARM and morphometric features.
We obtained the SPHARM coefficients by using the popular SPHARM-MAT toolbox71 that implements surface reconstruction and spherical parametrization using the CALD algorithm25. Then, we followed by the expansion of the object surface into a complete set of spherical harmonic basis functions of degree lâ=â13 (default setting). Finally, SHREC method72 was used to minimize the mean square distance between corresponding surface parts. SPHARM shape descriptors were computed as described by Ducroz et al.25 and used as feature vectors for classification.
Throughout, we employed the open-source Python package scikit-learn 0.17.067 to test a number of commonly used machine learning classification methods on derived feature vectors with default parameters for each method. Performance was compared using the L2OGO cross-validation scheme and the area under the receiver operating characteristic curve (AUC) as a performance metric. As shown in Table 1, 3D shape morphometric measures not only demonstrate comparable discriminative performance to SPHARM coefficients, but outperform them using all tested algorithms.
Fibroblast cell classification
The full collection of fibroblast masks for binary classification consists of total 965 nuclei (498 SS and 470 PROLIF) and 2,181 nucleoli (1,151 SS and 1,030 PROLIF). Figure 5A demonstrates the variability of the extracted morphometry measures in a t-SNE projection visualized in SOCRAT. Although there was a small degree of grouping, there was no clear separation between classes.
Fibroblast morphometric analysis: (A) SOCRAT visualization of t-SNE projection of morphometric feature space; (B) mean AUC for various cell set sizes; (C) top-10 features for classification by importance score (right, nucleolar feature names start with Avg, Min, Max or Var, feature names that were also reported in top-10 for PC3 cells are shown in blue font); and (D): SOCRAT visualization of interactions between top-3 features.
The best result by a single classifier was achieved using a stochastic gradient boosting classifier with 1,500 base learners, maximum tree depth 8, subsampling rate 0.5. Hyper-parameters were fine-tuned using a cross-validated grid search. To evaluate these classification results, we measured accuracy, precision, sensitivity, and AUC over L2OGO cross-validation, which are presented in Table 2 for single cell and 19-cell-set classifications. Figure 5B shows mean AUC values for set sizes from 3 to 19 cells. A 90% mean AUC was reached when classifying sets with 19 cells and 92.5% for sets with 31 or more cells.
The gradient boosting classifier also computes and reports cross-validated feature importance (Fig. 5C). These allow us to evaluate which measures differ between two cell conditions, and potentially propose novel research hypotheses that can be tested using prospective data. Previous analysis has reported quantifiable changes in both nuclear size and shape under serum-starvation70. In our results, both nuclear (top-6, out of top-10) and nucleolar (4 of top-10) morphometric size and shape features are reported to be of high importance for distinguishing SS fibroblasts from PROLIF (Fig. 5C). We also visualized the relationship between top-3 features using SOCRAT, see Fig. 5D. Visualizations suggest the smaller variation of morphometric measures in SS fibroblast nuclei compared to their PROLIF counterparts. This result may provide insight in further downstream analysis of potential underlying mechanisms that lead to these morphometric changes. We made the fibroblast morphometry data publicly available within SOCRAT for further analysis and validation64.
PC3 EPI/EMT cell classification
The second collection contains images of human prostate cancer cells (PC3). Through the course of progression to metastasis, malignant cancer cells undergo a series of reversible transitions between intermediate phenotypic states bounded by pure epithelium and pure mesenchyme9. These transitions in prostate cancer are associated with quantifiable changes in both nuclear and nucleolar structure10,73. PC3 cells were cultured in: (1) epithelial (EPI), and (2) mesenchymal transition (EMT) phenotypic states. The collection includes 458 nuclear (310 EPI and 148 EMT) and 1,101 nucleolar (649 EPI and 452 EMT) 3D binary masks. Figure 6A demonstrates the variability of the extracted morphometry measures in a t-SNE projection visualized in SOCRAT. Similar to fibroblasts, the projection of the PC3 morphometric feature space does not demonstrate clear separation between classes.
PC3 morphometric analysis: (A) SOCRAT visualization of t-SNE projection of morphometric feature space; (B) mean AUC for various cell set sizes; (C) top-10 features for classification by importance score (right, nucleolar feature names start with Avg, Min, Max or Var, feature names that were also reported in top-10 for Fibroblast cells are shown in blue font); and (D): SOCRAT visualization of interactions between top-3 features.
In this case, the best classification performance by single classifier is the result of applying a random forest model (1,000 trees, maximum tree depth 12, maximum number of features for the best split 40%). Hyper-parameters fine-tuning, accuracy metrics, and cross-validation procedures are identical to the ones reported in the previous fibroblast experiment. Classification of sets of 19 cells achieves a mean AUC of 76.2%, Table 3. Figure 6B reports the AUC for different group sizes to show how the classification performance increases with the cell-set size and reaches 80% for sets of 27 or more cells. In this experiment, we also examined the classifier-reported feature importance, Fig. 6C. The top-10 important features in this classification included nuclear (4 of top-10, which were also in Fibroblast top-6) and nucleolar (top-3, 6 out of top-10) shape morphometry features. Top feature interactions visualized using SOCRAT demonstrate the important changes in distributions of nucleolar morphometric measures, Fig. 6D. For example, it seems that the EPI nucleoli tend to have more variability in minimal curvedness and average fractal dimension, compared to EMT nucleoli. Previously reported PC3 morphological analyses73 only used simple 2D nuclear form measures, such as diameter and the size of the bounding box. While we confirmed the importance of nuclear form in our results and suggested the need for further investigation of other highly ranked features, such as nucleolar curvedness, shape index, and fractal dimension, which may provide additional mechanistic insights. PC3 morphometry data are made publicly available within SOCRAT for further analysis and validation49.
Discussion
In this study, we proposed, implemented, and validated a solution for 3D modeling, morphological feature extraction, analysis, and classification of cells by treatment conditions. Compared to other studies using 2D projections, this approach operates natively in 3D space and takes advantage of extrinsic and intrinsic morphometric measures that are more representative of the real, underlying nuclear and nucleolar geometry and allow easy human interpretation. Given the limitations of using 3D voxels for accurate shape representation, we employed 3D surface models to extract more informative size and shape measures to improve the morphology classification performance. Robust surface reconstruction allows accurate approximation of 3D object boundaries that was validated on synthetic data. Suggested shape morphometric measures outperform another popular approach and demonstrated their universality across different cell types, conditions, and even domains31,32,33.
Our computational protocol implementation is highly parallel with throughput, limited only by the number of available computing nodes, and it can process thousands of objects simultaneously with minimal human intervention. This pipeline workflow integrates a number of open-source tools for different steps of data processing and analytics. Every module in our workflow represents an individual component that can be easily modified, removed, or replaced by an alternative. Such modular software platform architectures have been shown to enable high reusability and ease of modification49. This allows the user to use the same workflow or customize and expand it (e.g., specification of new datasets, swapping of specific atomic modules) for other purposes that require the analysis of a diverse array of cellular, nuclear, or other studies. The live demo available via the LONI Pipeline demonstrates the simplicity of use and high efficiency of parallel data processing. LONI also provides guest access (see Supplementary Information) and an opportunity to utilize a 4,500-core LONI cluster after applying for a collaboration account.
We tested our approach on the 3D Cell Nuclear Morphology Microscopy Imaging Dataset18, which includes a total of ~1,500 nuclear and ~2700 nucleolar masks. The classification results on these data comparing epithelial vs. mesenchymal human prostate cancer cell lines, and serum-starved vs. proliferating fibroblast cell lines, demonstrate the high accuracy of cell type prediction using 3D morphometry, especially when applied to sets of cells. Although different classification algorithms appear to be optimal for different experiments, we observed that both nuclear and nucleolar morphometric measures are important features for discriminating between treatment conditions or cell phenotypes. In the case of fibroblast classification, the results show the importance of nuclear morphometry, the number of nucleoli per nucleus, and various internal nucleolar morphometric measures. These observations confirm and extend previously reported results. For PC3 cells, the most important classification features are the moments of the distributions of various nucleolar morphometric measures, along with nuclear size and shape. Interestingly, there were 3 common morphometric features among the top-10 most important ones for both cell lines. This confirms previously reported observations73, suggests new important morphological characteristics, and demonstrates that our method extracts relevant information from cell forms to successfully classify cells using a combination of criteria. In addition, this also demonstrates the importance of sophisticated shape metrics, compared to volume and surface area, that alone, were not the most informative features for the classification results. The use of SOCRAT enables interactive interrogation of morphometric data in a visual manner, supported by analytical tools. This method of interactive âvisual analyticsâ provides insight into feature dependencies and interactions, and can be used for result interpretation. We also demonstrated the visualization of 3D volumetric images and derived meshed surface representations using the SOCR Dynamic Visualization Toolkit web application65.
Our computational approach is scalable and capable of processing complex big 3D imaging data, and is not limited to nuclear and nucleolar shapes. With some changes, it can be applied to other cellular and nuclear compartments of interest. More specifically, the robust smooth surface reconstruction algorithm can be directly applied to any 3D shapes, as long as their topology is sphere-like. Together with molecular level techniques, such as Hi-C, our 3D shape morphometry workflow can form a powerful combination for the investigation of DNA architecture in the spatial and temporal framework of the 4D nucleome6,74. One example of the many possible future applications of this workflow is to study asymmetric cell division. Stem and progenitor cells are characterized by their ability to self-renew and produce differentiated progeny. A balance between these processes is achieved through controlled asymmetric divisions and is necessary to generate cellular diversity during development and to maintain adult tissue homeostasis. Disruption of this balance may result in premature depletion of the stem/progenitor cell pool, or abnormal growth75,76. In many tissues, dysregulated asymmetric divisions are associated with cancer. Whether there is a causal relationship between asymmetric cell division defects and cancer initiation is unknown. We proposed that our shape analysis pipeline will be useful in studying the 4D nucleome topology of morphogenesis and cancer initiation.
As one of the approach limitations, we pointed out that other geometric measures can be used to characterize shapes of interest, such as intrinsic shape context, compactness, symmetry, smoothness, convexity, etc. In the current representation, analyzable shapes are limited to genus zero surfaces, which is a fair assumption when modeling objects like nuclei or nucleoli. However, it might be not trivial when considering other nuclear structures, for example, chromosome territories or interchromosomal loops, since their topologies may not be homeomorphic to a sphere, or may not appear to be genus zero under some imaging conditions and modalities. It is also conceivable, yet not very likely for the discretized LB, that 2 different shapes may have the same spectra. In this case, the algorithm may fail to detect the intrinsic differences between them due to false-negative error. Even though our workflow only requires little intervention (classifier selection and tuning), further improvements would involve adaptive implementations with even less manual intervention, as well as extraction of additional features. Another option is to use deep learning-based methods that alleviate the need to define features and allow to learn relevant patterns directly from data77,78. Recent applications have demonstrated the ability of deep neural networks to successfully perform classification, segmentation, and detection on limited amounts of biomedical imaging data79,80,81. For example, textural features could possibly increase discriminatory power of the method and provide more information on chromatin reorganization82. Since nuclear deformation serves as a proxy to underlying processes, the importance of particular features and the methodâs ability to classify nuclei does not provide direct insight into the fundamental biological mechanism driving the observed morphometric differences between cell phenotypes or environmental conditions. The computational results should be further tested and externally validated using other experimental conditions and prospective data.
Conclusions
Quantification of cell nuclear morphology enables more subtle characterization of cellular phenotypic traits, which can be associated with functional changes coupled to underlying biological processes. Using the new methodology described in this paper, we compared the morphology of serum-starved vs. proliferating fibroblast cells as a control, followed by a comparison of epithelial with mesenchymal human prostate cancer cell lines. In the case of fibroblast classification, our results show the importance of nuclear morphometric change, along with the number of detected nucleoli per nucleus, and various internal nucleolar morphometric measures. Results for PC3 cells demonstrate that the changes in nucleolar morphology are the most informative. However, in both cell lines, both nuclear and nucleolar morphometric measures contribute to the discriminative power of the classification algorithms. To the best of our knowledge, this study is the first where a 3D morphometric assay could easily distinguish between the epithelial and mesenchymal cell nuclei. These findings suggest that further investigation of highly ranked features that were not previously reported, such as nucleolar curvedness, shape index, and fractal dimension, may provide interesting mechanistic insights.
The ability to automate the processes of specimen collection, image acquisition, data pre-processing, computation of derived biomarkers, modeling, classification, and analysis can significantly impact clinical decision-making and fundamental investigation of cell deformation. To our knowledge, this is the first attempt to combine 3D cell nuclear shape modeling by robust smooth surface reconstruction and extraction of shape morphometry measures into a highly parallel pipeline workflow protocol for morphological analysis of thousands of nuclei and nucleoli in 3D. This approach allows efficient and informative evaluation of cell shapes in the imaging data and represents a reproducible technique that can be validated, modified, and repurposed by the biomedical community. This facilitates result reproducibility, collaborative method validation, and broad knowledge dissemination.
Availability of Materials and Data
The documentation supporting the conclusions of this article together with the pipeline workflows and underlying source code are made available online on the project webpage SOCR 3D Cell Morphometry Project, http://socr.umich.edu/projects/3d-cell-morphometry54.
Change history
26 October 2018
A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.
References
Jevtic, P., Edens, L. J., Vukovic, L. D. & Levy, D. L. Sizing and shaping the nucleus: mechanisms and significance. Curr Opin Cell Biol 28, 16â27, https://doi.org/10.1016/j.ceb.2014.01.003 (2014).
Bustin, M. & Misteli, T. Nongenetic functions of the genome. Science 352, aad6933, https://doi.org/10.1126/science.aad6933 (2016).
Dundr, M. & Misteli, T. Biogenesis of nuclear bodies. Cold Spring Harb Perspect Biol 2, a000711, https://doi.org/10.1101/cshperspect.a000711 (2010).
Fedorchak, G. R., Kaminski, A. & Lammerding, J. Cellular mechanosensing: getting to the nucleus of it all. Prog Biophys Mol Biol 115, 76â92, https://doi.org/10.1016/j.pbiomolbio.2014.06.009 (2014).
Uhler, C. & Shivashankar, G. V. Regulation of genome organization and gene expression by nuclear mechanotransduction. Nat Rev Mol Cell Biol 18, 717â727, https://doi.org/10.1038/nrm.2017.101 (2017).
Higgins, G. A., Allyn-Feuer, A., Handelman, S., Sadee, W. & Athey, B. D. The epigenome, 4D nucleome and next-generation neuropsychiatric pharmacogenomics. Pharmacogenomics 16, 1649â1669, https://doi.org/10.2217/pgs.15.111 (2015).
Meaburn, K. J. Spatial Genome Organization and Its Emerging Role as a Potential Diagnosis Tool. Front Genet 7, 134, https://doi.org/10.3389/fgene.2016.00134 (2016).
Zink, D., Fischer, A. H. & Nickerson, J. A. Nuclear structure in cancer cells. Nature reviews. Cancer 4, 677â687, https://doi.org/10.1038/nrc1430 (2004).
Veltri, R. W. & Christudass, C. S. Nuclear morphometry, epigenetic changes, and clinical relevance in prostate cancer. Advances in experimental medicine and biology 773, 77â99, https://doi.org/10.1007/978-1-4899-8032-8_4 (2014).
Montanaro, L., Trere, D. & Derenzini, M. Nucleolus, ribosomes, and cancer. The American journal of pathology 173, 301â310, https://doi.org/10.2353/ajpath.2008.070752 (2008).
Pincus, Z. & Theriot, J. A. Comparison of quantitative methods for cell-shape analysis. Journal of microscopy 227, 140â156, https://doi.org/10.1111/j.1365-2818.2007.01799.x (2007).
Huang, H. et al. Cancer diagnosis by nuclear morphometry using spatial information. Pattern Recognit Lett 42, 115â121, https://doi.org/10.1016/j.patrec.2014.02.008 (2014).
Caserta, F. et al. Determination of fractal dimension of physiologically characterized neurons in two and three dimensions. J Neurosci Methods 56, 133â144 (1995).
Choi, H. J. & Choi, H. K. Grading of renal cell carcinoma by 3D morphological analysis of cell nuclei. Comput Biol Med 37, 1334â1341, https://doi.org/10.1016/j.compbiomed.2006.12.008 (2007).
Meyer, M. G. et al. Automated cell analysis in 2D and 3D: A comparative study. Pattern Recogn 42, 141â146 (2009).
Dufour, A. C. et al. Signal Processing Challenges in Quantitative 3-D Cell Morphology: More than meets the eye. IEEE Signal Processing Magazine 32, 30â40, https://doi.org/10.1109/MSP.2014.2359131 (2015).
Pegoraro, G. & Misteli, T. High-throughput Imaging as a versatile and unbiased discovery tool. Methods (San Diego, Calif.) 96, 1â2, https://doi.org/10.1016/j.ymeth.2016.01.003 (2016).
Kalinin, A. A. et al. 3D Cell Nuclear Morphology: Microscopy Imaging Dataset and Voxel-Based Morphometry Classification Results In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition Workshops, 2272â2280 (2018).
Brechbuhler, C., Gerig, G. & Kubler, O. Parametrization of Closed Surfaces for 3-D Shape-Description. Comput Vis Image Und 61, 154â170 (1995).
Antoine, J. P. & Vandergheynst, P. Wavelets on the 2-sphere: A group-theoretical approach. Appl Comput Harmon A 7, 262â291 (1999).
Khairy, K., Foo, J. & Howard, J. Shapes of Red Blood Cells: Comparison of 3D Confocal Images with the Bilayer-Couple Model. Cell Mol Bioeng 1, 173â181, https://doi.org/10.1007/s12195-008-0019-5 (2010).
Singh, S. et al. Non-parametric population analysis of cellular phenotypes. Med Image Comput Comput Assist Interv 14, 343â351 (2011).
Ducroz, C., Olivo-Marin, J. C. & Dufour, A. Characterization of cell shape and deformation in 3D using Spherical Harmonics In 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI), 848â851 (IEEE, 2012).
Du, C. J., Hawkins, P. T., Stephens, L. R. & Bretschneider, T. 3D time series analysis of cell shape using Laplacian approaches. BMC bioinformatics 14 (2013).
Shen, L. & Makedon, F. Spherical mapping for processing of 3D closed surfaces. Image Vision Comput 24, 743â761 (2006).
Shi, Y. G. et al. Robust Surface Reconstruction via Laplace-Beltrami Eigen-Projection and Boundary Deformation. IEEE transactions on medical imaging 29, 2009â2022 (2010).
Levy, B. Laplace-Beltrami Eigenfunctions Towards an Algorithm That âUnderstandsâ Geometry In IEEE International Conference on Shape Modeling and Applications 2006 (SMIâ06), 13-13 (2006).
Han, X., Xu, C. Y. & Prince, J. L. A topology preserving level set method for geometric deformable models. IEEE Transactions on Pattern Analysis and Machine Intelligence 25, 755â768 (2003).
Batchelor, P. G. et al. Measures of folding applied to the development of the human fetal brain. IEEE transactions on medical imaging 21, 953â965, https://doi.org/10.1109/TMI.2002.803108 (2002).
Koenderink, J. J. & Vandoorn, A. J. Surface Shape and Curvature Scales. Image Vision Comput 10, 565â565 (1992).
Dinov, I. D. et al. Efficient, Distributed and Interactive Neuroimaging Data Analysis Using the LONI Pipeline. Front Neuroinform 3, 22, https://doi.org/10.3389/neuro.11.022.2009 (2009).
Fani, N. et al. FKBP5 and attention bias for threat: associations with hippocampal function and shape. JAMA Psychiatry 70, 392â400, https://doi.org/10.1001/2013.jamapsychiatry.210 (2013).
Moon, S. W. et al. Structural Brain Changes in Early-Onset Alzheimerâs Disease Subjects Using the LONI Pipeline Environment. Journal of neuroimaging: official journal of the American Society of Neuroimaging 25, 728â737, https://doi.org/10.1111/jon.12252 (2015).
Ramo, P., Sacher, R., Snijder, B., Begemann, B. & Pelkmans, L. CellClassifier: supervised learning of cellular phenotypes. Bioinformatics (Oxford, England) 25, 3028â3030, https://doi.org/10.1093/bioinformatics/btp524 (2009).
Held, M. et al. CellCognition: time-resolved phenotype annotation in high-throughput live cell imaging. Nature methods 7, 747â754, https://doi.org/10.1038/nmeth.1486 (2010).
Pau, G., Fuchs, F., Sklyar, O., Boutros, M. & Huber, W. EBImageâan R package for image processing with applications to cellular phenotypes. Bioinformatics (Oxford, England) 26, 979â981, https://doi.org/10.1093/bioinformatics/btq046 (2010).
Kamentsky, L. et al. Improved structure, function and compatibility for CellProfiler: modular high-throughput image analysis software. Bioinformatics (Oxford, England) 27, 1179â1180, https://doi.org/10.1093/bioinformatics/btr095 (2011).
Chiang, M. et al. Analysis of in vivo single cell behavior by high throughput, human-in-the-loop segmentation of three-dimensional images. BMC bioinformatics 16, 397, https://doi.org/10.1186/s12859-015-0814-7 (2015).
Peng, T. & Murphy, R. F. Image-derived, three-dimensional generative models of cellular organization. Cytometry. Part A: the journal of the International Society for Analytical Cytology 79, 383â391, https://doi.org/10.1002/cyto.a.21066 (2011).
Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nature methods 9, 676â682, https://doi.org/10.1038/nmeth.2019 (2012).
de Chaumont, F. et al. Icy: an open bioimage informatics platform for extended reproducible research. Nature methods 9, 690â696, https://doi.org/10.1038/nmeth.2075 (2012).
Kankaanpaa, P. et al. BioImageXD: an open, general-purpose and high-throughput image-processing platform. Nature methods 9, 683â689, https://doi.org/10.1038/nmeth.2047 (2012).
Ollion, J., Cochennec, J., Loll, F., Escude, C. & Boudier, T. TANGO: a generic tool for high-throughput 3D image analysis for studying nuclear organization. Bioinformatics (Oxford, England) 29, 1840â1841, https://doi.org/10.1093/bioinformatics/btt276 (2013).
Eliceiri, K. W. et al. Biological imaging software tools. Nature methods 9, 697â710, https://doi.org/10.1038/nmeth.2084 (2012).
Peng, H., Bria, A., Zhou, Z., Iannello, G. & Long, F. Extensible visualization and analysis for multidimensional images using Vaa3D. Nature protocols 9, 193â208, https://doi.org/10.1038/nprot.2014.011 (2014).
Li, L., Zhou, Q., Voss, T. C., Quick, K. L. & LaBarbera, D. V. High-throughput imaging: Focusing in on drug discovery in 3D. Methods (San Diego, Calif.) 96, 97â102, https://doi.org/10.1016/j.ymeth.2015.11.013 (2016).
Ince, D. C., Hatton, L. & Graham-Cumming, J. The case for open computer programs. Nature 482, 485â488, https://doi.org/10.1038/nature10836 (2012).
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nature methods 9, 671â675 (2012).
Kalinin, A. A., Palanimalai, S. & Dinov, I. D. SOCRAT Platform Design: A Web Architecture for Interactive Visual Analytics Applications In Proceedings of the 2nd Workshop on Human-In-the-Loop Data Analytics, HILDAâ17 1â6 (ACM, 2017).
Dinov, I. D. et al. Neuroimaging study designs, computational analyses and data provenance using the LONI pipeline. PloS one 5, https://doi.org/10.1371/journal.pone.0013070 (2010).
Dinov, I. D. et al. Applications of the pipeline environment for visual informatics and genomics computations. BMC bioinformatics 12, 304, https://doi.org/10.1186/1471-2105-12-304 (2011).
Thompson, P. M., Schwartz, C., Lin, R. T., Khan, A. A. & Toga, A. W. Three-dimensional statistical analysis of sulcal variability in the human brain. The Journal of neuroscience: the official journal of the Society for Neuroscience 16, 4261â4274 (1996).
Meyer, M., Desbrun, M., Schröder, P. & Barr, A. H. Discrete Differential-Geometry Operators for Triangulated 2-Manifolds In Visualization and Mathematics III, (eds Hans-Christian, H. & Konrad, P.) 35â57 (Springer Berlin Heidelberg, 2003).
SOCR 3D Cell Morphometry Project, http://socr.umich.edu/projects/3d-cell-morphometry (2018).
Ferri, M. & Gagliardi, C. The only genus zero n-manifold is S n. Proc. Amer. Math. Soc. 85, 638â642, https://doi.org/10.1090/S0002-9939-1982-0660620-5 (1982).
Niethammer, M. et al. Global medical shape analysis using the Laplace-Beltrami spectrum. Med Image Comput Comput Assist Interv 10, 850â857 (2007).
Terzopoulos, D. The computation of visible-surface representations. IEEE Transactions on Pattern Analysis and Machine Intelligence 10, 417â438, https://doi.org/10.1109/34.3908 (1988).
Larson, R. & Edwards, B. Calculus. 10th edn, (Cengage Learning, 2009).
Santaló, L. A. Integral geometry and geometric probability. 2nd edn, (Cambridge University Press, 2004).
Mandelbrot, B. B. The fractal geometry of nature. 1st edn, Vol. 173 (W. H. Freeman and Company, 1982).
Jagannathan, A. Segmentation and Recognition of 3D Point Clouds within Graph-theoretic and thermodynamic frameworks: a thesis, Northeastern University, 2005 (2005).
Husain, S. S., Kalinin, A., Truong, A. & Dinov, I. D. SOCR data dashboard: an integrated big data archive mashing medicare, labor, census and econometric information. J Big Data 2, https://doi.org/10.1186/s40537-015-0018-z (2015).
van der Maaten, L. & Hinton, G. Visualizing Data using t-SNE. J Mach Learn Res 9, 2579â2605 (2008).
SOCRAT Web Demo, http://socr.umich.edu/HTML5/SOCRAT/ (2018).
SOCR Dynamic Visualization Toolkit http://socr.umich.edu/HTML5/DViewer/ (2018).
Cheplygina, V., Tax, D. M. J. & Loog, M. On classification with bags, groups and sets. Pattern Recognition Letters 59, 11â17, https://doi.org/10.1016/j.patrec.2015.03.008 (2015).
Pedregosa, F. et al. Scikit-learn: Machine learning in Python. The Journal of Machine Learning Research 12, 2825â2830 (2011).
van der Walt, S. et al. scikit-image: image processing in Python. PeerJ 2, e453, https://doi.org/10.7717/peerj.453 (2014).
Langan, T. J. & Chou, R. C. Synchronization of mammalian cell cultures by serum deprivation. Methods Mol Biol 761, 75â83, https://doi.org/10.1007/978-1-61779-182-6_5 (2011).
Seaman, L., Meixner, W., Snyder, J. & Rajapakse, I. Periodicity of nuclear morphology in human fibroblasts. Nucleus 6, 408â416, https://doi.org/10.1080/19491034.2015.1095432 (2015).
SPHARM-MAT, http://www.iu.edu/~spharm/ (2018).
Shen, L., Farid, H. & McPeek, M. A. Modeling Three-Dimensional Morphological Structures Using Spherical Harmonics. Evolution 63, 1003â1016 (2009).
Verdone, J. E., Parsana, P., Veltri, R. W. & Pienta, K. J. Epithelial-mesenchymal transition in prostate cancer is associated with quantifiable changes in nuclear structure. The Prostate 75, 218â224, https://doi.org/10.1002/pros.22908 (2015).
Zheng, G. et al. Hypothesis: Caco-2 cell rotational 3D mechanogenomic Turing patterns has clinical implications to colon crypts. Journal of Cellular and Molecular Medicine, https://doi.org/10.1111/jcmm.13853 (2018).
Morrison, S. J. & Kimble, J. Asymmetric and symmetric stem-cell divisions in development and cancer. Nature 441, 1068â1074, https://doi.org/10.1038/nature04956 (2006).
Gomez-Lopez, S., Lerner, R. G. & Petritsch, C. Asymmetric cell division of stem and progenitor cells during homeostasis and cancer. Cell Mol Life Sci 71, 575â597, https://doi.org/10.1007/s00018-013-1386-1 (2014).
Ching, T. et al. Opportunities and obstacles for deep learning in biology and medicine. Journal of The Royal Society Interface 15, https://doi.org/10.1098/rsif.2017.0387 (2018).
Kalinin, A. A. et al. Deep learning in pharmacogenomics: from gene regulation to patient stratification. Pharmacogenomics 19, 629â650, https://doi.org/10.2217/pgs-2018-0008 (2018).
Rakhlin, A., Shvets, A., Iglovikov, V. & Kalinin, A. A. Deep Convolutional Neural Networks for Breast Cancer Histology Image Analysis In Image Analysis and Recognition, Lecture Notes in Computer Science. Chapter 83, 737â744 (2018).
Iglovikov, V. I., Rakhlin, A., Kalinin, A. A. & Shvets, A. A. Pediatric Bone Age Assessment Using Deep Convolutional Neural Networks. In 4th Deep Learning in Medical Image Analysis (DLMIA) Workshop, (Springer, 2018).
Shvets, A., Rakhlin, A., Kalinin, A. & Iglovikov, V. Automatic Instrument Segmentation in Robot-Assisted Surgery Using Deep Learning. bioRxiv, https://doi.org/10.1101/275867 (2018).
Kalinin, A. A., Athey, B. D. & Dinov, I. D. Evaluation of Methods for Cell Nuclear Structure Analysis from Microscopy Data In Supplementary Proceedings of the 7th International Conference on Analysis of Images, Social Networks and Texts (AIST-SUP 2018), (CEUR-WS, 2018).
Acknowledgements
Many colleagues at the Michigan Institute for Data Science (MIDAS), the Department of Computational Medicine and Bioinformatics, and the Statistics Online Computational Resource (SOCR) provided contributions including ideas, pilot testing, improvement suggestions and other assistance in the development and validation of these methods. We also would like to thank colleagues at the Laboratory of Neuro Imaging (LONI), at the Keck School of Medicine, University of Southern California, for providing technical support for the LONI Pipeline environment. This work was partially supported by the National Science Foundation grants 1734853, 1636840, 1023115, 1022560, 1022636, 0089377, 9652870, 0442992, 0442630, 0333672, 0716055, 1734853; the National Institutes of Health grants P20 NR015331, U54 EB020406, P50 NS091856, P30 DK089503, P30 AG053760, T32 GM070449, UL1TR002240; the Department of Computational Medicine and Bioinformatics Research Account, University of Michigan; the James Buchanan Brady Urological Institute, John Hopkins University, within the contract âPromoting Scientific Progress through Biomedical Research, Biomedical Informatics and the Development Of Ontological and Biomedical Informatics Tools that enable Collaborative Biomedical Researchâ (Sponsor ID 420900), and the Elsie Andresen Fiske Research Fund.
Author information
Authors and Affiliations
Contributions
A.A. and D.D. designed and built modules for data pre-processing and curation. A.A.K., A.A., D.D., A.A.F., G.V.F., W.M., J.R.d.W., G.A.H., G.Z., A.C., J.W.W. iteratively evaluated various image analysis approaches. A.A.K., A.A.F., S.S.H. and I.D.D. designed and built modules for 3D shape modeling, morphometric feature extraction, and visual analytics. A.A.K., A.A.F., A.A., and I.D.D. designed, built, and executed high-throughput Pipeline workflows. AAK implemented and performed statistical analysis and result interpretation. G.V.F., W.M., J.R.d.W., G.A.H., G.Z., A.C., J.W.W., J.E.V., R.W.V., K.J.P. and D.S.C. contributed to conceptualization of study. I.D.D. and B.D.A. conceived the study. A.A.K., A.A.F., A.A., I.D.D. and B.D.A. wrote the manuscript with input from other authors. All authors participated in numerous project discussions and decision making. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kalinin, A.A., Allyn-Feuer, A., Ade, A. et al. 3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification. Sci Rep 8, 13658 (2018). https://doi.org/10.1038/s41598-018-31924-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-31924-2