Comparative proteomic analysis of the effect of temperature and fertilizer on gliadin and glutenin accumulation in the developing endosperm and flour from Triticum aestivum L. cv. Butte 86

Background Flour quality is largely determined by the gluten proteins, a complex mixture of proteins consisting of high molecular weight-glutenin subunits (HMW-GS), low molecular weight-glutenin subunits (LMW-GS), and α-, γ-, and ω-gliadins. Detailed proteomic analyses of the effects of fertilizer and high temperature on individual gliadin and glutenin protein levels are needed to determine how these environmental factors influence flour quality. Results Wheat plants (Triticum aestivum L. cv. Butte 86) were grown in greenhouses under moderate and high temperature regimens with and without post-anthesis fertilizer. Quantitative two-dimensional gel electrophoresis was used to construct accumulation profiles in developing endosperm for the entire complement of gluten proteins identified previously by tandem mass spectrometry. Amounts of individual gliadins and glutenins were also determined in flour produced under each of the regimens. Under all environmental regimens, most HMW-GS, LMW-GS, γ- and ω-gliadins accumulated rapidly during early stages of grain development and leveled off during middle stages of development. A subset of LMW-GS showed a second distinct profile, accumulating throughout development, while α-gliadins showed a variety of accumulation profiles. In flour, fourteen distinct gluten proteins responded similarly to fertilizer, high temperature, and high temperature plus fertilizer. The majority of HMW-GS and ω-gliadins and some α-gliadins increased while two LMW-GS and a minor γ-gliadin decreased. Fertilizer did not influence gluten protein accumulation under high temperature conditions. Additionally, the effects of fertilizer and high temperature were not additive; very few changes were observed when plants that received fertilizer were subjected to high temperature. Conclusions Although post-anthesis temperature and fertilizer have very different effects on grain development and yield, the two treatments elicit surprisingly similar effects on the accumulation of gluten proteins. The similarity of the responses to the different treatments is likely due to source-sink activities of nitrogen reserves in the wheat plant. Because each protein that showed a response in this study is linked to a gene sequence, the work sets the stage for transgenic studies that will better elucidate the roles of specific proteins in flour quality and in the response to the environment.


