MSight: an image analysis software for liquid chromatography-mass spectrometry.

Images obtained from high-throughput mass spectrometry (MS) contain information that remains hidden when looking at a single spectrum at a time. Image processing of liquid chromatography-MS datasets can be extremely useful for quality control, experimental monitoring and knowledge extraction. The importance of imaging in differential analysis of proteomic experiments has already been established through two-dimensional gels and can now be foreseen with MS images. We present MSight, a new software designed to construct and manipulate MS images, as well as to facilitate their analysis and comparison.

A study of the cambial zone and conductive phloem of common beech (Fagus sylvatica L.) using an image analysis method. I. Influence of tree age on the structure

Quantification is a major problem when using histology to study the influence of ecological factors on tree structure. This paper presents a method to prepare and to analyse transverse sections of cambial zone and of conductive phloem in bark samples. The following paper (II) presents the automated measurement procedure. Part I here describes and discusses the preparation method, and the influence of tree age on the observed structure. Highly contrasted images of samples extracted at breast height during dormancy were analysed with an automatic image analyser. Between three young (38 years) and three old (147 years) trees, age-related differences were identified by size and shape parameters, at both cell and tissue levels. In the cambial zone, older trees had larger and more rectangular fusiform initials. In the phloem, sieve tubes were also larger, but their shape did not change and the area for sap conduction was similar in both categories. Nevertheless, alterations were limited, and demanded statistical analysis to be identified and ascertained. The physiological implications of the structural changes are discussed.

A 2D/3D image analysis system to track fluorescently labeled structures in rod-shaped cells: application to measure spindle pole asymmetry during mitosis.

BACKGROUND: The yeast Schizosaccharomyces pombe is frequently used as a model for studying the cell cycle. The cells are rod-shaped and divide by medial fission. The process of cell division, or cytokinesis, is controlled by a network of signaling proteins called the Septation Initiation Network (SIN); SIN proteins associate with the SPBs during nuclear division (mitosis). Some SIN proteins associate with both SPBs early in mitosis, and then display strongly asymmetric signal intensity at the SPBs in late mitosis, just before cytokinesis. This asymmetry is thought to be important for correct regulation of SIN signaling, and coordination of cytokinesis and mitosis. In order to study the dynamics of organelles or large protein complexes such as the spindle pole body (SPB), which have been labeled with a fluorescent protein tag in living cells, a number of the image analysis problems must be solved; the cell outline must be detected automatically, and the position and signal intensity associated with the structures of interest within the cell must be determined. RESULTS: We present a new 2D and 3D image analysis system that permits versatile and robust analysis of motile, fluorescently labeled structures in rod-shaped cells. We have designed an image analysis system that we have implemented as a user-friendly software package allowing the fast and robust image-analysis of large numbers of rod-shaped cells. We have developed new robust algorithms, which we combined with existing methodologies to facilitate fast and accurate analysis. Our software permits the detection and segmentation of rod-shaped cells in either static or dynamic (i.e. time lapse) multi-channel images. It enables tracking of two structures (for example SPBs) in two different image channels. For 2D or 3D static images, the locations of the structures are identified, and then intensity values are extracted together with several quantitative parameters, such as length, width, cell orientation, background fluorescence and the distance between the structures of interest. Furthermore, two kinds of kymographs of the tracked structures can be established, one representing the migration with respect to their relative position, the other representing their individual trajectories inside the cell. This software package, called "RodCellJ", allowed us to analyze a large number of S. pombe cells to understand the rules that govern SIN protein asymmetry. CONCLUSIONS: "RodCell" is freely available to the community as a package of several ImageJ plugins to simultaneously analyze the behavior of a large number of rod-shaped cells in an extensive manner. The integration of different image-processing techniques in a single package, as well as the development of novel algorithms does not only allow to speed up the analysis with respect to the usage of existing tools, but also accounts for higher accuracy. Its utility was demonstrated on both 2D and 3D static and dynamic images to study the septation initiation network of the yeast Schizosaccharomyces pombe. More generally, it can be used in any kind of biological context where fluorescent-protein labeled structures need to be analyzed in rod-shaped cells. AVAILABILITY: RodCellJ is freely available under http://bigwww.epfl.ch/algorithms.html, (after acceptance of the publication).

