Multi-Environment Quantitative Trait Loci Mapping for Grain Iron and Zinc Content Using Bi-parental Recombinant Inbred Line Mapping Population in Pearl Millet

Pearl millet is a climate-resilient, nutritious crop with low input requirements that could provide economic returns in marginal agro-ecologies. In this study, we report quantitative trait loci (QTLs) for iron (Fe) and zinc (Zn) content from three distinct production environments. We generated a genetic linkage map using 210 F₆ recombinant inbred line (RIL) population derived from the (PPMI 683 × PPMI 627) cross using genome-wide simple sequence repeats (SSRs). The molecular linkage map (seven linkage groups) of 151 loci was 3,273.1 cM length (Kosambi). The content of grain Fe in the RIL population ranged between 36 and 114 mg/Kg, and that of Zn from 20 to 106 mg/Kg across the 3 years (2014–2016) at over the three locations (Delhi, Dharwad, and Jodhpur). QTL analysis revealed a total of 22 QTLs for grain Fe and Zn, of which 14 were for Fe and eight were for Zn on three consecutive years at all locations. The observed phenotypic variance (R²) explained by different QTLs for grain Fe and Zn content ranged from 2.85 (QGFe.E3.2014–2016_Q3) to 19.66% (QGFe.E1.2014–2016_Q3) and from 2.93 (QGZn.E3.2014–2016_Q3) to 25. 95% (QGZn.E1.2014–2016_Q1), respectively. Two constitutive expressing QTLs for both Fe and Zn co-mapped in this population, one on LG 2 and second one on LG 3. Inside the QTLs candidate genes such as Ferritin gene, Al³⁺ Transporter, K⁺ Transporters, Zn²⁺ transporters and Mg²⁺ transporters were identified using bioinformatics approaches. The identified QTLs and candidate genes could be useful in pearl millet population improvement programs, seed, restorer parents, and marker-assisted selection programs.

Introduction

Pearl millet [Pennisetum glaucum (L). R. Br.] is a climate-resilient and health-promoting nutritious crop of the semi-arid tropics of Africa and Asia (Satyavathi et al., 2015). It is a staple food crop for approximately 90 million people that live in Asia and Africa and practice low-input subsistence farming and livestock production systems (Gangashetty et al., 2016). Pearl millet is an ideal crop in the context of global climate change due to its inherent tolerance to heat stress, salinity, and drought (Shivhare and Lata, 2016). Deficiencies or imbalanced intakes of energy or nutrients, particularly vitamins and minerals, cause a number of dysfunctions and diseases in humans, that are collectively referred to as “micronutrient malnutrition or hidden hunger.” Currently more than 2 billion people in the world are in risk of micronutrient malnutrition (Athar et al., 2020). Iron (Fe) and Zinc (Zn) deficiency in humans is caused by lack of purchasing power for highly nutritious food, reduced dietary intake and less bioavailability and bioutilization, especially among resource poor women, infants, and children in the developing countries (Shariatipour and Heidari, 2020). Nutrition supplementation, dietary diversification, and food fortification are among the targeted strategies available to address micronutrient deficiencies. Among, the above strategies, improving the nutrient profile of agricultural crops through genetic means is more pleasing as it is a one-time investment and becomes self-sustaining over a long period. Breeding of pearl millet varieties/hybrids with improved nutritional quality is one of the priority areas for providing nutritional security in the developing countries. Many breeders (Govindaraj et al., 2011, 2016; Rai et al., 2015; Satyavathi et al., 2015; Kumar et al., 2016) have identified huge variability for the trait and utilizing them in development of a handful of varieties and hybrids rich in micronutrients, which now playing a major role in nutritional security in dry lands. However, hybrid released so far in India is having as mean Fe content at par or less than HarvestPlus target level of 77 mg kg^–1 (HarvestPlus, 2014). Underutilization of germplasm material with higher nutritional value may attribute due to poor agronomic qualities. Hence there is a great need to make constant effort to breed more and more parental lines which can combine both for yield and micronutrient content along with resistance/tolerance to drought, downy mildew and blast diseases, good keeping quality, etc., and showing a stable performance under different agro-geographical regions.

Genetic maps provide important information for detailed genetic analysis of quantitative traits and have been proven to be significant method for crop improvement (Mohan et al., 1997; Doerge, 2002). Identifying molecular markers which are closely associated with the nutritional trait such as Fe and Zn enable rapid introgression of such traits into elite backgrounds. Among the different marker systems, simple sequence repeat (SSRs) is considered as a powerful and practically useful marker system for marker-assisted breeding (MAB) because of its cheap and more user-friendly along with co-dominant nature, abundance, reproducibility, and variability (Wang et al., 2018). In addition, the size and type of the mapping population determine the gene effect for an economically important trait (Kumar et al., 2018). Recombinant Inbred Lines (RILs) are the best choice among different types of mapping population as they are homozygous as well as have undergone multiple cycles of recombination which helps in mapping of tightly linked markers. Further use of codominant markers along with large RIL population helps in the construction accurate high-resolution map.

There are so many examples that SSR markers being employed for marker assisted biofortification of crop varieties such as QPM and pro-vitamin rich maize hybrids (Gupta et al., 2009; Muthusamy et al., 2014) and development of high yielding basmati varieties with disease resistance (Singh et al., 2012).

There are two different strategies that the plant can able to uptake iron from the rhizosphere: Strategy I and Strategy II (Connorton et al., 2017). Most of the cereals like rice, maize, and wheat follow Strategy II through secretion of phytosiderophores (PSs) (Roberts et al., 2004). TOM1 (Transporter of Mugineic acid family phytosiderophores)/ZIFL4 is a member of major facilitator superfamily (MFS) which exports Fe³⁺-PS chelates in rice and barley (Nozoye et al., 2011). In addition to Strategy II, rice and barley both have a functional Iron-Regulated Transporter 1 (IRT1) homolog that allows direct Fe²⁺ uptake from the rhizosphere. Metal transporters, such as the P1B-ATPase family, zinc-regulated transporter (ZRT), iron-regulated transporter (IRT)-like protein (ZIP), natural resistance-associated macrophage protein (NRAMP) family, and cation diffusion facilitator (CDF) family, have been reported in plants, where these two metals play critical roles (Colangelo and Guerinot, 2006). Several genes/gene families involved in Fe and Zn homeostasis such as yellow stripe-like (YSL), ZRT, ZIP, NRAMP, nicotianamine synthase (NAS), nicotianamine aminotransferase (NAAT), heavy metal ATPases (HMA), zinc-induced facilitator-like (ZIFL), metal tolerant protein (MTP), and others have been characterized in cereal crops such as rice (Anuradha et al., 2012), sorghum (Kotla et al., 2019). In pearl millet, PglYSL, PglZIP, and PglNRAMP and a single gene for PglFER and PglNAS were identified (Mahendrakar et al., 2020). Hence, the aim of this study was to identify the QTL associated with grain iron and zinc content through construct of linkage map using SSRs markers in a large-sized RIL mapping population.

Materials and Methods

Plant Materials

The experiment materials comprised of 210 F_6:7 RILs mapping population derived from an intra specific cross between PPMI 683 (P₁, Female) and PPMI 627 (P₂, Male) made during rainy season, 2009. The RIL population advanced by single seed descent method from F₂ to F₆ from single plants selected during F₂, and was segregating for grain Fe and Zn contents, respectively. Parent (P₁), i.e., PPMI 683 having high grain iron and zinc content and parent (P₂), i.e., PPMI 627 having low grain iron and zinc content with checks (ICMB 98222 for high grain Fe and Zn content, ICMB 92111 for low grain Fe and Zn content) (Table 1). Both checks were planted after every 20^th row of the population.

TABLE 1

Table 1. Pedigree of the parental lines used to develop F_6:7 mapping populations along with the high and low iron checks.

The average grain yield of five randomly selected plants was taken as grain yield per plant and weight of randomly counted 1,000 seeds in each genotype was taken as 1,000 Seed Weight. This study is mainly focused on grain Fe and Zn content. Hence, only the two most important traits viz., seed yield per plant and 1,000 seed weight were considered for the study.

Environments

A total of 210 RILs along with both parents and checks were evaluated in an alpha lattice design with two replicates (Patterson and Williams, 1976), under field conditions at three different geographical areas, showing all three pearl millet growing agro-climatic zones of India, (i) E1: ICAR–Indian Agricultural Research Institute, New Delhi (28°382′N, 77°802′E) representing Zone A with annual rainfall of more than 400 mm, (ii) E2: ICAR-IARI Regional Centre farm, Dharwad (15°212′N, 75°052′E) from zone B (covering the southern peninsular India), and (iii) E3: Agricultural Research Station, All India coordinated pearl millet project, Mandor, Jodhpur (26°252′N, 72°992′E) falling in zone A₁ with annual rainfall <400 mm. The trials were conducted at all locations during three consecutive (south-west monsoon season, 2014, 2015, and 2016) growing seasons. Each RIL was sown on a four-rows plot of 4 m long with an inter-row spacing of 65 cm and plant-to-plant spacing of 10 cm. For a normal, healthy crop, standard pearl millet cultivation practices were strictly followed. Table 2 shows the weather and soil parameters at each location throughout the year. All climatic parameters except rainfall are presented as means over the crop-growing season (June–October). Rainfall is cumulative rainfall received during the crop-growing season. The relative humidity was calculated as the average of measurements taken in every morning and afternoon on each day.

TABLE 2

Table 2. Geographical location, climatic, and edaphic factors over different years and locations.

Sample Preparation and Extraction of Mineral Elements

At physiological maturity, panicles from five plants per plot were randomly selected from the central two rows and harvested separately and threshed using wooden mallet. Seeds from individual plants within the plot were mixed and washed with deionized double-distilled water to remove debris and other possible contaminations in order to obtain a representative sample for trait measurements. The grain samples were sealed in individual paper bags and placed in an oven at 75°C overnight. About 15 g of air-dried sample from each genotype was further used. Mineral (Fe and Zn) content in the grains during initial year (2014) was measured with EDXRF instrument (a non-destructive, bench-top, energy-dispersive X-ray fluorescence spectrometry, Oxford Instruments X-Supreme 8000) at all three locations (Anuradha et al., 2017). Later, in 2015–2016, a more precise method, di-acid digestion of ground floor samples was used (Dhyan et al., 2005), followed by readings taken on an atomic absorption spectrometer (AAS, ZEE nit 700 tech Analytikjena). To determine whether higher grain mineral content is due to contamination with dust, soil, or is inherent, ICPMS (Nex ION 300X, Perkin Elmer, United States) was used to analyze genotypes having high iron content (>90 mg kg^–1) (Pfeiffer and McClafferty, 2007).

Statistical Analysis

Descriptive statistics and Analysis of Variance (ANOVA) was carried out using PBTools v1.4. Histograms and correlations between pairs of traits were estimated through Pearson correlation co-efficient using R software. The model used for ANOVA was:

$$z_{ijklm}=\mathrm{\mu }+y_i+e_j+y_{eij}+r_{ijk}+b_{ijkl}+g_m$$

