# N Partitioning between Gluten Fractions in a Set of Italian Durum Wheat Cultivars: Role of the Grain N Content

Marina Mefleh, Rosella Motzo, Marie-Françoise Samson, Marie-Hélène Morel, Francesco Giunta

<figure><img src="image_1.png" /></figure>

## To cite this version:

Marina Mefleh, Rosella Motzo, Marie-Françoise Samson, Marie-Hélène Morel, Francesco Giunta. N Partitioning between Gluten Fractions in a Set of Italian Durum Wheat Cultivars: Role of the Grain N Content. Foods, 2020, 9 (11), pp.1684. 10.3390/foods9111684. hal-03127951

HAL Id: hal-03127951

https://hal.inrae.fr/hal-03127951

Submitted on 1 Feb 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# N Partitioning between Gluten Fractions in a Set of Italian Durum Wheat Cultivars: Role of the Grain N Content

Marina Mefleh ¹D, Rosella Motzo ¹,*, Marie-Françoise Samson ²D, Marie-Hélène Morel ² and Francesco Giunta ¹

¹ Dipartimento di Agraria, Sezione Agronomia, Coltivazioni erbacee e Genetica, Università degli Studi di Sassari, Via De Nicola, 07100 Sassari, Italy; mmefleh@uniss.it (M.M.); giunta@uniss.it (F.G.)

² IATE, Univ Montpellier, CIRAD, INRAE, Institut Agro, 34060 Montpellier, France; marie-francoise.samson@inrae.fr (M.-F.S.); marie-helene.morel@inrae.fr (M.-H.M.)

* Correspondence: motzo@uniss.it

Received: 28 September 2020; Accepted: 13 November 2020; Published: 18 November 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Grain protein content constitutes a key quality trait for durum wheat end-products and may also impact grain protein composition. A total of sixteen durum wheat cultivars were analyzed in a field trial during two seasons at two nitrogen (N) levels to evaluate whether and to what extent the variation in total grain N was associated with variation in the quantity of the various protein fractions and grain quality parameters. Genotypic variation in grain N content correlated with the variation in the content of all three protein fractions, although the strength of the correlation with gliadin and albumin-globulin was higher than that with glutenins. Genotypic variation in gliadin and glutenin content was more tightly correlated with the variation in the sulfur (S)-rich protein groups than with the S-poor protein groups and subunits. The variation in the percentage of unextractable polymeric proteins (UPP%) among genotypes was independent of their glutenin allelic composition. The significant genotypic differences in UPP% and in the ratios between protein groups and subunits were not influenced by the corresponding variation in grain N content. The final grain N content can only account for part of the variation in quality parameters and in the partitioning of total grain N between protein fractions since genotypic differences other than grain N content also contribute to these variations.

**Keywords:** durum wheat; grain N partitioning; grain N content; gliadin/glutenin; S-rich/S-poor

## 1. Introduction

Durum wheat (*Triticum turgidum* L. spp. *durum* Desf.) is mainly cultivated in the Mediterranean region, where it is primarily used as the preferred ingredient to make pasta, as well as a number of traditional and specialty breads [1]. A grain with a minimum protein content of 12–15% is required to manufacture pasta as it produces a uniform particle sized semolina, an elastic resilient non-sticky dough, and an al dente texture of the cooked pasta [2,3]. However, when the overall protein content is low to intermediate, a strong and tenacious gluten becomes necessary. Wheat gluten proteins endow wheat with its unique properties that are exploited in pasta and bread making. They form up to 80% of the total protein in grain, being accumulated and stored in the grain’s starchy endosperm. The mixing of durum wheat semolina with water leads to the formation of a viscoelastic network of gluten proteins, which endows the dough with its essential rheological properties [4,5]. Gluten proteins are composed of monomeric gliadins, which are categorized into the following groups: alpha/beta, omega, and gamma gliadins, and the polymeric glutenins, which comprise high molecular weight
(HMW) and low molecular weight (LMW) glutenin subunits (GS) [4]. These protein subunits can also be categorized as either sulfur (S)-poor (HMW-GS and omega gliadins) or S-rich (LMW-GS and alpha/beta and gamma gliadins) [6], a classification which highlights the fundamental role of di-sulfide bonds for the formation of high molecular mass polymers [7]. Total protein content, the amount and ratios of the various storage proteins (mainly gliadins/glutenins, HMW/LMW-GS, and S-rich/S-poor proteins), and the percentage of grain unextractable polymeric proteins (UPP%, the proportion of large sodium dodecyl sulfate insoluble polymeric proteins) partly determine the levels of dough strength, elasticity, extensibility, and viscosity, all of which contribute to the biophysical and functional properties of dough and pasta quality [5].

The total grain protein content is not independent of the amounts and ratios of the various storage proteins. Indeed, although the overall partitioning of grain N between gliadins and glutenins in bread wheat at maturity has been demonstrated to be unmodified by growing conditions, it was found to vary according to differences in the total quantity of N accumulated in the grain [8,9]. However, these conclusions were reached by considering the pattern of accumulation of the various protein fractions during grain development for just one or very few bread wheat cultivars under different thermal and nutritional conditions. Hence, Triboi et al. [8] referred to the observed variation in grain N content as a consequence of grain growth, but they also implied that the protein composition of the grains of some bread wheat genotypes at harvest can simply be deduced by the total quantity of N per grain. At the same time, in the case of cultivars with different genotypes at the loci regulating the synthesis of different gluten fractions, quantitative ratios between fractions may differ despite the same total grain N content. In the case of Italian durum wheat cultivars, for example, a higher glutenin to gliadin ratio has been observed in modern compared with old cultivars, deriving from the expression of B-type LMW-GS, which was on average two times higher in the modern than in the old cultivars [10]. The importance of transcriptional control on the accumulation of gliadins and glutenins has also been noticed in bread wheat, and several putative loci identified that might regulate the differential expression of the different grain storage proteins [11].

Information on the existence of such relationships between total grain N and protein components and on the extent of environmental and genetic deviations from grain N allocation relationships is lacking with respect to durum wheat. To fill this gap, we analyzed the protein content and composition of a set of 16 durum wheat cultivars from different eras of breeding, hypothesized to retain a high level of genotypic variation in relation to grain protein content and composition. The set specifically contained a large number of old cultivars since the level of genetic variability in gluten characteristics present in old durum wheat cultivars from the Mediterranean basin is known to be much larger than that of modern cultivars [12]. At the same time, substantial differences exist in protein content and in grain quality traits between old and modern cultivars, since modern cultivars generally have lower protein percentages [13,14] but stronger gluten as a consequence of the introgression of more favorable alleles for both HMW and LMW-GS. The aim of our study was to find out whether, and to what extent, the genotypic variation in total grain N at maturity is also responsible for variations in the amounts and proportions of the different grain protein fractions and subunits analyzed.

## 2. Materials and Methods

### 2.1. *Site, Treatments, and Experimental Design*

The experiment was carried out during seasons 2015/16 (‘2016’) and 2016/17 (‘2017’) at the experimental station belonging to the University of Sassari, located in Ottava (41° N; 8° E; 80 m above sea level). The location is a typically Mediterranean environment, with a long-term mean annual rainfall of 539 mm—the majority of which occurs between October and April.

A total of 16 cultivars of durum wheat (*Triticum turgidum* L. ssp. *durum* Desf.), from different eras of wheat breeding and with different gluten allelic composition, were compared (Table 1). Most of them were tall cultivars, which were diffuse across Italy in the first half of the 20th century.
Two cultivars, Ichnusa and Maristella—deriving from crosses between *mediterraneum* and *syriacum* types [15], and widely grown in the 1960s—are generally shorter and earlier than the old cultivars, but still prone to lodging. Two semi-dwarf varieties, Svevo and Aureo, were also investigated, which were chosen from the many available semi-dwarf Italian durum wheat cultivars because of their reputation for producing good quality pasta [10].

