Analysis of scaling methods to minimize experimental variations in two-dimensional electrophoresis quantitative data: application to the comparison of maize inbred lines

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.

2D-DIGE: comparative proteomics of cellular signalling pathways

Two-dimensional (2-D) gel electrophoresis concerted with protein identification by mass spectrometry (MS) is an extremely powerful method for comparative expression profiling of complex protein samples such as cell lysates. The highly resolutive 2-D electrophoresis allows the separation of heterogeneous protein samples on the basis of isoelectric point (pI), molecular mass (Mr), solubility, and relative abundance ((1) J Biol Chem 250: 4007-4021, 1975; (2) Electrophoresis 14: 1067-1073, 1993). Consequently, it provides a comprehensive view of a proteome state ((3) Electrophoresis 21: 1037-1053, 2000), where variations in protein expression levels, isoforms, or post-translational modifications (e.g. phosphorylation) can be highlighted and investigated ((4) Electrophoresis 21: 2196-2208, 2000). Furthermore, this allows the identification of biological markers that characterize a specific physiological or pathological background of a cell or a tissue ((5) Proteomics 1: 397-408, 2001; (6) J Bacteriol 179: 7595-7599, 1997). In this way one can compare the effects of a stimulus or drug on cells or tissue, or more importantly, analyse the effects of disease on the expression level of proteins. Relatively recently, conventional 2-D gel electrophoresis has been combined with protein labelling strategies using up to three different fluorescent dyes to allow comparative analysis of different protein samples within a single 2-D gel platform. In this technique, termed differential in-gel electrophoresis (DIGE), samples are labelled separately then combined and run on the same 2D gel minimizing experimental variation and greatly facilitating spot matching. When three CyDyes (Cy2, Cy3, and Cy5) have been used, three images of the gel are captured then superposed to localize the differentially regulated spots on the 2-D gel using image analysis software. This is an extremely powerful tool in comparative proteomics as these dyes provide a linear response to protein concentration up to five orders of magnitude and great sensitivity with detection down to 125 pg of a single protein, which is less than needed for MS identification. In this chapter, we describe the basic methods for protein labelling, optimization of the isoelectrofocusing parameters for the first dimension (where proteins are separated according to their isoelectric point (pI)), sodium-dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) separation for the second dimension (based on molecular weight (MW)), and different post-staining protocols of the 2-D gel and protein preparation for mass spectrometry identification.

The holm oak leaf proteome: analytical and biological variability in the protein expression level assessed by 2-DE and protein identification tandem mass spectrometry de novo sequencing and sequence similarity searching

As a first approach in establishing the holm oak leaf proteome, we have optimised a protocol for this plant and tissue which includes the following steps: trichloroacetic acid-acetone extraction, two-dimensional gel electrophoresis (2-DE) on pH 5 to 8 linear gradient immobilised pH gradient strips as the first dimension, and sodium dodecyl sulfate-polyacrylamide gel electrophoresis on 13% polyacrylamide gels as the second one. Proteins were detected by Coomassie staining. Gel images were recorded and digitalized, and the protein spots quantified by using a linear regression equation of protein quantity on spot volume obtained against standard proteins. Analytical variance was calculated for one-hundred protein spots from three replicate 2-DE gels of the same protein extract. Biological variance was determined for the same protein spots from independent tissue extracts corresponding to leaves from different trees, or the same tree at different orientations or sampling times during a day. Values of 26% for the analytical variance and 58.6% for the biological variance among independent trees were obtained. These values provide a quantified and statistical basis for the evaluation of protein expression changes in comparative proteomic investigations with this species. A representative set of the major proteins, covering the isoelectric point range of 5 to 8 and the relative molecular mass(r) range of 14 to 78 kDa, were subjected to liquid chromatography-tandem mass spectrometry analysis. Due to the absence of Quercus DNA or protein sequence databases, a method based on the procedure reported by Liska and Shevchenko including de novo sequencing and BLAST similarity searching against other plant species databases was used for protein identification. Out of 43 analysed spots, 35 were positively identified. The identified proteins mainly corresponded to enzymes involved in photosynthesis and energetic metabolism, with a significant number corresponding to RubisCO.

Maximising sensitivity for detecting changes in protein expression: Experimental design using minimal CyDyes