$$+{(yg)}_{im}+{(eg)}_{jm}+{(yeg)}_{ijm}+{}^{\prime }{}^{\prime }\mathrm{\Sigma }_{ijklm}$$

Where μ is the grand mean; yi is the fixed effect of year i; ej is the fixed effect of location j; ye_ij is the fixed effect of interaction between year i and location j; gm is the random effect of genotype m and is ∼NID (0, σ²g); r_ijk is the random effect of replication in location j and year i and is ∼NID (0, σ²r); b_ijkl is the random effect of block l nested with replication k in location j and year i and is ∼NID (0, σ²_b); (yg)_im is the random effect of the interaction between genotype m and year i and is ∼NID (0, σ²_yg); (eg)_jm is the random effect of the interaction between accession m in location j and ∼NID (0, 0, σ²_eg); (yeg)_ijm is the random effect of the interaction effect of the genotype m in year i and location j and ∼NID (0, 0, σ²_yg); and ε_ijklm is the random residual effect and ∼NID (0, 0, σ²_ε).

Analysis of variance was also conducted using data from each environment for grain Fe and Zn content.

Heritability (H², broad sense) at individual environment was estimated from analysis of variance. The formula used was-

$$H^2=\frac{{\mathrm{\sigma }}_g^2}{{\mathrm{\sigma }}_g^2+\frac{{\mathrm{\sigma }}_{\mathrm{\epsilon }}^2}{r}}$$

Whereas Heritability (H²) estimates across environments were estimated by the formula-

$$H^2=\frac{{\mathrm{\sigma }}_g^2}{{\mathrm{\sigma }}_g^2+\frac{{\mathrm{\sigma }}_{yg}^2}{y}+\frac{{\mathrm{\sigma }}_{eg}^2}{l}+\frac{{\mathrm{\sigma }}_{yeg}^2}{yl}+\frac{{\mathrm{\sigma }}_{\mathrm{\epsilon }}^2}{ryl}}$$

Where r, y, l denotes the number of replicates, years and environments, respectively.

The REML model also produced best linear unbiased predictors (BLUPs) of each genotype, thereby adjusting the influence of the neighboring rows. These BLUPs were used for downstream analysis.

Genomic DNA Extraction and Genotyping

The genomic DNA was extracted from fresh and young leaf tissue of each RIL and the parents by using modified cetyltrimethyl-ammonium bromide (CTAB) method (Aboul-Maaty and Oraby, 2019). Polymorphism was assessed using 372 SSR markers, EST-SSRs (ICMP and IPES series), gSSRs (PSMP and CTM series), and gene-based markers. SSRs were amplified in 10 μl PCR reaction mixture containing 25–30 ng genomic DNA, 100 μM of each of the four dNTPs, 1.5 mM MgCl₂, 5 pmol/μl of forward and reverse primer, 1 U of Taq DNA polymerase (Banglore-genei, Merck) and 1× PCR buffer in Biorad My Cycler Thermal Cycler, United States. PCR conditions were used as described by Anuradha et al. (2017). The amplified product was checked using 1.2% agarose gel and PCR products were separated by 10% polyacrylamide gel electrophoresis at 120 V (C.B.S. Scientific vertical electrophoresis unit). The SSR marker alleles were scored manually by comparing the position of bands with that of 100 bp ladder from bottom to top (smallest to largest). Based on the position, the allele size in base pairs (bp) was recorded for all polymorphic SSRs in all 210 RILs studied. The parent P₁ (PPMI 683) was scored as “A” and parent P₂ (PPMI 627) was scored as “B” and both type band present on single locus were designed as “H”. Diffused or ambiguous bands were recorded as missing bands with “0” notation. The chi-square method was conducted to test whether the inherited alleles of the mapping population were consistent with Mendelian segregation ratios. The segregation ratio across the mapping population was tested against a 1:1 ratio. The segregation pattern that did not fit either ratio (p < 0.05) were treated as distorted.

Linkage Analysis and QTL Mapping

The linkage map was constructed using 151 polymorphic SSR markers. The markers that were not taken for the linkage analysis were those with: a minor allele frequency of less than 5%, more than 20% missing data, and those that were monomorphic between the parents. The MAPMAKER/EXP3.0 software was used to develop a genetic linkage map (Lander et al., 1987). The maximum recombination fraction (θ) was set at 0.49 and the critical logarithm of odds (LOD) score for the test of marker pair independence was set at 3. Markers with high segregation distortion (p ≤ 0.01) were removed in a preliminary analysis of the data. Subsequently, the linkage analysis of genotyping data for all the markers including those with distorted segregation was used. Using MapChart software, a graphical map was generated (Voorrips, 2002). The consensus map developed by Supriya et al. (2011) and Rajaram et al. (2013) was used as a reference map in the current study for grouping markers and constructing the map. QTL analysis was carried out using the Inclusive Composite Interval Mapping method (ICIM). The ICIM approach uses a strategy in which a stepwise regression is firstly performed, so markers with significant effect on QTL are selected. ICIM-ADD method was used with the multi-environmental model builtin QTL IciMapping (Li et al., 2015). Significant LOD thresholds were set at 5% tail of the null distribution in a 1,000 permutations test (Luciano Da Costa et al., 2012).

Quantitative Trait Loci Validation

All QTL intervals identified for grain Fe and Zn content were analyzed in-silico with reference genome, Tift23D₂B₁–P₁–P₅ (PROJECT ID: PRJNA294988.) in order to navigate the presence of candidate gene associated with Fe and Zn content in grain. The physical position of flanked markers of identified QTLs were obtained from NCBI-BLASTn¹. The sequence similarities searched against related cereals (Setaria italica) by Ab initio gene prediction method by using Fgenesh software for the identification of potential protein coding regions (Salamov and Solovyev, 2000)². For each predicted exon similarity searches were run using the BLAST program on protein and EST databases (Altschul et al., 1997).

Results

Phenotypic Evaluations

The mean performance of parents and range of trait values of RILs during three consecutive years (2014, 2015, and 2016) at E1, E2, and E3 are summarized in Table 3. Both traits (Fe and Zn) were approximately normally distributed when analyzed across the years at their respective locations (Figure 1). The levels of grain Fe in RILs ranged from 43.27 to 108.33 mg kg^–1, 33.34–118.52 and 32.88–115.47 mg kg^–1, while the range of grain Zn content was 22.4–94, 20.97–120.87, and 17.33–105.17 mg kg^–1 across the years at E1, E2, and E3, respectively. The mean BLUP value of grain Fe and Zn content in both parents and RILs were higher across the described years at E1 when compared to E2 and E3, whereas mean BLUPs of grain Fe and Zn for both parents was lowest at E3 across the years. For grain Fe and Zn studied in different environments, both parents and RILs showed a wider range of variation.

TABLE 3

Table 3. Grain Fe and Zn content (mg Kg^–1) ±SE of the parents along with RILs of F_6:7 mapping populations over 3 years and three diverse locations.

FIGURE 1

Figure 1. Histograms of grain Fe and Zn in mapping populations of recombinant inbred lines evaluated across the year (2014–2016) at Delhi, Dharwad, and Jodhpur location.

The seed yield per plant and 1,000 seed weight of the 210 RILs along with both parents varied greatly. At E1 during three consecutive years (2014–2016), 1,000 seed weight ranged from 6.75 to 12.17 g with a mean of 9.05 g, and seed yield per plant ranged from 14.09 to 36.21 g with the mean of 24.05 g. And at E2 from 2014 to 2016, 1,000 seed weight ranged from 6.77 to 12.33 g with a mean of 9.09 g, and seed yield per plant ranged from 19.05 to 35.37 g with a mean of 26.81 g. From 2014 to 2016, 1,000 seed weights in E3 ranged from 6.48 to 11.41 g, with an average of 8.38 g, and seed yield per plant was 18.20–39.8 g, with an average of 26.75 g.

The ANOVA indicated a highly significant G × Y interaction [p < 0.0001, F_(66.45) = 1.112 for Fe] and [p < 0.0001, F_(62.43) = 1.112 for Zn]. The pooled mean, range, coefficient of variation (CV), genotypic variance (σ²_g), residual variance, broad sense heritability (H²), least square difference (LSD) and G × E variance in RILs evaluated during 2014–2016 cropping season at E1, E2, and E3 has shown in Table 4. In the RIL population, the genotypic components of variance (σ²_g) for both traits were significant, and the operational heritability (H²) estimates ranged from 0.83 (Zn) to 0.88 (Fe) (Table 4). The high heritability shows that much of the phenotypic variance in the population is genetically controlled and QTL can be mapped with a high degree of reliability. Broad-sense heritability was usually high enough to allow effective QTL mapping, indicating that the traits studied have moderately high to high proportions of genetic variance. Correlation coefficients (r) for grain Fe and Zn content were found statistically significant (p < 0.001) for each year of evaluations for each location.

TABLE 4

Table 4. Estimates of mean, variance components (2014, 2015, and 2016), broad-sense heritability and LSD at Delhi, Dharwad, and Jodhpur for grain Fe and Zn content in RILs.

Homogeneity of variance tests based on Bartlett’s test indicated homogeneous error variance for both the traits which allowed for a combined analysis across the years for all locations. When performing a combined ANOVA, the year was considered as a random and genotypes were considered as fixed effects. Combined ANOVA showed significant Genotype × Environment Interaction (GEI) (p < 0.001), exhibiting the influence of changes in environment on the grain micronutrient contents of genotypes. It was observed that Genotypes (G) and GEI effects were highly significant (p < 0.01) and the contribution of year is very less for both the traits (Table 5).

TABLE 5

Table 5. Mean sum of square of Fe and Zn at all locations based on combined analysis over three years for RILs and parents.

Correlation of Grain Fe, Zn, 1,000 Seed Weight and Yield

Grain Fe and Zn content were highly (0.711) and significantly (p < 0.01) correlated (Table 6). Thousand seed weight was also observed to be positively correlated with seed yield per plant (0.579, p < 0.01). For successful incorporation of high levels of both grain Fe and Zn content, association between these micronutrients and their association with thousand seed weight and yield is important to plan a successful breeding approach so as to develop new improved lines or varieties having higher grain yield along with enhanced levels of grain micronutrients.

TABLE 6

Table 6. Correlation of grain iron and zinc content with thousand seed weight and grain yield.

Linkage Map Construction

For the polymorphism survey of both parental lines, 366 SSR primer pairs (IPES, PSMP, and ICMP series) and six gene-based primers were used. Of them, 151 (40.59%) SSRs detected polymorphism between the two parents. Polymorphic SSR markers included 94 IPES, 42 PSMP, and 11 ICMP series primer pairs and four newly synthesized primer pair named as Ppmsb (Pusa pearl millet SSR biofortification) (Figure 2). These polymorphic primers were further used for linkage map construction with LOD score 3.0.

FIGURE 2