<table><caption>**Table 1.** The investigated genotypes, their relative year of release, geographic or genetic origin, and details of the allelic composition of gluten genes.</caption><thead><tr><th>Name</th><th>Geographic or Genetic Origin</th><th>Year of Release</th><th>Glu-A1</th><th>HMW-GS-B1</th><th>LMW-GS-B3</th><th>Gliadin-A1</th></tr></thead><tbody><tr><td>Calabria</td><td>Calabria</td><td>-</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Dauno</td><td>Apulia</td><td>1914</td><td>Null</td><td>6 + 8</td><td>2</td><td>45</td></tr><tr><td>Russello</td><td>Selected from a Sicilian landrace</td><td>1936</td><td>Null</td><td>13 + 16</td><td>2</td><td>45</td></tr><tr><td>Saragolla <sup>(1)</sup></td><td>Apulia</td><td>-</td><td>Null</td><td>6 + 8; 20</td><td>2</td><td>42</td></tr><tr><td>Scorzonera</td><td>Indigenous landraces from Sicily</td><td>-</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Senatore Cappelli</td><td>Nord-african landrace Jean Retifah</td><td>1915</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Taganrog</td><td>Russia</td><td>-</td><td>Null</td><td>20*</td><td>2</td><td>45</td></tr><tr><td>Timilia</td><td>Sicily</td><td>-</td><td>2*</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Trigu arrubiu</td><td>Sardinia</td><td>-</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Trigu biancu</td><td>Sardinia</td><td>-</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Trigu murru</td><td>Sardinia</td><td>-</td><td>2<sup>+</sup></td><td>20*</td><td>2</td><td>45</td></tr><tr><td>Triminia</td><td>Sicily</td><td>-</td><td>2<sup>+</sup></td><td>7</td><td>2</td><td>45</td></tr><tr><td>Ichnusa</td><td>Biancale × Capeiti 8</td><td>1968</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Maristella</td><td>Dauno-III × Capeiti-8</td><td>1968</td><td>Null</td><td>20</td><td>2</td><td>45</td></tr><tr><td>Svevo</td><td>Sel. CIMMYT × Zenit sib</td><td>1996</td><td>Null</td><td>7 + 8</td><td>2</td><td>45</td></tr><tr><td>Aureo</td><td>Kofa × Svevo</td><td>2009</td><td>Null</td><td>6 + 8</td><td>2</td><td>45</td></tr></tbody></table>

Cultivars were sown on 10 December 2015 (2016 season), and 15 November 2016 (2017 season), at a sowing rate of 250 viable seeds m⁻², and at two fertilization rates: a low N level ('N46') with a single nitrogen application at the time of sowing (46 kg ha⁻¹ in the form of urea); a high N level ('N86'), which benefited from a second application of 40 kg ha⁻¹ at the onset of stem elongation in the form of ammonium nitrate.

The soils, just 0.4–0.5 m deep due to an underlying layer of limestone, have an organic matter content of 1.4 ± 0.3%, 45 kg ha⁻¹ of mineral nitrogen, 8.4 ± 0.5 ppm of available phosphorus and high total CaCO₃ (40 ± 4.4%). The available water amounted to 61 mm. In both seasons, the preceding crop was faba bean, and the sowing bed was prepared by ploughing to a depth of 0.25 m followed by surface cultivation. Weeds, pests, and diseases were chemically controlled.

Each plot consisted of six 1.2 m rows with an inter-row distance of 0.18 m. Plots were separated from each other by 50 cm all around to ensure the absence of interference between neighboring plots. The plots were set out as a split-plot design with three replications. Cultivars were assigned to the main plots and nitrogen rates to the sub-plots.

## 2.2. *Measurements and Analysis*

Anthesis date and physiological maturity were recorded when 50% of plants had reached these phenological stages, as ascertained from periodical inspections of the plots. Grain filling duration was estimated as the difference between maturity and anthesis.

Grain moisture content and grain weight were obtained from four 250 grain sub-samples per plot. Nitrogen percentage was determined on each sub-sample by means of a Carbon/Hydrogen/Nitrogen Analyzer (628 Series, LECO Corporation, St. Joseph, MI, USA). Nitrogen data were used to calculate grain protein percentage, on a dry basis, as N percentage × 5.7 and the amount of nitrogen (in µg) per grain (total grain N) as grain weight at dry basis × N percentage.

<footnote>(1) Saragolla seeds were heterogeneous regarding the HMW-GS pattern.</footnote>
### 2.2.1. Sequential Extraction of Gliadin and Glutenin for Reverse Phase-High Performance Liquid Chromatography (RP-HPLC) Analysis

Albumins-globulins, gliadins and glutenins were sequentially extracted from 50 mg of ground wholemeal grain as described by Wieser and Seilmeier [16] and separated at 50 °C on an ACE C18 column (250 × 2.1 mm, 5 µm, 300 Å) using an Alliance HPLC system (Waters). The elution system consisted of (A) trifluoroacetic acid (TFA, 0.1%, *v*/*v*) and (B) acetonitrile/TFA (99.9/0.1%, *v*/*v*). Linear gradient: 0 min 20% B, 20 min 60% B (albumins-globulins); 0 min 28% B, 30 min 56% B (gliadins, glutenin subunits). Flow rate was 0.2 mL min⁻¹; detection was UV absorbance at 210 nm. Following RP-HPLC separation, the area under the protein curve was calculated, corresponding to the amount of protein extracted in each fraction. For each grain sample, the protein extracted was then expressed as an amount per grain. Gluten amount was obtained by summing the gliadin and glutenin peaks. Total protein content was estimated by adding the albumin-globulin, gliadin, and glutenin peaks together. Areas corresponding to gliadin groups (alpha/beta, omega, and gamma) were calculated by cutting the whole gliadin chromatogram up as described by Samson et al. [17] with some adjustments: 7 to 14 min (omega); 14 to 22 min (alpha/beta); 22 min until the last peak (gamma). The same procedure was applied to glutenin chromatograms: 7 to 17 min (HMW-GS); 17 min to the last glutenin peak (LMW-GS) [18].

The relative percentage of each protein fraction (albumin-globulin, gliadin, and glutenin) was calculated from total protein, and the relative percentage of each gliadin group (i.e., alpha/beta, omega, and gamma) and the glutenin subunits (HMW and LMW) was calculated from total gluten.

Gliadin over glutenin ratio was calculated as the ratio of the gliadin peak over glutenin peak and the S-rich over S-poor ratio was calculated as the ratio of the sum of alpha/beta, gamma, and LMW-GS peaks over the sum of omega and HMW-GS peaks.

### 2.2.2. Extraction of Extractable and Unextractable Polymeric Proteins for Size Exclusion-High Performance Liquid Chromatography (SE-HPLC) Analysis

