Next Article in Journal
Physiological and Full-Length Transcriptome Analyses Reveal the Dwarfing Regulation in Trifoliate Orange (Poncirus trifoliata L.)
Next Article in Special Issue
Genetic Diversity and Population Structure of Cowpea (Vigna unguiculata (L.) Walp.) Landraces from Portugal and Mozambique
Previous Article in Journal
Applications of Metabolomics for the Elucidation of Abiotic Stress Tolerance in Plants: A Special Focus on Osmotic Stress and Heavy Metal Toxicity
Previous Article in Special Issue
Genetic and Pomological Determination of the Trueness-to-Type of Sweet Cherry Cultivars in the German National Fruit Genebank
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variation of Zinc and Iron Concentration in Normal, Provitamin A and Quality Protein Maize under Stress and Non-Stress Conditions

by
Nakai Goredema-Matongera
1,2,3,*,
Thokozile Ndhlela
2,
Angeline van Biljon
3,
Casper N. Kamutando
4,
Jill E. Cairns
2,
Frederic Baudron
2 and
Maryke Labuschagne
3
1
Scientific and Industrial Research and Development Centre, 1574 Alpes Road, Harare 00263, Zimbabwe
2
International Maize and Wheat Improvement Centre, P.O. Box MP163, Mt Pleasant, Harare 00263, Zimbabwe
3
Department of Plant Sciences, University of the Free State, P.O. Box 339, Bloemfontein 9300, South Africa
4
Department of Plant Production Sciences and Technologies, University of Zimbabwe, P.O. Box MP167, Mt Pleasant, Harare 00263, Zimbabwe
*
Author to whom correspondence should be addressed.
Plants 2023, 12(2), 270; https://doi.org/10.3390/plants12020270
Submission received: 8 November 2022 / Revised: 22 December 2022 / Accepted: 22 December 2022 / Published: 6 January 2023

Abstract

:
The negative impacts of zinc (Zn) and iron (Fe) deficiency due to over-reliance on monotonous cereal-based diets are well-documented. Increasing micronutrient densities in maize is currently among top breeders’ priorities. Here, 77 single-cross Zn-enhanced hybrids with normal, provitamin A and quality protein maize genetic backgrounds were evaluated together with seven checks for grain Zn and Fe concentration and agronomic traits under optimum, low nitrogen (N) and managed drought conditions. Results showed a fairly wide variability for grain Zn (10.7–57.8 mg kg−1) and Fe (7.1–58.4 mg kg−1) concentration amongst the hybrids, across management conditions. Notable differences in Zn concentration were observed between the Zn-enhanced quality protein maize (QPM) (31.5 mg kg−1), Zn-enhanced provitamin A maize (28.5 mg kg−1), Zn-enhanced normal maize (26.0 mg kg−1) and checks (22.9 mg kg−1). Although checks showed the lowest micronutrient concentration, they were superior in grain yield (GY) performance, followed by Zn-enhanced normal hybrids. Genotypes grown optimally had higher micronutrient concentrations than those grown under stress. Genotype × environment interaction (G × E) was significant (p ≤ 0.01) for GY, grain Zn and Fe concentration, hence micronutrient-rich varieties could be developed for specific environments. Furthermore, correlation between grain Zn and Fe was positive and highly significant (r = 0.97; p ≤ 0.01) suggesting the possibility of improving these traits simultaneously. However, the negative correlation between GY and grain Zn (r = −0.44; p ≤ 0.01) and between GY and grain Fe concentration (r = −0.43; p ≤ 0.01) was significant but of moderate magnitude, suggesting slight dilution effects. Therefore, development of high yielding and micronutrient-dense maize cultivars is possible, which could reduce the highly prevalent micronutrient deficiency in sub-Saharan Africa (SSA).

1. Introduction

Grain production for major staple cereals such as maize, wheat and rice has been increasing due to advances in research and technology, with the ultimate goal of meeting the food demand for the continuously growing population [1]. While this improves food self-sufficiency at household and international levels, low densities of micronutrients such as Zn, Fe and vitamin A in these crops remain a global challenge and contribute significantly to micronutrient deficiency [2,3]. Micronutrient deficiency, also referred to as “hidden hunger” is detrimental to health, causing an array of ailments in populations whose diets are characterized by low micronutrient intakes [4,5]. Global statistics estimate that over two billion people suffer from micronutrient deficiency worldwide, although the highest prevalence rate occurs in low income countries where at least half of the population is affected [6,7]. Among the micronutrients, Zn and Fe deficiencies are currently considered as major human health risk factors and cause of stunting, mental illness, anaemia, diarrhoea and reduced ability to do physical work in both children and adults [8,9]. Several studies reported that more than 60% of all infant mortalities are associated with inadequate dietary intake of Zn and Fe [10,11]. All this evidence shows that though required in minute quantities, micronutrient deficiency can cause huge socio-economic losses through reduced national production and increased health care costs [3]. Despite a remarkable reduction in prevalence of various micronutrient deficiencies in most developed countries over the past decades, Zn and Fe deficiencies are still widespread in the developing regions such as in sub-Sahara Africa (SSA) and South Asia. Populations in these regions are worst affected due to overreliance on cereals as a food source, with limited nutritional supplements as well as ignorance of good nutritional practices [12,13].
One way of curbing micronutrient deficiencies in developing maize-based regions is to develop micronutrient-dense maize varieties through genetic and agronomic approaches [11,14]. Improving the nutritional quality of maize has always been a priority for breeders, and this strategy complements other micronutrient deficiency mitigation measures, such as dietary diversification, clinical supplementation and industrial fortification [2,15]. To date, several biofortified maize varieties with improved protein quality and high provitamin A and Zn content have been developed and commercialized worldwide, through concerted efforts from both national and international research organizations [2,16]. However, breeding for maize varieties with high levels of grain Zn and Fe concentrations is still lagging behind. To our knowledge, Zn and Fe enhanced maize varieties are still very few, following the recent commercialization of Zn-enhanced ICTA HB-18, ICTA HB-15 and BIO-MZN01 hybrids in Guatemala and Colombia, respectively [17].
Breeding for genotypes with high levels of grain Zn and Fe concentration depends on the existence of sufficient genetic variability for such traits in germplasm. Significant genetic variation in grain Zn and Fe concentration has been reported [4,18]. In addition, some studies reported positive association between grain Zn and Fe concentration and between grain yield and micronutrient concentration, suggesting the possibility of simultaneous improvement of such traits [4,19]. Therefore, knowledge of the nature of phenotypic correlations between traits is important in plant breeding and facilitates indirect selection of such traits. Moreover, the ability to accumulate significant concentrations of Zn and Fe in edible crop parts differs among maize genotypes, and sometimes depends on the root microbiome and soil physicochemical properties such as total soil Zn and Fe, extractable Zn and Fe, soil pH, organic matter and moisture content [20,21]. The soil must have sufficient plant-available Zn and Fe concentration for absorption by the maize roots. Among the soil chemical factors, soil pH plays the most important role in Zn and Fe solubility in soil solution. Previous studies reported about 30 to 50% decrease in Zn concentration in soil solution for each unit increase in pH [7,14]. Hence, a slight increase of soil pH from 5.5 can increase the risk of Zn deficiency in maize. Other studies reported sufficiently high Zn2+ in soil solution when pH is maintained at 5.0 [14]. In addition, soil organic matter increases solubility and transport of Zn to plant roots [14]. Low grain Zn and Fe concentrations in maize can be caused by moisture stress. Transport of these micronutrients to the root surface in soils occurs predominantly via diffusion, and hence adequate soil moisture is critical for providing suitable medium for that process [20].
Whilst several studies have focused on assessing the genetic variability of maize inbred lines and hybrids [1,4,13], very few studies have focused on how maize genotypes with different nutritional background, soil physicochemical properties, and crop management conditions can affect the accumulation of Zn and Fe in maize grain. This information is very useful in optimizing both agronomic and genetic biofortification as interventions to combat micronutrient deficiencies in affected regions. Therefore, the objectives of this study were: (i) to determine the variation of total soil extractable Zn and Fe and other physicochemical properties in maize experimental sites, (ii) to determine the variation of grain Zn and Fe concentration and other agronomic traits and its association with soil extractable Zn and Fe, (iii) to examine correlation between grain Zn and Fe and agronomic traits using Zn-enhanced normal (NML), provitamin A (ProA) and QPM genotypes grown under optimum, low N and managed drought conditions.

2. Results

2.1. Soil Chemistry

