Bread Wheat Quality under Limiting Environmental Conditions: I-Molecular Properties of Storage Proteins and Starch Constituents in Mature Grains

: To support the introduction of local soft wheat varieties, the Lebanese state has implemented, during recent years, an improvement program to select plant material with good productivity in semi-arid conditions that can lead to national production that can, at the same time, meet quality needs expressed by Lebanese processors. In the present study, the main biochemical components of grains conventionally associated with bread-making properties (storage proteins, starch polysaccharides) were physiochemically characterized for a panel of local soft wheat varieties selected. During the two consecutive years of cultivation, the post-ﬂowering thermal constraints signiﬁcantly modiﬁed the kinetics of accumulation of the major constituents by mainly limiting the grain ﬁlling time. The level of polymerization/aggregation (i.e., gliadin/glutenin ratio, unextractable polymeric protein (UPP) amount, molecular weight average of glutenins) of prolamins was very high. The reduction in the amount of starch was accompanied by a signiﬁcant change in the amylopectin/amylose ratio. Finally, the genotypes studied were characterized by signiﬁcantly different distributions of starch granules; the percentage of the volume occupied by A-type and B-type starch granules varied between genotypes for the two cropping years. All these observations must be considered because of their determining role in the technological aptitude of the ﬂours generated.


Introduction
Wheat count among the most important raw materials for human consumption. It is well known for its nutritional qualities. In Lebanon, this cereal is grown mainly in the Beqaa plain, considered as the country's agricultural basket. Annual production in 2019 reached a level of 140,000 tons [1], with 75% made up of durum wheat (Triticum durum), mainly due to its better adaptation to environmental conditions, resulting in acceptable levels of yield appreciated by farmers. At the same time, annual imports in 2019 amounted to 1,575,000 tons [1], mainly constituting soft wheat suitable for breading. Today, in order to limit this level of imports, the Lebanese state wishes to support the introduction of local soft wheat varieties. Thus, an improvement program has been put in place to select plant material adapted to environmental conditions (i.e., good productivity in semi-arid conditions) that can lead to national production of a good quality (i.e., respect the needs expressed by Lebanese processors (i.e., millers and bakers)).
From these works, the variation of the baking properties of wheat flour clearly depends on the ability of its storage protein fractions (i.e., prolamins) to establish, during mixing, a three-dimensional structure with viscoelasticity characteristics [16]. Consequently, the storage proteins found in common wheat flour show a strong tendency to auto-aggregate into a three-dimensional plastic network during the kneading of dough. Since the 2000s, particular attention has been paid to glutenin proteins because a close relationship/correlation has been established between baking properties (i.e., extensibility, mixing time, and loaf volume) and the molecular weight distribution (MWD) of the aggregates of polymeric proteins [17].
Beyond these important protein parameters, it is essential to consider modifications of the content and/or characteristics of starch components contained in the grains [18][19][20][21]. Indeed, several authors have previously pointed to the influence of different amylaceous characteristics (i.e., amylose/amylopectin ratio [22], distribution of starch granule sizes [23], grain surface-to-volume ratio [24]) on grain processing aptitudes.
Our work is part of a general study on the assessment of processing skills (mainly bread-making performance) of local soft wheat varieties selected by LARI's (Lebanese Agricultural Research Institute) Grain and Vegetable Research Department in collaboration with the International Center for Agricultural Research in the Dry Areas (ICARDA) and the Arab Center for the Studies of Arid Zones and Drylands (ACSAD) for their agronomic behavior to handle environmental constraints. The first part of this general study (part I: present paper) is dedicated to the characterization of the main biochemical components selected (i.e., molecular properties of storage proteins and starch of matured grains) and a second part, which is in press [9], is dedicated to the characterization of the rheological behavior of flours and doughs produced.

Plant Materials
Four different bread wheat genotypes, selected for their potential tolerance/resistance to hydric stress and rust [25][26][27] by two agricultural research centers operating in dry zones (ICARDA and ACSAD, Bekaa Valley, Lebanon), are delivered by LARI (Bekaa Valley, Lebanon) to be the material of this study. These cultivar pedigrees are shown in Table 1. Table 1. Parentage and pedigree of the bread wheat genotypes studied.