Extraction of polymeric protein and separation by SE-HPLC was conducted as described by Dachkevitch et al. [19] and Morel et al. [20] with minor modifications. Briefly, SDS-soluble proteins were extracted from 160 mg wholemeal flour by the addition of 20 mL 1% (*w*/*v*) SDS in 0.1 M sodium phosphate buffer (pH 6.9) and incubated on a rotary shaker (60 rpm at 60 °C for 80 min). Following centrifugation (37,000× g, 30 min, 20 °C), the supernatant was set aside for further SE-HPLC analysis. The unextractable proteins (or SDS insoluble proteins) were obtained by adding 5 mL SDS-phosphate buffer to the pellet, vortexing, and sonicating for 3 min at 7.5 watts, then centrifuging at 37,000× g at 20 °C for 30 min. The supernatant was subjected to SE-HPLC analysis and the pellet discarded. SE-HPLC separations were performed on an Alliance system (Waters) equipped with a TOSOBIOSCIENCE TSKgel G4000SWXL column (8 × 300 mm, 8 µm, 300 Å) protected by a guard column TOSOBIOSCIENCE TSKgel SWXL guard (6 × 40 mm, 7 µm mm). Proteins were eluted at ambient temperature with 0.1 M sodium phosphate buffer (pH 6.9) containing 0.1% (*w*/*v*) SDS at a constant flow of 0.7 mL min⁻¹ and absorbance was recorded at 214 nm.

The total area under each chromatogram obtained from SDS extractable and unextractable protein extracts was expressed as a percentage (EP% and UP%, respectively) of the sum of the total area of both chromatograms. The unextractable polymeric protein fraction (UPP%) was calculated as a percentage of the total polymeric proteins (UPP% = UP%/(sum of glutenins % of EP) + UP%).

Grain nitrogen data and the percentages of protein fractions calculated by RP-HPLC were used to estimate the µg of nitrogen of each protein fraction, considering that we can express proteins in terms of nitrogen and vice versa (grain protein = grain N × 5.7) and that the proportion of non-protein nitrogen in wheat grain at maturity and in the semolina is very low to negligible [21–23]. The nitrogen contents of EP and UP were calculated using their percentages (calculated by the SE-HPLC) multiplied by grain nitrogen (total grain N). The nitrogen contents of albumin-globulins, gliadins, and glutenins were calculated using their percentages calculated by RP-HPLC multiplied by the nitrogen content of EP.
The nitrogen contents of the gliadin groups and the glutenin subunits (alpha/beta, omega, and gamma and HMW and LMW-GS) were calculated using their percentages calculated by RP-HPLC multiplied by the nitrogen content of gluten. From here on, the nitrogen content of a fraction or subunit will simply be referred to as ‘content’; for example, the µg of nitrogen of gliadin will be referred to as ‘gliadin content’.

### 2.2.3. Electrophoresis

Gliadins were analyzed by A-PAGE, and the gliadin extraction was carried out on 30 mg durum wheat flour as described by Morel [24]. HMW-GS were separated by SDS-PAGE, and protein extraction was performed on 20 mg durum wheat flour as described by Singh et al. [25]. The identification of high-molecular-weight glutenin subunit (HMW-GS) alleles was based on the classification proposed by Payne and Lawrence [26].

### 2.3. Statistical Analysis

Agronomic data were analyzed by combined ANOVA, using R software [27], according to a split-split plot design with three replications, where years were assigned to the main plot, cultivars to the sub plots, and nitrogen treatment to the sub-sub plots. Cultivar means were separated by means of a multiple *t*-test once the level of statistical significance of the cultivar effect had been ascertained. The Pearson correlation coefficient was used to evaluate the existence of causal relationships between pairs of traits.

## 3. Results

### 3.1. Weather Conditions and Phenology

The rainfall recorded between October and May amounted to 286 mm in 2017 and to 363 mm in 2016, i.e., 60% and 76%, respectively, of the long-term mean (40 years) (Figure 1).

Rainfall was particularly scarce in the springs of both seasons; a total of only 8 mm of rain fell in April and May in 2017, and 25 mm in 2016. 2017 was also characterized by higher maximum temperatures in May (24.1 °C compared with 21.2 °C in 2016), when most of the grain filling occurs, as well as by a greater number of days in which the maximum temperature exceeded 25 °C (17 vs. 4 days in 2016).

Phenology was not affected by the N treatment, whereas the effect of Cultivar and Cultivar × Year interaction were both significant as shown by ANOVA (Table 2).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Rainfall (blue bars), reference evapotranspiration (ETo) (orange bars), and rainfall deficit (rainfall–ETo) (grey bars) between October and May for the years 2016 and 2017, plus the long-term mean (40 years; 1970–2010).</figcaption></figure>
<table><caption>**Table 2.** Anthesis date and grain filling duration. Least Square Means and LSD<sub>0.05</sub> for the comparison of cultivar × year means at the same or different levels of year.</caption><thead><tr><th rowspan="2">Cultivar</th><th colspan="2">Anthesis (Doy) <sup>(1)</sup></th><th colspan="2">Grain Filling Duration (Day)</th></tr><tr><th>2016</th><th>2017</th><th>2016</th><th>2017</th></tr></thead><tbody><tr><td>Calabria</td><td>124</td><td>121</td><td>32</td><td>30</td></tr><tr><td>Cappelli</td><td>125</td><td>122</td><td>31</td><td>30</td></tr><tr><td>Dauno</td><td>127</td><td>124</td><td>30</td><td>27</td></tr><tr><td>Russello</td><td>125</td><td>123</td><td>32</td><td>29</td></tr><tr><td>Saragolla</td><td>126</td><td>125</td><td>32</td><td>26</td></tr><tr><td>Scorzonera</td><td>128</td><td>127</td><td>29</td><td>25</td></tr><tr><td>Taganrog</td><td>126</td><td>124</td><td>31</td><td>27</td></tr><tr><td>Timilia</td><td>116</td><td>114</td><td>33</td><td>34</td></tr><tr><td>Trigu arrubiu</td><td>123</td><td>123</td><td>33</td><td>28</td></tr><tr><td>Trigu biancu</td><td>125</td><td>124</td><td>30</td><td>28</td></tr><tr><td>Trigu murru</td><td>123</td><td>123</td><td>32</td><td>28</td></tr><tr><td>Triminia</td><td>123</td><td>125</td><td>33</td><td>27</td></tr><tr><td>Ichnusa</td><td>109</td><td>109</td><td>39</td><td>38</td></tr><tr><td>Maristella</td><td>110</td><td>110</td><td>40</td><td>37</td></tr><tr><td>Svevo</td><td>107</td><td>110</td><td>40</td><td>35</td></tr><tr><td>Aureo</td><td>107</td><td>110</td><td>42</td><td>36</td></tr><tr><td>Mean</td><td>120</td><td>120</td><td>34</td><td>30</td></tr><tr><td>LSD<sub>0.05</sub></td><td colspan="2">2</td><td colspan="2">2</td></tr></tbody></table>

In relation to anthesis date, the cultivars could be divided into two groups: (i) those with a mean anthesis date that fell within the first week of May (in the 2017 season, this resulted in less than one month for grain filling in the case of some cultivars), which included almost all the tall and old cultivars; (ii) those with a mean anthesis date between April 17 and 20, which enabled a longer grain filling period compared with the first group, and comprised the four more recent constitutions. These differences in anthesis date resulted in different thermal conditions during grain filling: the first group of cultivars (made up of the old, later varieties) was exposed to higher maximum temperatures compared with the earlier flowering group (24.5 vs. 22.3 °C in 2017, and 21.7 vs. 20.4 °C in 2016) and had more days in which the maximum temperature exceeded 25 °C, especially in 2017 (14 vs. 10 days in 2016).

## 3.2. Characterization of Gliadin and Glutenin Alleles

Based on their mobility on the acid and SDS-PAGE, six HMW-GS patterns were revealed: 20, 20*, 7, 13 + 16, 6 + 8, and 7 + 8. Particularly, one variety, Saragolla, produced a heterogeneous HMW-GS pattern (6 + 8; 20) that was different to the one reported by De Santis et al. [10] This is probably because landraces tend to be mixtures of different pure lines. All the genotypes studied had the superior LMW2 associated with gliadin γ-45, except for Saragolla, which had the unusual combination of LMW2/gliadin γ-42 (Table 1).

