- Open Access
Differential abundance analysis of mesocarp protein from high- and low-yielding oil palms associates non-oil biosynthetic enzymes to lipid biosynthesis
Proteome Science volume 13, Article number: 28 (2015)
The oil palm Elaeis guineensis Jacq. which produces the highest yield per unit land area of the oil crops is the most important commercial oil crop in South East Asia. The fleshy mesocarp of oil palm fruit, where oil is mostly derived from, contains up to 90 % dry weight of oil (one of the most concentrated in plant tissues). Hence, there is attention given to gain insights into the processes of oil deposition in this oil rich tissue. For that purpose, two-dimensional differential gel electrophoresis (DIGE) coupled with western assays, were used here to analyze differential protein levels in genetically-related high-and low-yielding oil palm mesocarps.
From the DIGE comparative analysis in combination with western analysis, 41 unique differentially accumulated proteins were discovered. Functional categorization of these proteins placed them in the metabolisms of lipid, carbohydrate, amino acids, energy, structural proteins, as well as in other functions. In particular, higher abundance of fructose-1,6-biphosphate aldolase combined with reduced level of triosephosphate isomerase and glyceraldehyde-3-phosphate dehydrogenase may be indicative of important flux balance changes in glycolysis, while amino acid metabolism also appeared to be closely linked with oil yield.
Forty-one proteins in several important biological pathways were identified as exhibiting differential in abundance at critical oil production stages. These confirm that oil yield is a complex trait involving the regulation of genes in multiple biological pathways. The results also provide insights into key control points of lipid biosynthesis in oil palm and can assist in the development of genetic markers for use in oil palm breeding programmes.
The oil palm Elaeis guineensis Jacq. produces the highest yield per unit land area of the oil crops  and has become the most important commercial oil crop in tropical South East Asia. The oil palm fruit is a drupe with bunches containing 1000–3000 fruitlets and most of the commercially extracted oil is derived from its fleshly mesocarp. The mesocarp contains up to 90 % dry weight of oil, one of the most concentrated found in plant tissues [2, 3]. Yield remains the main focus in oil palm cultivation due to the limitation of available arable land and growing demand for edible and non-edible oils. Current standard commercial oil palm material can produce approximately 4 ton/hectare/year of oil at a planting density of 148 palms/hectare. However, some plants have attained yields greater than 10 ton/hectare/year in trial plots giving an indication of potential future yields.
Recent advances in biochemical analysis have allowed researchers to develop omics techniques to investigate crop traits at a cellular level. In particular, proteomics has been used extensively to investigate the biochemical processes attribute to traits in non-oil crops [4–9], as well as processes leading to increased lipid production in oil crops [10–15]. The proteomes of two near-isogenic sunflower varieties differing in seed oil traits were comparatively analyzed by Hajduch et al. . Proteins associated with glycolytic pathway were found highly up-regulated in the high oil variety. In contrast, proteins associated to amino acid biosynthesis were up-regulated in the low oil variety. Another group of researchers, Troncoso-Ponce et al.  investigated the glycolytic initial metabolites and enzyme activities from developing seed of two different sunflower lines, of high and low oil content, during storage lipid accumulation. They discovered that the amount of sucrose produced and available for lipid synthesis was higher in high oil line. In addition, the enzymatic activities of both phosphoglycerate kinase and enolase were found notably higher in this line, allowing increased synthesis of phosphoenolpyruvate, an intermediate carbon donor for lipid biosynthesis. Another comparative analysis of the soluble proteins found in early developing germs from high-oil and normal inbred lines in maize , revealed that 3 enzymes closely related to lipid biosynthesis, namely enoyl-ACP reductase, stearoyl-ACP desaturase and acetyl-CoA C-acyl-transferase, had higher abundance in high-oil lines than in normal.
Recently, several studies that employed comparative transcriptomics and metabolomics to gain insight into the development of oil palm fruit have been reported [16–21]. Bourgis et al.  examined the differences in transcriptome and metabolome between oil palm and date palm during mesocarp development, in the hope to elucidate the mechanisms that lead to an extreme difference in carbon partitioning in them (the mesocarp of oil palm accumulates oil while the mesocarp of date palm accumulates sugars). They implicated that synthesis of fatty acids and supply of pyruvate in the plastid, rather than acyl assembly into triacylglcerol, as the main factor for the accumulation of oil in the mesocarp of oil palm. Tranbarger et al.  and Dussert et al.  investigated the transcriptional basis of lipid accumulation in the mesocarp of oil palm. A transcript, homologous to Arabidopsis seed oil transcription factor WRINKLED1, was identified to coordinate its transcript level with several fatty acids biosynthetic transcripts and high rates of lipid deposition, suggesting the mesocarp homolog is an important regulatory factor in oil biosynthesis.
The DIGE (Difference Gel Electrophoresis) analysis [22, 23], a two dimensional gel electrophoresis employing sensitive fluorescent labeling dyes, has been successfully employed in protein expression analyses in rice and sunflower [10, 24]. Although DIGE has been utilized in the studies of quantitative proteomics in these plants, but none in oil palm. In this study, we report the first comparative proteomics analysis of two distinct groups of oil palms (low-yielding, LY palms and high-yielding, HY palms) using DIGE analysis, followed by confirmation using western blot immunoassays, to identify proteins associated with yield in oil palm and provide insights into the regulation of lipid biosynthesis.
Results and discussion
Mesocarp tissue extraction protocol for 2D gel electrophoresis
Three methods were tested for protein extraction from oil palm mesocarp tissues. They were (i) the sucrose method published by He et al.  for pine needles, (ii) a modified sucrose method and (iii) the phenol/ammonium acetate method . The sucrose method was modified in this study to give a better 2D gels, evidently with more protein spots being resolved in the alkaline range (Fig. 1). Moreover, the mean number of protein spots detected in 2D gels (3 replicate gels) produced from the original sucrose method was 576 ± 18 spots whilst 738 ± 80 spots detected in 2D gels from the modified sucrose method. Successively, the modified sucrose method was compared to the phenol/ammonium acetate method. In terms of spot resolution and separation, the 2D gels produced from these 2 methods were comparable (Fig. 1) in that 729 ± 78 spots were detected from 2D gels of the third method. Consequently, the modified sucrose method was adopted to extract protein from oil palm mesocarp tissues for subsequent 2D gel electrophoresis study. Moreover, the latter method involves the usage of toxic phenol and procedurally lengthy.
Selection of plant materials and comparative studies
In order to investigate biosynthetic controls leading to significant increases in yield, the proteomes of genetically related oil palms grown in the same field and exhibiting a 2-fold difference in oil yield  were selected for this study. Mesocarp samples were collected at three key phases preceding and during lipid biosynthesis: 12 weeks after pollination (WAP) (preceding), 16 WAP (onset of lipid deposition) and 18 WAP [3, 17, 27–29]. The DIGE protein analysis was used to identify pairwise differential proteins between three high-yielding (HY) and three comparatively low-yielding (LY) plants. A total of nine analytical gels were produced (Fig. 2; see all DIGE gels in Additional file 1) from the eighteen samples and subsequently used for cross-gel Biological Variation Analysis (BVA) (DeCyder™ 2D software Version 6.5 by Amersham BioSciences): HY versus LY at 12, 16 and 18 WAP, in addition to temporal analysis of HY palms: 16 and 18 WAP versus 12 WAP (Additional file 2). Out of approximately 2000 peptide spots detected, 84 were found to be at least 1.5-fold different and have a p-value < 0.1. The p-value < 0.1 was taken to compensate the small sample size (only 3 biological replicates) in this DIGE experiment. Moreover, the DIGE experiment here functioned mainly as a screening tool for subsequent western analysis that involved bigger populations. Mass spectrometry identification using MALDI-TOF/TOF was conducted on 61 distinct spots with isoforms removed. Fifty-three protein spots were identified based on combined MS and MS/MS matching to proteins in the public domain of NCBI non-redundant database, of which 45 were unique (see Additional file 3 for MASCOT protein-peptide identification results). Subsequently, the DIGE results were cross-checked with western analysis by using a larger sample size of 8 biological replicates (Figs. 3 and 4, Additional file 4, Additional file 5, Additional file 6, Additional file 7).
The functions of the majority of the differential proteins of those having significant p-value < 0.05 from DIGE and western analyses, a total of 41 proteins, could be grouped into six main categories (Table 1) according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways database : lipid metabolism, carbohydrate metabolism, amino acid metabolism, energy metabolism, stress-related proteins and structure-related proteins. Oil yield is known to be a complex trait involving the regulation of many genes across multiple biological pathways . In this study, only 2 of the 41 proteins are known to be involved in lipid biosynthesis, while the remaining are involved in a range of biosynthetic pathways.
None of the key enzymes known to be directly involved in lipid biosynthesis or assembly was found to be differential between the HY and LY oil palm samples in this study, for example plastid fatty acid synthase and thioesterase, cytosolic acyl-CoA synthetase and the endoplasmic reticulum located enzymes of the Kennedy pathway [12, 32]. This could possibly due to the fact that the DIGE experiment here was limited to the detection of proteins within pH 4 to 7 only. The pH 4–7 range was chosen in this study because it is known that the majority of cellular proteins are found in this particular pH range. The pI of oil palm fatty acid synthases and thioesterases, and acyl-CoA synthetases, as well as the endoplastic reticulum located enzymes of the Kennedy pathway are not known. However, several Arabidopsis fatty acid synthases, fatty acid thioesterases and acyl-CoA synthetases were reported to have pI in the range above 7 [33–38]. Incidentally, the pI of some of the enzymes in the Kennedy pathway, viz the glycerol-3-phosphate O-acyltransferase, lysophosphatidic acid acyltransferase and diacylglycerol acyltransferase are reported to be in the alkaline range, respectively, pI 9.36 in flowering plant Echium pitardii , pI 9.43 in peanut  and pI 9.46 in microalga Myrmecia incisa .
Interestingly, the accumulation of phospholipase D at 16 WAP was observed to be significantly higher [1.6x, p-value < 0.05](Additional file 2) in LY palms compared to HY palms using DIGE (Table 1) and western analysis (Fig. 3). Phospholipase D cleaves the ester bond between the phosphatidyl moiety and the head group of phospholipids, but can also act as a transferase [42–44]. A study by Yu  found Arabidopsis phospholipase D knockout mutants exhibited decreased seed oil content compared to wild type, although it is possible that a reduction in the activity of this lipase at the period of highest lipid biosynthesis may be beneficial for yield.
Biotin carboxylase subunit was found to be up-regulated in the HY palms during lipid biosynthesis from 1.8-fold at 16 WAP to 3.0-fold at 18 WAP [p-value < 0.05] (Additional file 2) compared to 12 WAP when lipid biosynthesis is still relatively minimal (Fig. 4). Tranbarger et al.  also observed a higher abundance of biotin carboxylase transcripts during oil biosynthesis. Biotin carboxylase subunit is a component of acetyl-coenzyme A carboxylase . Acetyl-CoA carboxylase catalyzes the conversion of acetyl-CoA to malonyl-CoA as the initial step of de novo fatty acid synthesis occurring in the plastid . It has been demonstrated that over-expression of foxtail millet acetyl-CoA carboxylase gene in maize increases seed oil content by 25 to 65 % . Although, no differential abundance was observed between HY and LY palms in this study, the importance of biotin carboxylase to lipid production in oil palm mesocarp is clearly evident.
DIGE results identified three enzymes involved in glycolysis as being produced at a different levels between HY and LY palms; fructose-1,6-biphosphate aldolase (FBA), triosephosphate isomerase (TPI) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Results of western analysis for FBA and TPI were concordant with DIGE data (Fig. 3). During the last stages of glycolysis, FBA fragments fructose-1,6-biphosphate into glyceraldehyde-3-phosphate (GAP) and dihydroxyacetonephosphate (DHAP), while TPI acts as a housekeeping enzyme to catalyze the inter-conversion and maintain appropriate balance of GAP and DHAP concentrations. Subsequently, GAP is converted to 1,3-bisphosphoglycerate by GAPDH and then eventually leads to pyruvate. Meanwhile, glycerol-3-phosphate dehydrogenase converts DHAP to glycerol 3-phosphate, an important intermediate of the acylglycerol backbone. FBA was observed to be up-regulated at 16 and 18 WAP, concordant with lipid biosynthesis (Fig. 4), and more interestingly have higher abundance in HY palms at 16 WAP (Fig. 3). The up-regulation of FBA in HY palms is likely to contribute to increased flux through glycolysis into de novo fatty acid synthesis [10, 48]. Interestingly, DIGE analysis indicated that TPI and GAPDH were both down-regulated in HY palms at 16 WAP. These findings are consistent with observations in sunflower where TPI and GAPDH activities were found to be higher in the low-oil content sunflower lines . However, they also reported that FBA exhibited higher activity in low-oil lines. Studies of yeast  and canola  also demonstrated that reduced TPI activity causes a shift in glycolytic flux towards glycerol leading to increased yield. It is likely that a concurrent reduction in GAPDH activity would further alter the flux balance at this important branch point in the glycolysis pathway. Detailed metabolite comparisons of HY and LY oil palms provide further evidence of this flux balance change where glycerol-3-phosphate levels were observed to be 43 % higher in HY mesocarp at 16 WAP, while downstream of GAP, 3-phosphoglyceric acid levels were 26 % lower .
Receiver operating characteristic (ROC) curves  were generated using TPI protein levels detected in HY and LY palms at 12 and 16 WAP (Figs. 5). The generalized area under the curve for both time points was calculated to be 0.80 and 0.85, respectively, indicating the potential use of TPI as a yield biomarker.
Amino acid metabolism
Several enzymes involved in amino acid metabolism were found to have significantly higher levels during lipid biosynthesis (Fig. 4). Although, no significant differences were observed in their levels between HY and LY palms. A separate study by Teh et al.  has shown there are significantly higher concentrations of amino acids in the mesocarp of HY oil palms at the onset of lipid biosynthesis. Tranbarger et al.  has also reported higher amino acid concentrations in oil palm compared to its non-oil-producing relative, date palm. Taken together, these results suggest that the regulation and metabolism of amino acids before and during lipid biosynthesis may play a crucial role in fruit maturation and oil yield.
Our DIGE and western immunoassay results indicated that the enzyme caffeic acid O-methyltransferase (COMT) was more abundance in LY palms compared to HY palms at 12 WAP [2.1x, p < 0.05] and 16 WAP [3.3x, p < 0.05] (Table 1, Fig. 3 and Additional file 2). COMT is a key enzyme in the lignin biosynthetic pathway  and it is plausible that a reduction in lignin biosynthesis may be complimentary to higher lipid production in mesocarp tissue.
Abundance of several other proteins were found to be altered in the course of oil palm fruit development. One is the Vacuolar H + −ATPase (V-ATPase). It is important for secondary ion transport into the expanding vacuole and as a regulator of membrane trafficking in higher plants . In pear trees, V-ATPase activity is reported to be down-regulated during the flowering stage and markedly increases during early developmental stages and fruit growth [53, 54]. In oil palm mesocarp, the level of V-ATPase was observed to be higher in LY palms [1.9x, p < 0.05] (Additional file 2) at 16 WAP compared to HY palms using DIGE (Table 1) and western analysis (Fig. 3). The direct role of V-ATPase in contributing to the onset of oil biosynthesis in oil palm at 16 WAP is unknown. Our results here, however, suggested that this enzyme may associate with oil biosynthesis in oil palm.
Other proteins that showed differential abundance in HY palms have been reported to be associated with oil yield. Our DIGE analysis led to the identification of seven stress related proteins such as 17.6 kDa class I small heat shock protein , glutathione S-transferase  and glutathione peroxidase . Studies of maize and sunflower have also found higher level of stress response proteins in high-yielding lines [10, 12]. We also observed that catalase was higher in abundance in the mesocarp of HY palms at 18 WAP (Additional file 2). Catalase plays a role in the regulation of hydrogen peroxide and reactive oxygen species, and has been reported to be important during sunflower and maize seed development and oil biosynthesis [13, 58, 59] as well as in peach mesocarp development . Large subunit ribulose-1,5-bisphosphate caboxylase/oxygenase (Rubisco) is another protein found to have higher levels in HY palms at 18 WAP (Additional file 2). The abundance of Rubisco was also clearly correlated with lipid biosynthesis, increasing steadily from 12 to 18 WAP (Fig. 4). Rubisco is the most abundant protein found in plants and is known to play a crucial role in photosynthetic carbon assimilation. Studies on developing embryos of Brassica napus indicated that Rubisco increases the efficiency of carbon use for triacylglycerol production without the Calvin cycle [14, 61] and its expression has been linked with yield in maize [62, 63].
Several transcriptomics studies report the importance of transcription factor WRINKLED1 in oil biosynthesis in Arabidopsis, maize and oil palm [17, 18, 21, 64–67]. In oil palm mesocarp, Tranbarger et al. , Dussert et al.  and Wong et al.  all noted positive correlation between oil deposition during fruit maturity and WRINKLED1 expression; mRNA of WRINKLED1 increased from 12 WAP to 16 WAP and subsequently decreased from 18 WAP onwards. Conversely, the temporal changes of WRINKLED1 protein abundance during fruit development was not observed in the present DIGE study. Incidentally, alteration in protein level of any transcription factors was neither observed in DIGE analysis of maize differing in oil trait by Hadjuch et al. , nor in 2D gel comparative analyses of Brassica campestri oilseed development  or of Olea europaea olive fruit development . Transcription factors are generally expressed in very low amount in a cell, both as RNA and protein, because only small quantity are required for activation of cognate genes. Hence, it is not impossible that the mesocarp WRINKLED1 protein is being produced in very low concentration that it is beyond the detection limit of the DIGE technique deployed here. In addition, it is not unknown that proteins of low abundance are being masked by proteins of higher abundance in this gel-based proteomics analysis [69, 70], thus making their detection challenging.
In summary, 2D-DIGE techniques identified forty-one protein candidates as being differentially accumulated in two genetically related populations of oil palm displaying a 2-fold difference in yield. While only two of the proteins are known to be involved in oil metabolism, many are involved in other important pathways such as glycolysis, energy metabolism, amino acid metabolism and structural proteins. These results confirm that oil yield is a complex trait involving the regulation of genes in multiple biological pathways. In particular, two important glycolytic enzymes FBA and TPI exhibit differential levels in high-yielding palms, and suggest altered flux balance between DHAP and GAP may lead to increased yield. Results also showed that the regulation of amino acid metabolism appears to be closely linked to lipid biosynthesis. Better understanding of the biochemical differences leading to increased oil yield in Elaeis guineensis will aid the identification of key biosynthetic control points in addition to enabling the identification of markers for use in breeding programmes.
Sixteen sibling or half-sibling palms derived from crosses of Serdang Avenue dura and AVROS pisifera and grown in the same field at Carey Island Estate, Malaysia were selected for this study. Oil yield and fruit bunch analysis data were collected for individual plants over a 7 year period. Eight palms were identified as high-yielding with yields of 10–12 t/Ha/y, and 8 were relatively low-yielding with 4–7 t/Ha/y . Fruit bunches were harvested at 12, 16, and 18 weeks after pollination to represent key lipid biosynthesis phases [3, 17, 27–29]. Mesocarp tissues were collected, immediately frozen in liquid nitrogen then stored at−80 °C prior to use. The permission for sample collection was obtained from the owner of the estate, Sime Darby Plantation Malaysia. No specific permissions were required for the location and the field studies, because the project did not involve any protected or endangered plant species.
Mesocarp protein extractions for 2D gel electrophoresis
In order to find a suitable protein extraction method for 2D gel experiments using mesocarp protein, 3 methods were compared. They were (i) the sucrose method , (ii) a modified sucrose method and (iii) the phenol/ammonium acetate method .
The sucrose method
He et al  10 mL sucrose extraction buffer (5 % sucrose, 4 % SDS and 5 % β-mercaptoethanol) and 0.1 g PVPP were added to 1 g pulverized tissue and the mixture was swirled gently for 10 min at room temperature. After centrifugation at 10,000 g for 5 min, the supernatant was heated at 100 °C for 3 min and cooled to room temperature. Eight volumes of cold acetone were then added to the cooled supernatant. The mixture was incubated at−20 °C for at least 1 h, prior to centrifugation at 10,000 g for 10 min. The resulting pellet was re-suspended in 5 mL sucrose extraction buffer before being centrifuged at 10,000 g for 10 min. The protein pellet was washed twice in 80 % cold acetone and finally precipitated in 4 volumes of cold acetone before air drying.
The modified sucrose method
In the sucrose modified method, the steps found in the sucrose method  were followed from the start until the re-suspension of protein pellet in 5 mL sucrose extraction buffer. From here, the insoluble material was removed by centrifugation at 10,000 g for 5 min at room temperature. Sequentially, 4 volumes of cold acetone were added to the clear supernatant prior to incubation for 1 h at−20 °C. Following centrifuging for 10 min, the new pellet was washed twice in 80 % cold acetone and once in cold acetone before air drying. All centrifugal procedures in this method were carried out at 4 °C, unless stated otherwise.
The phenol/ammonium acetate method
Two mL cold 10 % TCA/acetone was added to 300 mg of pulverized mesocarp tissue, and the mixture was thoroughly mixed before centrifuging at 16,000 g for 3 min at 4 °C. The supernatant was discarded and 2 mL ammonium acetate (0.1 M) in 80 % methanol was added to the pellet and the mix was vortexed prior to being centrifuged at 16,000 g for 3 min at 4 °C. The supernatant was discarded and the pellet was washed in 80 % acetone by vortexing. The mixture was centrifuged at 16,000 g for 3 min at 4 °C and the supernatant removed. The resulting pellet was air-dried at room temperature or incubated at least 50 °C to remove residual methanol. Sequentially, the protein in the dried pellet was extracted by vortexing vigorously in 2 mL of 1:1 phenol pH 8.0/SDS buffer [SDS buffer: 0.1 M Tris–HCl (pH 8.0), 30 % sucrose, 2 % SDS, 5 % β-mercaptoethanol] for 5 min. The mixture was centrifuged at 16,000 g for 3 min at 4 °C and the phenol phase (top phase) was collected. The protein in the phenol phase was precipitated by adding 2 mL methanol containing 0.1 M ammonium acetate and stored overnight at−20 °C. The protein was pelleted by centrifuging at 16,000 g for 5 min at 4 °C before being washed in methanol by vortexing. After centrifuged at 16,000 g for 3 min at 4 °C, the pellet was washed in 80 % acetone. A final centrifugation at 16,000 x g for 3 min at 4 °C was carried out and the supernatant discarded. The final pellet was allowed to air dry.
The protein pellet obtained from each method aforementioned was then re-suspended in lysis buffer containing 30 mM Tris–HCl (pH 8.8), 7 M urea, 2 M thiourea and 4 % CHAPS. The mixture was sonicated for 30 min at room temperature and centrifuged at 14,000 rpm for 30 min at room temperature and the resulting supernatant was collected. Protein concentration of the supernatant was determined using the Bradford assay method .
The extracted proteins were used for the production of 2D gels. Briefly, for first dimension separation, 3 μg protein extract was passive-loaded into a 7 cm pH3-10 NL IPG strip (GE Healthcare), prior to being isoelectric focused for 9 kvh using the Ettan IPGphor III system (GE Healthcare). For the second dimension separation, the proteins from the strip were electrophoresed through a 12 % SDS polyacrylamide gel. The 2D gel was silver-stained  and it image was digitalized by Image Master Image Scanner III (GE Healthcare). Subsequently, spot analysis was performed on the digitalized image utilizing Image Master Platinum 7 software (GE Healthcare).
Two dimensional differential gel electrophoresis (DIGE)
Mesocarp samples of 3 HY and 3 LY oil palms (biological replicates) at 12, 16 and 18 WAP were selected for initial DIGE screening.
Setup of 2D-DIGE analytical gels
A total of 9 DIGE gels were produced for comparative analysis with sample pairing of HY and LY palm samples at each time point. An internal standard containing a pool of all 18 samples was included on the DIGE gels to aid peptide spot comparison between gels.
For each sample, 30 μg of protein was mixed with 1.0 μL of diluted CyDye, and kept on ice in the dark for 30 min. Samples from each pair of HY and LY were labeled with Cy3 and Cy5, respectively, while the internal standard was labeled with Cy2. The labeling reaction was quenched with 1.0 μL of 10 mM lysine followed by further incubation for 15 min. The HY, LY and internal standard samples were then mixed together. Two times 2-D sample buffer (8 M urea, 4 % CHAPS, 20 mg/mL DTT, 2 % Pharmalyte and trace amount of bromophenol blue), Destreak solution (100 μL) and rehydration buffer (7 M urea, 2 M thiourea, 4 % CHAPS, 20 mg/mL DTT, 1 % Pharmalyte and trace amount of bromophenol blue) were added to the labeling mix to make a total volume of 250 μL. The resulting solution was mixed thoroughly and centrifuged prior to loading onto immobilized pH 4–7 gradient gel (IPG) strips.
IEF (isoelectric focusing) and SDS-PAGE (polyacrylamide gel electrophoresis)
IEF was run following the protocol provided by Amersham BioSciences . Once IEF was completed, the IPG strips were incubated in freshly prepared equilibration buffer-1 (50 mM Tris–HCl, pH 8.8, containing 6 M urea, 30 % glycerol, 2 % SDS, trace amount of bromophenol blue and 10 mg/mL DTT) for 15 min with gentle shaking. The strips were then rinsed in freshly prepared equilibration buffer-2 (50 mM Tris–HCl, pH 8.8, containing 6 M urea, 30 % glycerol, 2 % SDS, trace amount of bromophenol blue and 45 mg/mL DTT) for 10 min with gentle shaking. Lastly, the IPG strips were rinsed in SDS-gel running buffer before transferring into 12 % SDS-gels and run at 15 °C.
Image scan and data analysis
Gel images were scanned immediately following SDS-PAGE using Typhoon TRIO (Amersham BioSciences). The scanned images were then analyzed using Image-Quant software (Version 6.0, Amersham BioSciences), followed by Biological Variation Analysis (BVA) using DeCyder™ 2D software Version 6.5 (Amersham BioSciences).
Protein Identification by MALDI TOF/TOF mass spectrometer
Three preparative gels were prepared from running non-labeled proteins (ca 200 μg each) and spots were picked up from the gels by the aid of automated Ettan Spot Picker (Amersham BioSciences). The gel spots were washed, dried, and rehydrated in trypsin digestion buffer containing modified porcine trypsin protease (Trypsin Gold, Promega) at 37 °C. The trypsinized peptides were desalted by C18 Zip-tip (Millipore), eluted in matrix solution (5 mg/ml α-cyano-4-hydroxycinnamic acid in 50 % acetonitrile, 0.1 % trifluoroacetic acid, 25 mM ammonium bicarbonate) and spotted into wells of a MALDI plate. MALDI-TOF and TOF/TOF tandem MS were performed on an ABI 4700 mass spectrometer (Applied Biosystems, Framingham, MA). The mass spectra of MALDI-TOF were acquired in reflectron positive ion mode (ca. 4000 laser shots per spectrum). Tandem MS/MS fragmentation spectra of TOF/TOF were acquired (ca. 4000 laser shots per fragmentation spectrum) from 10 most abundant ions present in each sample. Protein identities were determined by searching against the database of National Center for Biotechnology Information non-redundant (NCBInr) with the combined MS and MS/MS spectra using GPS Explorer software equipped with the MASCOT search engine (Matrix Science). Searches were performed with carbamido methylation of cysteine and oxidation of methionine residues, and with allowance of one missed cleavage of trypsin. Hits with either protein score C.I.% greater than 80 or Ion C.I.% greater than 95 were considered significant. Precursor tolerance was set to 150 ppm and MS/MS tolerance was set to 0.4 Da.
Validation of DIGE results using western analysis
Dot blot immunoassays were conducted on all 16 of the HY and LY samples collected at each of 12, 16 and 18 WAP (Additional file 4, Additional file 5, Additional file 6, Additional file 7). The specificity of antibodies used in the immunoassays was tested by 1D or 2D SDS-PAGE blots (Additional file 4; 1D and 2D SDS-PAGE methods are described in Additional file 5). Mesocarp proteins for dot blot immunoassays were extracted using the TCA-acetone method described by Carpentier et al. . The resulting protein pellet was re-dissolved and stored in urea buffer (9 M urea, 4 % CHAPS, 1 % DTT, 1 % ampholytes pH3-10, 35 mM Tris base). The concentration of the protein extract was determined using the Bradford method . Using a 386-pin replicator, protein samples were spotted onto nitrocellulose membranes. The membranes were allowed to dry overnight before incubation with a primary antibody. An alkaline phosphatase labeled secondary antibody was used for colorimetric detection with NBT/BCIP. Images of the membranes were scanned using Adobe Photoshop C84 Extended & Olympus Micro software to automatically capture and transform spot intensities into Microsoft Excel® data format. Comparisons between two populations were performed by using the Student t-test (Microsoft Excel® 2007). The receiver operating characteristic (ROC) curve was plotted using mROC (Version 1.0) by CIS Bio International.
Weeks after pollination
Murphy DJ. Oil palm: future prospects for yield and quality improvements. Lipid Technol. 2009;21:257–60.
Salas JJ, Sánchez J, Ramli US, Manaf AM, Williams M, Harwood JL. Biochemistry of lipid metabolism in olive and other oil fruits. Prog Lipid Res. 2000;39:151–80.
Sundram K, Sambanthamurthi R, Tan YA. Palm fruit chemistry and nutrition. Asia Pac J Clin Nutr. 2003;12:355–62.
Wang L, Wang X, Jin X, Jia R, Huang Q, Tan Y, et al. Comparative proteomics of Bt-transgenic and non-transgenic cotton leaves. Proteome Sci. 2015;13:1–15.
Ndimba B, Ngara R. Sorghum and sugarcane proteomics. In: Paterson AH, editor. Genomics of the Saccharinae, vol. 11. New York: Springer; 2013. p. 141–68. Plant Genetics and Genomics: Crops and Models].
Zhang XL, Qi MF, Xu T, Lu XJ, Li TL. Proteomics profiling of ethylene-induced tomato flower pedicel abscission. J Proteomics. 2015;121:67–87.
Li X, Xu W, Chowdhury M, Jin F. Comparative proteomic analysis of labellum and inner lateral petals in Cymbidium Ensifolium flowers. Int J Mol Sci. 2014;15:19877.
Jiang C, Cheng Z, Zhang C, Yu T, Zhong Q, Shen JQ, et al. Proteomic analysis of seed storage proteins in wild rice species of the Oryza genus. Proteome Sci. 2014;12:1–12.
Budak H, Akpinar BA, Unver T, Turktas M. Proteome changes in wild and modern wheat leaves upon drought stress by two-dimensional electrophoresis and nanoLC-ESI-MS/MS. Plant Mol Biol. 2013;83:89–103.
Hajduch M, Casteel JE, Tang S, Hearne LB, Knapp S, Thelen JJ. Proteomic analysis of near-isogenic sunflower varieties differing in seed oil traits. J Proteome Res. 2007;6:3232–41.
Troncoso-Ponce MA, Garcés R, Martínez-Force E. Glycolytic enzymatic activities in developing seeds involved in the differences between standard and low oil content sunflowers (Helianthus annuus L.). Plant Physiol Biochem. 2010;48:961–5.
Liu Z, Yang X, Fu Y, Zhang Y, Yan J, Song T, et al. Proteomic analysis of early germs with high-oil and normal inbred lines in maize. Mol Biol Rep. 2009;36:813–21.
Tsaftaris AS, Scandalios JG. Comparison of the glyoxysomes and the glyoxysomal enzymes in maize lines with high or low oil content. Plant Physiol. 1983;71:447–50.
Hajduch M, Casteel JE, Hurrelmeyer KE, Song Z, Agrawal GK, Thelen JJ. Proteomic analysis of seed filling in Brassica napus. Developmental characterization of metabolic isozymes using high-resolution two-dimensional gel electrophoresis. Plant Physiol. 2006;141:32–46.
Bianco L, Alagna F, Baldoni L, Finnie C, Svensson B, Perrotta G. Proteome regulation during Olea europaea fruit development. PLoS One. 2013;8:e53563.
Bourgis F, Kilaru A, Cao X, Ngando-Ebongue G-F, Drira N, Ohlrogge JB, et al. Comparative transcriptome and metabolite analysis of oil palm and date palm mesocarp that differ dramatically in carbon partitioning. Proc Natl Acad Sci. 2011;108:12527–32.
Tranbarger TJ, Dussert S, Joet T, Argout X, Summo M, Champion A, et al. Regulatory mechanisms underlying oil palm fruit mesocarp maturation, ripening, and functional specialization in lipid and carotenoid metabolism. Plant Physiol. 2011;156:564–84.
Wong Y, Kwong Q, Lee H, Ong C, Mayes S, Chew F, et al. Expression comparison of Oil biosynthesis genes in Oil palm mesocarp tissue using custom array. Microarrays. 2014;3:263.
Teh HF, Neoh BK, Hong PL, Low YS, Ng LM, Ithnin N, et al. Differential metabolite during fruit development in high-yielding oil palm mesocarp. PLoS One. 2013;8:e61344.
Neoh BK, Teh HF, Ng TL, Tiong SH, Thang YM, Ersad MA, et al. Profiling of metabolites in oil palm mesocarp at different stages of oil biosynthesis. J Agric Food Chem. 2013;61:1920–7.
Dussert S, Guerin C, Andersson M, Joët T, Tranbarger TJ, Pizot M, et al. Comparative transcriptome analysis of three Oil palm fruit and seed tissues that differ in oil content and fatty acid composition. Plant Physiol. 2013;162:1337–58.
Mackintosh JA, Choi H-Y, Bae S-H, Veal DA, Bell PJ, Ferrari BC, et al. A fluorescent natural product for ultra sensitive detection of proteins in one-dimensional and two-dimensional gel electrophoresis. Proteomics. 2003;3:2273–88.
Chevalier F. Highlights on the capacities of “Gel-based” proteomics. Proteome Sci. 2010;8:1–10.
Teshima R, Nakamura R, Satoh R, Nakamura R. 2D-DIGE analysis of rice proteins from different cultivars. Regul Toxicol Pharmacol. 2010;58:S30–5.
He C-y, J-g Z, A-g D, J-y Y, D-s Z. Comparison of methods for protein extraction from pine needles. Forestry Stud China. 2005;7:20–3.
Wang W, Vignani R, Scali M, Cresti M. A universal and rapid protocol for protein extraction from recalcitrant plant tissues for proteomic analysis. Electrophoresis. 2006;27:2782–6.
Othman A, Lazarus C, Fraser T, Stobart K. Cloning of a palmitoyl-acyl carrier protein thioesterase from oil palm. Biochem Soc Trans. 2000;28:619–22.
Asemota O, Shah FH. Detection of mesocarp oleoyl-thioesterase gene of the South American oil palm Elaeis oleifera by reverse transcriptase polymerase chain reaction. Afr J Biotechnol. 2004;3:595–8.
Shaarani SM, Cardenas-Blanco A, Amin MHG, Soon NG, Hall LD. Monitoring development and ripeness of oil palm fruit (Elaeis guineensis) by MRI and bulk NMR. Int J Agric Biol. 2010;12:101–5.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Dudley JW, Lambert RJ: 100 Generations of Selection for Oil and Protein in Corn. In Plant Breeding Reviews. Volume 24. Edited by Janick J. Hoboken, New Jersey: John Wiley & Sons, Inc.; 2010: 79–110.
Ramli US, Baker DS, Quant PA, Harwood JL. Use of control analysis to study the regulation of plant lipid biosynthesis. Biochem Soc Trans. 2002;30:1043–6.
Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, et al.: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Research 2012;40:D1202-12102, 40:D1202-1210.
The Arabidopsis Information Resource (TAIR) [https://www.arabidopsis.org/servlets/TairObject?id=137515&type=locus]
The Arabidopsis Information Resource (TAIR) [https://www.arabidopsis.org/servlets/TairObject?id=137248&type=locus]
The Arabidopsis Information Resource (TAIR) [https://www.arabidopsis.org/servlets/TairObject?id=38410&type=locus]
The Arabidopsis Information Resource (TAIR) [https://www.arabidopsis.org/servlets/TairObject?id=27844&type=locus]
The Arabidopsis Information Resource (TAIR) [https://www.arabidopsis.org/servlets/TairObject?id=37115&type=locus]
Mañas-Fernández A, Li-Beisson Y, Alonso D, García-Maroto F. Cloning and molecular characterization of a glycerol-3-phosphate O-acyltransferase (GPAT) gene from Echium (Boraginaceae) involved in the biosynthesis of cutin polyesters. Planta. 2010;232:987–97.
Chen S-L, Huang J-Q, Lei Y, Ren X-P, Wen Q-G, Chen Y-N, et al. Cloning and expression analysis of lysophosphatidic acid acyltransferase (LPAT) encoding gene in peanut. Acta Agron Sin. 2012;38:245–55.
Chen C-X, Sun Z, Cao H-S, Fang F-L, Ouyang L-L, Zhou Z-G. Identification and characterization of three genes encoding acyl-CoA: diacylglycerol acyltransferase (DGAT) from the microalga Myrmecia incisa Reisigl. Algal Res. 2015;12:280–8.
Heller M. Phospholipase D. Adv Lipid Res. 1978;16:267–326.
Wang X. Multiple forms of phospholipase D in plants: the gene family, catalytic and regulatory properties, and cellular functions. Prog Lipid Res. 2000;39:109–49.
Qin C, Wang X. The Arabidopsis phospholipase D family. Characterization of a calcium-independent and phosphatidylcholine-selective PLDζ1 with distinct regulatory domains. Plant Physiol. 2002;128:1057–68.
Yu Y. Investigation of potential roles of phospholipase D in Arabidopsis thaliana seed oil accumulation. Dept Botany: University of British Columbia; 2009.
Shorrosh BS, Roesler KR, Shintani D, van de Loo FJ, Ohlrogge JB. Structural analysis, plastid localization, and expression of the biotin carboxylase subunit of acetyl-coenzyme A carboxylase from tobacco. Plant Physiol. 1995;108:805–12.
Dong Z, Zhao H, He J, Huai J, Lin H, Zheng J, et al. Overexpression of a foxtail millet Acetyl-CoA carboxylase gene in maize increases sethoxydim resistance and oil content. Afr J Biotechnol. 2011;10:3986–95.
Hajduch M, Matusova R, Houston NL, Thelen JJ. Comparative proteomics of seed maturation in oilseeds reveals differences in intermediary metabolism. Proteomics. 2011;11:1619–29.
Capitanio D, Merico A, Ranzi BM, Compagno C. Effects of the loss of triose phosphate isomerase activity on carbon metabolism in Kluyveromyces lactis. Res Microbiol. 2002;153:593–8.
Swets J. Measuring the accuracy of diagnostic systems. Science. 1988;240:1285–93.
Sewalt V, Ni W, Blount JW, Jung HG, Masoud SA, Howles PA, et al. Reduced lignin content and altered lignin composition in transgenic tobacco down-regulated in expression of L-phenylalanine ammonia-lyase or cinnamate 4-hydroxylase. Plant Physiol. 1997;115:41–50.
Schumacher K, Krebs M. The V-ATPase: small cargo, large effects. Curr Opin Plant Biol. 2010;13:724–30.
Shiratake K, Kanayama Y, Maeshima M, Yamaki S. Changes in H(+)-pumps and a tonoplast intrinsic protein of vacuolar membranes during the development of pear fruit. Plant Cell Physiol. 1997;38:1039–45.
Shiratake K, Martinoia E. Transporters in fruit vacuoles. Plant Biotech. 2007;24:127–33.
Efeoglu B. Heat shock proteins and heat shock response in plants. Gazi Univ J Sci. 2009;22:67–75.
Jain M, Ghanashyam C, Bhattacharjee A. Comprehensive expression analysis suggests overlapping and specific roles of rice glutathione S-transferase genes during development and stress responses. BMC Genomics. 2010;11:73.
Rodriguez Milla MA, Maurer A, Huete AR, Gustafson JP. Glutathione peroxidase genes in Arabidopsis are ubiquitous and regulated by abiotic stresses through diverse signaling pathways. Plant J. 2003;36:602–15.
Kibinza S, Bazin J, Bailly C, Farrant JM, Corbineau F, El-Maarouf-Bouteau H. Catalase is a key enzyme in seed recovery from ageing during priming. Plant Sci. 2011;181:309–15.
Polidoros AN, Scandalios JG. Response of the maize catalases to light. Free Radic Biol Med. 1997;23:497–504.
Nilo R, Saffie C, Lilley K, Baeza-Yates R, Cambiazo V, Campos-Vargas R, et al. Proteomic analysis of peach fruit mesocarp softening and chilling injury using difference gel electrophoresis (DIGE). BMC Genomics. 2010;11:43.
Schwender J, Goffman F, Ohlrogge JB, Shachar-Hill Y. Rubisco without the Calvin cycle improves the carbon efficiency of developing green seeds. Nature. 2004;432:779–82.
Loza-Tavera H, Serrano B. J.D. M, Ortega-Delgado ML, Sanchez-de-Jimenez E: C02-fixation enzymes in maize (Zea mays) and grain yield. Can J Bot. 1987;65:607–11.
Loza-Tavera H, Martinez-Barajas E, Sanchez-de-Jimenez E. Regulation of ribulose-1,5-bisphosphate carboxylase expression in second leaves of maize seedlings from low and high yield populations. Plant Physiol. 1990;93:541–8.
Baud S, Mendoza MS, To A, Harscoët E, Lepiniec L, Dubreucq B. WRINKLED1 specifies the regulatory action of LEAFY COTYLEDON2 towards fatty acid metabolism during seed maturation in Arabidopsis. Plant J. 2007;50:825–38.
Cernac A, Benning C. WRINKLED1 encodes an AP2/EREB domain protein involved in the control of storage compound biosynthesis in Arabidopsis. Plant J. 2004;40:575–85.
Pouvreau B, Baud S, Vernoud V, Morin V, Py C, Gendrot G, et al. Duplicate maize Wrinkled1 transcription factors activate target genes involved in seed Oil biosynthesis. Plant Physiol. 2011;156:674–86.
Shen B, Allen WB, Zheng P, Li C, Glassman K, Ranch J, et al. Expression of ZmLEC1 and ZmWRI1 increases seed Oil production in maize. Plant Physiol. 2010;153:980–7.
Li W, Gao Y, Xu H, Zhang Y, Wang J. A proteomic analysis of seed development in Brassica campestri L. PLoS One. 2012;7:e50290.
Corthals GL, Wasinger VC, Hochstrasser DF, Sanchez JC. The dynamic range of protein expression: a challenge for proteomic research. Electrophoresis. 2000;21:1104–15.
Gygi SP, Corthals GL, Zhang Y, Rochon Y, Aebersold R. Evaluation of two-dimensional gel electrophoresis-based proteome analysis technology. Proc Natl Acad Sci U S A. 2000;97:9390–5.
Bradford MM. A dye-binding method for the determination of microgram quantities of protein. Anal Biochem. 1976;72:248–54.
Blum H, Beier H, Gross H. Improved silver staining of plant proteins, RNA and DNA in polyacrylamide gels. Electrophoresis. 1987;8:93–9.
GE Healthcare. 2-D electrophoresis: principles and methods. 2004.
Carpentier SC, Witters E, Laukens K, Deckers P, Swennen R, Panis B. Preparation of protein extracts from recalcitrant plant tissues: an evaluation of different methods for two-dimensional gel electrophoresis analysis. Proteomics. 2005;5:2497–507.
The authors acknowledge Sime Darby Plantation, Malaysia for funding the project. The authors offer thanks to their colleagues at Sime Darby R & D Centres, especially to Claude Bakoume and Mohaimi Mohamed for ascertainment and sample attainment of the low/high yielding palms. Applied Biomics, Inc. of San Francisco, California is recognized for their assistance in DIGE work.
The authors have declared no competing interests in this study.
TEKO, WCY, LDJD, BZN, FTC and HK participated in the conception and designing the study. TEKO, LDJD, WCY and BZN performed the 2D gel and DIGE experiments. WCY, BZN, FCL and AMO carried out the western immunoassays. TEKO, WCY, FCL, FTC and DRA help to draft the manuscript. All authors read and approved the final manuscript.
Two dimensional DIGE gels produced for spot analysis using 13 cm IEF strips of pH 4–7 with 20-kVh focusing. Gel images of DIGE experiments labeled with Cy3 and Cy5 and their overlays. (PDF 2904 kb)
Data of Biological Variation Analysis (BVA) on DIGE gels. Results of Biological Variation Analysis (BVA) on DIGE gels are given in a table. (PDF 95 kb)
Protein-peptide summary of MASCOT protein identification. Protein-peptide summary of MASCOT protein identification. (XLSX 396 kb)
Validation of DIGE results by western analysis. Results of western immunoblots carried out in this study. (PDF 193 kb)
Protocols for electrophoresis and immuno-blotting for western analysis, and antibodies used in this study. Methods for 1D and 2D-SDS PAGE for immunoblots analysis and the list of antibodies used in this study. (PDF 77 kb)
Analysis of dot intensity of immunoblots for comparison of protein levels in mesocarp of high yielding and low yielding oil palms at 12, 16 and 18 week after pollination. Analysis of dot intensity of immunoblots for comparison of protein levels in mesocarp of high yielding and low yielding oil palms at 12, 16 and 18 week after pollination. (PDF 120 kb)
Dot intensity analysis of immunoblots for time points comparison (12 vs 16 week after pollination and 12 vs 18 week after pollination) of mesocarp protein levels in high yielding oil palms. Dot intensity analysis of immunoblots for time points comparison (12 vs 16 week after pollination and 12 vs 18 week after pollination) of mesocarp protein levels in high yielding oil palms (PDF 74 kb)
About this article
Cite this article
Ooi, T.E.K., Yeap, W.C., Daim, L.D.J. et al. Differential abundance analysis of mesocarp protein from high- and low-yielding oil palms associates non-oil biosynthetic enzymes to lipid biosynthesis. Proteome Sci 13, 28 (2015). https://doi.org/10.1186/s12953-015-0085-2
- Oil palm
- Elaeis guineensis
- Lipid biosynthesis
- Plant proteome