Background
Wheat flour is used worldwide as an ingredient in a variety of baked products due to its unique protein properties. When flour is mixed with water, the proteins interact to form gluten, a continuous protein network that is responsible for the viscoelastic properties and gas holding capacity of flour dough. Gluten is a complex mixture of several hundred proteins with molecular weights ranging from 30,000 to 88,000 Da. The gluten proteins are synthesized in the endosperm of the developing wheat grain and comprise 60-80% of flour protein. Based on extraction and solubility properties, the gluten proteins can be separated into two major groups: the alcohol-soluble gliadins and the alcoholinsoluble glutenins [1]. The gliadins contribute mainly to the viscosity and extensibility of wheat dough and the glutenins to elasticity and dough strength. The monomeric gliadins are classified as α-, γ-, and ω-gliadins based on their mobility in acid PAGE gels. The glutenins are classified as high molecular weight-glutenin subunits (HMW-GS) and low molecular weight-glutenin subunits (LMW-GS), which link together via intermolecular disulfide bonds to form large insoluble polymers. Variations in flour quality are related to genetic differences in gliadin and glutenin composition among wheat cultivars and to effects of environmental conditions on their relative proportions [2,3].
Since consistent flour quality is important to end users, an increased understanding of the effects of environment on individual gliadins and glutenins and their relationship to quality is needed. Due to the difficulties in separating individual gluten proteins by chromatographic methods [4], quantitative two-dimensional gel electrophoresis (2-DE) has become a method of choice for identifying proteins that respond to environmental cues. To date, a number of studies have used 2-DE to discover gluten protein responses to high temperature, fertilizer, or drought [5][6][7][8][9][10][11] and have reported changes in as few as three [8] to as many as 107 2-DE spots [11]. Mass spectrometry (MS) was used to identify some of the proteins that showed responses to environmental cues. However, in most studies, identifications were based on very few peptides and it was generally not possible to distinguish individual gluten proteins within the major classes.
Recently, we developed improved methods for tandem mass spectrometry (MS/MS) identification of the gluten proteins. Proteins in 2-DE spots were digested individually with chymotrypsin, thermolysin, and trypsin to increase sequence coverage. In addition, search strategies were improved and databases that included cultivarspecific gluten protein sequences were constructed for analysis of spectral data [12]. This methodology made it possible to develop a comprehensive proteome map of flour from the US wheat Butte 86 in which 233 individual 2-DE spots were assigned to specific gene sequences [13]. Gluten proteins identified included 5 HMW-GS, 22 LMW-GS, 13 γ-gliadins, 23 α-gliadins and 4 ω-gliadins, many of which had sequence coverage greater than 50%. Among the proteins identified in this study were α-, γand ω-gliadins that contain extra cysteine residues and may be linked into the glutenin polymer. In addition, αgliadins containing specific epitopes involved in celiac disease were distinguished from those that do not contain these sequences. The proteome map made possible a detailed analysis of the effects of fertilizer on gliadin and glutenin levels in flour [14]. One-hundred twentytwo gluten protein spots representing 19 gluten protein gene sequences responded significantly to fertilizer: five HMW-GS, six α-gliadins, and four ω-gliadins increased and three LMW-GS and one γ-gliadin decreased.
In this paper, we compare the effects of high temperature and fertilizer singly and in combination on the entire complement of gliadins and glutenins in developing endosperm and in flour from mature grain. The study extends earlier studies in which we analyzed the effects of high temperature and fertilizer on a limited number of gliadins and glutenins in developing endosperm by 2-DE and evaluated flour composition by RP-HPLC [5,6]. Because the Butte 86 proteome map links gluten proteins to gene sequences, the goal of this study was to identify genes encoding proteins that respond to fertilizer and high temperature that can be used in transgenic approaches to determine the roles of specific proteins in flour quality and in the response of the grain to environmental cues.