Figure 2. Linkage map (linkage groups LG 1–LG 7) of the pearl millet recombinant inbred line (RIL) population based on the cross (PPMI 683 × PPMI 627). QTL identified in red color font indicates the grain Fe QTLs and QTL identified in blue color font indicates the grain Zn QTLs.

The total length of the linkage map was 3273.1cM (Kosambi). The individual linkage group ranged from 910 cM for LG 2 with the highest number of markers (31) to 135.3 cM for LG 3 with 14 markers, while the least number of markers were assigned on LG4 (10) followed by LG 3 (14) and LG 5 (14). The average interval size was 8.79 cM. The average LG length was 467.5 cM with an average of 21.57 loci. The average adjacent-marker interval lengths ranged from 3.07 cM (LG 3) to 14.91 cM (LG 2) followed by 11.34 cM (LG 1). Map distance between adjacent markers varied from 0.9 to 49.2 cM (LG 1) and 8.33% of the intervals (31 out of 372 intervals) were smaller than 10 cM. Details of all linkage group, total markers, skewed marker and their percentage, total map length (cM) and average distance (cM) were given in Table 7.

TABLE 7

Table 7. Details of the simple sequence repeats (SSRs)-based linkage map of the pearl millet RIL population based on the cross PPMI 683 × PPMI 627.

Quantitative Trait Loci Analysis

The multi-environmental QTL analysis exhibited the association of various genomic regions with grain Fe and Zn content at all the three locations during three consecutive years, 2014, 2015, and 2016 (Table 8). The analysis indicated the presence of four significant common QTL for grain Fe content, flanked by following markers Ppmsb001–Ppmsb002, IPES0142–IPES0180, IPES0157–IPES0093, and IPES0206–IPES0015, which explained the total phenotypic variance of 18.34, 17.98, and 15.98% for QTL1 on LG 2, 14.99, 16.26, and 14.8% for QTL 2 on LG 3, 19.66, 18.83, and 16.85% for QTL 3 on LG 5, 15.66, 16.61, and 14.56% for QTL 4 on LG 7 across three consecutive years (2014–2016) at E1, E2, and E3, respectively. These four QTLs were named as: QGFe.E1.2014–2016_Q1, QGFe.E1.2014–2016_Q2, QGFe.E1.2014–2016_Q3, and QGFe.E1.2014–2016_Q4 for Delhi location (E1), QGFe.E2.2014–2016_Q1, QGFe.E2.2014–2016_Q2, QGFe.E2.2014–2016_Q3, and QGFe.E2.2014–2016_Q4 for Dharwad location (E2), and QGFe.E3.2014–2016_Q1, QGFe.E3.2014–2016_Q2, QGFe.E3.2014–2016_Q3, and QGFe.E3.2014–2016_Q4 for Jodhpur location (E3) (Q = QTL; G = grain; Fe = iron). In addition, two more QTLs for GFe were also found at E3, flanked by IPES0101–IPES0139 and IPES0027–IPES0236 and having very less PVE of 2.92 and 2.85% on LG 1 and LG 2, respectively. The genome wide significant threshold was LOD = 3.5 for E1, E2, and E3.

TABLE 8

Table 8. Quantitative trait loci for grain Fe (GFe) and Zn (GZn) content in RIL over three locations.

For GZn, one common QTL was found in E1, E2, and E3 on LG 3 at significant threshold LOD value of 3.5 which was flanked by IPES0142–IPES0180 with total phenotypic variance of 14.39, 16.07, and 22.26% at E1, E2, and E3, respectively. And one more common QTL was also found at E1 and E2 on LG 2, flanked by Ppmsb001–Ppmsb002 having total PVE of 25.95 and 21.88%, respectively. Whereas, at E3, one QTL with flanked markers Ppmsb002–IPES0181was mapped on LG 2 with comparatively low total PVE of 11.18%. Two more un-repetitive QTLs between IPES0166–PSMP2249 (PVE 6.2%) on LG 3 and IPES0141–IPES0151 (PVE 2.93%) on LG 6 were detected in E1 and E3, respectively.

The proportion of the total phenotypic variance (PVE) explained by the QTLs for GFe ranged from 14.99 to 19.66%, 16.26–18.83, and 2.85–16.85% in E1, E2, and E3, respectively (Table 8). On the other hand, the PVE of the QTLs for GZn in E1, E2, and E3 ranged from 6.2 to 25.95%; 16.07–21.88 and 2.93–22.26% (Table 8). The QTLs that had the highest PVE for GFe in E1, E2, and E3 were 19.66% (IPES0157–IPES0093), 18.83% (IPES0157–IPES0093), and 16.85% (IPES0157–IPES0093), respectively. For GZn, the QTLs with the highest PVE were 25.95% (Ppmsb001–Ppmsb002), 21.88% (Ppmsb001–Ppmsb002), and 22.26% (IPES0142–IPES0180) at E1, E2, and E3, respectively.

We localized two genomic regions with pleiotropic effects on LG 2 and LG 3 for GFe and GZn at E1, E2, and E3. The first genomic region was QGFe.E1.2014–2016_Q1 to QGZn.E1.2014–2016_Q1 at E1, QGFe.E2.2014–2016_Q1 to QGZn.E2.2014–2016_Q1 at E2 and QGFe.E3.2014–2016_Q2 to QGZn.E3.2014–2016_Q1 at E3 linked to marker Ppmsb1–Ppmsb2 on LG 2 (Table 8). The second region was associated with QGFe.E1.2014–2016_Q2 to QGZn.E1.2014–2016_Q3, QGFe.E2.2014–2016_Q2 to QGZn.E2.2014–2016_Q2, and QGFe.E3.2014–2016_Q4 to QGZn.E3.2014–2016_Q2 linked to marker IPES0142–IPES0180, across the years, at all the three locations.

The QGFe.E1.2014–2016_Q3 region linked to marker IPES0157–IPES0093 having largest PVE for GFe (19.66%) on LG 5 over the locations, which was also displayed 3.6 additive effect of QTL (Table 8). Whereas, the QTL for GZn having largest PVE (25.95%) was QGZn.E1.2014–2016_Q1 on LG 2 over the locations (Table 8). Similarly, one more QTL for grain Zn (QGZn.E2.2014–2016_Q1) was also found with high PVE 21.88% on LG 2 at E2. The additive effect was detected for all the QTLs on which the alleles from PPMI 683 (donor parent) had positive additive effect.

Validation of Identified QTL in Diverse Set of Germplasm

Eleven high (PPMI 708, PPMI 1102, PPMFeZMP 199, PPMI 1278, 841B, PPMI 1232, PPMI 1107, PPMI 1225, PPMI 1089, PPMWGI 146, and PPMI 1231), nine low (PPMI 1011, PPMI 759, J 2467, PPMI 1087, PPMFeZMP 47, PPMFeZMP 34, J2405, and H77/833-2, 5540B) and four checks (ICMB 98222, ICMB 93222, ICMB 07999, and ICMB 92111) for grain iron expressing pearl millet genotypes (Table 9) were selected and their DNA was used for amplification with ten consistently associated primers, IPES0142, IPES0180, IPES 0160, IPES0093, IPES0206, IPES0015, IPES0166, PSMP2249, Ppmsb001, and Ppmsb002. The amplicons of these genotypes were obtained by PCR amplification of these primers. The similar regions of reference sequence obtained were thus compared with the amplified sequences of high and low grain iron of expressing pearl millet genotypes. These markers were compared with reference genome Tift₂₃D₂B_1.

TABLE 9

Table 9. List of pearl millet genotypes used to validate the flanked markers of QTLs for grain iron and zinc content.

In-silico Identification of Candidate Gene

Candidate genes identified for each QTL was shown in Table 10. In the study, we identified ten metal transporter genes. Out of the ten genes seven were identified in the genome of Seteria italica whereas remaining three genes were identified in the genome of Zea mays. Among the genes from the Seteria italica, Ferritin-1 was associated with the QTL QFe 2.1 and QZn 2.1, Aluminum-activated malate transporter 5 and Potassium transporter 25 was identified with QFe 3.1 and QZn 3.1, Probable magnesium transporter NIPA4, Metal transporter Nramp3 isoform X2, Probable cadmium/zinc-transporting ATPase HMA1, chloroplastic isoform X2 and Cadmium/zinc-transporting ATPase HMA2 were linked with QFe 5.1. However, among the genes from the Zea mays, Potassium transporter 3 was associated with QFe 5.1 and QFe 7.1.

TABLE 10

Table 10. Identification of putative candidate gene within QTL involved in enhancing Fe/Zn content.

Discussion

Variance Components and Heritability

Pearl millet has a great variation of mineral and agronomic traits. For both traits, the parents of this RIL population exhibited statistically significant divergent phenotypes. The BLUP means of RILs were significant for Fe and Zn content. A very wide range of variation was detected among the population for grain Fe and Zn content (Table 3).

Wider range of variation existed in the present mapping population for grain yield per plant (18.89–36.70 g) and thousand seed weight (6.84–11.79 g). Wide range of variation was earlier reported by many workers in pearl millet for yield and micronutrient traits (Govindaraj et al., 2011; Sankar et al., 2013; Vinodhana et al., 2013; Bashir et al., 2014; Sabiel et al., 2014; Anuradha et al., 2017). Greater amount of variation in the present population show that there is a lot of scope for selection of better lines to increase the grain yield in pearl millet.

The operational heritability (H²) estimates were very high and the genotypic components of variance (σ²g) for both traits were significant. Similar study was also done by Kumar et al. (2016, 2018). The trait heritability information is helpful to choose the selection procedure to be followed to improve the trait. Broad-sense heritability was high enough to allow effective QTL mapping, with moderately high to high proportions of genetic variance for the traits investigated. Higher estimates of heritability with genetic advance as a percent of mean were observed, suggesting the presence of additive gene action efficacy of selection for both traits. The traits which expressed high heritability and low genetic advance showed non additive gene action, hence heterosis breeding would be recommended for these traits. Sumathi et al. (2010) and Govindaraj et al. (2011) were also found similar result.

Frequency Distributions

The frequency distributions of the overall means for grain iron and zinc content in the mapping population are represented by histogram (Figure 1). Likewise, previously Kumar et al. (2016, 2018) were also concluded the normal frequency distribution for both traits in pearl millet. The normal Frequency distribution for both traits was also observed in other crop like wheat (Crespo-Herrera et al., 2017). Both traits in the RIL population have phenotypic normal distributions and transgressive segregations, indicating polygenic inheritance. The presence of transgressive segregation also suggests genetic recombination (Falconer, 1996), which means that both favorable and unfavorable alleles for the traits are dispersed between the parents.

Correlation Analysis

