Warping two-dimensional electrophoresis gel images to correct for geometric distortions of the spot pattern

Two-Dimensional Gel Electrophoresis Image Analysis

Gel-based proteomics is still quite widespread due to its high-resolution power; the experimental approach is based on differential analysis, where groups of samples (e.g., control vs diseased) are compared to identify panels of potential biomarkers. However, the reliability of the result of the differential analysis is deeply influenced by 2D-PAGE maps image analysis procedures. The analysis of 2D-PAGE images consists of several steps, such as image preprocessing, spot detection and quantitation, image warping and alignment, spot matching. Several approaches are present in literature, and classical or last-generation commercial software packages exploit different algorithms for each step of the analysis. Here, the most widespread approaches and a comparison of the different strategies are presented.

2D-DIGE and Fluorescence Image Analysis

2D-DIGE is still a very widespread technique in proteomics for the identification of panels of biomarkers, allowing to tackle with some important drawback of classical two-dimensional gel-electrophoresis. However, once 2D-gels are obtained, they must undergo a quite articulated multistep image analysis procedure before the final differential analysis via statistical mono- and multivariate methods. Here, the main steps of image analysis software are described and the most recent procedures reported in the literature are briefly presented.

GENOCOP Algorithm and Hierarchical Grid Transformation for Image Warping of Two-Dimensional Gel Electrophoretic Maps

Hierarchical grid transformation is a powerful hierarchical approach to 2-D map warping, able to model both global and local deformations. The algorithm can be stopped when a desired degree of accuracy in the images alignment is obtained. The deformed image is warped and aligned to the target image using a grid where the number of nodes increases in each step of the algorithm. The numerical optimization of the position of the nodes of the grid can be efficiently solved by genetic algorithms, ensuring the achievement of the optimal position of the nodes with a low computational cost with respect to other methods. Here, the optimization of the position of the nodes is carried out by GENOCOP (genetic algorithm for numerical optimization of constrained problems), refined by the following conjugate gradient optimization step. The modeling of the warped space is then achieved by a spline model where some constraints are introduced in the choice of the nodes that are moved. The whole procedure can be intended as an evolutionary method that models the deformation of the gel map at different levels of detail.

Image Pretreatment Tools II: Normalization Techniques for 2-DE and 2-D DIGE

Gel electrophoresis is usually applied to identify different protein expression profiles in biological samples (e.g., control vs. pathological, control vs. treated). Information about the effect to be investigated (a pathology, a drug, a ripening effect, etc.) is however generally confounded with experimental variability that is quite large in 2-DE and may arise from small variations in the sample preparation, reagents, sample loading, electrophoretic conditions, staining and image acquisition. Obtaining valid quantitative estimates of protein abundances in each map, before the differential analysis, is therefore fundamental to provide robust candidate biomarkers. Normalization procedures are applied to reduce experimental noise and make the images comparable, improving the accuracy of differential analysis. Certainly, they may deeply influence the final results, and to this respect they have to be applied with care. Here, the most widespread normalization procedures are described both for what regards the applications to 2-DE and 2D Difference Gel-electrophoresis (2-D DIGE) maps.

Spot Matching of 2-DE Images Using Distance, Intensity, and Pattern Information

The analysis of a large number of two-dimensional gel electrophoresis (2-DE) images requires developing automatic methods. In such analyses, spot matching plays a fundamental role, in particular for the identification of proteins. We describe a simple and accurate method which allows to automatically and accurately match spots in 2-DE images. The method consists of simultaneously exploiting the distance between the spots, their intensity, and the pattern formed by their spatial configuration.

Algorithms for Warping of 2-D PAGE Maps

Software-based image analysis of 2-D PAGE maps is an important step for the investigation of proteome. Warping algorithms, which are employed to register spots among gels, are able to overcome the difficulties due to the low reproducibility of this analytical technique. Over the years, the research of new matching and warping mathematical methods has allowed the development of several routine applications of easy-to-use software. This chapter describes common and basic spatial transformations used for the alignment of protein spots present in different gel maps; some recently new approaches are also presented.

Effects of increased CO2 on fish gill and plasma proteome