Pore space analysis of beech wood: The vessel network

Water transport in wood is vital for the survival of trees. With synchrotron radiation X-ray tomographic microscopy (SRXTM), it has become possible to characterize and quantify the three-dimensional (3D) network formed by vessels that are responsible for longitudinal transport. In the present study, the spatial size dependence of vessels and the organization inside single growth rings in terms of vessel-induced porosity was studied by SRXTM. Network characteristics, such as connectivity, were deduced by digital image analysis from the processed tomographic data and related to known complex network topologies.

Sexual dimorphism of the human mandible: demonstration by elliptical Fourier analysis.

A new quantitative approach of the mandibular sexual dimorphism, based on computer-aided image analysis and elliptical Fourier analysis of the mandibular outline in lateral view is presented. This method was applied to a series of 117 dentulous mandibles from 69 male and 48 female individuals native of Rhenish countries. Statistical discriminant analysis of the elliptical Fourier harmonics allowed the demonstration of a significant sexual dimorphism in 97.1% of males and 91.7% of females, i.e. in a higher proportion than in previous studies using classical metrical approaches. This original method opens interesting perspectives for increasing the accuracy of sex identification in current anthropological practice and in forensic procedures.

Deoxyribonucleic acid content as an indicator of progression of squamous cell carcinogenesis in the esophagus: comparative analysis on imprint-cytospin and tissue section preparation.

BACKGROUND: The purpose of this study was to explore the potential use of image analysis on tissue sections preparation as a predictive marker of early malignant changes during squamous cell (SC) carcinogenesis in the esophagus. Results of DNA ploidy quantification on formalin-fixed, paraffin-embedded tissue using two different techniques were compared: imprint-cytospin and 6 microm thick tissue sections preparation. METHODS: This retrospective study included 26 surgical specimens of squamous cell carcinoma (SCC) from patients who underwent surgery alone at the Department of Surgery in CHUV Hospital in Lausanne between January 1993 and December 2000. We analyzed 53 samples of healthy tissue, 43 tumors and 7 lymph node metastases. RESULTS: Diploid DNA histogram patterns were observed in all histologically healthy tissues, either distant or proximal to the lesion. Aneuploidy was observed in 34 (79%) of 43 carcinomas, namely 24 (75%) of 32 early squamous cell carcinomas and 10 (91%) of 11 advanced carcinomas. DNA content was similar in the different tumor stages, whether patients presented with single or multiple synchronous tumors. All lymph node metastases had similar DNA content as their primary tumor. CONCLUSIONS: Early malignant changes in the esophagus are associated with alteration in DNA content, and aneuploidy tends to correlate with progression of invasive SCC. A very good correlation between imprint-cytospin and tissue section analysis was observed. Although each method used here showed advantages and disadvantages; tissue sections preparation provided useful information on aberrant cell-cycle regulation and helped select the optimal treatment for the individual patient along with consideration of other clinical parameters.

On the convergence of EM-like algorithms for image segmentation using Markov random fields.

Inference of Markov random field images segmentation models is usually performed using iterative methods which adapt the well-known expectation-maximization (EM) algorithm for independent mixture models. However, some of these adaptations are ad hoc and may turn out numerically unstable. In this paper, we review three EM-like variants for Markov random field segmentation and compare their convergence properties both at the theoretical and practical levels. We specifically advocate a numerical scheme involving asynchronous voxel updating, for which general convergence results can be established. Our experiments on brain tissue classification in magnetic resonance images provide evidence that this algorithm may achieve significantly faster convergence than its competitors while yielding at least as good segmentation results.

Sparse regularization for fiber ODF reconstruction: from the suboptimality of ℓ2 and ℓ1 priors to ℓ0.