The summary of the descriptive statistical data for soil physiochemical properties analysed in soil samples that were taken in four experimental stations across two years is presented in Table 1. All the sites showed highly significant differences (p ≤ 0.05) for all the soil physicochemical properties except for inorganic carbon (IC). Soil pH ranged from moderately acidic (5.70–5.80), slightly acidic (6.10–6.40) to slightly alkaline (7.60) across sites and years (Table 1). The lowest and highest pH was observed at ART farm and Chisumbanje, respectively. Total carbon (TC) was the highest at CIMMYT under optimum management conditions and lowest at DR&SS under low N conditions. Nitrogen content (TN) was similar for all sites except for soil from DR&SS (low N site), which had <5% N content (Table 1). A similar trend was also observed for soil organic carbon content (SOC), which is the primary component of soil organic matter. Results showed significant differences (p ≤ 0.05) in available P (Olsen P extractable) content, with the highest value observed in soils from ART farm (25.9 mg kg−1) and lowest at CIMMYT (10.4 mg kg−1). Cation exchange capacity (CEC), a useful indicator of soil fertility ranged from 14.2–28.1 and 12.1–23.4 meq 100 g−1 in 2018/19 and 2019/20, respectively (Table 1). Extractable Zn concentration was lower than for Fe in all sites. Soils from ART farm had the highest soil Zn concentration with 5.10 and 5.70 mg kg−1 in 2018/19 and 2019/20, respectively, while Chisumbanje had the lowest values across the two years. In terms of soil texture, the highest clay content was observed in Chisumbanje, whilst DR&SS had proportionally higher sand content than other sites.

2.2. Agronomic and Micronutrient Performance of Zn-Enhanced Hybrids

The means for days to mid-anthesis (AD) for the hybrids did not differ much among the three management conditions (Table 2). However, the means for ASI showed notable differences. Under optimum management, means for anthesis-silking-interval (ASI) were ≤2.00 days across sites and years, while slightly higher mean values were observed under low N and managed drought conditions. Plants were shortest in E8 (161 cm) under managed drought conditions, and tallest in E2 (245 cm) under optimum conditions. Similarly, the mean grain yield performance of hybrids was highest in E6 (9.30 t ha−1) under optimum conditions, and lowest in E8 (1.10 t ha−1) under managed drought conditions. Combined analysis for GY indicated that most hybrids performed better under optimum than stress conditions (Figure 1).
The means for grain Zn concentration were higher than for Fe in most sites. The highest mean for grain Zn was observed in E1 (30.7 mg kg−1) followed by E6 (30.4 mg kg−1) which were optimum sites, whereas stress induced sites E3 and E8 had the lowest grain Zn means of 24.8 mg kg−1 and 25.1 mg kg−1, respectively. Grain Fe concentration was lowest in E3 (20.4 mg kg−1) and highest in E6 (30.8 mg kg−1). Genetic variability for grain Zn was highest in E6 (15.6–57.8 mg kg−1), followed by E2 (16.4–49.2 mg kg−1). For grain Fe concentration, highest and lowest genetic variability was observed in E5 (12.6–58.4 mg kg−1) and E3 (7.81–34.9 mg kg−1). Narrower genetic variability for grain Zn and Fe concentration was observed in E3, E4, E7 and E8 (drought and low N) than under optimum conditions (Table 2).

2.3. Genotypic Variance, G × E Interaction and Heritability

Highly significant differences for δ2g were evident for all the six traits, except for AD (E2) and ASI (E2 and E6) as shown in Table 3. Among the agronomic traits, plant height (PH) had the highest δ2g in most individual experimental sites, with the highest value in E2. A similar trend was observed across sites. GY also showed highly significant δ2g in both individual locations and across environments, with the highest value at E2 (Table 3). The results indicate that in all the sites across years, δ2g for GY was higher at optimum sites than at stress-induced sites. The δ2gy interactions were not significant for all the agronomic traits, whereas for δ2ge interactions, only AD and GY were significant. Comparing the two agronomic traits, δ2ge for AD was significant but with a small magnitude, whilst for GY, δ2ge was highly significant (Table 4). Although δ2gye interactions were significant for all the agronomic traits, its significance was less than for genotypic variance of these traits.
For micronutrient concentration, the δ2g for grain Zn and Fe concentration was highly significant (p ≤ 0.01) in all individual experimental sites (Table 3). Grain Zn concentration had the highest genotypic variance in E1 (CIMMYT 2018/19), and lowest in E7 (DR&SS 2019/20). For grain Fe concentration, the highest δ2g was observed in E5 (CIMMYT 2019/20), whereas E4 (Chisumbanje 2018/19) had the lowest genotypic variance. Similar to observations on GY, δ2g for grain micronutrient concentration was higher under optimum sites than under managed drought or low N conditions. While δ2g was significant across sites using pooled data, δ2gy interactions were not significant for both grain Zn and Fe concentrations (Table 4). However, both δ2ge and δ2gye interactions were highly significant (p ≤ 0.01). Apparently, δ2ge interactions for the micronutrients were higher than δ2g.
Most of the phenotypic traits were highly heritable (H2 > 60%) at individual experimental sites except for AD (E2) and ASI in several sites (Table 3). In stress and non-stress environments, PH, GY, grain Zn and Fe concentrations had the highest heritability. Grain yield was more heritable under optimum than stress sites. In contrast, heritability estimates in stress and non-stress environments were similar for grain micronutrients. Compared to agronomic traits, higher heritability was observed for grain micronutrients (Table 3). In the combined analysis, heritability ranged from 56–89%, with GY having the highest heritability, followed by PH (Table 4). Comparing the two micronutrients, grain Zn concentration was more heritable than grain Fe concentration across sites and years.

2.4. Grain Yield and Micronutrient Performance by Nutritional Type

The Analysis of Variance (ANOVA) showed significant differences (p ≤ 0.05) for GY performance among the four nutritional groups: Zn × NML, Zn × ProA, Zn × QPM and checks (ANOVA table not shown). The effects of management and management × nutritional group interactions were significant (p ≤ 0.05). For all the nutritional groups, GY performance was highest under optimum conditions and lowest under managed drought conditions. Commercial checks had the highest GY under optimum (8.40 t ha−1) and low N (2.40 t ha−1) conditions and across sites (Figure 2). However, the GY performance of the Zn × NML category was comparable to that of the checks. Despite the lower GY performance, Zn × ProA and Zn × QPM hybrids yielded at least 7.00 t ha−1 under optimum conditions. GY performance of these two nutritional groups was comparable for most management conditions. However, Zn × ProA hybrids had the highest GY performance under managed drought conditions, although there were no large differences among the four nutritional groups. ANOVA for grain Zn and Fe concentration showed highly significant differences (p ≤ 0.05) for the management and nutritional groups (data not shown). However, management × nutritional type interaction was non-significant for grain Zn concentration, whereas high significance was observed for Fe. For all the nutritional groups, higher grain Zn concentrations were observed under optimum than stress conditions (Figure 3). Zn × QPM hybrids accumulated the highest grain Zn concentration (34.0 mg kg−1) across all management levels (Figure 3), followed by Zn × ProA hybrids (31.3 mg kg−1). Commercial checks had the lowest grain Zn concentration (25.2 mg kg−1). The general observation across sites in terms of grain Zn accumulation was Zn × QPM > Zn × ProA > Zn × NML > check hybrids (Figure 4). The Zn × QPM and Zn × ProA hybrids accumulated higher grain Fe concentration than Zn × NML and checks (Figure 5). However, cross-over interactions were observed between nutritional groups and management levels for grain Fe concentration.

2.5. Correlations between Agronomic and Micronutrient Traits

Grain Zn and Fe concentration showed positive and significant correlations (Table 5) at all individual sites (E1 = 0.75, E2 = 0.72, E3 = 0.52, E4 = 0.67, E5 = 0.56, E6 = 0.70, E7 = 0.54 and E8 = 0.51; p ≤ 0.01). Across site correlation analysis, a similar trend was observed (r = 0.97; p ≤ 0.01). The association between GY and grain Zn concentration was negative and non-significant in all environments, except for two sites (E7 and E8), but in all cases the magnitude was small. However, across site analysis revealed that GY and grain Zn concentration were moderately and significantly correlated (r = −0.44; p ≤ 0.01). In contrast, no trend was observed for correlation between GY and grain Fe concentration. However, the correlation of these traits was moderate, negative and significant (r = −0.43; p ≤ 0.01) across sites (Table 5). The association between grain Zn concentration and either AD or ASI was not significant for all sites except for AD in E2. In contrast, across site analysis indicated positive and significant association between these traits. Highly positive and significant correlation was observed between PH and GY in individual sites and across sites (E1 = 0.58, E2 = 0.66, E3 = 0.57, E4 = 0.28, E5 = 0.74, E6 = 0.58, E7 = 0.61 and E8 = 0.83; Pooled = 0.99; p ≤ 0.01). AD was positively associated with GY at all optimum sites except for E6, while high and significant negative correlation was observed under stress conditions (E3, E4, E7 and E8). The association between ASI and GY was negative and stronger under stress than non-stress environments.