Ocean acidification and warming are both primarily caused by increased levels of atmospheric CO₂, and marine organisms are exposed to these two stressors simultaneously. Although the effects of temperature on fish have been investigated over the last century, the long-term effects of moderate CO₂ exposure and the combination of both stressors are almost entirely unknown. A proteomics approach was used to assess the adverse physiological and biochemical changes that may occur from the exposure to these two environmental stressors. We analysed gills and blood plasma of Atlantic halibut (Hippoglossus hippoglossus) exposed to temperatures of 12°C (control) and 18°C (impaired growth) in combination with control (400 µatm) or high-CO₂ water (1000 µatm) for 14 weeks. The proteomic analysis was performed using two-dimensional gel electrophoresis (2DE) followed by Nanoflow LC-MS/MS using a LTQ-Orbitrap. The high-CO₂ treatment induced the up-regulation of immune system-related proteins, as indicated by the up-regulation of the plasma proteins complement component C3 and fibrinogen β chain precursor in both temperature treatments. Changes in gill proteome in the high-CO₂ (18°C) group were mostly related to increased energy metabolism proteins (ATP synthase, malate dehydrogenase, malate dehydrogenase thermostable, and fructose-1,6-bisphosphate aldolase), possibly coupled to a higher energy demand. Gills from fish exposed to high-CO₂ at both temperature treatments showed changes in proteins associated with increased cellular turnover and apoptosis signalling (annexin 5, eukaryotic translation elongation factor 1γ, receptor for protein kinase C, and putative ribosomal protein S27). This study indicates that moderate CO₂-driven acidification, alone and combined with high temperature, can elicit biochemical changes that may affect fish health.

Effects of post-electrophoretic analysis on variance in gel-based proteomics

2D electrophoresis (2DE) is a prominent separation method for complex proteomes. Although recent advances have increased the utility of this method in quantitative proteomics studies, many sources of variance still exist. This review discusses the post-electrophoretic sources of variance in current 2DE analysis. The essential improvements in protein visualization and software algorithms that have made 2DE a leading quantitative proteomics method are briefly reviewed. A number of shortcomings in the post-electrophoretic analysis of 2DE data that require further attention are highlighted. Topics discussed include protein visualization and image acquisition, internal standards and normalization methods, background subtraction algorithms, normality of distribution, and the need for standardized tests for the evaluation of 2DE analysis software packages.

GENOCOP algorithm and hierarchical grid transformation for image warping of two dimensional gel eletrophoretic maps

Hierarchical grid transformation is a powerful approach to SDS 2DPAGE maps warping. The hierarchy of the warping transformation is able to model both global and local deformations of the gels and the algorithm can be stopped when a certain degree of accuracy in the image alignment is obtained. The numerical optimization of the position of the nodes of the grid that are responsible for the image warping is a multivariate task that can be solved efficiently using Genetic Algorithms. The use of Genetic Algorithms ensures that an optimal position of the nodes can be defined with a low computational cost with respect to other methods. The optimal positions of the nodes of the grid can be successfully used for defining a good warping of the gels.

Image analysis tools and emerging algorithms for expression proteomics

Since their origins in academic endeavours in the 1970s, computational analysis tools have matured into a number of established commercial packages that underpin research in expression proteomics. In this paper we describe the image analysis pipeline for the established 2-DE technique of protein separation, and by first covering signal analysis for MS, we also explain the current image analysis workflow for the emerging high-throughput 'shotgun' proteomics platform of LC coupled to MS (LC/MS). The bioinformatics challenges for both methods are illustrated and compared, whereas existing commercial and academic packages and their workflows are described from both a user's and a technical perspective. Attention is given to the importance of sound statistical treatment of the resultant quantifications in the search for differential expression. Despite wide availability of proteomics software, a number of challenges have yet to be overcome regarding algorithm accuracy, objectivity and automation, generally due to deterministic spot-centric approaches that discard information early in the pipeline, propagating errors. We review recent advances in signal and image analysis algorithms in 2-DE, MS, LC/MS and Imaging MS. Particular attention is given to wavelet techniques, automated image-based alignment and differential analysis in 2-DE, Bayesian peak mixture models, and functional mixed modelling in MS, and group-wise consensus alignment methods for LC/MS.