The understanding of the association between nutritional and agronomic traits will support the breeders to select suitable selection/breeding program for genetic improvement of associated traits such as mineral content. Several crops have been studied to determine the correlation between grain Fe and Zn contents. Iron and zinc didn’t record any significant association with grain yield and thousand seed weight albeit they in turn are correlated. In many previous studies there was negative association of Fe and Zn with yield (Kanatti et al., 2014), which hampers the effective selection of all these traits. This study is in accordance with Anuradha et al. (2017), where they didn’t notice any associations. Hence, there is a possibility for selection of high Fe and Zn content along with high yield. In current study, positive correlation between iron and zinc was observed as also reported in many previous studies (Rai et al., 2012; Govindaraj et al., 2013; Kanatti et al., 2014; Kumar et al., 2016, 2018; Anuradha et al., 2017, 2018; Velu et al., 2017). This could be due to common molecular mechanisms that regulate mineral uptake and metabolism in grains, or to common transporters controlling mineral movement within plants (Vreugdenhil et al., 2004; Ghandilyan et al., 2006; Magallanes-López et al., 2017). The strong association between these minerals in populations may be due to the co-segregation of genes for these traits. The direction and intensity of the association suggest that co-transferring superior alleles regulating these traits into the genetic backgrounds of elite lines could provide good opportunities for simultaneous genetic improvement of both micronutrients (Allouis et al., 2001). Genotypic correlations were found to be higher for all traits, implying that there is less interaction between the genetic makeup of the traits and environment.

Linkage Map

Linkage mapping is essential for QTL mapping and marker-assisted breeding programs. The most preferred mapping population is RIL, but the accuracy of the mapping resolution is dependent on population size. The linkage maps have been predominantly constructed using F₂ populations and up to certain extent RILs in pearl millet (Ferreira et al., 2006). Only few previous studies were done on linkage map and QTLs in pearl millet like Supriya et al. (2011) had constructed the linkage map based on 321 loci (258 DArTs and 63 SSRs) for 140 RILs. Similarly, Kumar et al. (2016, 2018) had also used DArTs and SSRs for construction of linkage map and used a smaller number of SSRs for genotyping of RILs and map development. In our study, we found 151 polymorphic SSRs out of 366 and used for linkage map development for 210 RILs. The mapping population size and number of SSRs was high compared to the Supriya et al. (2011) and Kumar et al. (2016, 2018). The present study recorded lesser polymorphism for SSRs compared to that of Senthilvel et al. (2008) who reported 74%. It may be because the parents used here for the development of RILs are closely related than that of Senthilvel et al. (2008). The distribution of SSRs across the linkage group was uniform in the current study. In terms of SSR marker positioning on the genetic map, all SSRs mapped on the similar linkage group (LG) as previously reported pearl millet maps by various researchers (Yadav et al., 2004; Senthilvel et al., 2008; Supriya et al., 2011; Kumar et al., 2016, 2018). One hundred fifty-one polymorphic SSR markers were assigned to seven previously established pearl millet linkage groups (Figure 2). The present map had a larger length than the one reported by Yadav et al. (2004); Supriya et al. (2011) and Kumar et al. (2016, 2018).

Recombinant inbred lines usually exhibit segregation distortion, because, many recessive lethal genes become homozygous during the RIL development process (Kumar et al., 2016). Segregation distortion was found in 31.57% of the loci on the current linkage map. Whereas Kumar et al. (2018) found 60% segregation distortion of DArts and SSRs. Pollen abortion is most common in pearl millet than abnormalities in female gametes, leading to a comparatively greater loss of male alleles and the subsequent skewness toward the female parent (Kumar et al., 2018). Distortion from expected Mendelian segregation has been observed previously in maize (Wendel et al., 1987; Lu et al., 2002), barley (Graner et al., 1991; Devaux et al., 1995), rice (Causse et al., 1994; Xu et al., 1997), wheat (Blanco et al., 2004; Quarrie et al., 2005), and pearl millet (Supriya et al., 2011). The protogynous nature of pearl millet also contributes to segregation distortion (Liu et al., 1994). According to Cloutier et al. (1997) and Livingstone et al. (1999), residual heterozygosity and inbreeding depression during inbred development may also contribute to segregation distortion. The residual heterozygosity in some RIL may be advantageous at the same time, because deleterious genetic combinations in the form of reduced fitness or lethality can be avoided. The segregation of nearly all loci of LG 3 was distorted in comparison to earlier reports, which may explain part of the increase in LG 3 map length.

The presence of large gaps between the centromere and telomere is a characteristic feature of pearl millet linkage maps. Senthilvel et al. (2008) found a large gap in LG 4. The previously built framework map for the cross ICMB 841-P3 × 863B-P2 had big gaps in LG 2 and LG 7 (Senthilvel et al., 2008; Yadav et al., 2011). This new population also revealed a large gap (>25 cM) in LG 1, LG 2, LG 4, and LG 6, possibly due to extreme recombination localization at the ends of LG. The interval distance between the markers in LG 4 was high suggesting additional polymorphic markers are needed to fill-up the gap. According to many researchers like Devos et al. (2000); Varshney et al. (2006), and Senthilvel et al. (2008), large gaps in the distal regions reflect regions of high recombination, rather than a lack of markers in these regions. However, these linkage groups are still incomplete, and need the addition of new markers that are located on the distal regions of the linkage groups. The number of large gaps has tried to decrease in the present study, although there is still the possibility to map more markers to fill these gaps.

Quantitative Trait Loci for Grain Iron (Fe) and Zinc (Zn) Content and Environment Interaction

The use of QTL analysis to map populations is useful not only for identifying genomic regions associated with traits of interest, but also for using the associated marker information in breeding programs to integrate particular loci in elite germplasm. The application of QTL by environment interaction with composite interval mapping was used in the present study (Li et al., 2015), it was possible to locate several genomic regions associated with GFe and GZn during the three consecutive years over three diverse locations. Our findings support other authors’ findings that GFe and GZn are quantitative traits (Kumar et al., 2016; Anuradha et al., 2017). These previous QTL studies have also mapped QTL or genomic region for GFe and GZn in different linkage group of pearl millet, including LG 3 and LG 5, with PVE 13.6 and 10.0%, respectively, for iron content (Kumar et al., 2016). on the other hand, Anuradha et al. (2017), identified different markers on different linkage groups for both high grain iron and zinc content, namely PSMP 2261 (LG 5), IPES0096 (LG 7), IPES 0180 (LG3), and Xsinramp 6 (unmapped). Kumar et al. (2016) reported the co-localization of high grain iron and zinc content alleles/QTLs in pearl millet. It may happen because starting from uptake to final deposition into grain of iron and zinc may share some common pathways (Grotz and Guerinot, 2006). In present study we found four repetitive QTL on LG 2, LG 3, LG 5, and LG 7 with PVE ranging 16.0–18.6, 14.83–16.26, 16.88–19.63, and 14.59–16.61% in E1, E2, and E3, respectively, and the largest PVE (19.63%) was displayed by QGFe.E1.2014–2016_Q3 in LG 5. Kumar et al. (2016) has mapped one SSR marker (IPES142) on LG 3 and Anuradha et al. (2017) has mapped IPES180 on LG3. In the present study one QTL was mapped between IPES142–IPES180 on LG 3 which is in accordance with that of Kumar et al. (2016) and Anuradha et al. (2017). Further, Anuradha et al. (2017) mapped a marker ICMP 3092 on LG 7 for Fe. Similarly, we also identified one repetitive QTL between IPES0206–IPES0015 on LG 7 and interestingly the nearest marker to this QTL is ICMP 3092 as per Rajaram et al. (2013). However, the genomic regions of this study and other studies are not comparable since Kumar et al. (2016) used DArT markers while this study used SSR markers. Similarly, the mapping method used by Anuradha et al. (2017) was association mapping while it is QTL in the present study.

The genotype by environment interaction, particularly the cross-over type interactions, has significant implications for crop performance and breeding. The ideal case for pearl millet biofortification is to obtain stable pearl millet genotypes that perform well without cross-over interaction when tested in different environments or years in a particular geographical area. The analysis of the phenotypic data in our study revealed the presence of a significant G × Y interaction. QTL analysis revealed that most of the LOD scores for the additive average effect were higher than the LOD score for the interaction (Table 7), showing that QTL with higher LOD (Add) are more stable than those with higher LOD (G Y) (Li et al., 2015). The additive effect was detected for all the QTLs for grain Fe and Zn, on which the alleles from donor parent had positive additive effect. Gaoh et al. (2020) found the additive × additive gene effects had the most important effects for grain Fe content, while additive × dominance gene effects were significant for grain Zn content in pearl millet. Similarly, Garcia-Oliveira et al. (2018) and Kumar et al. (2018) reported that epistasis plays a role in regulating grain iron and zinc content in pearl millet and cereals, respectively. Lu et al. (2008) reported 28 genome-wide additive × additive interactions for mineral elements in rice. Anuradha et al. (2012) also reported epistatic interactions between loci on chromosomes 1 and 5 for grain Zn concentration. Jin et al. (2013) reported the additive effect for Zn and partial dominant for Fe in maize.

Validation of Identified QTLs in a Diverse Set of Germplasm

Validation of marker trait associations can reveal their worthiness for further utilization in breeding. Hence, DNA from eleven high, nine low grain iron and four check (two for high and two for low) expressing genotypes were amplified with ten SSR markers which were consistently associated with both the traits. Flanked SSR markers of identified QTLs were compared with the reference Tift₂₃D₂B₁ genome.

In silico Identification of Candidate Gene

The availability of the pearl millet genome sequence has allowed researchers to search for candidate genes involved in Fe and Zn accumulation in grains. The presence of candidate genes was investigated in genomic regions encompassing all QTLs contributing to Fe and Zn in pearl millet. QTLs identified for Fe and Zn content using composite interval mapping were analyzed In-silico for the presence of putative candidate genes. In the study, we identified ten metal related genes within the obtained QTLs.

Foods with higher iron content can be produced using modern genetic and molecular technologies. To increase the amount of iron and zinc in wheat, ferritin is over-expressed (Liu et al., 2016). Other mineral concentrations were also affected by ferritin expression, as reported by several studies (Drakakaki et al., 2005; Qu et al., 2005). The promoter selected to control ferritin expression is critical to iron accumulation in specific tissues. The Fe and Zn concentrations in transgenic indica rice grain improved by endosperm-specific glutelin gene promoter of soybean ferritin (Vasconcelos et al., 2003). According to Liu et al. (2016), over-expression of sickle alfalfa ferritin, which is controlled by the seed-storage protein glutelin GluB-1 gene promoter, increases grain Fe and Zn concentrations but also affects mineral homeostasis in transgenic wheat grain. Wheat TaALMT1 (ALMT, for Al-activated Malate Transporter) encoding a malate transporter involved in Al tolerance (Raman et al., 2005). Over expression of TaALMT1 in wheat, barley, and tobacco-cell suspension increases the efflux of Al-activated malate and increases the tolerance to Al stress (Pereira et al., 2010). Root exudation of organic acid such as malate/citrate may play an important role in providing iron to plants (Brown, 1983). Citrate efflux is crucial for iron translocation, and this is mediated by the efflux transporter FRD3 in Arabidopsis (Green and Rogers, 2004) and its orthologue FRDL1 in rice (Yokosho et al., 2016). Barak and Chen (1984) reported a positive effect of K on iron absorption which was associated with acidification of the rhizosphere. Sakaguchi et al. (1999), observed secretion of mugineic acid family phytosiderophores (MAs) from the use of potassium gradient, and content of potassium in barley roots increased with iron-deficiency. Grotz et al. (1998) identified the zinc transporter genes in Arabidopsis thaliana. Yang et al. (2020) studied on activity of zinc transporter (OsZIP9) mediates zinc uptake in rice. Plants acquire Mg from the environment and distribute in the ionic form via Mg²⁺-permeable transporters. In Arabidopsis, MGT6, (plasma membrane-localized magnesium transporter) mediates Mg²⁺ uptake under Mg-limited conditions (Yan et al., 2018). The Magnesium/proton exchanger protein is a vacuolar exchanger of protons with cytosolic Mg²⁺ and Zn²⁺ in Arabidopsis. Hence, it may be possible that Mg and Zn share a common pathway of transport in plants.

