Adaptation and Agronomic Performance of Domesticated Moroccan Oat ( Avena magna ssp. domestica ) Lines under Subsistence Farming Conditions at Multiple Locations in Morocco

: Common hexaploid oat ( Avena sativa L.) is an important global cereal crop. A Moroccan tetraploid sister species, A. magna Murphy et Terrel, was exclusively a wild species until recently. The goal of domestication was to exploit its superior groat-protein content and climatic tolerances. We set up replicated trials of 41 domesticated A. magna lines on eight Moroccan farms during the 2017–2018 and 2018–2019 growing seasons. Twenty traits were measured and analyses of variance detected signiﬁcant differences among lines. The highest grain yield was at Berrechid in 2017–2018 (63.56 q/ha), with an average annual yield across sites of 43.50 q/ha, the site factor explaining 82% and the genotype-environment interaction explaining 15% of the variability. In the second year, El Kebab recorded the highest yield at 20.03 q/ha over the annual average of 14.78 q/ha. In this second year, the site factor was highly signiﬁcant, explaining 42.25% of the variation, with the genotype-environment interaction explaining 26.61% of the variability. An additional main effect and multiplicative interaction analysis of the eight two-year trials identiﬁed several accessions with good yield stability. Twelve lines exhibited a ASVs ≤ 1.50, with ﬁve accessions (A34, A40, A23, A05, A04) exceeding the overall average yield of 29.53 and A34 having the greatest mean grain yield and stability. The versatility and stability of A. magna can provide a sustainable protein source and an economic resource for farmers seeking products that are resilient to climatic instability. A. magna as a climate change-resistant crop for human and animal nutrition in vulnerable areas of the subtropical developing world and especially in well-drained soils.


Introduction
The oat genus Avena L. (x = 7) includes the seventh most important cereal worldwide: common or white oat (A. sativa L. 2n = 6x = 42; AACCDD subgenomes). Other domesticated taxa include red oat A. byzantina C. Koch (AACCDD), a fodder oat grown in mild winterproduction conditions; Ethiopian endemic oat A. abyssinica Hochst (AABB); the diploid lopsided or sand oat complex (A. strigosa Schreb. A s A s ); and the recently domesticated A. magna Murphy et Terrel ssp. domestica Ladizinsky [1]. The oat genus' center of origin is the Maghreb and southern Iberia, though the most likely tetraploid progenitor of the hexaploid species-A. insularis Ladizinsky (CCDD)-is currently found only in Tunisia and Sicily in the central Mediterranean [2,3]. The hexaploid cultivated forms appear to have been domesticated from weedy A. sterilis L. on at least two occasions in the ancient Near East [4,5], with a secondary center of diversification of the hulless or naked oat in Northwest China and Mongolia [6].
Avena magna (Moroccan oat) was recognized by Ladizinsky [1] as having approximately double the seed protein content of common oat. Consequently, Ladizinsky embarked upon a project to transfer some components of the domestication syndrome-reduced awn length, glabrous lemma and palea, non-shattering seed, and white hull color-from hexaploid common oat to tetraploid Moroccan oat through sexual hybridization using two cycles of backcrossing. One of his stable backcross lines displaying the domestication syndrome was Ba13-13. Oliver et al. [7] produced a linkage map for A. magna from a Ba13-13 × wild #169 recombinant inbred line (RIL) population containing 1013 molecular, three morphological, and one cytological marker arranged in 14 linkage groups. They also found that the domestication syndrome components displayed by Ba13-13 were associated with inheritance of genomic material from hexaploid oat and were distinguished from wild A. magna and similar-phenotype RILs by lack of a cytological knob on chromosome 5C. The tight linkage in the coupling phase with genes controlling shattering (Ba) and awn formation (Awn) at the terminus of linkage group 13 (chromosome 5C), with the third major gene for lemma pubescence (Lp) near the terminus of linkage group 11, suggested that these traits could be easy to retain in further crosses to develop a broad, diverse gene pool for cultivated Moroccan oat breeding. Unfortunately, while the domestication syndrome transfer from hexaploid oat into A. magna represented an improvement over wild A. magna, introduction of other traits important for broad cultivation, such as awn deletion, complete resistance to shattering, semi-dwarfism, lodging resistance, erect growth habit, resistance to seed dormancy, phototropism (reduced time to flower and maturity), reduced groat length, and ease of dehulling proved difficult while retaining high seed protein content [8,9].
Years of further work, beginning with a wild-phenotype F 3 plant from the Ba13-13 × #169 population as described by Jackson [9], resulted in a series of genetically diverse A. magna lines having improved domestication traits without relying on the hexaploid 5C chromosomal material while retaining a protein content exceeding 25%. Moreover, this new domesticated species should be well adapted to changing climatic conditions in its native Morocco, where models forecasting 2-3 • C temperature increases and 10-20% precipitation decreases by 2050 have prompted policy changes away from optimizing and towards stabilizing crop production [10]. Additionally, the increased macro-and micrnutrient density of Moroccan oat [9] could help address stunting and other malnutrition problems that persist in the country, especially in rural areas [11]. The present study was designed to assess genetic diversity, yield stability, and overall adaptive behavior of 41 A. magna domesticated lines through eight agro-morphological characterization trials at seven geographically and climatically diverse sites in Morocco spanning the 2017-2018 and 2018-2019 growing seasons.