Wheat Genotypes
Line/Pedigree Source Bread wheat varieties were cultivated for 2 successive years (2015 and 2016) at Tal Amara Lebanon in LARI experimental farm, located at 905 m altitude, 33 • 28 N latitude, and 36 • 30 E longitude. The four cultivars were grown in a sandy loam clay soil with three replications according to a randomized complete block design. Experiments were performed without irrigation intervention. Seeding took place in early December, whereas harvesting was done in late June. Central Beqaa climate is semi-arid, wherein rain is concentrated between December and April.

Wheat Grain Quality Assessment
Studied traits of wheat homogeneous grain samples were based on the determination of the thousand kernels weight (TKW), the protein total content (N × 5.7) by Dumas method (AACC 76.13.01) on LECO FP-528 (Perkin Elmer, Villebon sur Yvette, France), and starch content [28].

Wheat Flour Preparation
Wheat grains were dampened at 16% (w/w) humidity and tempered for 24 h before being milled on an experimental Bühler miller (type MLU-202, Bühler, Villepinte, France) according to (AACC 26-21.02) [28]. Resulting white flour, having a mean extraction rate of 65% and a mean ash content of 0.68%, was subjected to different analysis with 3 replicates of each sample.

Quantification of Grain Proteins by Size Exclusion-High Performance Liquid Chromatography (SE-HPLC)
For the extraction of Sodium Dodecyl Sulphate (SDS)-extractable proteins, 10 mg of white flour samples were mixed with sodium phosphate buffer solution (1 mL, pH 6.90, 0.05 M, SDS (2%, w/v)). After 2 h of incubation at 60 • C with continuous stirring, mixtures were centrifuged at 12,500× g for 30 min at 20 • C. Supernatants containing the SDS-soluble proteins were filtrated through a syringe filter of regenerated cellulose (porosity = 0.45 µm), and 20 µL was then injected into the SE-HPLC system. The pellets were resuspended in of the same sodium phosphate solution (1 mL, pH 6.90, 0.05 M, SDS (2%, w/v)) and sonicated for 20 s using a 3 mm microtip probe (Sonics Materials, Thermo Fisher Scientific, Les Ulis, France, model 75038) and centrifuged for 30 min (12,500× g at 20 • C). The supernatants containing the SDS-unextractable proteins were filtered through 0.45 µm followed by injection of 20 µL.
For the total protein extraction, 10 mg of white flour was incubated at 60 • C for 2 h with a continuous stirring in the presence of sodium phosphate buffer solution (1 mL, pH 6.90, 0.05 M, SDS (2%, w/v)). Sonication was then performed for 20 s using a 3 mm microtip probe (Sonics Materials, Thermo Fisher Scientific, Les Ulis, France model 75038). After this, samples were centrifuged at 12,500× g for 30 min at 20 • C. As above, the supernatants were filtered through 0.45 µm, and 20 µL was then injected into the column.

Determination of the Molecular Distribution of Flour Protein by Asymmetrical Flow Field Flow Fractionation (A4F)
Wheat protein molecular weight distribution was performed as reported previously [17]. Flour (30 mg) was mixed with sodium phosphate solution (1 mL, pH 6.90, 0.05 M, SDS (2%, w/v)) and incubated for 2 h at 60 • C with vortexing every 10 min. Samples were then sonicated (20 s at power setting ≈30%) with a 3 mm microtip probe before being centrifuged (12,500× g) for 15 min at 25 • C and filtered through 0.45 µm and then injected (30 µL) in A4F/MALLS system. The A4F machine used was an Eclipse3 F System (Wyatt Technology, Santa Barbara, CA, USA) combined with a multi-angle light scattering (MALLS) detector (Dawn multi-angle Heleos TM, Wyatt Technology, Santa Barbara, CA, USA) and an Optilab T-rEX refractive index detector (Wyatt Technology, Santa Barbara, CA, USA). In addition, The Agilent HPLC 1200 Series (Agilent Technologies, Les Ulis, France) was used in tandem with the AF4 system. The Trapezoidal channel had 286 mm of length and the used spacer was 350 µm. The ultrafiltration wall consisted of regenerated cellulose membrane with a cutoff of 10 kDa. The mobile phase consisting of a sodium phosphate buffer solution (pH 6.90, 0.05 M, SDS (0.1%, w/v)) passed continually through a 0.1 µm regenerated cellulose filter (Merck Millipore, Guyancourt, France). Absorbance was recorded at 214 nm. The cross flow used for fractionation was focused at 0.5 min at a constant flow rate of 2 mL·min −1 . Injection step was realized during 1.0 min at a flow rate of 0.2 mL·min −1 , and a relaxation phase was adjusted to 0.5 min. These steps were followed by elution at flow detector of 1.0 mL·min −1 and cross-flow rate decreasing from 3.0 to 0.0 mL·min −1 . After 14 min of elution, cross-flow rate was maintained at 0.0 mL·min −1 for 9 min. Determination of number-average molar mass (Mn), weight-average (Mw) molar mass, radius of gyration (Rz), and hydrodynamic radius of gyration (Rh) was performed using ASTRA 7.1.2 software (Wyatt Technology, Santa Barbara, CA, USA).