3. Discussion

The current study aimed to, first, determine the soil extractable Zn and Fe concentration in major maize experimental sites in Zimbabwe, second, evaluate the genetic variability of grain Zn and Fe concentration in newly developed Zn-enhanced normal, provitamin A and QPM testcrosses and link this with soil physicochemical properties and crop management; third, examine correlation between grain Zn and Fe concentration and agronomic traits, in an effort to develop high yielding Zn-enhanced cultivars with substantial grain Fe concentration.
Moderate but significant diversity at the experimental sites in terms of soil extractable Zn and Fe concentration, pH, total soil carbon, soil organic carbon and texture was observed. All these soil properties, including the root microbiome are critical in determining the extent of soil Zn and Fe solubility and subsequent accumulation in edible crop parts [14,22]. The soil pH for sites ranged from slightly acidic to slightly alkaline where maize normally grows optimally [23,24]. Having soil pH of such range implies that the pH had minimum inhibitory effects to grain Zn and Fe accumulation, considering that bioavailability of these micronutrients is greatly affected under extremes of both acidity and alkalinity [25]. Experimental sites differed in terms of soil extractable Zn and Fe concentration, with ART farm and CIMMYT having the highest soil extractable Zn and Fe concentration, respectively, while Chisumbanje had the lowest concentration for the two micronutrients. Although experimental sites differed in terms of soil extractable Zn and Fe, all the sites had sufficient concentrations of these micronutrients required to sustain maize growth and development based on critical limits of 1.00 and 5.00 mg kg−1 for soil Zn and Fe, respectively [26]. Previous studies have demonstrated that grain Zn and Fe concentration in maize sharply increases with an increase in soil extractable Zn and Fe concentration [7,27]. Whilst genotypic differences observed in the present study for grain Zn and Fe concentration between sites could have been attributed to differences in soil extractable Zn and Fe concentration, the impact was small. This is because the low N sites (E3 and E7) had comparably higher soil extractable micronutrients but lower means for grain Zn and Fe concentration.
Differences in means and extent of genetic variability for grain Zn and Fe concentration observed among the three crop management practices suggest high significance of crop management in determining grain micronutrient concentration. It is well documented that nutrient uptake decreases drastically under moisture stress conditions [28,29] and this could explain the low grain Zn and Fe concentration at all the managed drought sites. In addition, the combination of both optimal N fertilizer rates applied and higher inherent soil N content at CIMMYT 2018/19 (E1) and ART farm 2019/20 (E5) could be responsible for higher grain Zn and Fe in those sites as compared to low N sites. Previously, a 14.6–26.9% increase in grain Zn concentration in maize compared to trials with no N application was reported [30]. In the present study, optimum crop management increased the grain Zn concentration by 21.5% as compared to the average of all low N sites. Similar findings have also been reported by other authors [31,32]. Contrary to this, a different study [33] reported negative correlations between N fertilizer application and grain Zn concentration. As expected, higher N sites in the present study were associated with higher soil organic carbon, a major component of soil organic matter [34]. Direct and indirect influence of soil organic matter on soil micronutrient availability and uptake by plants was reported in different studies [35,36]. All these results demonstrate that the accumulation of micronutrients in maize grain is highly dependent of several environmental factors [14]. Although results from this study demonstrate significant differences in terms of soil physicochemical properties that might influence grain Zn and Fe concentration in maize, the experimental sites used were not fully representative of target environments, since most smallholder farmers in Zimbabwe grow maize in highly acidic sandy soils with very low organic matter [37]. Future studies should consider conducting such experiments in smallholder farms that are often characterized by low productivity potential. In addition to that, future studies should explore the correlation between extractable Zn and Fe concentration in soil solution and the micronutrient concentration accumulated in maize grains.
The mean performance of genotypes in both individual and across environments shows that sufficient genetic variability for both agronomic and micronutrient traits exist. This suggests that improvement of grain micronutrients and agronomic traits in locally adapted normal, provitamin A and QPM germplasm is possible. Genotypic variability in the gene pool available to the breeder is crucial, as this facilitates selection of genotypes with superior performance in traits in question [38]. Genotypic variances observed for most traits in the present study were highly significant coupled with high heritability, indicating great potential to improve such traits (Table 3). Plant height had the highest genotypic variance in almost all sites, indicating strong genetic control of this trait and minimum environmental influence [39]. The traits phenotyped in this study are quantitative and are influenced by several genes interacting with various environmental factors, and this justifies highly significant δ2ge for these traits. G × E interaction reduces the broad-sense heritability across sites and this ultimately reduces the selection accuracy [40]. Comparing the two micronutrients, grain Fe concentration had higher δ2ge than grain Zn concentration and this indicates stronger influence of the environment for Fe than for Zn. However, δ2gy was non-significant, implying that growing seasons were similar across the years. In contrast, δ2gye was highly significant and of greater magnitude compared to δ2g, for all traits except for grain Zn and Fe concentration. This shows that the environment highly influenced agronomic traits compared to micronutrient concentration. A similar low contribution of the environment to grain Fe and Zn concentration was reported [41,42] in pearl millets. Contrary to these findings, another study [19] reported a bigger contribution of the environment to the total micronutrient variation in sorghum. Comparing the two micronutrients, δ2gye for grain Fe was higher than for grain Zn concentrations and these results corroborate previous findings [19] in sorghum, but are contrary to reports in maize [43]. Greater proportion of genotype × environment × year effect for grain Fe suggests that this micronutrient is more sensitive to environmental changes than grain Zn concentration. This further supports the well-documented evidence that the accumulation of grain micronutrients is highly influenced by several factors including soil moisture, pH, soil type and fertility, genotype and interactions of among nutrients in soil [14,43]. The presence of large genotypic variance and moderate G × E interactions for grain Zn and Fe concentration suggests the possibility to increase the mineral concentration through both modern and conventional breeding approaches.
The optimum environments were more favourable to most Zn-enhanced genotypes compared to stress environments. Both drought and low N stress increases kernel abortion, delays silking, reduces the number of ears per plant, which all have negative impacts to grain yield [44]. Furthermore, reduced soil moisture could have contributed to the low micronutrient absorption and uptake. Several studies have demonstrated that grain Zn concentration in maize [7,30], and wheat [45] increases with N supply. Therefore, the N fertilizer limitation in DR&SS site might have decreased grain Zn concentration. Evidence on whether N supply has direct influence on grain Fe accumulation in maize is still limited, although few studies have reported significant improvement of grain Fe with N application in cowpea [46].
The nutritional backgrounds (types) of the evaluated hybrids were also significant in determining grain yield and micronutrient concentration. The highest grain yield observed in non-biofortified commercial checks, followed by Zn × NML hybrids concur with several studies reporting dilution effects [4,17]. However, the differences in grain yield performance between Zn-enhanced hybrids and checks were not large and this suggests the possibility of developing high yielding and nutritionally superior maize genotypes. In terms of micronutrient performance, the QPM-based hybrids had comparably better performance than other nutritional backgrounds. The superiority of QPM germplasm for grain Zn concentration has been reported by several authors [1,8,13,47]. However, high grain Zn in normal maize has also been reported [48] and this is quite encouraging to breeders. The positive and highly significant correlation between grain Zn and Fe concentration suggests that the genes governing accumulation of these micronutrients could be linked, implying that simultaneous improvement is possible [4]. Correlation between grain yield and either grain Zn or Fe concentration was negative and moderate, although it was highly significant. While other studies have reported similar findings [4], there have also been contrasting results [13,49]. Based on the results, it is apparent that the dilution effects were significant but only to a limited extent. This shows that genotypes with moderate to high yield potential coupled with high levels of grain micronutrients can be developed. Application of foliar Zn-containing fertilizers can also increase the grain Zn concentration in case of high yielding cultivars with moderate ability to accumulate grain Zn in soil [49].

4. Materials and Methods

4.1. Plant Materials

Seventy-seven Zn-enhanced normal, provitamin A and QPM testcross hybrids were developed by crossing 11 Zn donors derived from CIMMYT-Mexico and the International Institute of Tropical Agriculture (IITA) in West Africa and 7 testers of normal, provitamin A and QPM backgrounds. Testers were developed by CIMMYT in Southern Africa. Testcrossing was done during the 2017/18 cropping season at CIMMYT’s Harare Station (Latitude = 17°43′ S; Longitude = (31°05′ E); Altitude = 1480 masl; Average annual rainfall = 850 mm). Planting followed an 11 × 7 North Carolina II design, with Zn donors used as males while normal, provitamin A and QPM line testers were used as females. The testers were locally adapted tropical germplasm that has been widely used in evaluating newly developed inbred lines. Baseline Zn and Fe concentrations for the potential parental lines were measured (Table 6) and final selections were based on considerable micronutrient densities and tolerance to biotic and abiotic stresses such as drought and disease (maize streak virus and grey leaf spot). In addition, the provitamin A and QPM testers were previously analysed and selected for high levels of β-carotene and tryptophan content. All the generated 77 testcross hybrids had enough seed quantities to be evaluated in trials together with seven commercial checks.

