Grain Quality of Maize Cultivars as a Function of Planting Dates, Irrigation and Nitrogen Stress: A Case Study from Semiarid Conditions of Iran

: Maize grain is an important source of human and animal feed, and its quality can be affected by management practices and climatic conditions. This study aimed to evaluate the concentration and composition of starch, protein and oil in grain of maize cultivars in response to different planting dates (20 June and 21 July), irrigation (12-day and 6-day intervals) and nitrogen rates (0 and 184 kg N ha − 1 ). The ﬁrst two principal components (PCs) accounted for 84.5% of the total variation. High N fertilization increased protein (by 6.0 and 10.9 g kg − 1 ) and total nonessential amino acids (by 3.4 and 2.4 g kg − 1 ) during 2018 and 2019, respectively. With the high irrigation rate, the high N rate increased oil, total unsaturated fatty acids, and starch and amylopectin, whereas with the low irrigation rate, there was no effect of the N rate. With earlier planting, total saturated fatty acids were higher. The ﬁndings highlight the complicated relationship between the different factors and how they affect quality characteristics of maize grain. There was a large impact of year, which to a great extent cannot be controlled, even in this environment where water supply was controlled and rainfall did not affect the results.


Introduction
Dent maize (Zea mays L. var. indentata) is cultivated because of its various end uses including livestock nutrition, human consumption and ethanol production [1].
The mature grain of dent maize is composed of storage components including 60% to 72% starch [2], 8 to 11% protein [3] and 4 to 6% oil [1]. Amylose and amylopectin are glucose storage polymers that can alter physicochemical properties of starch [4]. The optimal amylose-amylopectin ratio differs depending on the purpose of use. While highamylose maize grain is a good source to supply resistant starch for food industries [5], highamylopectin grain can be appropriate for livestock feeding because it is easier to digest [6]. For poultry diets, amylose content is significantly correlated with lower digestibility or higher resistant starch content because of a very compact physical structure [7].
Protein of cereals is critically important with regards to nutritional quality. Protein is relevant in terms of its overall concentration, and the essential and nonessential amino acids and their ratios [8]. Diets with essential amino acids as the only N source are used less efficiently than diets with an optimal ratio of essential to nonessential amino acids [9], which affects the amount of N losses from poultry excrement [10]. Saturated and unsaturated fatty acids and their distribution in triacylglycerol molecules are also important compositional factors of crude fat or oil to estimate metabolizable energy values in animal diets and oil nutritional value for human consumption [11].
The abovementioned compositional factors, as substantial dimensions of cereal grain quality are strongly affected by the genetic potential, the growing environment and agricultural practices [12,13]. Among agronomic practices, key strategies include cultivar selection [14]; water regimes [15]; nutritional status, especially N utilization [8]; and edaphic and climatic conditions, in particular, temperature changes due to year, location and planting date [16]. The effect of these factors on grain quality is more complex than yield. As an example, apart from interspecies differences, reports have shown that higher N rates reduced rice cooking quality by decreasing amylose [17], and in wheat grain, higher N doses caused an increase in different fractions of proteins and amino acids [18]. Lower temperatures resulted in smaller starch granule formation and decreased amylose content of maize grain starch [19]. Concentrations of amino acids increased with canopy warming in winter wheat [20]. Soil water deficit decreased the contents of amylose and starch in wheat grain [21] and of oil and linolenic acid in maize grain [22]. Water stress caused an increase in amino acid concentrations of maize grain [23]. An increase in starch [11] and a decrease in protein concentration [24] were reported for sorghum grain grown under high water-deficit conditions. Interactions between factors can also affect grain quality. For example, variations in nitrogen-based compounds like protein and amino acids can be a consequence of the dilution effect-decreasing the concentration of elements in plant tissue due to changes in environmental conditions-in response to irrigation and N interactions [25]. As a result, cultivar, water regime, N rate and planting date impact the qualitative characteristics of cereal grains; however, more knowledge is needed on how these factors work in an interactive way.
Thus, in this study, we assessed the interactive effects of planting date, irrigation timing and N fertilization on various measures of maize grain quality. The aim of the research was to provide a basis of knowledge for defining future research and for developing guidelines for improving maize grain quality.