Results
High temperature and fertilizer alter fresh weight and gluten protein levels during grain development The plant growth regimens used and the comparisons made in this study are summarized in Table 1. All plants were grown at 24/17°C (day/night maximum temperatures) with fertilizer until anthesis, when pots were placed into  Figure 1A). Grain fresh weight was 79 mg at 37 dpa compared to 69 mg without fertilizer. The 37/28°C regimens shortened the duration of grain fill from 37 to 21 dpa and grain fresh weight was substantially less. Fresh weight/grain at 21 dpa under the high temperature regimens was similar with or without fertilizer. Although mature grain produced under the high temperature regimens was much smaller, the germination rate was 100% (data not shown). Total gluten protein amount/grain, determined by the method of Lowry et al. [15,16] for the KCl-insoluble protein fraction, was highest under the 24/17°C regimen with fertilizer ( Figure 1B). At 37 dpa under the 24/17°C regimens, gluten protein amount was somewhat higher with fertilizer than without. At 21 dpa under the 37/28°C regimens, gluten protein amount was similar with or without fertilizer, but substantially less than under the 24/17°C regimens. The high temperature regimens also increased the percentage of gluten protein per grain. With fertilizer, gluten protein was 1.8% of grain fresh weight under the 24/17°C regimen compared to 2.6% under the 37/28°C regimen. Without fertilizer, gluten protein was 1.6% of fresh weight under the 24/17°C regimen compared to 2.2% under the 37/28°C regimen. The increase in the proportion of the gluten proteins per grain under the 37/28°C regimens is due to a decrease in starch levels under the high temperature regimens [17,18].  Figure 2A). The three replicate gels for each time point clustered together, demonstrating excellent technical reproducibility of the 2-D gel patterns. PCA separated the gel patterns into four clusters that were aligned from the earliest to the latest time points along the x-axis. Cluster a contained the 8 dpa gel patterns for the 24/17°C regimens with or without fertilizer. Cluster b contained the 6 dpa gel patterns for the 37/28°C regimen with fertilizer. Clusters c and d contained the gel patterns for the remaining time points. Within cluster c, the gel patterns for the 24/17°C regimens with or without fertilizer at 13 dpa were related to those for the 37/28°C regimens at 8 dpa with fertilizer and 6 and 8 dpa without fertilizer. Within cluster d, the gel patterns for the 24/17°C regimen without fertilizer at 18, 23, 28, 33 and 37 dpa were located below the origin of the y-axis and the gel patterns for the 24/17°C regimen with fertilizer at 18, 23, 28, 33 and 37 dpa and the 37/28°C regimens with and without fertilizer at 13, 18 and 21 dpa were located above the origin. This distribution suggests that the 24/17°C regimen without fertilizer was different from the other regimens and that it was a suitable treatment for determining gluten protein responses to fertilizer and high temperature.  Table 1.
Comparison of the 2-D gel patterns of proteins in the developing endosperm provided a partial explanation for the temporal distribution of the clusters generated by PCA ( Figures 2B-E). Under the 24/17°C regimen at 8 dpa, the most abundant proteins were ATP synthase A (AtpA), calreticulin (CRT), 70 kD heat shock protein (Hsp70), protein disulfide isomerase (PDI), and stress responsive protein R40 (R40) ( Figure 2B). PDI and Hsp70 were identified previously [13,19] and ATP synthase A (AtpA), calreticulin (CRT), and stress responsive protein R40 (R40) were identified in this study (Additional file 2 contains the MS/MS identification data). In keeping with the onset of gluten protein biosynthesis, CRT, Hsp70, and PDI, which reside in the ER, facilitate the folding and maturation of proteins. CRT also modulates calcium homeostasis and Hsp70 assists in transport of precursor   Table 1. Boxes in E outline positions of the HMW-GS, the ω-gliadins, and the LMW-GS plus αand γ-gliadins. The proteome map and mass spectrometry data for identifications of individual proteins were reported previously [13]. Additional proteins identified in this study are indicated in panel B: AtpA, ATP synthase subunit alpha; CRT, calreticulin; Hsp70, 70 kD heat shock-induced protein; PDI, protein disulfide isomerase; R40, ABA-and stress-induced protein (Additional file 2 contains identification data). 13   proteins into organelles and targets damaged proteins for degradation [20]. PDI regulates the formation, reduction, and isomerization of disulfide bonds associated with protein folding and AtpA is a key enzyme in bioenergetics. The specific role of the R40 proteins remains unclear. It has been proposed that these ABA-and osmotic stress-responsive proteins [21,22] have both protein-and carbohydrate-binding domains enabling them to cross-link proteins and glycoconjugates [23]. The non-gluten proteins appear to decrease throughout development, but this is due to a concomitant large increase in gluten protein accumulation that masks the presence of most non-gluten proteins. In an earlier study [19], a fraction in which gluten proteins were depleted was essential for analysis of non-gluten protein accumulation during Butte 86 grain development. Under the 24/ 17°C regimens, gluten proteins were present at very low levels at 8 dpa ( Figure 2B), but comprised a major proportion of proteins from 13 to 37 dpa ( Figure 2D and E, Additional file 1). Under the 37/28°C regimens, gluten protein levels were more abundant at 6 dpa ( Figure 2C) than at 8 dpa under the 24/17°C regimens and comprised a major proportion of proteins from 8 to 21 dpa (Additional file 1).