Ferritin Gene (Ppmsb001–Ppmsb002)

Plant ferritins are important iron storage proteins that serve as both an acceptor and a donor of iron in metabolic processes (Briat et al., 2010). It shares structurally and functionally similarities with animal ferritins. Plant ferritin is generally observed in cells which are inactive in photosynthesis such as seed, root, root nodules, and others. Plant ferritins are nuclear encoded and located primarily in plastids (Zielińska-Dawidziak, 2015). This protein also controls the cellular concentration of transition metals, metal ions apart from iron in the mineral core (including cadmium, beryllium, zinc, aluminum, and lead) (Kumar and Prasad, 1999). Many gene families like FER, ZIP, NRAMP, YSL, and NAS, were involved for mineral translocation from soil to the grain filling site. These gene families have been identified in plants such as Arabidopsis thaliana (Weber et al., 2004), Setaria italica (Alagarasan et al., 2017), and Oryza sativa (Neeraja et al., 2018). In the current report, the ferritin gene was reported by in silico method, which is similar to the Ferritin-like (PglFER1) gene family previously reported by Mahendrakar et al. (2020). With reference to previous studies it is clear that ferritin gene regulates the iron metabolism in plants.

Al³⁺ Transporter in Plants (IPES0142–IPES0180)

Aluminum is a potentially phytotoxic metal and Al tolerant plants excrete organic acids such as malate, citrate, and oxalate from root tips, depending on the plant (Ryan et al., 2011). The organic acid exudation is vital for plants tolerating metal and nutritional stress at the root-soil interface. Aluminum interferes the uptake or transport of nutrients such as Ca, B, Fe, Mn, P, Mg, and Cu or K (Keltjens and Tan, 1993; Keltjens, 1995; Lukaszewski and Blevins, 1996; Ślaski et al., 1996). Root elongation and nutrient (B, Fe, Mg, Ca, and P) absorption may also be inhibited by Al or low pH soils. Root exudation of organic acid such as malate/citrate may play a key role in supplying Fe to plants (Brown, 1983). Hence, with the consideration of previous studies, we assumed that the aluminum activated malate transporter gene is responsible to cease the aluminum absorption by plant root and also induce the iron absorption.

K⁺ Transporters in Plants (IPES0142–IPES0180, IPES0157–IPES0093, and IPES0206–IPES0015)

Oertli and Opoku (1974), as well as Barak and Chen (1984), identified a positive effect of K on Fe nutrition, which was linked to rhizosphere acidification (excessive K⁺ uptake and consequent release of H⁺ ions by roots to maintain a cation/anion balance). The activity of H+ ATPase pumps located at the plasma membrane is responsible for rhizosphere acidification, which is essential for nutrient acquisition. The liberation of proton is in the favor of uptake of Fe and other macro and micronutrients, especially under the deficiency conditions (Houmani et al., 2015). Hughes et al. (1992) found a significant effect of K⁺ on proton efflux and ferric reduction mechanisms (Strategy I of iron uptaking by roots). They also identified that K⁺ play specific physiological functions in the biosynthetic pathway of mugineic acid production and in the transport of Fe³⁺- phytosiderophore complex (Strategy II). So, we can assume that K⁺ transporter also regulate the iron absorption in plants under iron stress condition.

Zn²⁺ Transporters in Plants (XIPES0157–XIPES0093)

In present study, we have notified one zinc transporter gene of ZIP superfamily that was earlier reported in pearl millet by Mahendrakar et al. (2020). In higher plants, members of the ZIP family are involved in metal uptake, transport, and accumulation of iron and zinc in plant cells. The first identified zinc transporter genes were ZIP1, ZIP2, and ZIP3 in Arabidopsis thaliana (Grotz et al., 1998). Several ZIPs, including OsIRT1, OsIRT2, OsZIP1, OsZIP3, OsZIP4, OsZIP5, OsZIP7, and OsZIP8, have been identified as being responsible for Zn uptake from soil, translocation within the root and from the root to the shoot, and storage in rice seeds (Yang et al., 2009; Lee et al., 2010).

Mg²⁺ Transporters in Plants (IPES0157–IPES0093)

Mg²⁺ transporters are used to maintain Mg in plant organs. Among the proteins potentially involved in Mg²⁺ transport, Magnesium transporter (MGT) family Mg²⁺ transporters plant was the most well-investigated protein (Kobayashi and Tanoi, 2015). Meanwhile, the other transporters of Mg²⁺ are also possible such as Non-selective cation channels (NSCCs) (Guo et al., 2003). According to Demidchik and Maathuis (2007), Voltage-independent NSCC (VI-NSCC) reported the uptake of several cations including Mg²⁺, Ca²⁺, and Zn²⁺ at the resting membrane potentials. The Arabidopsis Magnesium/proton exchanger (MHX) protein is a vacuolar exchanger of protons with cytosolic Mg²⁺ and Zn²⁺. Hence, it may be possible that Mg and Zn share a common pathway of transport in plants.

We have identified one of the gene in ferritin gene family that was already reported in pearl millet via insilico studies and validated by qRT-PCR gene expression (Mahendrakar et al., 2020). The transporter genes like Al, Mg, and K transporter are reporting for the first time in pearl millet. Whereas these transporter gene family were found previously in many different crops such as Al transporters gene in wheat (Li T. et al., 2017) and Rice (Sonoda et al., 2003); Mg transporter gene in Maize (Li et al., 2016, 2018) and Arabidopsis (Li H. et al., 2017); K transporter gene in Maize (Zhang et al., 2012) and Barley (Santa-María et al., 1997); and Zn transporter gene in Arabidopsis (Grotz et al., 1998), rice (Ishimaru et al., 2005; Lee et al., 2010; Neeraja et al., 2018) and Maize (Li et al., 2013).

Conclusion

Enhancing the level of grain Fe and Zn content among elite cultivars of pearl millet become the priority area in answering the nutritional imbalance among poor inhabitance of the semi-arid tropics, especially women and children. Unravelling the genetic principles involved in inheritance of complex traits like grain micronutrient content in pearl millet via, trait-specific mapping enables rapid genetic gains. Huge genetic variation among RIL population observed under the present study can be further used to develop high yielding micronutrient rich cultivar through heterosis breeding. The two co-localized QTLs on LG 2 and LG 3, for both the minerals found to promising as it could simultaneously transferred to the parental lines of elite pearl millet hybrids through marker-assisted breeding programs.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

CS, SPS, and AK designed and supervised the overall research and contributed in the preparation of the manuscript. CS and SPS provided technical guidance for creating liaison among multi-environment locations, and edited the manuscript for final submission. TS, SMS, JB, and MM executed the field experiments. TS and SMS performed the phenotyping, carried out statistical analysis, and prepared the manuscript draft. CB and NA reviewed the manuscript. All authors contributed to the article and approved the submitted version.

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.

Funding

Financial support of ICAR through the project entitled “CRP on biofotification pearl millet [Pennisetum glaucum (L.) Br.]” (ICAR-Sanction number: CS.11/6/2014-IA-IV).

Acknowledgments

The authors are thankful to A. K. Singh, Director, IARI, New Delhi, Division of Genetics, IARI, New Delhi, ICAR-IARI for providing the necessary facilities for smooth conductance of research.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.659789/full#supplementary-material

Supplementary Figure 1 | A to D shown the QTLs for grain Fe content on LG 2,3,5, and 7 at Delhi location.

Supplementary Figure 2 | A and B shown the QTLs for grain Zn content at LG 2 and 3 at Delhi location.

Supplementary Figure 3 | A to D shown the QTLs for grain Fe content on LG 2,3,5, and 7 at Dharwad location.

Supplementary Figure 4 | A and B shown the QTLs for grain Zn content at LG 2 and 3 at Dharwad location.

Supplementary Figure 5 | A to D shown the QTLs for grain Fe content on LG 2,3,5, and 7 at Jodhpur location.

Supplementary Figure 6 | A and B shown the QTLs for grain Zn content at LG 2 and 3 at Jodhpur location.

Footnotes

1. ^ https://blast.ncbi.nlm.nih.gov/
2. ^ http://www.softberry.com/berry.phtml?topic=fgenesh&group=programs&subgroup=gfind

References

Aboul-Maaty, N. A. F., and Oraby, H. A. S. (2019). Extraction of high-quality genomic DNA from different plant orders applying a modified CTAB-based method. Bull. Natl. Res. Cent. 43:25. doi: 10.1186/s42269-019-0066-1

CrossRef Full Text | Google Scholar

Alagarasan, G., Dubey, M., Aswathy, K. S., and Chandel, G. (2017). Genome wide identification of orthologous ZIP genes associated with zinc and Iron translocation in Setaria italica. Front. Plant Sci. 8:775. doi: 10.3389/fpls.2017.00775

PubMed Abstract | CrossRef Full Text | Google Scholar

Allouis, S., Qi, X., Lindup, S., Gale, M. D., and Devos, K. M. (2001). Construction of a BAC library of pearl millet, Pennisetum glaucum. Theor. Appl. Genet. 102, 1200–1205. doi: 10.1007/s001220100559

CrossRef Full Text | Google Scholar

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | CrossRef Full Text | Google Scholar

Anuradha, K., Agarwal, S., Rao, Y. V., Rao, K. V., Viraktamath, B. C., and Sarla, N. (2012). Mapping QTLs and candidate genes for iron and zinc concentrations in unpolished rice of Madhukar× Swarna RILs. Gene 508, 233–240. doi: 10.1016/j.gene.2012.07.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Anuradha, N., Satyavathi, C. T., Bharadwaj, C., Nepolean, T., Sankar, S. M., Singh, S. P., et al. (2017). Deciphering genomic regions for high grain iron and zinc content using association mapping in pearl millet. Front. Plant Sci. 8:412. doi: 10.3389/fpls.2017.00412

PubMed Abstract | CrossRef Full Text | Google Scholar

Anuradha, N., Satyavathi, C. T., Bharadwaj, C., Sankar, M., and Pathy, L. (2018). Association of agronomic traits and micronutrients in pearl millet. Int. J. Chem. Stud. 6, 181–184.

Google Scholar

