Statistical challenges in the analysis of two-dimensional difference gel electrophoresis experiments using DeCyder

MOTIVATION: The DeCyder software (GE Healthcare) is the current state-of-the-art commercial product for the analysis of two-dimensional difference gel electrophoresis (2D DIGE) experiments. Analyses complementing DeCyder are suggested by incorporating recent advances from the microarray data analysis literature. A case study on the effect of smallpox vaccination is used to compare the results obtained from DeCyder with the results obtained by applying moderated t-tests adjusted for multiple comparisons to DeCyder output data that was additionally normalized. RESULTS: Application of the more stringent statistical tests applied to the normalized 2D DIGE data decreased the number of potentially differentially expressed proteins from the number obtained from DeCyder and increased the confidence in detecting differential expression in human clinical studies.     

Fluorescence two-dimensional difference gel electrophoresis and mass spectrometry based proteomic analysis of Escherichia coli

Separation and relative quantitation of complex protein mixtures remain two of the most challenging aspects of proteomics. Here an advanced technique called fluorescence difference 2-D gel electrophoresis technology (2D-DIGE) has been applied to a model system study of the Escherichia coli proteome after benzoic acid treatment. The molecular weight and charge matched cyanine dyes enable pre-electrophoretic labelling of control and treated samples which are then mixed and run in the same gel. Pooled control and treated samples labelled with Cy trade mark 3 were used as an internal standard for both Cy5 labelled control and treated E. coli samples. Together with DeCyder trade mark imaging analysis software, more accurate quantitative analysis than conventional two-dimensional polyacrylamide gel electrophoresis was achieved. Using matrix-assisted laser desorption/ionization-time of flight and quadrupole-time of flight mass spectrometry a total of 179 differentially expressed protein spots were identified. These included enzymes, stress related and substrate (e.g. amino acids, maltose, ribose and TRP repressor) binding proteins. Of the spots analysed, 77% contained only one protein species per spot, hence the change in protein expression measured was solely attributed to the identified protein. Many membrane proteins and protein isoforms were identified indicating both adequate solubilization of E. coli samples and potential post-translational modification. The results indicate that the regulatory mechanisms following benzoic acid treatment of E. coli are far more complicated than hitherto expected.     

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.     

Semiquantitative proteomic analysis of the human spliceosome via a novel two-dimensional gel electrophoresis method

More than 200 proteins associate with human spliceosomes, but little is known about their relative abundances in a given spliceosomal complex. Here we describe a novel two-dimensional (2D) electrophoresis method that allows separation of high-molecular-mass proteins without in-gel precipitation and thus without loss of protein. Using this system coupled with mass spectrometry, we identified 171 proteins altogether on 2D maps of stage-specific spliceosomal complexes. By staining with a fluorescent dye with a wide linear intensity range, we could quantitate and categorize proteins as present in high, moderate, or low abundance. Affinity-purified human B, B(act), and C complexes contained 69, 63, and 72 highly/moderately abundant proteins, respectively. The recruitment and release of spliceosomal proteins were followed based on their abundances in A, B, B(act), and C spliceosomal complexes. Staining with a phospho-specific dye revealed that approximately one-third of the proteins detected in human spliceosomal complexes by 2D gel analyses are phosphorylated. The 2D gel electrophoresis system described here allows for the first time an objective view of the relative abundances of proteins present in a particular spliceosomal complex and also sheds additional light on the spliceosome's compositional dynamics and the phosphorylation status of spliceosomal proteins at specific stages of splicing.     

Proteome analysis of human colon cancer by two-dimensional difference gel electrophoresis and mass spectrometry

Two-dimensional difference gel electrophoresis (2-D DIGE) coupled with mass spectrometry (MS) was used to investigate tumor-specific changes in the proteome of human colorectal cancers and adjacent normal mucosa. For each of six patients with different stages of colon cancer, Cy5-labeled proteins isolated from tumor tissue were combined with Cy3-labeled proteins isolated from neighboring normal mucosa and separated on the same 2-D gel along with a Cy2-labeled mixture of all 12 normal/tumor samples as an internal standard. Over 1500 protein spot-features were analyzed in each paired normal/tumor comparison, and using DIGE technology with the mixed-sample internal standard, statistically significant quantitative comparisons of each protein abundance change could be made across multiple samples simultaneously without interference due to gel-to-gel variation. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) and tandem (TOF/TOF) MS provided sensitive and accurate mass spectral data for database interrogation, resulting in the identification of 52 unique proteins (including redundancies due to proteolysis and post-translationally modified isoforms) that were changing in abundance across the cohort. Without the benefit of the Cy2-labeled 12 sample mixture internal standard, 42 of these proteins would have been overlooked due to the large degree of variation inherent between normal and tumor samples.     