High temperature and fertilizer alter individual gliadin and glutenin accumulation profiles during grain development
Using protein identifications of 2-DE spots determined by MS/MS in Dupont et al. [13], protein profiles were grouped into the HMW-GS, LMW-GS, and the α-, γ-, and ω-gliadins. These were examined by hierarchical cluster analysis to provide a comprehensive image of the accumulation patterns of the gluten proteins during endosperm development under the four environmental regimens ( Figure 3; original data in Additional file 3). The resulting display contains 8,208 data points: 24 experimental time points, 114 gluten protein spots, and spot volumes from triplicate 2-D gels. The spot volumes for the replicates clustered together at each time point, demonstrating excellent technical reproducibility of the 2-D gel patterns. The hierarchical cluster analysis ordered the accumulation profiles for each protein class by protein amount, which had a dynamic range of six orders of magnitude. Due to this wide dynamic range, a color scale of 200,000 to 7,000,000 was chosen to display the spot volume data. Depending on gluten protein class, 0.6-5.4% of the protein profiles with spot volumes less than 200,000 are off-scale and 7.7-16% with spot volumes greater 7,000,000 are off-scale. The dendrograms in Figure 3 show the relationships between the protein profiles within each of the gluten protein classes. The identities of the gluten proteins corresponding to the profiles depicted in Figure 3 are located in Table 2 and are listed in Additional file 3 in the order determined by the hierarchical cluster analysis. Profiles of gluten protein accumulation under the 24/17°C regimens appeared similar with or without fertilizer. Most HMW-GS, LMW-GS, and α-, γ-, and ω-gliadins increased rapidly between 8 and 13 dpa and then increased gradually from 13 to 37 dpa. The 37/28°C regimens compressed the developmental profiles; proteins increased rapidly between 6 and 8 dpa and attained maximal levels from 13 to 21 dpa.
Because individual profiles were difficult to distinguish in Figure 3, correlation analysis was used to investigate accumulation profiles in more detail. Protein profiles were grouped into HMW-GS, LMW-GS, and the α-, γ-, and ωgliadins, spot volumes log normalized, and accumulation patterns analyzed ( Figure 4). The profiles include the three replicate spot volumes for each time point and are color coded to match the spot numbers in Figure 3. This analysis showed, with the exception of a few outliers, that the profiles for the HMW-GS and the γand ω-gliadins formed single clusters under each of the environmental treatments while the LMW-GS formed two clusters. The α-gliadins exhibited greater variability and only five of the 21 profiles clustered together.
Under the 24/17°C regimen without fertilizer, the gluten proteins exhibited four developmental profiles (Figure 4, MN). The HMW-GS and α-gliadins reached maximum levels at 18 dpa and remained at these levels until 37 dpa.
The ω-gliadins attained maximum levels at 13-18 dpa and declined slightly by 37 dpa. The majority of LMW-GS ( Figure 4, LMW-GS A ) had profiles similar to the γgliadins and attained maximum levels at 28 dpa and declined slightly by 37 dpa. A small group of LMW-GS (spots 140, 145, 153, 170, 315) increased gradually throughout development and attained maximum levels at (See figure on previous page.) Figure 3 Accumulation patterns of gluten proteins during wheat endosperm development. Accumulation profiles were analyzed by hierarchical cluster analysis and visualized using the TIGR Multiexperiment Viewer. Each square represents the normalized volume of a single 2-D gel spot. The color scale indicates relative protein amount. The dendrograms show the relationships between the profiles within each of the gluten protein classes. Additional file 3 contains the data set of normalized spot volumes. Spot numbers are color coded to match the subset of proteins selected by correlation analysis shown in Figure 4. HMW-GS N56 and the LMW-GS N61 were identified in developing endosperm and not detected in flour; Additional file 1, Figure MN-13 shows the map location and Additional file 2 contains the mass spectrometry data for these proteins. MN, MF, HN, and HF are defined in Table 1.   [13]. 2 LMW-GS are classified as m-, s-or i-type based on their N-terminal amino acid sequence. 3 Denotes gliadins containing odd numbers of cysteines. 4 Denotes α-gliadins that do not contain Glia-α, Glia-α2, Glia-α-9 or Glia-α-20 celiac disease epitopes. 5 Spots containing same gene sequences. 6 Ave and SD refer to average and standard deviation of normalized volumes from all spots with the same identification. Additional file 4 contains the data set of normalized spot volumes. 7 MN, 24/17°C minus fertilizer. 8 MF, 24/17°C plus fertilizer. 9 HN, 37/28°C minus fertilizer. 10 HF, 37/28°C plus fertilizer. 11 Increases of 20% or more are indicated in bold type and decreases of 20% or more are indicated in bold italics type. † Proteins that responded significantly to fertilizer in Altenbach et al. [14]. ‡ Spot numbers of proteins comprising LMW-GS B , a group of LMW-GS sharing a unique developmental profile. 37 dpa (Figure 4, LMW-GS B , indicated by ‡ in the Spot Numbers column in Table 2). These spots were identified as the m-type LMW-GS Bu-1, Bu-7, Bu-11 and Bu-18 and the i-type LMW-GS Bu-4 [13]. Under the 24/17°C regimen with fertilizer, the gluten proteins exhibited only two developmental profiles (Figure 4, MF). The HMW-GS, LMW-GS A , and the α-, γ-, and ω-gliadins peaked at 28 and the LMW-GS B were maximal at 37 dpa. Compared to the 24/17°C regimen without fertilizer, the HMW-GS, αgliadins, and ω-gliadins were higher from 28-37 dpa, the LMW-GS A were lower from 28-37 dpa, the LMW-GS B were lower from 18-37 dpa, and the γ-gliadins were somewhat lower at 33 dpa (Figure 4, compare MN and MF). At 37 dpa, the most notable changes in gluten protein levels were increases in the αand ω-gliadins and decreases in the LMW-GS B .
Under the 37/28°C regimens, the gluten proteins exhibited three developmental profiles without and with post-anthesis fertilizer (Figure 4, HN and HF). The HMW-GS and LMW-GS A attained maximum levels at 13 dpa and remained at these levels until 21 dpa. The α-, γ-, and ω-gliadins also peaked at 13 dpa, but declined slightly by 21 dpa. The LMW-GS B increased gradually throughout development and attained maximum levels at 21 dpa. The most notable changes in gluten protein levels at 21 dpa under the two high temperature regimens were increases in the ω-gliadins.