DIGE is a powerful tool for measuring changes in protein expression between samples. Here we assess the assumptions of normality and heterogeneity of variance that underlie the univariate statistical tests routinely used to detect proteins with expression changes. Furthermore, the technical variance experienced in a multigel experiment is assessed here and found to be reproducible within- and across-sample types. Utilising the technical variance measured, a power study is completed for several "typical" fold changes in expression commonly used as thresholds by researchers. Based on this study using DeCyder, guidance is given on the number of gel replicates that are needed for the experiment to have sufficient sensitivity to detect expression changes. A two-dye system based on utilising just Cy3 and Cy5 was found to be more reproducible than the three-dye system. A power and cost-benefit analysis performed here suggests that the traditional three-dye system would use fewer resources in studies where multiple samples are compared. Technical variance was shown to encompass both experimental and analytical noise and thus is dependent on the analytical software utilised. Data is provided as a resource to the community to assess alternative software and upgrades.

Statistical exploration of variation in quantitative two-dimensional gel electrophoresis data

Two-dimensional gel electrophoresis is a major technique in global analysis at the protein level. This paper presents an examination of spot volume data from three gel sets with radioactively labeled yeast Saccharomyces cerevisiae proteins. A strong variance versus mean dependence in data was found to be stabilized by applying a shifted logarithmic transformation. However, transformed data showed a remaining substantial variance heterogeneity for different proteins. Furthermore, examination of studentized residuals revealed that transformed data were approximately normally distributed and that there were spatial correlations among the measurement errors in the gel.

Water deficits affect caffeate O-methyltransferase, lignification, and related enzymes in maize leaves. A proteomic investigation

Drought is a major abiotic stress affecting all levels of plant organization and, in particular, leaf elongation. Several experiments were designed to study the effect of water deficits on maize (Zea mays) leaves at the protein level by taking into account the reduction of leaf elongation. Proteomic analyses of growing maize leaves allowed us to show that two isoforms of caffeic acid/5-hydroxyferulic 3-O-methyltransferase (COMT) accumulated mostly at 10 to 20 cm from the leaf point of insertion and that drought resulted in a shift of this region of maximal accumulation toward basal regions. We showed that this shift was due to the combined effect of reductions in growth and in total amounts of COMT. Several other enzymes involved in lignin and/or flavonoid synthesis (caffeoyl-CoA 3-O-methyltransferase, phenylalanine ammonia lyase, methylenetetrahydrofolate reductase, and several isoforms of S-adenosyl-l-methionine synthase and methionine synthase) were highly correlated with COMT, reinforcing the hypothesis that the zone of maximal accumulation corresponds to a zone of lignification. According to the accumulation profiles of the enzymes, lignification increases in leaves of control plants when their growth decreases before reaching their final size. Lignin levels analyzed by thioacidolysis confirmed that lignin is synthesized in the region where we observed the maximal accumulation of these enzymes. Consistent with the levels of these enzymes, we found that the lignin level was lower in leaves of plants subjected to water deficit than in those of well-watered plants.

Deciphering genetic variations of proteome responses to water deficit in maize leaves

The proteome of the basal part of growing Zea mays leaves was analyzed from 4 to 14 d after stopping watering and in well watered controls. The relative quantity of 46 proteins was found to increase in leaves of plants submitted to water deficit. Different types of responses were observed, some proteins showing a constant increase during water deficit, while others showed stabilization after a first increase or a transient increase. Isoforms encoded by the same gene showed different responses. The response to water deficit showed genetic variation. Some increased proteins were induced specifically in one of the two studied genotypes (e.g. ASR1) while others were significantly induced in both genotypes but to a different level or with different kinetics. Analyses of relations between protein quantities, relative water content (RWC) and abscisic acid (ABA) concentration allowed us to show that the quantitative variation of some proteins (e.g. ABA45 and OSR40 proteins) was linked to differences in ABA accumulation between the genotypes. Other proteins showed genetic variations that were not related to differences in water status or ABA concentration (e.g. a cystatin). Data obtained from these experiments, together with data from other experiments, contribute to the characterization of maize proteome response to drought in different conditions and in different genotypes. This characterization allows the search for candidate proteins, i.e. for protein whose genetic variation of expression could be partly responsible for the variability of plant responses to drought.

The UV-B stimulon of the terrestrial cyanobacterium Nostoc commune comprises early shock proteins and late acclimation proteins