Proteome analysis of human hepatocellular carcinoma tissues by two-dimensional difference gel electrophoresis and mass spectrometry

Proteome analysis of human hepatocellular carcinoma tissues was conducted using two-dimensional difference gel electrophoresis coupled with mass spectrometry. Paired samples from the normal and tumor region of resected human liver were labeled with Cy3 and Cy5, respectively while the pooled standard sample was labeled with Cy2. After analysis by the DeCyder software, protein spots that exhibited at least a two-fold difference in intensity were excised for in-gel tryptic digestion and matrix-assisted laser desorption/ionization-time of flight mass spectrometry. A total of 6 and 42 proteins were successfully identified from the well- and poorly-differentiated samples, respectively. The majority of these proteins are related to detoxification/oxidative stress and metabolism. Three down-regulated metabolic enzymes, methionine adenosyltransferase, glycine N-methyltransferase, and betaine-homocysteine S-methyltransferase that are involved in the methylation cycle in the liver are of special interest. Their expression levels, especially, methionine adenosyltransferase, seemed to have a major influence on the level of S-adenosylmethionine (AdoMet), a vital intermediate metabolite required for the proper functioning of the liver. Recent work has shown that chronic deficiency in AdoMet in the liver results in spontaneous development of steatohepatitis and hepatocellular carcinoma, and hence the down-regulation of hepatic methionine adenosyltransferase in our hepatocellular carcinoma samples is in line with this observation. Moreover, when a comparison is made between the differentially expressed proteins from our human hepatocellular carcinoma samples and from the liver tissues of knockout mice deficient in methionine adenosyltransferase, there is a fairly good correlation between them.     

A novel experimental design for comparative two-dimensional gel analysis: two-dimensional difference gel electrophoresis incorporating a pooled internal standard

The comparison of two-dimensional (2-D) gel images from different samples is an established method used to study differences in protein expression. Conventional methods rely on comparing images from at least 2 different gels. Due to the high variation between gels, detection and quantification of protein differences can be problematic. Two-dimensional difference gel electrophoresis (Ettan trade mark DIGE) is an emerging technique for comparative proteomics, which improves the reproducibility and reliability of differential protein expression analysis between samples. In the application of DIGE different samples are labelled with mass and charge matched spectrally resolvable fluorescent dyes and are then separated on the same 2-D gel. Using an Escherichia coli lysate "spiked" with varying amounts of four different known proteins, we have tested a novel experimental design that exploits the sample multiplexing capabilities of DIGE, by including a standard sample in each gel. The standard sample comprises equal amounts of each sample to be compared and was found to improve the accuracy of protein quantification between samples from different gels allowing accurate detection of small differences in protein levels between samples.     

Comparison of proteomes in various human plasma preparations by two-dimensional gel electrophoresis

Serum or plasma can be utilized in a variety of studies targeted toward the discovery of disease biomarkers. In this study, the proteome profiles of plasma samples prepared using various anticoagulants (EDTA, heparin or citrate), were compared with those of serum using two-dimensional electrophoresis (2-DE). Proteins which evidenced different levels in the plasma and serum were screened and identified using ESI-Q-TOF MS/MS. The proteins which became detectable after the removal of fibrinogen from serum were identified as pigment epithelial differentiating factor (four spots), fetuin-like protein, and the hemopexin precursor. In particular, three proteins, pre-serum amyloid P component, plasma glutathione peroxidase precursor, and tetranectin, evidenced increased volume intensity only in the plasma samples prepared with EDTA.     

Comparative proteomic analysis for hCTLA4Ig production in transgenic rice suspension cultures using two-dimensional difference gel electrophoresis

The new technology, two-dimensional difference gel electrophoresis (2D DIGE), uses fluorescent dyes to simplify the process of detecting and matching proteins between multiple gels by allowing for the separation of up to three separate protein samples within the same gel. In this study, recombinant human cytotoxic T lymphocyte-associated antigen 4-immunoglobulin (hCTLA4lg) was produced in transgenic rice suspension cell cultures and the intracellular proteins were analyzed by 2D DIGE. The highest level of hCTLA4Ig (25.4 mg/L) was obtained five days after induction. The intracellular proteins expressed at both the growth and induction culture stages were separated and analyzed using DeCyder software. At least 2,218 spots were detected with two-fold thresholds and 95% confidence. We found that 29 spots increased and 20 spots decreased in their intensities during the production of recombinant hCTLA4Ig. In addition, the 2D Western blot of hCTLA4Ig revealed that this fusion protein was expressed in a variety of isoforms.     