Experiment Site
Field experiments were performed in 2018 and 2019 in Pakdasht city (35.4669 • N, 51.6861 • E), Tehran province, Iran. The geographical location of Pakdasht city in Iran is shown in Figure 1. This region has a semiarid climate with relatively cold winters and hot summers. Annual precipitation in Pakdasht is approximately 160 mm and concentrated in late autumn and winter. Irrigated summer maize in this region is typically cultivated after winter cereals like wheat and barley. Summaries of the climatic parameters during the summer maize growing season for two years of the experiment are presented in Table 1. The overall means of daily temperatures were 25.3 and 24.9 • C for 2018 and 2019, respectively, and the total precipitation was less than 1.5 mm in each growing season. The highest temperatures occurred in July for both years. Soil sampling was conducted before the first cultivation year from two depths (0-30 and 30-60 cm). The samples were air dried; passed through a 2-mm sieve; and tested for organic carbon using the Walkley-Black method [26], total nitrogen (N) using the Kjeldahl method [27], for available phosphorus (P) by the Olsen procedure [28] and for available Potassium (K) by using a flame photometer [29]. The soil is characterized as clay-loam, and the initial physicochemical characteristics in 0−30 cm soil depth were as follows: organic carbon, 1.09%; total N, 0.12%; P, 90.17 mg kg −1 ; and K, 453 mg kg −1 . Moisture release curve was used to determine soil water [30,31]. Water content at field capacity and permanent wilting point were 0.346 and 0.115 g·cm −3 , respectively.

Experimental Design
In mid-June 2018, experimental strips and plots were designated at a farm where no crops had been cultivated for the previous three years. The field was divided into two strips (24-m length) to apply different irrigation timings (12-day and 6-day intervals). Two planting dates (20 June and 21 July) were randomized within the irrigation treatments. Two ditches crossed the strips, partitioning each of them into three plots. Thus, the experiment comprised twelve main plots. The four combinations of cultivars (KSC704 and KSC260) and N rates (0 and 184 kg ha −1 ) were randomized to subplots within each plot. Thus, a total of 12 main plots containing a total of 48 subplots were obtained from three replications in each cultivation year. Individual subplots were 6 m in length and consisted of 6 rows sown at a density of six and eight plants per square meter for KSC704 and KSC260, respectively. Treatments were applied to the same plots in each year.
KSC704 is a locally popular high-yield late maturity hybrid which has been cultivated in Iran since 1980 and accounts for 80% of the maize growing areas of Iran [32]. KSC260, an early maturity maize cultivar, was introduced in 2008 and has shown good