Plant Materials
Avena magna ssp. domestica lines in this study were developed from a single Ba13-13 × wild #169 recombinant inbred line (RIL) [7] expressing the wild-type growth habit and resistance to field races of crown rust (Puccinia coronata) in Baton Rouge, LA [9]. Seed from this line was subjected to six cycles of selfing and selection for improvements in agronomic and seed domestication, resulting in six foundational lines displaying various combinations of traits for domestication and cultivation including awn reduction, shattering resistance, semi-dwarfism, lodging resistance, erect growth habit, resistance to seed dormancy, phototropism, reduced groat length, and ease of dehulling [9]. A set of 114 RILs, including the 41 experimental lines in this study (Table 1), were produced by intercrossing the foundational lines as instructed by a virtual pedigree produced using a genotype/phenotype model with the JMP Genomics software package (SAS Institute Cary, Cary, NC, USA). The 41 lines were compared to the internal control 'Avery' (A40), a commercial A. magna ssp. domestica variety produced by General Mills in 2013 [9] and released commercially in Morocco late 2019.

Experimental Sites
Experimental site locations are mapped in Figure 1, shown photographically in Figure 2, and described in Table 2. Experimental trials of 2017-18 were conducted at four locations: a commercial potato production farm 14 km southeast of Berrechid (central Atlantic coastal plain); a modern production farm 15 km south of Meknes; a subsistence farm in Lahri (low-altitude, Khenifra region); and a subsistence farm 4 km southwest of El Kebab (high-altitude, Khenifra region). Trials were sown at Berrechid on 6 December 2017; at Lahri and El Kebab on 8 December; and at Meknes on 27 December. Experimental trials of 2018-19 were conducted in the fall-winter-spring growing season on an organic production farm 16 km west of Tiflet (northern Atlantic coastal plain); at the Royal Agricultural Domain Ain Hamra Farm in Seba Ayoun 20 km east of Meknes (north-central interior lowlands); at the El Kebab farm; and at a cooperative subsistence farm 18 km east of Youssoufia in Bouchane (semiarid Phosphate Plateau) with sowing dates of November 27, 28, 29, and 30, respectively. Locations were purposely selected to test the adaptational range of A. magna as a grain crop and, consequently, included a range of technological sophistication from commercial to subsistence; climate zones ranging from warm to cool and dry to humid; and traditional cultural contexts including both Amazigh (Lahri, El Kebab) and Arab (i.e., Bouchane).
consequently, included a range of technological sophistication from commercial to subsistence; climate zones ranging from warm to cool and dry to humid; and traditional cultural contexts including both Amazigh (Lahri, El Kebab) and Arab (i.e., Bouchane).   consequently, included a range of technological sophistication from commercial to subsistence; climate zones ranging from warm to cool and dry to humid; and traditional cultural contexts including both Amazigh (Lahri, El Kebab) and Arab (i.e., Bouchane).   Soil and climatic data at each of the seven experimental locations are provided in Table 2. The Bouchane and Tiflet experiments were placed on sites that had unworked fallow fields during the previous year. Tillage was performed using two moldboard crosspassages for the Berrechid, Bouchane and Meknes sites and using conventional tillage at the Tiflet and El Kebab sites. The sowing dose was determined for each line as a function of the germination rate of the seed. Weed control was done by hand, which took place over two periods: at the tillering and flowering stages. The number of weedings was based on the degree of weed infestation. For example, the Tiflet site was repeatedly weeded in accordance with the organic management of the farm. However, the Bouchane site experienced minimal weed infestation due to the exceptionally dry conditions. At El Kebab, nitrogen fertilization was carried out at the early tillering stage in the form of ammonium nitrate. Weather data from nearby stations and presented in Table 2 reveal that December was very dry in both years, but especially following planting in 2018-2019, and that sufficient rainfall to support abundant crop growth did not occur until very late in the growing season in March-April. Additionally, it should be noted that average temperatures were low throughout the winter of 2017-2018 in the Middle Atlas region encompassing the El Kebab site.

Experimental Design
The experimental set-up at the 2017-2018 locations was in randomized complete blocks (RCB) with two repetitions. Planting was done in paired three-meter rows, spaced 20 cm apart, with a seeding rate of 5 g per paired lines. The spacing between lines was 60 cm. Each site was characterized by peculiarities or variations in the planting plan in accordance with terrain constraints. The experiments in 2018-2019 were set up in randomized complete blocks (RCB) with 1.83 m rows in three blocks at Meknes and Tiflet; 1 m rows with three blocks at El Kebab; and 2 m rows with three blocks at Bouchane.