Athar, T., Khan, M. K., Pandey, A., Yilmaz, F. G., Hamurcu, M., Hakki, E. E., et al. (2020). Biofortification and the involved modern approaches. J. Elem. 25, 717–731.

Google Scholar

Barak, P., and Chen, Y. (1984). The effect of potassium on iron chlorosis in calcareous soils. J. Plant Nutr. 7, 125–133. doi: 10.1080/01904168409363180

CrossRef Full Text | Google Scholar

Bashir, E. M., Ali, A. M., Ali, A. M., Melchinger, A. E., Parzies, H. K., and Haussmann, B. I. (2014). Characterization of Sudanese pearl millet germplasm for agro-morphological traits and grain nutritional values. Plant Genet. Resour. 12:35. doi: 10.1017/s1479262113000233

CrossRef Full Text | Google Scholar

Blanco, A., Simeone, R., Cenci, A., Gadaleta, A., Tanzarella, O. A., Porceddu, E., et al. (2004). Extension of the Messapia x dicoccoides linkage map of Triticum turgidum (L.) Thell. Cell Mol. Biol. Lett. 9, 529–542.

Google Scholar

Briat, J. F., Duc, C., Ravet, K., and Gaymard, F. (2010). Ferritins and iron storage in plants. Biochim. Biophys. Acta 1800, 806–814. doi: 10.1016/j.bbagen.2009.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, D. J. A. (1983). Effect of calcium and aluminum concentrations on the survival of brown trout (Salmo trutta) at low pH. Bull. Environ. Contam. Toxicol. 30, 582–587. doi: 10.1007/bf01610178

PubMed Abstract | CrossRef Full Text | Google Scholar

Causse, M. A., Fulton, T. M., Cho, Y. G., Ahn, S. N., Chunwongse, J., Wu, K., et al. (1994). Saturated molecular map of the rice genome based on an interspecific backcross population. Genetics 138, 1251–1274. doi: 10.1093/genetics/138.4.1251

CrossRef Full Text | Google Scholar

Cloutier, S., Cappadocia, M., and Landry, B. S. (1997). Analysis of RFLP mapping inaccuracy in Brassica napus L. Theor. Appl. Genet. 95, 83–91. doi: 10.1007/s001220050535

CrossRef Full Text | Google Scholar

Colangelo, E. P., and Guerinot, M. L. (2006). Put the metal to the petal: metal uptake and transport throughout plants. Curr. Opin. Plant Biol. 9, 322–330. doi: 10.1016/j.pbi.2006.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Connorton, J. M., Balk, J., and Rodríguez-Celma, J. (2017). Iron homeostasis in plants–a brief overview. Metallomics 9, 813–823. doi: 10.1039/c7mt00136c

PubMed Abstract | CrossRef Full Text | Google Scholar

Crespo-Herrera, L. A., Govindan, V., Stangoulis, J., Hao, Y., and Singh, R. P. (2017). QTL mapping of grain Zn and Fe concentrations in two hexaploid wheat RIL populations with ample transgressive segregation. Front. Plant Sci. 8:1800. doi: 10.3389/fpls.2017.01800

PubMed Abstract | CrossRef Full Text | Google Scholar

Demidchik, V., and Maathuis, F. J. (2007). Physiological roles of nonselective cation channels in plants: from salt stress to signalling and development. New Phytol. 175, 387–404. doi: 10.1111/j.1469-8137.2007.02128.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Devaux, P., Kilian, A., and Kleinhofs, A. (1995). Comparative mapping of the barley genome with male and female recombination-derived, doubled haploid populations. Mol. Gen. Genet. 249, 600–608. doi: 10.1007/bf00418029

PubMed Abstract | CrossRef Full Text | Google Scholar

Devos, K. M., Pittaway, T. S., Reynolds, A., and Gale, M. D. (2000). Comparative mapping reveals a complex relationship between the pearl millet genome and those of foxtail millet and rice. Theor. Appl. Genet. 100, 190–198. doi: 10.1007/s001220050026

CrossRef Full Text | Google Scholar

Dhyan, S., Chhonkar, P. K., and Dwivedi, B. S. (2005). Manual on Soil, Plant and Water Analysis. New Delhi: Westville Publishing House.

Google Scholar

Doerge, R. W. (2002). Mapping and analysis of quantitative trait loci in experimental populations. Nat. Rev. Genet. 3, 43–52. doi: 10.1038/nrg703

PubMed Abstract | CrossRef Full Text | Google Scholar

Drakakaki, G., Marcel, S., Glahn, R. P., Lund, E. K., Pariagh, S., Fischer, R., et al. (2005). Endosperm-specific co-expression of recombinant soybean ferritin and Aspergillus phytase in maize results in significant increases in the levels of bioavailable iron. Plant Mol. Biol. 59, 869–880. doi: 10.1007/s11103-005-1537-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Falconer, D. S. (1996). Introduction to Quantitative Genetics. Noida: Pearson Education India.

Google Scholar

Ferreira, A., Silva, M. F. D., and Cruz, C. D. (2006). Estimating the effects of population size and type on the accuracy of genetic maps. Genet. Mol. Biol. 29, 187–192. doi: 10.1590/s1415-47572006000100033

CrossRef Full Text | Google Scholar

Gangashetty, P. I., Motagi, B. N., Pavan, R., and Roodagi, M. B. (2016). “Breeding crop plants for improved human nutrition through biofortification: progress and prospects,” in Advances in Plant Breeding Strategies: Agronomic, Abiotic and Biotic Stress Traits, eds J. Al-Khayri, S. Jain, and D. Johnson (Cham: Springer), 35–76. doi: 10.1007/978-3-319-22518-0_2

CrossRef Full Text | Google Scholar

Gaoh, B. S. B., Gangashetty, P. I., Mohammed, R., Dzidzienyo, D. K., and Tongoona, P. (2020). Generation mean analysis of pearl millet [Pennisetum glaucum (L.) R. Br.] grain iron and zinc contents and agronomic traits in West Africa. J. Cereal Sci. 96:103066. doi: 10.1016/j.jcs.2020.103066

CrossRef Full Text | Google Scholar

Garcia-Oliveira, A. L., Chander, S., Ortiz, R., Menkir, A., and Gedil, M. (2018). Genetic basis and breeding perspectives of grain iron and zinc enrichment in cereals. Front. Plant Sci. 9:937. doi: 10.3389/fpls.2018.00937

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghandilyan, A., Vreugdenhil, D., and Aarts, M. G. (2006). Progress in the genetic understanding of plant iron and zinc nutrition. Physiol. Plant. 126, 407–417. doi: 10.1111/j.1399-3054.2006.00646.x

CrossRef Full Text | Google Scholar

Govindaraj, M., Rai, K. N., and Shanmugasundaram, P. (2016). Intra-population genetic variance for grain iron and zinc contents and agronomic traits in pearl millet. Crop J. 4, 48–54. doi: 10.1016/j.cj.2015.11.002

CrossRef Full Text | Google Scholar

Govindaraj, M., Rai, K. N., Shanmugasundaram, P., Dwivedi, S. L., Sahrawat, K. L., Muthaiah, A. R., et al. (2013). Combining ability and heterosis for grain iron and zinc densities in pearl millet. Crop Sci. 53, 507–517. doi: 10.2135/cropsci2012.08.0477

CrossRef Full Text | Google Scholar

Govindaraj, M., Selvi, B., Rajarathinam, S., and Sumathi, P. (2011). Genetic variability and heritability of grain yield components and grain mineral concentration in India’s pearl millet (Pennisetum glaucum (L) R. Br.) accessions. Afr. J. Food Agric. Nutr. Dev. 11:4758.

Google Scholar

Graner, A., Jahoor, A., Schondelmaier, J., Siedler, H., Pillen, K., Fischbeck, G., et al. (1991). Construction of an RFLP map of barley. Theor. Appl. Genet. 83, 250–256. doi: 10.1007/bf00226259

PubMed Abstract | CrossRef Full Text | Google Scholar

Green, L. S., and Rogers, E. E. (2004). FRD3 controls iron localization in Arabidopsis. Plant Physiol. 136, 2523–2531. doi: 10.1104/pp.104.045633

PubMed Abstract | CrossRef Full Text | Google Scholar

Grotz, N., Fox, T., Connolly, E., Park, W., Guerinot, M. L., and Eide, D. (1998). Identification of a family of zinc transporter genes from Arabidopsis that respond to zinc deficiency. Proc. Natl. Acad. Sci. U.S.A. 95, 7220–7224. doi: 10.1073/pnas.95.12.7220

PubMed Abstract | CrossRef Full Text | Google Scholar

Grotz, N., and Guerinot, M. L. (2006). Molecular aspects of Cu, Fe and Zn homeostasis in plants. Biochim. Biophys. Acta Mol. Cell Res. 1763, 595–608. doi: 10.1016/j.bbamcr.2006.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, T. R., Zhang, G. P., Lu, W. Y., Wu, H. P., Chen, J. X., Wu, F. B., et al. (2003). Effect of Al on dry matter accumulation and Al and nutrition contents of barleys differing in Al tolerance. Plant Nutr. Fertil. Sci. 9, 324–330.

Google Scholar

Gupta, H. S., Agrawal, P. K., Mahajan, V., Bisht, G. S., Kumar, A., Verma, P., et al. (2009). Quality protein maize for nutritional security: rapid development of short duration hybrids through molecular marker assisted breeding. Curr. Sci. 96, 230–237.

Google Scholar

HarvestPlus (2014). Biofortification Progress Briefs. Available online at: www.HarvestPlus.org (accessed August, 2014)

Google Scholar

Houmani, H., Rabhi, M., Abdelly, C., and Debez, A. (2015). “Implication of rhizosphere acidification in nutrient uptake by plants: cases of potassium (K), phosphorus (P), and iron (Fe),” in Crop Production and Global Environmental Issues, ed. K. Hakeem (Cham: Springer), 103–122. doi: 10.1007/978-3-319-23162-4_4

CrossRef Full Text | Google Scholar

Hughes, D. F., Jolley, V. D., and Brown, J. C. (1992). Roles for potassium in the iron-stress response mechanisms of strategy I and strategy II plants. J. Plant Nutr. 15, 1821–1839. doi: 10.1080/01904169209364442

CrossRef Full Text | Google Scholar

Ishimaru, Y., Suzuki, M., Kobayashi, T., Takahashi, M., Nakanishi, H., Mori, S., et al. (2005). OsZIP4, a novel zinc-regulated zinc transporter in rice. J. Exp. Bot. 56, 3207–3214. doi: 10.1093/jxb/eri317

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, T., Zhou, J., Chen, J., Zhu, L., Zhao, Y., and Huang, Y. (2013). The genetic architecture of zinc and iron content in maize grains as revealed by QTL mapping and meta-analysis. Breed. Sci. 63, 317–324. doi: 10.1270/jsbbs.63.317

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanatti, A., Rai, K. N., Radhika, K., Govindaraj, M., Sahrawat, K. L., Srinivasu, K., et al. (2014). Relationship of grain iron and zinc content with grain yield in pearl millet hybrids. Crop. Improv. 41, 91–96.

