A statistical approach to computer processing of cryo-electron microscope images: virion classification and 3-D reconstruction

The scattering density of the virus is represented as a truncated weighted sum of orthonormal basis functions in spherical coordinates, where the angular dependence of each basis function has icosahedral symmetry. A statistical model of the image formation process is proposed and the maximum likelihood estimation method computed by an expectation-maximization algorithm is used to estimate the weights in the sum and thereby compute a 3-D reconstruction of the virus particle. If multiple types of virus particle are represented in the boxed images then multiple 3-D reconstructions are computed simultaneously without first requiring that the type of particle shown in each boxed image be determined. Examples of the procedure are described for viruses with known structure: (1). 3-D reconstruction of Flockhouse Virus from experimental images, (2). 3-D reconstruction of the capsid of Nudaurelia Omega Capensis Virus from synthetic images, and (3). 3-D reconstruction of both the capsid and the procapsid of Nudaurelia Omega Capensis Virus from a mixture of unclassified synthetic images.     

A visual analytics approach for models of heterogeneous cell populations

In recent years, cell population models have become increasingly common. In contrast to classic single cell models, population models allow for the study of cell-to-cell variability, a crucial phenomenon in most populations of primary cells, cancer cells, and stem cells. Unfortunately, tools for in-depth analysis of population models are still missing. This problem originates from the complexity of population models. Particularly important are methods to determine the source of heterogeneity (e.g., genetics or epigenetic differences) and to select potential (bio-)markers. We propose an analysis based on visual analytics to tackle this problem. Our approach combines parallel-coordinates plots, used for a visual assessment of the high-dimensional dependencies, and nonlinear support vector machines, for the quantification of effects. The method can be employed to study qualitative and quantitative differences among cells. To illustrate the different components, we perform a case study using the proapoptotic signal transduction pathway involved in cellular apoptosis.     

A coarse-to-fine approach to prostate boundary segmentation in ultrasound images

Background  

In this paper a novel method for prostate segmentation in transrectal ultrasound images is presented.     

3D surface analysis and classification in neuroimaging segmentation

This work emphasizes new algorithms for 3D edge and corner detection used in surface extraction and new concept of image segmentation in neuroimaging based on multidimensional shape analysis and classification. We propose using of NifTI standard for describing input data which enables interoperability and enhancement of existing computing tools used widely in neuroimaging research. In methods section we present our newly developed algorithm for 3D edge and corner detection, together with the algorithm for estimating local 3D shape. Surface of estimated shape is analyzed and segmented according to kernel shapes.     

Automated image analysis methods for 3-D quantification of the neurovascular unit from multichannel confocal microscope images.

BACKGROUND: There is a need for integrative and quantitative methods to investigate the structural and functional relations among elements of complex systems, such as the neurovascular unit (NVU), that involve multiple cell types, microvasculatures, and various genomic/proteomic/ionic functional entities. METHODS: Vascular casting and selective labeling enabled simultaneous three-dimensional imaging of the microvasculature, cell nuclei, and cytoplasmic stains. Multidimensional segmentation was achieved by (i) bleed-through removal and attenuation correction; (ii) independent segmentation and morphometry for each corrected channel; and (iii) spatially associative feature computation across channels. The combined measurements enabled cell classification based on nuclear morphometry, cytoplasmic signals, and distance from vascular elements. Specific spatial relations among the NVU elements could be quantified. RESULTS: A software system combining nuclear and vessel segmentation codes and associative features was constructed and validated. Biological variability contributed to misidentified nuclei (9.3%), undersegmentation of nuclei (3.7%), hypersegmentation of nuclei (14%), and missed nuclei (4.7%). Microvessel segmentation errors occurred rarely, mainly due to nonuniform lumen staining. CONCLUSIONS: Associative features across fluorescence channels, in combination with standard features, enable integrative structural and functional analysis of the NVU. By labeling additional structural and functional entities, this method can be scaled up to larger-scale systems biology studies that integrate spatial and molecular information.     

A collaborative framework for 3D alignment and classification of heterogeneous subvolumes in cryo-electron tomography

The limitation of using low electron doses in non-destructive cryo-electron tomography of biological specimens can be partially offset via averaging of aligned and structurally homogeneous subsets present in tomograms. This type of sub-volume averaging is especially challenging when multiple species are present. Here, we tackle the problem of conformational separation and alignment with a “collaborative” approach designed to reduce the effect of the “curse of dimensionality” encountered in standard pair-wise comparisons. Our new approach is based on using the nuclear norm as a collaborative similarity measure for alignment of sub-volumes, and by exploiting the presence of symmetry early in the processing. We provide a strict validation of this method by analyzing mixtures of intact simian immunodeficiency viruses SIV mac239 and SIV CP-MAC. Electron microscopic images of these two virus preparations are indistinguishable except for subtle differences in conformation of the envelope glycoproteins displayed on the surface of each virus particle. By using the nuclear norm-based, collaborative alignment method presented here, we demonstrate that the genetic identity of each virus particle present in the mixture can be assigned based solely on the structural information derived from single envelope glycoproteins displayed on the virus surface.     