Measurements and Observations
Agronomic and morpho-physiological characteristics were measured to estimate yield potential through grain and straw productivity; yield stability across environments; tolerance to the three main oat diseases (crown rust (CR), barley yellow dwarf virus (BYDV), powdery mildew); and to determine the most optimally and stably productive lines in terms of yield. Measurements were taken at three different growth stages in order to carry out the disease and agro-morphological measurements.
At the flowering stage, rust resistance was assessed using a scale based on the percentage of leaf cover by the pathogen (severity of the disease). Tolerance to both powdery Agronomy 2021, 11, 1037 6 of 23 mildew and BYDV were assessed based on "presence" or "absence", though powdery mildew was only scored in 2017-2018. The evaluation of the three diseases was made for each line at the four sites in 2017-2018, though powdery mildew presence was erratic and was therefore omitted from the final analyses.
Observations were made on a sample of four plants at the seedling stage and on three samples from each plot at harvest. In each plot, we counted the stand and measured the agro-morphological parameters (stem height, length of the roots and number of tillers). For the individual plants, we separated the different organs (seeds, stems, roots) apart to determine their dry weight. For the rest of the plots, we also separated the seeds from the remaining biomass. All of the harvested tissue was then weighed after oven-drying at 70 • C for 48 h for stems, roots, and biomass and at 35 • C for 48 h for seeds hulled after threshing. From these measurements the harvest index was then calculated using Equation (1), after estimating the grain yield and straw yield. where: • HI: harvest index • GYH: grain yield (q/ha) • DYH: dry matter yield (q/ha).
For the estimation of grain and straw yields, we extrapolated the value of the average grain and straw yield of three random samples of 0.5 m row length. The thousand-seed weight was calculated by counting and weighing the threshed and oven-dried seeds. This allowed us to calculate the number of grains per individual plant using the following formula: where: • N grains : number of seeds per plant • P grains : grain weight per plant • TSW: thousand-seed weight At maturity, a sample of three individuals per line, block and site; and a sample from three 0.5 m row length per line, block and site constituted the material on which we carried out measurements and estimated the agro-morphological variables.

Statistical Analyses
The collected data were used to carry out statistical analyses. We assessed intra-locality variability by comparing the lines with each other at each agro-climatic site. On the one hand, the measured parameters were the subject of a descriptive analysis and an analysis of variance (ANOVA) with two sources of variation (lines and blocks). Before any analysis, the normality of the variables was tested for each indicator through a Kolmogorov-Smirnov (KS) test at a significance level of 5%. If the null hypothesis was rejected, a test to compare the Student-Newman-Keuls means (SNK) was used to distinguish the different homogeneous groups. For experimental sites where the variables associated with productivity (grain yield, dry matter yield, harvest index) presented significant differences due to the effect of the genotype, we carried out a principal component analysis (PCA).
Lastly, interlocality and growing season analyses for data were performed using an ANOVA with three factors (lines, environments and blocks), a principal component analysis (PCA), and an ascending hierarchical classification (AHC). Then, the evaluation of the genotype × environment (G×E) interaction on the basis of all the parameters evaluated was carried out using an additive main effects and multiplicative interaction (AMMI) analysis. The AMMI is a hybrid procedure that provides visual inspection and interpretation of the components of the yield G×E interaction [12]. The AMMI model separates the main additive effects from genotypes and environments using ANOVA, and then analyzes the effect of interaction using the multiplicative model provided by the PCA. This analysis mainly provides information on the productivity and stability of the genotypes. Consequently, the AMMI stability value (ASV) was used to quantify genotype yield stability as a function of the first two axes of the main components of the interaction, and was calculated according to the formula developed by Purchase [13][14][15]: where:

Assessments of Disease Tolerances
The observation of diseases was carried out in the middle of the cycle for all crops, specifically on the 92nd, 111th, 118th, and 156th days after sowing at the Meknes, Berrechid, Lahri, and El Kebab sites in 2017-18, respectively. The three most important diseases attacking oats in Morocco are crown rust (Puccinia coronata f. sp. avenae), BYDV, and powdery mildew (Blumeria graminis f. sp. avenae). Crown rust developed by the middle of the cycle at the Berrechid site but was not observed at the other locations until very late in the growth cycle. We therefore assessed crown rust severity only at Berrechid via visual scoring using the modified Cobb scale [17]. Following this we performed a two-factor ANOVA, with the summary table showing only significant differences presented in Table 3. Powdery mildew was only observed at Berrechid and was therefore not scored. As for BYDV, it was observed at all sites, though infestations within a given site did not appear to be uniform.
In the 2018-2019 experiments, observations of diseases were made on the 152nd day after sowing at the Meknes and El Kebab sites and on the 154th and 161st days for the Bouchane and Tiflet sites, respectively. The assessment focused on two diseases: crown rust and BYDV.
For the assessment of crown rust tolerance, the lines and blocks terms showed no significant differences among the 41 lines (Table 3) in 2017-2018. On the other hand, the lines*blocks interaction showed a highly significant difference, indicating that there was likely a non-uniform infestation across the plots-something that should not be surprising given that oat is a rare crop in these growing areas of central Morocco. During a 209th day post-sowing visit at Lahri, controls T1 and T2, as well as lines A01, and A02 displayed crown rust symptoms. At Meknes and El Kebab, we noticed some symptoms of crown rust on several lines at harvest time. Table 3. Summary two-factor ANOVA table for significant disease traits. Crown rust (Puccinia coronata f. sp. avenae) infestation was sufficiently uniform to allow for assessment at Berrechid in 2017-2018 and at all sites in 2018-2019. NS = not significant; * significant at p = 0.05; ** significant at p = 0.01; *** significant at p = 0.001.  (Table 3), significant differences were found among lines, while the non-significant differences for crown rust infestation detected by the line*block ANOVA interaction indicated there was uniform infestation throughout the plots. At Tiflet, the infestation was assessed on only the first block; however, from the ANOVA it is clear that there were significant differences among the lines and there was likely an exogenous source of the inoculum. Because of its relative proximity to the Atlantic Ocean, the Tiflet site experiences chronic crown rust infestations [18]. The narrow gradient of the distribution can be explained by the dissemination of the pathogenic agent into the experimental plot, which was cleared of riparian vegetation and cultivated for the first time for this particular study.

