Avenanthramides: Unique Bioactive Substances of Oat Grain in the Context of Cultivar, Cropping System, Weather Conditions and Other Grain Parameters

Our study was focused on the evaluation of the content of a wider spectrum of eight avenanthramides (AVNs) as unique components of oat grain under the effects of four selected factors (cultivar, locality, cropping system, and year). The weather effects on changes in the AVN content and their relationship to other important parameters of oat grain were further evaluated in more detail. A sensitive UHPLC system coupled with a QExactive Orbitrap mass spectrometer was used for AVN quantification. AVNs confirmed a high variability (RDS = 72.7–113.5%), which was dominantly influenced by the locality and year factors. While most AVN types confirmed mutually high correlations (r = 0.7–0.9), their correlations with the other 10 grain parameters were lower (r ≤ 0.48). Their significant correlations (0.27–0.46) with β-D-glucan could be used in perspective in breeding programs for the synergetic increase of both parameters. PCA analysis and Spearman correlations based on individual cultivars confirmed a significant effect of June and July precipitation on the increase of Σ AVNs. However, the results also indicated that higher precipitation can generate favorable conditions for related factors, such as preharvest sprouting evoking a direct increase of AVNs synthesis in oat grain.


Introduction
Oat (Avena sativa L.), with a current production share of about 2.3% of all cereals in the European Union, can be considered as a minor cereal compared to wheat or barley production [1]. The genus Avena includes about 70 species. The most important cultivated oat belongs to the hexaploid species Avena sativa L. (2n = 6x = 42; AACCDD). Within this species, two different types of cultivars were cultivated: oat with hulled and oat with naked grain [2]. The commercial oat cultivars with naked grain are referred to in some studies [3,4] as a separate species of Avena nuda L. Nevertheless, the species A. nuda belongs taxonomically to diploid oats which are not commonly cultivated [2].
The nutritionally high-value composition of oat grain has found its application in both livestock fattening and human nutrition. Oat consumption in the human diet has been increased because of health benefits associated with advantageous composition in macronutrients: (i) lipids with a high degree of unsaturation, including oleic and linoleic acids (about 40% and 36% of total fatty acids, respectively), (ii) proteins with a favorable  The three major AVNs (AVN-2p, AVN-2f, and AVN-2c) accounted for about 2/3 of the share of all eight AVNs detected. On average, AVN-2f (47.2 mg/kg of dw) had the highest values in the tested group. In contrast, the lowest mean value was confirmed for AVN-3p (1.0 mg/kg of dw).
Similar relative variability (RSD = 74.8%) as for the Σ AVN was confirmed only in the content of immunoreactive avenin peptides G12. The achieved average value of the tested set (11.2 mg/kg of dw) was significantly below the permitted value of 20 mg/kg for coeliac patients. However, in some individual cases (e.g., 'Patrik' 2018-locality CB, 'Patrik' 2019-locality PR, and 'Seldon' 2019-locality PR), this limit was slightly exceeded (25-40 mg/kg of dw). The total average values of other monitored parameters, including grain yield and its qualitative parameters (TGW, content of crude proteins, fat, starch, β-D-glucan, and ash), varied to a much lesser extent (from RSD = 4.3% for starch to RSD = 40.3% for yield; see Table 1). In summary, the values for most parameters were typical of oats. On average, higher values were mainly detected in the protein content (CPmax = 20.1% and both protein fractions, AVEmax = 3.5%; GLUmax = 3.5%). A summary of all detected parameters depending on the cultivar, locality, growing system, and year is further presented in Supplementary Tables S1 and S2.

