Segmentation and Modelling of the Nuclear Envelope of HeLa Cells Imaged with Serial Block Face Scanning Electron Microscopy

• Published in: Journal of imaging 2019-09, Vol.5 (9), p.75
• Main Authors: Karabağ, Cefa, Jones, Martin L., Peddie, Christopher J., Weston, Anne E., Collinson, Lucy M., Reyes-Aldasoro, Constantino Carlos
• Format: Article
• Language: English
• Subjects: active contours; Algorithms; Automation; Cervical cancer; Datasets; Hausdorff distance; HeLa cells; Image segmentation; Jaccard index; Medical research; Metric space; Morphology; Pixels; Scanning electron microscopy; segmentation; Similarity; Two dimensional models
• ISSN: 2313-433X; 2313-433X
• DOI: 10.3390/jimaging5090075

This paper describes an unsupervised algorithm, which segments the nuclear envelope of HeLa cells imaged by Serial Block Face Scanning Electron Microscopy. The algorithm exploits the variations of pixel intensity in different cellular regions by calculating edges, which are then used to generate superpixels. The superpixels are morphologically processed and those that correspond to the nuclear region are selected through the analysis of size, position, and correspondence with regions detected in neighbouring slices. The nuclear envelope is segmented from the nuclear region. The three-dimensional segmented nuclear envelope is then modelled against a spheroid to create a two-dimensional (2D) surface. The 2D surface summarises the complex 3D shape of the nuclear envelope and allows the extraction of metrics that may be relevant to characterise the nature of cells. The algorithm was developed and validated on a single cell and tested in six separate cells, each with 300 slices of 2000 × 2000 pixels. Ground truth was available for two of these cells, i.e., 600 hand-segmented slices. The accuracy of the algorithm was evaluated with two similarity metrics: Jaccard Similarity Index and Mean Hausdorff distance. Jaccard values of the first/second segmentation were 93%/90% for the whole cell, and 98%/94% between slices 75 and 225, as the central slices of the nucleus are more regular than those on the extremes. Mean Hausdorff distances were 9/17 pixels for the whole cells and 4/13 pixels for central slices. One slice was processed in approximately 8 s and a whole cell in 40 min. The algorithm outperformed active contours in both accuracy and time.