Site and Season
Since humid conditions favor the development of the Puccinia coronata fungus, drought conditions during winter 2018-2019 at Meknes and El Kebab likely explained the absence of crown rust infestations at these study sites. This can also possibly be explained by local applications of fungicides on nearby cereal fields, especially at Meknes where the plots were located on a large commercial farm within the Domaine Royaume. This hypothesis is supported by the results from the prior year's study when little crown rust appeared, in spite of more abundant precipitation.
The means comparison for degree of BYDV infestation at El Kebab in 2018-2019 enabled us to distinguish three homogeneous groups, with lines A01 and A03 being the most susceptible (data not shown). The average infestation at El Kebab was the lowest in both growing seasons.
For BYDV, the disease was present at all sites, with by far the highest incidence being at Tiflet. This observation is not surprising, given the relative abundance and diversity of aphid vectors of this disease [18,19], combined with A. magna's well-known susceptibility to BYDV [1]. However, the lowest incidence of BYDV in both growing seasons was at El Kebab, suggesting that the relatively lower temperatures, higher elevation, and distance from oceanic humidity at this site were unfavorable for the aphid vectors of this disease. Future A. magna germplasm enhancement will have to put more effort into BYDV tolerance, but it will likely rely on variation induction because of the lack of resistance genes in the wild species [1].

Analyses of Agro-Morphological Traits
Significant agro-morphological effects at the eight environments are highlighted via descriptive statistics provided in Table 4 and two-way ANOVA Table 5. Variables studied included the measured parameters, the calculated parameters, the thousand-seed weight (TSW) and the harvest index (HI). Explanations: PH = plant height (cm); RL = root length (cm); NFT = fertile tiller number; RW = root weight (g); DYP = dry matter weight per plant (g); NGP = number of grain per plant; GWP = grain yield per plant (g); TSW = thousand-seed weight (g); GYH = grain yield per hectare (q); DYH = dry matter yield per hectare (q), HI = harvest index.   Descriptive statistics of the agro-morphological traits are presented in Table 4. Based on the internal analyses at each site, PH and TSW were the most stable variables as they present low coefficients of variation. The highest CVs were observed for DYP, NGP, and GYP. In the first growing season of 2017-2018, plant development was much greater than the second season. The PH values measured between 110.24 and 176.86 cm in the first year and 45.15 and 79.14 cm the second year, and the average grain yield between 36.79 and 14.78 q/ha in the two years, respectively.
In the 2018-2019 season, Bouchane had the shortest plant size (45.15 cm) in comparison to the other sites, and the average RL of the lines was more proportionate to the PL than at the other sites. This demonstrated the capacity of A. magna ssp. domestica roots to penetrate deeply into sandy-loamy soils under harsh growing conditions (Figure 2).
The accessions in Tiflet presented the lowest values for the agro-morphological parameters (Table 4). We suspect this was mainly due to the poorer growing conditions of the trial of the new managed loamy-sandy soil of an organic farm. The plants were also exposed to severe competition from weeds, pests, and diseases, which strongly impacted the growth and development of the lines (Figure 2).
The two-factor ANOVA revealed significant differences among accessions. At Berrechid, the genotype effect was highly significant for the variables PH and GYH and moderately significant for RW and HI ( Table 5). The line*block interaction was moderate for the PH parameter and minor for RW, NGP, GYP, and GYH.
The Bouchane site ANOVA analysis showed highly significant genotypic effects for six of ten variables (NFT, RW, GYP, GYH, DYH, and HI), a moderate effect of PH, and a minor effect of RL ( Table 5). The block factor had significant effects on all the variables except RL. The line*block interaction effect was highly significant for PH, moderate for RW and GYP, and significant at the 5% level for the rest of the parameters.
At the Lahri site, the ANOVA (Table 5) showed no significant effects linked to the line factor for only RW. The block factor showed highly to moderately significant effects for all the traits. The line*block interaction showed highly significant effects for all parameters and slightly significant (at the 5% level) for RW and DYH.
At El Kebab the effect of the lines factor was highly significant for all the parameters in the 2018-2019 growing season but not significant for DYP, NGP, and GYP for the 2017-2018 season ( Table 5). There was no significant block effect for PH, RL, NGP, and GYP for the first year trials. Finally, the line*block interaction was highly significant for PH in the first season and for NGP, DYH, and HI in the second season.
At Meknes, the 2018-2019 trial revealed significant differences among lines for GYP and GYH, and highly significant differences for the rest of the traits ( Table 5). The block and line*block interaction factor effects were highly significant for PH, DYP, NGP, GYP, GYH, DYH, and HI, but the block factor was not significant for the NFT and the line*block interaction for RL.
At the Tiflet site, the test detected highly significant differences among the lines for all the variables except NFT ( Table 5). The block factor had only a slightly significant effect on DYP, DYH, and HI and no effect on all other variables. The line*block interaction was highly significant for PH, RL, NFT, DYP, DYH, and HI but not for RW, NGP, GYP, and GYH.