Experimental Design
In mid-June 2018, experimental strips and plots were designated at a farm where no crops had been cultivated for the previous three years. The field was divided into two strips (24-m length) to apply different irrigation timings (12-day and 6-day intervals). Two planting dates (20 June and 21 July) were randomized within the irrigation treatments. Two ditches crossed the strips, partitioning each of them into three plots. Thus, the experiment comprised twelve main plots. The four combinations of cultivars (KSC704 and KSC260) and N rates (0 and 184 kg ha −1 ) were randomized to subplots within each plot. Thus, a total of 12 main plots containing a total of 48 subplots were obtained from three replications in each cultivation year. Individual subplots were 6 m in length and consisted of 6 rows sown at a density of six and eight plants per square meter for KSC704 and KSC260, respectively. Treatments were applied to the same plots in each year.
KSC704 is a locally popular high-yield late maturity hybrid which has been cultivated in Iran since 1980 and accounts for 80% of the maize growing areas of Iran [32]. KSC260, an early maturity maize cultivar, was introduced in 2008 and has shown good performance in various experiments comparing new early-maturing cultivars [33]. Harvesting date, Agriculture 2021, 11, 11 4 of 16 maturity period and the total growing degree days (GDDs) of cultivars in the present study are shown in Table 2. GDD was calculated using the GDD calculator program [34]. For this study, 10 • C was set as the base temperature, 34 • C was set as the optimum temperature and 40 • C was set as the maximum temperature threshold, adapted from Cutforth and Shaykewich [35]. The two planting dates were chosen to create different temperature conditions during grain filling. The variability of daily average temperature during the reproductive stage (VT-R6) for both planting dates is shown in Figure 2. There were clear differences between the recorded temperatures for the planting dates in each year and differences between the two years. The second year was a warmer year during the reproductive period. The temperature thresholds used to calculate GDDs were 10 • C (base temperature), 34 • C (optimum) and 40 • C (maximum) (Cutforth and Shaykewich, 1990).
Agriculture 2021, 11, x FOR PEER REVIEW 4 of 17 performance in various experiments comparing new early-maturing cultivars [33]. Harvesting date, maturity period and the total growing degree days (GDDs) of cultivars in the present study are shown in Table 2. GDD was calculated using the GDD calculator program [34]. For this study, 10 °C was set as the base temperature, 34 °C was set as the optimum temperature and 40 °C was set as the maximum temperature threshold, adapted from Cutforth and Shaykewich [35]. The two planting dates were chosen to create different temperature conditions during grain filling. The variability of daily average temperature during the reproductive stage (VT-R6) for both planting dates is shown in Figure 2. There were clear differences between the recorded temperatures for the planting dates in each year and differences between the two years. The second year was a warmer year during the reproductive period. The temperature thresholds used to calculate GDDs were 10 °C (base temperature), 34 °C (optimum) and 40 °C (maximum) (Cutforth and Shaykewich, 1990). An irrigation regime with 6-day intervals and N rate of 184 kg ha −1 were specified as typical non-stressed growing conditions. Furrow irrigation was used for irrigation. The closed-end furrows were constructed with a ditcher. The ridge and furrow widths were 40 cm and 20 cm, respectively, and the depth of the furrow was 15 cm. Based on the flow rate per furrow (approximately 1.3 L/s), water was delivered until the furrows were completely full and the duration of each irrigation event was recorded. The known flow rate of the irrigation pump was then used to estimate irrigation volumes.
In total, the KSC704 cultivar received about 9200 m 3 ha −1 of water with 6-day irrigation timing and about 5800 m 3 ha −1 with 12-day irrigation timing. KSC260 received less water due to a shorter growing period and less irrigation events: 7800 m 3 ha −1 with 6-day An irrigation regime with 6-day intervals and N rate of 184 kg ha −1 were specified as typical non-stressed growing conditions. Furrow irrigation was used for irrigation. The closed-end furrows were constructed with a ditcher. The ridge and furrow widths were 40 cm and 20 cm, respectively, and the depth of the furrow was 15 cm. Based on the flow rate per furrow (approximately 1.3 L/s), water was delivered until the furrows were completely full and the duration of each irrigation event was recorded. The known flow rate of the irrigation pump was then used to estimate irrigation volumes.
In total, the KSC704 cultivar received about 9200 m 3 ha −1 of water with 6-day irrigation timing and about 5800 m 3 ha −1 with 12-day irrigation timing. KSC260 received Agriculture 2021, 11, 11 5 of 16 less water due to a shorter growing period and less irrigation events: 7800 m 3 ha −1 with 6-day irrigation timing and 4800 m 3 ha −1 with 12-day irrigation timing. The Food and Agriculture Organization (FAO) estimates the water requirement of late-maturing maize cultivars in arid and semiarid regions to be around 7000 to 8000 m 3 ha −1 [36]. Urea was used as the N source and top-dressed in equal proportions at two stages (pre-planting and V4-V6). Fresh irrigation water was used from the main ditch for each replication, and excess water was not reused.

Sample Preparation
At the R6 (physiological maturity) stage, based on the Hanway standard [37], 15 plants in an area of approximately 11 m 2 (4-m row length × 0.7-m row spacing × 4 rows) were randomly cut from each subplot. Ears were threshed and dried at 60 • C to a constant weight. Seeds were separated from ears and were weighed and ground. The milled samples were maintained at 5 ± 1 • C until the start of laboratory procedures.

Analysis of Starch Content and Composition
Total starch and its composition (amylose and amylopectin) were determined using a Spectrophotometer (V-M5 model, BEL Engineering, Monza, Italy) set at 510 nm according to the amylose/amylopectin Megazyme procedure [38].

Analysis of Protein Content and Composition
Crude protein, essential amino acids (composed of methionine, threonine, valine, isoleucine, leucine, phenylalanine, histidine, arginine, lysine and tryptophan) and nonessential amino acids (composed of cysteine, asparagine, serine, glutamine, proline, glycine and alanine) were estimated using Near-Infrared Reflectance Spectroscopy (NIRS) (NIRS-XDS model, Foss, Hilleroed, Denmark) in the 1100-2500 nm wavelength range at five-nm intervals. A dent maize library of global origin samples that previously had been assayed by wet chemistry or non-NIR methods was used for model calibration.

Analysis of Oil Concentration and Fatty Acid Composition
Oil extraction was conducted using a Soxhlet extractor for approximately 4 h with hexane as the solvent, with a solid to solvent ratio of 1:7 m v −1 . The fatty acid methyl esters (FAME) were extracted according to AOAC 996.06 protocol [39]. The synthesized FAME was injected into gas chromatography (GC) (CP-Sil 88 model, Varian, Walnut Creek, CA, USA) to detect fatty acid composition by curve and retention time. The GC was equipped with a flame ionization detector (FID) (Column: CP-Sil 88 (100 m * 250 µm * 0.2 µm) removable phase: Nitrogen, 28.8 min, heat injection chamber: 270 • C, heat detector: 260 • C).