High temperature and fertilizer alter levels of specific gliadins and glutenins in flour
Since flour quality depends on gluten protein composition, the effect of fertilizer and high temperature on  Figure 3 and may be used with Table 2 and Additional file 3 to determine the identities of the individual proteins within the gluten protein classes. MN, MF, HN, and HF are defined in Table 1.
gluten protein levels was examined in flour milled from mature grain harvested from plants used for the developmental studies (original data in Additional file 4). Because the same protein sequences were often associated with multiple 2-D spots, the volumes of spots containing proteins with the same identifications were summed for this analysis as was done by Altenbach et al. [14]. This resulted in 42 different gluten proteins ( Table 2). The percent change in volume relative to the 24/17°C regimen without fertilizer was calculated to determine specific gluten protein responses to fertilizer, high temperature, and high temperature plus fertilizer. In a previous study, 19 gluten proteins showed statistically significant changes of 20% or more in response to fertilizer under a 24/17°C regimen ( Table 2, indicated by † in the Predominant Protein column) [14]. In this study, it is notable that 14 proteins responded not only to fertilizer (MNvMF), but also to high temperature (MNvHN) and high temperature plus fertilizer (MNvHF). These included three HMW-GS (Ax2*, Dx5, By9), one m-type LMW-GS (Bu-11), five α-gliadins (Bu-5, Bu-11, Bu-12, Bu-14, Bu-BQ807130), one γ-gliadin (Bu-1) and one ω-5 gliadin and three ω-1,2 gliadins. The effects of fertilizer and high temperature were not additive. In fact, only three low abundance proteins differed 20% or more in amount when plants that received fertilizer were subjected to high temperature (MFvHF); an i-type LMW-GS (GenBank: AAB48469) increased and two α-gliadins (Bu-1, Bu-BQ806209) decreased. In addition, protein abundance changed little when plants subjected to high temperature received fertilizer (HNvHF). A comparison of percent change of the gluten protein groups ( Figure 5) in response to fertilizer, temperature, and high temperature plus fertilizer showed that the HMW-GS, α-gliadins, and ω-gliadins increased, the LMW-GS decreased, and the γ-gliadins showed little change in response to fertilizer, high temperature or fertilizer plus high temperature. The increases were greatest for the ω-gliadins and much less for the HMW-GS and α-gliadins. After accounting for gliadins containing extra cysteines (indicated in Table 2), ratios of gliadins to glutenins increased from 0.71 for the MN treatment to 0.86, 0.81 and 0.84 for the MF, HN and HF treatments, respectively. Ratios of HMW-GS to LMW-GS also increased from 0.44 for the MN treatment to 0.57, 0.58 and 0.59 for the MF, HN and HF treatments, respectively.