Comparison of proteomes in various human plasma preparations by two-dimensional gel electrophoresis

Serum or plasma can be utilized in a variety of studies targeted toward the discovery of disease biomarkers. In this study, the proteome profiles of plasma samples prepared using various anticoagulants (EDTA, heparin or citrate), were compared with those of serum using two-dimensional electrophoresis (2-DE). Proteins which evidenced different levels in the plasma and serum were screened and identified using ESI-Q-TOF MS/MS. The proteins which became detectable after the removal of fibrinogen from serum were identified as pigment epithelial differentiating factor (four spots), fetuin-like protein, and the hemopexin precursor. In particular, three proteins, pre-serum amyloid P component, plasma glutathione peroxidase precursor, and tetranectin, evidenced increased volume intensity only in the plasma samples prepared with EDTA.     

乙型肝炎病毒疫苗免疫不应答群体中蛋白质组表达差异性的初步研究

以往研究显示,人群中约有10%无法通过接种乙型肝炎病毒(HBV)疫苗而产生保护性水平的应答。为研究HBV疫苗免疫后不应答发生的机制,本研究对108例经正常免疫程序(大连汉信产HBV疫苗,10μg/支,共接种3针)后不同应答水平人群的血浆蛋白质组差异进行分析,将其分为无应答(抗体效价检测结果为阴性)组和高应答(抗体效价约1 000mIU/ml)组。用多重亲和去除系统(MARS)亲和柱去除14种高丰度蛋白后,采用双向荧光差异电泳(2D-DIGE)与基质辅助激光解析电离串联飞行时间质谱(MALDI-TOFMS)相结合的方法,鉴定获得11种差异表达蛋白质。在这些差异表达蛋白质中,α1微球蛋白、激肽原1的表达在无应答人群中较高应答人群上调,而α1抗胰凝乳蛋白酶、维生素D结合蛋白、afamin、抗凝血酶Ⅲ和玻连蛋白表达下调。其中,α1抗胰凝乳蛋白酶、激肽原1和α1微球蛋白的差异性表达进一步得到酶联免疫吸附试验(ELISA)或蛋白免疫印迹法验证。以上蛋白质的表达情况与HBV疫苗免疫应答状况相关,可能成为检测HBV疫苗免疫应答水平的分子标记,为进一步研究免疫不应答机制奠定基础。     

Genetic analysis of human lymphocyte proteins by two-dimensional gel electrophoresis

Summary We describe a genetic polymorphism of cytosol polypeptide with mol. wt. of 38,000 detected in phytohemagglutinin (PHA)-stimulated peripheral blood lymphocytes by two-dimensional gel electrophoresis. Three different electrophoretic phenotypes (type 1-1, 2-1, 2-2) of the polypeptide have been identified in a Japanese population. Family and population studies indicate that three phenotypes are determined by two common alleles at a single autosomal locus. Since the polypeptide is mainly present in cytosol of cells, we propose that the polypeptide be temporarily designated as cytosol polypeptide with mol. wt. of 38,000 (CP 38) and that the gene for CP 38 be designated as CP 38. The gene frequencies of two common alleles (CP 38 ¹ and CP 38 ²) are 0.899 and 0.101, respectively, in a Japanese population. The data on gel filtration of cytosol proteins on a Sephadex G-100 column suggest that CP 38 exists as a dimer in the cytosol. CP 38 was observed in the wide range of different cells, including B-lymphoblastoid cells, adult skin fibroblasts, HeLa cells, and erythrocytes. In 11 out of 72 individuals, the phenotypes of CP 38 were different from those of adenosine deaminase which is similar to CP 38 in subunit size, cell distribution, and allele frequencies. These data indicate that CP 38 is a new polymorphic polypeptide encoded by an autosomal locus.     

Genetic analysis of human lymphocyte proteins by two-dimensional gel electrophoresis