Statistical Analysis
Principal component analysis (PCA) scores were derived using Minitab, v.19 [40] after standardizing the variables by using the correlation matrix. Results from the two cultivation years were analyzed separately, using mixed model procedures with PROC MIXED in SAS, v.9.4 software [41]. The mixed-effects model included fixed effects of irrigation, planting date, cultivar, N, their interactions, and the random effects of irrigation × planting date and irrigation × planting date × block interactions. Because of the experimental design, irrigation effects have to be interpreted with caution-possible meaningful effects of irrigation and its interactions can be either caused by differences between the planned irrigation treatments or by differences in the two sections of the field where the irrigation treatments were applied. Tukey's statistic was used to test differences (p ≤ 0.05) among means.

Correlations between Variables and Treatments
The PCA comprising the first two principal components accounted for 84.5% of the total variance. The loading plot ( Figure 3A) shows the eigenvectors for nine variables: starch, amylose, amylopectin, protein, total essential amino acids (ΣEAA), total nonessential amino acids (ΣNEAA), oil, total saturated fatty acids (ΣSFA) and total unsaturated fatty acids (ΣUSFA).

Correlations between Variables and Treatments
The PCA comprising the first two principal components accounted for 84.5% of the total variance. The loading plot ( Figure 3A) shows the eigenvectors for nine variables: starch, amylose, amylopectin, protein, total essential amino acids (ΣEAA), total nonessential amino acids (ΣNEAA), oil, total saturated fatty acids (ΣSFA) and total unsaturated fatty acids (ΣUSFA). year and irrigation treatment. Each marker represents the mean of replicates in the field. Abbreviations; I1, Irrigation at 12-day intervals; I2, Irrigation at 6 day-intervals; ΣEAA, total essential amino acids; ΣNEAA, total non-essential amino acids; ΣSFA, total saturated fatty acids; ΣUSFA, total unsaturated fatty acids.
The cosine of the angle between two vectors estimates the correlation between them; therefore, clustered points are highly correlated with each other. There are two clusters of variables that are strongly correlated with each. The first cluster includes protein concentration, protein composition (ΣEAA and ΣNEAA) and starch composition (amylose concentration) of the grain. The eigenvector for amylopectin points in the opposite direction; thus, the first cluster is highly negatively correlated with amylopectin (e.g., the correlation between ΣNEAA and amylopectin is −0.917). Running perpendicular to the first cluster, the second cluster of highly correlated variables includes oil, starch and ΣUSFA. This cluster is highly negatively correlated with ΣSFA (e.g., the correlation between starch and ΣSFA is −0.663).
The score plot means are shown in Figure 3B. To illustrate the most noticeable trends, points are labelled according to combinations of year and irrigation treatment. The first component separates the first year of the experiment from the second year-markers for the first year (2018) are located in the left of the score plot, whereas markers for the second year (2019) are in the right. The second component separates well-irrigated samples from those under potential water-stress. Markers for well-watered samples are distributed toward the top, while markers from water-limited plots are mostly located toward the bottom. However, the differences of irrigation timing are more clearly separated in the first year of the experiment than second. shows the means, grouped by year and irrigation treatment. Each marker represents the mean of replicates in the field. Abbreviations; I1, Irrigation at 12-day intervals; I2, Irrigation at 6 day-intervals; ΣEAA, total essential amino acids; ΣNEAA, total non-essential amino acids; ΣSFA, total saturated fatty acids; ΣUSFA, total unsaturated fatty acids.
The cosine of the angle between two vectors estimates the correlation between them; therefore, clustered points are highly correlated with each other. There are two clusters of variables that are strongly correlated with each. The first cluster includes protein concentration, protein composition (ΣEAA and ΣNEAA) and starch composition (amylose concentration) of the grain. The eigenvector for amylopectin points in the opposite direction; thus, the first cluster is highly negatively correlated with amylopectin (e.g., the correlation between ΣNEAA and amylopectin is −0.917). Running perpendicular to the first cluster, the second cluster of highly correlated variables includes oil, starch and ΣUSFA. This cluster is highly negatively correlated with ΣSFA (e.g., the correlation between starch and ΣSFA is −0.663).
The score plot means are shown in Figure 3B. However, the differences of irrigation timing are more clearly separated in the first year of the experiment than second.
The data points presented in score plots coincide with the directions of change in maize grain compositional variables in the loading plots. The two clusters of eigenvectors and their negative correlations define two axes, making an x shape. The first axis (from the southwest to northeast quadrants) separates 2018 from 2019. i.e., 2019 samples were higher in amylose, ΣEAA, ΣNEAA and protein. The second axis (from the northwest to southeast quadrants) separates the irrigation treatments. Irrigated plots were higher in starch, oil and ΣUSFA.