4.2. Experimental Sites

The study was conducted in eight environments in Zimbabwe (Table 7), during the 2018/19 and 2019/20 cropping seasons. These sites are widely used for evaluating agronomic performance of breeding materials under development. The testcross hybrids were grown under two optimum, one low N and one managed drought conditions for two years. Both optimum and low N trials were planted in the main (summer) cropping season, while the managed drought trials were grown under irrigation during the winter, rain-free period. Maize plants were exposed to drought stress by withholding irrigation two weeks before flowering up to 21 days after flowering, so that drought stress would coincide with flowering, which is the most sensitive growth stage. Similarly, the low N trials were grown in low N screening sites that were developed by continuously depleting N to less than 7 ppm, which is estimated to cause 30% maize yield reduction [50].
In Zimbabwe, agro-ecological zones, also known as natural farming regions (NFRI–V), are defined based on the annual rainfall amount received in a particular area, and soil quality and type of vegetation, among other factors [51]. NFR1 receives the highest annual rainfall amount of >1000 mm and NFRV receives the least rainfall amount of <450 mm7. A representation of high and low potential environments with contrasting agro-ecologies was included in this study (Table 7).

4.3. Trial Layout, Management and Data Collection

The 84 genotypes were laid out in the field using an alpha (0.1) lattice design with two replications. Briefly, each replicate accommodated 12 incomplete blocks with a block size of 7 genotypes. Genotypes were planted in single row plots of 4 m long with 17 planting stations. Two seeds were sown in each planting hole and thinned two weeks after emergence to achieve a uniform crop stand. Spacing was 0.75 m (inter-row) and 0.25 m (in-row) to achieve a plant density of 53,000 ha–1. A total of four border rows were planted on each side of the main experimental plot.
Standard agronomic and cultural practices were applied at all sites during crop growth and development. Weed control was done using herbicides and in some cases hand weeding was applied. Supplementary irrigation was applied when necessary in all the trials during the early vegetative stages. Two random plants from each plot were selfed to produce seed for the micronutrient analysis and all other plants were allowed to open-pollinate for phenotypic traits data collection, primarily for grain yield (GY), days to mid-anthesis (AD), anthesis-silking-interval (ASI) and plant height (PH). GY recorded using plants in the net plot area as the two border plants close to the alley were discarded. AD was recorded as number of days from planting to when 50% of the plants had tassels shedding pollen. Plant height (PH) was recorded using a laser distance meter as the average of all plants in the plot. Lastly, ASI was determined as the difference between days to silking and anthesis.

4.4. Soil Micronutrient Analysis

Soil analysis was done to determine the soil micronutrient concentration and other physicochemical properties for each experimental site used. Soil samples were taken using a 20 cm hand operated auger in an M-shaped pattern after scraping away surface litter [52]. The edges of the field were excluded from sampling to avoid contamination from residual fertilizers, since bags are normally placed at the periphery of the field during application. A total of 15 subsamples were collected at each experimental site to make a composite soil sample for analysis. Soil samples were sent to Rothamsted Research Soil Science Laboratory, United Kingdom. The primary micronutrients, including Zn, Fe, magnesium (Mg), and copper (Cu) were simultaneously extracted using the Mehlich 3 method [53]. Extractable micronutrients were determined by an inductively coupled plasma-optical emission spectrophotometer (ICP-OES; model: ICAP 6000 series, Thermo Fisher Scientific, Dreieich, Germany). The conventional Olsen P method [54] was used for determining available P. Extraction of exchangeable metal cations such as calcium (Ca), potassium (K), Mg, and sodium (Na) was done using ammonium chloride [55]. Total nitrogen (TN) and total carbon (TC) were determined using the Dumas combustion procedure using the Leco CNS 2000 analyzer [56]. Inorganic carbon (IC) was determined using finely ground (0.15 mm) soils in the Primacs SNC analyser (Skalar Inc., Buford, Georgia). Soil pH was measured in a 1:2.5 soil/water mixture using an electrode pH meter.

4.5. Grain Micronutrient Analysis

Ears from the selfed maize plants were manually harvested and shelled to avoid any metal contamination from mechanical harvesting. Grain of all selfed plants in a plot was bulked to form a 100 g composite sample. Fifty kernels per sample were ground using a Foss mini flour industrial maize grinder. Milling time was 60 s at 15 Hz. The flour was sieved through a 1 mm stainless steel sieve and stored at 4 °C in screw-top polycarbonate vials. Two grams of flour was dry-ashed in porcelain crucibles, placed in a muffle furnace at 550 °C for 3 h and digested using concentrated HNO3 acid [57]. Grain micronutrient concentration (Zn and Fe) was determined using the atomic absorption spectrometry (AAS) [58].

4.6. Statistical Analysis

Analysis of variance (ANOVA) for the phenotypic traits of testcross hybrids across environments was carried out using the Best Linear Unbiased Prediction (BLUP) approach in Multi-Environmental Trial Analysis (META-R) v6.0 [59]. The variance components (δ2g, δ2gy, δ2ge, δ2gyl and δ2ε) were estimated in the linear mixed model, where the environment and year were considered as fixed factors, while genotypes, replicates, blocks and genotype interactions with year and environment were considered as random effects. Grain yield was adjusted to 12.5% moisture content before analysis. Pooling data for combined analysis was done after conducting a Bartlett’s test for testing for homogeneity of variances between environments and years [60]. GenStat software 17th edition [61] was also used for generating descriptive statistics.
Genotypic correlations among micronutrient and agronomic traits were performed using BLUPs of individual and across sites in META-R. The combined analysis, taking into account the δ2gye effect was performed using the following model:
Yjklmn = µ + Yj + Ek + YEjk + rjkl + Bjklm + Gn + (YG)jn + (EG)kn + (YEG)jkn + εjklmn
where µ is the grand mean; Yj is the fixed effect of year j; Ek is the fixed effect of environment k; YEjk is the fixed effect of interaction between year j and environment k; rjkl is the random effect of replication l in environment k and year j; Bjklm is the random effect of block m nested in replication l in environment k and year j; Gn is the random effect of genotype n; (YG)jn is the random effect of the interaction between genotype n and year j; (EG)kn is the random effect of the interaction between genotype n and environment k; (YEG)jkn is the random effect of the interaction effect of the genotype n in environment k and year j; εjklmn is the residual variance.
Broad sense heritability (H2) for traits was estimated in META-R using the following formula:
H 2 = δ g 2 δ g   2 + δ ε / r   2  
Where δ2g is the genotypic variance, r and δ2ε denotes the number of replicates and residual variance, respectively.

5. Conclusions

The current study showed significant differences (genotypic and between nutritional type) in terms of both micronutrient and agronomic traits under optimum, low N and managed drought conditions. Experimental sites with high soil N, total carbon and organic matter had the highest grain Zn and Fe concentration and therefore agronomic biofortification could be useful in increasing these micronutrients in maize grain. Further studies using large populations are required for validation of the genotypic differences in micronutrient densities between normal, provitamin A and QPM germplasm. Considering limited differences in soil physicochemical properties in experimental sites observed in the present study, further studies should explore differences in smallholder farms, which are normally characterized by low clay to sand ratio. High and significant correlation between grain Zn and Fe concentration indicate the possibility of simultaneous improvement of these traits, which is quite cost-effective. Despite showing lower grain yields, micronutrient content of Zn-enhanced hybrids was higher than the control. Although, grain yield was negatively correlated with both grain micronutrients, the magnitude of the correlation coefficient was moderate and this gives hope in developing high yielding and micronutrient-dense maize cultivars, which could reduce the highly prevalent micronutrient deficiency in SSA.

Author Contributions

Conceptualization, N.G.-M. and M.L.; investigation, N.G.-M.; resources, T.N. and M.L.; Writing—Original draft preparation, N.G.-M.; Writing—Review and editing, N.G.-M., T.N., J.E.C., C.N.K., A.v.B. and M.L.; data analysis; N.G.-M., C.N.K., F.B. and J.E.C., supervision, M.L., T.N. and A.v.B.; funding acquisition, T.N. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bill & Melinda Gates Foundation and USAID funded project “Accelerating Genetic Gains” (grant number INV-003439, OPP1134248), the National Research Foundation of South Africa through the Research Chairs Initiative (UID 84647) and UK Global Challenges Research Fund administered by the Biotechnology and Biological Sciences Research Council for the project “Addressing malnutrition with biofortified maize in Zimbabwe: from crop management to policy and consumers” (IATI Identifier: GB-GOV-13-FUND–GCRF-BB_T009047_1).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data presented in this study are available on request from the corresponding author.