Summary We describe a genetic polymorphism of cytosol polypeptide with mol.wt. of 20,000 detected in lymphocytes the arythrocytes by two-dimensional gel electrophoresis. Three different electrophoretic phenotypes (type 1-1, 2-1, and 2-2) of the polypeptide have been identified in a Japanese population. Family studies indicate that the phenotypes are determined by two common alleles at a single autosomal locus. The polypeptide is present in the cytosol of various kinds of cells and is abundant in erythrocytes. The data on a gel filtration of the erythrocyte cytosol proteins on a Sephadex G-100 column suggest that the polypeptide exists as a dimer in cells. In nine out of 79 individuals, the phenotypes of the polypeptide were different from those of glyoxalase 1 (GLO1) which has similar properties in subunit size, cell distribution, and allele frequencies. These date indicate that the polypeptide with mol. wt. of 20,000 is a new polymorphic cellular polypeptide. We propose that the polypeptide be temporarily designated as cytosol polypeptide with mol. wt. of 20,000 (CP20) and that the gene for CP20 be designated as CP20. The gene frequencies of two common alleles (CP20 ¹ and CP20 ²) are 0.955 and 0.045, respectively, in a Japanese population.     

Multiplex proteomic analysis by two-dimensional differential in-gel electrophoresis

This paper describes the use of fluorescence two-dimensional differential in-gel electrophoresis in a multiplex analysis of two distinct proteomes. As a model system, cerebral cortex tissues were analyzed from neurokinin1 receptor knockout (NK(1)R-/-) and wild type (NK(1)R+/+) mice in an attempt to identify molecular pathways involved in the function of this protein. Paired NK(1)R-/- and NK(1)R+/+ samples were labeled with fluorescent Cy3 and Cy5 dyes and electrophoresed on the same two-dimensional gels. Scanning the gels at wavelengths specific for each dye revealed the two different proteomes which were overlaid and the differences in abundance of specific protein spots were determined by the Amersham Biosciences DeCyder Differential In-gel Analysis software. A Cy2-labeled sample pool was co-electrophoresed with all Cy3- and Cy5-labeled sample pairs as an internal standard providing a link for inter-gel comparisons and for more robust statistical analysis of the data. Eight spots were found to be upregulated and two downregulated in the NK(1)R-/- mice compared to NK(1)R+/+ controls. Matrix assisted laser desorption/ionisation-time of flight (MALDI-TOF) mass fingerprinting was used to identify the proteins. The results illustrate the power of this multiplex proteomics technology and illustrate how proteomics can be used to understand gene function.     

The proteomic analysis of an adipocyte differentiated from human mesenchymal stem cells using two-dimensional gel electrophoresis

Adipose tissues play a crucial endocrine role in the control of whole body glucose homeostasis and insulin sensitivity. Considering the current substantial rise in obesity and obesity-related diseases, including diabetes, it is important to understand the molecular basis of adipocyte differentiation and its control. In this study, we have analyzed the protein expression inherent to adipogenic differentiation, by 2-DE, MALDI-TOF, and RT-PCR. This study focused on proteins that were differentially expressed by the differentiation of human mesenchymal stem cells (hMSCs) to adipocytes. We conducted 2-DE for each set of proteins in the cytosol of adipocytes that had differentiated from hMSC, in a pH range from 3-10. Thirty-two protein spots were shown to have different expression levels. Among these, eight up-regulated proteins were identified by MALDI-TOF/MS, as the following: syntaxin binding protein 3, OSBP-related protein 3, phosphodiesterase, glycophorin, immunoglobulin kappa chain variable region, peroxisome proliferative activated receptor gamma (PPAR-gamma), bA528A10.3.1 (novel protein similar to KIAA01616, isoform 1), and T cell receptor V-beta 4. Four proteins: syntaxin-3, OSBP-related protein 3, PPAR-gamma and glycophorin were associated with adipogenesis.     

Quantitative analysis of mitochondrial protein expression in methylmalonic acidemia by two-dimensional difference gel electrophoresis

Isolated methylmalonic acidemia (MMA) is a rare metabolic disease due to the deficient activity of L-methylmalonyl-CoA mutase (MCM). This mitochondrial enzyme converts L-methylmalonyl-CoA to succinyl-CoA using adenosylcobalamin (Adocbl) as cofactor. Isolated MMA is subdivided into five forms: mut MMA associated with MCM deficiency, three different defects related to mitochondrial Adocbl formation (cblA, cblB, and cblH), and cblD variant 2. We performed proteomic analysis on mitochondria from an individual with cblH/cblD disorder using 2-D DIGE to identify differentially expressed proteins in this disease. Comparative analysis of control/patient mitochondrial proteome allowed us to identify differential expression of 10 proteins. The most notable groups included proteins involved in apoptosis (cytochrome c), oxidative stress (manganese superoxide dismutase) and cell metabolism (succinyl-CoA ligase (GDP forming) and mitochondrial glycerophosphate dehydrogenase). Immunoblot analysis further validated 2-D DIGE results of two of these proteins in multiple MMA patients, suggesting that the differences in expression are a general effect in this disorder. It is feasible that the differential proteins identified in this study have a biological significance and might be related to the pathophysiology of MMA.     