Quantitation of High Molecular Weight Glutenin Subunits (HMW-GS) by Lab-on-a-Chip
Extraction and quantification of glutenin subunits was achieved as reported in our previous study [30]. Wheat flour (30 mg) was suspended for 20 min at 22 • C in a buffer solution (1 mL, 0.08 M Tris-HCL, pH 7.50, propanol-1 (50%, v/v)). Samples were centrifuged for 15 min (15,900× g). Monomeric proteins (i.e., gliadins, albumins, and globulins) contained in the supernatant were discarded. Pellets containing polymeric proteins were mixed with 0.6 mL of a Tris-HCl solution with SDS (2%, w/v) and DTT (1%, w/v) and sonicated for 20 s (30%) with a microtip probe of 3 mm. The mixes were incubated for 30 min at 60 • C before being centrifuged for 15 min (12,500× g at 22 • C). The separation of glutenin subunits was performed on a Lab-on-a-Chip apparatus (LabChip GXII Touch HT Protein Characterization System, PerkinElmer, Waltham, MA, USA), where the chip channel was primed with a polymer solution of polydimethyl methacrylate of high molecular weight (pDMA) in a solution of Tris-Tricine containing SDS and a non-covalent coloration of 0.27% (w/v).

Extraction of Starch Granules of the Wheat Flours and Determination of Their Size
Starch granule extraction was based on the developed INRA (Institut National de la Recherche Agronomique) protocol by Bancel et al., 2010 [31]. A sample of 500 mg of flour was added to 10 mL of ultrapure water, mixed and filtered (100 µm), and then centrifuged at 4000× g for 5 min at 22 • C. The precipitated material was twice treated in the same way. A total of 5 mL of Tris-HCL buffer (55 mM, pH 6.80, SDS 2.3% (w/v), 10% glycerol (v/v), and 1% thiothreitol (w/v)) were then added to the precipitate and vigorously stirred for 30 min. Samples were then sonicated (20 s, 30% at power max) (3 mm microtip probe) and centrifuged (4000× g) for 5-10 min at 20 • C. Pellets were washed for 3 consecutive times using ultra-pure water (10 mL) and then washed twice with acetone (10 mL). Centrifugation for 5-10 min (4000× g, 20 • C) was performed after each washing step. Dried starch was kept at 4 • C.
Starch granule size distribution was performed through laser diffraction using Malvern Masterizer 2000 (Malvern, Palaiseau, France). Distribution of starch granules is expressed by the percentage of the volume occupied by granules due to their size. It is a particle size analyzer that sizes starch particles through a wet dispersion unit. The measurements were conducted using a MU hydro dispersion unit. The stirring energy ensured the homogenization of the sample suspension. The chosen pump speed was 2500 revolution per minute (rpm) in order to obtain a maximum homogenization of the suspension in the beaker by elimination of air bubbles. A particular detector registers the intensity of laser light in the optical cell that is converted to particles size distribution.

The Molecular Weight Distribution of Starch Polysaccharides by Asymmetrical Flow Field Flow Fractionation (A4F)
The molar mass distribution of starch polysaccharides was performed according to the developed procedure reported by Chiaramonte et al. [32]. Purified starch (10 mg) was dissolved in dimethyl sulfoxide (1 mL, DMSO 95% (w/v)) at 100 • C for 1 h, then precipitated in 5 mL of ethanol and centrifuged (20,000× g) for 20 min (20 • C). The precipitate was added to 4 mL of NaOH (20 mM) and solubilized for 8 min at 135 • C in a microwave (Discover CEM, Saclay, France), then filtered through a 0.45 µm syringe filter (Merck Millipore, Guyancourt, France) and injected (100 µL) in the A4F/MALLS system. The A4F machine used was an Eclipse3 F System (Wyatt Technology, Toulouse, France), which can be similarly described as previously (Section 2.5). However, the used channel was the short channel that had a 195 mm length and the used spacer was of 0.35 mm. The system was normalized and checked with pullulan 110 kDa and 800 kDa, respectively. Detection of molecules was obtained by a differential refractometer (DRI) and through multiangular light scattering MALLS.