Analysis of the Influence of the Main Factors on the Monitored Grain Parameters
The significance of the four individual factors (cultivar, cropping system, locality, and year) was calculated using a four-way analysis of variance (ANOVA). The variability of all AVNs and most other parameters was significantly affected by all four factors (Tables 2  and 3). The only exceptions were the three parameters-ash, G12, and β-D-glucan, where the cropping system was not statistically significant, and in the case of β-D-glucan, the locality was also insignificant (Table 3).    Statistical analysis of the effects of the main factors on the content of AVNs (ANOVA) and also a subsequent graph describing their mutual interactions are illustrated in Table 2 and Figure 1. The analysis confirmed a higher influence of locality and year on AVN variability compared to cultivar and the cropping system. In particular, the CB locality generally showed higher contents of Σ AVN compared to PR (158.6 mg/kg of dw vs. 60.8 mg/kg of dw). In the conventional conditions of 2019 and 2020, the increase of Σ AVNs compared to the PR locality was even 3-4 times. On the contrary, an example of a grain parameter with a high influence of genotypes (cultivar) on its variability was evident in the content of β-D-glucan, presented in Figure  2.  On the contrary, an example of a grain parameter with a high influence of genotypes (cultivar) on its variability was evident in the content of β-D-glucan, presented in Figure 2. The three oat cultivars confirmed similar average values of Σ A of dw), which were significantly higher compared with ´Korok´ ( 'Raven' (106.3 mg/kg of dw) cultivars. Significant differences betwe AVNs were also found in the above three cultivars with the highest AVNs ( Table 2). The results further showed a significantly higher c the case of conventional cultivation (Table 3). This difference was m tected high content of Σ AVN in cultivars at the conventional CB loc This is especially true for the 'Seldon' cultivar, which in 2020 show increase in the content of AVNs (407 mg/kg) in conventional conditi A direct quantitative comparison of the percentage influence of the variability of the tested parameters (the share of the main signific sum of squares) is illustrated in Figure 3. The results confirmed tha vidual AVNs were most affected by the locality (26-38%) and then The only exception was AVN-2fd, with a lower total influence of The three oat cultivars confirmed similar average values of Σ AVNs (116-123 mg/kg of dw), which were significantly higher compared with Korok (84.9 mg/kg of dw) or 'Raven' (106.3 mg/kg of dw) cultivars. Significant differences between individual types of AVNs were also found in the above three cultivars with the highest content of the sum of AVNs ( Table 2). The results further showed a significantly higher content of Σ AVNs in the case of conventional cultivation (Table 3). This difference was mainly due to the detected high content of Σ AVN in cultivars at the conventional CB locality in 2019 and 2020. This is especially true for the 'Seldon' cultivar, which in 2020 showed an extremely high increase in the content of AVNs (407 mg/kg) in conventional conditions of the CB locality.
A direct quantitative comparison of the percentage influence of significant factors on the variability of the tested parameters (the share of the main significant factors in the total sum of squares) is illustrated in Figure 3. The results confirmed that the contents of individual AVNs were most affected by the locality (26-38%) and then by the year (9-24%). The only exception was AVN-2fd, with a lower total influence of all four main factors (only 31.5%) and the highest influence of the year (14.8%). The influence of the other two factors on the variability of individual AVNs was significantly lower, with ranges of 1-11% for the cultivar and 1-3% for the cropping system, respectively. For Σ AVNs, the effects of the cultivar and cropping system were only 3 and 2%, respectively. As Table 2 has already indicated, the percentage effects of the main significant factors for the other parameters were different ( Figure 3). In the case of β-D-glucan, TGW, and G12, the main dominant factor for their variability was cultivar (40.6%, 28.5%, and 6.0%).
The year as the most important factor was found in the CP (54.1%) and ASH (59.8%) parameters. The cropping system was then most significantly reflected in the variability of yield (18.1%), protein content (8.9%), and starch (12.8%). The low percentage effect of all major significant factors (9.1%) and a higher proportion of the sum of significant interactions (57.7%) are evident in the content of immunoreactive avenins (G12). Nevertheless, the effect of the sum of all factors (including their interactions) on the variability of this parameter was significantly lower compared to the others ( Figure 3).

Mutual correlations between AVNs and other grain parameters of oat
The contents of most AVNs were closely correlated with each other (r = 0.68-0.99; Figure 4). Only AVN-2fd showed a specific low insignificant correlation (r ≤ 0.25). Pearson's correlations of most AVNs to other monitored grain components were low to moderately strong (0.20 ≤ |r| ≤ 0.48). Significant negative correlations of most AVNs (excluding AVN-2fd) were recorded only for yield, protein content, and protein fractions (−0.48 ≤ |r| ≤ −0.26). Furthermore, the group of four AVNs (AVN 2f, AVN 5f, AVN 3p, and AVN 3f) significantly (positively) correlated with β-D-glucan (0.27 ≤ |r| ≤ 0.43) ( Figure 5). Only the mentioned AVN-2fd showed a very specific (negative) correlation to β-D-glucan (r = −0.43). Cluster analysis then graphically classified the tested parameters into three more separate clusters C1-C3 according to mutual correlation coefficients ( Figure 5). The cluster C1 grouped Σ AVN with a group of seven individual types of AVNs (except AVN-2fd). The cluster C2 was composed of a set of five following parameters (β-GLU, AVN-2fd, ST, TGW, and YLD). The more separate cluster C3 contained highly correlated parameters of As Table 2 has already indicated, the percentage effects of the main significant factors for the other parameters were different ( Figure 3). In the case of β-D-glucan, TGW, and G12, the main dominant factor for their variability was cultivar (40.6%, 28.5%, and 6.0%). The year as the most important factor was found in the CP (54.1%) and ASH (59.8%) parameters. The cropping system was then most significantly reflected in the variability of yield (18.1%), protein content (8.9%), and starch (12.8%). The low percentage effect of all major significant factors (9.1%) and a higher proportion of the sum of significant interactions (57.7%) are evident in the content of immunoreactive avenins (G12). Nevertheless, the effect of the sum of all factors (including their interactions) on the variability of this parameter was significantly lower compared to the others ( Figure 3).

Mutual Correlations between AVNs and Other Grain Parameters of Oat
The contents of most AVNs were closely correlated with each other (r = 0.68-0.99; Figure 4). Only AVN-2fd showed a specific low insignificant correlation (r ≤ 0.25). Pearson's correlations of most AVNs to other monitored grain components were low to moderately strong (0.20 ≤ |r| ≤ 0.48). Significant negative correlations of most AVNs (excluding AVN-2fd) were recorded only for yield, protein content, and protein fractions (−0.48 ≤ |r| ≤ −0.26). Furthermore, the group of four AVNs (AVN 2f, AVN 5f, AVN 3p, and AVN 3f) significantly (positively) correlated with β-D-glucan (0.27 ≤ |r| ≤ 0.43) ( Figure 5). Only the mentioned AVN-2fd showed a very specific (negative) correlation to β-D-glucan (r = −0.43). Cluster analysis then graphically classified the tested parameters into three more separate clusters C1-C3 according to mutual correlation coefficients ( Figure 5). The cluster C1 grouped Σ AVN with a group of seven individual types of AVNs (except AVN-2fd). The cluster C2 was composed of a set of five following parameters (β-GLU, AVN-2fd, ST, TGW, and YLD). The more separate cluster C3 contained highly correlated parameters of CP and protein fractions avenins (AVE) and glutelins (GLU), to which the parameters ASH, G12, and fat content (FT) were assigned. Despite the high variability of AVNs due to external (non-genetic) factors, their negative correlations to the increasing content of CP and conversely a positive correlation to β-D-glucan (except minor AVN-2fd) is evident.
Plants 2021, 10, x FOR PEER REVIEW 9 of 19 ASH, G12, and fat content (FT) were assigned. Despite the high variability of AVNs due to external (non-genetic) factors, their negative correlations to the increasing content of CP and conversely a positive correlation to β-D-glucan (except minor AVN-2fd) is evident.  It is also interesting to compare the low but significant correlations of both storage protein fractions (AVE and GLU) to the content of immunoreactive avenins (G12). Although the reactive (homologous) components of gluten are detected in the alcohol-soluble avenin fraction using the test kit AgraQuant Gluten G12 (see the manufacturer's protocol), a higher significant correlation to G12 content was confirmed with the content of soluble glutelins (rGLU = 0.37 vs. rAVE = 0.28) (Figure 4).

Effect of Weather Conditions on the Variability of AVNs
Principal component analysis (PCA) and Spearman's correlations ( Figure 5) were used to estimate and illustrate the relationships between Σ AVNs and selected weather parameters on the background of all tested oat cultivars, both cultivation systems, different localities, and three years. Both principal components explained together 68.8% of the total variability (the first: 41.78%, the second: 27.02%).
Principal component analysis (PCA) and Spearman's correlations ( Figure 5A,B) were used to estimate and illustrate the relationships between Σ AVNs and selected weather parameters. In the case of PCA analysis (5A), the mutual relationships are summarized on the background of all tested oat cultivars, both culture systems, different localities, and three years of evaluation. Spearman's correlation further describes these relations on the background of five individual oat cultivars ( Figure 5B).
Both principal components of PCA explained together 73.18% of the total variability (the first: 46.33%, the second: 26.85%). Closer positive relations to the variable Σ AVNs were mainly confirmed by the sum of precipitation in May (V_P) and June (VI_P). In contrast, the average July temperatures (VII_T) showed an antagonistic relationship to the Σ AVN contents.
Subsequent calculations of Spearman's correlation coefficients (rs) between Σ AVNs and weather parameters performed for individual cultivars ( Figure 5B) confirmed positive, strong, and statistically significant correlations between the sum of precipitation in May (V_P) and the growth of Σ AVNs (0.61 ≤ |r s | ≤ 0.83). Positive medium to strong correlations, which were even statistically significant in the case of Seldon, Kertag, and Korok cultivars, were also confirmed by the relationships between the sum of precipitation in June (VI_P) and Σ AVNs (0.47 ≤ |r s | ≤ 0.81). It is also possible to mention the trend of antagonistic relations between the average temperatures in June and July (VI_T and VII_T) and Σ AVNs. In the case of the Seldon cultivar, these correlations were even statistically significant −0.65-≤ |r s | ≤ −0.59).

Discussion
More than 40 distinct AVNs in oat seeds have been published so far [8,19]. However, the levels of these various AVNs in the samples are usually very different [8]. Their levels can be in intervals of more than three orders of magnitude. Therefore, we focused only to assess the eight most commonly abundant AVNs in our samples, whose levels were reported to be usually higher than 1 mg/kg dry weight (dw).
The three main detected AVN concentration ranges (AVN-2P-AVN AVN-2f and AVN-2c) were found in oats [4]. Our previous study [8] revealed in 10 selected varieties of oats both a high proportion of three main AVNs (65-70%) and a range of another 33 AVNs, including the structures of 10 novel AVNs. That study also identified nine additional quantitatively significant AVNs (≥1 mg/kg of dw), which represented approximately 25-30% of Σ AVNs [8]. In the context of these results, it can be estimated that the detected number of eight of the most quantitatively significant AVNs in this study represented the majority of all free AVNs of the tested grain. Chromatographic separation of these eight AVNs is shown in Figure S2.
Immunoreactive avenins (G12), which as a single parameter showed a similar level of variability as AVNs, indicated in some cases a slight exceedance of the concentration limit of gluten (20 mg/kg of dw) safe for celiac patients [5]. In these cases, the possible individual contamination of samples with residual gluten from other cereal species during harvesting or post-harvest operations of the samples cannot be completely ruled out. However, some inaccuracies may be related to the unclear suitability of the diagnostic immunochemical kit for monitoring gluten epitopes in oats. On the other hand, the Spanish authors have successfully used this system to identify oat genotypes with different contents of gluten epitopes. At the same time, the results of the genotypes correlated well with the clinical analysis of blood samples in patients with celiac disease [20]. Vice versa, Gilissen et al. [5] disagreed with the interpretation of the immunogenic responses of the G12 monoclonal antibodies (Anbs) in the case of oat. According to these authors, it could be caused by crossreactivity with some homological sequences in avenins, which do not have to correspond with clinically verified toxic sequences. Currently, only two adenine sequences, which are perhaps recognized by G12 Anbs, are mentioned as resistant to trypsin and chymotrypsin digestion [21], but these most likely occur in every oat cultivar [5].
The total average values of other monitored oat parameters (Table 1), including grain yield and its qualitative parameters (TGW, content of crude proteins, fat, starch, β-Dglucan, and ash) corresponded to published values [22,23]. A 2-3% higher content was on average detected in the avenin fraction-AVE (about 15% in crude protein) compared to the results of Van den Broeck et al. [24]. Concentrations of extractable glutelins (GLU) were similar to those of AVE with a higher RSD value (20.9%). Additionally, in this case, these values were higher compared to the review publication [25] mentioning the proportion of extractable GLUs below 10% in total protein. The main cause of higher contents of the mentioned protein fractions as well as crude protein can probably be found in connection with significantly higher vegetation temperatures in all monitored years in comparison with long-term averages in both localities ( Figure 6). As previously mentioned by Capouchová et al. [26], the effect of higher temperatures is manifested by increased plant respiration, which reduces the amount of carbohydrate assimilated and thus increases the percentage of proteins.
individual contamination of samples with residual gluten from other cereal species ing harvesting or post-harvest operations of the samples cannot be completely ruled However, some inaccuracies may be related to the unclear suitability of the diagn immunochemical kit for monitoring gluten epitopes in oats. On the other hand, the S ish authors have successfully used this system to identify oat genotypes with diffe contents of gluten epitopes. At the same time, the results of the genotypes correlated with the clinical analysis of blood samples in patients with celiac disease [20]. Vice v Gilissen et al. [5] disagreed with the interpretation of the immunogenic responses o G12 monoclonal antibodies (Anbs) in the case of oat. According to these authors, it c be caused by cross-reactivity with some homological sequences in avenins, which do have to correspond with clinically verified toxic sequences. Currently, only two ade sequences, which are perhaps recognized by G12 Anbs, are mentioned as resistant to t sin and chymotrypsin digestion [21], but these most likely occur in every oat cultivar The total average values of other monitored oat parameters (Table 1), including g yield and its qualitative parameters (TGW, content of crude proteins, fat, starch, β-Dcan, and ash) corresponded to published values [22,23]. A 2-3% higher content wa average detected in the avenin fraction-AVE (about 15% in crude protein) compare the results of Van den Broeck et al. [24]. Concentrations of extractable glutelins (G were similar to those of AVE with a higher RSD value (20.9%). Additionally, in this these values were higher compared to the review publication [25] mentioning the pro tion of extractable GLUs below 10% in total protein. The main cause of higher conten the mentioned protein fractions as well as crude protein can probably be found in con tion with significantly higher vegetation temperatures in all monitored years in comp son with long-term averages in both localities ( Figure 6). As previously mentioned by pouchová et al. [26], the effect of higher temperatures is manifested by increased p respiration, which reduces the amount of carbohydrate assimilated and thus increase percentage of proteins.  [4,15,16]. In particular, our results corresponded well with a study by Oraby and Ahmed [27] which demonstrated that harvesting years and planting locations affected by higher levels of biotic and/or abiotic stresses induce higher levels of AVNs, nevertheless in a genotype-dependent manner.
One aspect that is very interesting is the mentioned significant increase in the content of AVNs only at the conventional locality in CB in 2019 and 2020. It can be assumed that generally higher plant density in a conventional cropping system, combined with a higher sum of precipitation in July and August 2019 and 2020 (Figure 6), could generate more favorable microclimatic conditions for potential stressors, causing an increase in AVNs. Tybursky et al. [28] also mentioned in this context that the strict adherence to good agricultural practices in organic farming (such as good sowing practices, organic fertilization increasing the diversity of soil micro-organisms) should significantly protect crops from pathogen attacks. Since no other pesticide treatment other than herbicides has been applied in the conventional cropping system, a higher pressure of pathogens in the conventional cropping system is associated with an increased response of AVNs as protective phytoalexins can generally be expected. On the other hand, there was no visible difference in the health status of the plants kept under conventional or organic conditions. At the same time, the trend of AVN content in oat cultivars in conventional and organic cropping systems at the PR locality was rather the opposite ( Figure 1).
The above-detected effects of the four main factors on other parameters (Table 3) corresponded to a number of scientific studies. For example, in accordance with other scientific studies [29,30], the contribution of mineral nitrogen in the conventional cropping system had a positive effect on the growth of CP, storage protein fractions (AVE and GLU), TGW, and yield of grain (YLD). On the contrary, the starch content decreased. The significant effect of the cultivar as well as significant annual differences in β-D-glucan content were also confirmed [31]. These authors also mention the significant effect of mineral nitrogen fertilization on the increase of β-D-glucan. However, this was not confirmed by our results.
The evaluation of the percentage influence of significant factors (Figure 3) confirmed the high and-at the same time-specific level of the influence of individual main factors (including their interactions) on the variability of the monitored parameters. The mentioned exception was only the content of immunoreactive avenins (G12), indicating a certain level of undefined variability. This fact can be caused by a lower degree of epitope specificity of the oat avenin kit [5]. Another source of uncertainty in the G12 content is the lower but significant positive correlation of the glutelin fraction with the G12 content ( Figure 4). The presence of a variable proportion of the glutelin fraction with reactive groups of peptides in the isolation of avenines using the commercial ELISA test kit AgraQuant cannot be completely ruled out so far. Despite the very low celiac reactivity of oat protein, which is generally safe for celiacs, more attention should be paid to this glutelin fraction.
According to our findings, studies correlating a wider range of AVNs to selected nutritional and technological parameters of oat grain are not yet freely available. Studies explaining the causes of significant variability in the content of AVNs in multi-year field experiments are also very limited. Due to the large influence of non-genetic factors (especially weather conditions), the lower mutual correlations between AVN and other grain parameters (Figure 4) are not surprising. It is also evident that the analysis of a larger number of oat cultivars (genotypes), in combination with greater environmental variability, will be necessary to further refine these correlations. Despite these facts, the significant effect of cultivars on the content of AVNs and the detected positive significant correlation of some AVNs to genotype-highly dependent β-D-glucan [23] could be one of the breeding strategies for obtaining oats with improved nutritional quality.
PCA analysis and Spearman correlations based on individual cultivars confirmed that higher sums of precipitation, especially in May and June, have a significant effect on the increase of Σ AVN in hulled and naked oats ( Figure 5A,B). Due to the highly above-average temperatures in all three monitored years, higher precipitation amounts could have an essential effect on plant development and the beginning of grain formation, including the synthesis of AVNs. Some major differences in the content of Σ AVN between the organic and conventional cropping system at the identical locality and year (e.g., CB in years 2019 and 2020, see Figure 1) also indicated that precipitation probably evoked other effects (stressors), which increase Σ AVNs further accelerated. The scientific literature mentions two main factors is a possible high increase in the content of Σ AVNs. The first factor (stressor), to which more humid conditions suit and which has been shown to increase contents of AVNs, are fungal diseases, e.g., oat rust [32]. However, a higher incidence of fungal diseases in the CB locality in the conventional variant compared to the organic system was not recorded.
The high sum of precipitation in the final stage of grain ripening, often in combination with lodging, can cause the pre-harvest sprouting of grains. Thus, it can be another factor that could significantly increase the Σ AVN content. Significant increases of Σ AVNs in germinated grain have been confirmed from a number of studies [33,34]. Although the individual Spearman correlation coefficients between Σ AVNs and the sum of precipitation in July and August (VII_P and VIII_ P) showed trends of positive moderately strong correlations ( Figure 5B), in our experiments, the visual control of grains did not reveal potential latent pre-harvest sprouting. On the other hand, the above-mentioned higher density (interconnection) of plants in the conventional crop system could, due to aboveaverage precipitation, generate a wetter microclimate, evoking latent preharvest sprouting more easily. Preharvest sprouting of grain, which significantly damages starch granules and other important technological parameters of cereals, is not commonly monitored in oats in comparison with wheat and rye [35]. The quantification of damaged starch using the Megazyme starch damage kit was evaluated as part of a nutritional analysis of oats [24]. However, the detected ranges of 1.8-4.0% were not related to the definition of any technological limit or possible changes of bioactive compounds. Despite as yet probably unpublished relationships between the level of damaged starch and the AVN content, a close positive correlation can be expected.
Based on these results, it can be concluded that a basic accelerator of AVN growth in oat grain will be above-average rainfall during grain ripening. Nevertheless, precipitation will probably only generate the basic conditions for the development of other factors (abiotic and biotic stressors), which further trigger the direct synthesis of AVN in oat grain.

Field Experiments
The exact field experiments with 5 selected oat varieties were carried out in the Czech Republic in two different localities in two cropping systems during three vegetation periods in 2018-2020. We obtained 60 combinations of samples (cultivar-5, locality-2, cropping system-2, and year-3) which were measured in 4 replicates. In total, 240 oat samples were evaluated in this study.
The After harvest, the yield of grain reported at 14% moisture and TGW (Thousand Grain Weight) were determined. Then, the oat grains of the hulled varieties were dehulled using a laboratory dehuller Codema LH 5095 (Maple Grove, MN, USA) and prepared for future analyses.

Weather Conditions
Weather vegetation conditions (April-August) at both localities in individual years (2018-2020) are shown in Figure 6. The course of mean monthly temperatures in both localities was similar and all 3 years (2018-2020) in the Czech Republic can be described as very warm compared to the long-term temperature average. However, significant differences between the two localities can be found in the monthly sum of precipitation for vegetation, which was significantly higher at the locality in CB. In particular, the sums of precipitation in May exceeded the sum of precipitation at the PR locality in all years, as well as the long-term average of the locality in the CB. The sum of temperatures during the vegetation period showed that the locality in the PR was warmer by an average of 1 • C only in 2018. In the following years, the sum of vegetation temperatures was comparable between the two localities. However, the achieved ranges of total precipitations 176 mm (2018), 251 mm (2019), and 249 mm (2020) for PR locality lagged significantly behind the CB locality with total precipitations 302 mm (2018), 320 mm (2019), and 488 mm (2020).

Basic Quality Parameters of Oat Grain
Analyses included quantification of following basic quality parameters of the oat grain: content of crude protein (CP) according to Kjeldahl method [36] on Kjeltec KT 200, FOSS, Sweden, and the protein content was calculated by multiplying the nitrogen content by 6.25, content of starch (ST) according to Ewers' polarimetric method [37] on Polamat A, Carl Zeis Jena, Germany, content of fat (FT) according to Randall method [38] on Randall Hot Extraction apparatus E6, Behr Labor Technik GmbH, Düsseldorf, Germany, and ash content (ASH) after burning at 750 • C for 4 h [39]. The dry matter content of seed samples (5 g) was further dried in an electric hot-air drier at 105 • C for 4 h, according to the standard method [40].

Analyses of Specific Nutritional Parameters of Oat Grain
Analyses of 8 selected most abundant avenathramides defined according to Dimberg's nomenclature (AVN 2c, AVN 2f, AVN 2p, AVN 2fd, AVN 2pd, AVN 3p, AVN 3f, and AVN 5f) were carried out according to Jágr et al. [8]. Briefly, oat flour was extracted with 80% MeOH solvent, then it was filtered and analyzed using UHPLC system (Dionex, UltiMate 3000 UHPLC system, Dionex Softron GmbH, Germany) coupled with QExactive Orbitrap mass spectrometer (Thermo Fisher Scientific, Rockwood, TN, USA). Chromatographic separation of AVNs was realized using gradient elution with 0.1% formic acid in water as a solvent A and methanol with 0.1% formic acid as a solvent B. LC gradient: 0 min: 75% of solvent A + 25% of solvent B, then linear gradient in 11.0 min: 25% A + 75% B. Column was flushed with 100% of B in 12.0 min and then column was equilibrated back with 75% A and 25% B in 15.0 min. The column was maintained at 40 • C at a flow rate of 0.35 mL/min. MS detection was performed using electrospray ionization in positive mode. Identification of AVN compounds in samples was based on their retention times and on mass spectral data obtained by UHPLC-MS/MS, which were compared with those described in our previous study [8]. The amounts of AVNs in the samples were calculated using calibration curves generated by the analysis of three external standards (AVN 2c, AVN 2f, and AVN 2p) purchased from Sigma-Aldrich (St. Louis, MO, USA). These standards were also used for the semi-quantitative analysis of all the other AVNs (AVN 2fd, AVN 2pd, AVN 3f, AVN 3p, and AVN 5f), for which no commercial authentic standards were available.
Immunoreactive gluten (avenin) peptides (G12) were assessed using ELISA test kit AgraQuant Gluten G12 (Romer labs Diagnostics) according to the manufacturer's protocol. All ELISA measurements were in a separate room to avoid gluten contaminations. The absorbance was recorded at 450 nm using Sunrise (Tecan) microplate reader [42].
The separation of both protein fractions was determined by high-performance liquid chromatography (RP-HPLC) [44]. Protein extracts (10 µL) were injected into RP-HPLC using a Waters 2965 apparatus with UV detector and 300 SB-C8 Zorbax Poroshell ™ column (75 × 2.1 mm, 5 µm particles) linked to a Zorba 300SB-C8 cartridge guard column (Rockland Technologies, Inc., Newport, DE, USA). Elution of avenins and glutelins was monitored at the wavelength 210 nm. Wheat gliadin (Sigma Aldrich) was used as an external standard for final quantification.

Statistical Methods
All chemical and nutritional parameters were assessed in 4 repetitions. Basic statistics included the determination of the mean, standard deviation (SD), standard error (SE), minimum/maximum, and relative standard deviation (RSD) in all 19 tested parameters.
The next statistical methods further included 4-way analyses of variance (ANOVA) followed by Tukey's Honest Significance Test (Tukey HSD test) to determine significant differences between tested parameter averages. The percentage effect of significant factors was calculated based on the ANOVA as a percentage of the sum of squares (% SS) of the significant factor to the total sum of squares [46].
Pearson's correlation matrix of individual parameters with mutual hierarchical clustering (standardize Euclidean distance and UPGMA method as a clustering technique) were also used.
The applied calculations of Principal Component Analysis (PCA) [47] and Pearson's correlations were used to illustrate the relationship between AVNs as selected parameters of weather conditions. All above mentioned statistical methods were computed with the Statistica 7.1CZ software StatSoft, Inc. (2005) (StatSoft, Tulsa, USA).

Conclusions
The achieved three-year results confirmed the quantitative differences of a wider spectrum of three main and five minor AVNs between both oat cultivars and individual types of AVNs. Individual AVNs confirmed a very high variability (RDS = 72.7-113.5%), which was dominantly influenced by the locality and year factor. The effect of the cultivar and cropping system was also statistically significant, but with a significantly lower share in the total variability of AVN content. With the exception of AVN-2fd, the other AVN types confirmed high mutual correlations 0.7 ≤ |r| ≤ 0.9. AVNs showed only lower to medium correlations to the other 10 grain parameters. Despite this fact, detected significant correlations (0.27-0.46) between four AVNs (AVN 2f, AVN 5f, AVN 3p, and AVN 3f) and strongly genetically fixed β-D-glucan contents could be effectively used in breeding programs for the synergetic increase of both parameters. PCA analysis and Spearman correlations based on individual cultivars confirmed significant positive relationships between the sum of precipitation in June and July and the increase of Σ AVNs. Precipitation during this period at otherwise very above-average temperatures could favorably affect grain formation as well as AVN synthesis. However, the results also indicate that higher precipitation can also generate favorable conditions for related factors, such as preharvest sprouting, evoking a direct increase in the AVN content of the oat grain.