The UV-B and desiccation-tolerant terrestrial cyanobacterium Nostoc commune was grown under defined UV irradiation. Proteome changes were monitored in the membrane and the cytosolic and the extracellular fractions. Tools were developed to separate stress-triggered from growth stage-dependent changes. UV-B changed the relative cellular concentration of 493 out of 1,350 protein spots at least by a factor of three, rendering the UV-B stimulon of N. commune the most complex one described so far. It comprises two different parts: an early shock response influencing 214 proteins and a late acclimation response involving 279 proteins. The shock response comprised many membrane or membrane-associated proteins, whereas the acclimation response mainly changed cytosolic proteins. Most of the shock-induced changes were transient and did not overlap with the acclimation response. In the extracellular fraction, UV irradiation induced superoxide dismutase and the water stress protein. In total, 27 intracellular, UV-B-induced proteins were partially sequenced by electrospray ionization tandem mass spectrometry. Three functional classes were identified: proteins involved in lipid metabolism, in carbohydrate metabolism and in regulatory pathways. About 50% of the sequenced proteins were homologous to cyanobacterial database entries with un-known function. Interestingly, all of these proteins belong to the UV-B acclimation response. We conclude that the UV-B shock response and the UV-B acclimation response represent two completely different and remarkably complex strategies of N. commune to protect itself against UV-B radiation in its natural environment.

Compression wood-responsive proteins in developing xylem of maritime pine (Pinus pinaster ait.)

When a conifer shoot is displaced from its vertical position, compression wood (CW) is formed on the under side and can eventually return the shoot to its original position. Changes in cell wall structure and chemistry associated with CW are likely to result from differential gene/protein expression. Two-dimensional polyacrylamide gel electrophoresis of differentiating xylem proteins was combined with the physical characterization of wooden samples to identify and characterize CW-responsive proteins. Differentiating xylem was harvested from a 22-year-old crooked maritime pine (Pinus pinaster Ait.) tree. Protein extracted from different samples were revealed by high-resolution silver stained two-dimensional polyacrylamide gel electrophoresis and analyzed with a computer-assisted system for single spot quantification. Growth strain (GS) measurements allowed xylem samples to be classified quantitatively from normal wood to CW. Regression of lignin and cellulose content on GS showed that an increase in the percentage of lignin and a decrease of the percentage of cellulose corresponded to increasing GS values, i.e. CW. Of the 137 studied spots, 19% were significantly associated with GS effect. Up-regulated proteins included 1-aminocyclopropane-1-carboxylate oxidase (an ethylene forming enzyme), a putative transcription factor, two lignification genes (caffeic O-methyltransferase and caffeoyl CoA-O-methyltransferase), members of the S-adenosyl-L-methionine-synthase gene family, and enzymes involved in nitrogen and carbon assimilation (glutamine synthetase and fructokinase). A clustered correlation analysis was performed to study simultaneously protein expression along a gradient of gravistimulated stressed xylem tissue. Proteins were found to form "expression clusters" that could identify: (a) Gene product under similar control mechanisms, (b) partner proteins, or (c) functional groups corresponding to specialized pathways. The possibility of obtaining regulatory correlations and anticorrelations between proteins provide us with a new category of homology (regulatory homology) in tracing functional relationships.

Effects of renovascular hypertension on myocardial protein patterns: analysis by computer-assisted two-dimensional gel electrophoresis

Hypertensive heart disease caused by renovascular hypertension reflects the response of the heart to an increased afterload and neurohormonal stimulation. We hypothesized that in this condition the composition of the myocardial proteins of rats was altered. To identify yet unknown quantitative and qualitative differences in myocardial proteins in renovascular hypertensive heart disease, we analyzed protein patterns by computer-assisted two-dimensional polyacrylamide large gel electrophoresis. Renovascular hypertension was induced by placing a silver clip on the left renal artery in 9-week-old rat siblings. Sham-operated animals served as controls. Systolic blood pressure (197 +/- 19 mm Hg) and heart/body weight ratios (0.36 +/- 0.04%) were significantly increased in the hypertensive animals. Twenty protein patterns from the left ventricle of five hypertensive and five control rats were compared. The molecular mass and isoelectric point (pI) of proteins spots ranging from 13 to 100 kDa and from 4.5 to 8.5, respectively, were determined using marker proteins. In total, 761 +/- 88 protein spots were resolved in all twenty gels. For the quantitative data analysis a univariate (Mann-Whitney test) as well as a multivariate statistical approach (correspondence analysis) were applied. Only one myocardial protein spot (molecular mass = 41.3 kDa; pI = 6.3) was decreased by more than twofold (p < 0.05) in renovascular hypertension. The vast majority of spots did not indicate a significant alteration of intensity. Left ventricular hypertrophy in early renovascular hypertension induces a form of myocardial hypertrophy that conserves the naturally occurring protein expression pattern.