## 3.3. Grain Nitrogen and Its Partitioning

No effect of year was detected on the nitrogen content of the grain in spite of the more severe terminal drought and higher temperatures characterizing the grain filling period in 2017, whereas both grain nitrogen content and protein percentage increased at the higher N rate (Table 3).

<footnote>(1) Doy: day of the year.</footnote>
<table><caption>**Table 3.** Grain protein percentage, total grain N, N in gliadins, glutenins and albumins-globulins (µg of N in each component), the ratios gliadins to glutenins (GLI/GLU), high molecular to low molecular weight glutenins subunits (HMW/LMW) and sulfur rich to sulfur poor (Srich/Spoor), and unextractable polymeric proteins percentage (UPP%). Least Square Means and results of the analysis of variance.</caption><thead><tr><th></th><th>Grain Protein (%)</th><th>Total N</th><th>N in Gliadins</th><th>N in Glutenins</th><th>N in ALBGLOB</th><th>GLI/GLU</th><th>HMW/LMW</th><th>Srich/Spoor</th><th>UPP%</th></tr><tr><th></th><th colspan="9" style="text-align: center;">(µg grain<sup>-1</sup>)</th></tr></thead><tbody><tr><td>Year (Y)</td><td>ns</td><td>ns</td><td>**</td><td>ns</td><td>ns</td><td>**</td><td>ns</td><td>*</td><td>***</td></tr><tr><td>2016</td><td>14.6</td><td>1110</td><td>383 b</td><td>468</td><td>179</td><td>0.85 b</td><td>0.19</td><td>6.16 b</td><td>28.5 a</td></tr><tr><td>2017</td><td>14.8</td><td>1146</td><td>437 a</td><td>458</td><td>171</td><td>1.03 a</td><td>0.18</td><td>6.43 a</td><td>27.5 b</td></tr><tr><td>Nitrogen (N)</td><td>***</td><td>**</td><td>*</td><td>***</td><td>**</td><td>**</td><td>ns</td><td>***</td><td>***</td></tr><tr><td>N46</td><td>14.3 b</td><td>1093 b</td><td>433 a</td><td>434 b</td><td>166 b</td><td>1.04 a</td><td>0.18</td><td>6.62 a</td><td>28.4 a</td></tr><tr><td>N86</td><td>15.2 a</td><td>1138 a</td><td>387 b</td><td>492 a</td><td>184 a</td><td>0.84 b</td><td>0.19</td><td>5.98 b</td><td>27.6 b</td></tr><tr><td>Cultivar (C)</td><td>***</td><td>***</td><td>***</td><td>***</td><td>***</td><td>***</td><td>***</td><td>***</td><td>***</td></tr><tr><td>Calabria</td><td>14.3 df</td><td>1158 bc</td><td>405 bc</td><td>493 b</td><td>165 be</td><td>0.83 cf</td><td>0.15 gi</td><td>7.67 ab</td><td>30.2 ab</td></tr><tr><td>Cappelli</td><td>14.6 df</td><td>1190 bc</td><td>420 ac</td><td>498 b</td><td>170 be</td><td>0.85 cf</td><td>0.15 hi</td><td>7.52 bc</td><td>30.4 ab</td></tr><tr><td>Dauno</td><td>16.8 ab</td><td>1235 b</td><td>473 ab</td><td>498 b</td><td>185 bd</td><td>0.98 bd</td><td>0.24 c</td><td>3.91 h</td><td>25.5 e</td></tr><tr><td>Russello</td><td>15.0 ce</td><td>1090 ce</td><td>418 ac</td><td>450 bd</td><td>155 de</td><td>0.93 ce</td><td>0.20 d</td><td>6.45 eg</td><td>25.3 e</td></tr><tr><td>Saragolla</td><td>15.8 bc</td><td>1160 bc</td><td>458 ab</td><td>440 bd</td><td>185 bd</td><td>1.09 ac</td><td>0.26 b</td><td>3.53 h</td><td>25.7 e</td></tr><tr><td>Scorzonera</td><td>14.5 df</td><td>1000 df</td><td>423 ac</td><td>360 ef</td><td>158 ce</td><td>1.26 ab</td><td>0.21 d</td><td>6.67 ce</td><td>25.9 e</td></tr><tr><td>Taganrog</td><td>15.2 cd</td><td>1203 bc</td><td>458 ab</td><td>483 bc</td><td>183 bd</td><td>0.96 bd</td><td>0.17 ef</td><td>8.29 a</td><td>26.5 ce</td></tr><tr><td>Timilia</td><td>13.8 cd</td><td>1125 a</td><td>455 ab</td><td>385 de</td><td>200 ab</td><td>0.80 cf</td><td>0.18 e</td><td>6.27 cf</td><td>29.0 ac</td></tr><tr><td>Trigu arrubiu</td><td>14.8 cf</td><td>1245 b</td><td>378 bc</td><td>585 a</td><td>195 bc</td><td>0.65 ef</td><td>0.15 hi</td><td>7.39 bd</td><td>28.5 ab</td></tr><tr><td>Trigu biancu</td><td>15.2 f</td><td>1380 bd</td><td>448 ab</td><td>603 a</td><td>238 a</td><td>1.26 ab</td><td>0.16 fh</td><td>6.65 eg</td><td>28.1 bd</td></tr><tr><td>Trigu murru</td><td>17.3 a</td><td>1453 a</td><td>523 a</td><td>613 a</td><td>238 a</td><td>0.95 ce</td><td>0.33 a</td><td>3.62 h</td><td>23.7 e</td></tr><tr><td>Triminia</td><td>15.1 ce</td><td>1030 de</td><td>445 ab</td><td>353 ef</td><td>168 be</td><td>1.35 a</td><td>0.22 d</td><td>5.98 fg</td><td>26.5 de</td></tr><tr><td>Ichnusa</td><td>14.4 df</td><td>1158 bc</td><td>425 ac</td><td>478 bc</td><td>160 ce</td><td>0.91 cf</td><td>0.14 ij</td><td>7.39 bd</td><td>29.9 ab</td></tr><tr><td>Maristella</td><td>12.5 g</td><td>893 f</td><td>255 d</td><td>410 ce</td><td>155 de</td><td>0.63 f</td><td>0.14 ij</td><td>6.18 eg</td><td>30.9 a</td></tr><tr><td>Svevo</td><td>12.4 g</td><td>735 g</td><td>258 d</td><td>310 f</td><td>110 f</td><td>0.85 cf</td><td>0.16 fg</td><td>6.03 dg</td><td>28.8 ab</td></tr><tr><td>Aureo</td><td>14.0 ef</td><td>993 ef</td><td>323 cd</td><td>453 bd</td><td>138 ef</td><td>0.7 df</td><td>0.13 j</td><td>6.14 g</td><td>28.9 ab</td></tr><tr><td>Y × C</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td></tr><tr><td>N × C</td><td>ns</td><td>ns</td><td>*</td><td>*</td><td>ns</td><td>*</td><td>ns</td><td>*</td><td>ns</td></tr><tr><td>N × Y</td><td>***</td><td>***</td><td>***</td><td>ns</td><td>**</td><td>***</td><td>ns</td><td>ns</td><td>*</td></tr></tbody></table>

Level of significance: * *p* ≤ 0.05; ** *p* ≤ 0.01, *** *p* ≤ 0.001; ns: not significant. Means with the same letter are not statistically different as assessed using the Student Test for *p* ≤ 0.05.