Analyses of Productivity Traits
To assess the productivity of the A. magna oat lines, we analyzed the GYH, DYH, and HI of the eight experimental trials. Descriptive statistics for these productivity traits appear in Table 6 and the rankings of the 41 lines in Table 7. In comparing the two seasons, genotype variability as measured by intra-site CV's was much higher in the second year (2018-2019), 83.66% against 48.19% for GYH 70.21% against 40.47% for DYH, and 45.17% against 27.62% for HI (Table 6). For GYH, the greatest variability was found at Meknes (96.55%), which explains the high variability of HI in the same trial (0.66%). DYH was most variable at Bouchane (CV = 88.71%). Berrechid had the least variability for both productivity traits; GYH (40.69%), and DYH (33.98%). The highest HI variability (66.32%) was that of Meknes in the second trial season, 2018-2019.
In general, the productivity parameters for the 2017-2018 trials were higher than the 2018-2019 trials, with respective means of 43.48 and 14.68 q/ha for GYH and 161.92 and 39.85 q/ha for DYH ( Table 6). The HI values were lower in the first season with an average of 22.52% compared to the second season (28.76%). The overall mean grain yield was 29.08 q/ha, with the highest value at Berrechid (63.56 q/ha) followed by Meknes's first season (56.56 q/ha), and the lowest at Tifelt (5.10 q/ha). For DYH, the average overall yield was 100.89 q/ha with a maximum at Berrechid (255.08 q/ha) and a minimum at Tifelt (18.08 q/ha). The best HI was that of Bouchane (39.24%) and the lowest at Lahri (17.30%).
In 2017-2018, the Meknes site had the second-highest HI, with an average of 22.57%. The two lowest sites for HI in 2017-2018 were Berrechid and Lahri, 20.12 and 17.30%, respectively. In the second-year trial, Bouchane, followed by El Kebab, presented the highest HI (39.24 and 33.40%, respectively), compared to Meknes and Tifelt at 18.72 and 23.70%, respectively ( Table 6).
The two-way ANOVA detected highly significant differences for GYH among the 41 lines at all sites except Meknes in 2017-18 (Table 5). For DYH and HI, the differences were highly significant at six and eight sites, respectively, and were slightly significant (at the 5% level) at the Berrechid and Meknes first-year trials for DYH.
The overall ranking of lines for GYH, DYH, and HI are presented in Table 7. GYH ranged from a low of 16.96 q/ha for A22 to a high of 30.28 q/ha for A23; the A40 control line yield across the eight sites was 20.82 q/ha. DYH ranged from 64.44 to 102.83 q/ha with the respective accessions A22 and A02. For HI, A15 had the best value of 32.19% and A25 the lowest of 23.21%; the HI of the control, A40 was 24.36%.
In the first year of trials, A35, A27, and A14 ranked first in Berrechid for GYH, reaching yields of 108.73, 88.85, and 85.45 q/ha, respectively ( Table 8). The control (A40) yielded 73.18 q/ha while line A22 produced the lowest GYH (38.83 q/ha) at the site. Meknes the second-best site; A37 and A16 ranked first with GYH at 83.90 and 75.02 q/ha, respectively; line A40 produced 61.82 q/ha; and the last lines, A36 and A13, yielded 39.22 and 34.87 q/ha, respectively.  In the first-year trials at El Kebab, lines A04, A15, and A25 had the lowest GYH 15.09, 15.23, and 15.41 q/ha, respectively, the control yield was 20.08 q/ha, and A36 and A38 recorded the highest GYH, 44.80 and 38.34 q/ha, respectively. The most productive line in Lahri, A16, produced 48.57 q/ha and the lowest yielding accession, A42 produced only 9 q/ha.
In the second-year trial, the early inauspicious conditions for optimum plant growth (December-February) were suitable for genotype discrimination at Bouchane and Tiflet. At Bouchane, lines A20 and A33 were the highest yielding genotypes, 28.86 and 26.29 q/ha, respectively ( Table 8). The exceptionally dry conditions at the Bouchane site are best conveyed photographically; Figure 2 shows conditions on 3 May, approximately one month before harvest, where the A. magna plots provided seed when the wheat and barley crops elsewhere on the farm were a complete loss, being plowed under or opened up for livestock grazing. In terms of dry matter yield, line A02 was highest at 43.93 q/ha. At Tiflet, line A20 recorded the highest GYH, at 20.34 q/ha ( Table 8). The highest DYQ was in lines A33 and A41 with a biomass production of 32.93 and 31.14 q/ha, respectively. As with the agro-morphological traits, the lines at the Tiflet site were likely well below their productivity potential due to the poorly managed organic production conditions. At El Kebab, GYH and HI were higher compared to the other sites. Line A02 produced an average yield of 44.8 q/ha, followed by A41 with a yield of 38.34 q/ha (Table 8). Lines A43, A02, and A44 had the highest harvest indices at 40.52%, 39.23%, and 39.22%, respectively. The grain yield results at the El Kebab site were significantly higher in the second experimental campaign (2018-2019) compared to the first (2017-2018) because of the improved soil tillage during cultivation, fertilization, and irrigation along with the warmer temperatures and the absence of prolonged cold period.
The Meknes trial was characterized by relatively high biomass production. The DYH reached 179.3 q/ha for line A02 and 133.46 q/ha for A17. The results at Meknes in 2018-2019 were lower, though consistent with, those of the previous year in the same region.
The ranking of the genotypes changed from one environment to another. The genotypes yield fluctuation across the locations and years is an indication of the significant effects of the G×E interaction and a differential performance of the genotypes across the environments as well as yield instability between the experimental locations [12].