Statistical Analysis
Data statistical analysis was performed using XLSTAT 2020 (Addinsoft, Paris, France). Multiple mean comparison by ANOVA multifactor was calculated through a variance analysis-type ANOVA (p = 0.05). Multiple mean comparisons were performed using Tukey's test HSD (honestly significant difference) (p = 0.05).

Characteristics of the Grain Filling Period
Climatic changes impact on ecophysiology, and potentials of wheat cultures essentially depend on the temperature and water. In semi-arid zones, such as in Lebanon, mean day temperatures and maximal day temperature are of higher levels [33]. The frequency of day temperature level exceeding 30 • C is significantly elevated. Thus, in these particular zones, high temperatures can lead to higher evapotranspiration and the apparition of a hydric deficit, resulting in a lower photosynthesis activity. These environmental conditions can have important consequences if they occur after the anthesis during the phases of cell division and/or the cell enlargement in wheat kernels.
In our study, the grain filling period of the different wheat genotypes studied can be characterized by the sum of degree day (ΣDD) registered during two different seasons of culture (i.e., 2015 and 2016) from anthesis up to the harvest maturity of the grains. Figure 1 shows the data recorded in terms of (i) sum of degree day (ΣDD) but also in terms of (ii) grain filling period duration for the four different wheat genotypes studied. The sums of degree day observed were fairly homogeneous; they varied between 677 DD (mean value for 2015) and 679 DD (mean value for 2016). KATILA and ACSAD 1133 are two genotypes characterized by higher precociousness, and this is regardless of the cultivation year concerned. At the same time, the grain filling period length observed for the same wheat genotypes varied between 34 days after anthesis (DAA) (mean value for 2016) and 39 DAA (mean value for 2015). Whatever the wheat genotype studied, these values correspond to a severe limitation of the accumulation time length of the grain reserves (storage proteins and starch) classically observed (i.e., more than 10 or 20 days of limitation) [34]. Thus, the higher the daily temperature, the number of days of accumulation decreases. For the two cultivation seasons, the temperatures and the precipitation recorded dur- For the two cultivation seasons, the temperatures and the precipitation recorded during the same period are presented in Figure 2.   (1) Multiple mean comparisons were made using Tukey's test (HSD-honestly significant difference) and different letters indicate a statistically difference (p < 0.05). (2) Mean values of two consecutive cultivation years and three blocks.  These first physiological observations seem to confirm that the different soft wheat genotypes used during our study were subjected to environmental conditions that can be qualified as "limiting" physiological conditions. These environmental conditions have led to a natural limitation of the length of the grain filling phase that, in the absence of modification of the filling speed of the grains [2], has had, as is shown later, a significant impact on the intensity of the accumulation of total dry matter in the grains.
The general characteristics of wheat grains (i.e., TKW, protein content, starch content) harvested in 2015 and in 2016 are presented in Table 2. Thus, the TKW values were between 21.42 and 25.96 g, with a very low multi-year average value corresponding to 24.5 g (from 23.35 to 25.69 g for 2016 and 2015, respectively). In fact, TKW values observed in our study were significantly reduced (i.e., from −25% to −35%) compared to those previously obtained by El-Kareem and El-Saidy [35], who mentioned mean values of about 34.2 g for the genotypes parented of the wheat cultivars studied here in classical environmental conditions in Egypt. In our study, SHAM 8 was the only wheat genotype that differentiated statistically (p = 0.05), with an average value of 21.42 g. These results, which reflect a strong limitation of the dry matter accumulation in the grains for all the genotypes studied, seemed to be mainly related to a significant dysfunction of the post-anthesis carbon metabolism (i.e., the biosynthesis and the accumulation of the starch resulting from post-anthesis photosynthetic activity). Indeed, as shown in the results presented in Table 2, the grains of the different genotypes studied were characterized by very high levels of total proteins (average protein content of 17.43 g/100 g DM (dry matter) with annual averages between 14.78 and 20.29 for 2016 and 2015, respectively) and low levels of starch (average starch content of 55.16 g/100 g DM with annual averages between 53.57 and 56.75 for 2015 and 2016, respectively). In contrast to the metabolism of protein accumulation, which is mainly is a remobilization metabolism, the starch accumulation metabolism essentially corresponded to a de novo synthesis metabolism from the post-anthesis photosynthetic activity of the terminal leaf surfaces (i.e., mainly the last two foliage stages) [36]. In fact, any reduction in photosynthetic activity will result in a deficit of biosynthesis and starch accumulation in the grains and, consequently, a lack of dilution of the accumulated protein content, resulting in the consequent reduction of TKW [37]. These results are in total accordance with some previous observations made in different environmental limiting conditions [38,39] and, in the case of our study, these limiting conditions were even more marked in 2015 than in 2016.