The ranking of the cultivars in relation to the content and percentage of N in the grains was not influenced by the year, nor by the nitrogen treatment, as demonstrated by the lack of any significant interaction of N and/or year with cultivar. The genotypic variation expressed by the 16 cultivars on the average of the two years and two N treatments was large for both grain protein percentage (12.4–17.3%) and total grain N content (735–1453 µg N grain⁻¹), and was mirrored in a significant effect of cultivar for all the traits presented in Table 3. The differences in grain protein percentage between the cultivars was derived from the corresponding differences in grain N content, the two being associated traits, with an *r* = 0.81 ***.

A lower gliadin content was observed at the higher N rate (N86) (Table 3) due to a lower gamma gliadin content, while alfa/beta and omega groups remained unchanged (Table 4).

<table><caption>**Table 4.** Grain content of the gliadin groups (alfa/beta, omega, and gamma), and of the glutenin subunits (high- and low-molecular weight: HMW and LMW) (µg of N in each component). Results from Least Square Means and the analysis of variance ANOVA.</caption><thead><tr><th></th><th>N in Alfa/Beta Gliadins<br>(µg)</th><th>N in Omega Gliadins<br>(µg)</th><th>N in Gamma Gliadins<br>(µg)</th><th>N in HMW-GS<br>(µg)</th><th>N in LMW-GS<br>(µg)</th></tr></thead><tbody><tr><td>Year (Y)</td><td>**</td><td>ns</td><td>**</td><td>ns</td><td>ns</td></tr><tr><td>2016</td><td>202</td><td>51.6</td><td>129</td><td>73.9</td><td>394</td></tr><tr><td>2017</td><td>225</td><td>58.9</td><td>152</td><td>70.8</td><td>388</td></tr><tr><td>Nitrogen (N)</td><td>ns</td><td>ns</td><td>*</td><td>***</td><td>***</td></tr><tr><td>46</td><td>222</td><td>55.1</td><td>156</td><td>66.3</td><td>368</td></tr><tr><td>86</td><td>206</td><td>55.4</td><td>125</td><td>78.4</td><td>413</td></tr></tbody></table>

Level of significance: * *p* ≤ 0.05; ** *p* ≤ 0.01, *** *p* ≤ 0.001; ns: not significant.
However, this negative effect of the fertilization treatment on gliadins was only true in 2016 (464 µg grain⁻¹ at N46 vs. 302 µg grain⁻¹ at N86), whereas the opposite was true in the drier 2017 season. In 2017, the gliadin content at N86 was 14.6% higher than that at N46 (400 µg grain⁻¹ at N46 vs. 471 µg grain⁻¹ at N86) for all cultivars except for Trigu biancu, Trigu murru, and Trigu arrubiu. Consequently, significant N × Year and N × C interactions were generated (Table 3).

The glutenin plus albumin-globulin content was higher at N86 than at N46 (12% and 10%, respectively). However, the glutenin contents of Cappelli and Dauno were negatively affected by the additional N fertilization in both years, generating an N × C interaction (Table 3). Both glutenin subunits (HMW and LMW) varied, by the same degree, with N fertilization (Table 4). The described variation in gliadins and glutenins caused the gliadins to glutenins ratio (GLI/GLU) to be higher at N46 (1.04) than at N86 (0.84) (Table 3). The gliadins constituted the only protein fraction to be affected by year, with a higher gliadin content being observed in 2017—the year characterized by severe water and temperature stress (437 vs. 383 µg grain⁻¹ in 2016). The content of the different gliadin groups increased accordingly, with the exception of omega gliadins. As a consequence, GLI/GLU was higher in 2017. UPP% was lower at N86 and in 2017—the drier and hotter year. The negative effect of an additional N application on UPP% was stronger in 2017 than in 2016 for all the cultivars studied.

To evaluate whether the genotypic variation in total grain N also implied a variation in the various protein fractions, cultivar means were used to calculate the relationships between albumin-globulin, gliadins, and glutenins with total N content. Given the significant cultivar × N interaction for gliadins and glutenins, as shown by ANOVA (Table 3), we first regressed the two N treatments separately, but because their slopes and intercepts were not statistically different, only the results for a regression based on cultivar means across year and N treatments are discussed (Figure 2).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Relationships between the variation in the total grain nitrogen and the N content of the three protein fractions: glutenins (GLU) (squares), gliadins (GLI) (circles), and albumins-globulins (ALBGLOB) (triangles). Points are cultivar means across years and N treatments. *** *p* ≤ 0.001.</figcaption></figure>

These relationships provided information about the strength of the relationship between genotypic variation in protein fractions and the total grain N (coefficient of correlation and determination), and the extent of the variation in protein fractions for each unit of variation of total grain N, or partitioning coefficient (slope of the regression).

The albumins-globulins protein fraction was the one most tightly associated with total grain N, which accounted for 82% of its genotypic variation, compared with 78% for gliadins and 65% for glutenins. On the other hand, the albumins-globulins fraction was also the one with the lowest slope of regression (0.17 ± 0.02) compared with the storage proteins, gliadins (0.34 ± 0.07), and glutenins (0.44 ± 0.06); i.e., it was the fraction that varied the least in response to the genotypic variation in total
grain N. The relationship between total grain N and gliadin was stronger than that between total grain N and glutenins, and no significant differences were detected between either the slopes of these two relationships, or between their intercepts, which were not different from zero. A significant effect of cultivars on the GLI/GLU ratio was detected by ANOVA, reflecting the variation in GLI/GLU from 0.63 to 1.35 (Table 3), but this variation in GLI/GLU was not associated with the variation in total grain N (data not shown).

Within the gliadin groups, the sulfur(S)-rich alpha/beta and gamma were generally present in larger amounts (101–261 µg for the alpha/beta, and 79–166 µg for the gamma gliadins), than the S-poor omega gliadins (31–108 µg) (Figure 3). The relationship between total gliadin content and alpha/beta gliadins was stronger ($R^2$ = 0.89 **)) than the relationship with gamma gliadins ($R^2$ = 0.68 **)), whereas omega gliadin levels varied independently from the total gliadin content, in spite of the large genotypic variation in this component. The slope of regression of alpha/beta (0.56 ± 0.05) was significantly different and two times higher than that of gamma (0.27 ± 0.05).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Relationships between the variation in the µg of total gliadins per grain and the content of the three gliadin components: alpha/beta (squares), gamma (triangles), and omega (circles). Points are cultivars means across years and N treatments. ***, *p* ≤ 0.001; ns, not significant.</figcaption></figure>

The S-rich LMW-GS were present in larger amounts than the S-poor HMW-GS (Figure 4) and showed greater genotypic variation: 270–510 µg grain⁻¹ (coefficient of variation = 36 vs. 19% for HMW-GS). The genotypic variation in grain N allocated to glutenins was more tightly associated with the variation in LMW-GS ($R^2$ = 0.93 **)) than in HMW-GS ($R^2$ = 0.47 **)).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Relationships between the variation in the µg of total glutenins per grain and the amount of the two glutenin components: LMW-GS (triangles) and HMW-GS (circles). Points are cultivar means across years and N treatments. *** *p* < 0.001, ** *p* < 0.01.</figcaption></figure>
The slope of regression of the content of LMW-GS per grain (0.80 ± 0.06) was more than four times greater than that of the content of HMW-GS (0.20 ± 0.06). In spite of these linear relationships, the HMW/LMW ratio was significantly different between genotypes as shown by ANOVA, ranging from 0.13 to 0.33 (Table 3), but no relationship was found between that ratio and total N grain glutenin content.