Acknowledgments

We are grateful to the Department of Research and Specialist Services in Zimbabwe for providing some of the experimental sites used in this study. We thank Mr Chikwari at the Chemistry and Soils Research Institute for his kind support in milling grain samples. We thank Cosmos Magorokosho for guidance and advice. Thanks to Ridokunda Sivhada, Toverengwa Chitana and Passmore Kasere for their assistance in the execution of seed preparation, field and laboratory activities.

Conflicts of Interest

The authors declare no known competing interest.

References

  1. Hindu, V.; Palacios-Rojas, N.; Babu, R.; Suwarno, W.B.; Rashid, Z.; Usha, R.; Saykhedkar, G.R.; Nair, S.K. Identification and validation of genomic regions influencing kernel zinc and iron in maize. Theor. Appl. Genet. 2018, 131, 1443–1457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bouis, H.E.; Saltzman, A. Improving nutrition through biofortification: A review of evidence from HarvestPlus, 2003 through 2016. Glob. Food Secur. 2017, 12, 49–58. [Google Scholar] [CrossRef] [PubMed]
  3. Botoman, L.; Nalivata, P.C.; Chimungu, J.G.; Munthali, M.W.; Bailey, E.H.; Ander, E.L.; Lark, R.M.; Mossa, A.W.; Young, S.D.; Broadley, M.R. Increasing zinc concentration in maize grown under contrasting soil types in Malawi through agronomic biofortification: Trial protocol for a field experiment to detect small effect sizes. Plant Direct 2020, 4, e00277. [Google Scholar] [CrossRef] [PubMed]
  4. Menkir, A. Genetic variation for grain mineral content in tropical-adapted maize inbred lines. Food Chem. 2008, 110, 454–464. [Google Scholar] [CrossRef] [PubMed]
  5. Garcia-Oliveira, A.L.; Chander, S.; Ortiz, R.; Menkir, A.; Gedil, M. Genetic basis and breeding perspectives of grain iron and zinc enrichment in cereals. Front. Plant Sci. 2018, 9, 937. [Google Scholar] [CrossRef]
  6. Wessells, K.R.; Brown, K.H. Estimating the global prevalence of zinc deficiency: Results based on zinc availability in national food supplies and the prevalence of stunting. PLoS ONE 2012, 7, e50568. [Google Scholar] [CrossRef] [Green Version]
  7. Manzeke, G.M.; Mtambanengwe, F.; Watts, M.J.; Hamilton, E.M.; Lark, R.M.; Broadley, M.R.; Mapfumo, P. Fertilizer management and soil type influence grain zinc and iron concentration under contrasting smallholder cropping systems in Zimbabwe. Sci. Rep. 2019, 9, 6445. [Google Scholar] [CrossRef] [Green Version]
  8. Welch, R.M.; Graham, R.D. Breeding for micronutrients in staple food crops from a human nutrition perspective. J. Exp. Bot. 2004, 55, 353–364. [Google Scholar] [CrossRef] [Green Version]
  9. Hotz, C.; Brown, K.H. Assessment of the risk of zinc deficiency in populations and options for its control. Food Nutr. Bull. 2004, 25, S94–S203. [Google Scholar]
  10. Krebs, N.F.; Miller, L.V.; Hambridge, K.M. Zinc deficiency in infants and children: A review of its complex and synergistic interactions. Paediatr. Int. Child Health 2014, 34, 279–288. [Google Scholar] [CrossRef]
  11. Cakmak, I.; Kutman, U.B. Agronomic biofortification of cereals with zinc: A review. Eur. J. Soil Sci. 2018, 69, 172–180. [Google Scholar] [CrossRef]
  12. Nuss, E.T.; Tanumihardjo, S.A. Quality protein maize for Africa: Closing the protein inadequacy gap in vulnerable populations. Adv. Nutr. 2011, 2, 217–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mageto, E.K.; Lee, M.; Dhliwayo, T.; Palacios-Rojas, N.; Vicente, F.S.; Burgueño, J.; Hallauer, A.R. An evaluation of kernel zinc in hybrids of elite Quality Protein Maize (QPM) and non-QPM inbred lines adapted to the tropics based on a mating design. Agronomy 2020, 10, 695. [Google Scholar] [CrossRef]
  14. Cakmak, I. Enrichment of cereal grains with zinc: Agronomic or genetic biofortification? Plant Soil 2008, 302, 1–17. [Google Scholar] [CrossRef]
  15. Xia, H.; Kong, W.; Wang, L.; Xue, Y.; Liu, W.; Zhang, C.; Yang, S.; Li, C. Foliar Zn spraying simultaneously improved concentrations and bioavailability of Zn and Fe in maize grains irrespective of foliar sucrose supply. Agronomy 2019, 9, 386. [Google Scholar] [CrossRef] [Green Version]
  16. Vivek, B.S.; Krivanek, A.F.; Palacios-Rojas, N.; Twumasi-Afriyie, S.; Diallo, A.O. Breeding Quality Protein Maize (QPM): Protocols for Developing QPM Cultivars; CIMMYT: El Batán, Mexico, 2008. [Google Scholar]
  17. Maqbool, M.A.; Beshir, A.R. Zinc biofortification of maize (Zea mays L.): Status and challenges. Plant Breed. 2019, 138, 1–28. [Google Scholar] [CrossRef] [Green Version]
  18. Chakraborti, M.; Prasanna, B.M.; Hossain, F.; Mazumdar, S.; Singh, A.M.; Guleria, S.; Gupta, H.S. Identification of kernel iron- and zinc-rich maize inbreds and analysis of genetic diversity using microsatellite markers. J. Plant Biochem. Biotechnol. 2011, 20, 224–233. [Google Scholar] [CrossRef]
  19. Phuke, R.M.; Anuradha, K.; Radhika, K.; Jabeen, F.; Anuradha, G.; Ramesh, T.; Hariprasanna, K.; Mehtre, S.P.; Deshpande, S.P.; Anil, G.; et al. Genetic variability, genotype × environment interaction, correlation, and GGE biplot analysis for grain iron and zinc concentration and other agronomic traits in RIL population of sorghum (Sorghum bicolor L. Moench). Front. Plant Sci. 2017, 8, 712. [Google Scholar] [CrossRef] [Green Version]
  20. Bukvic, G.; Antunovic, M.; Popovic, S.; Rastija, M. Effect of P and Zn fertilisation on biomass yield and its uptake by maize lines (Zea mays L.). Plant Soil Environ. 2003, 49, 505–510. [Google Scholar] [CrossRef] [Green Version]
  21. Cheah, Z.X.; Kopittke, P.M.; Scheckel, K.G.; Noerpel, M.R.; Bell, M.J. Comparison of Zn accumulation and speciation in kernels of sweetcorn and maize differing in maturity. Ann. Bot. 2020, 125, 185–193. [Google Scholar] [CrossRef]
  22. Rutkowska, B.; Szulc, W.; Bomze, K.; Gozdowski, D.; Spychaj-Fabisiak, E. Soil factors affecting solubility and mobility of zinc in contaminated soils. Int. J. Env. Sci. Technol. 2015, 12, 1687–1694. [Google Scholar] [CrossRef]
  23. Opala, P.A. Influence of lime and phosphorus application rates on growth of maize in an acid soil. Adv. Agric. 2017, 2017, 7083206. [Google Scholar] [CrossRef] [Green Version]
  24. Tandzi, L.N.; Mutengwa, C.S.; Ngonkeu, E.L.M.; Woïn, N.; Gracen, V. Breeding for quality protein maize (QPM) varieties: A Review. Agronomy 2017, 7, 80. [Google Scholar] [CrossRef] [Green Version]
  25. Sadeghzadeh, B. A review of zinc nutrition and plant breeding. J. Soil Sci. Plant Nutr. 2013, 13, 905–927. [Google Scholar] [CrossRef] [Green Version]
  26. Lindsay, W.L.; Norvell, W.A. Development of a DTPA soil test for zinc, iron, manganese and copper. Soil Sci. Soc. Am. J. 1978, 42, 421–448. [Google Scholar] [CrossRef]
  27. Aref, F. Manganese, iron and copper contents in leaves of maize plants (Zea mays L.) grown with different boron and zinc micronutrients. Afr. J. Biotechnol. 2012, 11, 896–903. [Google Scholar] [CrossRef]
  28. Hussain, H.A.; Men, S.; Hussain, S.; Chen, Y.; Ali, S.; Zhang, S.; Zhang, K.; Li, Y.; Xu, Q.; Liao, C.; et al. Interactive effects of drought and heat stresses on morphophysiological attributes, yield, nutrient uptake and oxidative status in maize hybrids. Sci. Rep. 2019, 9, 3890. [Google Scholar] [CrossRef] [Green Version]
  29. Abbas, M.; Abdel-Lattif, H.; Shahba, M. Ameliorative effects of calcium sprays on yield and grain nutritional composition of maize (Zea mays L.) cultivars under drought stress. Agriculture 2021, 11, 285. [Google Scholar] [CrossRef]
  30. Xue, Y.F.; Yue, S.C.; Liu, D.Y.; Zhang, W.; Chen, X.P.; Zou, C.Q. Dynamic zinc accumulation and contributions of pre-and/or post-silking zinc uptake to grain zinc of maize as affected by nitrogen supply. Front. Plant Sci. 2019, 10, 1203. [Google Scholar] [CrossRef] [Green Version]
  31. Manzeke, G.M.; Mtambanengwe, F.; Watts, M.J.; Broadley, M.R.; Lark, R.M.; Mapfumo, P. Nitrogen effect on zinc biofortification of maize and cowpea in Zimbabwean smallholder farms. Agron. J. 2020, 112, 2256–2274. [Google Scholar] [CrossRef]
  32. Grujcic, D.; Yazici, A.M.; Tutus, Y.; Cakmak, I.; Singh, B.R. Biofortification of silage maize with zinc, iron and selenium as affected by nitrogen fertilization. Plants 2021, 10, 391. [Google Scholar] [CrossRef]
  33. Bojtor, C.; Illès, A.; Mousavi, S.M.; Szèles, A.; Tóth, B.; Nagy, J.; Marton, C.L. Evaluation of the nutrient composition of maize in different NPK fertilizer levels based on multivariate method analysis. Int. J. Agron. 2021, 2021, 5537549. [Google Scholar] [CrossRef]
  34. Nanganoa, L.T.; Ngome, F.A.; Suh, C.; Basga, S.D. Assessing soil nutrients variability and adequacy for the cultivation of maize, cassava, and sorghum in selected agro-ecological zones of Cameroon. Int. J. Agron. 2020, 2020, 8887318. [Google Scholar] [CrossRef]
  35. Rengel, Z.; Batten, G.G.; Crowley, D.E. Agronomic approaches for improving the micronutrient density in edible portions of field crops. Field Crops Res. 1999, 60, 27–40. [Google Scholar] [CrossRef]
  36. Dhaliwal, S.S.; Naresh, R.K.; Mandal, A.; Singh, R.; Dhaliwal, M.K. Dynamics and transformations of micronutrients in agricultural soils as influenced by organic matter build-up: A review. Environ. Sustain. Indic. 2019, 1–2, 100007. [Google Scholar] [CrossRef]
  37. Gwenzi, W.; Munondo, R. Long-term impacts of pasture irrigation with treated sewage effluent on nutrient status of a sandy soil in Zimbabwe. Nutr. Cycl. Agroecosys. 2019, 82, 197–207. [Google Scholar] [CrossRef]
  38. Ortiz-Monasterio, J.I.; Palacios-Rojas, N.; Meng, E. Enhancing the mineral and vitamin content of wheat and maize through plant breeding. J. Cereal Sci. 2007, 46, 293–307. [Google Scholar] [CrossRef]
  39. Tucker, S.L.; Dohleman, F.G.; Grapov, D.; Flagel, L.; Yang, S.; Wegener, K.M.; Kosola, K.; Swarup, S.; Rapp, R.A.; Bedair, M.; et al. Evaluating maize phenotypic variance, heritability, and yield relationships at multiple biological scales across agronomically relevant environments. Plant Cell Environ. 2019, 43, 880–902. [Google Scholar] [CrossRef] [Green Version]
  40. Liu, X.; Wang, H.; Wang, H.; Guo, Z.; Xu, X.; Liu, J.; Wang, S.; Li, W.; Zou, C.; Prasanna, B.M.; et al. Factors affecting genomic selection revealed by empirical evidence in maize. Crop J. 2018, 6, 341–352. [Google Scholar] [CrossRef]
  41. Govindaraj, M.; Rai, K.N.; Kanatti, A.; Upadhyaya, H.D.; Shivade, H.; Rao, A.S. Exploring the genetic variability and diversity of pearl millet core collection germplasm for grain nutritional traits improvement. Sci. Rep. 2020, 10, 21177. [Google Scholar] [CrossRef]
  42. Pujar, M.; Govindaraj, M.; Gangaprasad, S.; Kanatti, A.; Shivade, H. Genetic variation and diversity for grain iron, zinc, protein and agronomic traits in advanced breeding lines of pearl millet [Pennisetum glaucum (L.) R. Br.] for biofortification breeding. Genet. Resour. Crop Evol. 2020, 67, 2009–2022. [Google Scholar] [CrossRef]
  43. Prasanna, B.M.; Mazumdar, S.; Chakraborti, M.; Hossain, F.; Manjaiah, K.M.; Agrawal, P.K. Genetic variability and genotype × environment interactions for kernel Fe and Zn concentrations in maize (Zea mays L.) genotypes. Indian J. Agric. Sci. 2011, 81, 704–711. [Google Scholar]
  44. Meseka, S.; Menkir, A.; Bossey, B.; Mengesha, W. Performance assessment of drought tolerant maize hybrids under combined drought and heat stress. Agronomy 2018, 8, 274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kutman, U.B.; Yildiz, B.; Ozturk, L.; Cakmak, I. Biofortification of durum wheat with zinc through soil and foliar applications of nitrogen. Cereal Chem. 2010, 87, 1–9. [Google Scholar] [CrossRef]
  46. Manzeke-Kangara, M.G.; Mtambanengwe, F.; Watts, M.J.; Broadley, M.R.; Lark, R.M.; Mapfumo, P. Can nitrogen fertilizer management improve grain iron concentration of agro-biofortified crops in Zimbabwe? Agronomy 2021, 11, 124. [Google Scholar] [CrossRef]
  47. Maqbool, M.A.; Issa, A.B.; Khokhar, E.S. Quality protein maize (QPM): Importance, genetics, timeline of different events, breeding strategies and varietal adoption. Plant Breed. 2021, 140, 375–399. [Google Scholar] [CrossRef]
  48. Queiroz, V.A.V.; Guimarães, P.E.O.; Queiroz, L.R.; Guedes, E.O.; Vasconcelos, V.D.B.; Guimarães, L.J.; Ribeiro, P.E.A.; Schaffert, R.E. Iron and zinc availability in maize lines. Ciência Tecn. Alim. 2011, 31, 577–583. [Google Scholar] [CrossRef] [Green Version]
  49. Liu, D.Y.; Zhang, W.; Liu, Y.M.; Chen, X.P.; Zou, C.Q. Soil application of zinc fertilizer increases maize yield by enhancing the kernel number and kernel weight of inferior grains. Front. Plant Sci. 2020, 11, 188. [Google Scholar] [CrossRef]
  50. Bänziger, M.; Edmeades, G.O.; Beck, D.; Bellon, M. Breeding for Drought and Nitrogen Stress Tolerance in Maize: From Theory to Practice; CIMMYT: El-Batán, Mexico, 2000. [Google Scholar]
  51. Mugandani, M.; Wuta, A.; Makarau, A.; Chipindu, B. Re-classification of agro-ecological regions of Zimbabwe in conformity with climate variability and change. Afr. Crop Sci. J. 2012, 20, 361–369. [Google Scholar]
  52. Motsara, M.R.; Roy, R.N. Guide to Laboratory Establishment for Plant Nutrient Analysis; FAO Fertilizer and Plant Nutrition Bulletin: Rome, Italy, 2008; Volume 19. [Google Scholar]
  53. Mehlich, A. Mehlich 3 soil test extractant: A modification of Mehlich 2 extractant. Commun. Soil Sci. Plant Anal. 1984, 15, 1409–1416. [Google Scholar] [CrossRef]
  54. Olsen, S.R.; Cole, C.V.; Watanabe, F.S.; Dean, L.A. Estimation of available phosphorus in water and sodium bicarbonate extracts from soils. Soil Sci. Am. Proc. 1954, 29, 677–678. [Google Scholar]
  55. Tucker, B.M. The determination of exchangeable calcium and magnesium in carbonate soils. Austr. J. Agric. Res. 1954, 5, 705–715. [Google Scholar] [CrossRef]
  56. Etheridge, R.D.; Pesti, G.M.; Foster, E.H. A comparison of nitrogen values obtained utilizing the Kjeldahl nitrogen and Dumas combustion methodologies (Leco CNS 2000) on samples typical of an animal nutrition analytical laboratory. Anim. Feed. Sci. Technol. 1998, 73, 21–28. [Google Scholar] [CrossRef]
  57. Galicia, L.; Nurit, E.; Rosales, A.; Palacios-Rojas, N. Laboratory Protocols. Maize Nutrition Quality and Plant Tissue Analysis Laboratory; CIMMYT: El Batán, Mexico, 2009. [Google Scholar]
  58. Zarcinas, B.A.; Cartwright, B.; Spouncer, L.R. Nitric acid digestion and multi element analysis of plant material by inductively coupled plasma spectrometry. Commun. Soil Sci. Plant Anal. 1987, 18, 131–146. [Google Scholar] [CrossRef]
  59. Alvarado, G.; López, M.; Vargas, M.; Pacheco, A.; Rodríguez, F.; Burgueño, J.; Crossa, J. Multi-Environment Trial Analysis with R (META-R) Version 6.04; CIMMYT: El Batán, Mexico, 2015. [Google Scholar]
  60. Gomez, K.A.; Gomez, A.A. Statistical Procedures for Agricultural Research; John Wiley and Sons: New York, NY, USA, 1984. [Google Scholar]
  61. VSN International. Genstat for Windows, 17th ed.; VSN International: Hemel Hempstead, UK, 2020. [Google Scholar]