Computational methods for analysis of two-dimensional gels

Two-dimensional gel electrophoresis (2D gels) is an essential quantitative proteomics technique that is frequently used to study differences between samples of clinical relevance. Although considered to have a low throughput, 2D gels can separate thousands of proteins in one gel, making it a good complementary method to MS-based protein quantification. The main drawback of the technique is the tendency of large and hydrophobic proteins such as membrane proteins to precipitate in the isoelectric focusing step. Furthermore, tests using different programs with distinct algorithms for 2D-gel analysis have shown inconsistent ratio values. The aim here is therefore to provide a discussion of algorithms described for the analysis of 2D gels.

Identifying calpain substrates in intact S2 cells of Drosophila

Calpains are cysteine proteases involved in a number of physiological and pathological processes, yet our knowledge of substrates cleaved in vivo, in intact cells, is scarce. In this work we made an attempt to develop a technique for finding calpain substrates in intact Drosophila Schneider S2 cells. The procedure consists in comparative 2D gelelectrophoresis: three identical samples were treated in different ways: A (control, no addition), B, activated (Ca(2+) and ionomycin added), C, inactivated (additions as in B+specific calpain inhibitor). 2D gel pattern were analyzed by densitometry. Spots showing density relation A>B<<C were identified by mass spectroscopy. In a typical run, 11 candidate substrates were recognized; out of these, four were randomly selected: all four were verified to be calpain substrates, by digestion of the recombinant protein with recombinant calpain.

A new method for 2D gel spot alignment: application to the analysis of large sample sets in clinical proteomics

BACKGROUND

In current comparative proteomics studies, the large number of images generated by 2D gels is currently compared using spot matching algorithms. Unfortunately, differences in gel migration and sample variability make efficient spot alignment very difficult to obtain, and, as consequence most of the software alignments return noisy gel matching which needs to be manually adjusted by the user.

RESULTS

We present Sili2DGel an algorithm for automatic spot alignment that uses data from recursive gel matching and returns meaningful Spot Alignment Positions (SAP) for a given set of gels. In the algorithm, the data are represented by a graph and SAP by specific subgraphs. The results are returned under various forms (clickable synthetic gel, text file, etc.). We have applied Sili2DGel to study the variability of the urinary proteome from 20 healthy subjects.

CONCLUSION

Sili2DGel performs noiseless automatic spot alignment for variability studies (as well as classical differential expression studies) of biological samples. It is very useful for typical clinical proteomic studies with large number of experiments.

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.

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.

Automated image alignment for 2D gel electrophoresis in a high-throughput proteomics pipeline

MOTIVATION

The quest for high-throughput proteomics has revealed a number of challenges in recent years. Whilst substantial improvements in automated protein separation with liquid chromatography and mass spectrometry (LC/MS), aka 'shotgun' proteomics, have been achieved, large-scale open initiatives such as the Human Proteome Organization (HUPO) Brain Proteome Project have shown that maximal proteome coverage is only possible when LC/MS is complemented by 2D gel electrophoresis (2-DE) studies. Moreover, both separation methods require automated alignment and differential analysis to relieve the bioinformatics bottleneck and so make high-throughput protein biomarker discovery a reality. The purpose of this article is to describe a fully automatic image alignment framework for the integration of 2-DE into a high-throughput differential expression proteomics pipeline.

RESULTS

The proposed method is based on robust automated image normalization (RAIN) to circumvent the drawbacks of traditional approaches. These use symbolic representation at the very early stages of the analysis, which introduces persistent errors due to inaccuracies in modelling and alignment. In RAIN, a third-order volume-invariant B-spline model is incorporated into a multi-resolution schema to correct for geometric and expression inhomogeneity at multiple scales. The normalized images can then be compared directly in the image domain for quantitative differential analysis. Through evaluation against an existing state-of-the-art method on real and synthetically warped 2D gels, the proposed analysis framework demonstrates substantial improvements in matching accuracy and differential sensitivity. High-throughput analysis is established through an accelerated GPGPU (general purpose computation on graphics cards) implementation.

AVAILABILITY