Diffusion MRI is a well established imaging modality providing a powerful way to probe the structure of the white matter non-invasively. Despite its potential, the intrinsic long scan times of these sequences have hampered their use in clinical practice. For this reason, a large variety of methods have been recently proposed to shorten the acquisition times. Among them, spherical deconvolution approaches have gained a lot of interest for their ability to reliably recover the intra-voxel fiber configuration with a relatively small number of data samples. To overcome the intrinsic instabilities of deconvolution, these methods use regularization schemes generally based on the assumption that the fiber orientation distribution (FOD) to be recovered in each voxel is sparse. The well known Constrained Spherical Deconvolution (CSD) approach resorts to Tikhonov regularization, based on an ℓ(2)-norm prior, which promotes a weak version of sparsity. Also, in the last few years compressed sensing has been advocated to further accelerate the acquisitions and ℓ(1)-norm minimization is generally employed as a means to promote sparsity in the recovered FODs. In this paper, we provide evidence that the use of an ℓ(1)-norm prior to regularize this class of problems is somewhat inconsistent with the fact that the fiber compartments all sum up to unity. To overcome this ℓ(1) inconsistency while simultaneously exploiting sparsity more optimally than through an ℓ(2) prior, we reformulate the reconstruction problem as a constrained formulation between a data term and a sparsity prior consisting in an explicit bound on the ℓ(0)norm of the FOD, i.e. on the number of fibers. The method has been tested both on synthetic and real data. Experimental results show that the proposed ℓ(0) formulation significantly reduces modeling errors compared to the state-of-the-art ℓ(2) and ℓ(1) regularization approaches.

Mitotic figure counts are significantly overestimated in resection specimens of invasive breast carcinomas.

Several authors have demonstrated an increased number of mitotic figures in breast cancer resection specimen when compared with biopsy material. This has been ascribed to a sampling artifact where biopsies are (i) either too small to allow formal mitotic figure counting or (ii) not necessarily taken form the proliferating tumor periphery. Herein, we propose a different explanation for this phenomenon. Biopsy and resection material of 52 invasive ductal carcinomas was studied. We counted mitotic figures in 10 representative high power fields and quantified MIB-1 immunohistochemistry by visual estimation, counting and image analysis. We found that mitotic figures were elevated by more than three-fold on average in resection specimen over biopsy material from the same tumors (20±6 vs 6±2 mitoses per 10 high power fields, P=0.008), and that this resulted in a relative diminution of post-metaphase figures (anaphase/telophase), which made up 7% of all mitotic figures in biopsies but only 3% in resection specimen (P<0.005). At the same time, the percentages of MIB-1 immunostained tumor cells among total tumor cells were comparable in biopsy and resection material, irrespective of the mode of MIB-1 quantification. Finally, we found no association between the size of the biopsy material and the relative increase of mitotic figures in resection specimen. We propose that the increase in mitotic figures in resection specimen and the significant shift towards metaphase figures is not due to a sampling artifact, but reflects ongoing cell cycle activity in the resected tumor tissue due to fixation delay. The dwindling energy supply will eventually arrest tumor cells in metaphase, where they are readily identified by the diagnostic pathologist. Taken together, we suggest that the rapidly fixed biopsy material better represents true tumor biology and should be privileged as predictive marker of putative response to cytotoxic chemotherapy.

Response characteristics of arsenic-sensitive bioreporters expressing the gfp reporter gene

This paper describes the development of an analytical technique for arsenic analyses that is based on genetically-modified bioreporter bacteria bearing a gene encoding for the production of a green fluorescent protein (gfp). Upon exposure to arsenic (in the aqueous form of arsenite), the bioreporter production of the fluorescent reporter molecule is monitored spectroscopically. We compared the response measured as a function of time and concentration by steady-state fluorimetry (SSF) to that measured by epi-fluorescent microscopy (EFM). SSF is a bulk technique; as such it inherently yields less information, whereas EFM monitors the response of many individual cells simultaneously and data can be processed in terms of population averages or subpopulations. For the bioreporter strain used here, as well as for the literature we cite, the two techniques exhibit similar performance characteristics. The results presented here show that the EFM technique can compete with SSF and shows substantially more promise for future improvement; it is a matter of research interest to develop optimized methods of EFM image analysis and statistical data treatment. EFM is a conduit for understanding the dynamics of individual cell response vs. population response, which is not only a matter of research interest, but is also promising in the practical terms of developing micro-scale analysis.