Prolamin Content of Wheat Grains and Their Molecular Weight Distribution
All the different results of these assays are compiled in Table 2. Although there were slight significant differences between the four genotypes studied (i.e., between TAL AMARA 2 and the other wheat genotypes), the polymeric prolamin contents were very close (from 29.92% to 36.51% for KATILA and TAL AMARA 2, respectively). At the same time, the content of gliadin (monomeric prolamins) represented between 50.20% and 55.47% of the total protein content of wheat grains, with a multi-year average value corresponding to 54.06%.
In view of these results, it appears that the four genotypes in the context of our study were characterized by relatively high gliadin contents. In fact, this fraction of prolamins represented here more than 54% of the total protein content, i.e., approximately +10 to +15% more than the protein content commonly encountered for soft wheat varieties [40]. Gliadins were generally synthesized and accumulated most rapidly between 10 and 30 DAA. On the contrary, glutenin subunits [Low molecular weight glutenin subunits (LMW-GS) and particularly high molecular weight glutenin subunits (HMW-GS)] accumulated more slowly but for a longer period during grain filling.
As was shown above, the environmental conditions that were recorded during our study led to a significant reduction in the filling time of the grains. As a consequence, these limiting environmental conditions caused a significant reduction in the synthesis and accumulation of polymeric prolamins in favor of the monomeric prolamins.
Since the 2000s, scientists agreed that the molecular weight distribution (MWD) of prolamin proteins are the main factors determining the technological properties of flours [41]. However, variations in MWD that may exist between wheat genotypes can theoretically result from modification of the ratio between monomeric proteins and polymeric proteins (gliadin to glutenin ratio) but also from modification of size distribution of polymeric proteins alone [42] determined by quantification of the non-extractable polymeric protein (UPP) and/or by analysis of the glutenin fraction using A4F.
In our study, whatever the wheat genotype studied, the gliadin/glutenin ratio (also called the prolamin polymerization index) was herein characterized by high values (i.e., >1.2-1.3) ( Table 2). Indeed, this ratio varied between 1.38 and 1.89 for TAL AMARA 2 and KATILA, respectively, with a general mean value of 1.72. TAL AMARA 2 is the only wheat genotype that differentiates statistically (p = 0.05) with an average value of 1.38. At the same time, UPP contents were quite comparable between genotypes, with values between 41.35% and 46.70% for TAL AMARA 2 and SHAM 8, respectively. On average (i.e., 45.00%), these high values corresponded to highly aggregated/polymerized glutenin contents. Table 3 presents the macromolecular features (Mn, Mw, Rz, and Rh) of polymeric prolamins of the wheat grains. As can be seen from the calculation of the polydispersity index (Ip = Mw/Mn), all the polymeric fractions present in the different wheat grains were characterized by a high heterogeneity (i.e., Ip > 5.0). Indeed, the values collected varied between 6.96 and 9.41 for TAL AMARA 2 and SHAM 8, respectively. The genotypes studied differed significantly (p = 0.05) in two statistically distinct groups. (1) Mn: molecular weight number-average, Mw: molecular weight weight-average, Ip: polydispersity index (Mw/Mn); (2) Rz: radii of gyration, and Rh: hydrodynamic radii of gyration; (3) multiple mean comparisons made using Tukey's test (HSD) and different letters indicate a statistically difference (p < 0.05); (4) mean values of two consecutive years and three blocks.
The number-average molar mass (Mn) and the weight-average molar mass (Mw) respectively varied from 0.55 × 10 6 g/mol to 0.61 × 10 6 g/mol, and from 3.92 × 10 6 g/mol to 5.13 × 10 6 g/mol. These two macromolecular features allowed us to discriminate some genotypes studied (mainly SHAM8 and TAL AMARA 2), with the calculated average molecular weights being characteristic of highly polymerized/aggregated prolamins, thus confirming the results obtained previously by quantifying the UPP fraction. In the same way, radius of gyration Rz and hydrodynamic radius Rh respectively ranged from 85.70 nm to 94.16 nm and from 8.01 nm to 8.23 nm. Table 4 shows, for the first time to our knowledge, the high molecular weight glutenin subunits (HMW-GS) composition of the polymeric prolamins of each of the selected genotypes and the HMW-GS/LMW-GS (Low molecular weight glutenin subunits) ratio. (1) Glutenin subunit content: relative percentage expressed in relation to the total content of HMW-GS (high molecular weight glutenin subunits); (2) multiple mean comparisons made using Tukey's test (HSD) and different letters indicate a statistically difference (p < 0.05); (3) mean values of two consecutive years and three blocks. LMW-low molecular weight.
The interactions that may exist between polymeric proteins were strongly related to the nature of the HMW-GS (especially HMW-GS pair 5 + 10 vs. HMW-GS pair 2 + 12 coded by Glu-D1) and the HMW-GS/LMW-GS ratio [41]. Thus, any modification of HMW-GS composition within the polymeric prolamins will cause a modification of the MWD of these polymers. As seen in SHAM and TAL AMARA 2, for the same allelic combination of HMW-GS, any improvement in the HMW-GS/LMW-GS ratio (Table 4) resulted in an increase in the weight-average molecular weight (Mw) ( Table 3) and vice versa.