Treatment Effects on Starch Content and Composition
After calculating p-values of all main effects and interactions (Supplementary Materials S1), the treatment main effect means are presented in Table 3 and significant interactions were plotted (unless no pairs of means were significantly different or the effect could not be clearly interpreted). Table 3. Least square means, significances and standard errors of maize grain starch, amylopectin, amylose, protein, total nonessential amino acids (ΣNEAA), total essential amino acids (ΣEAA), oil, total unsaturated fatty acids (ΣUFA) and total saturated fatty acids (ΣSFA) in response to treatment main effects (irrigation, planting date, cultivar and nitrogen rate) in 2018 and 2019. In 2019, average starch concentration in KSC704 was higher than KSC260 by 11 g kg −1 ( Table 3). In 2018, there was an interactive effect of irrigation and N rate on starch ( Figure 4A). With the high irrigation rate, the high N rate increased starch, whereas with low irrigation, the high N rate decreased starch. Grain starch was higher with the high irrigation rate, regardless of the N rate. In 2019, high N increased amylopectin concentration by 17 g kg −1 (Table 3). In 2018, there was an interactive effect of irrigation and N rate on amylopectin ( Figure 4B). With the high irrigation rate, the high N rate increased amylopectin (649 g kg −1 ), whereas with low irrigation, there was no N effect. Grain amylopectin was higher with the high irrigation rate, regardless of the N rate.
In 2018, there were interactive effects of irrigation and N rate ( Figure 4C); irrigation rate and cultivar ( Figure 4D); and N rate, planting date and cultivar ( Figure 4E) on amylose concentration. For either irrigation rate, high N decreased amylose ( Figure 4C). Amylose was higher with the high irrigation rate (with either N rate) than with the low irrigation rate ( Figure 4C). With the high irrigation rate, there was no cultivar effect, whereas with low irrigation, KSC260 (103 g kg −1 ) had higher amylose than KSC704 (97 g kg −1 ) ( Figure 4D). Grain amylose was higher with the high irrigation rate regardless of the cultivar ( Figure 4D). With the zero N rate, KSC704 planted late had higher amylose than that planted earlier.
With the high N rate, KSC260 planted late had higher amylose than that planted earlier. In general, amylose was higher with zero N fertilizer than with the high N rate ( Figure 4E).  In 2019, there were interactive effects of N rate, cultivar and irrigation rate ( Figure 4F); planting date and cultivar ( Figure 4G); and planting date and N rate ( Figure 4H) on amylose concentration. Low irrigation rate and zero N rate in KS704 had the greatest amylose concentration (205.15 g kg −1 ), whereas high irrigation and N rate in KSC260 decreased amylose to its lowest value (154.53 g kg −1 ) ( Figure 4F). With an early planting date, KSC704 had higher amylose than KSC260 (198.98 g kg −1 ), whereas with a late planting date, there was no cultivar effect ( Figure 4G). For both planting dates, N application decreased amylose concentration; however, the effect was greater with the late planting date ( Figure 4H).

Treatment Effects on Protein Content and Composition
Irrigation, planting date and cultivar did not influence protein concentration in either year. However, high N increased proteins by 6.0 g kg −1 in 2018 and 10.9 g kg −1 in 2019. Moreover, high N increased ΣNEAA concentrations by 3.4 g kg −1 in 2018 and 2.4 g kg −1 in 2019 (Table 3). In 2018, low irrigation rate increased ΣNEAA by 3.2 g kg −1 , and KSC260 had higher ΣNEAA than KSC704 by 1.2 g kg −1 (Table 3). In 2019, there was an interactive effect of planting date and cultivar on ΣNEAA ( Figure 5A). With the early planting date, KSC704 had a higher ΣNEAA than KSC260, whereas with late planting, KSC260 had a higher ΣNEAA than KSC704.
Agriculture 2021, 11, x. https://doi.org/10.3390/xxxxx www.mdpi.com/journal/agriculture on ΣEAA concentration. For either irrigation rate, a high N rate increased the ΣEAA concentration ( Figure 5B); however, the effect was greater with a low irrigation rate. For either planting date, the high N rate increased the ΣEAA concentration ( Figure 5C); however, high N had a greater effect on ΣEAA concentration (32.35 g kg −1 ) for the later planting date. In 2019, there was an interactive effect of planting date and cultivar on ΣEAA ( Figure  5D). With the early planting date, KSC704 had higher ΣEAA than KSC260, whereas with the late planting date, there was no cultivar effect.