Early functional differentiation in the chick embryonic disc: interaction between mechanical activity and extracellular matrix

The mechanical behaviour of ectodermal cells in the area opaca and the supracellular organization of fibronectin in the adjacent extracellular matrix were studied in whole chick blastoderms developing in vitro. The pattern of spontaneous mechanical activity and its modification by immunoglobulins against fibronectin were determined using a real-time image-analysis system. The pattern of fibronectin was studied using immunocytochemical techniques. It was found that the ectodermal cells in the area opaca actively develop a radially oriented contraction, which leads to a distension of the area pellucida from which the embryo develops. Abnormally increased tension resulted in perturbations of gastrulation and neurulation. An optimized mechanical equilibrium within the blastoderm seems to be necessary for normal development. Anti-fibronectin antibodies applied to the basal side of the blastoderm led rapidly and reversibly to an increase of tension in the contracted cells. This observation indicates that modifications of the extracellular matrix can be transmitted to cytoskeletal elements within adjacent cells. The extracellular matrix of the area opaca contains fibronectin arranged in radially oriented fibrils. This orientation corresponds to the direction of migration of the mesodermal cells. Interestingly, the radial pattern of fibronectin is found in the regions where the ectodermal cells are contracted and develop radially oriented forces. This observation suggests that the supracellular assembly of the extracellular materials could be influenced by the mechanical activity of adjacent cells. Possible modulations of the supracellular organization of extracellular matrix by other factors, e.g. diffusible metabolites, is also discussed. The presence of characteristically organized extracellular matrix components, of spatially differentiated cell activities and of reciprocal interactions between them makes the young chick blastoderm an excellent system for physiological studies of the coordinated cellular activities that lead to changes in form, complexity and function.

Pathological Reorganization of NMDA Receptors Subunits and Postsynaptic Protein PSD-95 Distribution in Alzheimer's Disease.

In Alzheimer's disease (AD), synaptic alterations play a major role and are often correlated with cognitive changes. In order to better understand synaptic modifications, we compared alterations in NMDA receptors and postsynaptic protein PSD-95 expression in the entorhinal cortex (EC) and frontal cortex (FC; area 9) of AD and control brains. We combined immunohistochemical and image analysis methods to quantify on consecutive sections the distribution of PSD-95 and NMDA receptors GluN1, GluN2A and GluN2B in EC and FC from 25 AD and control cases. The density of stained receptors was analyzed using multivariate statistical methods to assess the effect of neurodegeneration. In both regions, the number of neuronal profiles immunostained for GluN1 receptors subunit and PSD-95 protein was significantly increased in AD compared to controls (3-6 fold), while the number of neuronal profiles stained for GluN2A and GluN2B receptors subunits was on the contrary decreased (3-4 fold). The increase in marked neuronal profiles was more prominent in a cortical band corresponding to layers 3 to 5 with large pyramidal cells. Neurons positive for GluN1 or PSD-95 staining were often found in the same localization on consecutive sections and they were also reactive for the anti-tau antibody AD2, indicating a neurodegenerative process. Differences in the density of immunoreactive puncta representing neuropile were not statistically significant. Altogether these data indicate that GluN1 and PSD-95 accumulate in the neuronal perikarya, but this is not the case for GluN2A and GluN2B, while the neuropile compartment is less subject to modifications. Thus, important variations in the pattern of distribution of the NMDA receptors subunits and PSD-95 represent a marker in AD and by impairing the neuronal network, contribute to functional deterioration.