The S-rich to S-poor ratio did not vary with total grain N, but it was significantly different among genotypes and varied from 3.7 to 8.3 (Table 3). Similarly, UPP% did not vary with total grain N, although it was negatively correlated with HMW/LMW (r = 0.88 **)) and positively correlated with S-rich/S-poor (r = 0.93 **)).

4. Discussion

Grain protein percentage makes the largest contribution (40%) to the EU Quality index for durum wheat (European Commission Regulation No. 2237/2003, 23 December 2003), and in some countries, grain protein concentration influences the amount of money paid to wheat farmers due to the importance of this trait for the quality of both pasta and bread. In bread wheat, a higher protein content is generally equated with higher protein quality for bread making [28], consistent with the association between grain N content and protein fraction composition, as first highlighted by Triboi et al. [29] and then used by Martre et al. [30] to simulate grain N composition in the Sirius Quality model. An important finding obtained from this approach was that the water availability, the temperature during grain filling, and the nitrogen fertilization rate can only affect protein composition through a variation in total grain N. All of these relationships were investigated in the present paper in durum wheat, including testing for the effect of genotype, because not only did we expect genotype to affect the proportions of gluten protein fractions and components, we also expected it to interact with the environment.

4.1. *Nitrogen and Season Effect*

The effect of N fertilization on grain protein fractions across seasons and genotypes was mediated by an effect on the total grain N, as previously suggested by Triboi et al. [8], and consisted in a change in the GLI/GLU ratio, as glutenins increased and the gliadin fraction decreased. This result lies in contrast with some previous studies on bread wheat [31–33]. The significant Cultivar × N and N × Year interactions for gliadin and glutenin content in our study suggest that the effect of N fertilization on grain protein is strongly dependent on the specific cultivar considered, in agreement with Pechanek et al. [34] and Wieser and Seilmeier [16], and on the environmental conditions under which it is grown.

N fertilization also affected the quality parameters S-rich/S-poor and UPP%. Our results agree with those of Johansson et al. [35,36], who showed that N fertilization decreases the UPP%, and with other authors who found that N supply leads to a lower S-rich and higher S-poor protein component [37,38]. In bread wheat grown under controlled conditions, this change in S-rich/S-poor ratio was not due to a S deficiency, and the addition of post-anthesis S had no effect on protein composition [36]. On the other hand, our results are in contrast with those of Gagliardi et al. [39], who showed that N supply increases UPP% in durum wheat.

We found a similar total grain N and grain protein percentage between the two years in spite of
the more severe stress conditions during grain filling and the shorter grain filling duration encountered
in 2017. This shorter grain filling was the likely cause of the lower UPP% observed in 2017 since the
accumulation of SDS-insoluble large polymer proteins in the grain increases in the final stages of grain
growth and is related to continuous grain dehydration [40]. The only protein fraction affected by
year was the gliadins—as also noted by Daniel and Triboi [41]—specifically the S-rich gliadin groups,
which led to an increase in the S-rich/S-poor ratio. However, the interactions between year and N
supply in grain protein components and the UPP% make the interpretation of the effect of each factor
alone complicated, as also noted by Johansson et al. [36] and Malik et al. [42]. On the other hand,
no Cultivar × Year interaction was observed for any of the protein fractions considered in spite of the different environmental conditions experienced by the set of cultivars analyzed as a consequence of their different anthesis dates.

## 4.2. Genotypic Effect

The variation in grain protein percentage and absolute content was considerable for the set of cultivars investigated. This allowed us to assess the dependence of genotypic variation in grain protein percentage on the variation in total grain N, and to describe the relationships between total grain N and protein fractions for durum wheat. By expressing nitrogen on a per grain basis in µg rather than as a percentage, we were also able to exclude the confounding effect of the variation in carbohydrate accumulation—the processes of N and carbohydrate accumulation being largely independent of each other [43].

As total grain N varied among genotypes, the variation in the content of the storage proteins exceeded that of the metabolic proteins; i.e., each µg of total grain N difference between cultivars was associated with an increase in the µg of N in gliadins and glutenins, and this increase was more than twice the increase in µg of N observed in albumins-globulins. Indeed, metabolic protein accumulation is known to be sink limited [9,30], in contrast with storage proteins which are source-limited.

The genotypic variation in total grain N was more tightly associated with the variation in gliadins than in glutenins, at least in the set of cultivars and environmental conditions analyzed here. This meant that genotypic deviations for glutenins were large, a result that can be explained by the presence of polymorphisms in the HMW-B1 glutenin subunits (six different HMW-B1 types among our 16 cultivars), which also implied differential expression levels of the specific storage proteins [10].

Interestingly, the partitioning coefficients calculated in this experiment between N in gliadins and total grain N (0.34 ± 0.07), and N in glutenins and total grain N (0.44 ± 0.06), were similar to those calculated by Triboi et al. [8] for the second phase of the grain filling period of bread wheat (0.34 ± 0.01 for gliadins and 0.46 ± 0.01 for glutenins). This means that the relationships described by Triboi et al. [8] between the storage protein fractions and total N of developing grains of bread wheat grown under varying environmental conditions (post-anthesis temperature and water availability, and rate and timing of N fertilization) also hold for mature grains of durum wheat when the variation is induced by genotypic differences.

The proportions of the different gliadin groups in the total gliadins were the highest for alpha/beta compared to gamma or omega, which is in accordance with the literature on bread wheat [29,34,44], and the proportion of glutenin subunits in total glutenins was higher for LMW than HMW-GS, as already found in bread wheat [29].

The genotypic variation in both gliadins and glutenins was more strongly associated with their respective S-rich components (alfa-beta gliadins gamma gliadins, and LMW-GS) than with their S-poor components. The S-poor omega subunit did not vary with the variation in total µg of gliadins per grain, in accordance with Uthayakumaran et al. [45], and the coefficient of determination between total glutenins and LMW-GS was double the value obtained for HMW-GS. The S concentration of the grain is an important indicator of the size distribution of the gluten macropolymers that govern the rheological properties of dough [37,45–47]. Therefore, it is reasonable to suppose that durum wheat cultivars with high total grain N and gluten content would also produce a more extensible dough due to their high content of S-rich proteins [6].

The significant genotypic differences in the GLI/GLU, S-rich/S-poor, and HMW/LMW ratios and in UPP% were not influenced by the corresponding variation in grain N content. In contrast with our results, Blumenthal et al. [48], Jia et al. [33], and Triboi et al. [29] found a positive relationship between GLI/GLU and grain N content, and Pechanek et al. [34] found a positive relationship between HMW/LMW and grain N content. In accordance with our results, Triboi et al. [29] did not identify any relationship between the HMW/LMW ratio and grain N content, and Zhang et al. [49] agree that UPP% largely depends on genotype. The reasons for this discrepancy could be a consequence of the
different species analyzed (bread wheat was addressed in the cited papers, whereas durum wheat was the subject of the present study) and/or in the size of the sets of cultivars studied, which was smaller in most of the cited papers, whereas the present study included cultivars from different breeding eras, which therefore expressed large polymorphisms in their protein fractions.

According to Wrigley et al. [50], UPP% is greatly affected by the gluten allelic composition of the cultivars. However, our modern varieties, Aureo and Svevo, contain the improved versions of the glutenin HMW alleles at the Glu-B1 locus (6 + 8 and 7 + 8, respectively) and had high UPP% values that did not differ to those for Calabria and Cappelli, which have the HMW 20 genotype at the Glu-B1 locus that is considered weaker than both HMW 7 + 8 and 6 + 8 for producing good quality pasta and bread [51,52].

## 5. Conclusions