Discussion
Effect of fertilizer and high temperature on developmental profiles of gluten proteins  Table 2). These proteins had the same MS/ MS identities as other proteins in the LMW-GS A group, but had different pIs and M r s. Although MS/MS sequence coverage for these LMW-GS was quite high, generally about 50%, differences in the expression profiles suggest that very similar proteins in the LMW-GS A and LMW-GS B groups are encoded by different genes. Juhász et al. [24] recently found two expression profiles for LMW-GS genes. However, they reported that s-type and i-type LMW-GS had expression profiles similar to the LMW-GS B in the present study. Liu et al. [25] also reported multiple expression profiles for LMW-GS genes in a proteomic analysis of three Chinese wheat cultivars.
The developmental profiles demonstrate that gluten protein accumulation is a complex process that is subject to spatial and temporal regulation as well as environmental signaling. Individual proteins within each gluten protein class accumulated to different levels, suggesting that the corresponding genes had different basal levels of expression. Ultimately, the proportions of ω-gliadins and HMW-GS and a few α-gliadins and LMW-GS changed in flour in response to both postanthesis fertilizer and high temperature while those of most αand γ-gliadins and LMW-GS were unchanged. The data suggest that there are differences in the promoter sequences of the gluten protein genes, both between and within the major gene families. Promoter elements responsible for tissue specificity of the gluten proteins are the endosperm box, which is comprised of the prolamin box and a GCN4-like motif (GLM), and the AACA motif [26,27]. GLM, also known as the N motif, is thought to play a role in the response of the gluten protein genes to nitrogen [28,29]. Recently, Juhász et al. [24] compared regulatory elements in the 5 0 non-coding regions of LMW-GS genes available in public databases and developed a model for their transcriptional regulation based on an in silico analysis of expression profiles. Clearly, it will be important to confirm this model using experimental data and extend the model to include other classes of gluten protein genes. Comparisons of promoter regions from genes encoding gluten proteins characterized in this study should prove valuable for this type of analysis.