Treatment Effects on Oil Content and Composition
Planting date and cultivar did not affect oil and ΣUSFA concentration in either year. In both years, there was an interactive effect of irrigation and N rate on oil ( Figure 6A,B) and ΣUSFA ( Figure 6C,D). In 2018, with the high irrigation rate, there was no N effect on oil and ΣUSFA, whereas with low irrigation, the high N rate decreased oil and ΣUSFA ( Figure 6A,C). In 2019, with the high irrigation rate, the high N rate increased oil and ΣUSFA, whereas with low irrigation, there was no N effect ( Figure 6B,D). In both years, oil and ΣUSFA were higher in the high irrigated plots (with either N rate) than with the low irrigated ones. In 2019, ΣSFA concentration with the early planting date was higher In 2018, KSC260 had higher ΣEAA than KSC704 by 0.8 g kg −1 (Table 3). In 2019, high N increased the ΣEAA concentration by 1.4 g kg −1 (Table 3). In addition, in 2018, there was an interactive effect of irrigation and N ( Figure 5B), and planting date and N ( Figure 5C) on ΣEAA concentration. For either irrigation rate, a high N rate increased the ΣEAA concentration ( Figure 5B); however, the effect was greater with a low irrigation rate. For either planting date, the high N rate increased the ΣEAA concentration ( Figure 5C); however, high N had a greater effect on ΣEAA concentration (32.35 g kg −1 ) for the later planting date. In 2019, there was an interactive effect of planting date and cultivar on ΣEAA ( Figure 5D). With the early planting date, KSC704 had higher ΣEAA than KSC260, whereas with the late planting date, there was no cultivar effect.

Treatment Effects on Oil Content and Composition
Planting date and cultivar did not affect oil and ΣUSFA concentration in either year. In both years, there was an interactive effect of irrigation and N rate on oil ( Figure 6A,B) and ΣUSFA (Figure 6C,D). In 2018, with the high irrigation rate, there was no N effect on oil and ΣUSFA, whereas with low irrigation, the high N rate decreased oil and ΣUSFA ( Figure 6A,C). In 2019, with the high irrigation rate, the high N rate increased oil and ΣUSFA, whereas with low irrigation, there was no N effect ( Figure 6B,D). In both years, oil and ΣUSFA were higher in the high irrigated plots (with either N rate) than with the low irrigated ones. In 2019, ΣSFA concentration with the early planting date was higher (by 0.22 g kg −1 ) than with the late planting date (Table 3). In both years, there was an interactive effect of irrigation and N rate on ΣSFA ( Figure 6E,F). In 2018, with the high irrigation rate, there was no N effect on ΣSFA, whereas with low irrigation, the high N rate decreased ΣSFA ( Figure 6E). In 2019, with the high irrigation rate, the high N rate increased ΣSFA, whereas with low irrigation rate, the high N rate decreased ΣSFA ( Figure 6F).
(by 0.22 g kg −1 ) than with the late planting date (Table 3). In both years, there was an interactive effect of irrigation and N rate on ΣSFA (Figure 6E,F). In 2018, with the high irrigation rate, there was no N effect on ΣSFA, whereas with low irrigation, the high N rate decreased ΣSFA ( Figure 6E). In 2019, with the high irrigation rate, the high N rate increased ΣSFA, whereas with low irrigation rate, the high N rate decreased ΣSFA ( Figure  6F).

Correlations between Variables and Treatments
The results suggest that year and irrigation were the most influential effects in this experiment. The findings are well supported by Hammac et al. [42], who reported that temperature and water changes are more effective than soil nutrient status for changing rapeseed composition. In this study, there are two main potential sources for the differences among the trend of the years: (a) daily temperature differences among years and (b) performing the experiment after three years of fallow. The grain filling period is an enzyme-dependent stage of accumulating storage materials, primarily starch and protein and is sensitive to factors affecting photosynthesis, especially temperature and soil moisture, and nutritional status [43].
Many studies point to the negative impact of elevated temperatures and water deficit on oil, starch and dry matter accumulation in cereal grain [44,45]. Examples from the literature report 76.8% of the changes in oat oil content [46] and 52% of the changes in canola oil content [47] explained by climatic variation, especially water and temperature variability among years. However, Riccardi et al. [48] found that water stress induces the expression of proteins not specifically related to this stress but rather to reactions against cell damage. This may be the reason why the protein and its composition did not show a clear response to further irrigation compared to starch and oil, according to the PCA results. Castro et al. [49] reported that slight heat stress increased protein of wheat grain by shortening the grain filling period and by increasing the rate of N remobilization to grain. This is a potential reason for the long eigenvector of protein concentration in the second year.