Three-Factor ANOVA's and CV's
Descriptive statistics for the agro-morphological and productivity traits showed different degrees of variation between the parameters (data not presented). BYDV susceptibility rate varied between 0 and 100% with an average of 20% and had the highest variability coefficient (129.56%). The RW presented a high coefficient of variation (119%). The GYH, DYH, DYP, GYP, NGP, and NFT all showed variability ranging from 60-119%. The GYH varied between 0.34 and 164.80 q/ha with a mean of 24.38 q/ha. The thousand-seed weight (TSW) had the lowest CV (27.46%), ranging from 12.50 to 58.49 g with a mean of 35.19 g.
We performed the three-way ANOVAs for the site, line, and block factors. Significant sources of variation are displayed in Table 9. Very highly significant site, line, and their interaction for the 13 agro-morphological traits are shown except for the effect of the accession on plant vigor in the vegetative growth phase. The block factor-related effects were high for PH, RW, and vegetative growth vigor and were not significant for most traits. The accession-block interaction was significant for PH and DYQ (Table 9). Table 9. Three-way ANOVA for the agro-morphological and productivity traits. Explanations: PH = plant height (cm); RL = root length (cm); NFT = fertile tiller number; RW = root weight (g); SW = stem weight (g); BP = plant biomass (g); NGP = number of grains per plant; DYP= dry matter weight per plant (g); GYH = grain yield (q/ha); DYH = dry matter yield (q/ha); HI = harvest index; TSW = thousand-seed weight (g); VGV = growth vigor at vegetative stage. * Significant at p = 0.05; ** significant at p = 0.01; *** significant at p = 0.001.

Correlation Matrix
The Pearson correlation matrix values among variables are displayed in Table 10. Significant values are highlighted in bold type. The highest positive correlations of paired variables were observed between DYP and GYP (R = 0.97) and between NGP and GYP (R = 0.93). Other significant positive correlations were observed between GYH and DYH (R = 0.65) and between DYP and RW (R = 0.60). There were also significant correlations with R-values of 0.54 between DYP and GYP; 0.52 between DYP and SW; and 0.51 between DYH and VGV, SW, and DYP (Table 10). The low positive correlations ranged from 0.49 to 0.40 and were observed between NGP and GYP (R = 0.47); NGP and GYH (R = 0.45); TSW and PH (R = 0.41); and HI and GYP (R = 0.40). The lowest correlations were measured between NFT and RL (R = 0.36); NFT and TSW (R = 0.36); NFT and GYP (R = 0.35); NGP and DYP (R = 0.33); and GYH and DYP (R = 0.31). In addition, there were two negative correlations between parameters: the first was between HI and DYH (R = −0.43) and the second between GYP and BYDV sensitivity (R = −0.34).