Starch Content of Wheat Grains and Their Molecular Weight Distribution
Starch consists of two structurally different polymers: (i) amylose (AML), a linear polymer of (1,4) α-linked glucose units, and (ii) amylopectin (AMP), a highly branched structure of glucose units with (1,6) and (1,4) α-linkages [43]. Amylopectin and amylose are present in a ratio of about 3:1 in most plant starches, and analysis of the mature grains in our study showed about 34.22% amylose (Table 5). Our results can be considered important; however, this quantification is very dependent on the methodology used. As has already been demonstrated [32], the separation and quantification of starch polymers by A4F makes it possible to obtain more reliable and reproducible results. In any case, this content seems very stable (i.e., no significant differences between the different genotypes during the two years of culture).
The two constituent polymers of starch (AML and AMP) can be characterized by different macromolecular features (Mn, Mw, and Rz). As shown by the results of these measurements (Table 5), the amylopectin polymers (AMP), which were the most abundant polymers, were characterized by very high molecular weights: 16.1 × 10 6 g/mol and 58.8 × 10 6 g/mol for Mn the Mw, respectively. (1) Mn = molecular weight number-average (×10 6 g/mol), Mw = molecular weight-average (×10 6 g/mol), Rz = radii of gyration (nm), V = slope of the log-log plot of Rz vs. Mw; (2) multiple mean comparisons made using Tukey's test (HSD) and different letters indicate a statistically difference (p < 0.05); (3) mean values of two consecutive years and three blocks.
For the same AMP content, TAL AMARA 2 was the only genotype to differ in having significantly higher molecular weights than the other three genotypes studied (i.e., 24.0 × 10 6 g/mol and 66.2 × 10 6 g/mol for Mn and Mw, respectively). At the same time, amylose polymers (AML) were characterized by more limited molecular weights: 0.66 × 10 6 g/mol and 0.33 × 10 6 g/mol for Mn and Mw, respectively. As in the previous case, for the same AML content, TAL AMARA 2 was the only genotype to differ in having significantly higher molecular weights than the other three genotypes studied (i.e., 0.80 × 10 6 g/mol and 0.37 × 10 6 g/mol for Mn and Mw, respectively).
For macromolecular assemblies such as starch polymers, the shape or the branching rate was most reliably determined from the slope V of the log-log plot of Rz vs. Mw. The calculated values of V for the AMP polymers (≈0.32) corresponded to branched macromolecular assemblies. As can be seen from the results presented in Table 5 even if the calculated data for V were close, TAL AMARA 2 differentiated significantly (p = 0.05) from the other wheat genotypes retained, with a value of V = 0.29. Thus, under the conditions of our study, for a quantity of stable AMP comparable to that of the other genotypes (i.e., no statistical difference in AMP%), TAL AMARA 2 presented polymers of AMP whose branching rate was higher (i.e., significant decrease of V).