In the set of cultivars studied and environmental conditions explored, the linear relationships calculated between protein fractions and total grain N accounted for the majority of the variation in protein composition, but cultivar-dependent deviations with respect to these relationships were generally large enough as to significantly influence the ratios between fraction and, hence, presumably, gluten quality. The variation in protein fractions associated with a variation in total grain N content depends on the cause of the variation in µg N per grain. A greater N availability implies an increase in glutenin and a concomitant decrease in the average gliadin content, although this response was found to be genotype-dependent. On the contrary, when a cultivar-dependent increase in N content occurs, the relationship with total grain N is positive for both gliadins and glutenins, although the S-rich gluten components are more strongly associated with total grain N.

**Author Contributions:** Conceptualization, F.G. and M.M.; methodology, F.G., M.-F.S. and M.-H.M.; formal analysis, M.M. and R.M.; data curation, M.M.; writing—original draft preparation, M.M. and R.M.; writing—review and editing, R.M. and F.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** This research represents a chapter of a PhD project carried out by Marina Mefleh at the Doctoral School in Agriculture Sciences at the University of Sassari. We are extremely grateful to Antonio Canalis and to the technicians, A. Piredda, A. Ara and D. Carta, of the Unit of Agronomia, Coltivazioni erbacee e genetica at the University of Sassari for their assistance in field management and laboratory analysis.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Gallo, G.; Bianco, M.L.; Bognanni, R.; Saimbene, G.; Orlando, A.; Grillo, O.S.C.A.R.; Saccone, R.; Venora, G. Durum wheat bread: Old sicilian varieties and improved ones. In Proceedings of the 5th International Tech. Symp. Food Process. Monit. Technology Bioprocesses Food Qual. Management, Potsdam, Germany, 31 August–2 September 2009; Volume 4, pp. 1128–1138.

2. Peña, R.J.; Trethowan, R.; Pfeiffer, W.H.; Van Ginkel, M. Quality (End-Use) Improvement in Wheat. *J. Crop. Prod.* **2002**, *5*, 1–37. [CrossRef]

3. Dexter, J.E.; Matsuo, R.R. Influence of protein content on some durum wheat quality parameters. *Can. J. Plant Sci.* **1977**, *57*, 717–727. [CrossRef]

4. Wieser, H. Chemistry of gluten proteins. *Food Microbiol.* **2007**, *24*, 115–119. [CrossRef] [PubMed]

5. Mefleh, M.; Conte, P.; Fadda, C.; Giunta, F.; Motzo, R. From seed to bread: Variation in quality in a set of old durum wheat cultivars. *J. Sci. Food Agric.* **2019**, *100*, 4066–4074. [CrossRef]

6. Naeem, H.A.; MacRitchie, F. Sulphur nutrition and quality attributes of wheat. In *Sulphur in Plants*; Abrol, Y.P., Ahmad, A., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003; pp. 303–322.

7. Shewry, P.; Tatham, A. Disulphide Bonds in Wheat Gluten Proteins. *J. Cereal Sci.* **1997**, *25*, 207–227. [CrossRef]

8. Triboi, A.-M.; Martre, P.; Triboï-Blondel, A. Environmentally-induced changes in protein composition in developing grains of wheat are related to changes in total protein content. *J. Exp. Bot.* **2003**, *54*, 1731–1742. [CrossRef] [PubMed]
9. Martre, P.; Porter, J.R.; Jamieson, P.D.; Triboi, E. Modeling Grain Nitrogen Accumulation and Protein Composition to Understand the Sink/Source Regulations of Nitrogen Remobilization for Wheat. *Plant Physiol.* **2003**, *133*, 1959–1967. [CrossRef]

10. De Santis, M.A.; Giuliani, M.M.; Giuzio, L.; De Vita, P.; Lovegrove, A.; Shewry, P.R.; Flagella, Z. Differences in gluten protein composition between old and modern durum wheat genotypes in relation to 20th century breeding in Italy. *Eur. J. Agron.* **2017**, *87*, 19–29. [CrossRef]

11. Plessis, A.; Ravel, C.; Bordes, J.; Balfourier, F.; Martre, P. Association study of wheat grain protein composition reveals that gliadin and glutenin composition are trans-regulated by different chromosome regions. *J. Exp. Bot.* **2013**, *64*, 3627–3644. [CrossRef]

12. Nazco, R.; Pena, R.J.; Ammar, K.; Villegas, D.; Crossa, J.; Moragues, M. Durum wheat (Triticum durum Desf.) Mediterranean landraces as sources of variability for allelic combinations at Glu-1/Glu-3 loci affecting gluten strength and pasta cooking quality. *Genet. Resour. Crop Evol.* **2014**, *61*, 1219–1236. [CrossRef]

13. Giunta, F.; Pruneddu, G.; Motzo, R. Grain yield and grain protein of old and modern durum wheat cultivars grown under different cropping systems. *Field Crop. Res.* **2019**, *230*, 107–120. [CrossRef]

14. Mefleh, M.; Conte, P.; Fadda, C.; Giunta, F.; Piga, A.; Hassoun, G.; Motzo, R. From ancient to old and modern durum wheat varieties: Interaction among cultivar traits, management, and technological quality. *J. Sci. Food Agric.* **2018**, *99*, 2059–2067. [CrossRef] [PubMed]

15. Ali Dib, T.; Monneveux, P.; Araus, J.L. Adaptation à la sècheresse et notion d'ideotype chez le blé dur. II. Caractères physiologiques d'adaptation. *Agronomie* **1992**, *12*, 381–393. [CrossRef]

16. Wieser, H.; Seilmeier, W. The influence of nitrogen fertilization on quantities and proportions of different protein types in wheat flour. *J. Sci. Food Agric.* **1998**, *76*, 49–55. [CrossRef]

17. Samson, M.-F.; Mabille, F.; Chéret, R.; Abecassis, J.; Morel, M.-H. Mechanical and Physicochemical Characterization of Vitreous and Mealy Durum Wheat Endosperm. *Cereal Chem. J.* **2005**, *82*, 81–87. [CrossRef]

18. Wieser, H.; Antes, S.; Seilmeier, W. Quantitative Determination of Gluten Protein Types in Wheat Flour by Reversed-Phase High-Performance Liquid Chromatography. *Cereal Chem. J.* **1998**, *75*, 644–650. [CrossRef]

19. Dachkevitch, T.; Autran, J.-C. Prediction of baking quality of bread wheats in breeding programs by size-exclusion high-performance liquid chromatography. *Cereal Chem.* **1989**, *66*, 448–456.

20. Morel, M.-H.; Dehlon, P.; Autran, J.C.; Leygue, J.P.; Bar-L'Helgouac'H, C. Effects of Temperature, Sonication Time, and Power Settings on Size Distribution and Extractability of Total Wheat Flour Proteins as Determined by Size-Exclusion High-Performance Liquid Chromatography. *Cereal Chem. J.* **2000**, *77*, 685–691. [CrossRef]

21. Jennings, A.; Morton, R. Changes in Carbohydrate, Protein, and Non-Protein Nitrogenous Compounds of Developing Wheat Grain. *Aust. J. Biol. Sci.* **1963**, *16*, 318. [CrossRef]

22. Mariotti, F.; Tome, D.; Mirand, P. Converting Nitrogen into Protein -Beyond 6.25 and Jones' Factors. *Crit. Rev. Food. Sci. Nutr.* **2008**, *48*, 177–184. [CrossRef]

23. Bell, P.M. Non-protein-nitrogen in wheat flour. *JFSA* **2006**, *14*, 133–137. [CrossRef]

24. Morel, M.H. Acid-polyacrylamide gel electrophoresis of wheat glutenins: A new tool for the separation of high and low molecular weight subunits. *Cereal Chem.* **1994**, *71*, 238–242.