Principal Component Analysis
The principal component analyses (PCA) were conducted at three sites ( Table 2): Berrechid on the relatively humid central Atlantic Coastal Plain; Bouchane at a semiarid location on the Central Plateau; and El Kebab representing a relatively high-altitude environment in the Middle Atlas Mountains. In addition, we carried out a combined-site PCA.
The Berrechid experiment's first two principal axes' contribution to the total variability was 57.78%; for El Kebab and Bouchane these values were 57.03% and 56.14%, respectively (Table 11). The variables positively correlated to the first principal component in the three locations were GYP, NGP, NYP, and NFT. The GYH, DYH, and SW parameters contributed to the PC1 in two of the three locations. The PH and TSW contributed mainly to the 2nd and 3rd main axes. Unique significant contributions observed were for RW to the PC1 at Bouchane and HI to the PC2 at El Kebab. The variables of the eight experiments that contributed significantly to the first main axis were DYP, SW, GYP, and NGP (Table 11). The second axis correlated with TSW, DYH, VGV, and NGP, and the third main axis to HI, NFT, PH, and RL. The projection of individuals and variables onto the two first main axes using biplots revealed groups of lines with comparable field performance ( Figure 3). The Berrechid site biplot (a) shows a first cluster (green) linking lines that performed better in terms of GYH, DYH, NGP, NFT, and RW; these included lines A35, A27, A01, A30 A45, A44, and A14. The second group (blue) highlighted lines that were tall, having intermediate yields, and greater harvest indexes: A18, A05, A36, A03, A39, A43, and A28. The poorest performing lines in Berrechid clustered (yellow) in the opposite ends of the PC1; these were less productive, having low HI and PH. The last group (red)wase revealed at Berrechid to have significant VGV, producing intermediate biomass and low HI values; these included lines A40, A23, A16, A07, A25, and A17. The rest of the accessions performed intermediately well for most of the studied parameters.
In the semiarid and hot environment of Bouchane, the biplot delimited four clusters (Figure 3b). The first group grouped the highest-yielding lines: A20, A39, A33 A26, and A29. The second cluster in yellow grouped together drought-sensitive and poor-productivity lines, including A42, A40, A25, A11, A10, A12, A06, and A17. The two other clusters included lines that showed intermediate agronomical potential; the blue one on the top grouped lines A08, A16, A01, A02, A04, A23, A05, and A15 that had better root systems (RL and RW) and were fairly susceptible to BYDV. The other red group below the biplot center reassembles A35, A14, A18, A22, A30, A27, A38 A34, and A32 lines that have high HI, TSW, PH values and relative susceptibility to crown rust.
The factorial plane constructed from axes PC1 and PC2 (Figure 3c) revealed a strong grouping of lines towards the origin of the axis at the high-elevation site of El Kebab. Nevertheless, two distinct groups of experimental lines were evident. The group to the left of the origin (cluster in blue) consists of lines that were generally less productive, while those to the right of the plot were higher yielding. In the El Kebab first-year evaluation, genotypes with the highest grain yields, NGP, and NFT and HI values much higher than those of the experimental lines were A38, A36, A34, A39, A06, A17, A33, A17. Below this cluster, on the negative side of PC2, lines A10, A23, A40 produced more biomass and presented low harvest indices-factors indicative of good potential for ensilage. The factorial plane constructed from axes PC1 and PC2 (Figure 3c) revealed a strong grouping of lines towards the origin of the axis at the high-elevation site of El Kebab. Nevertheless, two distinct groups of experimental lines were evident. The group to the left of the origin (cluster in blue) consists of lines that were generally less productive, while those to the right of the plot were higher yielding. In the El Kebab first-year evaluation, genotypes with the highest grain yields, NGP, and NFT and HI values much higher than those of the experimental lines were A38, A36, A34, A39, A06, A17, A33, A17. Below this cluster, on the negative side of PC2, lines A10, A23, A40 produced more biomass and presented low harvest indices-factors indicative of good potential for ensilage.
The last biplot showing the eight combined experimental trials is in Figure 3d. Each line's potential corresponds to its average performance across testing environments. The green cluster on the positive side of PC1 and the negative side of the PC2 axis grouped lines having high yield potential. In the upper right side of axis 2, the red cluster highlights The last biplot showing the eight combined experimental trials is in Figure 3d. Each line's potential corresponds to its average performance across testing environments. The green cluster on the positive side of PC1 and the negative side of the PC2 axis grouped lines having high yield potential. In the upper right side of axis 2, the red cluster highlights lines that displayed good biomass potential through their vegetative vigor, high dry matter yield, and root and on-ground vegetation weight. The lines on the left side of PC1 had the lowest productivity and were relatively more sensitive to crown rust and BYDV.