Figure 1. Density plot for grain yield (t ha−1) performance of test hybrids under optimum, low nitrogen and managed drought conditions. The area under the curve represents the probability of getting hybrids with certain grain yield performance along the x-axis.
Figure 1. Density plot for grain yield (t ha−1) performance of test hybrids under optimum, low nitrogen and managed drought conditions. The area under the curve represents the probability of getting hybrids with certain grain yield performance along the x-axis.
Plants 12 00270 g001
Figure 2. Grain yield performance (t ha−1) for the Zn-enhanced normal (Zn × NML), provitamin A (Zn × ProA) and QPM (Zn × QPM) testcross hybrids and checks evaluated under stress and non-stress environments.
Figure 2. Grain yield performance (t ha−1) for the Zn-enhanced normal (Zn × NML), provitamin A (Zn × ProA) and QPM (Zn × QPM) testcross hybrids and checks evaluated under stress and non-stress environments.
Plants 12 00270 g002
Figure 3. Grain zinc concentration (mg kg−1) for the Zn-enhanced normal (Zn × NML), provitamin A (Zn × ProA) and QPM (Zn × QPM) testcross hybrids and checks evaluated under stress and non-stress environments.
Figure 3. Grain zinc concentration (mg kg−1) for the Zn-enhanced normal (Zn × NML), provitamin A (Zn × ProA) and QPM (Zn × QPM) testcross hybrids and checks evaluated under stress and non-stress environments.
Plants 12 00270 g003
Figure 4. Nutritional profile showing Zn concentration versus Zn yield of testcross hybrids and checks grown under low N and managed drought conditions where optimum management was the reference.
Figure 4. Nutritional profile showing Zn concentration versus Zn yield of testcross hybrids and checks grown under low N and managed drought conditions where optimum management was the reference.
Plants 12 00270 g004
Figure 5. Grain iron concentration (mg kg−1) for the Zn-enhanced normal (Zn × NML), provitamin A (Zn × ProA) and QPM (Zn × QPM) testcross hybrids and checks evaluated under stress and non-stress environments.
Figure 5. Grain iron concentration (mg kg−1) for the Zn-enhanced normal (Zn × NML), provitamin A (Zn × ProA) and QPM (Zn × QPM) testcross hybrids and checks evaluated under stress and non-stress environments.
Plants 12 00270 g005
Table 1. Means for extractable soil physicochemical properties of different experimental sites.
Table 1. Means for extractable soil physicochemical properties of different experimental sites.
2019
SiteManagementpHTC (%)IC (%)TN (%)SOC (%)Olsen P mg kg−1CEC
meq−100 g
Zn mg kg−1Fe mg
kg−1
% Clay%
Silt
% Sand
ART farmOptimum5.801.400.020.151.3025.914.25.1057.839.549.511.0
CIMMYTOptimum6.402.200.010.202.2010.424.52.7076.140.540.519.0
DR&SSLow N6.100.800.020.040.7017.816.44.7056.152.026.521.5
ChisumbanjeManaged drought7.601.700.030.151.7015.728.11.9041.173.026.01.00
F-test *** ******************
2020
ART farmOptimum5.701.400.020.151.4023.413.25.7052.337.548.514.0
CIMMYTOptimum6.302.300.010.182.3012.125.82.7077.941.539.519.0
DR&SSLow N6.200.600.020.050.6018.514.94.8057.651.526.022.5
ChisumbanjeManaged drought7.601.800.020.171.7016.029.12.0039.172.525.02.50
F-test **** ******************
F-test for significant differences between experimental sites. * p < 0.05, ** p < 0.01.
Table 2. Descriptive statistics of the performance of Zn-enhanced testcross hybrids grown under optimum, low N and managed drought conditions.
Table 2. Descriptive statistics of the performance of Zn-enhanced testcross hybrids grown under optimum, low N and managed drought conditions.
EnvironmentsManagementStatisticsTraits Ұ
ADASIPHGYZnFe
CIMMYT 2018/19 (E1)OptimumMean69.21.402186.8030.729.7
Range60.0–74.0−1.00–3.00163–2532.90–13.712.9–46.711.3–55.1
SD2.160.7420.82.106.398.14
ART farm 2018/19 (E2)OptimumMean68.61.402457.5029.930.1
Range58.0–74.0−3.00–5.00159–2881.90–18.816.4–49.212.2–54.9
SD3.061.0922.33.535.807.91
DR&SS 2018/19 (E3)Low NMean71.03.001861.8024.820.4
Range63.0–76.00.00–6.00148–2260.40–4.1013.9–39.57.81–34.9
SD2.811.6517.40.865.246.08
CHISUMBANJE 2018/19 (E4)Managed DroughtMean68.02.701741.8026.223.0
Range64.0–73.00.00–7.00140–2050.40–4.1010.7–39.47.10–37.6
SD1.881.5014.10.825.094.85
CIMMYT 2019/20 (E5)OptimumMean69.61.402078.5030.529.9
Range62.0–76.0−5.00–4.00150–2454.30–14.317.3–45.312.6–58.4
SD3.140.9022.32.075.878.06
ART farm 2019/20 (E6)OptimumMean72.62.002069.3030.430.8
Range65.0–84.0−5.00–6.00150–2406.50–13.015.6–57.812.5–54.4
SD3.751.2818.81.306.327.69
DR&SS 2019/20 (E7)Low NMean71.12.901822.0025.223.1
Range64.0–76.00.00–7.00145–2200.40–5.6012.4–36.69.20–37.1
SD2.451.5616.40.894.926.12
CHISUMBANJE 2019/20 (E8)Managed DroughtMean69.42.8.01611.1025.121.1
Range64.0–77.00.00–6.00116–2140.10–3.8012.1–36.87.90–42.8
SD2.711.5118.10.725.236.10
Ұ AD, Days to 50% anthesis; ASI, Anthesis silking interval (days); PH, Plant height (cm); GY, Grain yield (t ha−1); Zn, grain zinc concentration (mg kg−1); Fe, Grain iron concentration (mg kg−1), SD, Standard deviation.
Table 3. Genotypic variance (δ2g), standard error (SE) and broad-sense heritability (H2) for agronomic and micronutrient traits at eight experimental sites.
Table 3. Genotypic variance (δ2g), standard error (SE) and broad-sense heritability (H2) for agronomic and micronutrient traits at eight experimental sites.
Trait ҰCIMMYT 2018/19 (E1)ART Farm 2018/19 (E2)DR&SS 2018/19 (E3)Chisumbanje 2018/19 (E4)
δ2gSEH2δ2gSEH2δ2gSEH2δ2gSEH2
AD3.05 **0.170.791.790.240.556.75 **0.220.921.89 **0.150.72
ASI0.19 **0.060.510.160.080.532.07 **0.130.861.46 **0.120.78
PH20.4 **1.610.68305 **1.720.79202 *1.350.8183.2 **1.090.65
GY3.14 **0.160.859.57 **0.270.870.47 **0.070.780.38 **0.060.73
Zn38.9 **0.490.9727.9 **0.450.8924.5 **0.400.9422.8 **0.390.94
Fe49.8 **0.630.8758.3 **0.610.9734.6 **0.470.9620.6 **0.370.94
TraitCIMMYT 2019/20 (E5)ART Farm 2019/20 (E6)DR&SS 2019/20 (E7)Chisumbanje 2019/20 (E8)
δ2gSEH2δ2gSEH2δ2gSEH2δ2gSEH2
AD7.72 **0.240.9011.5 *0.290.923.83 *0.190.785.98 **0.210.89
ASI0.24 *0.250.560.15 ns0.100.511.41 *0.120.761.27 **0.120.73
PH274 **1.720.73217 **1.450.76142 **1.270.6986.4 *1.390.63
GY3.59 **0.160.921.19 **0.100.830.47 **0.070.750.29 **0.060.73
Zn32.0 **0.450.9634.7 **0.490.9321.3 **0.380.9323.4 **0.400.92
Fe59.8 **0.620.9652.9 **0.590.9435.4 **0.470.9734.1 **0.470.95
Ұ AD, Days to 50% anthesis; ASI, Anthesis silking interval (days); PH, Plant height (cm); GY, Grain yield (t ha−1); Zn, grain zinc concentration (mg kg−1); Fe, Grain iron concentration (mg kg−1). * p < 0.05, ** p < 0.01; ns, not significant at p < 0.05.
Table 4. Genotypic variance (δ2g), Genotype × Year (δ2gy), Genotype × Environment (δ2ge), Genotype × Year × Environment (δ2gye) interaction, standard error (SE) and broad-sense (H2) heritability for traits of Zn-enhanced hybrids.
Table 4. Genotypic variance (δ2g), Genotype × Year (δ2gy), Genotype × Environment (δ2ge), Genotype × Year × Environment (δ2gye) interaction, standard error (SE) and broad-sense (H2) heritability for traits of Zn-enhanced hybrids.
Trait ҰCombined Data
δ2gSEδ2gySEδ2geSEδ2gyeSEH2
AD0.94 **0.290.120.251.63 *0.412.59 *0.210.59
ASI0.08 *0.040.090.040.010.110.87 *0.150.56
PH64.0 **13.811.07.741.320.0174.7 **0.420.69
GY0.30 **0.110.010.101.23 **0.170.84 *0.170.89
Zn11.2 **2.011.280.6614.0 **1.313.99 **0.320.81
Fe9.36 **2.240.551.4929.3 **2.405.35 **0.050.66
Ұ AD, Days to 50% anthesis; ASI, Anthesis silking interval (days); PH, Plant height (cm); GY, Grain yield (t ha−1); Zn, grain zinc concentration (mg kg−1); Fe, grain iron concentration (mg kg−1). * p < 0.05, ** p < 0.01.
Table 5. Correlation of grain Zn and Fe concentrations with agronomic traits in Zn-enhanced hybrids.
Table 5. Correlation of grain Zn and Fe concentrations with agronomic traits in Zn-enhanced hybrids.
Traits ҰCIMMYT 2018/19 (E1)Art Farm 2018/19 (E2)DR&SS 2018/19 (E3)Chisumbanje 2018/19 (E4)
GYZnFeGYZnFeGYZnFeGYZnFe
AD0.42 **0.070.220.18 **0.31 **0.25−0.72 **0.11−0.01−0.48 **0.20−0.03
ASI0.50 **−0.020.29 *−0.35 **0.140.38 **−0.88 **0.08−0.00−0.43 **0.040.15
PH0.58 **−0.14−0.160.66 **−0.04−0.190.57 **−0.11−0.22 *0.28 *−0.30 *−0.11
GY−0.010.04−0.12−0.18−0.000.06−0.18−0.12
Fe0.75 **0.72 **0.52 **0.67 **
TraitsCIMMYT 2019/20 (E5)Art farm 2019/20 (E6)DR&SS 2019/20 (E7)Chisumbanje 2019/20 (E8)
GYZnFeGYZnFeGYZnFeGYZnFe
AD0.41 **0.05−0.09−0.210.14−0.02−0.52 **0.110.18−0.32 **0.160.23 *
ASI0.16−0.08−0.09−0.070.060.24−0.64 **0.080.19−0.67 **−0.16−0.07
PH0.74 **−0.25 *−0.29 *0.58 **−0.01−0.060.61 **−0.090.040.83 **−0.110.02
GY−0.01−0.14−0.14−0.160.050.030.110.12
Fe0.56 **0.70 **0.54 **0.51 **
Across environments
TraitsGYZnFe
AD−0.090.43 **0.25 *
ASI−0.57 **0.34 **0.60 **
PH0.99 **−0.39 **−0.51 **
GY−0.44 **−0.43 **
Fe0.97 **
Ұ AD, Days to 50% anthesis; ASI, Anthesis silking interval (days); PH, Plant height (cm); GY, Grain yield (t ha−1); Zn, grain zinc concentration (mg kg−1); Fe, Grain iron concentration (mg kg−1). * p < 0.05, ** p < 0.01.
Table 6. Description of the selected plant materials used for making Zn-enhanced testcrosses and checks.
Table 6. Description of the selected plant materials used for making Zn-enhanced testcrosses and checks.
NO.CodeParent TypeNutritional TypeOriginZn (mg kg−1)Fe (mg kg−1)
1D2MaleZinc donorCIMMYT-Mexico33.8528.27
2D3MaleZinc donorCIMMYT-Mexico33.7228.98
3D5MaleZinc donorCIMMYT-Mexico30.3635.00
4D6MaleZinc donorCIMMYT-Mexico28.6831.52
5D7MaleZinc donorCIMMYT-Mexico32.1828.06
6D8MaleZinc donorIITA30.5234.40
7D9MaleZinc donorIITA30.2526.17
8D10MaleZinc donorIITA34.2928.77
9D11MaleZinc donorIITA30.0932.14
10D12MaleZinc donorIITA30.0226.72
11D13MaleZinc donorIITA27.2528.39
12NML1FemaleNormalCIMMYT-SARO34.3934.39
13NML3FemaleNormalCIMMYT-SARO28.3428.34
14NML5FemaleNormalCIMMYT-SARO30.1130.11
15PROA1FemaleProvitamin ACIMMYT-SARO30.8228.78
16PROA3FemaleProvitamin ACIMMYT-SARO28.1034.53
17QPM4FemaleQPMCIMMYT-SARO35.4833.65
18QPM6FemaleQPMCIMMYT-SARO29.1930.67
19C1CheckProvitamin ASeed companyNDND
20C2CheckProvitamin ASeed companyNDND
21C3CheckQPMSeed companyNDND
22C4CheckNormalSeed companyNDND
23C5CheckNormalSeed companyNDND
24C6CheckNormalSeed companyNDND
25C7CheckNormalSeed companyNDND
ND—Not yet determined; CIMMYT-SARO = CIMMYT Southern Africa; IITA = International Institute of Tropical Agriculture.
Table 7. Description of experimental sites, planting dates and crop management used.
Table 7. Description of experimental sites, planting dates and crop management used.
Site/YearAgro-EcologyLatitudeLongitudeAltitude
(masl)
Annual Rainfall (mm)ManagementEntriesPlanting TimeSoil Type
CIMMYT 2018/19 (E1)IIa17°48′ S31°03′ E1483850Optimum84 × 2 repsNovemberFerralsol
ART farm 2018/19 (E2)IIa17°42′ S31°5′ E1556850Optimum84 × 2 repsNovemberLixisol
DR&SS 2018/19 (E3)IIa17°13′ S31°03′ E1506850Low N84 × 2 repsNovemberFerralsol
Chisumbanje 2018/19 (E4)V20°47′ S32°13′ E480450Managed drought84 × 2 repsMayVertisol
CIMMYT 2019/20 (E5)IIa17°48′ S31°03′ E1483850Optimum84 × 2 repsNovemberFerralsol
ART farm 2019/20 (E6)IIa17°42′ S31° 5′ E1556850Optimum84 × 2 repsNovemberLixisol
DR&SS 2019/20 (E7)IIa17°13′ S31°03′ E1506850Low N84 × 2 repsNovemberFerralsol
Chisumbanje 2019/20 (E8)V20°47′ S32°13′ E480450Managed drought84 × 2 repsMayVertisol
Masl = metres above sea level; low N = managed low nitrogen conditions.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Goredema-Matongera, N.; Ndhlela, T.; van Biljon, A.; Kamutando, C.N.; Cairns, J.E.; Baudron, F.; Labuschagne, M. Genetic Variation of Zinc and Iron Concentration in Normal, Provitamin A and Quality Protein Maize under Stress and Non-Stress Conditions. Plants 2023, 12, 270. https://doi.org/10.3390/plants12020270

AMA Style

Goredema-Matongera N, Ndhlela T, van Biljon A, Kamutando CN, Cairns JE, Baudron F, Labuschagne M. Genetic Variation of Zinc and Iron Concentration in Normal, Provitamin A and Quality Protein Maize under Stress and Non-Stress Conditions. Plants. 2023; 12(2):270. https://doi.org/10.3390/plants12020270

Chicago/Turabian Style

Goredema-Matongera, Nakai, Thokozile Ndhlela, Angeline van Biljon, Casper N. Kamutando, Jill E. Cairns, Frederic Baudron, and Maryke Labuschagne. 2023. "Genetic Variation of Zinc and Iron Concentration in Normal, Provitamin A and Quality Protein Maize under Stress and Non-Stress Conditions" Plants 12, no. 2: 270. https://doi.org/10.3390/plants12020270

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop