An improved pixel-based approach for analyzing images in two-dimensional gel electrophoresis

An improved pixel-based approach for analyzing 2-DE images is presented. The key feature of the method is to create a mask based on all gels in the experiment using image morphology, followed by multivariate analysis on the pixel level. The method reduces the impact of noise and background by identifying regions in the image where protein spots are present, but make no assumption on individual spot boundaries for isolated spots. This makes it possible to detect significant changes in complex regions, and visualize these changes over multiple gels in an easy way. False missing values and spot volumes caused by imposing erroneous spot boundaries are thus circumvented. The approach presented gives improved pixel-based information from the gels, and is also an alternative to existing methods for data-reduction, significance testing and visualization of 2-DE data. Results are compared with software using a common spot boundary approach on an experiment consisting of 35 full size gel images. Gel alignment is required before analysis.     

A review of the CLIP system for the quantitative analysis of two-dimensional electrophoresis gels

This paper reviews the CLIP image processing system for the complete analysis of two-dimensional electrophoresis images. The analysis problem for two-dimensional gel images can be broken down into three issues: segmentation of individual gel images, alignment and comparison of pairs of gel images, and information storage and retrieval. This paper describes these problems and reviews how the CLIP system handles each of them. Segmentation is the location and isolation of each protein spot on an individual gel image and also the extraction of individual spot data such as position, area and volume. There are three basic stages: background field correction, noise filtering, spot detection and information extraction. Alignment and comparison of gel images involves matching protein spots between two gels. This can be quite difficult because there is not a simple relationship which can transform one gel image onto another. The database issues concern storing all the information which has been obtained from the above operations such that retrieval of this information can be readily performed. The advantage of the CLIP system over others is speed of processing. CLIP series computers use one processor for every pixel of the camera image such that image processing algorithms run in parallel. The main disadvantage is in the cost of these machines. With the declining trend in the cost of parallel processors, these machines will become more and more viable alternatives. This papers reviews the algorithms for the analysis of two-dimensional gels. It is shown that CLIP is flexible enough to perform more than one type of algorithm for a particular operation.     

Adaptive contrast enhancement of two-dimensional electrophoretic protein gel images facilitates visualization, orientation and alignment

2-DE is a powerful technique to discriminate post-translationally modified protein isoforms. However, all steps of 2-DE preparation and gel-staining may introduce unwanted artefacts, including inconsistent variation of background intensity over the entire 2-DE gel image. Background intensity variations limit the accuracy of gel orientation, overlay alignment and spot detection methods. We present a compact and efficient denoising algorithm that adaptively enhances the image contrast and then, through thresholding and median filtering, removes the gray-scale range covering the background. Applicability of the algorithm is demonstrated on immunoblots, isotope-labeled gels, and protein-stained gels. Validation is performed in contexts of (i) automatic gel orientation based on Hough transformation, (ii) overlay alignment based on cross correlation and (iii) spot detection. In gel stains with low background variability, e.g. Sypro Ruby, denoising will lower the spot detection sensitivity. In gel regions with high background levels denoising enhances spot detection. We propose that the denoising algorithm prepares images with high background for further automatic analysis, without requiring manual input on a gel-to-gel basis.     

A new method for assigning common spot boundaries for multiple gels in two-dimensional gel electrophoresis

The benefits of defining common spot boundaries when several gels from 2-DE are compared and analyzed have lately been stressed by both commercial software producers and users of this software. Though the importance of common spot boundaries is clearly stated, few reports exist that target this issue explicitly. In this study a method for defining common spots boundaries is developed, called the spot density method. The method consists of the following steps: segmentation and spot identification on each individual gel, transferring the spot-center coordinates for all gels onto a single new gel, collecting spot centers clustered together in the new gel and finally assigning pixels and new spot boundaries based on the spots in each cluster. The method is compared to a synthetic gel approach, and validated by visual inspection of three representative areas in the gels. The gel images need to be aligned prior to segmentation and spot identification, but the method can be used regardless of the choice of segmentation procedure. This makes the method an easy extension to existing methods for spot identification and matching. Conclusions based on the visual inspection are that the spot density method identifies partly overlapping spots and low-intensity spots better than the synthetic gel approach.     