AMMI Analysis
The objective of the AMMI analysis was to accurately characterize the genotype and environmental effects on the productivity and stability of the A. magna germplasm across environments. The combined analysis of variance for grain yield of the 41 Avena magna lines across the eight environments is in Table 12. The effects of the genotypes, the environments, and their interaction were highly significant (p ≤ 0.0001) on the total variance. The environment's main effect accounted for 82.89%, whereas genotype and G × E interaction effects accounted for 1.81% and 15.30% of the total variation, respectively. The excessive value of the sum of squares of the environment reflected the large differences among testing locations and years (Table 12). Consequently, the A. magna lines performed differently across the testing environments. The high percentage of the experimental location effects indicates that the major factor influencing yield was the environment. Even if the genotype factor showed minimal contribution to the total variance for grain yield, it exhibited highly significant genetic differences among domesticated oat lines for grain yield.
In Table 12, the AMMI multiplicative component partitioned the G × E interaction into three interaction principal component axes (IPCAs). The first two axes showed a significant contribution to the G × E in the AMMI model. The first interaction component explained 47.7% of the SCE of the G × E interaction for grain yield. The second axis explained an additional 31.7% of the SCE of the interaction. The AMMI biplot in Figure 4 accounted for 79.4% of the G × E interaction, providing the interaction principal component scores with 90 degrees of freedom.  Genotypes clustering closer to the intersection of the two PC axes tended to be more stable, while those that plotted farther apart from the origin had lower yield stability. According to Purchase [13] and Adugna and Labuschagne [20], the greater the IPCA score (positive or negative), the more adapted a genotype is to a specific environment; con- Genotypes clustering closer to the intersection of the two PC axes tended to be more stable, while those that plotted farther apart from the origin had lower yield stability. According to Purchase [13] and Adugna and Labuschagne [20], the greater the IPCA score (positive or negative), the more adapted a genotype is to a specific environment; conversely, the more the IPCA score approaches zero, the more stable it is. Accordingly, oat lines A35, A13, A27, and A14 were relatively unstable, while lines A06, A08, and A34 displayed high yield stability across environments.
The projection point of a genotype close to an environmental vector indicated specific interaction between both factors. Genotypes A13, A02, A39, A36, and A08 interacted positively at the El Kebab_19 environment. The genotype with the highest positive interaction with Berrechid_18 was A35; A13 interacted positively with El Kebab_19, while A42 had high interaction with Bouchane_19 and A16 was the best genotype in Meknes_18 (Figure 4).
When looking across environments, it is clear that there is tremendous variation among the different experimental locations. Berrechid_18, Meknes_18, and El Kebab_19 showed good potential for discriminating among genotypes, as indicated by their distance from the biplot origin ( Figure 4). However, because of their IPCA2 score, the data may not accurately represent their agronomic potential at a specific location. A closer relationship was observed among Meknes_19, Bouchane_19, Tiflet_19, and El Kebab_2018.
Examination of the factorial plane also allowed for the classification of environments according to their average productivity ( Figure 5) Figure 5. Plot depicting the relationship between AMMI stability (IPC1) and grain yield. IPCA1 (y-axis) represents genotype x environment (G × E) interaction plotted against grain yields in q/ha. Values closest to 0 on the y-axis represent lines with the greatest stability for yield.
As evidenced by their orientation on the right side of the factor plane, the set of lines including A35, A45, A40 and A37 presented the highest yields of the A. magna collection. In contrast, lines A42 and A22, located on the far-left side, were the least productive. However, some of the genotypes with the highest yields were the least stable, including lines Figure 5. Plot depicting the relationship between AMMI stability (IPC1) and grain yield. IPCA1 (y-axis) represents genotype × environment (G × E) interaction plotted against grain yields in q/ha. Values closest to 0 on the y-axis represent lines with the greatest stability for yield.
As evidenced by their orientation on the right side of the factor plane, the set of lines including A35, A45, A40 and A37 presented the highest yields of the A. magna collection. In contrast, lines A42 and A22, located on the far-left side, were the least productive. However, some of the genotypes with the highest yields were the least stable, including lines A35, A27, A13, and A14.
Although the common oat (A. sativa) would out-yield the set of A. magna experimental accessions, it is important to note that the 108.7 q/ha seed yield of A35 at Berrechid likely represents a substantially higher seed protein yield in comparison with any common oat variety, given the previously reported protein content in the Avery check (line A40) of 25.9% versus A. sativa at 12.9% [9]. In addition, considering that common oats are known to require more water for grain production than wheat or barley [21] it is remarkable that the A. magna accessions at Bouchane included lines approaching 28.86 q/ha (line A20) under a mere 201.5 mm (7.93 in) of ambient rainfall in well-drained soils ( Figure 2, Table 2) and without supplementary irrigation-conditions under which surrounding wheat and barley fields did not provide measurable seed yields. Varietal registration trials of check 'Avery' performed in 2019-20 recorded a grain yield of 32.43 q/ha with two supplemental irrigations at El Kebab (unpublished). These results warrant continued efforts to develop A. magna as a climate change-resistant crop for human and animal nutrition in The graph differentiates four groups with good performance and higher-than-average yield. Group one was the most stable (ASVs < 1) and included lines A32 (32.44 q/ha) and A34 (32.02 q/ha). Group two was moderately stable (1 < ASV < 2) and consisted of five lines: A05 (30.45 q/ha), A40 (32.85 q/ha), A33 (31.58 q/ha), A23 (30.56 q/ha), and A04 (30.01 q/ha). Group three was relatively unstable with ASVs between 2 and 4; it comprised the following 12 lines: A45 (33.93), A37 (32.77), A44 (32.70), A38 (32.34), A16 (32.09), A29 (31.98 q/ha), A02 (31.79 q/ha), A31 (30.72 q/ha), A01 (30.67 q/ha), A41 (30.64 q/ha), A18 (30.09 q/ha), and A14 (29.65 q/ha). Group four consisted of two lines having very high yields but with low stability: A35 (38.57 q/ha), and A27 (31.38 q/ha).
Although the common oat (A. sativa) would out-yield the set of A. magna experimental accessions, it is important to note that the 108.7 q/ha seed yield of A35 at Berrechid likely represents a substantially higher seed protein yield in comparison with any common oat variety, given the previously reported protein content in the Avery check (line A40) of 25.9% versus A. sativa at 12.9% [9]. In addition, considering that common oats are known to require more water for grain production than wheat or barley [21] it is remarkable that the A. magna accessions at Bouchane included lines approaching 28.86 q/ha (line A20) under a mere 201.5 mm (7.93 in) of ambient rainfall in well-drained soils ( Figure 2, Table 2)