25. Singh, N.K.; Shepherd, K.W.; Cornish, G.B. Rapid communication a simplified SDS-PAGE procedure for separating LMW subunits of glutenin. *J. Cereal Sci.* **1991**, *14*, 203–208. [CrossRef]

26. Payne, P.I.; Lawrence, G.J. Catalogue of alleles for the eomplex gene loci, GLU-Al, GLU-BI, and GLU-DI which code for high-molecular-weight subunits of glutenin in hexaploid wheat. *Cereal Res. Commun.* **1983**, *11*, 29–35.

27. R Core Team. *A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2014.

28. Weegels, P.L.; Hamer, R.J.; Schofield, J.D. Critical review: Functional properties of wheat glutenin. *J. Cereal Sci.* **1996**, *23*, 1–18. [CrossRef]

29. Triboi, E.; Abad, A.; Michelena, A.; Lloveras, J.; Ollier, J.; Daniel, C. Environmental effects on the quality of two wheat genotypes: 1. quantitative and qualitative variation of storage proteins. *Eur. J. Agron.* **2000**, *13*, 47–64. [CrossRef]

30. Martre, P.; Jamieson, P.D.; Semenov, M.A.; Zyskowski, R.F.; Porter, J.R.; Triboi, E. Modelling protein content and composition in relation to crop nitrogen dynamics for wheat. *Eur. J. Agron.* **2006**, *25*, 138–154. [CrossRef]
31. Doekes, G.J.; Wennekes, M.J. Effect of nitrogen fertilization on quantity and composition of wheat flour protein. *Cereal Chem.* **1982**, *59*, 276–278.

32. Gupta, R.B.; Batey, I.L.; MacRitchie, F. Relationships between protein and functional properties of wheat flours. *Cereal Chem.* **1992**, *69*, 125–131.

33. Jia, Y.Q.; Masbou, V.; Aussenac, T.; Fabre, J.L.; Debaeke, P. Effects of nitrogen fertilization on maturation conditions on protein aggregates and on the breadmaking quality of Soissons, a common wheat cultivar. *Cereal Chem.* **1996**, *73*, 123–130.

34. Pechanek, U.; Karger, A.; Gröger, S.; Charvat, B.; Schöggl, G.; Lelley, T. Effect of Nitrogen Fertilization on Quantity of Flour Protein Components, Dough Properties, and Breadmaking Quality of Wheat. *Cereal Chem. J.* **1997**, *74*, 800–805. [CrossRef]

35. Johansson, E.; Prieto-Linde, M.L.; Svensson, G. Influence of nitrogen application rate and timing on grain protein composition and gluten strength in Swedish wheat cultivars. *J. Plant Nutr. Soil Sci.* **2004**, *167*, 345–350. [CrossRef]

36. Johansson, E.; Kuktaite, R.; Andersson, A.; Prieto-Linde, M.L. Protein polymer build-up during wheat grain development: Influences of temperature and nitrogen timing. *J. Sci. Food Agric.* **2005**, *85*, 473–479. [CrossRef]

37. Randall, P.J.; Wrigley, C.W. Effects of sulfur supply on the yield, composition and quality of grain from cereals, oilseeds, and legumes. *Adv. Cereal Sci. Technol.* **1986**, *8*, 171–206.

38. Dupont, F.; Hurkman, W.; Vensel, W.; Chan, R.; Lopez, R.; Tanaka, C.; Altenbach, S. Differential accumulation of sulfur-rich and sulfur-poor wheat flour proteins is affected by temperature and mineral nutrition during grain development. *J. Cereal Sci.* **2006**, *44*, 101–112. [CrossRef]

39. Gagliardi, A.; Carucci, F.; Masci, S.; Flagella, Z.; Della Gatta, G.; Giuliani, M.M. Effects of Genotype, Growing Season and Nitrogen Level on Gluten Protein Assembly of Durum Wheat Grown under Mediterranean Conditions. *Agronomy* **2020**, *10*, 755. [CrossRef]

40. Ferrise, R.; Bindi, M.; Martre, P. Grain filling duration and glutenin polymerization under variable nitrogen supply and environmental conditions for durum wheat. *Field Crop. Res.* **2015**, *171*, 23–31. [CrossRef]

41. Daniel, C.; Triboï, E. Effects of temperature and nitrogen nutrition on the accumulation of gliadins analysed by RP-HPLC. *Funct. Plant Biol.* **2001**, *28*, 1197. [CrossRef]

42. Malik, A.H.; Kuktaite, R.; Andersson, A.; Johansson, E.; Prieto-Linde, M.L. Individual and interactive effects of cultivar maturation time, nitrogen regime and temperature level on accumulation of wheat grain proteins. *J. Sci. Food Agric.* **2011**, *91*, 2192–2200. [CrossRef]

43. Jenner, C.; Ugalde, T.; Aspinall, D. The Physiology of Starch and Protein Deposition in the Endosperm of Wheat. *Funct. Plant Biol.* **1991**, *18*, 211–226. [CrossRef]

44. Wieser, H.; Seilmeier, W.; Belitz, H.-D. Quantitative Determination of Gliadin Subgroups from Different Wheat Cultivars. *J. Cereal Sci.* **1994**, *19*, 149–155. [CrossRef]

45. Uthayakumaran, S.; Tömösközi, S.; Tatham, A.S.; Savage, A.W.J.; Gianibelli, M.C.; Stoddard, F.L.; Bekes, F. Effects of Gliadin Fractions on Functional Properties of Wheat Dough Depending on Molecular Size and Hydrophobicity. *Cereal Chem. J.* **2001**, *78*, 138–141. [CrossRef]

46. Zhao, F.; Hawkesford, M.; McGrath, S. Sulphur Assimilation and Effects on Yield and Quality of Wheat. *J. Cereal Sci.* **1999**, *30*, 1–17. [CrossRef]

47. Zorb, C.; Steinfurth, D.; Seling, S.; Langenkämper, G.; Koehler, P.; Wieser, H.; Lindhauer, M.G.; Muhling, K.H. Quantitative Protein Composition and Baking Quality of Winter Wheat as Affected by Late Sulfur Fertilization. *J. Agric. Food Chem.* **2009**, *57*, 3877–3885. [CrossRef] [PubMed]

48. Blumenthal, C.; Barlow, E.; Wrigley, C. Growth Environment and Wheat Quality: The Effect of Heat Stress on Dough Properties and Gluten Proteins. *J. Cereal Sci.* **1993**, *18*, 3–21. [CrossRef]

49. Zhang, P.; He, Z.H.; Zhang, Y.; Xiangchun, Z.; Chen, D.; Zhang, Y. Association between % SDS-Unextractable Polymeric Protein (%UPP) and End-Use Quality in Chinese Bread Wheat Cultivars. *Cereal Chem. J.* **2008**, *85*, 696–700. [CrossRef]

50. Wrigley, C.W.; Békés, F.; Bushuk, W. *Gliadin and Glutenin: The Unique Balance of Wheat Quality*; Chapter 1, Gluten: A Balance of Gliadin and Glutenin; AACC: Washington, DC, USA, 2006; pp. 3–32.
51. Ammar, K.; Kronstad, W.E.; Morris, C.F. Breadmaking Quality of Selected Durum Wheat Genotypes and Its Relationship with High Molecular Weight Glutenin Subunits Allelic Variation and Gluten Protein Polymeric Composition. *Cereal Chem. J.* **2000**, *77*, 230–236. [CrossRef]

52. Boggini, G.; Pogna, N. The breadmaking quality and storage protein composition of Italian durum wheat. *J. Cereal Sci.* **1989**, *9*, 131–138. [CrossRef]

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).