Figures
Abstract
Grain filling during the linear phase contributes most of the dry matter accumulated in the maize kernel, which in turn determines the final grain yield. Endosperms and embryos of three elite maize hybrids (Zhengdan 958, Nongda 108, and Pioneer 335) were sampled 17, 22, 25, and 28 days after pollination, during the linear phase of grain filling, for proteomic analysis to explore the regulatory factors critical for grain filling rate. In total, 39 and 43 protein spots that showed more than 2-fold changes in abundance at P<0.01 between any two sampling stages in the endosperm and embryo were analyzed by protein mass spectrometry. The changing patterns in expression index of these proteins in the endosperm were evenly distributed, whereas up-regulation patterns predominated (74%) in the embryo. Functional analysis revealed that metabolism was the largest category, represented by nine proteins in the endosperm and 12 proteins in the embryo, of the proteins that significantly changed in abundance. Glycolysis, a critical process both for glucose conversion into pyruvate and for release of free energy and reducing power, and proteins related to redox homeostasis were emphasized in the endosperm. Additionally, lipid, nitrogen, and inositol metabolism related to fatty acid biosynthesis and late embryogenesis abundant proteins were emphasized in the embryo. One protein related to cellular redox equilibrium, which showed a more than 50-fold change in abundance and was co-localized with a quantitative trait locus for grain yield on chromosome 1, was further investigated by transcriptional profile implying consistent expression pattern with protein accumulation. The present results provide a first step towards elucidation of the gene network responsible for regulation of grain filling in maize.
Citation: Jin X, Fu Z, Ding D, Li W, Liu Z, Tang J (2013) Proteomic Identification of Genes Associated with Maize Grain-Filling Rate. PLoS ONE 8(3): e59353. https://doi.org/10.1371/journal.pone.0059353
Editor: Ive De Smet, University of Nottingham, United Kingdom
Received: October 12, 2012; Accepted: February 13, 2013; Published: March 19, 2013
Copyright: © 2013 Jin et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the Major State Basic Research Development Program of China (2012CB723001,http://www.nsfc.gov.cn.) and a grant from the National High Technology Research and Development Program of China (no 2012AA10A305 http://www.863.gov.cn/) and Key Base Research Project of Henan Province. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Maize (Zea mays L.) is one of the most important cereal crops cultivated worldwide, and serves as a source of food, feed and fuel. Among the three key components of grain yield in maize, kernel weight is crucial. The embryo, which requires the coordinated development of multiple tissues and cell types, and the endosperm, a starch-rich tissue that supports the embryo during germination, are the main contributors to final kernel weight. Grain productivity is the result of grain-filling rate and duration, which contribute greatly to the accumulation and distribution of dry matter in the kernel [1]. During the effective grain-filling period, the grain-filling rate is positively correlated with grain weight [2]–[4]. Thus, improvement of the grain-filling rate is an important objective of breeders in order to achieve high yields, especially under stress conditions [5]–[7].
In maize, accumulation of kernel dry matter influenced by grain filling rate is a continuous process that can be divided into three characteristic stages [8]. The first stage occurs from fertilization to 12 days after pollination (DAP). During this lag phase, the cells expand rapidly in number and size, which determines the ‘storeroom’ for accumulation of photosynthate [9], [10]. The cell number might have a greater impact on final grain yield than cell size [11]–[15]. As the primary source of maize kernel dry matter [16], carbohydrates from leaf photosynthates are rapidly transported to the developing grain during the linear phase of grain filling. The linear phase ranges from 12 to 40 DAP with the most rapid grain filling occurring between 21 and 25 DAP. Notably, the conversion of imported sucrose and amino acids into starch and storage proteins during this stage can account for about 90% of kernel total dry matter [17], [18]. The final phase is the ‘desiccation–maturation’ stage, which occurs from 40 to 70 DAP. In this stage, a decreasing grain-filling rate culminates in physiological maturity and black-layer formation, and is related to seed dormancy [18], [19]. The duration of each stage is variable depending on genetic background and environmental temperature [4].
Given that the linear phase has important implications for final maize grain yield, its underlying molecular mechanisms have been studied at different levels. Liu et al. evaluated the grain-filling rate of a set of 203 recombinant inbred lines at two locations over two years using 217 simple sequence repeat (SSR) markers and identified six stable quantitative trait loci (QTL) for grain-filling rate [20]. Méchin et al. established a proteome reference map for an INRA F2 maize inbred line of the endosperm at 14 DAP by means of two-dimensional (2-D) gel electrophoresis and protein identification with liquid chromatography–mass spectrometry (MS) [17]. Prioul et al. further investigated the endosperm transcriptome of this F2 inbred line with cDNA libraries constructed at the lag phase (10 DAP), onset of carbohydrate storage (14 DAP), and maximum starch accumulation rate (21 DAP). The 10 DAP transcriptome and 14 DAP proteome showed good agreement [18]. Although some QTL and important functional protein classes were identified during the grain-filling stage, derivation of the endosperm proteome from a single inbred line is inadequate to understand the molecular mechanisms and crucial regulatory factors responsible for grain filling. Hybrid maize genotypes, which show a grain-filling rate distinct from that of inbred lines, predominate in terms of cultivated area and grain yield worldwide compared with inbred lines. Thus, dissection of the molecular mechanism of grain-filling rate in hybrids rather than inbred lines is more important for improvement of grain yield by genetic manipulation in maize. The embryo, which develops from the fertilized egg and thus is genomically distinct from the maternal tissue-derived endosperm, predominantly accumulates fatty acids and also contributes importantly to kernel weight. However, in cereals proteomic analysis of the embryo during the linear phase has been reported only in rice [21]. Thus, to elucidate the molecular mechanisms responsible for grain yield in maize, integration of genomic and proteomic assays for endosperm and embryo of hybrid maize should be undertaken synchronously.
As an important exploratory technique at the translational level, proteomic technology, which combines the resolution of 2-D gel electrophoresis with the sensitivity of MS, was used in the present study to dissect the specific proteins that accumulate in the endosperm and embryo during the linear phase of grain filling in three elite hybrids of maize (Zhengdan 958, Nongda 108, and Pioneer 335), which are widely cultivated in China. The main objectives were: (1) to obtain comparative information on proteins differentially accumulated between the endosperm and embryo during the linear phase of grain filling; and (2) to identify potential candidate genes that influence grain filling in maize.
Results
Analysis of Grain-filling Rate
During the entire period of dry matter accumulation (10 to 50 DAP), the three elite hybrids showed significant differences in final grain dry weight (P = 0.006) (Figure 1B) and grain-filling rate (P = 0.004) (Figure 1A). In addition, the grain-filling rate at 17 and 33 DAP (Figure 1A), and grain dry weight at 28 DAP (Figure 1B) also showed significant differences between the three hybrids. Grain-filling rate dynamics was different for the three tested hybrids but with the same peak at 25 DAP (Figure 1A). Zhengdan 958 showed a steep increase before 25 DAP and a more moderate decrease after 28 DAP in grain-filling rate, which was significantly distinguish from that of Nongda 108 and Pioneer 335 with distinct grain-filling rate from 25 to 40 DAP (Figure 1A). These results implied that dry matter accumulation might mainly derive from rapid grain filling before 25 DAP and a relatively uniform decline after 25 DAP in Zhengdan 958; in contrast, dry matter accumulation was mainly attributable to the comparatively high grain-filling rate after 25 DAP, especially between 25 and 40 DAP, in Pioneer 335.
Identification of Differentially Expressed Proteins
The endosperm and embryo of grains were sampled at 17, 22, 25, and 28 DAP for proteomic analysis within a linear gradient of pH 4–7 on the basis of the observed changes in grain-filling rate. At these respective sampling stages, 568, 655, 583, and 507 protein spots in the endosperm and 638, 523, 614, and 463 protein spots in the embryo were identified for Zhengdan 958; 565, 437, 527, and 498 protein spots in the endosperm and 691, 704, 699, and 345 protein spots in the embryo were identified for Nongda 108; and 439, 501, 426, and 383 protein spots in the endosperm and 600, 707, 675, and 802 protein spots in the embryo were identified for Pioneer 335 (Table 1). Compared with the embryo, a greater number of protein spots were identified at 25 DAP and fewer protein spots were identified at 17 DAP for each genotype in the endosperm. Furthermore, analysis of variance (ANOVA) and t-test analyses revealed that 281 and 183, 278 and 209, and 244 and 193 protein spots showed significantly differential expression in the endosperm and embryo of Zhengdan 958, Nongda 108, and Pioneer 335, respectively, at the four sampling stages (P<0.01) (Table 1). Among these spots, 39 spots in the endosperm (17 for Zhengdan 958, 10 for Nongda 108, and 12 for Pioneer 335) and 43 spots in the embryo (17 for Zhengdan 958, 14 for Nongda 108, and 12 for Pioneer 335), all of which showed greater than 2-fold changes in abundance at P<0.01, were excised and analyzed by MS for protein identification (Figure 2; Figure S1).
Maps from left to right represent expressed proteins at 17, 22, 25, and 28 days after pollination in the endosperm (A) and embryo (B) of Zhengdan 958, respectively. Differentially expressed protein spots at any two stages with more than 2-fold changes in expression are indicated by arrows with numeral.
Changing Pattern and Functional Analysis of Differentially Expressed Proteins
The differentially expressed proteins in the endosperm (39 protein spots) and embryo (43 protein spots) were revealed in seven groups (Figure S2). Compared with the endosperm, continuously up-regulated proteins formed the largest group (group I) represented by 13 out of 39 proteins in the embryo. Proteins with the accumulation peak at 25 DAP (Group III) were also frequent in the embryo, which was represented by 12 out of 39 proteins. In contrast, expression patterns in the endosperm were more uniformly distributed over the sampling period (Table S1).
On the basis of functional annotations, seven categories of proteins in the endosperm (Figure 3A) and eight categories in the embryo were screened (Figure 3B). The largest functional category was metabolism, which was represented by nine and twelve proteins in the endosperm and embryo, respectively. This category was further subdivided into five and seven subcategories in order to analyze in more detail dynamic changes in metabolism for the endosperm and embryo, respectively. Proteins involved in glycolysis (five proteins) in the endosperm, and proteins involved in lipid metabolism (three proteins) and nitrogen metabolism (three proteins) in the embryo were the predominant subcategories of the metabolism category. Defense proteins, the second-largest category, were represented by eight proteins in the endosperm and ten proteins in the embryo. In the endosperm, the remaining proteins were classified as protein synthesis (six proteins), protein destination (four proteins), cellular communication proteins (four proteins), energy proteins (three proteins), and ionic homeostasis proteins (two proteins). Proteins associated with development (seven proteins), cellular biogenesis (six proteins), transport facilitation proteins (three proteins), energy proteins (two proteins), cellular communication proteins (one protein) and ionic homeostasis proteins (one protein) were the remaining categories in the embryo.
A: the protein category in the endosperm; B: the protein category in the embryo.
In the endosperm, transcription proteins and eukaryotic translation initiation factor 5A (Table 2) were down-regulated. Elongation factor 1-delta was up-regulated in abundance during grain filling. The expression changes of transcription, translation initiation and elongation proteins during the linear phase of grain filling were in accordance with the protein synthesis process (Figure 4). Proteins involved in protein destination were up-regulated. Energy-related proteins that showed continuous down-regulation were identified only in Zhengdan 958. In the defense category, proteins responsive to environmental factors were up-regulated but the reverse applied to redox homeostasis proteins. In the largest metabolism category, two major glycolysis proteins were identified, namely triose-phosphate isomerase in Nongda 108 (spot no. 27) and Pioneer 335 (no. 39), showing a differential expression pattern, and glyceraldehyde-3-phosphate dehydrogenase in Zhengdan 958 (no. 7) and Nongda 108 (no. 25) showing an identical up-regulation expression pattern (Figure 4).
In the embryo, almost all proteins involved in development, defense, transport facilitation, cellular biogenesis, nitrogen metabolism, and lipid metabolism were up-regulated from 17 DAP to 25 or 28 DAP. Notably, most development-related proteins were late embryogenesis abundant (LEA) proteins and defense-related proteins comprised redox homeostasis regulation or antifungal proteins (Table 2).
Transcriptional Profile of Protein sdh1 (Sorbitol Dehydrogenase Homolog 1, Spot no. 11)
Among the protein spots identified by MS, protein sdh1 was continuously up-regulated and showed the greatest difference in accumulation levels of all proteins–a more than 50-fold increase in Zhengdan 958 from 17 to 28 DAP. The transcriptional profile of the gene encoding this protein was assayed to explore the potential molecular mechanism at the transcriptional level. A similar pattern, with the expression peak at 22 DAP, was observed in each of the hybrids (Figure 5). In addition, the target sites of miRNA168 were identified in the coding gene with PlantPAN (http://plantpan.mbc.nctu.edu.tw/gene_group/index.php) programme.
Integration with Genomic Information
To validate the potential implications of these functional annotations, the present results and those of a previous genetic analysis [20] were integrated for grain-filling rate (Table S3). Seven endosperm proteins, including one GenBank accession (gi|195618512) identified twice in the present investigation, were mapped to five QTL intervals on chromosomes (Chr.) 1, 2, 7, and 9. Nine distinct embryo proteins were mapped to three QTL intervals on Chr. 3, 8, and 9. Protein spot no. 11 was co-localized with a QTL for grain yield (qGY1b) [20] on Chr. 1 between the SSR markers umc1396 and umc1122.
Discussion
Key Period of Grain Filling
Grain yield is the result of accumulation of carbohydrates from leaf photosynthates during grain filling and non-structural carbohydrates in culms and leaf sheaths [22], [23]. Previous studies have demonstrated the importance of the former during the linear phase of grain filling for final yield [16]. In cereal crops, the effective and compensatory behavior of the duration and rate of grain filling influence the accumulation in the grain of carbohydrates from leaf photosynthates. Different from grain-filling duration, the grain-filling rate is positively correlated with, and critical for, the final grain yield [3], [4], [20], [23]. Therefore, the rate of carbohydrate transportation during the linear phase of grain filling has a fundamental effect on final grain yield based on the particular grain volume formed during the lag phase. Regardless of the genetic background, the peak grain-filling rate in maize inbred lines and hybrids occurs between 21 and 25 DAP [4], [20]. In the present study, the peak grain-filling rate was observed at 25 DAP for the three tested hybrids (Figure 1A). On the basis of changes in the grain-filling rate for the tested hybrids and data from a previous study [20], the critical period that influences dry matter accumulation was determined to be from 17 to 28 DAP. Four sampling stages were assayed both for the endosperm and embryo of each hybrid by 2-D gel electrophoresis in conjunction with MS identification, which identified some key factors that influence the carbohydrate transportation to the sink tissue (Table 2). Among the identified proteins, many were annotated as the ones involved in metabolic pathways, which implied to be important contributors to the grain filling rate and final yield. Exploration of the expression pattern and function of these proteins during the linear phase of grain filling should shed light on the genetic mechanism of kernel dry matter accumulation and contribute to understanding the regulation of grain development. The expression profiling for gene encoding sorbitol dehydrogenase homolog 1 implied that a specific signal might induce the high expression of this gene at 22 DAP (Figure 5).
The Remobilization of Nutrients during Grain Filling Course
Two phases are involved in maize grain filling: 1) the production of carbon and nitrogen from carbohydrates, fats and proteins derived from photosynthetic and/or non-photosynthetic organs; 2) the synthesis and accumulation of the three nutrients-carbohydrates, fats and proteins in endosperm and embryo. This remobilization of nutrients would directly influence the grain filling rate and dry weight in different hybrids involves. Distinct metabolic processes should be integrated into the remobilization of nutrients, which requires a coordinating protein network [24].
In maize kernels, carbohydrates such as glucose, fructose and glutamine, as sources of carbon and nitrogen [25], are used by the endosperm and embryo to generate and accumulate different biomass products, which requires that these two tissues have different patterns of metabolism. Transportation of carbohydrates is a complicated process, and the synthesis of biomass products both in the endosperm and embryo requires several metabolic pathways to produce ATP and reducing power. Glycolysis/Gluconeogenesis pathway located in the cytoplasm of the cell is important for glucose conversion into pyruvate/the regeneration of a carbon source. The former is accompanied by the release of free energy that is used to form the high-energy compounds ATP and NADH. The enzymes involved in glycolysis pathway determine the dynamic balance of metabolites. Triose-phosphate isomerase plays a central role in glycolysis and influences efficient energy production, and has been identified in numerous plant studies [17]–[19], [26]–[29]. Triose-phosphate isomerase was identified in the endosperm of Zhengdan 958 (spot no. 17), Nongda 108 (no. 27), and Pioneer 335 (no. 39), which had different expression patterns. Glyceraldehyde-3-phosphate dehydrogenase, which functions in the sixth step of glycolysis to break down glucose for energy and carbon molecules, was also identified in the endosperm of Zhengdan 958 (spot no. 7) and Nongda 108 (no. 25), showing continuously up-regulated expression during grain filling course. This protein was identified only in the embryo of Zhengdan 958 with the least accumulation level at 22 DAP (Table 2, Table S1). Malate dehydrogenase involved in gluconeogenesis [30] is important for the generation of malate which is used to get oxaloacetate for tricarboxylic acid (TCA) cycle in the mitochondria. Its importance for the oxidation of malate to oxaloacetate might partially explain the identification in the endosperm and embryo. The location of the above mentioned proteins in the Glycolysis/Gluconeogenesis pathway (KEGG pathway No. zma00010) might reflect the importance of differential accumulation of Glycolysis/Gluconeogenesis related proteins during endosperm and embryo development. And the fact that these proteins were detected from different hybrids suggested that the differential regulation of Glycolysis/Gluconeogenesis pathway depended on the demand level for various metabolic intermediates of Glycolysis/Gluconeogenesis pathway during the course of kernel development in different hybrids. Pyruvate produced during glycolysis is used to provide precursors for biosynthesis of certain amino acids and the reducing agent NADH for the oxidative phosphorylation pathway to produce ATP energy in the TCA cycle. The TCA cycle is important for energy generation through the oxidation of the three nutrients-carbohydrates, fats and proteins and for precursors production of certain amino acids and reducing agent NADH, but only malate dehydrogenase, an enzyme involved in the regeneration of oxaloacetate, was identified with continuous down-regulation in the endosperm of Pioneer 335 (spot no. 32) and in the embryo of Nongda 108 with peak at 22 DAP (spot no. 59) in this study (Table 2). The mitochondrial localization of the reactions might be one possible reason for the small proportion of identified enzymatic proteins in the TCA cycle [17], [18]. The extreme importance of TCA cycle for the whole development of each organ might be another reason for the small proportion of differentially expressed proteins in this reaction during the sampling stages. The down-regulation expression of this protein was consistent with a great decrease in TCA cycle activity at later stages of development in the embryo and endosperm of this study, the embryo of maize [31], and the seed of Arabidopsis [32].
Based on the sources of carbon and nitrogen, carbohydrates, fats and proteins are re-synthesized as storage substances both in sink tissues, such as the endosperm, and in viable tissues, such as the embryo, and grain weight is added. Proteins control almost every reaction within a cell, provide structure and serve as signals to other cells. Thus, protein synthesis, including initiation, elongation, and termination steps, should have an important regulatory function in seed development. During the development process, a translation initiation factor is required for the correct assembly of an initiation complex and an elongation factor is required for effective elongation. Translation initiation factor 5A (spot no. 36) and translation elongation factor (spot no. 28) were identified in the endosperm and displayed down- and up-regulation expression trends, respectively, over the course of the grain-filling process in the present study (Table 2, Table S1, Figure 4). Translation initiation factors and elongation factors are the actual cellular effectors of protein synthesis, and their abundance, in turn, depends on the balance between protein translation and protein degradation. Protein-binding proteins that function to ensure correct protein folding and degradation were also identified in the endosperm of Zhengdan 958 (spot no. 2), Nongda 108 (spot nos. 19 and 21), and Pioneer 335 (spot no. 29) and were continuously up-regulated or peaked at 25 DAP. However, none of the protein synthesis-related proteins was identified in the embryo of the tested hybrids. The difference might be attributable to the distinct developmental requirement of the endosperm and embryo during the linear phase of grain filling. This conclusion is supported by the accumulation of storage proteins in the embryo (spot nos. 52, 78, and 79). The storage proteins globulin 1 (spot no. 78) and 2 (spot nos. 52 and 79), which are destinations of remobilized nitrogen and sulfur [33], were up-regulated during embryo development, which suggested globulin modulated efficient partitioning of nutrients during the linear phase of grain filling. The level of storage proteins, a destination of accumulated nutrients, may contribute to grain-filling rate and grain weight. Proteins involved in lipid metabolism were represented by different protein spots (nos. 50, 65, and 74) but showed a similar up-regulated expression pattern (Table S1) in the three hybrids. Meanwhile, proteins involved in nitrogen metabolism were detected with different expression patterns in Nongda 108 (spot nos. 62 and 67) and Pioneer 335 (spot no. 75) (Table S1). Phytic acid is a phosphorylated derivative of myo-inositol, functions as the major storage form of phosphorus in plant seeds, and plays diverse roles in plants as a signal transduction molecule, osmoprotectant, and cell wall constituent. The inositol-3-phosphate synthase (spot no. 71) identified in the present study catalyzes the first step in de novo synthesis of myo-inositol. This protein showed continuous down-regulation, opposite to the expression of proteins involved in lipid metabolism and nitrogen metabolism, in Pioneer 335 (spot no. 71), which supported its negative regulation of grain-filling rate during the linear phase.
Specific Grain Filling Related Proteins in Embryo and Endosperm
Embryo development represents a critical stage in the sporophytic life cycle. The embryo becomes desiccated and dormant at late stages of development [34]–[36]. A group of hydrophilic proteins, known as LEA proteins, accumulate to high levels during the final stage of seed maturation [37]. These proteins are neither enzymes nor true storage proteins, but probably serve to protect embryos from desiccation stress and stabilize other proteins and membranes that protect the cell [38]–[40]. The genes that encode LEA proteins belong to a multigene family [40], [41]. The protein LEA 34 identified in the present study (spot nos. 51, 64, and 76) belongs to LEA group 5, of which transcripts accumulate during late stages of seed development and in response to drought [42], [43]. Generally, LEA proteins accumulate in embryos soon after storage protein accumulation and their synthesis in embryos coincides with desiccation of the seed. LEA 34 was identified in the embryo of the three hybrids and showed similar up-regulated expression during seed development. An additional embryo-abundant protein (spot no. 68) and embryo-specific protein (spot no. 80) identified in this study also showed up-regulated expression during seed development. These proteins may represent markers for embryo development stages [44].
The endosperm, which is an important nutritional supply tissue for embryo development during seed germination, should accumulate abundant starch during the linear phase of grain filling. However, none of the enzymes involved in the starch synthesis pathway was identified in the present study. Three possible reasons might explain this result: 1) differential expression of enzymes involved in starch synthesis was less than 2-fold during the linear phase of grain filling; or 2) starch synthesis, in which ADP glucose is converted to amylose or amylopectin by starch synthases and starch-branching enzyme, does not have a strong and direct regulatory relationship with sucrose transportation in the cytosol controlled by grain-filling rate during the linear phase; or 3) the quantitative sensitivity limitations of the technology used in this study might also be a possible reason. Sorbitol dehydrogenase is an enzyme in carbohydrate metabolism converting sorbitol, the sugar alcohol form of glucose, into fructose [45]. The spot no. 11 identified in this study is a sorbitol dehydrogenase homolog 1 protein, which was continuously up-regulated with the greatest difference in accumulation levels of all proteins and expressed a peak at 22 DAP at the transcriptional level for all hybrids tested here (Figure 5). Integration with genomic information verified that its coding gene was co-localized with a QTL for grain yield (qGY1b) on Chr. 1 [20]. In addition, the target sites of miRNA168 were also identified in this gene. These results implied that sorbitol dehydrogenase gene regulated by miRNA might be an important regulator of maize grain filling, especially for hybrid Zhengdan 958.
Materials and Methods
Plant Material
Three elite maize hybrids that show significant differences in grain-filling rate and grain-filling duration were used in this study. Nongda 108 (Xu 178×Huang C) was the most widely planted hybrid (in terms of planting area) between 2002 and 2004 in China, Zhengdan 958 (Zheng 58×Chang 7-2) was the most widely planted hybrid between 2005 and 2012 in China, and Pioneer 335 had the third-largest planting area between 2008 and 2012 in China. To obtain a similar flowering time in the field, Nongda 108 was planted on 1 May 2011, and Zhengdan 958 and Pioneer 335 were planted on 5 May 2011, which resulted in a flowering time between 4 and 9 July. For each genotype, two replication plots were planted on the research farm of Henan Agricultural University, Zhengzhou, China (113°42′ E, 34°48′ N). The plots consisted of eight rows of 4 m length, with inter-row spacing of 75 cm and within-row spacing of 25 cm. The hybrids were self-pollinated by hand on 6 July and ears were collected from each plot at 10, 17, 22, 25, 28, 33, 40, and 50 DAP. Only the middle rows were sampled to avoid edge effects. Kernels from the middle part of 90 ears were separated into two groups for each plot. Three one-hundred kernels per replication plot were oven-dried at 70°C for 24 h and the dry weight was recorded. The final dry weight of each replication plot was calculated by averaging the three 100-kernel dry weight values. The grain-filling rate was calculated by dividing the dry matter accumulated by the number of days between two sampling stages. The remainder of the kernels was dissected to remove the pericarp, and the endosperm and embryo were dissected for three biological replicates and stored at –80°C until use for protein and RNA extraction, respectively. Data analysis was performed using the SAS 8.0 statistical software package with the PROC MIXED procedure. On the basis of the change in grain-filling rate during the sampling stages (Figure 1A), endosperm and embryo samples harvested at 17, 22, 25, and 28 DAP for each genotype were subjected to further proteomic analysis.
Protein Extraction
The three biological replicates per tissue sample on each sampling date per hybrid were prepared. Approximately 1 g for each sample was ground in liquid nitrogen. Ten milliliters of chilled trichloroacetate buffer containing 10% trichloroacetate and 0.07% β-mercaptoethanol (in acetone) was added and mixed in tube. After 2 h incubation at –20°C, the extracts were centrifuged at 18,000×g for 30 min at 4°C and the supernatant was discarded. Ten milliliters of chilled 80% acetone containing 0.07% β-mercaptoethanol was added to each sample, mixed well and incubated for 1 h at –20°C before centrifuging at 18,000×g for 15 min at 4°C. This step was repeated twice and the pellet was freeze-dried under vacuum for embryo samples. The pellet for endosperm samples was added the following steps for the starchy characteristics. The pellet was resuspended in extraction buffer (30% sucrose, 0.1 M Tris-HCl [pH 8.0], 1% 1,4-dithio-DL-threitol [DTT], 100 mM KCl, 5 mM EDTA, 2 mM phenylmethylsulfonyl fluoride). The mixture was centrifuged at 18,000×g for 10 min at 4°C. The final combined supernatant was transferred to a new tube, and an equal volume of Tris-phenol was added and mixed. After centrifugation at 18,000×g for 30 min at 4°C, the phenol phase and a 5-fold volume of chilled 0.1 M ammonium acetate (in methanol) were incubated at –20°C for 30 min. The supernatant was discarded, and the protein pellet was collected by centrifugation at 12,000×g at 4°C for 30 min. After two thorough washes with cooled 0.1 M ammonium acetate (in methanol) and one wash with 80% cooled acetone, the resulting proteins were freeze-dried under vacuum. A solution containing 9 M urea, 2 M thiourea, 4% (w/v) CHAPS, and 1% DTT was added to the dried protein pellets at a ratio of 20 µl mg−1. After incubation for 1 h at 28°C with frequent shaking, the insoluble fraction was removed via centrifugation at 18,000×g for 40 min at 25°C. The concentration of the supernatant containing the soluble protein fraction was determined by the Bradford method with bovine serum albumin used as a standard [46].
Two-dimensional Gel Electrophoresis
The final extracted proteins were subjected immediately to 2-D gel electrophoresis or stored in aliquots at −80°C. Isoelectric focusing (IEF) of the proteins for each biological replicate of each genotype was performed with three technical replicates. For each genotype, samples of the four different developmental stages were analyzed simultaneously for convenience of comparison. For each replication, equal amounts of total protein extract (1 mg) were assayed using 24 cm immobilized dry strips (Immobiline DryStrips, Bio-Rad, Hercules, CA, USA) with a linear pH gradient of pH 4–7. After rehydration for 11 h at 50 V, IEF was performed at 18°C with the following conditions: 200 V for 1 h, 500 V for 1 h, 1000 V for 1 h, 5000 V for 1 h, and 10,000 V for 10 h. Current was limited to 0.05 mA per immobilized pH gradient gel strip. Equilibration of the strips was performed immediately with 10 mL of two types of sodium dodecyl sulfate (SDS) equilibration buffer for 15 min each. Buffer 1 contained 0.375 M Tris-HCl (pH 8.8), 6 M urea, 20% glycerol, 4% SDS, and 2% DTT, and buffer 2 contained 0.375 M Tris-HCl (pH 8.8), 6 M urea, 20% glycerol, 4% SDS, and 2.5% iodoacetamide. The immobilized pH gradient gel strips with the proteins were embedded in the top of a 2-D gel in 1% agarose solution after equilibration. Proteins were separated on the basis of their Mr on a 12% SDS-PAGE gel at 16°C under the following conditions: 1 W of constant power for 30 min, then increased to 6 W, which was maintained until the electrophoresis was finished. Preparative gels were stained with Coomassie Brilliant Blue R250 and imaged by a laser scanner (AlphaImager HP, Alpha Innotech, Santa Clara, CA, USA).
Identification of Protein Spots by Mass Spectrometry
Only protein spots that appeared consistently in nine technical replicates of the three biological replicates for each genotype were counted. Their expression was analyzed using Imagemaster 2-D Platinum software version 5.0 (Bio Rad). Only protein spots with more than 2-fold changes between the maximum expression level and the lowest expression level (absolutely intense difference) at different DAP stages for each genotype, as determined by ANOVA and Student’s t-test (P<0.01) after normalization, were excised manually and processed for MS. Spots from the 2-D gel were transferred individually to a microcentrifuge tube (Eppendorff) and destained with 25 mM NH4HCO3 in 50% (v/v) acetonitrile for three periods of 15 min each at room temperature. Gel spots were dried using a vacuum centrifuge at room temperature and incubated in 50 µl digestion solution, which consisted of 25 mM NH4HCO3 in 0.1% acetic acid and 12.5 ng mL−1 sequencing-grade trypsin at 37°C overnight (Promega, Madison, WI, USA). After digestion, a 0.3 µl sample was spotted onto a matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF) sample plate with the same volume of matrix (10 mg mL−1 α-cyano-4-hydroxycinnamic acid in 50% acetonitrile, 0.1% trifluoroacetic acid). Peptide mass spectra were obtained on a MALDI-TOF/TOF mass spectrometer (4700 Proteomics Analyzer, Applied Biosystems) in the positive ion reflector mode. The MS spectra were calibrated internally with the Mass Standards kit for the 4700 Proteomics Analyzer. Proteins were identified by automated peptide mass fingerprinting using the Global Proteome Server Explorer™ software version 3.5 (Applied Biosystems) against maize protein FASTA sequences (94,845 sequences) in the NCBI database. The unknown and unnamed protein spots were further searched in the KEGG database with blast function. Peak lists (S/N >10) were extracted from the raw data for data processing, and positive identifications were accepted up to the 95% confidence level using the following criteria: 1) precursor ion mass tolerance ±0.1 Da; 2) fragment ion mass tolerance ±0.3 Da; 3) maximum of one missed cleavage per peptide allowed; and 4) fixed modifications by cysteine carboxyamidomethylation and variable modifications by methionine oxidation.
Transcript Profiling
Endosperm tissue at 17, 22, 25, and 28 DAP were assayed for transcript profiling in accordance with the protein sampling date. Three samples were prepared on each sampling date per hybrid. Total RNA was isolated using the Plant Total RNA Extraction kit (Bioteke Corporation, China) in accordance with the manufacturer’s protocol. Absence of genomic DNA contamination was confirmed by PCR, using primers designed to flank the intron sequences of a control gene (Actin: 5′-TGGCATTGTCAACAACTGG-3′, 5′-CTCCTTGCTCATACGATCGG-3′). First-strand cDNA was synthesized using oligo (dT)15 as a primer and M-MLV Reverse Transcriptase (Promega, USA) in a 25-µl reaction volume. Quantitative real-time PCR was performed using SYBR® Premix Ex Taq™ (Takara, Dalian, China) on a DNA Engine Option®2 Continuous Fluorescence detection system (Bio-Rad, Waltham, USA). Melting curve analysis was performed to verify primer specificity. The relative quantity of the transcripts was calculated using the comparative threshold cycle method. An Actin control was used to normalize data across samples. Each sample was analyzed in triplicate for each gene and Actin, and the average value for each was used in analyses. The sequences of Actin primers for real-time PCR were 5′-CGATTGAGCATGGCATTGTCA-3′ and 5′-CCCACTAGCGTACAACGAA-3′.
Supporting Information
Figure S1.
2-D maps of the maize hybrids Nongda 108 and Pioneer 335. From left to right, 17, 22, 25, and 28 DAP are presented. Differentially expressed protein spots at any two stages with more than 2-fold changes in expression are indicated by arrows with numeral. (A) 2-D maps for endosperm of Nongda 108. (B) 2-D maps for endosperm of Pioneer 335. (C) 2-D maps for embryo of Nongda 108. (D) 2-D maps for embryo of Pioneer 335.
https://doi.org/10.1371/journal.pone.0059353.s001
(TIF)
Figure S2.
Changing pattern classification of differentially expressed proteins identified by MS at sampling stages. I, Proteins whose abundance increased linearly from 17 to 28 DAP; II, proteins whose abundance decreased linearly from 17 to 28 DAP; III, proteins up-regulated at 25 DAP; IV, proteins down-regulated at 25 DAP; V, proteins up-regulated at 22 DAP; VI, proteins down-regulated at 22 DAP; VII, proteins that changed irregularly.
https://doi.org/10.1371/journal.pone.0059353.s002
(TIF)
Table S1.
Expression levels of differential proteins identified by MS at sampling stages. I, Proteins whose abundance increased linearly from 17 to 28 DAP; II, proteins whose abundance decreased linearly from 17 to 28 DAP; III, proteins up-regulated at 25 DAP; IV, proteins down-regulated at 25 DAP; V, proteins up-regulated at 22 DAP; VI, proteins down-regulated at 22 DAP; VII, proteins that changed irregularly.
https://doi.org/10.1371/journal.pone.0059353.s003
(XLS)
Table S2.
Functional annotation of differentially expressed proteins identified by MS. Repeated peptides were omitted.
https://doi.org/10.1371/journal.pone.0059353.s004
(XLS)
Table S3.
Integration of QTL intervals and differentially expressed proteins identified by MS for grain-filling rate.
https://doi.org/10.1371/journal.pone.0059353.s005
(XLS)
Acknowledgments
Protein sequencing and MS identification of this work were performed by the Beijing Genomics Institute.
Author Contributions
Conceived and designed the experiments: JHT ZHL ZYF. Performed the experiments: XNJ ZYF WL. Analyzed the data: ZYF XNJ DD. Contributed reagents/materials/analysis tools: XNJ JHT. Wrote the paper: ZYF.
References
- 1. Yang W, Peng S, Dionisio-Sese ML, Laza RC, Visperas RM (2007) Grain filling duration a crucial determinant of genotypic variation of grain yield in fieldgrown tropical irrigated rice. Field Crop Res 105: 211–227.
- 2. Hanft JM, Jones RJ, Stumme AB (1986) Dry matter accumulation and carbohydrate concentration patterns of field grown and in vitro cultured maize kernels from the tip and middle ear positions. Crop Sci 26: 568–572.
- 3. Wang GL, Kang MS, Moreno O (1999) Genetic analyses of grain-filling rate and duration in maize. Field Crops Res 61: 211–222.
- 4.
Gasura E, Setimela P, Tarekegne A, Okori P, Edema R, et al.. (2010) Grain filling patterns of CIMMYT early maize inbred lines. Scond RUFORUM Biennial Meeting 20–24 september, Entebbe, Uganda Research Application Summary.
- 5. Saini HS, Westgate ME (2000) Reproductive development in grain crops during drought. Adv Agron 68: 59–95.
- 6. Zahedi M, Jenner CF (2003) Analysis of effects in wheat of high temperature on grain filling attributes estimated from mathematical models of grain filling. J Agr Sci 141: 203–212.
- 7. Yang JC, Zhang JH (2006) Grain filling of cereals under soil drying. New Phytol 169: 223–236.
- 8.
Kling JG, Edmeades G (1997) Morphology and growth of maize IITA/CIMMYT. Research Guide 9 Training programme. International Institute of Tropical Agriculture IITA Ibadan. Nigeria 1997. pp: 36.
- 9. Olsen OA (2004) Nuclear endosperm development in cereals and Arabidopsis thaliana. Plant Cell 16: S214–S227.
- 10. Sabelli PA, Larkins BA (2009) The development of endosperm in grasses. Plant Physiol 149: 14–26.
- 11. Mizukami Y, Fischer RL (2000) Plant organ size control AINTEGUMENTA regulates growth and cell numbers during organogenesis. Proc Natl Acad Sci USA 97: 942–947.
- 12. Mizukami YA (2001) A matter of size-developmental control of organ size in plants. Curr Opin Plant Biol 4: 533–539.
- 13. Hu Y, Xie Q, Chua NH (2003) The Arabidopsis auxin-inducible gene ARGOS controls lateral organ size. Plant Cell 15: 1951–1961.
- 14. Pedro RC, Roberto HR, Fernando HA (2007) Ecophysiological yield components of maize hybrids with contrasting maturity. American society of Agronomy 99: 1111–1118.
- 15. Guo M, Rupe MA, Dieter JA, Zou JJ, Spielbaver D, et al. (2010) Cell Number Regulator 1 affects plant and organ size in maizeimplications for crop yield enhancement and heterosis. Plant Cell 22: 1057–1073.
- 16. Jenner CF, Rathjen AJ (1972) Factors limiting the supply of sucrose to the developing wheat grain. Ann Bot 36: 729–741.
- 17. Méchin V, Balliau T, Château-Joubert S, Davanture M, Langella O, et al. (2004) A two-dimensional proteome map of maize endosperm. Phytochemistry 65: 1609–1618.
- 18. Prioul JL, Méchin V, Lessard P, Thévenot C, Grimmer M, et al. (2008) A joint transcriptomic proteomic and metabolic analysis of maize endosperm development and starch filling. Plant Biotechnolog Journal 6: 855–869.
- 19. Finnie C, Melchior S, Roepstorff P, Svensson B (2002) Proteome analysis of grain filling and seed maturation in barley. Plant Physiol 129: 1308–1319.
- 20. Liu ZH, Ji HQ, Cui ZT, Wu X, Duan LJ, et al. (2011) QTL detected for grain-filling rate in maize using a RIL population. Mol Breeding 27: 25–35.
- 21. Xu H, Zhang WP, Gao Y, Zhao Y, Guo L, et al. (2012) Proteomic analysis of embryo development in rice (Oryza sativa). Planta 235: 687–701.
- 22. Schnyder H (1993) The role of carbohydrate storage and redistribution in the source–sink relations of wheat and barley during grain filling. New Phytol 123: 233–245.
- 23. Takai T. Fukuta Y, Shiraiwa T, Horie T (2005) Time-related mapping of quantitative trait loci controlling grain-filling in rice (Oryza sativa L.) J Exp Bot. 56: 2107–2118.
- 24. Xu SB, Yu HT, Yan LF, Wang T (2010) Integrated proteomic and cytological study of rice endosperms at the storage phase. J Proteome Res 9: 4906–4918.
- 25. Alonso AP, Dale VL, Shachar-Hill Y (2010) Understanding fatty acid synthesis in developing maize embryos using flux analysis. Metab Eng 12: 488–497.
- 26. Gallardo K, Job C, Groot SP, Puype M, Demol H, et al. (2001) Proteomic analysis of Arabidopsis seed germination and priming. Plant Physiol 126: 835–848.
- 27. Skylas DJ, Copeland L, Rathmell WG, Wrigley CW (2001) The wheat-grain proteome as a basis for more efficient cultivar identification. Proteomics 1: 1542–1546.
- 28. Koller A, Washburn MP, Lange BM, Andon NL, Deciu C, et al. (2002) Proteomic survey of metabolic pathways in rice. Proc Natl Acad Sci USA 99: 11969–11974.
- 29. Ostergaard O, Melchior S, Roepstorff P, Svensson B (2002) Initial proteome analysis of mature barley seeds and malt. Proteomics 2: 733–739.
- 30. Hung G-C, Brown CR, Wolfe AB, Liu J, Chiang H-L (2004) Degradation of the gluconeogenic enzymes fructose-16-bisphosphatase and malate dehydrogenase is mediated by distinct proteolytic pathways and signaling events. J Biol Chem 279: 49138–49150.
- 31. Lee JM, Williams ME, Tingey SV, Rafalski JA (2002) DNA array profiling of gene expression changes during maize embryo development. Funct Integr Genomics 2: 13–27.
- 32. Fait A, Angelovici R, Less H, Ohad I, Urbanczyk-Wochniak E, et al. (2006) Arabidopsis seed development and germination is associated with temporally distinct metabolic switches Plant Physiol. 142: 839–854.
- 33. Kriz AL (1989) Characterization of embryo globulins encoded by maize Glb genes. Biochem Genet 27: 239–251.
- 34. Mayer U, Torres-Ruiz RA, Berleth T, Miseéra S, Jürgens G (1991) Mutations affecting body organization in the Arabidopsis embryo. Nature 353: 402–407.
- 35. Laux T, Würschum T, Breuninger H (2004) Genetic regulation of embryonic pattern formation. Plant Cell 16: S190–S202.
- 36. De Smet I, Lau S, Mayer U, Jürgens G (2010) Embryogenesis-the humble beginnings of plant life. Plant J 61: 959–970.
- 37. Tunnacliffe A, Wise M J (2007) The continuing conundrum of the LEA proteins. Naturwissenschaften 94: 791–812.
- 38. Cordova-Tellez L, Burris JS (2002) Embryo drying rates during the acquisition of desiccation tolerance in maize seed. Crop Sci 42: 1989–1995.
- 39. Wise MJ, Tunnacliffe A (2004) POPP the question what do LEA proteins do? Trends Plant Sci 9: 13–17.
- 40. Battaglia M, Olvera-Carrillo Y, Garciarrubio A, Campos F, Covarrubias AA (2008) The enigmatic LEA proteins and other hydrophilins. Plant Physiol 148: 6–24.
- 41. Bies-Ethève N, Gaubier-Comella P, Debures A, Lasserre E, Jobet E, et al. (2008) Inventory evolution and expression profiling diversity of the LEA late embryogenesis abundant protein gene family in Arabidopsis thaliana. Plant Mol Biol 67: 107–124.
- 42. Zegzouti H, Jones B, Marty C, Lelíèvre JM, Latchè A, et al. (1997) Er5 a tomato cDNA encoding an ethylene-responsive LEA-like protein characterization and expression in response to drought ABA and wounding. Plant Mol Biol 35: 847–854.
- 43. Kim HS, Lee JH, Kim JJ, Kim CH, Jun SS, et al. (2005) Molecular and functional characterization of CaLEA6 the gene for a hydrophobic LEA protein from Capsicum annuum. Gene 344: 115–123.
- 44. Joosen R, Cordewener J, Supena ED, Vorst O, Lammers M, et al. (2007) Combined transcriptome and proteome analysis identifies pathways and markers associated with the establishment of rapeseed microspore-derived embryo development. Plant Physiol 144: 155–172.
- 45. E1-Kabbani O, Darmanin C, Chung RP (2004) Sorbitol dehydrogenase: structure, function and ligand design. Curr Med Chem 11(4): 465–476.
- 46. Bradford MM (1976) A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72: 248–254.