Mapping the proteome of Leishmania Viannia parasites using two-dimensional polyacrylamide gel electrophoresis and associated technologies.

In this study we have demonstrated the potential of two-dimensional electrophoresis (2DE)-based technologies as tools for characterization of the Leishmania proteome (the expressed protein complement of the genome). Standardized neutral range (pH 5-7) proteome maps of Leishmania (Viannia) guyanensis and Leishmania (Viannia) panamensis promastigotes were reproducibly generated by 2DE of soluble parasite extracts, which were prepared using lysis buffer containing urea and nonidet P-40 detergent. The Coomassie blue and silver nitrate staining systems both yielded good resolution and representation of protein spots, enabling the detection of approximately 800 and 1,500 distinct proteins, respectively. Several reference protein spots common to the proteomes of all parasite species/strains studied were isolated and identified by peptide mass spectrometry (LC-ES-MS/MS), and bioinformatics approaches as members of the heat shock protein family, ribosomal protein S12, kinetoplast membrane protein 11 and a hypothetical Leishmania-specific 13 kDa protein of unknown function. Immunoblotting of Leishmania protein maps using a monoclonal antibody resulted in the specific detection of the 81.4 kDa and 77.5 kDa subunits of paraflagellar rod proteins 1 and 2, respectively. Moreover, differences in protein expression profiles between distinct parasite clones were reproducibly detected through comparative proteome analyses of paired maps using image analysis software. These data illustrate the resolving power of 2DE-based proteome analysis. The production and basic characterization of good quality Leishmania proteome maps provides an essential first step towards comparative protein expression studies aimed at identifying the molecular determinants of parasite drug resistance and virulence, as well as discovering new drug and vaccine targets.

3D Automated Lung Nodule Segmentation in HRCT

A fully-automated 3D image analysis method is proposed to segment lung nodules in HRCT. A specific gray-level mathematical morphology operator, the SMDC-connection cost, acting in the 3D space of the thorax volume is defined in order to discriminate lung nodules from other dense (vascular) structures. Applied to clinical data concerning patients with pulmonary carcinoma, the proposed method detects isolated, juxtavascular and peripheral nodules with sizes ranging from 2 to 20 mm diameter. The segmentation accuracy was objectively evaluated on real and simulated nodules. The method showed a sensitivity and a specificity ranging from 85% to 97% and from 90% to 98%, respectively.

Do we truly see what we think we see? The role of cognitive bias in pathological interpretation.

In the histomorphological grading of prostate carcinoma, pathologists have regularly assigned comparable scores for the architectural Gleason and the now-obsolete nuclear World Health Organization (WHO) grading systems. Although both systems demonstrate good correspondence between grade and survival, they are based on fundamentally different biological criteria. We tested the hypothesis that this apparent concurrence between the two grading systems originates from an interpretation bias in the minds of diagnostic pathologists, rather than reflecting a biological reality. Three pathologists graded 178 prostatectomy specimens, assigning Gleason and WHO scores on glass slides and on digital images of nuclei isolated out of their architectural context. The results were analysed with respect to interdependencies among the grading systems, to tumour recurrence (PSA relapse > 0.1 ng/ml at 48 months) and robust nuclear morphometry, as assessed by computer-assisted image analysis. WHO and Gleason grades were strongly correlated (r = 0.82) and demonstrated identical prognostic power. However, WHO grades correlated poorly with nuclear morphology (r = 0.19). Grading of nuclei isolated out of their architectural context significantly improved accuracy for nuclear morphology (r = 0.55), but the prognostic power was virtually lost. In conclusion, the architectural organization of a tumour, which the pathologist cannot avoid noticing during initial slide viewing at low magnification, unwittingly influences the subsequent nuclear grade assignment. In our study, the prognostic power of the WHO grading system was dependent on visual assessment of tumour growth pattern. We demonstrate for the first time the influence a cognitive bias can have in the generation of an error in diagnostic pathology and highlight a considerable problem in histopathological tumour grading.