Protein changes in response to progressive water deficit in maize . Quantitative variation and polypeptide identification

Three-week-old plants of two unrelated lines of maize (Zea mays L.) and their hybrid were submitted to progressive water stress for 10 d. Changes induced in leaf proteins were studied by two-dimensional electrophoresis and quantitatively analyzed using image analysis. Seventy-eight proteins out of a total of 413 showed a significant quantitative variation (increase or decrease), with 38 of them exhibiting a different expression in the two genotypes. Eleven proteins that increased by a factor of 1.3 to 5 in stressed plants and 8 proteins detected only in stressed plants were selected for internal amino acid microsequencing, and by similarity search 16 were found to be closely related to previously reported proteins. In addition to proteins already known to be involved in the response to water stress (e.g. RAB17 [Responsive to ABA]), several enzymes involved in basic metabolic cellular pathways such as glycolysis and the Krebs cycle (e.g. enolase and triose phosphate isomerase) were identified, as well as several others, including caffeate O-methyltransferase, the induction of which could be related to lignification.

Identification of procaryotic developmental stages by statistical analyses of two-dimensional gel patterns

Multivariate statistical comparisons of two-dimensional protein (2-D) gel patterns were used for the first time to define stages of a biological developmental system. The differentiating procaryote, Streptomyces coelicolor, was radiolabeled in liquid cultures at 16 intervals during development, and radioactive proteins were separated and quantified on 2-D gels. Cluster, principal component, and correlation analyses classified these gel patterns into four distinct groups, each reflecting a pattern of gene expression specific for a stage of development. These studies focused our attention on a phase of arrested growth as a key regulatory transition leading to secondary metabolism and a phase of renewed growth. Proteins whose synthesis was switched on or off during the "transitional" phase (some 21 and 18, respectively) were identified and will be the focus of future studies designed to identify their physiological or regulatory function.

Proteome research: complementarity and limitations with respect to the RNA and DNA worlds

A methodological overview of proteome analysis is provided along with details of efforts to achieve high-throughput screening (HTS) of protein samples derived from two-dimensional electrophoresis gels. For both previously sequenced organisms and those lacking significant DNA sequence information, mass spectrometry has a key role to play in achieving HTS. Prototype robotics designed to conduct appropriate chemistries and deliver 700-1000 protein (genes) per day to batteries of mass spectrometers or liquid chromatography (LC)-based analyses are well advanced, as are efforts to produce high density gridded arrays containing > 1000 proteins on a single matrix assisted laser desorption ionisation/time-of-flight (MALDI-TOF) sample stage. High sensitivity HTS of proteins is proposed by employing principally mass spectrometry in an hierarchical manner: (i) MALDI-TOF-mass spectrometry (MS) on at least 1000 proteins per day; (ii) electrospray ionisation (ESI)/MS/MS for analysis of peptides with respect to predicted fragmentation patterns or by sequence tagging; and (iii) ESI/MS/MS for peptide sequencing. Genomic sequences when complemented with information derived from hybridisation assays and proteome analysis may herald in a new era of holistic cellular biology. The current preoccupation with the absolute quantity of gene-product (RNA and/or protein) should move backstage with respect to more molecularly relevant parameters, such as: molecular half-life; synthesis rate; functional competence (presence or absence of mutations); reaction kinetics; the influence of individual gene-products on biochemical flux; the influence of the environment, cell-cycle, stress and disease on gene-products; and the collective roles of multigenic and epigenetic phenomena governing cellular processes. Proteome analysis is demonstrated as being capable of proceeding independently of DNA sequence information and aiding in genomic annotation. Its ability to confirm the existence of gene-products predicted from DNA sequence is a major contribution to genomic science. The workings of software engines necessary to achieve large-scale proteome analysis are outlined, along with trends towards miniaturisation, analyte concentration and protein detection independent of staining technologies. A challenge for proteome analysis into the future will be to reduce its dependence on two-dimensional (2-D) gel electrophoresis as the preferred method of separating complex mixtures of cellular proteins. Nonetheless, proteome analysis already represents a means of efficiently complementing differential display, high density expression arrays, expressed sequence tags, direct or subtractive hybridisation, chromosomal linkage studies and nucleic acid sequencing as a problem solving tool in molecular biology.