Google Scholar

Keltjens, W. G. (1995). Magnesium uptake by Al-stressed maize plants with special emphasis on cation interactions at root exchange sites. Plant Soil 171, 141–146. doi: 10.1007/bf00009577

CrossRef Full Text | Google Scholar

Keltjens, W. G., and Tan, K. (1993). “Interactions between aluminium, magnesium and calcium with different monocotyledonous and dicotyledonous plant species,” in Plant Nutrition–from Genetic Engineering to Field Practice, ed. N. J. Barrow (Dordrecht: Springer), 719–722. doi: 10.1007/978-94-011-1880-4_158

CrossRef Full Text | Google Scholar

Kobayashi, N. I., and Tanoi, K. (2015). Critical issues in the study of magnesium transport systems and magnesium deficiency symptoms in plants. Int. J. Mol. Sci. 16, 23076–23093. doi: 10.3390/ijms160923076

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotla, A., Phuke, R., Hariprasanna, K., Mehtre, S. P., Rathore, A., Gorthy, S., et al. (2019). Identification of QTLs and candidate genes for high grain Fe and Zn concentration in sorghum [Sorghum bicolor (L.) Moench]. J. Cereal Sci. 90, 102850. doi: 10.1016/j.jcs.2019.102850

CrossRef Full Text | Google Scholar

Kumar, S., Hash, C. T., Nepolean, T., Mahendrakar, M. D., Satyavathi, C. T., Singh, G., et al. (2018). Mapping grain iron and zinc content quantitative trait loci in an iniadi-derived immortal population of pearl millet. Genes 9:248. doi: 10.3390/genes9050248

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Hash, C. T., Thirunavukkarasu, N., Singh, G., Rajaram, V., Rathore, A., et al. (2016). Mapping quantitative trait loci controlling high iron and zinc content in self and open pollinated grains of pearl millet [Pennisetum glaucum (L.) R. Br.]. Front. Plant Sci. 7:1636. doi: 10.3389/fpls.2016.01636

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, T. R., and Prasad, M. N. V. (1999). Ferritin induction by iron mediated oxidative stress and ABA in Vigna mungo (L.) Hepper seedlings: role of antioxidants and free radical scavengers. J. Plant Physiol. 155, 652–655. doi: 10.1016/s0176-1617(99)80068-6

CrossRef Full Text | Google Scholar

Lander, E. S., Green, P., Abrahamson, J., Barlow, A., Daly, M. J., Lincoln, S. E., et al. (1987). MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics 1, 174–181. doi: 10.1016/0888-7543(87)90010-3

CrossRef Full Text | Google Scholar

Lee, S., Jeong, H. J., Kim, S. A., Lee, J., Guerinot, M. L., and An, G. (2010). OsZIP5 is a plasma membrane zinc transporter in rice. Plant Mol. Biol. 73, 507–517. doi: 10.1007/s11103-010-9637-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Du, H., Huang, K., Chen, X., Liu, T., Gao, S., et al. (2016). Identification, and functional and expression analyses of the CorA/MRS2/MGT-type magnesium transporter family in maize. Plant Cell Physiol. 57, 1153–1168. doi: 10.1093/pcp/pcw064

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Liu, C., Zhou, L., Zhao, Z., Li, Y., Qu, M., et al. (2018). Molecular and functional characterization of the magnesium transporter gene ZmMGT12 in maize. Gene 665, 167–173. doi: 10.1016/j.gene.2018.04.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Wang, N., Ding, J., Liu, C., Du, H., Huang, K., et al. (2017). The maize CorA/MRS2/MGT-type Mg transporter, ZmMGT10, responses to magnesium deficiency and confers low magnesium tolerance in transgenic Arabidopsis. Plant Mol. Biol. 95, 269–278. doi: 10.1007/s11103-017-0645-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Wang, J., and Zhang, L. (2015). Inclusive composite interval mapping of QTL by environment interactions in biparental populations. PLoS One 10:e0132414. doi: 10.1371/journal.pone.0132414

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Zhou, X., Huang, Y., Zhu, L., Zhang, S., Zhao, Y., et al. (2013). Identification and characterization of the zinc-regulated transporters, iron-regulated transporter-like protein (ZIP) gene family in maize. BMC Plant Biol. 13:114. doi: 10.1186/1471-2229-13-114

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Liao, K., Xu, X., Gao, Y., Wang, Z., Zhu, X., et al. (2017). Wheat ammonium transporter (AMT) gene family: diversity and possible role in host–pathogen interaction with stem rust. Front. Plant Sci. 8:1637. doi: 10.3389/fpls.2017.01637

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C. J., Witcombe, J. R., Pittaway, T. S., Nash, M., Hash, C. T., Busso, C. S., et al. (1994). An RFLP-based genetic map of pearl millet (Pennisetum glaucum). Theor. Appl. Genet. 89, 481–487. doi: 10.1007/bf00225384

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D. J., Wang, Y. B., Guo, C. H., Cong, Q., Gong, X., and Zhang, H. J. (2016). Enhanced iron and zinc accumulation in genetically engineered wheat plants using sickle alfalfa (Medicago falcata L.) ferritin gene. Cereal Res. Commun. 44, 24–34. doi: 10.1556/0806.43.2015.039

CrossRef Full Text | Google Scholar

Livingstone, K. D., Lackney, V. K., Blauth, J. R., Van Wijk, R. I. K., and Jahn, M. K. (1999). Genome mapping in Capsicum and the evolution of genome structure in the Solanaceae. Genetics 152, 1183–1202.

Google Scholar

Lu, H., Romero-Severson, J., and Bernardo, R. (2002). Chromosomal regions associated with segregation distortion in maize. Theor. Appl. Genet. 105, 622–628. doi: 10.1007/s00122-002-0970-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, K., Li, L., Zheng, X., Zhang, Z., Mou, T., and Hu, Z. (2008). Quantitative trait loci controlling Cu, Ca, Zn, Mn and Fe content in rice grains. J. Genet. 87:305. doi: 10.1007/s12041-008-0049-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Luciano Da Costa, E. S., Wang, S., and Zeng, Z. B. (2012). Composite interval mapping and multiple interval mapping: procedures and guidelines for using Windows QTL Cartographer. Methods Mol. Biol. 871, 75–119. doi: 10.1007/978-1-61779-785-9_6

CrossRef Full Text | Google Scholar

Lukaszewski, K. M., and Blevins, D. G. (1996). Root growth inhibition in boron-deficient or aluminum-stressed squash may be a result of impaired ascorbate metabolism. Plant Physiol. 112, 1135–1140. doi: 10.1104/pp.112.3.1135

PubMed Abstract | CrossRef Full Text | Google Scholar

Magallanes-López, A. M., Hernandez-Espinosa, N., Velu, G., Posadas-Romano, G., Ordoñez-Villegas, V. M. G., Crossa, J., et al. (2017). Variability in iron, zinc and phytic acid content in a worldwide collection of commercial durum wheat cultivars and the effect of reduced irrigation on these traits. Food Chem. 237, 499–505. doi: 10.1016/j.foodchem.2017.05.110

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahendrakar, M. D., Parveda, M., Kishor, P. K., and Srivastava, R. K. (2020). Discovery and validation of candidate genes for grain iron and zinc metabolism in pearl millet [Pennisetum glaucum (L.) R. Br.]. Sci. Rep. 10:16562.

Google Scholar

Mohan, M., Nair, S., Bhagwat, A., Krishna, T. G., Yano, M., Bhatia, C. R., et al. (1997). Genome mapping, molecular markers and marker-assisted selection in crop plants. Mol. Breed. 3, 87–103.

Google Scholar

Muthusamy, V., Hossain, F., Thirunavukkarasu, N., Choudhary, M., Saha, S., Bhat, J. S., et al. (2014). Development of β-carotene rich maize hybrids through marker-assisted introgression of β-carotene hydroxylase allele. PLoS One 9:e113583. doi: 10.1371/journal.pone.0113583

PubMed Abstract | CrossRef Full Text | Google Scholar

Neeraja, C. N., Kulkarni, K. S., Madhu Babu, P., Sanjeeva Rao, D., Surekha, K., and Ravindra Babu, V. (2018). Transporter genes identified in landraces associated with high zinc in polished rice through panicle transcriptome for biofortification. PLoS One 13:e0192362. doi: 10.1371/journal.pone.0192362

PubMed Abstract | CrossRef Full Text | Google Scholar

Nozoye, T., Nagasaka, S., Kobayashi, T., Takahashi, M., Sato, Y., Sato, Y., et al. (2011). Phytosiderophore efflux transporters are crucial for iron acquisition in graminaceous plants. J. Biol. Chem. 286, 5446–5454. doi: 10.1074/jbc.m110.180026

PubMed Abstract | CrossRef Full Text | Google Scholar

Oertli, J. J., and Opoku, A. A. (1974). Interaction of potassium in the availability and uptake of iron from ferric hydroxide. Soil Sci. Soc. Am. J. 38, 451–454. doi: 10.2136/sssaj1974.03615995003800030023x

CrossRef Full Text | Google Scholar

Patterson, H. D., and Williams, E. R. (1976). A new class of resolvable incomplete block designs. Biometrika 63, 83–92. doi: 10.1093/biomet/63.1.83

CrossRef Full Text | Google Scholar

Pereira, J. F., Zhou, G., Delhaize, E., Richardson, T., Zhou, M., and Ryan, P. R. (2010). Engineering greater aluminium resistance in wheat by over-expressing TaALMT1. Ann. Bot. 106, 205–214. doi: 10.1093/aob/mcq058

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfeiffer, W. H., and McClafferty, B. (2007). HarvestPlus: breeding crops for better nutrition. Crop Sci. 47, S-88–S-105.

Google Scholar

Qu, L., Yoshihara, T., Ooyama, A., Goto, F., and Takaiwa, F. (2005). Iron accumulation does not parallel the high expression level of ferritin in transgenic rice seeds. Planta 222, 225–233. doi: 10.1007/s00425-005-1530-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Quarrie, S. A., Steed, A., Calestani, C., Semikhodskii, A., Lebreton, C., Chinoy, C., et al. (2005). A high-density genetic map of hexaploid wheat (Triticum aestivum L.) from the cross Chinese Spring× SQ1 and its use to compare QTLs for grain yield across a range of environments. Theor. Appl. Genet. 110, 865–880. doi: 10.1007/s00122-004-1902-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Rai, K. N., Govindaraj, M., and Rao, A. S. (2012). Genetic enhancement of grain iron and zinc content in pearl millet. Qual. Assur. Saf. Crop. Foods 4, 119–125. doi: 10.1111/j.1757-837x.2012.00135.x

CrossRef Full Text | Google Scholar

Rai, K. N., Velu, G., Govindaraj, M., Upadhyaya, H. D., Rao, A. S., Shivade, H., et al. (2015). Iniadi pearl millet germplasm as a valuable genetic resource for high grain iron and zinc densities. Plant Genet. Resour. 13, 1–8.