Supplementary material, software and images used in the validation are available at http://www.proteomegrid.org/rain/.

Comparative proteomics of industrial lager yeast reveals differential expression of the cerevisiae and non-cerevisiae parts of their genomes

The proteomes of three industrial lager beer strains, CMBS33, OG2252 and A15, were analysed under standardised laboratory growth conditions. Protein spots in the 2-DE pattern of the lager strains were subjected to MS/MS to identify protein variants. We found the protein composition of the three lager strains to be qualitatively rather similar, while being substantially different from the Saccharomyces cerevisiae strain BY4742. Database searches using several fully sequenced genomes from the Saccharomyces genera indicated that the non-cerevisiae proteins in the 2-D pattern of lager strains were most closely related to S. bayanus. For many proteins the regulation of the bayanus-like protein and its cerevisiae counterpart varied in a strain-dependent manner, e.g. the bayanus-like form of Tdh3p was roughly eight-fold more abundant than the cerevisiae form in the OG2252 strain. We also found differential regulation of cerevisiae- and bayanus-like proteins during various stress conditions like low temperature growth, and adaptation to high temperatures or high salinity, e.g. for Arg1p, Sti1p and Pdc1p. Our data on the differential regulation of the two genomes in these hybrid strains may have important industrial implications for strain improvement and strain protection.

The state of the art in the analysis of two-dimensional gel electrophoresis images

Software-based image analysis is a crucial step in the biological interpretation of two-dimensional gel electrophoresis experiments. Recent significant advances in image processing methods combined with powerful computing hardware have enabled the routine analysis of large experiments. We cover the process starting with the imaging of 2-D gels, quantitation of spots, creation of expression profiles to statistical expression analysis followed by the presentation of results. Challenges for analysis software as well as good practices are highlighted. We emphasize image warping and related methods that are able to overcome the difficulties that are due to varying migration positions of spots between gels. Spot detection, quantitation, normalization, and the creation of expression profiles are described in detail. The recent development of consensus spot patterns and complete expression profiles enables one to take full advantage of statistical methods for expression analysis that are well established for the analysis of DNA microarray experiments. We close with an overview of visualization and presentation methods (proteome maps) and current challenges in the field.

Statistics for proteomics: Experimental design and 2-DE differential analysis

Proteomics relies on the separation of complex protein mixtures using bidimensional electrophoresis. This approach is largely used to detect the expression variations of proteins prepared from two or more samples. Recently, attention was drawn on the reliability of the results published in literature. Among the critical points identified were experimental design, differential analysis and the problem of missing data, all problems where statistics can be of help. Using examples and terms understandable by biologists, we describe how a collaboration between biologists and statisticians can improve reliability of results and confidence in conclusions.

Aligning families of two-dimensional gels by a combined multiresolution forward-inverse transformation approach

A new method for aligning families of two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) images arising in proteomics studies is presented. Forward piecewise bilinear transformations are used to determine an ideal gel and to obtain an initial alignment of the family of gels to this ideal gel. Both the ideal landmarks and the coefficients defining the transformations are obtained by solving a quadratic programming problem. The alignment is then improved by using inverse transformations on finer grids. Numerical results for a family of 123 gels are reported.

Examination of 2-DE in the Human Proteome Organisation Brain Proteome Project pilot studies with the new RAIN gel matching technique

The Human Proteome Organisation (HUPO) Brain Proteome Project (BPP) pilot studies have generated over 200 2-D gels from eight participating laboratories. This data includes 67 single-channel and 60 DIGE gels comparing 30 whole frozen C57/BL6 female mouse brains, ten each at embryonic day 16, postnatal day 7 (juvenile) and postnatal day 54-56 (adult); and ten single-channel and three DIGE gels comparing human epilepsy surgery of the temporal front lobe with a corresponding post-mortem specimen. The samples were generated centrally and distributed to the participating laboratories, but otherwise no restrictions were placed on sample preparation, running and staining protocols, nor on the 2-D gel analysis packages used. Spots were characterised by MS and the annotated gel images published on a ProteinScape web server. In order to examine the resultant differential expression and protein identifications, we have reprocessed a large subset of the gels using the newly developed RAIN (Robust Automated Image Normalisation) 2-D gel matching algorithm. Traditional approaches use symbolic representation of spots at the very early stages of the analysis, which introduces persistent errors due to inaccuracies in spot modelling and matching. With RAIN, image intensity distributions, rather than selected features, are used, where smooth geometric deformation and expression bias are modelled using multi-resolution image registration and bias-field correction. The method includes a new approach of volume-invariant warping which ensures the volume of protein expression under transformation is preserved. An image-based statistical expression analysis phase is then proposed, where small insignificant expression changes over one gel pair can be revealed when reinforced by the same consistent changes in others. Results of the proposed method as applied to the HUPO BPP data show significant intra-laboratory improvements in matching accuracy over a previous state-of-the-art technique, Multi-resolution Image Registration (MIR), and the commercial Progenesis PG240 package.