Potential of two-dimensional electrophoresis in routine identification of closely related durum wheat lines

Four closely related durum wheat varieties were compared by computer-assisted analysis of two-dimensional electrophoretic maps of leaf proteins. A low inter-varietal polymorphism was revealed and seven reliable qualitatively varying proteins allowed rapid visual identification of genotypes. For numerous spots, presence/absence or quantitative variations were greatly affected by a batch effect. Several criteria that should be used to discard unreliable spots or gels a priori were reviewed. Nevertheless, it was shown that, provided that the experimental design allows the integration of the batch effect, screening for discriminant markers as well as computing distances based on protein quantity variations are possible and allow variety identification. Euclidean and Mahalanobis distances allowed variety discrimination and single gel classification with a minimum risk of error, not only by taking into account the quantitative variations in discriminant proteins selected by analysis of variance, but also by taking into account all reproducible spots. The possible applications of two-dimensional electrophoresis in variety identification are discussed.

Amino acid analysis of proteins separated by two-dimensional electrophoresis in maize: isoform detection and function identification

The possibility of using experimentally determined amino acid composition to assess relatedness between 75 proteins separated by two-dimensional electrophoresis (2-DE) and to identify them was tested on maize. Two independent parameters, the relative Euclidean distance and the correlation coefficient between the amino acid compositions, were evaluated and used. Previous sequence information for 31 out of the 75 proteins made it possible to evaluate the method for the detection of isoforms and for identification. However, the extension of the interrogation beyond maize to all plant sequences raised the problem of false positives that could nevertheless be limited by replications and by using additional information. The efficiency of the method to assess relatedness between proteins should make amino acid composition analysis a valuable tool in large protein characterization programs based on 2-DE, by facilitating the transfer of information from one well-documented organ/tissue or genotype to another.

Distance indices in a comparison between the A, D, I and R genomes of the Triticeae tribe

Two-dimensional gel electrophoresis of shoot proteins was used to study the relationships between Triticum, Secale, and Hordeum. A high level of polymorphism was found among the 1275 spots scored: only 198 spots were found common to all. But, under the hypothesis that only allelic variations were observed, the mean number of alleles per locus was only two. Phenograms were built from different distance indices. All of them showed Triticum genomes A and D close to each other, Hordeum far from the Triticum cluster, and Secale at an intermediate position. A discussion on the use of various distance indices is presented.

Protein mutations revealed by two-dimensional electrophoresis

High-resolution two-dimensional electrophoresis (2DE) can resolve many hundreds of proteins present in complex mixtures depending on the method of detection. These proteins can be characterised qualitatively, with respect to their electrophoretic mobilities (i.e. charge and apparent molecular mass) and quantitatively, using densitometry, to determine their amounts. There has been a widespread application of 2DE in the analysis and characterisation of protein mutations for a range of organisms. This review presents examples of the use of 2DE to study naturally occurring protein mutations and polymorphisms as well as the characterisation of induced protein mutations in prokaryotes and eukaryotes. Examples are presented to illustrate the use of 2DE to detect mutations affecting the electrophoretic mobility and biosynthesis of individual proteins as well as mutations leading to global alterations in cellular protein synthesis. The advantages and disadvantages of 2DE in the detection of protein mutations are discussed.

Quantification of silver-stained proteins resolved by two-dimensional electrophoresis: genetic variability as related to abundance and solubility in two maize lines

Relative abundance and solubility of proteins from etiolated coleoptiles of maize were investigated using two-dimensional electrophoresis (2-D PAGE). Automatic quantification of silver-stained polypeptides on replicate 2-D gels made it possible to test the linearity of the relationship between spot integrated optical density and protein amount in the range from 15 micrograms to 135 micrograms per gel, in two inbred lines. A linear response was found for more than 60% of the spots in each genotype. When a linear response was found in both lines for a given spot, the slope values were similar in 94% of the cases, indicating the reliability of silver staining for polypeptide quantification. The parameters of the curves allowed the definition of protein classes of different abundances that could be compared for genetic variability between the two lines. From a comparison between standard 2-D patterns (trichloroacetic acid-acetone extracted proteins) and patterns obtained from Tris buffer extracted proteins, it appeared that 92% of the proteins visualized in the standard gels were soluble. No difference in genetic variability, either qualitative or quantitative, was evidenced between the various classes of abundance, or between soluble and insoluble proteins.