Correlations between Variables and Treatments
The results suggest that year and irrigation were the most influential effects in this experiment. The findings are well supported by Hammac et al. [42], who reported that temperature and water changes are more effective than soil nutrient status for changing rapeseed composition. In this study, there are two main potential sources for the differences among the trend of the years: (a) daily temperature differences among years and (b) performing the experiment after three years of fallow. The grain filling period is an enzymedependent stage of accumulating storage materials, primarily starch and protein and is sensitive to factors affecting photosynthesis, especially temperature and soil moisture, and nutritional status [43].
Many studies point to the negative impact of elevated temperatures and water deficit on oil, starch and dry matter accumulation in cereal grain [44,45]. Examples from the literature report 76.8% of the changes in oat oil content [46] and 52% of the changes in canola oil content [47] explained by climatic variation, especially water and temperature variability among years. However, Riccardi et al. [48] found that water stress induces the expression of proteins not specifically related to this stress but rather to reactions against cell damage. This may be the reason why the protein and its composition did not show a clear response to further irrigation compared to starch and oil, according to the PCA results. Castro et al. [49] reported that slight heat stress increased protein of wheat grain by shortening the grain filling period and by increasing the rate of N remobilization to grain. This is a potential reason for the long eigenvector of protein concentration in the second year.
The great impact of irrigation amount on the starch, protein and oil content of maize grain was reported by Kresović et al. [15]. Moreover, other studies point to the large effect of water availability [50] and genotype-environment interactions [51] on compositional attributes of maize and wheat grain. On the potential effect of fallow, fallows improve soil fertility, organic matter and physical properties to supply essential nutrients needed for assimilation and changes in grain composition. In contrast, continuous cultivation may result in some elements being deficient in grain, if they are not provided through fertilization [52].

Treatment Effects on Starch Content and Composition
The results showed that the average starch concentration for KSC704 was higher than for KSC260. One potential reason for this result is genetic differences between the two cultivars [53]. There is a positive correlation between grain weight and starch content in cereals [54,55]. Considering that the 100-grain weight in KSC704, as a late maturity hybrid, is higher than KSC260 [53], a higher starch concentration is expected. The highest concentration of starch was achieved by interactive effects of higher N and irrigation rate. Starch accumulation in grains is a physiological process of transportation and conversion of photosynthetic assimilates into starch and is expected to increase under irrigation and adequate N supply [56]. Moreover, the results showed that high N application with low irrigation decreased starch concentration. One potential reason for this result is the inverse starch and protein relationship in grain with high N application [57].
In the first year, the highest value of amylopectin was achieved by the interactive effect of high irrigation and N rate. In the second year, a high N rate increased amylopectin concentration. These results are in agreement with Jiyun et al. [58] and Kaplan et al. [59], who reported that amylopectin contents of maize hybrids increased under high irrigation and N rate in a similar trend to starch. The effect of N on amylose in the present experiment was that high N application was associated with reduced grain amylose. This is in close agreement with Kaufman et al. [60], who reported that high N application reduces type A granules in sorghum grain compared to type B or C, [61], suggesting that high N application reduces amylose by decreasing type A granules.
In addition, the results of the first year suggest that amylose was higher with the high irrigation rate regardless of N rate. The results are not in agreement with Kaplan et al. [59], who reported that the amylose content of maize grain decreased with an increase in irrigation rate, and it increased with an increase in N rate. Potential reasons for the discrepancy between results are different N and irrigation rates, and climatic conditions. The second-year results suggest that planting KSC704 earlier significantly increased amylose concentration. Few studies have investigated the simultaneous effects of maturity group and planting date on amylose content of grain. However, the increase in amylose in response to earlier planting [19] and lengthening maize kernel maturation [62,63] has been suggested.