A multiresolution approach to automated classification of protein subcellular location images

Background  

Fluorescence microscopy is widely used to determine the subcellular location of proteins. Efforts to determine location on a proteome-wide basis create a need for automated methods to analyze the resulting images. Over the past ten years, the feasibility of using machine learning methods to recognize all major subcellular location patterns has been convincingly demonstrated, using diverse feature sets and classifiers. On a well-studied data set of 2D HeLa single-cell images, the best performance to date, 91.5%, was obtained by including a set of multiresolution features. This demonstrates the value of multiresolution approaches to this important problem.     

Localization of telomeres in plant interphase nuclei by in situ hybridization and 3D confocal microscopy

We investigated the three-dimensional (3D) arrangement of telomeres in structurally well preserved, interphase nuclei of Pisum stativum and Vicia faba root tips using in situ hybridization of a probe to telomeric sequences. The probe was labelled with either digoxygenin or biotin and hybridized sequences were detected by immunofluorescence. Three-dimensional data sets were collected by confocal optical microscopy or using a cooled CCD camera. Twelve stacks of optical sections of P. sativum nuclei and nine of V. faba nuclei were studied in detail. Projections through the stacks of optical sections revealed that, in both species, most of the telomeres were adjacent to the nuclear envelope except for a small number next to the nucleolar periphery. In V. faba nuclei, the telomeres were clearly clustered at one pole while in P. sativum there was only a slight tendency for clustering. In V. faba, clusters were found at opposite poles in pairs of sister nuclei rather than at adjacent poles as would be expected if the arrangement at telophase were maintained into interphase.by D. Bazett-Jones     

Automatic segmentation and supervised learning-based selection of nuclei in cancer tissue images

Analysis of preferential localization of certain genes within the cell nuclei is emerging as a new technique for the diagnosis of breast cancer. Quantitation requires accurate segmentation of 100-200 cell nuclei in each tissue section to draw a statistically significant result. Thus, for large-scale analysis, manual processing is too time consuming and subjective. Fortuitously, acquired images generally contain many more nuclei than are needed for analysis. Therefore, we developed an integrated workflow that selects, following automatic segmentation, a subpopulation of accurately delineated nuclei for positioning of fluorescence in situ hybridization-labeled genes of interest. Segmentation was performed by a multistage watershed-based algorithm and screening by an artificial neural network-based pattern recognition engine. The performance of the workflow was quantified in terms of the fraction of automatically selected nuclei that were visually confirmed as well segmented and by the boundary accuracy of the well-segmented nuclei relative to a 2D dynamic programming-based reference segmentation method. Application of the method was demonstrated for discriminating normal and cancerous breast tissue sections based on the differential positioning of the HES5 gene. Automatic results agreed with manual analysis in 11 out of 14 cancers, all four normal cases, and all five noncancerous breast disease cases, thus showing the accuracy and robustness of the proposed approach. ? Published 2012 Wiley Periodicals, Inc.     

Three-dimensional component labeling of digital confocal microscope images enumerates replication centers in BrdUrd labeled fibroblasts.

An algorithm (3Dlabel-1) for labeling three dimensional binary data sets has been developed. Serial optical sections are acquired using a confocal laser scanning microscope. After filtering and thresholding operations, contiguous elements in the three dimensional data sets are identified and labeled. The results are used to calculate the number and size of objects. Additionally, the labeled data are displayed by applying an algorithm (3Ddisp-1) to generate stereo pairs, in which label numbers are color coded and depth is cued by intensity. These procedures have been applied in investigations of the temporal and spatial distribution of replication centers throughout S phase in BrdUrd pulse labeled mouse fibroblasts.     

A cytokinetic model for heterogeneous mammalian cell populations. I. Cell growth and cell death

A general model is proposed for describing the growth behavior of mammalian cell populations, which features:(a) a cell cycle time distribution function with properties such that mean and variance increase with increasing population size; (b) maturation age and maturation rate functions which constrain the maturational pathways of individual cells; and (c) a death rate function, where cell death is construed as irreparable damage to a cell's reproductive apparatus. The biological implications of the model are discussed, and methods for relating the model to real cell systems by means of commonly used experimental techniques are described. The model is compared with earlier models.     

Semi-automatic 3D segmentation of carotid lumen in contrast-enhanced computed tomography angiography images