Protein Image Alignment via Piecewise Affine Transformations

We present a new approach for aligning families of 2D gels. Instead of choosing one of the gels as reference and performing a pairwise alignment, we construct an ideal gel that is representative of the entire family and obtain a set of piecewise affine transformations that optimally align each gel of the family to the ideal gel. The coefficients defining the transformations as well as the ideal landmarks are obtained as the solution of a large-scale quadratic programming problem that can be solved efficiently by interior-point methods.

Analysis of electrophoretic patterns of arbitrarily primed PCR profiling

We present a mathematical algorithm for the analysis of electrophoretic patterns resulting from arbitrarily primed PCR profiling. The algorithm is based on the established mathematical procedures applied to the analysis of digital images of gel patterns. The algorithm includes (a) transformation of the image into a matrix form, (b) identification of every electrophoretic lane as a set of matrix columns that are further mathematically processed, (c) averaging of matrix columns corresponding to electrophoretic lanes that define lane representatives, (d) elimination of "smiling" bands, (e) solving the problem of a lane offset, and (f) removal of the background. Representation of individual electrophoretic lanes in the form of functions allows interlane comparisons and further mathematical analysis. Direct comparison of selected lanes was obtained by employing correlation analysis. Gel images were those obtained after arbitrarily primed PCR analysis of DNA that underwent damage induced by gamma radiation from a (60)Co source. The applied method proved to be useful for elimination of subjectivity of visual inspection. It offers the possibility to avoid overlooking important differences in case of suboptimal electrophoretic resolution. In addition, higher precision is achieved in the assessment of quantitative differences due to better insight into experimental artifacts. These simple mathematical methods offer an open-type algorithm, i.e., this algorithm enables easy implementation of different parameters that may be useful for other analytical needs.

Analysing proteomic data

The rapid growth of proteomics has been made possible by the development of reproducible 2D gels and biological mass spectrometry. However, despite technical improvements 2D gels are still less than perfectly reproducible and gels have to be aligned so spots for identical proteins appear in the same place. Gels can be warped by a variety of techniques to make them concordant. When gels are manipulated to improve registration, information is lost, so direct methods for gel registration which make use of all available data for spot matching are preferable to indirect ones. In order to identify proteins from gel spots a property or combination of properties that are unique to that protein are required. These can then be used to search databases for possible matches. Molecular mass, pI, amino acid composition and short sequence tags can all be used in database searches. Currently the method of choice for protein identification is mass spectrometry. Proteins are eluted from the gels and cleaved with specific endoproteases to produce a series of peptides of different molecular mass. In peptide mass fingerprinting, the peptide profile of the unknown protein is compared with theoretical peptide libraries generated from sequences in the different databases. Tandem mass spectroscopy (MS/MS) generates short amino acid sequence tags for the individual peptides. These partial sequences combined with the original peptide masses are then used for database searching, greatly improving specificity. Increasingly protein identification from MS/MS data is being fully or partially automated. When working with organisms, which do not have sequenced genomes (the case with most helminths), protein identification by database searching becomes problematical. A number of approaches to cross species protein identification have been suggested, but if the organism being studied is only distantly related to any organism with a sequenced genome then the likelihood of protein identification remains small. The dynamic nature of the proteome means that there really is no such thing as a single representative proteome and a complete set of metadata (data about the data) is going to be required if the full potential of database mining is to be realised in the future.

