1. Introduction
Barley (Hordeum vulgare L) is considered one of the founding crops of old-world agriculture. Its cultivation began 10,000 years ago in the “Fertile Crescent” (1). It is one of the most adapted crops that can be grown under diverse climatic and edaphic conditions due to its ability to escape late drought and its salinity tolerance compared to other cereal crops (2).
Barley was a staple food crop for several thousand years, extending from Ethiopia and Egypt across the southern coast of the Mediterranean Sea to Morocco (3). Barley has played a major role in providing food security in Morocco throughout history; large grain storage has been maintained to ensure security by the ruling dynasties since the beginning of the second millennium BC (4).
The proven health benefits of barley have been the main driver for the renewed interest in the use of barley as food. Its health claims include the potential to reduce the risk of cardiovascular disease by lowering blood cholesterol (5). Barley has been preferred not only for its nutritional profile but also for its high concentrations of several nutraceutical compounds that have been shown to play an important role in functional foods (6). Barley is also a great source of many phytochemicals, including phenolics, phytates, and tocols, whose activity reduces oxidative degenerative reactions in the human body as evidenced by clinical studies (7).
The major components of barley grain are starch, dietary fiber, and protein, the contents of which are affected by genotypic and environmental factors as well as their interactions (8). It has been reported that among cereals, barley has a greater amount of soluble dietary fiber, particularly mixed linkage β-D-glucan, known for its health benefits (8). β-D-glucans are polysaccharides found principally in the cell walls of the aleurone layer and endosperm in both barley and oat kernels; they are more concentrated in the endosperm in barley, while in oats they are more concentrated in the aleurone layer (6). Changes in β-glucan levels have been attributed most of the time to environmental conditions, but the genotypic effect is also important (6). Increased levels of β-D-glucan have been correlated with dry and hot environmental growing conditions and the opposite when the environmental growing conditions are moist, such as rains during ripening (6). The variation in β-D-glucan is genetically associated with starch and with hulless cultivars known to have higher β-glucan content (9). A great diversity is reported in barley kernel colors at maturity, white, yellow, blue, purple, green, brown, and black pigmentation, which is usually due to the phenolic compound and anthocyanin content that give a different pigment in the pericarp and aleurone layer (10).
The use of barley in the starch industry generates proteins as a by-product that could be used in the food industry (11). Barley kernels are composed of 10–20% of proteins, mainly prolamin hordeins, which represent 87% of storage proteins in the endosperm layer, and glutenins, which represent 23% (12). Barley proteins are considered a good source of essential and non-essential amino acids (12). The protein content is under genetic control, but its expression is affected by the environment (13). The protein content is increased when nitrogen is supplied and under the effects of drought and heat stress (13).
In addition to its high-soluble fiber content, barley has other constituents with nutritional benefits, including vitamins (B1, B3,…) and micronutrients (calcium, iron, magnesium, phosphorus, potassium, zinc, and copper) (14). Indeed, zinc and iron deficiencies along with iodine, selenium, and vitamin A deficiencies are among the “big five” micronutrient deficiencies in humans throughout the world, responsible for over 65% of childhood deaths (15). Two billion people throughout the world are affected by micronutrient malnutrition also known as “hidden hunger,” which is common in resource-poor parts of the world where alternative food is rare to supplement the predominant cereal-based food (14). Furthermore, the cereal grown in mineral-deficient soils is very poor in both concentration and bioavailability of iron and zinc (16). Abraha et al. (17) reported in Ethiopian barley that the uptake of micronutrients and their translocation into grain are heritable traits. Mamo et al. (14) have reported a high heritability of Fe, Zn, and Se contents under greenhouse conditions; however, Joshi et al. (18) reported a low heritability under field conditions. For wheat, the major source of phenotypic variation for micronutrients is explained by the genotype * environment interaction (19).
To overcome micronutrient deficiencies, many strategies have been implemented, including supplementation, dietary diversification, fortification, and biofortification of crop plants (20). Compared to the other techniques, biofortification is considered more sustainable as it aims to increase genetically the bioavailable mineral content of food crops (21). Tools have been used to develop micronutrient-enriched cereals using genetic and genomic approaches (22).
Biofortification is either achieved through an agronomic approach, or genetic selection using conventional plant breeding methods, or biotechnology techniques (23). The agronomic approach is achieved by the application of foliar micronutrient fertilizers to cereal crops, but other strategies could be applied, such as crop rotation, correcting soil alkalinity, and the introduction of beneficial soil microorganisms (24–26). Barley biofortification through fertilization was applied to increase its bioavailable zinc and iron contents in the grains (27).
Genetic resources of cereal crops, including their wild relatives stored in gene banks, are a valuable source of alleles used in the improvement of crops, including their nutritional attributes (21). Screening of landraces and wild relative species is essential to identifying sources of micronutrients and other nutritional compounds. In this regard, the accessions of Hordeum spontaneum are a great source of high contents of β-glucans and micronutrients (28). Many genomic regions have been associated with richness in β-glucans and micronutrients, most of which are linked to wild alleles (29).
In crop improvement, the efficiency of selection is hindered by the magnitude of genotype–environment interactions. Breeders are using multi-environment testing to select varieties with good stability in the expression of the desired trait. There are several statistical methods used to characterize the stability and adaptation of varieties (30, 31), including the additive main effects and multiplicative interaction (AMMI) analysis, which is used for the stability analysis of zinc and iron in maize (32, 33). However, the genotype main effect plus genotype-by-environment interaction (GGE) biplot is increasingly used as it allows for better visualization of the performance and stability of genotypes, the discrimination ability and representativeness of environments, and the correlations between environments (34).
The aims of this study are: (i) to capture the effects of genotypes, environments, and their interactions on nutritional and quality attributes of some Moroccan barley varieties and elite lines, (ii) to identify barley varieties with higher nutritional value for human consumption, and (iii) to recommend testing environments for efficient breeding for quality and nutritional attributes in barley.
2. Materials and methods
2.1. Plant material and field trial
Five Moroccan barley varieties (Hulless and hulled, six-row and two-row), along with two breeding lines (green and dark brown aleurone), were obtained from the National Barley Breeding Program of the National Agricultural Research Institute of Morocco (INRA) (Table 1).
Table 1. Characteristics of barley varieties and breeding lines used.
These entries were grown in two replications in a randomized complete block design at each of the four Moroccan locations (Marchouch, Annoceur, Jamaat Shaim, and Tassaout) during the 2016–2017 cropping season and at Marchouch, Annoceur, Jamaat Shaim, and Sidi El Aydi during the 2017–2018 cropping season. Table 2 summarizes the geographic coordinates, altitude, and the agro-climatic zones of the locations. The trials were planted between mid-November and mid-December except at Tassaout, where planting was done on January 5, and was conducted based on the recommended agronomic packages, including the use of fertilizers and weeding using both herbicides and manually. The grains were harvested in the four central rows (5 m2) of each experimental unit, consisting of six rows of 5 m length and 0.25 m row spacing.
Table 2. Coordinates and eco-geographic zones of different testing locations.
The soil characteristics across the studied sites exhibit distinct patterns. Jamaat Shaim and Sidi El Aydi sites are characterized by clay soil texture, while Marchouch has heavy clay soil. Conversely, the Annoceur and Tassaout sites exhibited a medium texture with sandy clay loam soils. The presence of a calcium content risk is evident in all soils except at the Marchouch site. In terms of pH levels, all soils tend to be slightly basic to basic, except for the soil at Jamaat Shaim, where the pH remains neutral. Importantly, none of the soils are saline. The soils at Jamaat Shaim and Marchouch sites have medium organic matter content, while those at the other stations are poorer. The soils at the Marchouch and Sidi El Aydi sites exhibit high contents of available phosphorus, while Annoceur, Jamaat Shaim, and Tassaout sites display medium contents. Notably, all sites feature high available potassium content (more than 300 ppm), except for the soil at the Marchouch site (110 ppm).
Meteorological data for the experimental locations have been collected from the National Institute for Agriculture Research (Figure 1). Cumulative annual precipitation ranged from 286.6 mm in the first growing season to 497.4 mm in the second growing season for Marchouch station, from 204.8 mm (from January 2017 to August 2017) to 712 mm (from September 2017 to June 2018) for Annoceur station, and from 249 mm in the first growing season to 240 mm during the second growing season for Jemaa Shaim station. Higher temperatures were observed between May and September at the three locations. Cumulative precipitation in the other experimental locations was 262.8 mm and 325 mm, respectively, for Tassaout in 2016–2017 and Sidi El Aydi in 2017–2018. Only the trials at Tassaout received two supplemental irrigations of 20 mm each at sowing and at late tillering.
Figure 1. Monthly mean, maximum, minimum temperatures, and sums of precipitation of the testing locations during the two growing seasons 2016–2018.
2.2. Grain sample preparation
Dried and cleaned grains were ground using a laboratory cyclone mill (Retsch Cyclone Mill Twister, Retsch, Haan, Germany) with a sieve of 1 mm. The grinding mill was cleaned before and after milling each sample to avoid cross-contamination. The moisture content of each sample was determined following the American Association of Cereal Chemists (AACC) approved method [44–40]-44–40.01 Moisture Modified Vacuum-Oven Method (35).
2.3. Beta glucan determination
The β-glucan content was determined according to a fluorometric method using a continuous flow analyzer (San ++ SKALAR, Breda, Netherlands). Prior to this determination, an acid extraction was performed using the procedure recommended by the European Brewery Convention (36). The extract was pumped and injected in the β-glucan module alongside the calcofluor reagent prepared previously and mixed, and then the mixture passed through the flow cell of the fluorimeter for the lecture. The signal was fed to a computer equipped with a specific software analysis system (FlowAccess 3.4.0, Skalar, Breda, Netherlands), and quantification was based on the calibration curve. The calibration curve is generated using known standards with known β-glucan concentrations. The software uses this curve to convert the measured signal to a corresponding β-glucan concentration in the sample.
2.4. Minerals’ determination
The determination of microelement content was performed according to a standard protocol (37). The elemental determination was carried out with a simultaneous multi-element inductively coupled plasma-optical emission spectroscopy (ICP-OES) (ICAP-OES Duo, ICAP 7600 Series. Thermo Fisher Scientific, Waltham, USA), which detects ~30 elements simultaneously, including iron, zinc, and selenium.
2.5. Thousand kernel weight and protein determination
The thousand kernel weight was evaluated by weighing samples of 1000 kernels from the harvest of each plot, and the protein content was determined on the whole grain using near-infrared reflectance (NIR) (NIRS DS 2500, FOSS, Hillerod, Denmark) spectroscopy according to the AACC approved method [39–10] (35).
2.6. Statistical analysis
The analyses of variance for combined environments (location × year) and for each of the eight environments were conducted. The descriptive statistics, including the mean, LSD (5%), and CV (%), and tests of normal distribution were conducted for each trait in each environment. Data were also subject to GGE analysis, which uses principal component analysis to explain the maximum amount of variance in the original data by a few linearly transformed uncorrelated components. The resulting scatter plot based on PC1 and PC2 allows the visualization of mega-environments, the mean and stability of traits for each genotype and its comparison with the ideal genotype, the discriminating ability and representativeness of the environments, and the correlations among genotypes and environments (34). For each trait, similar genotypes are positioned closely in the biplot; genotypes far from the origin have large genotype plus interaction effects; and genotypes close to the origin have average performance in all environments. Closely associated environments suggest redundant locations, while large gaps indicate the need to extend the test environments to other locations. The plot ranking environments allow the identification of the most discriminant environment, which is nearest to the smallest circle representing the ideal environment. The same analysis can be conducted to compare the genotypes evaluated with the ideal genotype. The cosine of the angle between the vectors of genotypes or environments to the origin approximates the correlations between them; the acute angle shows a strong positive correlation, the obtuse angle shows a negative correlation, and the right angle indicates no correlation. Pearson’s correlations were performed to assess relationships between traits. All statistical analyses were conducted using the GenStat software (19th Edition GenStat, VSNi, Hemel Hempstead, England) (38).
3. Results
3.1. Analysis of variance
The combined analysis of variance showed highly significant differences among entries for all measured traits except for zinc content (Table 3). Furthermore, the effect of the environment was highly significant for all traits except for β-glucan content. The interactions between entry × environment were highly significant for β-glucan, TKW, iron, zinc, and selenium and significant for protein content.
Table 3. Analysis of variance for quality attributes and TKW of barley varieties and breeding lines grown in eight environments.
Because the interactions between entry × environment were significant to highly significant, an analysis of variance was conducted for each trait in each environment, which allowed comparison among the entries and the derivation of descriptive statistics.
For β-glucan, the highest content of 11.58% was obtained for the variety Chifaa at the Sidi El Aydi station in 2017–2018. This variety has also given the highest contents in all environments except Tassaout and Jamaat Shaim in 2016–2017, where the variety Assiya has given similar contents (Table 4). At Annoceur, in 2016–2017, the β-glucan content varied from 2.78% for Adrar to 7.37% for Chifaa, with 5.94% for the breeding line BFH129; during the 2017–2018 season, Adrar showed 6.02% of β-glucan, while Chifaa showed 9.04%. At Marchouch, the ranges were 1.95 to 8.78% in 2016–2017 and from 3.94 to 10.84% in 2017–2018, respectively, for the varieties Amalou and Chifaa. At Jamaat Shaim, the same ranking of these two varieties was observed during the two consecutive seasons, with Amalou having the lowest content and Chifaa the highest. At Tassaout, during 2016–2017, Assiya and Chifaa showed the highest content, while the lowest values were given by the varieties Amalou and Adrar.
Table 4. β-glucan content mean, CV, and LSD for barley entries in different environments.
The protein content varied from 9.28% for the landrace Rabat 071 at Annoceur in 2017–2018 to 20.14 for the breeding line M9V5 at Marchouch in 2016–2017 (Table 5). The highest protein concentrations were registered at Marchouch, Jamaat Shaim, and Tassaout in 2016–2017 and at Marchouch and Sidi El Aydi in 2017–2018. Annoceur in 2017–2018 showed the lowest protein percentage. The breeding line M9V5 had the highest value in each of the eight environments and ranged from 12.38% to 20.14%, followed by the variety Assiya (10.10 to 18.07%) and the breeding line BFH129 (11.15 to 16.49%). In most environments, the varieties Adrar and Amalou had low protein content, while Rabat 071 and Chifaa, Assiya, and the breeding line BFH129 had intermediate protein content.
Table 5. Protein content mean, CV, and LSD for barley entries in different environments.
For the thousand kernel weight, the highest values were recorded in the varieties Rabat 071 (41.13 to 59.30 g) and Adrar (42.46 to 54.31 g), while the lowest values were given by the breeding line M9V5 (30.15 to 39.88 g) and the variety Assiya (27.87–41.02 g) over the eight environments (Table 6). The remaining entries have a medium grain weight. The highest TKW was recorded at Annoceur in 2017–2018, followed by Tassaout in 2016–2017 and Marchouch in 2017–2018, while the lowest values were recorded at Marchouch in 2016–2017 and Jamaat Shaim and Sidi El Aydi in 2017–2018.
Table 6. TKW mean, CV, and LSD for barley entries in different environments.
Regarding the microelements, large variations were observed between entries and between environments. For iron content, it ranged from 17.49 ppm at Sidi El Aydi in 2017–2018 to 51.41 ppm at Marchouch in 2016–2017 for the same variety, Amalou (Table 7). It appeared that the iron contents were higher during the 2016–2017 season than in the 2017–2018 season. Among the entries, the varieties Adrar and Assiya had the lowest Fe content in most environments, while the breeding line M9V5 had higher values in the three environments in the 2017–2018 season. The breeding line BFH129, the landrace Rabat 071, and the varieties Amalou and Chifaa showed large variations with high to medium contents in most environments. For zinc content, the entries did not show significant differences when averaged over all environments. However, their content differed in some environments (Table 7). At the Marchouch station, the ranges were 29.15 ppm for Adrar to 42.96 ppm for Rabat 071 in the 2016–2017 season and from 25.08 ppm for Amalou to 34.28 ppm for Rabat 071 in 2017–2018. At Jamaat Shaim station, only the Adrar variety showed the lowest content of 21.86 ppm compared to the breeding line BHF129 (30.10 ppm) in 2016–2017, while no significant differences were declared among entries in 2017–2018. At Annoceur station, no significant differences were found among entries in 2016–2017, while in 2017–2018, the variety Adrar has given the lowest content of 23.18 ppm, while the rest of the entries have shown contents ranging between 38.21 and 42.82 ppm. For Sidi El Aydi in 2017–2018, no significant differences were found among entries. At Tassaout in 2016–2017, higher Zn contents were shown by the breeding line M9V5 (46.53 ppm) and the variety Chifaa (43.15 ppm) when compared with the variety Amalou (33.24 ppm). For selenium content, no significant differences were found between entries at Annoceur in 2016–2017, while the variety Amalou (0.09 ppm) showed significantly less than Chifaa (0.21 ppm) in 2017–2018 at the same station (Table 7). At the Marchouch site, the varieties Assiya, Chifaa, and Amalou outperformed the varieties Rabat 071 and Adrar, while the two breeding lines were intermediate in 2016–2017. In 2017–2018, the varieties Rabat 071 and Assiya showed by far the lowest values of 0.05 and 0.07 ppm, respectively, followed by the breeding line M9V5 (0.13 ppm) and Adrar (0.18 ppm), Chifaa (0.20 ppm), and Amalou (0.21 ppm), while the highest content was shown by the breeding line BFH129. At Jamaat Shaim, no significant differences were declared among entries in 2016–2017; however, large differences were found in 2017–2018, with a range of 0.07 ppm for the variety Assiya to the highest content given by the breeding line BFH129 (0.23 ppm), followed by the varieties Adrar (0.21 ppm) and Chifaa (0.20 ppm). At Tassaout in 2016–2017, the lowest content was registered for the variety Adrar and the highest by Rabat 071, while the remaining entries ranged from 0.15 to 0.17 ppm. At the Sidi El Aydi site in 2017–2018, the lowest content was shown by the breeding line M9V5 (0.06 ppm), followed by BFH129 (0.13 ppm), while the three varieties Adrar, Amalou, and Rabat 071 had a content of 0.21 ppm, and the hulless varieties Assiya and Chifaa had 016 ppm.
Table 7. Iron, zinc, and selenium contents mean, CV, and LSD for barley entries in different environments.
3.2. Pearson’s correlations
Pearson’s correlations between traits were only highly significant between TKW and protein content (-0.544) and significant between iron and β-glucan contents with a negative correlation of (−0.192) and between iron and protein with a positive correlation of (0.278) (Table 8). The correlations among the micronutrients were not significant.
Table 8. Correlation between quality and nutritional traits measured for seven barley entries over eight environments in Morocco.
3.3. GGE biplot analysis
The GGE biplot analysis was conducted only for β-glucan and the micronutrient contents, and four different plots were selected for the purpose of this study: (i) the scatter plot of the “which-wins-where” pattern, which allows to identify different sectors by joining the winner genotype in each sector to identify different mega-environments and the best genotypes in each mega-environment; (ii) the ranking plot, which allows to evaluate the genotype performance and stability through the projection of the genotype performance on average environment coordinate (AEC) abscissa. The shorter the projection from the AEC is, the more stable the genotype is; (iii) the comparison plot, which allows to identify the environment with the most discriminating power among test environments; and (iv) the joint biplot, which allows to compare the variety Chifaa (recommended for human food) with the landrace Rabat 071, the widely grown landrace in Morocco.
The two principal components of the biplot analysis for β-glucan accounted for 87.79% of the total variation (Figure 2). The scatter plot identified two mega-environments, with all of Annoceur, Marchouch, Jamaat Shaim, and Sidi El Aydi in one and only Tassaout 2016–2017 environments in the other. The ranking biplot showed that the variety Chifaa has the highest β-glucan content in all environments with relatively good stability. The breeding line BFH129 had a higher content than the mean and a stable performance over the environments, while all the rest of the entries performed less than the average, with the varieties Assiya and Amalou being the least stable. The comparison plot showed that Marchouch’s location was the most discriminant site, as both related environments were within the inner circle, representing the ideal environment. The joint biplot shows that Chifaa outperformed Rabat 071.
Figure 2. Scatter plot, ranking biplot, joint biplot, and comparison biplot for β-glucan of barley entries evaluated in eight environments (ANN, Annoceur; JSH, Jamaat Shaim; TST, Tassaout; MCH, Marchouch; and SEA, Sidi El Aydi).
For iron content, the two first principal components explained 74.46% of the variation (Figure 3). The scatter plot identified three mega-environments with a vertex of four sectors, with Marchouch 2016–2017 and Jamaat Shaim 2016–2017 forming the first mega-environments, Marchouch 2017–2018 the second mega-environment, and the remaining environments forming the third mega-environment. The ranking plot showed that the varieties Adrar and Assiya performed less than the average genotype and were unstable, the variety Amalou had average content but was the most unstable, and the remaining entries performed more than the average genotype with more stable performances. The comparison biplot identified only Tassaout 2016–2017 as the most discriminant environment, while both Marchouch environments were the least discriminant for this trait. The variety Chifaa has a slightly higher content of iron than Rabat 071, which appeared more stable.
Figure 3. Scatter plot, ranking biplot, joint biplot, and comparison biplot for iron of barley entries evaluated in eight environments (ANN, Annoceur; JSH, Jamaat Shaim; TST, Tassaout; MCH, Marchouch; and SEA, Sidi El Aydi).
For zinc content, 75.74% of the variation was accounted for by the two first principal components (Figure 4). Two mega-environments were identified, one having only Annoceur 2017–2018 and the other having the rest of the environments. Among the entries, the varieties Adrar and Amalou had the lowest zinc content but were the most unstable, while the remaining entries had above-average genotypes, with only the breeding line M9V5 being the most unstable. Marchouch in 2016–2017 allowed us to discriminate better among the entries for this trait. Chifaa slightly outperformed the landrace Rabat 071.
Figure 4. Scatter plot, ranking biplot, joint biplot, and comparison biplot for zinc of barley entries evaluated in eight environments (ANN, Annoceur; JSH, Jamaat Shaim; TST, Tassaout; MCH, Marchouch; and SEA, Sidi El Aydi).
For selenium content, the two first principal components explained 71.77% of the variation (Figure 5). The scatter plot allowed us to identify five mega-environments; one for each of the environments Annoceur 2017–2018, Jamaat Shaim 2016–2017, and Tassaout 2016–2017; one mega-environment including Jamaat Shaim and Marchouch in 2017–2018; and fifth mega-environment including the rest of the environments as shown by the scatter biplot. The ranking biplot classified the entries M9V5, Assiya, and Rabat 071 as having below-average content of selenium, with M9V5 being the most unstable entry overall. The remaining entries had above-average content and showed less stability of performance; Chifaa had higher content than Rabat 071.
Figure 5. Scatter plot, ranking biplot, joint biplot, and comparison biplot for selenium of barley entries evaluated in eight environments (ANN, Annoceur; JSH, Jamaat Shaim; TST, Tassaout; MCH, Marchouch; and SEA, Sidi El Aydi).
4. Discussion
Micronutrient malnutrition, referred to as “hidden hunger,” affects over two billion people in the world, including those living under harsh conditions. Iron and zinc deficiencies are among the most common and are posing serious health challenges in developing countries, mainly for resource-poor farming communities. Morocco is known for its highest per capita consumption of barley as food in the world, mainly by the rural communities in the mountainous and arid areas. Due to the predominant cereal-based diets and the limitations of diversifying the diets, biofortification through breeding micronutrient-rich varieties has become a research priority in many crop breeding programs. Several studies have promoted genetic biofortification as a more cost-effective and sustainable approach compared to agronomic biofortification, food diversification, or supplements (39, 40). The HarvestPlus and Biofortification Challenge research programs in collaboration with CGIAR (The Consultative Group on International Agricultural Research) centers and other national partners were able to develop zinc-enriched varieties of maize, rice, and wheat; iron-enriched varieties of bean and pearl millet; and vitamin A-enriched varieties of cassava, maize, and sweet potato (https://www.harvestplus.org/home/crops). ICARDA and many countries, including Morocco, are working on the development of barley with high β-glucan, which allowed the release by INRA-Morocco of two hulless barley varieties, Assiya and Chifaa, with the latter having high β-glucan.
The results of the present study on five released barley varieties in Morocco and two breeding lines showed large genotypic differences in protein content, thousand kernel weight, β-glucan, and Fe and Se micronutrient contents, significant effects of environments except for β-glucan, and significant genotype × environment interactions for all traits, indicating the potential for selecting nutrient-enriched germplasm and the need for multi-environment testing for selecting more stable genotypes. Many authors have reported similar findings, highlighting the greater significance of genotypic effects over environmental effects on protein content (9, 41–43). A non-significant effect of environment is reported for Ethiopian barley for β-glucan and protein contents (17). A negative correlation between β-glucan content and accumulated temperature to 30°C and total precipitation during grain filling (13), as well as between β-glucan and protein contents in drought-stressed environments, was reported (17). Zhang et al. (13) reported that the variations in β-glucan and protein contents are mainly due to environmental factors. The Moroccan variety Chifaa has the highest levels of β-glucan in all environments, with a percentage higher than what is reported in other studies (44).
Several studies have shown that the protein content is mainly under genetic control, with a slicing effect of the environment (13) and a high genotype × environment interaction (45). In our study, the interactions between genotype × environment were significant for all traits, including protein content. Andersson and Börjesdotter (46) have reported that the protein content has significant correlations with rainfall and degree days.
In this study, among micronutrients, Zn appears to be influenced more by the environment, while Se showed a high influence of the environment and a high genotype × environment interaction. Significant genotypic and environmental effects and their interactions were also reported for micronutrient contents in wheat (47). In other studies, the variation of Fe and Zn contents in barley genotypes depended on the variation of environmental conditions (48), while they were under genetic control for Ethiopian barley (17). Xu et al. (32) reported that the concentrations of Fe and Zn in grains of maize under optimum conditions were higher compared to under low-N conditions. Among the micronutrients, Se was the least studied in barley despite its health benefits on thyroid function, and our study showed a high contribution of genotype × environment interactions to its variation, making the selection of this trait more difficult.
Thousand kernel weight is an important yield component and quality attribute of barley that is most appreciated by farmers. This study showed the predominance of genetic effects and high genotype × environment interactions. Similar conclusions were reported by Boudiar et al. (49) and Bocianowski et al. (50); however, Kumar et al. (51) did not report significant genotype × environment interactions for barley grown in India. Among all the traits evaluated, our results showed the highest negative correlation between TKW and protein content, while Fe showed the highest positive correlation with β-glucan and a negative correlation with protein content, but no significant correlations were found between the micronutrient contents. Other studies have reported significant correlations between Zn and Fe in Syrian barley (52), in wheat (53), and with protein content in wheat (54). β-glucan was reported to be correlated with TKW in oats (46). Based on these results, selection should be conducted independently for the micronutrients and for β-glucan, and the negative correlation between TKW and protein content should be considered during the selection. The identification of molecular markers could facilitate the selection for quality and nutritional traits if associated QTLs and markers are identified, as reported in barley for protein and β-glucan contents (55, 56).
Morocco lies within the Mediterranean region, known for high climate variability within and between the seasons. Therefore, a comprehensive evaluation of GE interactions requires more advanced statistical methods than the standard ANOVA. Since the genotype × environment interactions were significant to highly significant for all traits measured, selection and identification of lines with higher β-glucan and higher micronutrient contents will be facilitated by the identification of locations with the most discriminating power for plant selection during segregating populations and by the evaluation of the performance stability of the derived fixed lines using multi-environment testing. GGE biplot was used since it offers many features not shown in other stability analysis approaches, GGE biplot provides many visual interpretations, including crossover G×E interaction, comparing tested genotypes with ideal genotypes, and identifying environments with the most discriminating ability (57). In the GGE biplot, the first component is highly correlated with genotype performance, while the second component allows for assessing the variation of the trait due to genotype × environment. For β-glucan, most of the environments could be used to identify genotypes with high content, but the Marchouch station appears more appropriate for efficient selection for this trait in early and later generations. Additional locations will be needed to fit the two additional sectors identified through the scatter plot. Among the tested genotypes, the hulless barley variety Chifaa has the highest β-glucan content and can express its high β-glucan when grown in all environments. For iron content, except for Tassaout in 2016–2017, none of the environments has high discriminating power, and Marchouch station has the least discriminating ability. Among the genotypes tested, the varieties Assiya and Adrar showed the lowest mean content and highest instability and therefore cannot be used as biofortified varieties for this trait. However, the variety Chifaa and the hulless breeding line M9V5 have above-average iron content. For zinc content, all environments were grouped in one sector except for Annoceur 2017–2018, and none of the environments covered the remaining two sectors. Based on the ranking biplot, only the varieties Adrar and Amalou have low zinc content, while Chifaa and the breeding line M9V5 have above-average content but less stable content. Among the locations, Marchouch appears to be the most appropriate for the selection and evaluation of this trait. For selenium, all the genotypes tested showed high instability, but selection for high content could be done at Marchouch station.
Based on the results of the analysis of variance and stability analysis using the GGE biplot, the hulless barley varieties Chifaa and Assiya can be recommended for human consumption, with the former one having the advantage of being biofortified for β-glucan with acceptable levels of Fe, Zn, and Se contents. The two hulless breeding lines BFH129 and M9V5 can also be promoted as varieties for human consumption or used as sources of micronutrients to be crossed with Chifaa for the development of varieties combining high β-glucan and high content of micronutrients. For micronutrients, a special search for parental germplasm with high contents is needed to develop biofortified varieties. For the remaining varieties, the landrace Rabat 071 and the newly released varieties Amalou and Adrar can be considered for multiple purposes, but their use as human food will require a dehulling process which could reduce the micronutrients. Among the experimental stations used, Marchouch appears to be the most appropriate for selecting β-glucan glucans and for Se and Zn but not for Fe. As a breeding implication, this station can be used to select within segregating populations, but multilocational trials are needed for most advanced lines to select the ones combining high-quality and more stable varieties for quality and nutritional traits. For commercial production of biofortified varieties, it will be good to consider environments that allow the best expression of these nutritional traits. More research is needed on the stability and bioavailability of these nutrients during processing and consumption.
5. Conclusion
Barley can be further promoted as a human food to partially replace the imported wheat grain. This could allow consumers to benefit from the nutritional value as well as the health claims attributed to barley. In this study, we characterized five barley varieties released in Morocco and two hulless breeding lines for their protein, β-glucan, Fe, Zn, and Se contents. We could identify genotypes suitable for human consumption and recommend the hulless variety Chifaa as a β-glucan biofortified variety for cultivation in Morocco and as a parental germplasm for breeding germplasm combining high β-glucan with high micronutrient contents using Marchouch as selection site for early generations, followed by a multi-environment evaluation of the fixed lines.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
AA, MI, and FE: conceptualization. FE, AJ, and ZK: methodology. FE and ZK: software. AA and MI: validation. FE, ZK, and AA: statistical analysis. FE, AJ, and AA: fieldwork. FE and AE-B: laboratory analysis. FE and AA: writing draft manuscript. AA, MI, and ZK: editing and finalizing. All authors have read and agreed to the submitted version of the manuscript.
Funding
This research was partially funded by the Deutsche Gesellschaft für Internationale Zusammenarbeit (GIZ) Genebank Uplift Funding to ICARDA.
Acknowledgments
The authors would like to acknowledge the support of Dr. Khaled El-Shamaa for the statistical advice and help with the analysis of data and thank the heads of the INRA experimental stations of Annoceur, Jamaat Shaim, Marchouch, Sidi El Aydi, and Tassaout for helping with fieldwork.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Badr A, Sch R, Rabey HE, Effgen S, Ibrahim HH, Pozzi C, et al. On the origin and domestication history of barley (Hordeum vulgare). Mol Biol Evol. (2000) 17:499–510. doi: 10.1093/oxfordjournals.molbev.a026330
2. Grando S, Macpherson HG. Food barley: importance, uses and local knowledge. In: Proceedings of the International Workshop on Food Barley Improvement, Hammamet, Tunisia, 14–17 January. (2005).
3. Newman RK, Newman CW. Barley for Food and Health: Science, Technology, and Products. New York, NY: John Wiley and Sons (2008). doi: 10.1002/9780470369333
4. Saidi S, Lemtouni A, Amri A, Moudden M. Use of barley grain for food in Morocco. Food Barley: Importance, Uses and Local Knowledge. 17 (2005).
5. Behall KM, Scholfield DJ, Hallfrisch J. Lipids significantly reduced by diets containing barley in moderately hypercholesterolemic men. J Am Coll Nutr. (2004) 23:55–62. doi: 10.1080/07315724.2004.10719343
6. Ehrenbergerová J, Brezinová Belcredi N, Psota V, Hrstková P, Cerkal R, Newman CW, et al. Changes caused by genotype and environmental conditions in beta-glucan content of spring barley for dietetically beneficial human nutrition. Plant Foods Hum Nutr. (2008) 63:111. doi: 10.1007/s11130-008-0079-7
8. Åman P, Newman CW. Chemical composition of some different types of barley grown in Montana, USA. J Cereal Sci. (1986) 4:133–41. doi: 10.1016/S0733-5210(86)80016-9
9. Zhang G, Junmei W, Jinxin C. Analysis of β-glucan content in barley cultivars from different locations of China. Food Chem. (2002) 79:251–4. doi: 10.1016/S0308-8146(02)00127-9
10. Glagoleva AY, Shmakov NA, Shoeva OY, Vasiliev GV, Shatskaya NV, Börner A, et al. Metabolic pathways and genes identified by RNA-seq analysis of barley near-isogenic lines differing by allelic state of the Black lemma and pericarp (Blp) gene. BMC Plant Biol. (2017) 17:182. doi: 10.1186/s12870-017-1124-1
11. Yalçin E, Çelik S, Ibanoglu E. Foaming properties of barley protein isolates and hydrolysates. Eur Food Res Technol. (2008) 226:967. doi: 10.1007/s00217-007-0618-8
12. Sullivan P, Arendt E, Gallagher E. The increasing use of barley and barley by-products in the production of healthier baked goods. Trends Food Sci Technol. (2013) 29:124–34. doi: 10.1016/j.tifs.2012.10.005
13. Zhang G, Chen J, Wang J, Ding S. Cultivar and environmental effects on (1 → 3, 1 → 4)-β-D-glucan and protein content in malting barley. J Cer Sci. (2001) 34:295–301. doi: 10.1006/jcrs.2001.0414
14. Mamo BE, Barber BL, Steffenson BJ. Genome-wide association mapping of zinc and iron concentration in barley landraces from Ethiopia and Eritrea. J Cereal Sci. (2014) 60:497–506. doi: 10.1016/j.jcs.2014.08.007
16. Cakmak I, Ozkan H, Braun HJ, Welch RM, Romheld V. Zinc and iron concentrations in seeds of wild, primitive, and modern wheats. Food Nutr Bull. (2000) 21:401–3. doi: 10.1177/156482650002100411
17. Abraha A, Uhlen AK, Abay F, Sahlstrøm S, Bjørnstad Å. Evaluation of Ethiopian Barleys Genetic Variation for Agronomic and Nutritional Quality Traits in Northern Ethiopia (2018).
18. Joshi AK, Crossa J, Arun B, Chand R, Trethowan R, Vargas M, et al. Genotype × environment interaction for zinc and iron concentration of wheat grain in eastern Gangetic plains of India. Field Crops Res. (2010) 116:268–77. doi: 10.1016/j.fcr.2010.01.004
19. Ficco DBM, Riefolo C, Nicastro G, Simone D, Di Gesù V, Beleggia AM, et al. Phytate and mineral elements concentration in a collection of Italian durum wheat cultivars. Field Crops Res. (2009) 111:235–42. doi: 10.1016/j.fcr.2008.12.010
20. Rawat N, Neelam K, Tiwari VK, Dhaliwal HS. Biofortification of cereals to overcome hidden hunger. Plant Breeding. (2013) 132:437–45. doi: 10.1111/pbr.12040
21. Sakellariou M, Mylona PV. New uses for traditional crops: the case of barley biofortification. Agronomy. (2020) 10:1964. doi: 10.3390/agronomy10121964
22. Uauy C, Distelfeld A, Fahima T, Blechl A, Dubcovsky J. A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science. (2006) 314:1298–301. doi: 10.1126/science.1133649
24. Rana A, Joshi M, Prasanna R, Shivay YS, Nain L. Biofortification of wheat through inoculation of plant growth promoting rhizobacteria and cyanobacteria. Eur J Soil Biol. (2012) 50:118–26. doi: 10.1016/j.ejsobi.2012.01.005
25. Venkatesh MS, Hazra KK, Ghosh PK, Khuswah BL, Ganeshamurthy AN, Ali M, et al. Long–term effect of crop rotation and nutrient management on soil–plant nutrient cycling and nutrient budgeting in Indo–Gangetic plains of India. Arch Agron Soil Sci. (2017) 63:2007–22. doi: 10.1080/03650340.2017.1320392
26. Khoshgoftarmanesh AH, Norouzi M, Afyuni M, Schulin R. Zinc biofortification of wheat through preceding crop residue incorporation into the soil. Eur J Agronomy. (2017) 89:131–9. doi: 10.1016/j.eja.2017.05.006
27. Palmgren MG, Clemens S, Williams LE, Krämer U, Borg S, Schjørring JK, et al. Zinc biofortification of cereals: problems and solutions. Trends Plant Sci. (2008) 13:464–73. doi: 10.1016/j.tplants.2008.06.005
28. Elouadi F, Amri A, El-baouchi A, Kehel Z, Salih G, Jilal A. Evaluation of a Set of Hordeum vulgare subsp. spontaneum accessions for β-glucans and microelement contents. Agriculture. (2021) 11:950. doi: 10.3390/agriculture11100950
29. Wiegmann M, Thomas WT, Bull HJ, Flavell AJ, Zeyner A, Peiter E, et al. Wild barley serves as a source for biofortification of barley grains. Plant Sci. (2019) 283:83–94. doi: 10.1016/j.plantsci.2018.12.030
30. Annicchiarico P. Defining Adaptation Strategies and Yield-Stability Targets in Breeding Programmes. Quantitative Genetics, Genomics and Plant Breeding. Wallingford: CABI Publishing Wallingford UK (2002), 365–83.
31. Mohammadi R, Amri A. Analysis of genotype × environment interaction in rain-fed durum wheat of Iran using GGE-biplot and non-parametric methods. Can J Plant Sci. (2012) 92:757–70. doi: 10.4141/cjps2011-133
32. Xu J, Wang X, Zhu H, Yu F. Maize genotypes with different zinc efficiency in response to low zinc stress and heterogeneous zinc supply. Front Plant Sci. (2021) 12:736658. doi: 10.3389/fpls.2021.736658
33. Singhal T, Satyavathi CT, Singh SP, Kumar A, Sankar SM, Bhardwaj C, et al. Multi-environment quantitative trait loci mapping for grain iron and zinc content using bi-parental recombinant inbred line mapping population in pearl millet. Front Plant Sci. (2021) 12:659789. doi: 10.3389/fpls.2021.659789
34. Yan W. Mega-environment analysis and test location evaluation based on unbalanced multiyear data. Crop Sci. (2015) 55:113–22. doi: 10.2135/cropsci2014.03.0203
35. AACC. Approved Methods of the American Association of Cereal Chemists, 10th Edn. St. Paul, MN: Cereals and Grains Assn (2000).
36. Manzanares P, Sendra JM. Determination of total (1 → 3), (1 → 4)-β-D-glucan in barley and malt flour samples. J Cereal Sci. (1996) 23:293–6. doi: 10.1006/jcrs.1996.0030
37. Pequerul A, Pérez C, Madero P, Val J, Monge E. A rapid wet digestion method for plant analysis. Optim Plant Nutr. (1993) 31:3–6. doi: 10.1007/978-94-017-2496-8_1
39. Garg M, Sharma N, Sharma S, Kapoor P, Kumar A, Chunduri V, et al. Biofortified crops generated by breeding, agronomy, and transgenic approaches are improving lives of millions of people around the world. Front Nutr. (2018) 5:12. doi: 10.3389/fnut.2018.00012
40. Goredema-Matongera N, Ndhlela T, Magorokosho C, Kamutando CN, van Biljon A, Labuschagne M, et al. Multinutrient biofortification of maize (Zea mays L) in Africa: current status, opportunities and limitations. Nutrients. (2021) 10:1039. doi: 10.3390/nu13031039
41. Holtekjølen AK, Uhlen AK, Bråthen E, Sahlstrøm S, Knutsen SH. Contents of starch and non-starch polysaccharides in barley varieties of different origin. Food Chem. (2006) 94:348–58. doi: 10.1016/j.foodchem.2004.11.022
42. Hang A, Obert D, Gironella AIN, Burton CS. Barley amylose and β-glucan: their relationships to protein, agronomic traits, and environmental factors. Crop Sci. (2007) 47:1754–60. doi: 10.2135/cropsci2006.06.0429
43. Griffey C, Brooks W, Kurantz M, Thomason W, Taylor F, Obert D, et al. Grain composition of Virginia winter barley and implications for use in feed, food, and biofuels production. J Cereal Sci. (2010) 51:41–9. doi: 10.1016/j.jcs.2009.09.004
44. Nair S, Ullrich SE, Blake TK, Cooper B, Griffey CA, Hayes PM, et al. Variation in kernel hardness and associated traits in US barley breeding lines. Cereal Chem. (2010) 87:461–6. doi: 10.1094/CCHEM-12-09-0163
45. Jeuffroy MH, Casadebaig P, Debaeke P, Loyce C, Meynard JM. Agronomic model uses to predict cultivar performance in various environments and cropping systems. Agron Sust Dev. (2014) 34:121–37. doi: 10.1007/s13593-013-0170-9
46. Andersson AA, Börjesdotter D. Effects of environment and variety on content and molecular weight of β-glucan in oats. J Cereal Sci. (2011) 54:122–8. doi: 10.1016/j.jcs.2011.03.003
47. Alomari DZ, Eggert K, Von Wirén N, Polley A, Plieske J, Ganal MW, et al. Whole-genome association mapping and genomic prediction for iron concentration in wheat grains. Int J Mol Sci. (2018) 20:76. doi: 10.3390/ijms20010076
48. Ghandilyan A, Ilk N, Hanhart C, Mbengue M, Barboza L, Schat H, et al. A strong effect of growth medium and organ type on the identification of QTLs for phytate and mineral concentrations in three Arabidopsis thaliana RIL populations. J Exp Bot. (2009) 60:1409–25. doi: 10.1093/jxb/erp084
49. Boudiar R, Casas AM, Cantalapiedra CP, Gracia MP, Igartua E. Identification of quantitative trait loci for agronomic traits contributed by a barley (Hordeum vulgare) Mediterranean landrace. Crop Past Sci. (2016) 67:37–46. doi: 10.1071/CP15149
50. Bocianowski J, Warzecha T, Nowosad K, Bathelt R. Genotype by environment interaction using AMMI model and estimation of additive and epistasis gene effects for 1000-kernel weight in spring barley (Hordeum vulgare L.). J Appl Genetics. (2019) 60:127–35. doi: 10.1007/s13353-019-00490-2
51. Kumar D, Narwal S, Verma RPS, Kharub AS. Genotypic and growing location effect on grain β-glucan content of barley under sub-tropical climates. Indian J Genet Pl Br. (2017) 77:235–41. doi: 10.5958/0975-6906.2017.00031.1
52. El-Haramein FJ, Grando S. Determination of iron and zinc content in food barley. In: Proceedings of the 10th International Barley Genetics Symposium, Alexandria: Egypt. (2010).
53. Genç Y, Verbyla AP, Torun AA, Cakmak I, Willsmore K, Wallwork H, et al. Quantitative trait loci analysis of zinc efficiency and grain zinc concentration in wheat using whole genome average interval mapping. Plant Soil. (2009) 314:49–66. doi: 10.1007/s11104-008-9704-3
54. Zhao FJ, Su YH, Dunham SJ, Rakszegi M, Bedo Z, McGrath SP, et al. Variation in mineral micronutrient concentrations in grain of wheat lines of diverse origin. J Cereal Sci. (2009) 49:290–5. doi: 10.1016/j.jcs.2008.11.007
55. Kling JG, Hayes PM, Ullrich SE. Barley: genetics and breeding. In:Wrigley C, Corke H, Walker C, , editors. Encyclopedia of Grain Science. Elsevier Limited (2004). p. 27–38.
56. Emebiri LC, Moody DB, Horsley R, Panozzo J, Read BJ. The genetic control of grain protein content variation in a doubled haploid population derived from a cross between Australian and North American two-rowed barley lines. J Cereal Sci. (2005) 41:107–14. doi: 10.1016/j.jcs.2004.08.012