Effect of fertilizer and high temperature on gluten protein composition in flour
Dupont et al. [6] used RP-HPLC to quantify the effects of fertilizer and high temperature on the accumulation of Butte 86 flour proteins and reported changes in the levels of ω-gliadins and HMW-GS. However, changes in specific LMW-GS, αand γ-gliadins were difficult to detect because of the complexity of these groups of proteins. In this study, 2-DE analysis revealed a number of changes in distinct gluten proteins within these classes in response to fertilizer and temperature. Of the 42 uniquely identified gluten proteins, 16 increased with fertilizer under moderate temperatures (four HMW-GS, eight α-gliadins, and four ω-gliadins) and three decreased (two LMW-GS and one γ-gliadin). The response to the 37/28°C regimen without fertilizer (MNvHN) was comparable to that of 24/17°C regimen with fertilizer (MNvMF); thirteen of the same proteins increased (four HMW-GS, five α-gliadins, and four ω-gliadins) and two of the same proteins decreased (one LMW-GS and one γ-gliadin). Most of these proteins also responded to high temperature plus fertilizer (MNvHF). It is notable that the effects of fertilizer and high temperature were not additive and that very few significant changes were observed when plants receiving fertilizer were subjected to high temperature (MFvHF). Likewise, the application of fertilizer under a high temperature regimen had little effect on gluten protein accumulation (HNvHF). As noted previously [5,14,30,31], the response to high temperature and fertilizer is related to the proportion of the S-containing amino acids in the different gluten protein classes.
Generally, the proportions of low S-proteins (ω-gliadins) and low to medium S-proteins (HMW-GS and α-gliadins) increased and those of the high S-proteins (LMW-GS and γ-gliadins) remained the same or decreased.
In agreement with previous studies [6,14,32], the ratios of gliadins to glutenins increased in response to fertilizer and high temperature. The higher ratios reflect increases in the proportions of ω-gliadins and certain α-gliadins. The ratios of the HMW-GS to LMW-GS also increased in response to both fertilizer and high temperature, largely due to increases in HMW-GS and decreases in certain LMW-GS. Of the gliadins containing an odd number of cysteines and likely to function as LMW-GS, only the cys-type ω-gliadins showed responses to fertilizer or high temperature. Others have reported that the ratio of the HMW-GS to LMW-GS in a glutenin polymer fraction of SDS unextractable protein is significantly correlated with dough strength [33][34][35][36][37].
Because individual proteins in this analysis were associated with specific gene sequences, changes in the immunogenic potential of the flour that result from the growth conditions of the plant can now be assessed. It is notable that the Gli-B3 ω-gliadins (also referred to as ω-5 gliadins) are associated with the food allergy wheatdependent exercise-induced anaphylaxis. These proteins increased significantly in response to both fertilizer and temperature. Among the gluten proteins, the α-gliadins are considered most immunogenic in celiac disease. Several of the α-gliadins that increased in response to both fertilizer and temperature contain epitopes important in celiac disease (Bu-5, Bu-11, Bu-14) while others do not contain these sequences (Bu-12) ( Table 2) [13,38]. It is interesting that proteins that do not contain celiac epitopes comprised 31-32% of the total α-gliadins under all treatments, suggesting that neither fertilizer nor temperature influenced the immunogenic potential of the flour with respect to celiac disease. The accu-  Table 1.
mulation of other immunogenic proteins were not assessed in this study. However, Altenbach et al. [14] analyzed a fraction of total SDS-extractable protein and reported changes in the proportions of other proteins with fertilizer that are likely food allergens. These included increases in serpins and decreases in lipid transfer protein (LTP), chitinase, the α-amylase/trypsin inhibitor CM3 and globulin-2. Hurkman et al. [19] analyzed an albumin/globulin fraction and reported that the same proteins had opposite responses to high temperature when plants were supplied with fertilizer. They observed increases in LTP, chitinase and globulin-2 and decreases in the proportions of several serpins.

Source-sink activities and gluten protein responses to fertilizer and high temperature
The question arises as to why post-anthesis high temperature and the post-anthesis addition of fertilizer alter gluten protein composition in the same manner. The similarity in responses of the gliadins and glutenins to high temperature and fertilizer might be explained by sourcesink activities in the wheat plant [39][40][41][42][43][44][45]. All plants were grown with fertilizer and adequate water under moderate temperatures until anthesis so the amount of nitrogen and carbon reserves that could be mobilized to the developing grain for the synthesis of protein and starch were comparable at the beginning of the treatments. The addition of post-anthesis fertilizer enhances the remobilization of reserves during grain fill and the rate of protein accumulation in the grain, resulting in increases in flour protein percentage and specific changes in gliadin and glutenin composition in flour [5,6,14]. On the other hand, high temperature conditions elicit premature leaf senescence and shorten the period of grain fill [39]. Because remobilization of nitrogen reserves to the developing grain occurs over a significantly shorter period of time, the temperature treatment may essentially mimic high nitrogen conditions, resulting in changes in the same complement of gluten proteins. The lack of response to fertilizer under high temperature conditions suggests that the plants are unable to utilize fertilizer under the temperatures used in this study. This would also explain why post-anthesis fertilizer and temperature do not elicit an additive response.