Elastic registration of biological images using vector-spline regularization

We present an elastic registration algorithm for the alignment of biological images. Our method combines and extends some of the best techniques available in the context of medical imaging. We express the deformation field as a B-spline model, which allows us to deal with a rich variety of deformations. We solve the registration problem by minimizing a pixelwise mean-square distance measure between the target image and the warped source. The problem is further constrained by way of a vector-spline regularization which provides some control over two independent quantities that are intrinsic to the deformation: its divergence, and its curl. Our algorithm is also able to handle soft landmark constraints, which is particularly useful when parts of the images contain very little information or when its repartition is uneven. We provide an optimal analytical solution in the case when only landmarks and smoothness considerations are taken into account. We have applied our approach to perform the elastic registration of images such as electrophoretic gels and fly embryos. The validation of the results by experts has been favorable in all cases.

ProteomeGRID: towards a high-throughput proteomics pipeline through opportunistic cluster image computing for two-dimensional gel electrophoresis

The quest for high-throughput proteomics has revealed a number of critical issues. Whilst improved two-dimensional gel electrophoresis (2-DE) sample preparation, staining and imaging issues are being actively pursued by industry, reliable high-throughput spot matching and quantification remains a significant bottleneck in the bioinformatics pipeline, thus restricting the flow of data to mass spectrometry through robotic spot excision and protein digestion. To this end, it is important to establish a full multi-site Grid infrastructure for the processing, archival, standardisation and retrieval of proteomic data and metadata. Particular emphasis needs to be placed on large-scale image mining and statistical cross-validation for reliable, fully automated differential expression analysis, and the development of a statistical 2-DE object model and ontology that underpins the emerging HUPO PSI GPS (Human Proteome Organization Proteomics Standards Initiative General Proteomics Standards). The first step towards this goal is to overcome the computational and communications burden entailed by the image analysis of 2-DE gels with Grid enabled cluster computing. This paper presents the proTurbo framework as part of the ProteomeGRID, which utilises Condor cluster management combined with CORBA communications and JPEG-LS lossless image compression for task farming. A novel probabilistic eager scheduler has been developed to minimise make-span, where tasks are duplicated in response to the likelihood of the Condor machines' owners evicting them. A 60 gel experiment was pair-wise image registered (3540 tasks) on a 40 machine Linux cluster. Real-world performance and network overhead was gauged, and Poisson distributed worker evictions were simulated. Our results show a 4:1 lossless and 9:1 near lossless image compression ratio and so network overhead did not affect other users. With 40 workers a 32x speed-up was seen (80% resource efficiency), and the eager scheduler reduced the impact of evictions by 58%.

Multivariate approaches in plant science

The objective of proteomics is to get an overview of the proteins expressed at a given point in time in a given tissue and to identify the connection to the biochemical status of that tissue. Therefore sample throughput and analysis time are important issues in proteomics. The concept of proteomics is to encircle the identity of proteins of interest. However, the overall relation between proteins must also be explained. Classical proteomics consist of separation and characterization, based on two-dimensional electrophoresis, trypsin digestion, mass spectrometry and database searching. Characterization includes labor intensive work in order to manage, handle and analyze data. The field of classical proteomics should therefore be extended to also include handling of large datasets in an objective way. The separation obtained by two-dimensional electrophoresis and mass spectrometry gives rise to huge amount of data. We present a multivariate approach to the handling of data in proteomics with the advantage that protein patterns can be spotted at an early stage and consequently the proteins selected for sequencing can be selected intelligently. These methods can also be applied to other data generating protein analysis methods like mass spectrometry and near infrared spectroscopy and examples of application to these techniques are also presented. Multivariate data analysis can unravel complicated data structures and may thereby relieve the characterization phase in classical proteomics. Traditionally statistical methods are not suitable for analysis of the huge amounts of data, where the number of variables exceed the number of objects. Multivariate data analysis, on the other hand, may uncover the hidden structures present in these data. This study takes its starting point in the field of classical proteomics and shows how multivariate data analysis can lead to faster ways of finding interesting proteins. Multivariate analysis has shown interesting results as a supplement to classical proteomics and added a new dimension to the field of proteomics.