The atherosclerosis disease is one of the major causes of the death in the world. Atherosclerosis refers to the hardening and narrowing of the arteries by plaques. Carotid stenosis is a narrowing or constriction of carotid artery lumen usually caused by atherosclerosis. Carotid artery stenosis can increase risk of brain stroke. Contrast-enhanced Computed Tomography Angiography (CTA) is a minimally invasive method for imaging and quantification of the carotid plaques. Manual segmentation of carotid lumen in CTA images is a tedious and time consuming procedure which is subjected to observer variability. As a result, there is a strong and growing demand for developing computer-aided carotid segmentation procedures. In this study, a novel method is presented for carotid artery lumen segmentation in CTA data. First, the mean shift smoothing is used for uniformity enhancement of gray levels. Then with the help of three seed points, the centerlines of the arteries are extracted by a 3D Hessian based fast marching shortest path algorithm. Finally, a 3D Level set function is performed for segmentation. Results on 14 CTA volumes data show 85% of Dice similarity and 0.42 mm of mean absolute surface distance measures. Evaluation shows that the proposed method requires minimal user intervention, low dependence to gray levels changes in artery path, resistance to extreme changes in carotid diameter and carotid branch locations. The proposed method has high accuracy and can be used in qualitative and quantitative evaluation.     

A nonparametric approach to the discrimination of two cell populations

An approach to the statistical testing of differences between two cell populations the elements of which are characterized by multivariate data is described. The approach is based on the Fisher discriminant and the Kruskal-Wallis test. The method, which is sufficient but not necessary, makes no assumptions about the normality of the data or about the equality of the covariance matrices for each population.     

Segmentation of center brains and optic lobes in 3D confocal images of adult fruit fly brains

Automatic alignment (registration) of 3D images of adult fruit fly brains is often influenced by the significant displacement of the relative locations of the two optic lobes (OLs) and the center brain (CB). In one of our ongoing efforts to produce a better image alignment pipeline of adult fruit fly brains, we consider separating CB and OLs and align them independently. This paper reports our automatic method to segregate CB and OLs, in particular under conditions where the signal to noise ratio (SNR) is low, the variation of the image intensity is big, and the relative displacement of OLs and CB is substantial.We design an algorithm to find a minimum-cost 3D surface in a 3D image stack to best separate an OL (of one side, either left or right) from CB. This surface is defined as an aggregation of the respective minimum-cost curves detected in each individual 2D image slice. Each curve is defined by a list of control points that best segregate OL and CB. To obtain the locations of these control points, we derive an energy function that includes an image energy term defined by local pixel intensities and two internal energy terms that constrain the curve’s smoothness and length. Gradient descent method is used to optimize this energy function. To improve both the speed and robustness of the method, for each stack, the locations of optimized control points in a slice are taken as the initialization prior for the next slice. We have tested this approach on simulated and real 3D fly brain image stacks and demonstrated that this method can reasonably segregate OLs from CBs despite the aforementioned difficulties.     

A graphical model approach to automated classification of protein subcellular location patterns in multi-cell images

Background  

Knowledge of the subcellular location of a protein is critical to understanding how that protein works in a cell. This location is frequently determined by the interpretation of fluorescence microscope images. In recent years, automated systems have been developed for consistent and objective interpretation of such images so that the protein pattern in a single cell can be assigned to a known location category. While these systems perform with nearly perfect accuracy for single cell images of all major subcellular structures, their ability to distinguish subpatterns of an organelle (such as two Golgi proteins) is not perfect. Our goal in the work described here was to improve the ability of an automated system to decide which of two similar patterns is present in a field of cells by considering more than one cell at a time. Since cells displaying the same location pattern are often clustered together, considering multiple cells may be expected to improve discrimination between similar patterns.     

Detection of nuclei in 4D Nomarski DIC microscope images of early Caenorhabditis elegans embryos using local image entropy and object tracking

Background  

The ability to detect nuclei in embryos is essential for studying the development of multicellular organisms. A system of automated nuclear detection has already been tested on a set of four-dimensional (4D) Nomarski differential interference contrast (DIC) microscope images of Caenorhabditis elegans embryos. However, the system needed laborious hand-tuning of its parameters every time a new image set was used. It could not detect nuclei in the process of cell division, and could detect nuclei only from the two- to eight-cell stages.     

A novel approach to the study of kinetically changing cell populations

A method is described for estimating changes in cell cycle times during periods of rapid change in proliferation rate. This method, which depends upon the interpretation of pre- and post-velocity sedimentation fractionation continuous thymidine labelling patterns, exploits the relationship between sedimentation rate and cell cycle location. By this means, cycle times can be estimated under conditions that are difficult (if not impossible) to analyse by FLM methods.