Treatment Effects on Protein Content and Composition
Higher protein content in plots with higher N rates was not surprising. These results are similar to the results obtained by Saint Pierre et al. [64], Yang et al. [65], and Cao et al. [17], who reported there is a positive correlation between N application and protein content of cereal grains. This is because N stimulates the activity of panel enzymes involved in protein biosynthesis [66]. Increasing the crude protein concentration in grain can be achieved via two scenarios: (a) increasing N utilization and (b) sustaining higher partitioning of N to grain (nitrogen harvest index-NHI) [67]. Due to the nonsignificant differences between cultivars in terms of protein, it can be concluded that NHI did not differ between cultivars.
Our findings on the significant effect of N, irrigation and their interaction on ΣNEAA and ΣEAA are in agreement with Zhang et al. [8], who reported that ΣEAA in wheat grain increased in response to high N rate and low frequency irrigation. They also reported an interactive effect of irrigation and N rate on wheat grain amino acids but did not have a consistent response during the three years of their experiment. Since N is one of the basic element of amino acids and protein compounds in grain [68], higher amino acids in N-contained treatments is expected. In the first year, ΣNEAA was lower in plots with more frequent irrigation. One potential reason for a decrease in amino acids in high irrigation conditions is due to yield dilution effects on N-containing compounds in the grain [8].
The first-year results showed that KSC260 had higher ΣNEAA and ΣEAA than KSC704. In the second year, amino acids differed between cultivars and changing planting dates. In general, the inconsistency of cultivars amino acids among years and planting dates could suggest a large impact of exogenous factors in addition to the maturity group and genetic potentials on grain quality [69], which can further be associated with the complexities of genotype-environment-management interactions. The results are in close agreement with Huang et al. [70] in rice, who found that, with earlier planting dates, early maturity cultivars enter the reproductive phase earlier than late maturity ones, potentially causing flowering to coincide with high summer temperatures and consequently reducing amino acids in early maturity cultivars. Accordingly, we infer that the absence of ΣEAA differences between maturity classes in the late planting date may be because of lower temperatures during grain filling.

Treatment Effects on Oil Content and Composition
The results demonstrated that, for any year of cultivation, oil concentrations were affected by interaction effects of irrigation and N rate. This observation gives weight to the results from Aguirrezábal et al. [71] and Kaplan et al. [22], who reported that the interaction effects of irrigation and N rate are very influential on grain oil. In any year of cultivation, grain oil was higher with the high irrigation rate, irrespective of N rate. This implies that water was more critical than N for increasing oil concentration and quality. This is because water deficiency reduces grain oil by decreasing N uptake [72] and germ growth and by reducing the enzyme activity responsible for lipid biosynthesis [73].
The results suggest that, for any year of cultivation, ΣUSFA was higher with high irrigation, regardless of N rate. Our findings for irrigation effects on fatty acids were somewhat different with Kaplan et al. [22]. These authors reported that high irrigation decreased linoleic acid, as the most abundant unsaturated fatty acids in maize grain. Higher irrigation regimes, lower air temperature and higher precipitation during grain filling are potential reasons for decreasing linoleic acid under higher irrigation rates in the study conducted by Kaplan et al. [22].
The results showed a negative impact of high N rate on ΣSFA under the low irrigation rate. This is in close agreement with the study of Ali and Ullah [74], who reported that a high N rate (225 kg ha −1 ) decreased palmitic acid and stearic acid in sunflower hybrids. In the second year, ΣSFA concentration was higher with earlier planting. Similar results were reported by Obeng et al. [75] in camelina cultivars. A possible reason for the difference between planting dates in the second year is a higher air temperature during the reproductive stage of maize growth. Some studies proposed that, as the average daily temperature rises during grain filling, the crop tends to produce more saturated fatty acids in sunflower [76] and oilseed crops [77]. This may be related to high-temperature impacts on lipid profiles by destabilizing enzymes effective in unsaturated fatty acid synthesis; as a result, saturated fatty acids increase in the grain [78].
There are some commonalities between trends from PCA and ANOVA results; however, one possible reason for some discrepancy between ANOVA and PCA trends is that the PCA plots only present the data in two dimensions (principle components 1 versus 2), whereas the correlation matrix values take into account all dimensions [79].

Conclusions
This study aimed to provide insights into understanding the relationship among quality characteristics of maize grain. Applying principal components analysis, the first two PCs accounted for 84.5% of the total variation. Year and irrigation had the greatest effect on yield and quality. Plots with high irrigation were associated with higher starch, oil and total unsaturated fatty acids (ΣUSFA), and data points from 2019, a warmer year, were associated with higher amylose, protein and amino acids. Analysis of variance results revealed more details on the effects of other factors and their interactions on maize grain components. In any year of cultivation, N application significantly increased protein and ΣNEAA values. A combination of high irrigation and N rate often increased oil and fatty acids values, whereas with the low irrigation rate, increased N had no effect. The cultivar KSC704 had a higher starch concentration, and KSC260 had a higher amino acids concentration. With earlier planting, ΣSFA was higher. The study was limited to two cultivars, and although there were clear differences between them, further studies that include additional cultivars would provide more confidence in the results. However, year-to-year variations in the effects of factors on amylose, amylopectin and amino acids suggest that the response of cultivars to the environment plays an important role in the final composition of starch and protein. The findings highlight the complicated relationship between the experimental factors and the large impacts of growing season conditions on quality attributes of maize grain.