Starch Granule Distribution
Starch presents as granules in the wheat grain endosperm. These granules have been reported to have trimodal size distributions [44]. The biosynthesis of the A-type starch granules (generally with diameters > 10 nm) begins in the first days of endosperm formation. The formation of B-type starch granules (with diameters < 10 nm) starts from about 10 DAA. Finally, the synthesis of C-type starch granules (less than 5 nm) initiates from approximately 21 DAA [45]. Table 5 shows the distribution of starch granules in the grains of the four wheat genotypes studied. As shown in our results, A-type were present in a very large majority (76.3%), and B-type and C-type represented 18.8% and 4.0%, respectively. Our observations were in total agreement with previous results obtained from the majority of common wheat genotypes. Under the physiological conditions of our study (i.e., "limiting" conditions, §3.1), certain genotypes differed statistically (p = 0.05). Thus, SHAM 8 was characterized by a population of granules poorer in A-type starch granules (72.8% vs. 76.3% on average) and richer in B-type starch granules (22.8% vs. 18.8% on average). Conversely, KATILA was characterized by a starch granule distribution with a large majority of A-type starch granules (79.2% vs. 76.3% on average) and a lower proportion in B-type starch granules (16.8% vs. 18.8% on average).
The modifications of the size and/or the number of the different types of starch granules in grains can be a characteristic of wheat genotype studied, or it can be caused by the modification of environmental conditions (particularly growth temperature) during the grain-filling period. The volume percentage of A-type and B-type starch granules can be modified when growth temperature increases from about 15 to 40 • C during the grain storage accumulation period [46]. These environmental effects could induce changes in the grain-filling pattern of starch granules (A-and B-type granules specifically). A-type granules are formed in the amyloplast from approximately 5 DAA and continue to grow until achieving a maximum diameter at physiological maturity [43]. However, the final number of A-amyloplasts is reached at about 7 DAA, coinciding with the cessation of cell division. On the other hand, B-type granules are initiated at about 11 DAA and continue to increase until 21 DAA, and up to a maximum diameter at maturity (35 DAA). Therefore, considering the behavior of type A and type B granules, we argue that any increase in temperature during the cell enlargement phase could induce a reduction in the activity of starch synthase but also a reduction in the grain filling phase (Figure 1 and Table 2), resulting in significant changes in the size and number of different types of starch granules in the endosperm. All of these elements must be considered because, due to their significant differences in terms of composition, molecular structure, granule swelling, gelatinization properties, and pasting/rheological behavior, proportion of the different type of starch granules (particularly A-and B-type) in flour has an impact on the properties of traditional wheat-based products [46].

Conclusions
Within the general framework of the implementation of a program to improve the production of local common wheat by the Lebanese State, we characterized the plant material selected for its adaptation to environmental conditions (i.e., good productivity in semi-arid conditions) by initially retaining all the biochemical components of the grains that are able to guarantee the desired processing performance (i.e., mainly baking performance).
Despite the selection of the plant material used in this study, our results demonstrate that environmental conditions induce limitations in its physiological functioning during the important phases of grain development and maturation. Thus, the post-flowering thermal constraints characteristic of semi-arid environments significantly modify the kinetics of accumulation of the major constituents of the grain by mainly limiting the grain filling time. The wheat grains thus formed were characterized by reduced TKW resulting from a deficit in starch accumulation.
At the same time, the synthesized and accumulated prolamins can be characterized by their particular molecular distributions. Whatever plant material selected, the level of polymerization/aggregation of polymeric prolamins was very high. The reduction in the amounts of synthesized and accumulated starch polysaccharides was accompanied by a significant change in the amylopectin/amylose ratio, with amylose content being greater than normal (i.e., >34%). Finally, the different genotypes studied were characterized by significantly different distributions of starch granules, with the percentage of the volume occupied by A-type and B-type starch granules varying between genotypes for the two cropping years. Here, again, we confirm with our results that the specific thermal regime during the synthesis and accumulation phases of prolamins and starch can explain the polymerization/aggregation changes of prolamins but also the modifications in the distribution of granules starch (A-type vs. B-type).
All these observations (i.e., molecular distribution of prolamins and starch, distribution of starch granules) must be considered because of their determining role for the definition of the technological aptitudes of the flours generated. These relationships between composition and rheological properties will be discussed in a future publication.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.