Genetic analysis of human lymphocyte proteins by two-dimensional gel electrophoresis:

Summary Approximately 250 phytohemagglutinin (PHA)-stimulated peripheral blood lymphocyte polypeptides from three unrelated healthy males were compared by high-resolution two-dimensional gel electrophoresis and double-label autoradiography. Comparisons by all possible pairwise combinations of [¹⁴C]leucine-labeled proteins from an individual and [³H]leucine-labeled proteins from another revealed that only three polypeptides differed qualitatively among the three individuals. The degree of variation in lymphocyte polypeptides between different individuals was similar to that in fibroblast polypeptides reported previously. Among the three variant polypeptides, two polypeptides with mol.wt. 64,000 and mol. wt. 37,000 coexisted with a polypeptide with the same molecular weight, and they showed the behavior expected of two allelic gene products separated in the isoelectric focusing dimension by charge differences. Analysis of [¹⁴C]leucine labeled peripheral blood lymphocyte proteints, from the parents of each individual, by two-dimensional gel electrophoresis indicated that the variant polypeptides with mol. wt. 64,000 and mol. wt. 37,000 in the propositus were inherited from one of his parents. The data indicate that genetic analysis of PHA-stimulated peripheral blood lymphocyte proteins is feasible by high-resolution two-dimensional gel electrophoresis in combination with double-label autoradiography and pedigree analysis.     

Identification of human hepatocellular carcinoma-related biomarkers by two-dimensional difference gel electrophoresis and mass spectrometry

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related death throughout the world. Although hepatitis B or C viral infections are main risk factors for HCC, the molecular mechanisms leading to HCC formation have not been clarified. To reduce the mortality and improve the effectiveness of therapy, it is important to search for changes in tumor-specific biomarkers whose function may involve in disease progression and which may be useful as potential therapeutic targets. In this study, we employed two-dimensional difference gel electrophoresis (2D-DIGE) combined with nano flow liquid chromatography tandem mass spectrometry (nanoLC-MS/MS) to investigate differentially expressed proteins in HCC. For each of eight HCC patients, Cy3-labeled proteins isolated from tumor tissue were combined with Cy5-labeled proteins isolated from the surrounding nontumor tissue and separated by 2D gel electrophoresis along with a Cy2-labeled mixture of all tumor and nontumor samples as an internal standard. Thirty-four protein spots corresponding to 30 different proteins were identified by nanoLC-MS/MS as showing significant change (paired t-test, p < 0.05) in the level of expression between tumor and nontumor tissues. Sixteen proteins were up-regulated and 14 were down-regulated in HCC; they seem to play important roles in a variety of pathways including glycolysis, fatty acid transport and trafficking, amino acid metabolism, iron and xenobiotic metabolism, ethanol metabolism, cell cycle regulation, cytoskeleton, and stress. A remarkable finding is the up-regulation of 14-3-3gamma protein in HCC. 14-3-3 isoforms had been linked to carcinogenesis because they are involved in various cellular processes such as cell cycle regulation, apoptosis, proliferation, and differentiation. In conclusion, 2D-DIGE is an efficient strategy that enables us to identify differentially expressed proteins in HCC. Identification of potential biomarkers, such as the pinpointing of 14-3-3gamma in our findings, may provide further useful insights into the pathogenesis of HCC.     

Genetic analysis of human lymphocyte proteins by two-dimensional gel electrophoresis:

Summary Three different electrophoretic types (1-1, 2-1 and 2-2) of a human cellular polypeptide with molecular weight of 31000 have been identified by the analysis of PHA-stimulated peripheral blood lymphocyte proteins using high resolution two-dimensional gel electrophoresis. Family and population studies indicate that the three phenotypes of the polypeptide are determined by two common alleles at a single autosomal locus. The 31k polypeptide appears to be present as a monomer in the cytosol in a wide range of different cell types, including permanent lymphoblastoid cell lines, fibroblasts and HeLa cells. In an individual with the 31k polypeptide type 2-2, the phenotypes of adenosine deaminase and uridine monophosphate kinase were both type 1. These data indicate that the 31K polypeptide is a new polymorphic protein encoded by a new autosomal locus. It is proposed that the polypeptide and its locus be temporarily designated cytosol 31k polypeptide (C31k polypeptide) and C31P, respectively. In a Japanese population, the gene frequencies of C31P ¹ and C31P ² were 0.940 and 0.060, respectively. The C31k polypeptide type 2-2 appears to be a molecular weight variant as well as a charge variant.