Conclusions
This study provides insight into the manner in which specific growth conditions during wheat grain development influence the accumulation of the individual gluten proteins that determine flour quality. While many gluten proteins respond to fertilizer and high temperature, most of the changes are small in magnitude and surprisingly similar between the two treatments. It is notable that both fertilizer and high temperature increase the proportions of most HMW-GS and ω-gliadins and some of the same αgliadins while decreasing the proportions of several LMW-GS and a minor γ-gliadin. A critical question is how the observed changes in gluten protein composition influence flour quality. Many different proteins contribute to flour quality. While the roles of HMW-GS have received considerable attention, the functions of individual LMW-GS and gliadins in flour quality remain poorly defined. The identification of genes within these classes that respond to fertilizer or temperature now makes it possible to utilize transgenic approaches to determine the importance of specific proteins in flour quality and in the response of the grain to the environment. For example, RNA interference was used recently to create transgenic Butte 86 plants in which genes encoding the ω-5 gliadins were silenced [46]. Since these proteins show some of the largest responses to fertilizer and temperature, analyses of transgenic lines should elucidate how the proteome responds to fertilizer and high temperature in their absence and reveal the contribution of the ω-5 gliadins to flour quality.

Protein preparation
A KCl-insoluble fraction enriched in gliadins and glutenin subunits was prepared from endosperm collected at each time point as described previously [48]. Fifty mg of endosperm was ground to a fine powder in liquid nitrogen using a mortar and pestle and suspended in 200 μl of cold (4°C) KCl buffer containing protease inhibitors (1 tablet per 10 ml of Mini Complete Protease Inhibitor Cocktail, Roche Applied Science, Indianapolis, IN). The suspensions were incubated on ice for 5 min with intermittent mixing and centrifuged at 14,000 g for 15 min at 4°C. The pellet or KCl-insoluble fraction was suspended in 800 μl of SDS buffer (2% SDS, 10% glycerol, 50 mM DTT, 40 mM Tris-Cl, pH 6.8), incubated for 1 h at room temperature, and insoluble material removed by centrifugation at 14,000 g for 10 min at room temperature (Eppendorf 5415C; Brinkman Instruments, Inc., Westbury, NY). Proteins were then precipitated by addition of 4 vol of cold (−20°C) acetone and incubation overnight at −20°C to remove SDS, which interferes with protein determination and prevents separation of proteins by IEF. Samples were centrifuged at 14,000 g for 10 min at room temperature, the pellet was rinsed by pipetting cold acetone onto the pellet, centrifuging at 14,000 rpm for 10 min at room temperature, and pipetting the acetone off of the pellet. The pellet was dried at room temperature and solubilized in urea buffer (9 M urea, 4% NP-40, 1% DTT, and 2% ampholytes) to a final concentration of 3 μg protein/μl. For protein determinations, triplicate 5 μl samples were removed from the SDS extracts, precipitated with acetone, and, following addition of 0.1 N NaOH, protein amount assayed by the procedure of Lowry et al. [16] as described by Hurkman and Tanaka [49].

Statistical analysis of 2-D gel images
2-DE gels were digitized with a calibrated scanner (UMAX Powerlook III; Dallas, TX) at 300 dpi using the same settings for all gels. Flour proteins were identified previously by MS/MS [13] and additional proteins in the endosperm were identified by the same methods. Progenesis SameSpots Ver. 4.0 (Nonlinear Dynamics Limited, Newcastle upon Tyne, UK) was used to detect spots, match spot patterns, normalize and quantify spot volumes, and create protein accumulation profiles. The means and standard deviations for individual spot volumes were calculated and found to be quite similar among the replicate gels for each time point. All of the gluten proteins had Anova p-values less than 0.05 and Power greater than 0.8. Progenesis SameSpots was also used to analyze 2-D gel patterns by principal component analysis and match protein accumulation patterns by correlation analysis. The spot volume data were analyzed by hierarchical cluster analysis using the TIGR MultiExperimentViewer (MeV) v. 4.0 ([50, http://www.tm4.org/mev/]).