Exploratory data analysis groupware for qualitative and quantitative electrophoretic gel analysis over the Internet-WebGel

Many scientists use quantitative measurements to compare the presence and amount, of various proteins and nucleotides among series of one- and two-dimensional (1-D and 2-D) electrophoretic gels. These gels are often scanned into digital image files. Gel spots are then quantified using stand-alone analysis software. However, as more research collaborations take place over the Internet, it has become useful to share intermediate quantitative data between researchers. This allows research group members to investigate their data and share their work in progress. We developed a World Wide Web group-accessible software system, WebGel, for interactively exploring qualitative and quantitative differences between electrophoretic gels. Such Internet databases are useful for publishing quantitative data and allow other researchers to explore the data with respect to their own research. Because intermediate results of one user may be shared with their collaborators using WebGel, this form of active data-sharing constitutes a groupware method for enhancing collaborative research. Quantitative and image gel data from a stand-alone gel image processing system are copied to a database accessible on the WebGel Web server. These data are then available for analysis by the WebGel database program residing on that server. Visualization is critical for better understanding of the data. WebGel helps organize labeled gel images into montages of corresponding spots as seen in these different gels. Various views of multiple gel images, including sets of spots, normalization spots, labeled spots, segmented gels, etc. may also be displayed. These displays are active and may be used for performing database operations directly on individual protein spots by simply clicking on them. Corresponding regions between sets of gels may be visually analyzed using Flicker-comparison (Electrophoresis 1997, 18, 122-140) as one of the WebGel methods for qualitative analysis. Quantitative exploratory data analysis can be performed by comparing protein concentration values between corresponding spots for multiple samples run in separate gels. These data are then used to generate reports on statistical differences between sets of gels (e.g., between different disease states such as benign or metastatic cancers, etc.). Using combined visual and quantitative methods, WebGel can help bridge the analysis of dissimilar gels which are difficult to analyze with stand-alone systems and can serve as a collaborative Internet tool in a groupware setting.     

The inter‐ and intra‐operator variability in manual spot segmentation and its effect on spot quantitation in two‐dimensional electrophoresis analysis

Separation of complex mixtures of proteins by 2‐DE is a fundamental component of current proteomic technology. Quantitative analysis of the images generated by digitization of such gels is critical for identifying alterations in protein expression within a given biological system. Software packages are designed for this purpose. The accurate definition of protein spot boundaries, using a suitable method of image segmentation, is a key requirement for image analysis. It is often necessary for operators to intervene manually to correct mistakes in spot segmentation; therefore operator subjectivity and differences in ability can weaken the analysis. We estimated the error in spot quantification after manual spot segmentation, which was performed by different operators, using two different software packages. Our results clearly show that this operation was associated with significant inter‐ and intra‐variability and an overestimation of subsequent spot intensity, especially when spots were weak. For comparative studies, we suggest separately analysing spots which have been manually segmented by imposing a requirement for at least a threefold difference in spot intensity in addition to use of statistical tests.     

An image analysis suite for spot detection and spot matching in two-dimensional electrophoresis gels

We propose a suite of novel algorithms for image analysis of protein expression images obtained from 2-D electrophoresis. These algorithms are a segmentation algorithm for protein spot identification, and an algorithm for matching protein spots from two corresponding images for differential expression study. The proposed segmentation algorithm employs the watershed transformation, k-means analysis, and distance transform to locate the centroids and to extract the regions of the proteins spots. The proposed spot matching algorithm is an integration of the hierarchical-based and optimization-based methods. The hierarchical method is first used to find corresponding pairs of protein spots satisfying the local cross-correlation and overlapping constraints. The matching energy function based on local structure similarity, image similarity, and spatial constraints is then formulated and optimized. Our new algorithm suite has been extensively tested on synthetic and actual 2-D gel images from various biological experiments, and in quantitative comparisons with ImageMaster2D Platinum the proposed algorithms exhibit better spot detection and spot matching.     

Optimization of large gel 2D electrophoresis for proteomic studies of skeletal muscle