Google Scholar

Rajaram, V., Nepolean, T., Senthilvel, S., Varshney, R. K., Vadez, V., Srivastava, R. K., et al. (2013). Pearl millet [Pennisetum glaucum (L.) R. Br.] consensus linkage map constructed using four RIL mapping populations and newly developed EST-SSRs. BMC Genomics 14:159. doi: 10.1186/1471-2164-14-159

PubMed Abstract | CrossRef Full Text | Google Scholar

Raman, H., Zhang, K., Cakir, M., Appels, R., Garvin, D. F., Maron, L. G., et al. (2005). Molecular characterization and mapping of ALMT1, the aluminium-tolerance gene of bread wheat (Triticum aestivum L.). Genome 48, 781–791. doi: 10.1139/g05-054

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts, L. A., Pierson, A. J., Panaviene, Z., and Walker, E. L. (2004). Yellow stripe1. Expanded roles for the maize iron-phytosiderophore transporter. Plant Physiol. 135, 112–120. doi: 10.1104/pp.103.037572

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryan, P. R., Tyerman, S. D., Sasaki, T., Furuichi, T., Yamamoto, Y., Zhang, W. H., et al. (2011). The identification of aluminium-resistance genes provides opportunities for enhancing crop production on acid soils. J. Exp. Bot. 62, 9–20. doi: 10.1093/jxb/erq272

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabiel, S. A., Ismail, M. I., Abdalla, E., Osman, K. A., and Ali, A. M. (2014). Genetic variation among pearl millet genotypes for yield and its components in semi-arid zone Sudan. IJACS 7, 822–826.

Google Scholar

Sakaguchi, T., Nishizawa, N. K., Nakanishi, H., Yoshimura, E., and Mori, S. (1999). The role of potassium in the secretion of mugineic acids family phytosiderophores from iron-deficient barley roots. Plant Soil 215, 221–227.

Google Scholar

Salamov, A. A., and Solovyev, V. V. (2000). Ab initio gene finding in Drosophila genomic DNA. Genome Res. 10, 516–522. doi: 10.1101/gr.10.4.516

PubMed Abstract | CrossRef Full Text | Google Scholar

Sankar, S. M., Satyavathi, C. T., Singh, M. P., Bharadwaj, C., Singh, S. P., and Barthakur, S. (2013). Genetic variability and association studies in pearl millet for grain yield and high temperature stress tolerance. Indian J. Dryland Agric. Res. Dev. 28, 71–76.

Google Scholar

Santa-María, G. E., Rubio, F., Dubcovsky, J., and Rodríguez-Navarro, A. (1997). The HAK1 gene of barley is a member of a large gene family and encodes a high-affinity potassium transporter. Plant Cell 9, 2281–2289. doi: 10.2307/3870585

CrossRef Full Text | Google Scholar

Satyavathi, C. T., Sankar, S. M., Singh, S. P., Bhowmick, P., Bhat, J., Singh, O., et al. (2015). Stability analysis of grain iron and zinc content in pearl millet (Pennisetum glaucum (L.) R. Br). J. Trop. Agric. 33, 1387–1394.

Google Scholar

Senthilvel, S., Jayashree, B., Mahalakshmi, V., Kumar, P. S., Nakka, S., Nepolean, T., et al. (2008). Development and mapping of simple sequence repeat markers for pearl millet from data mining of expressed sequence tags. BMC Plant Biol. 8:119. doi: 10.1186/1471-2229-8-119

PubMed Abstract | CrossRef Full Text | Google Scholar

Shariatipour, N., and Heidari, B. (2020). “Genetic-based biofortification of staple food crops to meet zinc and iron deficiency-related challenges,” in Plant Micronutrients, eds T. Aftab and K. R. Hakeem (Cham: Springer), 173–223. doi: 10.1007/978-3-030-49856-6_8

CrossRef Full Text | Google Scholar

Shivhare, R., and Lata, C. (2016). Selection of suitable reference genes for assessing gene expression in pearl millet under different abiotic stresses and their combinations. Sci. Rep. 6:23036.

Google Scholar

Singh, V. K., Singh, A., Singh, S. P., Ellur, R. K., Choudhary, V., Sarkel, S., et al. (2012). Incorporation of blast resistance into “PRR78”, an elite Basmati rice restorer line, through marker assisted backcross breeding. Field Crops Res. 128, 8–16. doi: 10.1016/j.fcr.2011.12.003

CrossRef Full Text | Google Scholar

Ślaski, J. J., Zhang, G., Basu, U., Stephens, J. L., and Taylor, G. J. (1996). Aluminum resistance in wheat (Triticum aestivum) is associated with rapid, Al-induced changes in activities of glucose-6-phosphate dehydrogenase and 6-phosphogluconate dehydrogenase in root apices. Physiol. Plant. 98, 477–484. doi: 10.1034/j.1399-3054.1996.980307.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sonoda, Y., Ikeda, A., Saiki, S., Yamaya, T., and Yamaguchi, J. (2003). Feedback regulation of the ammonium transporter gene family AMT1 by glutamine in rice. Plant Cell Physiol 44, 1396–1402. doi: 10.1093/pcp/pcg169

PubMed Abstract | CrossRef Full Text | Google Scholar

Sumathi, P., Madineni, S., and Veerabadhiran, P. (2010). Genetic variability for different biometrical traits in pearl millet genotypes (Pennisetum glaucum LR BR.). Electron. J. Plant Breed. 1, 437–440.

Google Scholar

Supriya, A., Senthilvel, S., Nepolean, T., Eshwar, K., Rajaram, V., Shaw, R., et al. (2011). Development of a molecular linkage map of pearl millet integrating DArT and SSR markers. Theor. Appl. Genet. 123, 239–250. doi: 10.1007/s00122-011-1580-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Varshney, R. K., Hoisington, D. A., and Tyagi, A. K. (2006). Advances in cereal genomics and applications in crop breeding. Trends Biotechnol. 24, 490–499. doi: 10.1016/j.tibtech.2006.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasconcelos, M., Datta, K., Oliva, N., Khalekuzzaman, M., Torrizo, L., Krishnan, S., et al. (2003). Enhanced iron and zinc accumulation in transgenic rice with the ferritin gene. Plant Sci. 164, 371–378. doi: 10.1016/s0168-9452(02)00421-1

CrossRef Full Text | Google Scholar

Velu, G., Tutus, Y., Gomez-Becerra, H. F., Hao, Y., Demir, L., Kara, R., et al. (2017). QTL mapping for grain zinc and iron concentrations and zinc efficiency in a tetraploid and hexaploid wheat mapping populations. Plant Soil. 411, 81–99. doi: 10.1007/s11104-016-3025-8

CrossRef Full Text | Google Scholar

Vinodhana, N. K., Sumathi, P., and Sathya, M. (2013). Genetic variability and inter-relationship among morpho-economic traits of pearl millet (Pennisetum glaucum (L.) r. Br.) and their implications in selection. IJPAES 3, 145–149.

Google Scholar

Voorrips, R. E. (2002). MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77–78. doi: 10.1093/jhered/93.1.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Vreugdenhil, D., Aarts, M. G. M., Koornneef, M., Nelissen, H., and Ernst, W. H. O. (2004). Natural variation and QTL analysis for cationic mineral content in seeds of Arabidopsis thaliana. Plant Cell Environ. 27, 828–839. doi: 10.1111/j.1365-3040.2004.01189.x

CrossRef Full Text | Google Scholar

Wang, X., Yang, S., Chen, Y., Zhang, S., Zhao, Q., Li, M., et al. (2018). Comparative genome-wide characterization leading to simple sequence repeat marker development for Nicotiana. BMC Genomics 19:500. doi: 10.1186/s12864-018-4878-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, M., Harada, E., Vess, C., Roepenack-Lahaye, E. V., and Clemens, S. (2004). Comparative microarray analysis of Arabidopsis thaliana and Arabidopsis halleri roots identifies nicotianamine synthase, a ZIP transporter and other genes as potential metal hyperaccumulation factors. Plant J. 37, 269–281. doi: 10.1046/j.1365-313x.2003.01960.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wendel, J. F., Edwards, M. D., and Stuber, C. W. (1987). Evidence for multilocus genetic control of preferential fertilisation in maize. Heredity 58, 297–301. doi: 10.1038/hdy.1987.44

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Zhu, L., Xiao, J., Huang, N., and McCouch, S. R. (1997). Chromosomal regions associated with segregation distortion of molecular markers in F 2, backcross, doubled haploid, and recombinant inbred populations in rice (Oryza sativa L.). Mol. Gen. Genet. 253, 535–545. doi: 10.1007/s004380050355

PubMed Abstract | CrossRef Full Text | Google Scholar

Yadav, R. S., Hash, C. T., Bidinger, F. R., Devos, K. M., and Howarth, C. J. (2004). Genomic regions associated with grain yield and aspects of post-flowering drought tolerance in pearl millet across stress environments and tester background. Euphytica 136, 265–277. doi: 10.1023/b:euph.0000032711.34599.3a

CrossRef Full Text | Google Scholar

Yadav, R. S., Sehgal, D., and Vadez, V. (2011). Using genetic mapping and genomics approaches in understanding and improving drought tolerance in pearl millet. J. Exp. Bot. 62, 397–408. doi: 10.1093/jxb/erq265

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Y. W., Mao, D. D., Yang, L., Qi, J. L., Zhang, X. X., Tang, Q. L., et al. (2018). Magnesium transporter MGT6 plays an essential role in maintaining magnesium homeostasis and regulating high magnesium tolerance in Arabidopsis. Front. Plant Sci. 9:274. doi: 10.3389/fpls.2018.00274

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, M., Li, Y., Liu, Z., Tian, J., Liang, L., Qiu, Y., et al. (2020). A high activity zinc transporter OsZIP9 mediates zinc uptake in rice. Plant J. 103, 1695–1709. doi: 10.1111/tpj.14855

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Wu, Y., Li, Y., Ling, H. Q., and Chu, C. (2009). OsMT1a, a type 1 metallothionein, plays the pivotal role in zinc homeostasis and drought tolerance in rice. Plant Mol. Biol. 70, 219–229. doi: 10.1007/s11103-009-9466-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Yokosho, K., Yamaji, N., and Ma, J. F. (2016). OsFRDL1 expressed in nodes is required for distribution of iron to grains in rice. J. Exp. Bot. 67, 5485–5494. doi: 10.1093/jxb/erw314

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Zhang, J., Chen, Y., Li, R., Wang, H., and Wei, J. (2012). Genome-wide analysis and identification of HAK potassium transporter gene family in maize (Zeamays L.). Mol. Biol. Rep. 39, 8465–8473. doi: 10.1007/s11033-012-1700-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zielińska-Dawidziak, M. (2015). Plant ferritin–a source of iron to prevent its deficiency. Nutrients 7, 1184–1201. doi: 10.3390/nu7021184

PubMed Abstract | CrossRef Full Text | Google Scholar