- 1Department of Plant Breeding, Institute for Sustainable Agriculture (IAS-CSIC), Córdoba, Spain
- 2Department of Crop and Forest Sciences, University of Lleida, Lleida, Spain
Temperature and nitrogen (N) availability are two important environmental factors that may produce important changes in grain composition during grain filling of bread wheat. In this study, four wheat lines with the down-regulation of gliadins by means of RNA interference (RNAi) have been characterized to determine the effect of thermal stress and N availability on grain weight and quality; with focus on gliadin and glutenin protein fractions. Grain weight was reduced with heat stress (HS) in all RNAi lines, whereas gliadin content was increased in the wild-types. With respect to gliadin content, RNAi lines responded to HS and N availability differently from their respective wild-types, except for ω-gliadin content, indicating a very clear stability of silencing under different environmental conditions. In a context of increased temperature and HS events, and in environments with different N availability, the RNAi lines with down-regulated gliadins seem well suited for the production of wheat grain with low gliadin content.
Introduction
Grain yield and quality are critical for wheat breeding and management. Both traits are determined during the grain-filling phase. Understanding the processes affecting grain weight and quality during grain filling is important for improving breeding and management strategies. Grain filling is commonly partitioned into three phases: the lag phase, the effective grain filling period, and the maturation drying phase (Egli, 1998). The lag phase is a period of active cell division, characterized by a rapid increase in water content with almost no dry matter accumulation. Grain dry weight then increases rapidly during the effective grain filling period until the maximum dry weight is attained, after which it remains approximately stable while the grain dries. During the effective grain filling period, starch and proteins are deposited in the endosperm (Jenner et al., 1991). It has been shown that the rate of their deposition is controlled by the source-sink balance (Fischer et al., 1977). Approximately 80% of total proteins in wheat grain are gluten (also termed prolamins) whereas the remaining 20% is composed of non-gluten proteins (NGPs) – mainly albumins and globulins (van den Broeck et al., 2009; Wen et al., 2012). Wheat gluten is able to form a network responsible for the viscoelastic properties of wheat flour since it allows the retention of carbon dioxide released during fermentation (Shewry, 2009). Gluten proteins can be further divided into two fractions: glutenins and gliadins (Lafiandra and Kasarda, 1985; Shewry, 2019). The glutenins form polymers linked by inter-chain disulfide bonds, they are insoluble in alcohol solutions, and can be divided according to their mobility in sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) into high molecular weight (HMW) and low molecular weight (LMW) glutenins. The gliadins are monomeric proteins, soluble in alcohol, and divided into three groups according to their mobility by electrophoresis in polyacrylamide gels at acidic pH (A-PAGE); ω, α and γ-gliadins. The glutenins are responsible for wheat dough elasticity and strength, while gliadins are important for viscosity and extensibility (Shewry and Halford, 2002; Shewry et al., 2003). Gluten proteins, particularly the gliadin fraction, are the primary factors responsible for triggering celiac disease (CD), since they contain the most immunogenic CD epitopes (Arentz-Hansen et al., 2000, 2002; Shan et al., 2002; Molberg et al., 2003; Ludvigsson et al., 2013; Box 1).
BOX 1. Celiac disease (CD) is a chronic enteropathy that results from the ingestion of gluten proteins present in wheat, and other similar proteins in barley and rye (Trier, 1998; Sollid, 2002). After ingestion of gluten, lesions form in the small intestine, characterized by flattening of the microvilli, hyperplasia of crypt cells, and infiltration of leukocytes (Sollid, 2002). As a result, symptoms such as diarrhea and malabsorption of food appear among others, since the spectrum of symptoms can be very broad. The immune response is triggered by the activation of CD4 T cells when they recognize the gluten peptides presented by serotypes HLA-DQ2 and HLA-DQ8. The presence of gluten peptides can be detected by the activity of the tissue transglutaminase 2 enzyme from the intestinal mucosa (Sollid, 2002; Sollid et al., 2012; Gayathri and Rashmi, 2014). CD is present throughout the world and the prevalence in the United States is around 1%, as in Europe, with the highest estimates in Finland and Sweden, and the lowest in Germany (Catassi et al., 2014). Gluten is present in many food products as the main element or as an additive. So far, the only possible treatment for CD is to follow a gluten-free diet for life (Sollid, 2002). The increase in the incidence of the disease was associated with the duration of exposure to gluten (Ventura et al., 1999), which increases the need to obtain wheat with a reduced content of proteins immunogenic for celiac sufferers.
In addition to CD, there are other pathologies related to wheat: (i) allergies as wheat-dependent exercise-induced anaphylaxis (WDEIA) – induced by the ω-5 gliadins and the HMW- (Morita et al., 2007; Morita et al., 2009), or baker’s asthma associated with non-specific lipid transfer proteins (Brant, 2007; Palacin et al., 2007); and (ii) non-celiac wheat sensitivity (NCWS) (Gibson et al., 2017), with an estimated prevalence ranging from 0.6 to 13% of global population (Aziz et al., 2016). Most of the allergens and proteins related to wheat pathologies have been mapped to the bread wheat Chinese Spring reference genome (RefSeq v1.0, International Wheat Genome Sequencing Consortium) (Appels et al., 2018) contributing to the knowledge of these diseases (Juhász et al., 2018). Moreover, there is a broad study on wheat allergens and CD peptides that allows their identification and composition for diagnostic assays by liquid chromatography-tandem mass spectrometry (Lexhaller et al., 2019).
RNA interference-based (RNAi) techniques are ideal for the down-regulation of specific protein fractions related to CD. Using this technology, γ-gliadins were silenced in two lines of bread wheat, providing reductions of up to 80% in this gliadin fraction (Gil-Humanes et al., 2008). Subsequently, the same workers used chimeric interference RNAs capable of silencing the genes from all the three groups of ω, γ and α-gliadins, to obtain several lines of two wheat genotypes with major reductions (in some cases up to 90%) in total gliadin content (Gil-Humanes et al., 2010). The crossing of the silenced lines with commercial varieties of wheat has allowed the obtention of new lines that effectively express the fragment of silencing in different genetic backgrounds, both to silence γ-gliadins (Gil-Humanes et al., 2012) and total gliadin fractions.
Climate model projections suggest that higher temperatures and heat stress events will become commonplace in most regions where grain crops are produced (Meehl and Tebaldi, 2004). Deleterious effects of high temperature on crop yield and quality are well documented in the literature (e.g. Barnabás et al., 2008). It is also well known that temperate species, such as wheat, maximize their vegetative growth during the period of colder temperatures, and the grain develops as temperature rises. Much works have studied the effects of very high temperature – moderate and short periods – during grain filling in wheat (Wardlaw et al., 2002; Farooq et al., 2011; Nuttall et al., 2018). Typically, high temperature during the grain filling linear phase results in the reduction of grain weight, mainly due to the decrease of soluble starch synthase activity under heat stress (Hawker and Jenner, 1993), reducing starch accumulation (e.g. Bhullar and Jenner, 1986). Simultaneously, heat stress (HS) not only increases protein percentage (Stone, 2001; Wardlaw et al., 2002), but also affects the synthesis of the different prolamin fractions and their ratios, which are responsible of the bread quality (Blumenthal et al., 1993; Stone, 2001).
Nitrogen (N) fertilization is one of the most common management practices used by farmers to improve yields. Consequently, there have been many studies analyzing yield (Foulkes et al., 1998; Barraclough et al., 2014) and protein content (Fischer et al., 1993) in response to soil N availability in wheat. In addition, several studies reported the effects of N on the types of proteins being synthesized during grain filling (Pechanek et al., 1997; Daniel and Triboi, 2000; Johansson et al., 2013), indicating that the synthesis of proteins in cereals is clearly influenced by temperature and N condition under which grain filling proceeds. However, genotypic variability can be found in all these responses (Saint Pierre et al., 2008; Elbashir et al., 2017); and even in the response to interactions between heat and N (Elía et al., 2018; Slafer and Savin, 2018). Therefore, determining the effects of heat and N availability on wheat genotypes with contrasting protein composition are of particular interest in the understanding grain protein distribution and its influence on grain weight and quality.
RNAi lines with low gluten content were subjected to various N and sulfur treatments to study the stability of the gliadin silencing under different fertilization conditions (García-Molina and Barro, 2017). In relation to N, this study showed that the RNAi lines had consistently lower gliadin levels than the wild-type across different N-fertilization regimes, but also that the level of gliadins in RNAi lines was sometimes significantly increased when N availability increased. In that study, N was applied when it would strongly affect grain number and yield and, therefore, may have diluted the availability of N-compounds during grain filling (i.e. more N available for absorption had to be shared between much higher grain numbers). As late N fertilization can be used to maintain green tissues during grain filling and to increase overall N content of the grains (Blandino et al., 2015), it may be relevant to explore whether the response of the RNAi lines would be even more marked. As mentioned before, grain filling is significantly affected by HS which would also favor the synthesis of proteins compared with that of starch (Barlow et al., 2015). It would thus be of interest to determine whether the synthesis of proteins in general, and gliadins in particular, in these RNAi lines is affected by the combinations of high temperatures and N availability during grain filling.
In this context, the objective of the present work was to determine the effects of contrasting temperature and N availability conditions on the silencing of gliadins. Thus, grain weight, total protein content and gluten protein distribution were studied in a set of RNAi lines and their respective wild-types. The aim was to determine to what degree the silencing of the synthesis of gliadins depends on environmental conditions, which is important for progressing in the development of low-gliadin wheat varieties suitable for new dietary approaches for gluten-related disorders.
Materials and Methods
Plant Material, Chamber Experiment and Treatments
Six lines of bread wheat were used: BW208, D770, D793, Gazul, J631, and M959. BW208 is a line from CIMMYT and Gazul is a commercial variety, and both were used as wild-types. D770 and D793 are lines derived from BW208 with RNA interference (RNAi) silencing of all gliadin fractions (Gil-Humanes et al., 2010). J631 and M959 are derived from crossing the Gazul genotype and line D770. Lines J631 and M959 were backcrossed at least four times with Gazul, always selecting the silencing character and the high and low molecular weight glutenins of Gazul, so that both silenced lines maintain the glutenin profile of this parent line.
We carried out a chamber experiment involving six wheat lines (two wild-type cultivars and four RNAi lines), two temperature treatments (control and heat stress, HS) during the linear phase of grain filling period, and two nitrogen (N) availabilities with three replicates, each replicate was composed of 6 plants (all in all there were 18 plants per genotype × N × HS; i.e. 216 plants per chamber). Two seeds were sown in pots (270 cm3) filled with a mixture of 30% peat and 70% soil. After emergence, one plant was left in each pot.
Plants were grown outdoors until heading when all pots were placed in a growth chamber set at 20/15°C. Different temperature treatments were imposed from 10 days after anthesis (DAA) during 10 consecutive days (Supplementary Figure 1). The control was set at 25/18°C in a chamber and the HS treatment to 40/18°C in another chamber. Minimum and maximum temperatures of 18 and 25 or 40°C were maintained for 8 and 6 h, respectively (Supplementary Figure 1). After the 10 days of treatment, temperatures were set to 25/18°C until maturity.
Chambers were set under long-day conditions (16 h). Pots inside the chambers were rearranged approximately once a week to minimize the effects of possible differences in microenvironment at different positions within each chamber. Pots were watered regularly to avoid water stress. N (9 mg per pot) was applied as urea diluted in all pots at flag leaf appearance (DC 4.5, Zadoks et al., 1974). At heading (DC 5.9), half of the plants received another dose of N (21 mg per pot).
Grain Weight and Total Protein Determination
At maturity, 18 plants per treatment were sampled. Mature grain weight was determined as the average of all grains from the main spikes harvested. Total grain protein content was determined by Dumas methodology (Dumas, 1831).
Prolamins Quantification by RP-HPLC
For gliadin and glutenin extraction, two grains from three different plants of each line and treatment were weighed and ground using a ball mill, and sequentially extracted following a previous protocol (Pistón et al., 2011) adapted to small samples. Briefly, gliadins were extracted stepwise three times with up to 400 μL of 60% (v/v) ethanol. Samples were centrifuged, and the supernatants collected, mixed together and filtered. The insoluble pellet was re-suspended in 50% (v/v) 1-propanol, 2 M urea, 0.05 M Tris–HCl (pH 7.5), and 2% (w/v) DTT for glutenin extraction, incubated for 30 min at 60°C and centrifuged stepwise three times. For each sample, the three collected supernatants were mixed together and filtered. The protein extracts were used for gliadin and glutenin quantification by Reverse-Phase High-Performance Liquid Chromatography (RP-HPLC, 1200 Series Quaternary LC System liquid chromatography from Agilent Technologies) with a DAD UV-V detector at 210 nm. A 25 cm long column LiChrospher® 100 RP8 (5 μm) (Merck) was used at 50°C and a sample volume equivalent to 2 mg of flour was injected. The flow rate was 0.5 mL⋅min–1. Mobile phase consisted in a mixture of Acetonitrile (ACN 0.1% TFA) and 0.1% aqueous TFA in a linear gradient (0 min 26% ACN, 60 min 54% ACN). The absolute amount of protein was calculated using bovine serum albumin protein as standard (BSA; BSA ≥ 98%, fraction V. Sigma-Aldrich, St Louis, MO, United States cat. no. A3294) (Supplementary Figure 2). The intervals of retention time used for the separation of prolamin fractions peaks are indicated in Supplementary Figure 3 according to Wieser et al. (1998). The integration of the peaks was performed automatically by RP-HPLC software with minor modifications if necessary.
Non-gluten Proteins (NGPs) Determination
The NGPs content was calculated by the difference between the total protein and prolamin content (glutenins and gliadins) for each line. The total protein content (μg/mg) was calculated from the percentage of N obtained by Dumas using the wheat conversion factor (5.83) (Merrill and Watt, 1973).
Data Processing
The retention time (min) and area (mAU) output of the RP-HPLC software was imported into a house developed software made in Python v2.71 to obtain the average values from the transformed technical repeats using the following formulas, that processes the hundreds of output files in a single run. The output of the software is a file with Microsoft Excel format.
The integration of the profiles, to obtain the area of each peak, and the subsequent transformation using the formula described, allow estimating the amount of protein for the samples. The arithmetic mean of the three technical repetitions was used for the variance analysis.
Statistical Analysis
To determine the effect of HS and N multifactorial univariates ANOVA were performed. Two variants of this model were tested: in the first, genotype, temperature, nitrogen and their interactions were independent variables, while grain weight and protein fractions were dependent variables. It was used to determine the general effect of the treatments on all genotypes. The second, has the same factors and variables, but it was performed for wild-types and RNAi lines separately to determine the effect of the treatments on each of these groups. Principal Components Analysis, PCA, was carried out with grain weight, total gliadin and its fractions, total glutenin and its fractions, and total prolamin as variables to evaluate their contribution to the model variance. The software used for the statistical analysis was R v 3.5.1 (R Core Team, 2018).
Results
Heat Stress and Nitrogen Treatment Effects on Grain Weight and Total Protein
Grain weight was significantly decreased in both wild-types and RNAi lines by heat stress (HS) (Figure 1A, Table 1, and Supplementary Table S1). Additional applications of nitrogen (N) had no significant effect on the RNA interference (RNAi) lines or wild-types (Figure 1B and Table 1). No significant differences were found in total grain protein content among all genotypes (Table 1). HS for a short period did not significantly modify the total protein content for both wild-types and RNAi lines (Figure 1C), but the late application of N (N1) resulted in a significant increase of the total protein for both the wild-types and the RNAi lines (Figure 1D).
Figure 1. Grain weight and total protein content for wild-types and RNAi lines under control and heat stress temperature treatments (A,C) and nitrogen availability (B,D). N0: no N application after heading, N1: N application after heading; control: 25/18°C during whole grain filling period, Heat stress (HS): 40/18°C for 10 days during grain filling period. The black line represents the median value. ∗ above the bars indicates significant difference (∗P ≤ 0.05; ∗∗∗P ≤ 0.001) between treatments according to the variance analysis.
Table 1. Significance of the variance of effects of genotype (6 lines: BW208, Gazul, D770, D793, J631 and M959), temperature (2 levels: Control temperature and heat stress), nitrogen (2 levels: No and N1) and their interactions for each variable studied (grain weight, protein fractions and protein ratios).
Heat Stress and Nitrogen Treatment Effects on Gliadins and Glutenins
We confirmed that total gliadin content was significantly higher in the wild-types than in the RNAi lines (Figure 2A and Table 1). Among RNAi lines, D793 had lower content of gliadins than that of the rest of RNAi lines (Supplementary Table S1). Both HS and the late application of N (N1) resulted in a significant increase in the total gliadin content for the wild-type lines, whereas no significant variation was observed for the RNAi lines (Figures 2B,C).
Figure 2. Content of total gliadin and total glutenin content for different genotypes with comparisons between RNAi lines and their wild-type by Dunnett’s test (A), wild-types and RNAi lines under different temperature treatments (B) and nitrogen availability (C). N0: no N application after heading, N1: application after heading; control: 25/18°C during whole grain filling period, Heat stress HS: 40/18°C for 10 days during grain filling period. The black line represents the median value. ∗ above the bars indicates significant difference (∗P ≤ 0.05; ∗∗P ≤ 0.01) between treatments according to the variance analysis.
The degree of silencing of ω-gliadins was lower than that of the rest of the gliadin fractions (Supplementary Figure 4A and Supplementary Table S1). A significant increase in the content of ω-gliadins due to supplementary N was found on both wild-types and RNAi lines, while HS only had an effect on the wild-types (Supplementary Figures 4B,C and Table 1). An overall effect of N level on the α-gliadin content (Supplementary Figure 4C and Table 1), as well as on the total gliadin content (Figure 2C), was observed in the wild-types, but the RNAi lines did not show this effect. HS treatment had no effect on the α-gliadin content for wild-types and RNAi lines (Supplementary Figure 4B). Conversely, γ-gliadin content was not affected by N availability, but a decrease in the amount of this fraction under HS was seen in wild-types, but not in the RNAi lines (Supplementary Figure 4B and Supplementary Table S1). Overall, the RNAi line D793 showed the highest reduction in α- and γ-gliadins (Supplementary Figure 4A and Supplementary Table S1).
The content of glutenins at grain maturity (Supplementary Table S1) was not significant higher for most RNAi genotypes than that of the wild-type lines (Figure 2A). The content of HMW was not statistically different between wild-types and RNAi lines (Supplementary Figure 5A and Table 1). HS had no effect on the HMW fraction, and for N treatment, RNAi lines tend to have a higher amount of HMW. In contrast, the LMW fraction was decreased in most of the RNAi lines in comparison to that of the wild-types (Supplementary Figure 5A, Table 1, and Supplementary Table S1). In addition, LMW content was affected by HS and N availability; in the wild-types LMW proteins decreased under HS; in RNAi lines LMW fractions increased at higher N availability (Supplementary Figures 5B,C and Supplementary Table S1).
The total gliadin/total protein ratio (GLI/TP) confirms that RNAi lines have lower gliadin content than the wild-types, particularly line D793 whatever experimental conditions. However, a lower total glutenin/total protein ratio (GLU/TP), was seen in wild-types than in RNAi lines without additional N supply, and in BW208 under HS with N1 than in RNAi lines (Figure 3 and Table 1).
Figure 3. Total gliadin/total grain protein ratio (GLI/TP), total glutenin/total grain protein ratio (GLU/TP) and NGPs/total grain protein ratio (NGPs/TP) between treatments and genotypes. The ratios were obtained with mean values for protein content. Low N: no N application after heading, High N: N application after heading; control: 25/18°C during whole grain filling period, Heat stress HS: 40/18°C for 10 days during grain filling period.
HS and N availability modified the GLI/TP ratio in the wild-type lines in different ways: it increased with N1 and HS in BW208 and decreased in Gazul (Figure 3). In contrast, minor effects were found in the GLI/TP ratio for the RNAi lines (Figure 3). Regarding the GLU/TP ratio, it was also modified by HS and N availability; it was strongly decreased in Gazul wild-type and RNAi derived lines under HS treatment and N0, while only a minor effect on lines with BW208 background was observed under those conditions; GLI/TP ratio was increased in BW208 RNAi lines under HS and N1, with only but minor effects in Gazul RNAi lines (Figure 3).
A Principal Component Analysis (PCA) was carried out, considering the effect of genotypes, temperature and N availability treatments on the variation of the protein fractions and grain weight (Figure 4A). Among the gliadin fractions, ω-gliadin and α-gliadin fractions contribute less and more, respectively, to the variance of the model. HMW proteins were the glutenin fraction that contributed most (Figure 4A). Gliadin and glutenin contents varied in perpendicular directions, indicating an independent behavior of both families of proteins (Figure 4A). The wild-types are separated from the RNAi lines (Figure 4B). The ellipses of 95% confidence level of each genotype indicated that there was a strong association between the variation of the glutenins and the silenced genotypes, and on the other hand, between the variation in the prolamins and the wild-types. D770 and J631 tend toward the direction of variation of the LMW proteins, and D793 and M959 toward that of the total glutenins (Figure 4). It is interesting to note that the variation in grain weight was independent of variations in the different protein fractions.
Figure 4. Principal Components Analysis (PCA). Effect of genotypes, temperature and N availability treatments on the variation of the protein fractions and grain weight (A). The high values in the color scale indicates a high contribution to the PCA. The direction and the size of the vectors indicate the relationship between all variables and their contribution to each axis. (B) Individuals are represented on the PCA axes with the 95% confidence ellipses showed for each genotype. The largest point for each genotype indicates the intersection of ellipse axes.
Discussion
Grain weight, total grain protein, and prolamin content under control temperature and N0 (control conditions) resulted in values similar to those previously reported for gliadin down-regulated lines (Gil-Humanes et al., 2010; Pistón et al., 2013). The decrease of the gliadin/total protein (GLI/TP) ratio in the RNA interference (RNAi) lines could be explained by the increase of the NGPs, as the glutenin/total protein (GLU/TP) ratio was higher in RNAi lines. In previous studies with these and other RNAi lines, protein compensation was observed (Altenbach et al., 2014; García-Molina and Barro, 2017) with increments in non-gluten proteins (NGPs) such as serpins, triticins and globulins (Gil-Humanes et al., 2011; Barro et al., 2016).
Brief heat stress (HS) events during the grain filling period generally result in a decrease in grain weight (Bhullar and Jenner, 1986; Savin et al., 1999). The range of variation depends on the genotype, timing and severity of HS (Balla et al., 2019). In the present study, we found a reduction in grain weight of 30% on average under HS and different availability of nitrogen (N). Grain weight was reduced by 35% for lines with BW208 genetic background, and about 24% for Gazul and its RNAi derived lines. This reduction could be mainly due to an extreme dependence on the temperature of starch synthesis, and an irreversible effect of HS on starch production after only a few days or even a few hours per day under control (Daniel and Triboi, 2000; Triboï et al., 2003; Spiertz et al., 2006; DuPont et al., 2006b; Liu et al., 2011; Hurkman et al., 2013) or field conditions (Savin et al., 1996; García et al., 2016; Elía et al., 2018). The percentage of grain protein generally increases under moderately high and very high temperatures (Stone, 2001; Wardlaw et al., 2002), either by a reduction of starch greater than the accumulation of protein, or by a reduction in starch without no change in protein accumulation. However, this response may not always occur for different genetic and environmental backgrounds (Graybosch et al., 1995). Interestingly, in the present study, HS resulted in a higher proportion of gliadin fraction in the wild-types. Other studies reported no effect of HS on total prolamins or even found a decrease with HS (DuPont et al., 2006b; Hurkman et al., 2013). However, Daniel and Triboi (2000) studied each fraction of gliadins and found that the proportion of ω- and α-gliadins increased with HS while γ-gliadins decreased, as found in the wild-types in the present study with the exception of α-gliadins. Also, in this work, grain protein content was increased under HS when post-anthesis N availability was higher (N1), whereas the total prolamin content of the wild-types differed in the response to HS. The RNAi lines, regardless of their genetic background, did not respond to temperature treatments for the total gliadin content, and for the gliadin fractions. This is an indication that these lines have robust gliadin silencing, independently of the temperature environment. Several authors have described a slight increase in LMW and HMW glutenin fractions with HS (DuPont et al., 2006a, b). Evidence has been also reported that the effect of HS is to cause a reduction in the size of glutenin polymers (e.g. Naeem et al., 2012). Nevertheless, in the present work, total glutenin and their fractions content were not modified significantly with HS in any of the genotypes, except for LMW in wild-types, but there is a non-significant decrease in all fractions in wild-types and RNAi lines.
Under higher N availability (N1), both wild-type and RNAi lines increase total grain protein, as previously described for other wheat genotypes (Daniel and Triboi, 2000; Triboï et al., 2003). Moreover, when increasing and splitting N doses, grain protein concentration increases and protein composition changes by increasing glutenin fractions (Xue et al., 2016). The response to N1 in wild-type lines in greater and in RNAi lines in lesser proportion, also confirmed that gliadin to glutenin ratio increase with N1 (Triboï et al., 2003). The content of ω-gliadins under N1 was increased, in comparison to N0, in the wild-types and RNAi lines, whereas the α-gliadin fraction was increased only in the wild-types as previously reported in D793 and other RNAi lines (Gil-Humanes et al., 2010; García-Molina and Barro, 2017). Total gliadin content did not increase in RNAi lines when additional N was supplied (N1), and this increase of ω-gliadins under N1 in RNAi lines has to be considered in further designing silencing constructs to improve their effectiveness since the ω-5 gliadins are related to wheat-dependent exercise-induced anaphylaxis (WDEIA) (Inomata, 2009; Morita et al., 2009) and ω-1,2 gliadins to CD (Tye-Din et al., 2010). However, the α-gliadins are reported as the major immunogenic complex in wheat, they contain three major celiac disease (CD) immunogenic peptides (Ozuna et al., 2015), and active peptides from this gliadin fraction were responsible for most of the immune response in patients with CD after eating wheat (Tye-Din et al., 2010). The α-gliadins were strongly reduced in the RNAi lines, and this was not affected either by HS or N application. Although some authors have indicated that LMW proteins decrease with high N availability at moderate temperatures (DuPont et al., 2006b; Hurkman et al., 2013), we found that LMW content increased in RNAi lines when additional N was supplied (N1). In contrast, the response in the HMW fraction and total glutenin content under N1, was not statistically significant in any of the genotypes.
Conclusion
Wheat grain proteins are important for the breadmaking quality of wheat, but they are also related to human pathologies as celiac disease (CD) and other gluten intolerances. RNA interference (RNAi) technology has provided wheat lines with all the gliadin fractions strongly down-regulated. Heat stress (HS) and nitrogen (N) availability could affect the synthesis and deposition of proteins during grain filling. Wild-types and RNAi lines studied in this work responded similarly for total grain protein and the content of ω-gliadins to additional N supply, as well as for the grain weight under HS. While the wild-types increase their total gliadin content under HS or high N availability the RNAi lines did not. Interestingly, the α-gliadin content, the most CD immunogenic fraction, is unaffected in the RNAi lines under additional N supply, but it was increased in wild-types. Therefore, under the specific scenario of brief events of temperature increase or additional application of N, studied in this work, the RNAi lines demonstrated a high stability of down-regulation of gliadins. However, further evaluations under field conditions will be necessary to confirm that the silencing of gliadin fractions in RNAi lines can be maintained under different abiotic stress environments.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
FB and RS designed the work. MM-S, MG, and RS carried out the work. All authors wrote the manuscript and approved the manuscript.
Funding
The Spanish Ministry of Economy and Competitiveness (Project AGL2016-80566-P) and the European Regional Development Fund (FEDER) supported this research.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The technical assistance of Ana García and the comments on the revised version of the manuscript by Gustavo A. Slafer are gratefully acknowledged. We also thank Dr. Paul A. Lazzeri (Agrasys SL) for his help in English review of the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00314/full#supplementary-material
Footnotes
References
Altenbach, S. B., Tanaka, C. K., and Seabourn, B. W. (2014). Silencing of omega-5 gliadins in transgenic wheat eliminates a major source of environmental variability and improves dough mixing properties of flour. BMC Plant Biol. 14:393. doi: 10.1186/s12870-014-0393-1
Appels, R., Eversole, K., Feuillet, C., Keller, B., Rogers, J., Stein, N., et al. (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. doi: 10.1126/science.aar7191
Arentz-Hansen, E., McAdam, S., Molberg, Ø., Kristiansen, C., and Sollid, L. (2000). Production of a panel of recombinant gliadins for the characterisation of T cell reactivity in coeliac disease. Gut 46, 46–51. doi: 10.1136/gut.46.1.46
Arentz-Hansen, H., Mcadam, S. N., Molberg, Ø., Fleckenstein, B., Lundin, K. E., Jørgensen, T. J., et al. (2002). Celiac lesion T cells recognize epitopes that cluster in regions of gliadins rich in proline residues. Gastroenterology 123, 803–809. doi: 10.1053/gast.2002.35381
Aziz, I., Dwivedi, K., and Sanders, D. S. (2016). From coeliac disease to noncoeliac gluten sensitivity; should everyone be gluten free? Curr. Opin. Gastroenterol. 32, 120–127. doi: 10.1097/MOG.0000000000000248
Balla, K., Karsai, I., Bónis, P., Kiss, T., Berki, Z., Horváth, Á., et al. (2019). Heat stress responses in a large set of winter wheat cultivars (Triticum aestivum L.) depend on the timing and duration of stress. PLoS One 14:e0222639. doi: 10.1371/journal.pone.0222639
Barlow, K., Christy, B., O’leary, G., Riffkin, P., and Nuttall, J. (2015). Simulating the impact of extreme heat and frost events on wheat crop production: a review. Field Crops Res. 171, 109–119. doi: 10.1016/j.fcr.2014.11.010
Barnabás, B., Jäger, K., and Fehér, A. (2008). The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ. 31, 11–38.
Barraclough, P. B., López-Bellido, R., and Hawkesford, M. J. (2014). Genetic variation in the uptake, partitioning and remobilisation of nitrogen during grain-filling in wheat. Field Crops Res. 156, 242–248. doi: 10.1016/j.fcr.2013.10.004
Barro, F., Iehisa, J. C., Giménez, M. J., García-Molina, M. D., Ozuna, C. V., Comino, I., et al. (2016). Targeting of prolamins by RNAi in bread wheat: effectiveness of seven silencing-fragment combinations for obtaining lines devoid of coeliac disease epitopes from highly immunogenic gliadins. Plant Biotechnol. J. 14, 986–996. doi: 10.1111/pbi.12455
Bhullar, S., and Jenner, C. (1986). Effects of temperature on the conversion of sucrose to starch in the developing wheat endosperm. Funct. Plant Biol. 13, 605–615.
Blandino, M., Vaccino, P., and Reyneri, A. (2015). Late-season nitrogen increases improver common and durum wheat quality. Agron. J. 107, 680–690. doi: 10.2134/agronj14.0405
Blumenthal, C., Barlow, E., and Wrigley, C. (1993). Growth environment and wheat quality: the effect of heat stress on dough properties and gluten proteins. J. Cereal Sci. 18, 3–21. doi: 10.1006/jcrs.1993.1030
Catassi, C., Gatti, S., and Fasano, A. (2014). The new epidemiology of celiac disease. J. Pediatr. Gastroenterol. Nutr. 59, S7–S9. doi: 10.1097/01.mpg.0000450393.23156.59
Daniel, C., and Triboi, E. (2000). Effects of temperature and nitrogen nutrition on the grain composition of winter wheat: effects on gliadin content and composition. J. Cereal Sci. 32, 45–56. doi: 10.1006/jcrs.2000.0313
DuPont, F. M., Hurkman, W. J., Vensel, W. H., Chan, R., Lopez, R., Tanaka, C., et al. (2006a). Differential accumulation of sulfur-rich and sulfur-poor wheat flour proteins is affected by temperature and mineral nutrition during grain development. J. Cereal Sci. 44, 101–112. doi: 10.1016/j.jcs.2006.04.003
DuPont, F. M., Hurkman, W. J., Vensel, W. H., Tanaka, C., Kothari, K. M., Chung, O. K., et al. (2006b). Protein accumulation and composition in wheat grains: effects of mineral nutrients and high temperature. Eur. J. Agron. 25, 96–107. doi: 10.1016/j.eja.2006.04.003
Egli, D. B. (1998). Seed Biology and the Yield of Grain Crops. New York, NY: CAB International, 178.
Elbashir, A. A., Gorafi, Y. S., Tahir, I. S., Elhashimi, A. M., Abdalla, M. G., and Tsujimoto, H. (2017). Genetic variation in heat tolerance-related traits in a population of wheat multiple synthetic derivatives. Breed. Sci. 67, 483–492. doi: 10.1270/jsbbs.17048
Elía, M., Slafer, G. A., and Savin, R. (2018). Yield and grain weight responses to post-anthesis increases in maximum temperature under field grown wheat as modified by nitrogen supply. Field Crops Res. 221, 228–237. doi: 10.1016/j.fcr.2018.02.030
Farooq, M., Bramley, H., Palta, J. A., and Siddique, K. H. (2011). Heat stress in wheat during reproductive and grain-filling phases. Crit. Rev. Plant Sci. 30, 491–507. doi: 10.1080/07352689.2011.615687
Fischer, R., Aguilar, I., and Laing, D. (1977). Post-anthesis sink size in a high-yielding dwarf wheat: yield response to grain number. Aust. J. Agric. Res. 28, 165–175.
Fischer, R. A., Howe, G. N., and Ibrahim, Z. (1993). Irrigated spring wheat and timing and amount of nitrogen fertilizer. I. Grain yield and protein content. Field Crops Res. 33, 37–56. doi: 10.1016/0378-4290(93)90093-3
Foulkes, M. J., Sylvester-Bradley, R., and Scott, R. K. (1998). Evidence for differences between winter wheat cultivars in acquisition of soil mineral nitrogen and uptake and utilization of applied fertilizer nitrogen. J. Agric. Sci. 130, 29–44. doi: 10.1017/s0021859697005029
García, G. A., Serrago, R. A., Dreccer, M. F., and Miralles, D. J. (2016). Post-anthesis warm nights reduce grain weight in field-grown wheat and barley. Field Crops Res. 195, 50–59. doi: 10.1016/j.fcr.2016.06.002
García-Molina, M. D., and Barro, F. (2017). Characterization of changes in gluten proteins in low-gliadin transgenic wheat lines in response to application of different nitrogen regimes. Front. Plant Sci. 8:257. doi: 10.3389/fpls.2017.00257
Gayathri, D., and Rashmi, B. (2014). Development of Celiac disease; pathogenesis and strategies to control: a molecular approach. J. Nutr. Food Sci. 4, 1–10.
Gibson, P. R., Skodje, G. I., and Lundin, K. E. (2017). Non-coeliac gluten sensitivity. J. Gastroenterol. Hepatol. 32, 86–89. doi: 10.1111/jgh.13705
Gil-Humanes, J., Pistón, F., Giménez, M. J., Martín, A., and Barro, F. (2012). The introgression of RNAi silencing of γ-gliadins into commercial lines of bread wheat changes the mixing and technological properties of the dough. PLoS One 7:e45937. doi: 10.1371/journal.pone.0045937
Gil-Humanes, J., Pistón, F., Hernando, A., Alvarez, J. B., Shewry, P. R., and Barro, F. (2008). Silencing of γ-gliadins by RNA interference (RNAi) in bread wheat. J. Cereal Sci. 48, 565–568.
Gil-Humanes, J., Pistón, F., Shewry, P. R., Tosi, P., and Barro, F. (2011). Suppression of gliadins results in altered protein body morphology in wheat. J. Exp. Bot. 62, 4203–4213. doi: 10.1093/jxb/err119
Gil-Humanes, J., Pistón, F., Tollefsen, S., Sollid, L. M., and Barro, F. (2010). Effective shutdown in the expression of celiac disease-related wheat gliadin T-cell epitopes by RNA interference. Proc. Natl. Acad. Sci. U.S.A. 107, 17023–17028. doi: 10.1073/pnas.1007773107
Graybosch, R. A., Peterson, C., Baenziger, P. S., and Shelton, D. (1995). Environmental modification of hard red winter wheat flour protein composition. J. Cereal Sci. 22, 45–51. doi: 10.1016/s0733-5210(05)80006-2
Hawker, J., and Jenner, C. (1993). High temperature affects the activity of enzymes in the committed pathway of starch synthesis in developing wheat endosperm. Funct. Plant Biol. 20, 197–209.
Hurkman, W. J., Tanaka, C. K., Vensel, W. H., Thilmony, R., and Altenbach, S. B. (2013). 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. Proteome Sci. 11:8. doi: 10.1186/1477-5956-11-8
Inomata, N. (2009). Wheat allergy. Curr. Opin. Allergy Clin. Immunol. 9, 238–243. doi: 10.1097/ACI.0b013e32832aa5bc
Jenner, C., Ugalde, T., and Aspinall, D. (1991). The physiology of starch and protein deposition in the endosperm of wheat. Funct. Plant Biol. 18, 211–226.
Johansson, E., Malik, A. H., Hussain, A., Rasheed, F., Newson, W. R., Plivelic, T., et al. (2013). Wheat gluten polymer structures: the impact of genotype, environment, and processing on their functionality in various applications. Cereal Chem. 90, 367–376. doi: 10.1094/cchem-08-12-0105-fi
Juhász, A., Belova, T., Florides, C. G., Maulis, C., Fischer, I., Gell, G., et al. (2018). Genome mapping of seed-borne allergens and immunoresponsive proteins in wheat. Sci. Adv. 4:eaar8602. doi: 10.1126/sciadv.aar8602
Lafiandra, D., and Kasarda, D. D. (1985). One-and two-dimensional (two-pH) polyacrylamide gel electrophoresis in a single gel: separation of wheat proteins. Cereal Chem. 62, 314–319.
Lexhaller, B., Colgrave, M. L., and Scherf, K. A. (2019). Characterization and relative quantitation of wheat, rye, and barley gluten protein types by liquid chromatography–tandem mass spectrometry. Front. Plant Sci. 10:1530. doi: 10.3389/fpls.2019.01530
Liu, P., Guo, W., Jiang, Z., Pu, H., Feng, C., Zhu, X., et al. (2011). Effects of high temperature after anthesis on starch granules in grains of wheat (Triticum aestivum L.). J. Agric. Sci. 149, 159–169. doi: 10.1017/s0021859610001024
Ludvigsson, J. F., Leffler, D. A., Bai, J. C., Biagi, F., Fasano, A., Green, P. H., et al. (2013). The Oslo definitions for coeliac disease and related terms. Gut 62, 43–52. doi: 10.1136/gutjnl-2011-301346
Meehl, G. A., and Tebaldi, C. (2004). More intense, more frequent, and longer lasting heat waves in the 21st century. Science 305, 994–997. doi: 10.1126/science.1098704
Merrill, A., and Watt, B. (1973). Energy Value of Foods: Basis and Derivation. Washington, DC: U.S. Department of Agriculture.
Molberg, Ø., Solheim Flaete, N., Jensen, T., Lundin, K. E., Arentz-Hansen, H., Anderson, O. D., et al. (2003). Intestinal T-cell responses to high-molecular-weight glutenins in celiac disease. Gastroenterology 125, 337–344. doi: 10.1016/s0016-5085(03)00890-4
Morita, E., Kunie, K., and Matsuo, H. (2007). Food-dependent exercise-induced anaphylaxis. J. Dermatol. Sci. 47, 109–117.
Morita, E., Matsuo, H., Chinuki, Y., Takahashi, H., and Dahlström, J. (2009). Food-dependent exercise-induced anaphylaxis—importance of omega-5 gliadin and HMW-glutenin as causative antigens for wheat-dependent exercise-induced anaphylaxis. Allergol. Int. 58, 493–498. doi: 10.2332/allergolint.09-RAI-0125
Naeem, H., Paulon, D., Irmak, S., and MacRitchie, F. (2012). Developmental and environmental effects on the assembly of glutenin polymers and the impact on grain quality of wheat. J. Cereal Sci. 56, 51–57. doi: 10.1016/j.jcs.2011.10.014
Nuttall, J. G., Barlow, K. M., Delahunty, A. J., Christy, B. P., and O’Leary, G. J. (2018). Acute high temperature response in wheat. Agron. J. 110, 1296–1308. doi: 10.2134/agronj2017.07.0392
Ozuna, C. V., Iehisa, J. C., Giménez, M. J., Alvarez, J. B., Sousa, C., and Barro, F. (2015). Diversification of the celiac disease α-gliadin complex in wheat: a 33-mer peptide with six overlapping epitopes, evolved following polyploidization. Plant J. 82, 794–805. doi: 10.1111/tpj.12851
Palacin, A., Quirce, S., Armentia, A., Fernández-Nieto, M., Pacios, L. F., Asensio, T., et al. (2007). Wheat lipid transfer protein is a major allergen associated with baker’s asthma. J. Allergy Clin. Immunol. 120, 1132–1138. doi: 10.1016/j.jaci.2007.07.008
Pechanek, U., Karger, A., Gröger, S., Charvat, B., Schöggl, G., and Lelley, T. (1997). Effect of nitrogen fertilization on quantity of flour protein components, dough properties, and breadmaking quality of wheat. Cereal Chem. 74, 800–805. doi: 10.1094/cchem.1997.74.6.800
Pistón, F., Gil-Humanes, J., and Barro, F. (2013). Integration of promoters, inverted repeat sequences and proteomic data into a model for high silencing efficiency of coeliac disease related gliadins in bread wheat. BMC Plant Biol. 13:136. doi: 10.1186/1471-2229-13-136
Pistón, F., Gil-Humanes, J., Rodríguez-Quijano, M., and Barro, F. (2011). Down-regulating γ-gliadins in bread wheat leads to non-specific increases in other gluten proteins and has no major effect on dough gluten strength. PLoS One 6:e24754. doi: 10.1371/journal.pone.0024754
R Core Team (2018). A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Saint Pierre, C., Peterson, C. J., Ross, A. S., Ohm, J. B., Verhoeven, M. C., Larson, M., et al. (2008). Winter wheat genotypes under different levels of nitrogen and water stress: changes in grain protein composition. J. Cereal Sci. 47, 407–416. doi: 10.1016/j.jcs.2007.05.007
Savin, R., Calderini, D. F., Slafer, G. A., and Abeledo, L. G. (1999). Final grain weight in wheat as affected by short periods of high temperature during pre-and post-anthesis under field conditions. Funct. Plant Biol. 26, 453–458.
Savin, R., Stone, P. J., and Nicolas, M. E. (1996). Responses of grain growth and malting quality of barley to short periods of high temperature in field studies using portable chambers. Aust. J. Agric. Res. 47, 465–477.
Shan, L., Molberg, Ø., Parrot, I., Hausch, F., Filiz, F., Gray, G. M., et al. (2002). Structural basis for gluten intolerance in celiac sprue. Science 297, 2275–2279. doi: 10.1126/science.1074129
Shewry, P. R. (2019). What is gluten-why is it special? Front. Nutr. 6:101. doi: 10.3389/fnut.2019.00101
Shewry, P. R., and Halford, N. G. (2002). Cereal seed storage proteins: structures, properties and role in grain utilization. J. Exp. Bot. 53, 947–958. doi: 10.1093/jexbot/53.370.947
Shewry, P. R., Halford, N. G., and Lafiandra, D. (2003). Genetics of wheat gluten proteins. Adv. Genet. 49, 111–184. doi: 10.1016/s0065-2660(03)01003-4
Slafer, G. A., and Savin, R. (2018). Can N management affect the magnitude of yield loss due to heat waves in wheat and maize? Curr. Opin. Plant Biol. 45, 276–283. doi: 10.1016/j.pbi.2018.07.009
Sollid, L. M. (2002). Coeliac disease: dissecting a complex inflammatory disorder. Nat. Rev. Immunol. 2, 647–655. doi: 10.1038/nri885
Sollid, L. M., Qiao, S.-W., Anderson, R. P., Gianfrani, C., and Koning, F. (2012). Nomenclature and listing of celiac disease relevant gluten T-cell epitopes restricted by HLA-DQ molecules. Immunogenetics 64, 455–460. doi: 10.1007/s00251-012-0599-z
Spiertz, J., Hamer, R., Xu, H., Primo-Martin, C., Don, C., and Van Der Putten, P. (2006). Heat stress in wheat (Triticum aestivum L.): effects on grain growth and quality traits. Eur. J. Agron. 25, 89–95. doi: 10.1016/j.eja.2006.04.012
Stone, P. (2001). “The effects of heat stress on cereal yield and quality,” in Crop Responses and Adaptations to Temperature Stress, ed. A. S. Basra (New York, NY: Food Products Press), 243–291.
Triboï, E., Martre, P., and Triboï-Blondel, A. M. (2003). Environmentally-induced changes in protein composition in developing grains of wheat are related to changes in total protein content. J. Exp. Bot. 54, 1731–1742. doi: 10.1093/jxb/erg183
Trier, J. S. (1998). Diagnosis of celiac sprue. Gastroenterology 115, 211–216. doi: 10.1016/s0016-5085(98)70383-x
Tye-Din, J. A., Stewart, J. A., Dromey, J. A., Beissbarth, T., van Heel, D. A., Tatham, A., et al. (2010). Comprehensive, quantitative mapping of T cell epitopes in gluten in celiac disease. Sci. Transl. Med. 2:41ra51. doi: 10.1126/scitranslmed.3001012
van den Broeck, H. C., America, A. H., Smulders, M. J., Bosch, D., Hamer, R. J., Gilissen, L. J., et al. (2009). A modified extraction protocol enables detection and quantification of celiac disease-related gluten proteins from wheat. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 877, 975–982. doi: 10.1016/j.jchromb.2009.02.035
Ventura, A., Magazzù, G., and Greco, L. (1999). Duration of exposure to gluten and risk for autoimmune disorders in patients with celiac disease. Gastroenterology 117, 297–303. doi: 10.1053/gast.1999.0029900297
Wardlaw, I. F., Blumenthal, C., Larroque, O., and Wrigley, C. W. (2002). Contrasting effects of chronic heat stress and heat shock on kernel weight and flour quality in wheat. Funct. Plant Biol. 29, 25–34.
Wen, S., Wen, N., Pang, J., Langen, G., Brew-Appiah, R. A., Mejias, J. H., et al. (2012). Structural genes of wheat and barley 5-methylcytosine DNA glycosylases and their potential applications for human health. Proc. Natl. Acad. Sci. U.S.A. 109, 20543–20548. doi: 10.1073/pnas.1217927109
Wieser, H., Antes, S., and Seilmeier, W. (1998). Quantitative determination of gluten protein types in wheat flour by reversed-phase high-performance liquid chromatography. Cereal Chem. 75, 644–650. doi: 10.1371/journal.pone.0192804
Xue, C., Rossmann, A., Schuster, R., Koehler, P., and Mühling, K.-H. (2016). Split nitrogen application improves wheat baking quality by influencing protein composition rather than concentration. Front. Plant Sci. 7:738. doi: 10.3389/fpls.2016.00738
Keywords: gluten proteins, heat stress, transgenic lines, celiac disease, Triticum aestivum
Citation: Marín-Sanz M, Giménez MJ, Barro F and Savin R (2020) Prolamin Content and Grain Weight in RNAi Silenced Wheat Lines Under Different Conditions of Temperature and Nitrogen Availability. Front. Plant Sci. 11:314. doi: 10.3389/fpls.2020.00314
Received: 20 November 2019; Accepted: 03 March 2020;
Published: 20 March 2020.
Edited by:
Katharina Anne Scherf, Karlsruhe Institute of Technology (KIT), GermanyReviewed by:
Ludovit Skultety, Institute of Virology (SAS), SlovakiaRobert Winkler, National Polytechnic Institute, Mexico
Eva Johansson, Swedish University of Agricultural Sciences, Sweden
Copyright © 2020 Marín-Sanz, Giménez, Barro and Savin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roxana Savin, c2F2aW5AcHZjZi51ZGwuY2F0