We describe improved methods for large format, two-dimensional gel electrophoresis (2DE) that improve protein solubility and recovery, minimize proteolysis, and reduce the loss of resolution due to contaminants and manipulations of the gels, and thus enhance quantitative analysis of protein spots. Key modifications are: (i) the use of 7 M urea and 2 M thiourea, instead of 9 M urea, in sample preparation and in the tops of the gel tubes; (ii) standardized deionization of all solutions containing urea with a mixed bed ion exchange resin and removal of urea from the electrode solutions; and (iii) use of a new gel tank and cooling device that eliminate the need to run two separating gels in the SDS dimension. These changes make 2DE analysis more reproducible and sensitive, with minimal artifacts. Application of this method to the soluble fraction of muscle tissues reliably resolves ~1800 protein spots in adult human skeletal muscle and over 2800 spots in myotubes.     

Superimposing two-dimensional gels to study genetic variation in malaria parasites.

Two-dimensional polyacrylamide gel electrophoresis is a valuable tool for studying genetic variation in the human malaria parasite, Plasmodium falciparum. It involves examining the position of protein spots in gel produced from different isolates. Some spots have been seen to vary, while others have had a constant position in all isolates so far examined. These invariant spots provide a reference frame to compare variations in other spots. This paper discusses the usefulness of digital image handling, warping and superimposition in a personal computer environment. Rather than produce a fully automatic interpretation system, we show how the computer may be used as a tool for manipulating gel images, although interpretation of the gels' features remains with the human expert. Autoradiographs are scanned on a desktop scanner, and the images in digital form can be displayed on a monitor attached to a personal computer. The coordinates of the invariant spots on each of several gels are identified by the user. Each of the gels is then warped so that the invariant spots of all the gels coincide as closely as possible. The variable spots are then examined. We have used both affine warping transformations, which match the invariant spots as closely as possible, and thin plate spline transformations, which match them exactly. Colour superimposition proved a useful way of examining the gels.     

Determination of protein spots separated by two-dimensional polyacrylamide gel electrophoresis

A simple method for quantitating proteins in the spots on two-dimensional polyacrylamide gel electropherograms is described. The system consists in three steps: (1) O'Farrell's two-dimensional gel electrophoresis of the proteins to be analysed; (2) staining of the gels with Coomassie brilliant blue; and (3) determination of the area and integrated density of the stained spots by the Joyce Loebl Magiscan-1 image analysis system. The method can be used for the determination of proteins in the range 0.5-100 micrograms/cm2; the amount of protein involved in most spots detected by the staining method actually falls within this range. As the minimum spot diameter that can easily be handled by the method is about 2 mm, as much as 30 ng of protein in such a spot can be determined. The method can also be applied to autoradiograms.     

Functional characterization of two-dimensional gel-separated proteins using sequential staining

Proteins separated by two-dimensional (2-D) gel electrophoresis can be visualized using various protein staining methods. This is followed by downstream procedures, such as image analysis, gel spot cutting, protein digestion, and mass spectrometry (MS), to characterize protein expression profiles within cells, tissues, organisms, or body fluids. Characterizing specific post-translational modifications on proteins using MS of peptide fragments is difficult and labor-intensive. Recently, specific staining methods have been developed and merged into the 2-D gel platform so that not only general protein patterns but also patterns of phosphorylated and glycosylated proteins can be obtained. We used the new Pro-Q Diamond phosphoprotein dye technology for the fluorescent detection of phosphoproteins directly in 2-D gels of mouse leukocyte proteins, and Pro-Q Emerald 488 glycoprotein dye to detect glycoproteins. These two fluorescent stains are compatible with general protein stains, such as SYPRO Ruby stain. We devised a sequential procedure using Pro-Q Diamond (phosphoprotein), followed by Pro-Q Emerald 488 (glycoprotein), followed by SYPRO Ruby stain (general protein stain), and finally silver stain for total protein profile. This multiple staining of the proteins in a single gel provided parallel determination of protein expression and preliminary characterization of post-translational modifications of proteins in individual spots on 2-D gels. Although this method does not provide the same degree of certainty as traditional MS methods of characterizing post-translational modifications, it is much simpler, faster, and does not require sophisticated equipment and expertise in MS.     

Microscale solution IEF combined with 2-D DIGE substantially enhances analysis depth of complex proteomes such as mammalian cell and tissue extracts

Current gel-based protein profiling methods such as 2-DE and fluorescent 2-D difference in gel electrophoresis (DIGE) evaluate small portions of complex proteomes. Hence, sample prefractionation is essential for more comprehensive proteome coverage and detection of low-abundant proteins. In this study, we describe the combination of DIGE labeling with microscale solution IEF (MicroSol-IEF) fractionation and subsequent analysis on slightly overlapping narrow pH range 2-D gels. By fluorescently tagging and mixing samples and controls prior to prefractionation, complications resulting from minor run-to-run variations during MicroSol-IEF separations of multiple samples are avoided. This greatly improves the reliability of quantitative comparisons. To illustrate its utility, this 3-D DIGE strategy was applied to analysis of human melanoma cells and mouse lung tissue extracts. Approximately 1000 reproducible spots can be obtained from narrow range 2-D gels of individual MicroSol-IEF fractions, and approximately 6000 spots can be obtained from entire proteomes. Quantitative changes in closely related samples could be more reliably detected and the method has a greatly increased capacity to distinguish between closely related protein isoforms. Thus the 3-D DIGE strategy produces a powerful method for more comprehensive and more reliable quantitative comparisons of protein profiles of very complex proteomes.     

The scaled volume as an image analysis variable for detecting changes in protein expression levels by silver stain.

A new variable for measuring the relative intensities of silver stained protein spots on two-dimensional gels is described. The scaled volume (SV) more accurately measures the intensity of protein spots and accounts for differences frequently encountered when trying to compare two gels than other variables such as relative volume ratio, optical density, or relative optical density. The SV scales the signal of interest by the noise (gel background) with secondary signals removed (spots not of interest, e.g., technical artifacts). The SV of spot intensities offers a better dynamic response to protein amount for the model proteins studied here. Depending on the quantity of protein loaded onto gels, we have observed a coefficient of variation range of 0.2 to 1.3. We also observe that the SV silver stain response follows a characteristic exponential profile for different proteins.     

中国小型猪心肌梗死模型的比较蛋白质组学研究

利用二维电泳(2DE)分离中国小型猪心肌梗死模型的正常与梗死心肌组织的蛋白提取液, 采用 PDQuest 软件对比分析了两种心肌组织在pH=5─8范围内的2DE谱图. 正常心肌组织检出851个蛋白点, 梗死组织检出1 032个蛋白点. 发现13个蛋白质点只在小型猪的正常心肌组织中表达, 而有14个蛋白质点只在梗死心肌组织中表达. 另外, 还有49个蛋白点在两种组织中表达量上有显著性变化(P<0.05), 选择进行质谱分析其中11个蛋白点, 成功地鉴定出7种蛋白, 蛋白功能分析结果表明, 这些蛋白的差异表达与心肌梗死过程相关.     

Identification of proteins in human cerebrospinal fluid using liquid-phase isoelectric focusing as a prefractionation step followed by two-dimensional gel electrophoresis and matrix-assisted laser desorption/ionisation mass spectrometry

Cerebrospinal fluid (CSF) is in close proximity to the brain and changes in the protein composition of CSF may be indicative of altered brain protein expression in neurodegenerative disorders. Analysis of brain-specific proteins in CSF is complicated by the fact that most CSF proteins are derived from the plasma and tend to obscure less abundant proteins. By adopting a prefractionation step prior to two-dimensional gel electrophoresis (2-DE), less abundant proteins are enriched and can be detected in complex proteomes such as CSF. We have developed a method in which liquid-phase isoelectric focusing (IEF) is used to prefractionate individual CSF samples; selected IEF fractions are then analysed on SYPRO-Ruby-stained 2-D gels, with final protein identification by matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOFMS). To optimise the focusing of the protein spots on the 2-D gel, the ampholyte concentration in liquid-phase IEF was minimised and the focusing time in the first dimension was increased. When comparing 2-D gels from individual prefractionated and unfractionated CSF samples it is evident that individual protein spots are larger and contain more protein after prefractionation of CSF. Generally, more protein spots were also detected in the 2-D gels from prefractionated CSF compared with direct 2-DE separations of CSF. Several proteins, including cystatin C, IgM-kappa, hemopexin, acetyl-coenzyme A carboxylase-alpha, and alpha-1-acid glycoprotein, were identified in prefractionated CSF but not in unfractionated CSF. Low abundant forms of posttranslationally modified proteins, e.g. alpha-1-acid glycoprotein and alpha-2-HS glycoprotein, can be enriched, thus better resolved and detected on the 2-D gel. Liquid-phase IEF, as a prefractionation step prior to 2-DE, reduce sample complexity, facilitate detection of less abundant protein components, increases the protein loads and the protein amount in each gel spot for MALDI-MS analysis.     

Assignment of human plasma polypeptides on a nondenaturing 2-D gel using MALDI-MS and PMF and comparisons with the results of intact protein mapping

Human plasma proteins were separated by 2-DE under nondenaturing conditions followed by the assignment of the CBB-stained spots using MALDI-MS and PMF, aiming to correlate the information of intact proteins with that of constituent polypeptides. A microgel system was employed to facilitate the analysis. Totally 157 spots on a nondenaturing micro-2-DE gel were numbered, the spots were excised, the proteins in the gel pieces were subjected to in-gel digestion with trypsin followed by polypeptide analysis using MALDI-MS and PMF. Two PMF algorithms, MASCOT (with Swiss-Prot database) and ProFound (with NCBInr database) were employed. A total of 153 spots out of the 157 provided significant match (p <0.05) with polypeptides in databases. Eighty spots were assigned to contain multiple (2-4) polypeptides, suggesting (i) noncovalent interaction between proteins/polypeptides, (ii) disulfide bonding of polypeptides, or (iii) overlapping of the protein locations on the gel. The results of polypeptide assignment coincided very well with the results of protein mapping previously reported, in which 33 plasma proteins were identified using blotting-immunochemical staining (Manabe, T., Takahashi, Y., Higuchi, N., Okuyama, T., Electrophoresis 1985, 6, 462-467). Further, 19 polypeptides in 25 spots were newly assigned. These results demonstrate that the techniques of MALDI-MS and PMF can be applied for analysis of proteins separated on nondenaturing 2-DE gels, providing information on their polypeptide structure. The integrated information on proteins and polypeptides would help the comprehensive understanding on the functions of complex protein systems.     

Quantitative and qualitative measure of intralaboratory two-dimensional protein gel reproducibility and the effects of sample preparation, sample load, and image analysis

We investigate one approach to assess the quantitative variability in two-dimensional gel electrophoresis (2-DE) separations based on gel-to-gel variability, sample preparation variability, sample load differences, and the effect of automation on image analysis. We observe that 95% of spots present in three out of four replicate gels exhibit less than a 0.52 coefficient of variation (CV) in fluorescent stain intensity (% volume) for a single sample run on multiple gels. When four parallel sample preparations are performed, this value increases to 0.57. We do not observe any significant change in quantitative value for an increase or decrease in sample load of 30% when using appropriate image analysis variables. Increasing use of automation, while necessary in modern 2-DE experiments, does change the observed level of quantitative and qualitative variability among replicate gels. The number of spots that change qualitatively for a single sample run in parallel varies from a CV = 0.03 for fully manual analysis to CV = 0.20 for a fully automated analysis. We present a systematic method by which a single laboratory can measure gel-to-gel variability using only three gel runs.     

Fast automatic registration of images using the phase of a complex wavelet transform: application to proteome gels

Image registration describes the process of manipulating a distorted version of an image such that its pixels overlay the equivalent pixels in a clean, master or reference image. The need for it has assumed particular prominence in the analysis of images of electrophoretic gels used in the analysis of protein expression levels in living cells, but also has fundamental applications in most other areas of image analysis. Much of the positional information of a data feature is carried in the phase of a complex transform, so a complex transform allows explicit specification of the phase, and hence of the position of features in the image. Registration of a test gel to a reference gel is achieved by using a multiresolution movement map derived from the phase of a complex wavelet transform (the Q-shift wavelet transform) to dictate the warping directly via movement of the nodes of a Delaunay-triangulated mesh of points. This warping map is then applied to the original untransformed image such that the absolute magnitude of the spots remains unchanged. The technique is general to any type of image. Results are presented for a simple computer simulated gel, a simple real gel registration between similar "clean" gels with local warping vectors distributed about one main direction, a hard problem between a reference gel and a "dirty" test gel with multi-directional warping vectors and many artifacts, and some typical gels of present interest in post-genomic biology. The method compares favourably with others, since it is computationally rapid, effective and entirely automatic.     

Quantitative comparison and evaluation of two commercially available,two-dimensional electrophoresis image analysis software packages,Z3 and Melanie

While a variety of software packages are available for analyzing two-dimensional electrophoresis (2-DE) gel images, no comparisons between these packages have been published, making it difficult for end users to determine which package would best meet their needs. The goal here was to develop a set of tests to quantitatively evaluate and then compare two software packages, Melanie 3.0 and Z3, in three of the fundamental steps involved in 2-DE image analysis: (i) spot detection, (ii) gel matching, and (iii) spot quantitation. To test spot detection capability, automatically detected protein spots were compared to manually counted, "real" protein spots. Spot matching efficiency was determined by comparing distorted (both geometrically and nongeometrically) gel images with undistorted original images, and quantitation tests were performed on artificial gels with spots of varying Gaussian volumes. In spot detection tests, Z3 performed better than Melanie 3.0 and required minimal user intervention to detect approximately 89% of the actual protein spots and relatively few extraneous spots. Results from gel matching tests depended on the type of image distortion used. For geometric distortions, Z3 performed better than Melanie 3.0, matching 99% of the spots, even for extreme distortions. For nongeometrical distortions, both Z3 and Melanie 3.0 required user intervention and performed comparably, matching 95% of the spots. In spot quantitation tests, both Z3 and Melanie 3.0 predicted spot volumes relatively well for spot ratios less than 1:6. For higher ratios, Melanie 3.0 did much better. In summary, results suggest Z3 requires less user intervention than Melanie 3.0, thus simplifying differential comparison of 2-DE gel images. Melanie 3.0, however, offers many more optional tools for image editing, spot detection, data reporting and statistical analysis than Z3. All image files used for these tests and updated information on the software are available on the internet (http://www.umbc.edu/proteome), allowing similar testing of other 2-DE image analysis software packages.     

Differential proteome analysis of tonsils from children with chronic tonsillitis or with hyperplasia reveals disease-associated protein expression differences

A proteomic approach has been used to establish a proteome map and differentiate between the protein composition of tonsils from patients with chronic tonsillitis (CT) and that of tonsils with hyperplasia (HPL). Two-dimensional gel analysis was performed with material from four patients with HPL and five patients with CT. An average of approximately 600 spots were detected in each gel. A total of 127 different proteins were identified in 158 spots analyzed by mass spectrometry. Our study revealed disease-associated differences between protein abundance for two protein spots, an HSP27 isoform and UMP-CMP kinase. Both protein spots were more abundant in the CT group. HSP27 ELISA was performed for 32 patients, 12 belonging to the HPL group and 20 to the CT group. ELISA could not be used to differentiate HSP27 isoforms nor to distinguish CT from HPL. HSP27 was found to migrate to two further protein spots in the 2D gels. The differently expressed HSP27 isoform migrated as the most acidic of all the HSP27 isoforms detected, indicating the highest degree of phosphorylation. The sum of all three HSP27 abundances in the gels from the CT group was not different from that of the HPL group, consistent with the ELISA results. Our results suggest that phosphorylation differences caused the observed migration differences of HSP27. Together with the UMP-CMP kinase abundance differences, we conclude that kinase and/or phosphatase activity are different in CT and HPL. This paper was presented at the 38th Annual Meeting of the German Society for Mass Spectrometry (DGMS) held in March 2005 in Rostock, Germany.