Genetic Correlation Between Fe and Zn Biofortification and Yield Components in a Common Bean (Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) is the most important legume for direct human consumption worldwide. It is a rich and relatively inexpensive source of proteins and micronutrients, especially iron and zinc. Bean is a target for biofortification to develop new cultivars with high Fe/Zn levels that help to ameliorate malnutrition mainly in developing countries. A strong negative phenotypic correlation between Fe/Zn concentration and yield is usually reported, posing a significant challenge for breeders. The objective of this study was to investigate the genetic relationship between Fe/Zn. We used Quantitative Trait Loci (QTLs) mapping and Genome-Wide Association Studies (GWAS) analysis in three bi-parental populations that included biofortified parents, identifying genomic regions associated with yield and micromineral accumulation. Significant negative correlations were observed between agronomic traits (pod harvest index, PHI; pod number, PdN; seed number, SdN; 100 seed weight, 100SdW; and seed per pod, Sd/Pd) and micronutrient concentration traits (SdFe and SdZn), especially between pod harvest index (PHI) and SdFe and SdZn. PHI presented a higher correlation with SdN than PdN. Seventy-nine QTLs were identified for the three populations: 14 for SdFe, 12 for SdZn, 13 for PHI, 11 for SdN, 14 for PdN, 6 for 100SdW, and 9 for Sd/Pd. Twenty-three hotspot regions were identified in which several QTLs were co-located, of which 13 hotpots displayed QTL of opposite effect for yield components and Fe/Zn accumulation. In contrast, eight QTLs for SdFe and six QTLs for SdZn were observed that segregated independently of QTL of yield components. The selection of these QTLs will enable enhanced levels of Fe/Zn and will not affect the yield performance of new cultivars focused on biofortification.

Introduction

Common bean (Phaseolus vulgaris L.) is one of the most important legumes for human consumption and is cultivated worldwide, particularly in tropical and subtropical countries of Africa and America (Broughton et al., 2003). Cultivated common bean is organized in two major gene pools, Mesoamerican and Andean, which resulted from independent domestication events, generating landraces with diverse seed quality, nutritional traits, resistance to diseases, tolerance to abiotic stress, morphologic, and agronomic traits (Schmutz et al., 2014). Bean has been documented as a highly valuable food for the human diet, an accessible and relatively inexpensive source of proteins and micronutrients, such as iron and zinc (Beebe, 2012).

Micronutrient deficiencies are among the most common global nutritional problems, affecting half of the population of the world, especially women and children in developing countries (Bouis and Welch, 2010). Fe and Zn deficiencies are among the main common causes of micromineral malnutrition worldwide. In humans, Fe is essential for several cellular mechanisms (Steinbicker and Muckenthaler, 2013). Iron Deficiency Anemia (IDA), a pathology associated with scarce iron intake and absorption, affects a third of the population of the world and has been associated with several chronic diseases (Lopez et al., 2016). Zn has a key role in regulating growth and in the functioning of the immune system (Chasapis et al., 2012), and its deficiency can contribute to infections (Hambidge, 2000).

Wide phenotypic variability has been reported in bean in terms of Fe and Zn accumulation, with seed Fe concentration from 35 to 90 mg/kg and seed Zn concentration from 20 to 54 mg/kg (Beebe et al., 2000), these ranges are higher than those of other main staple crops, such as rice, wheat, and maize (Petry et al., 2015). For this reason, common bean is one of the crops targeted for biofortification. A cup of beans supplies 25% of Fe and 15% of Zn to the recommended daily allowance in the human diet, but breeding programs have the potential to develop genotypes with two to three times more Fe and Zn concentrations (Cichy et al., 2009).

Bouis (2003) reported that biofortification breeding has multiple positive effects, such as more efficiency in nutrients, fertilizer, and irrigation demands, and can increase disease resistance and yield. However, some difficulties have been described in bean breeding for micronutrient concentration, mainly a strong negative correlation between Fe and Zn concentration and yield. Source genotypes with high iron contents (>100 mg/kg) have been identified, but most of these lines have poor agronomic performance (Raatz, 2018). For this reason, a major breeding objective is overcoming the negative correlation between the micronutrient concentrations with yield components and finding the optimum balance between these traits.

While the phenotypic variability for Fe and Zn in common bean has been utilized in breeding programs with some success, the development of genomic tools will likely be helpful to reach the goals of biofortification (Beaver and Osorno, 2009). Currently, the genetic architecture of Fe and Zn accumulation in seed in common bean has been the object of several studies (Izquierdo et al., 2018), usually using Quantitative Trait Loci (QTL) analysis to find genomic regions related to the genetic control of these traits. Nevertheless, molecular markers have not been implemented in bean breeding programs focused on mineral concentration, as parental materials employed for populations are often distant from elite breeding materials, and the QTLs found with major effects are scarce and cannot be extrapolated to other populations (Raatz, 2018).

Another factor to consider in biofortification breeding is the accompanying agronomic traits, such as drought tolerance. Hummel et al. (2018) predicted a large part of current common bean growing areas to become unsuitable for cultivation by 2050, generating reductions in yield for current varieties and potentially affecting the nutritional quality, especially the iron levels of the crop, if current climate change trends continue. Drought is the main cause of yield loss worldwide, and its incidence and duration are expected to increase due to climate change (Lauer et al., 2012). Drought stress is a principal abiotic constraint that affects beans, causing losses close to 60% in production (Beebe et al., 2008).

The objective of this study was to investigate the correlation of genetic architecture of Fe and Zn accumulation and yield components in a common bean under irrigated and drought conditions using QTL and Genome-Wide Association Studies (GWAS) analysis in three biofortified bi-parental Mesoamerican populations and to identify genomic regions associated with yield components, micromineral accumulation, and their possible relationship.

Materials and Methods

Plant Material

In this study we used three Mesoamerican bi-parental populations of common bean (1) SCR16 × SMC40, (2) SMC44 × SCR9, and (3) SMC33 × SCR16. Parents SCR16 and SCR9 are improved Mesoamerican lines from the International Center of Agricultural Tropical (CIAT) with traits related to drought tolerance and resistance to bean common mosaic virus (BCMV) and bean golden yellow mosaic virus (BGYMV). The parents SMC40, SMC44, and SMC33 are improved Mesoamerican lines with a high accumulation of Fe and Zn in the seed. Each population consists of 100 recombinant inbred lines (RILs). Populations were advanced by single seed descent, and families from which individual plants were drawn were evaluated for iron and zinc concentration, and yield components in the F3, F4, F5, F6, F7, and F8 generations (Table 1).

TABLE 1

Table 1. Description of trials used in the experiment.

Field Trials Design and Phenotyping

The three bi-parental populations were planted at the CIAT experimental station in Palmira, Colombia (altitude of 1,000 m.a.s.l., latitude 3° 32′ N and longitude 76° 18′ W). The experimental design used for this study was an augmented design with three blocks in which five different check lines (SER118, SER16, Tio Canela 75, SEN56, and DOR390), and each parent for each population was evenly distributed across the field. The experimental units in these trials were row plots of 2.22 m² for each line. The experiments were conducted in 5 years (2011, 2012, 2013, 2014, and 2015) in four irrigated trials and four drought trials, for a total of eight trials, with data taken on families as the RILs were advanced from F2.3 to F6.8 (Table 1). For more information about the climatic conditions see Supplementary Figure 1.

The seed was harvested manually by plot upon maturity (75–80 days after sowing). This was estimated when 50% of plants of the plot had at least one pod losing its green pigmentation. At the time of the harvest, 0.5 m per plot were collected separately to measure the number of pods per m² (pod number, PdN), number of seeds per m² (seed number, SdN), 100 seed weight (100SdW) was obtained estimating the weight of the seeds per plot and corrected to the moisture content in the seed of 14%, number of seeds per pod (Sd/Pd), and pod harvest index (PHI, %) defined as the ratio between seed weight and total pod weight, such as seed (Assefa et al., 2013). In 2011, only PHI was measured. Additional information on the description of phenotypic traits can be found in ‘‘Traits Dictionaries for Fieldbook Development.’’¹.

Fe and Zn Quantification

The samples to evaluate iron and zinc concentrations in the seed (SdFe and SdZn, mg/kg) were prepared according to the method described by Stangoulis and Sison (2008). In brief, we collected 5 g of seed and cleaned using a cloth dampened with high-purity water. Later, seeds were dried to 7–8% moisture content in an oven at 60°C. Finally, samples were ground with a grinding mill Retsch Mixer Mill MM 400 using grinding jars and zirconium balls (Retsch GmbH & Co KG, Germany) to avoid iron and zinc contamination. To quantify SdFe and SdZn a Energy Dispersive X-ray Fluorescence machine (EDXRF) Bruker S2 PUMA (Bruker Corporation, Billerica, MA, United States) was used. The quantification was developed according to the method described by Guild et al. (2017a). In summary, the EDXRF analysis conditions were as follows. Atmosphere: air, X-ray tube: palladium (50 W), voltage: 40 kV, current: 240 μA, peak detected: Kα, acquisition time: 30 s, tube filter: Al 500 nm, and detector: silicon drift detector. Analyses were conducted in supplied sample cups prepared as reported previously by Guild et al. (2017a). Sample cups were cleaned and prepared before each analysis to minimize cross-contamination between samples. For the calibration, we used 10 custom-made 40 mm diameter glass disks (FLUXANA GmbH & Co., Germany) with a range of nominal Fe and Zn levels to establish a non-matrix-matched calibration. SdFe and SdZn were measured in all trials except in the 2011.

Data Analysis

Statistical analysis of the phenotypic data was conducted with statistical software R (v3.3.2.). Best linear unbiased estimators and predictors (BLUEs/BLUPs) were calculated for each trait and each trial using the “lme4” Package (Bates et al., 2015). The data were modeled using the following equation:

${\mathit{\text{y}}}_{\mathit{mij}}=\mathit{\mu }+{\mathit{\text{G}}}_{\mathit{m}}+{\mathit{\text{R}}}_{\mathit{i}}+{\mathit{\epsilon }}_{\mathit{mij}}$

where y is a vector with the phenotypic data, μ is the overall intercept, G_m is the effect of the m-th genotype, R_i is the effect of the i-th replication, and ε_mij is the error term corresponding to y_mij. In this model, G_m term effects were treated either as fixed (to calculated BLUEs) or random (to obtain an estimate of the genetic variance and calculated BLUPs). We assumed that every random term μ and the residual ε adjust to a normal distribution with mean 0 and independent variances μ∼ N (0, σ²_uI) and ε∼ N (0, σ²_ε I).

To determine the proportion of the genetic variance controlling for evaluated traits, broad-sense heritability (H²) estimates were calculated using the method proposed by Cullis et al. (2006), using the equation below:

${\mathit{\text{H}}}^{\mathit{2}}=\mathbf{1}-\frac{{\mathit{\upsilon }}_{\mathit{BLUP}}}{\mathbf{2}{\mathit{\sigma }}_{\mathit{genotype}}^{\mathit{2}}}$

where υ_BLUP is the mean-variance of a difference of two BLUPs of genotypic effects, and $2{\mathrm{\sigma }}_{genotype}^2$ is the genetic variance. The phenotypic correlation between traits measured was expressed as Pearson’s correlation coefficients among BLUEs for individual trials and cross-trials means, and their significance was tested using a two-tailed t-test.

Genotyping

DNA extraction of the RILs and parents was developed on the F6.8 generation. For DNA extraction, seeds of each genotype were selected and germinated using a Calcium Sulfate 0.5-mM solution, the DNA was extracted from young leaves using liquid N² and the urea buffer-based extraction miniprep protocol. The genotyping was performed using the ApeKI-based genotyping-by-sequencing (GBS) protocol (Elshire et al., 2011), and the GBS libraries were sequenced at the HudsonAlpha Institute for Biotechnology.² Some lines were re-genotyped where seed quality was poor.

The mapping and variant calling processes for GBS reads are described in detail by Perea et al. (2016). In brief, the GBS reads were demultiplexed using the Next Generation Sequencing Eclipse Plugin (NGSEP; v3.3.2) (Tello et al., 2019). Adapters and low-quality bases from the raw sequencing data were trimmed using Trimmomatic (v0.36) (Bolger et al., 2014), and reads were aligned to the reference genome of P. vulgaris G19833 v2.1. (Schmutz et al., 2014) using Bowtie2 (v2.2.30) (Langmead and Salzberg, 2012) with default parameters. The variant calling process of RILs was performed using NGSEP following recommended parameters for GBS data (Perea et al., 2016), using the genotype calls of the parents (SMC44, SMC40, SMC33, SCR16, and SCR9) as reference. The resulting genotype matrix was filtered to remove genotype calls with quality below 40, remove markers with more than 6% heterozygote calls, minor allele frequency (MAF) below 0.02, and remove markers in repetitive regions of the reference genome (Lobaton et al., 2018). The process to construct a genotype matrix was carried out for each population separately and grouping the populations for GWAS analysis. The final matrices contained 1,650 markers in the SCR16 × SMC40 population, 3,033 markers in the SMC44 × SCR9 population, and 481 markers in the SMC33 × SCR16 population.

Genetic Map Construction, Quantitative Trait Loci Analysis, and Genome-Wide Association Studies Analysis

For the construction of the genetic maps, the integrated tools of the IciMapping Software v4.1 (Meng et al., 2015) were used, employing the matrices above described. The genotypic data were grouped and anchored on their physical position. The default logarithm of odds (LOD) score (3.0) and rfs (0.35) were used as a maximum threshold value for grouping markers. The Kosambi genetic mapping function was chosen for genetic interval estimation using recombination frequency. The orders of markers per chromosome were set by their physical position, and distortion of segregated markers was exhibited according to segregation ratios. QTL analysis for the populations was conducted using the genetic maps of each population and the phenotypic BLUEs for each RIL. Detection of QTLs and estimation of genetic parameters for SdFe, SdZn, PHI, PdN, SdN, 100SdW, and Sd/Pd were performed using the inclusive composite interval mapping for additive effects (ICIM-ADD) procedure of the software IciMapping (v4.1) (Meng et al., 2015) with 10 cm windows and a sliding parameter of 1 cm. Significant QTL was considered by defining the LOD score at a genome-wide type I error rate of α = 0.05 after 1,000 permutation tests for each trait. The detected QTLs having high phenotypic variation explained (PVE) (>10%) were justified as major QTL (Amanullah et al., 2021). QTL hotspots were defined as significant loci genetically linked at <5 cm between QTL.

Genome-Wide Association Studies analysis was carried out using the R package GAPIT (v3.0) (Tang et al., 2016). The association analysis was conducted using a Compression Mixed Linear Model (CMLM) approach. This model accounts for population structure using the top four principal components as fixed effects. It also accounts for random polygenic effects with a kinship matrix as variance-covariance structure, calculated using the method proposed by VanRaden (2008) implemented in the GAPIT package. Significant associations were defined when the p-value was equal to or smaller than the Bonferroni-corrected threshold (1.79 × 10^–6) (Johnson et al., 2010). Manhattan and quantile–quantile (Q-Q) plot graphics were made using the qqman R package (Turner, 2018).

Results

Phenotyping Variation for Traits Evaluated

Five elite lines of the CIAT bean-breeding program were used to develop three Mesoamerican bi-parental populations: (1) SCR16 × SMC40, (2) SMC33 × SCR16, and (3) SMC44 × SCR9. Families in several generations and F6.8 RILs were evaluated for a range of yield components and micronutrient concentrations of iron and zinc (SdFe and SdZn) in trials managed under irrigated and drought stress conditions. The phenotypic data show high and significant phenotypic variability among RIL in all traits evaluated (Supplementary Tables 1, 2 and Supplementary Figure 2). Transgressive segregation was also evident in the three populations with lower and higher values compared to parents in all traits. An effect of drought stress was not observed on SdFe, SdZn, PHI, and Sd/Pd in the populations evaluated. However, drought reduced PdN, SdN, and 100SdW by 29.5, 39.6, and 31.4%, respectively in all populations. The SdZn performance of SMC44 × SCR9 was better in most trials compared to the other populations (27.2 mg/kg on average), while SCR16 × SMC40 presented higher SdFe than the other two populations (78.3 mg/kg on average) (Figures 1, 2 and Supplementary Table 1).

FIGURE 1

Figure 1. Phenotypic distribution of micronutrients content traits. Fe accumulation (SdFe) (A) and Zn accumulation (SdZn) (B) evaluated of 2012, 2013, 2014, and 2015 under drought and irrigated conditions for the SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations.

The trial 2015_Irr_A presented the highest Zn accumulation for the three populations (>33 mg/kg for each population). In all trials, the SMC parents (SMC40, SMC44, and SMC33) presented a better performance in SdFe and SdZn than SCR parents (SCR9 and SCR16) (Figure 1). Conversely, SCR parents presented higher PHI, 100SdW, and Sd/Pd than the SMC parents (Figure 2). The broad-sense heritabilities for PHI, 100SdW, and Sd/Pd were high for the three populations: 0.70, 0.65, and 0.73 for SCR16 × SMC40; 0.79, 0.72, and 0.67 for SMC33 × SCR16; and 0.84, 0.79, and 0.72 for SMC44 × SCR9. In contrast, the heritabilities for the SdN and PdN were lower and not stable across trials (<0.53 for PdN and <0.47 for SdN). Heritabilities for micronutrient concentrations were intermediate ranging between 0.35 and 0.65 for SdFe and 0.58 and 0.70 for SdZn. SMC44 × SCR9 presented higher heritabilities compared to the other populations in most traits evaluated except for SdZn (Supplementary Table 1). Overall, the phenotypic datasets from field trials present a good quality for subsequent analyses.

FIGURE 2

Figure 2. Phenotypic distribution of agronomics traits. Pod harvest index (PHI) (A), 100 seed weight (100SdW) (B), pod number (PdN) (C), seed number (SdN) (D), seed per pod (Sd/Pd), and (E) evaluated of 2012, 2013, 2014, and 2015 under drought and irrigated conditions for the SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations.

Pearson correlations between traits suggested two trait clusters according to the positive and generally significant correlations within clusters and negative correlations between them. On the one hand, the micronutrient concentrations of SdFe and SdZn form a group that is positively correlated among all populations (0.71–0.76). On the other hand, a second group is formed between PdN and SdN showing a strong positive correlation between them (0.75–0.82). However, PHI only shows a stable correlation with SdN and Sd/Pd (0.49–0.56 and 0.64–0.76, respectively). These agronomic traits (PHI, PdN, SdN, and Sd/Pd) show strong negative correlations with micronutrient concentrations, especially between PHI and SdFe (−0.63 to −0.69) and SdZn (−0.72 to −0.78). 100SdW only presented weak negative correlations with SdN and PdN (−0.25 and −0.27) (Figure 3 and Supplementary Table 3).

FIGURE 3

Figure 3. Pearson’s correlation coefficients between yield components and mineral concentration traits on three populations: (A) SCR16 × SMC40, (B) SMC33 × SCR16, and (C) SMC44 × SCR9. Significance of correlations indicated as ***p < 0.0001; **p < 0.001; *p < 0.01; ns, not significant. SdFe, seed iron; SdZn, seed zinc; PHI, pod harvest index; PdN, pod number; SdN, seed number; Sd/Pd, number of seeds per pod; 100SdW, 100 seed weight. Blue color indicates positive correlation and red color indicates negative correlation.

Quantitative Trait Loci Analysis

Marker-trait associations were evaluated by QTL mapping using ICIM-ADD. In total, 79 QTLs for seven traits were identified for the three populations for eight individual trials: 14 for SdFe, 12 for SdZn, 13 for PHI, 11 for SdN, 14 for PdN, 6 for 100SdW, and 9 for Sd/Pd. Sixty-four QTLs of the 77 QTLs identified in this study are classified as “major QTL” (>10% PVE). Twenty-three hotspot regions were identified in which several QTLs were co-located (Table 2 and Supplementary Table 4). In GWAS analysis, results show low levels of significance (Supplementary Figures 3, 4).

TABLE 2

Table 2. Five most important “major QTLs” identified by inclusive composite interval mapping for additive effects (ICIM-ADD), genetic and physical position, LOD, phenotypic variation explained and additivity, mapped for SdFe, seed iron; SdZn, seed zinc; PHI, pod harvest index; PdN, pod number; SdN, seed number; Sd/Pd, number of seeds per pod in SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations evaluated in 2011, 2012, 2013, 2014, and 2015.

Quantitative Trait Loci for Fe accumulation were identified on seven chromosomes (Pv03, Pv05, Pv06, Pv07, Pv08, Pv09, and Pv011). Ten major QTLs were identified; SdFe6.3 in the 2013_Dro_B trial was the QTL with the highest PVE (19.47%) with an additive positive effect of 5.04 mg/kg, while SdFe6.2 had an additive effect of 6.28 mg/kg. In all QTLs found for SdFe, the additive positive effects came from SMC parents (SMC44 and SMC33) except for SdFe7.1. SdFe3.1 was the only QTL expressed in more than one trial. QTLs for SdZn were located on Pv02, Pv03, Pv04, Pv05, Pv06, Pv07, Pv08, and Pv10, of which eight were major QTLs. SdZn5.1 was the QTL that explained the greatest phenotypic variance on the trait (18.73%), with an additive effect of 2.74 mg/kg. SdZn2.3 was the only QTL with an additive positive effect from an SCR parent (1.86 mg/kg), while the other QTL for SdZn with additive positive effect came from SMC44 or SMC33. SdZn8.2 presented a higher additive effect compared with other QTL (3.57 mg/kg) (Table 2 and Supplementary Table 4). SdZn10.1 was identified in two trials.

Twelve QTLs for PHI were found on Pv01, Pv02, Pv03, Pv05, Pv07, Pv08, Pv09, and Pv11. Seven of these QTLs were classified as major QTLs. PHI7.2 presented the highest PVE and highest positive additive effect (20.61% and 13.54%, respectively) in 2014_Dro_B under drought stress. All QTLs with positive additive effect for this trait originates in an SCR parent, except PHI7.3. On Pv07, three different loci controlling PHI were identified, and two of them were present in two trials. Eight major QTLs for SdN were identified with SdN8.1 as the QTL with the highest PVE (18.59%) and SdN7.1 with the highest additive effect (186.93 seeds). Eleven major QTLs for PdN were found. PdN7.1 was the QTL with the highest phenotypic variance of 17.55% and PdN4.1 had the highest additive effect of 65.33 pods. Most QTLs found for PdN with positive effect came from the SMC parent except for PdN6.1, Pdn7.2, and PdN8.1. Six QTLs for seed weight were identified on Pv06, Pv07, Pv08, and Pv11, 100SdW7.1 was the QTL with the highest PVE (16.62%) and 100SdW8.2 from SCR9 presented a high additive effect (2.44 g). The sources for all QTL with positive additive allelic effects for 100SdW were SCR parents except for 100SdW6.1. Seven major QTLs for seed per pod were identified. Sd/Pd2.1 presented the highest PVE (24.50%) and Sd/Pd7.2 showed the highest positive allelic effect (0.67 seeds per pod) (Table 2 and Supplementary Table 4).

Nine QTL hotspots for yield components (PHI, SdN, PdN, SdW, and Sd/Pd) were identified. Similarly, 14 QTL hotspots for yield components and micronutrient concentration were recognized, and only one QTL hotspot was observed controlling the concentration of both Fe and Zn independently from yield components (Supplementary Table 4). Fourteen QTLs related with SdFe were identified, of which seven QTLs are co-localized with QTL related with Sd/Pd, SdN, PdN, 100SdW, and PHI. Similarly, of 12 QTLs identified for SdZn, eight were co-localized with QTL related to the yield components (Table 2 and Supplementary Table 4). In most of these “hotpots” (Hotpots 4, 5, 6, 11, 12, 16, 18, 20, 21, 22, and 23), it was observed that the positive additive allelic effects for higher micronutrient concentrations were associated with the negative allelic effect for yield component traits. The high yield component QTL was usually donated by the SCR parent, while the high mineral QTL came from the SMC parent. These strong negative associations are especially visible between SdFe and SdZn with PHI and Sd/Pd, with two strong associations with 100SdW.

Discussion

The populations analyzed showed ample phenotypic variability in terms of average accumulation of micronutrients: SCR16 × SMC40 (53.7–108.3 mg/kg for Fe and 17.5–36.8 mg/kg for Zn), SMC33 × SCR16 (50.0–101.6 mg/kg for Fe and 16.9–36.8 mg/kg for Zn) and SMC44 × SCR9 (49.2–107.5 mg/kg for Fe and 18.7–42.4 mg/kg for Zn) (Supplementary Table 1). These ranges of micronutrient concentrations are wider than those reported for other populations: Andean biparental populations AND696 × G19833 (39–79 mg/kg for Fe and 16–29 mg/kg for Zn) (Cichy et al., 2009); G21242 × G21078 (28–95 mg/kg for Fe and 17–49 mg/kg for Zn) (Blair et al., 2011); Mesoamerican biparental population G14519 × G4825 (35–97 mg/kg for Fe and 17–49 mg/kg for Zn) (Blair et al., 2010); intergene pool biparental populations DOR364 × G19833 (40–84 mg/kg for Fe and 17–42 mg/kg for Zn) (Blair et al., 2009); Cerinza × G10022 (54–100 mg/kg for Fe and 23–38 mg/kg for Zn) (Blair and Izquierdo, 2012); an African collection (23.6–78.3 mg/kg for Fe and 19–56.1 mg/kg for Zn) (Tryphone and Nchimbi-Msolla, 2010); an Andean collection (54.5–99.3 mg/kg for Fe and 21.6–39.7 mg/kg for Zn) (Katuuramu et al., 2018); a European collection (38.4–93.7 mg/kg for Fe and 18.9–43.6 mg/kg for Zn) (Caproni et al., 2020); and a multiparent advanced generation intercross (MAGIC) population (37.9–87.6 mg/kg for Fe and 18.5–39 mg/kg for Zn) (Diaz et al., 2020). Common bean has significant biofortification potential due to the wide phenotypic variability in terms of Fe and Zn concentration, producing improved lines that can accumulate >100 mg/kg of Fe and >45 mg/kg of Zn. This high accumulation capacity of micronutrients exceeds that of other legume species, such as chickpea (36–91 mg/kg of Fe), pea (26–91 mg/kg of Fe), and mung bean (35–87 mg/kg of Fe), and only surpassed by lentil (37–176 mg/kg of Fe) (Shiva et al., 2016; Upadhyaya et al., 2016; Ma et al., 2017; Jha and Warkentin, 2020; Kumar and Pandey, 2020; Sab et al., 2020).

Breeding programs need fast, accurate, and cheap methods for screening large numbers of breeding lines and germplasm for genotypes with high micronutrients contents. EDXRF is an alternative to inductively coupled plasma-optical emission spectrometry (ICP-OES), the most common tool used for measuring microelements (Guild et al., 2017b). The benefits of this method compared to ICP lie in minimal sample preparation before analysis, less labor time, and cost required [∼$0.15 per sample to X-ray diffraction (XRD) method]. In addition, it is a method with a high level of prediction. Guild et al. (2017a) reported values of prediction models of 96% for Fe and 93% for Zn compared with ICP reference concentration, with an SE of prediction of 2.5 mg/kg of Fe and 1.5 mg/kg of Zn with a maximum coefficient of variation of 5% for Fe and 3.5% for Zn.

The genetic basis underlying the micronutrient concentration have been widely studied in common bean using QTL analysis, and GWAS (Gonzalez et al., 2017; Izquierdo et al., 2018; Gunjača et al., 2021). QTL mapping is an excellent strategy for gene mapping in this case because only a bi-parental population and a small set of markers are necessary a bi-parental population and a small set of markers, which are utilized to identify the genomic regions that segregate with a trait. However, this strategy is usually low in resolution since only two alleles per locus are analyzed, and genetic recombination is limited (Islam et al., 2016). Traditional QTL mapping is highly dependent on the genetic diversity of the two parents, and the effects of the detected QTLs can vary between populations. For this reason, the validation is necessary in other populations. GWAS analysis did not show any significant association. This lacks correlation to the type of population used for the biofortification analysis, as biparental populations are more suited for QTL mapping (Diaz et al., 2021). The limitations of QTL analysis can be overcome using GWAS, which can narrow down the candidate regions using natural populations. However, its use in biparental populations presents problems by low events recombination typical of this type of populations and low allelic richness required for GWAS. On the other hand, GWAS is susceptible to report false-positive associations, and its results require validation. The number of markers used in the GWAS highly affects its results (Han et al., 2018). In this study, 14 QTLs for SdFe were identified, of which 9 QTLs have been previously reported. Astudillo et al. (2013) reported two QTLs co-located respectively to SdFe3.2 and SdFe11.2 in an inter-gene pool population; the positive allele for high iron in that study came from the Andean parental (G19833). In this study, SdFe6.3 showed the highest percentage of variance explained and was the most significant QTL identified for Fe. This QTL was identified previously in an Andean biparental population (Cichy et al., 2009), in an intergene pool population (Astudillo et al., 2013), in an Andean collection (Katuuramu et al., 2018), and a Meta-QTL analysis using different populations (Izquierdo et al., 2018). Other QTLs, such as SdFe6.2, SdFe6.4, SdFe7.1, SdFe8.1, and SdFe11.1, have been identified in different populations belonging to the Andean and Mesoamerican gene-pool (Blair et al., 2010, 2011; Blair and Izquierdo, 2012; Cichy et al., 2014; Izquierdo et al., 2018; Katuuramu et al., 2018; Diaz et al., 2020). SdFe6.1 was only reported previously by Diaz et al. (2020) in a Mesoamerican Multi-parent Population and the positive allele came from a biofortified genotype MIB778. It was suggested that this allele in MIB778 may have been derived from a Phaseolus dumosus parent. There are no reports about SdFe3.1 but in this study, the QTL proved to be highly important for the analyzed biofortified populations.

For SdZn, 12 QTLs were identified, of which 6 have been reported previously. SdZn6.1 has been reported in several populations of both gene pools (Cichy et al., 2009; Blair et al., 2010; Blair and Izquierdo, 2012; Izquierdo et al., 2018; Katuuramu et al., 2018). SdZn3.1 and SdZn10.1 were only identified in an inter-gene pool bi-parental population (Astudillo et al., 2013) and an Andean Diversity Panel (Katuuramu et al., 2018). Conversely, SdZn2.3 and SdZn8.2 only were reported in a population belonging to the Mesoamerican gene pool (Freyre et al., 1998; Blair et al., 2009; Izquierdo et al., 2018; Diaz et al., 2020). Mapping certain QTL in different gene pools demonstrates the complex control over micronutrient accumulation in the common bean that can vary between the Andean and Mesoamerican genotypes. It has been reported that Andean genotypes showed higher Fe concentrations than Mesoamerican genotypes and that Mesoamerican lines presented higher Zn contents than Andean lines (Islam et al., 2002). The high degree of correspondence with other reports for both Fe and Zn lends credibility to the current results.

One of the main challenges that common bean breeders working on biofortification face is the negative correlation between iron accumulation and yield (Raatz, 2018). Genotypes with high iron accumulation have been associated with poor yield and lower yield components (fewer pods, fewer seeds, fewer seeds per pod, or poorly filled seed) (Beebe, 2020). This same tendency is reflected in the correlations presented in Figure 3 and Supplementary Table 3 and has been reported previously (Diaz et al., 2020). Out of 14 loci for high iron, 6 were associated with low values of yield components (SdFe6.2, SdFe6.4, SdFe7.1, SdFe9.1, SdFe11.1, and SdFe11.2). These loci could contribute as much as 39% of the total additive value for high iron. Similarly, out of 12 loci for high Zn, 6 were associated with negative effects on yield components (SdZn2.1, SdZn2.2, SdZn2.3, SdZn6.1, SdZn7.1, and SdZn8.2), summing up to 49% total additive effect for Zn concentration. Several of these “false” loci were reported by other authors, so this has been a problem in the genetic analysis of mineral concentration for several years. On the other hand, the QTLs SdFe3.2, SdFe5.1, SdFe6.1, SdFe6.3, SdFe7.2, SdFe8.1, and SdFe11.3 were not associated with loci related with yield components, making them attractive for breeding through marker-assisted selection (MAS) because their selection would not adversely affect the yield. Another advantage of these identified QTLs is that some of these loci have been found in populations belonging to both gene pools, making them useful for improving Andean or Mesoamerican genotypes. The sum of additive effects not associated with poor yield could contribute as much as 40 mg/kg Fe. Even though it is likely not possible to maintain the full additive effects of these loci through genetic recombination, ample variability exists without affecting yield. Similarly, there are QTLs for Zn that were not co-located with QTLs related to yield components (SdZn3.1, SdZn4.1, SdZn5.1, SdZn5.2, SdZn8.1, and SdZn10.1). These QTLs can be implemented in the breeding program, which theoretically could increase Zn seed concentration by 10 mg/kg. Nonetheless, if limits exist to raising Fe and/or Zn concentration with the QTL studied here, additional sources of genetic variability for high minerals are being explored. Beebe (2020) reported that interspecific progeny derived from Phaseolus species that evolved on alkaline soils might be more receptive to raising iron levels. Initial phenotypic results with these interspecific crosses are promising but have not been studied yet to determine their genetic makeup.

The observed negative associations of seed mineral concentration with PHI and Sd/Pd, and in two cases with 100SdW, probably reflect simple factors of concentration of minerals. This may be due to a dilution effect of Fe and Zn caused by high carbohydrate translocation to the grain, a typical characteristic of genotypes with high yield and PHI. This has also been reported in other crops, such as chickpea (Diapari et al., 2014). Conversely, if seed mass is reduced, then as the iron in the pods is translocated to the seed, the concentration in the seed will be increased arithmetically. Furthermore, where selection for high iron has received priority in breeding, this will have unconsciously limited yield potential. These results reveal many of the causes behind the persistent negative correlation of yield and high minerals that have been the greatest challenge in breeding biofortified beans.

A strong phenotypic correlation between iron and zinc concentrations was observed here and has been reported previously (Blair et al., 2009; Cichy et al., 2009). Izquierdo et al. (2018) reported loci that contributed to high concentrations of both elements. They identified a MetaQTL (MQTL) denominated “MQTL_Fe&Zn_6.2,” containing eight QTLs identified in previous studies, such as sources of favorable alleles for Fe and Zn in both gene pools. MQTL_Fe&Zn_6.2 displayed higher PVE (27%) compared with other MQTLs identified in the same study. However, in this work, only one locus common to SdFe and SdZn was found, and this was associated with low values of PdN. This suggests that the correlation of Fe and Zn concentration that has been observed here and elsewhere may be driven in part by poor pod formation and/or poor pod and grain filling. The hotspot 12 is the only locus in which QTLs controlling the accumulation of Fe and Zn were co-located on Pv06 (∼28–29 Mpb).

The processes of iron and zinc uptake, transport, metabolism, and storage in sinks of interest are regulated by complex genetic mechanisms, of which several genes have been identified as possible candidates regulating these processes (Connorton et al., 2017). MQTL_Fe&Zn_6.2 mentioned above is located with Hotspot 12 (SdFe6.4/SdZn6.1), which is near the gene model Phvul.006G196600 that codifies a basic helix-loop-helix DNA-binding protein (bHLH). This gene is homologous to bHLH105/ILR3, an Arabidopsis gene codifying a transcriptional activator of responses to Fe deficiency, that also downregulates the expression of ferritin genes, implicated in the control of Fe homeostasis (Samira et al., 2018; Tissot et al., 2019; Gao et al., 2020). FER-LIKE REGULATOR OF IRON UPTAKE (FER) is required for induction of iron mobilization genes in Arabidopsis thaliana (Bauer et al., 2007). Phvul.005G430500 is homologous to FER and is located near SdFe5.1 (∼36.8 Mpb). Phvul.011G211900 is located near SdFe11.3 and is homologous to FRD3, which has a function of root xylem loading of iron chelator and efficient uptake out of the xylem into sink cells in Arabidopsis (Green and Rogers, 2004). AtNRAMP6 is a gene codifying a multispecific vacuolar metal transporter and plays an important role in the growth of Arabidopsis in Fe-deficient conditions (Li et al., 2019). Phvul.005G182000 (corresponding to SdFe5.1) and Phvul.009G127900 (corresponding to SdFe9.1) are homologous to AtNRAMP6. Phvul.003G143800 is homologous of AHA2, a proton pump H+-ATPASE2 controlling acidification-reduction-transport of iron uptake in Arabidopsis (Campos, 2020). The Zrt and Irt-like protein (ZIP) is well characterized for its role in Zn and Fe transporter (Astudillo et al., 2013); Phvul.002G184200 (SdZn2.2), Phvul.008G259200 (SdZn8.2), Phvul.010G059200 (SdZn10.1), and Phvul.011G058500 (SdFe11.1) have been identified as possible candidate genes belonging to ZIP metal ion transporter family.

In terms of conventional breeding, successful cases have already been reported using these populations evaluated for obtaining varieties with high accumulation of Fe and acceptable yield in different conditions: BIO-101 in Colombia (83 mg/kg of Fe and 44 mg/kg of Zn),³ SMR100 in Nicaragua (89 mg/kg of Fe and 35 mg/kg of Zn),⁴ and INTA Bio-Apante in Nicaragua⁵ were biofortified lines that used RILs of SMC40 × SCR16 population as parental materials.

An effect of drought was observed only on yield components, while Fe and Zn accumulation was not affected by the drought stress. The positive correlation between yield components has previously been reported under drought conditions (Polania et al., 2016a,b). Regarding the invariable Fe and Zn contents under stress conditions, Smith et al. (2019) explained that seeds have an evolutionary tendency to ensure adequate resources for germination and present buffering effects to maintain similar nutrient concentrations in drought and irrigated conditions, but the explicit mechanisms are not clear.

Pod harvest index expresses the ratio between the biomass in seed and pod and is a useful indicator of the remobilization of photosynthates to seed (Assefa et al., 2013). However, until now, it is not clear yet which of these two components have the key effect on PHI, and if there is a possible genetic relationship that explains the different behavior of PHI in each genotype. In this study, a stronger phenotypic correlation was observed between PHI with Sd/Pd and SdN than with PdN in all populations. However, the phenotypic correlation of PHI with 100SdW was low or not significant in this study (Figure 3), similar to results reported by Diaz et al. (2020) and Berny Mier Y Teran et al. (2019). Berny Mier Y Teran et al. (2019) suggested that sink strength due to the number of seeds is higher than sink strength based on the seed size. This could indicate that a key factor that affects PHI is the number of seeds rather than seed weight itself. In the genetic analysis, most QTLs identified for the number of seeds per pod (Sd/Pd) were co-localized with QTLs identified for SdN, and Sd/Pd has a high positive correlation with PHI. Moreover, all QTLs of PdN having a positive additive effect were associated with the SMC allele, which was parental with high micronutrient concentration but poor pod filling.

The main QTL found in this study for PHI (PHI 7.2), with the highest PVE and additive effect for the trait (Table 2), has been reported previously in Mesoamerican populations (Diaz et al., 2018; Berny Mier Y Teran et al., 2019) and in an Andean population (Sedlar et al., 2020), and could be used for breeding in both gene-pools. PHI1.1 and PHI8.1 have been identified previously in a Mesoamerican MAGIC population (Diaz et al., 2020). PHI1.2, PHI5.1, and PHI7.3 have been described in previous studies in Mesoamerican bi-parental populations (Asfaw et al., 2012; Diaz et al., 2018; Berny Mier Y Teran et al., 2019). PHI5.2 was reported by Sedlar et al. (2020) in an Andean bi-parental population. These QTLs seem to be independent of micronutrients content loci and could be used for MAS and theoretically could enhance PHI values by >20%, which discovers these markers very promising for breeding for drought tolerance and/or high yield. Other QTLs for other yield related-traits have also been reported, Berny Mier Y Teran et al. (2019) reported a QTL for seeds per pod near to Sd/Pd1.1. Similarly, Sedlar et al. (2020) reported a QTL near to Sd/Pd2.1, the QTL with the highest PVE. Other QTLs, such as PdN2.2, PdN4.1, SdN2.1, and SdN6.1, were reported previously in other Mesoamerican bi-parental populations (Diaz et al., 2017). Mukeshimana et al. (2014) found a QTL near to PdN8.2, this QTL was present in multiple trials.

One issue of PHI that has been discussed little is its role as a domestication trait associated with greater partitioning to grain, and whether the selection of genotypes with better pod filling was linked to other traits related to the domestication of common beans, such as pod indehiscence and determinate growth habit. PHI5.2 and PHI8.1 were co-located in the same genomic regions that controlled pod fiber or dehiscence for common bean (Parker et al., 2020). In the same way, PHI1.2 and PHI1.3 are near to the fin locus, a key region that controlled determinacy and photoperiod sensitivity in common bean (Kwak et al., 2012). This could indicate that PHI reflects traits that were enhanced in domestication and that loci associated with PHI could be related to domestication.

Understanding the segregation variation between the populations. Heterozygosity information in the populations is presented (Supplementary Figure 5). The utility of the early generation’s higher variability could be implemented as sources as a pedigree follow-up. In consequence, the analysis of the use of F₃ and F₄ populations was mandatory previously for further selection cycles.

Conclusion

This study evaluated the phenotypic and genetic correlation between Fe/Zn accumulation and yield components using QTL analysis in three biofortified populations developed in the CIAT bean-breeding program. The results in this work validate the negative phenotypic correlation previously described between yield and Fe/Zn accumulation in phenotypic terms. In addition, we are identifying several loci controlling both traits. We hypothesize this strong correlation is a dilution effect of Fe and Zn caused by high carbohydrate translocation to grain. Conversely, loci associated with high minerals and low values of yield components should be avoided. On the other hand, we observed QTL for Fe/Zn segregating independently of the other QTL for yield components, possibly mediated by transporters with affinity to minerals. The selection of these regions will enable enhanced levels of Fe/Zn and should not affect yield performance. QTLs observed in this and previous studies suggest similar genetic control in both gene pools at some loci and these QTLs show a high potential to be implemented in a breeding program focused on biofortification.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/LI31YG.

Author Contributions

SB conceived the study and designed the populations. CC and MG carried out the breeding populations advancement. JP did the field trial evaluation. SD did the phenotypic data analysis and linkage analysis. BR supervised the molecular laboratory. SD and DA-S performed the genotypic data analysis. SD and SB wrote the manuscript. All authors read and approved the final version of the manuscript.

Funding

This study was funded by Tropical Legumes III – Improving Livelihoods for Smallholder Farmers: Enhanced Grain Legume Productivity and Production in Sub-Saharan Africa and South Asia (OPP1114827) projects funded by the Bill & Melinda Gates Foundation. This project was supported by the CGIAR Research Program on Grain Legumes and we would like to thank USAID for their contributions, particularly in the form of the “Global Hunger and Food Security Research Strategy: Climate Resilience, Nutrition, and Policy – Feed the Future Innovation Lab for Climate Resilience in Beans (CRIB)” Project (prime award number AID-0AA-A-13-00077, sub-award number 4943-CIDT-AID-0077).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank the contributions and technical support made by the bean breeding and physiology team at the International Center for Tropical Agriculture (CIAT) for phenotyping work and population development. We are also thankful to HarvestPlus for providing support on the micromineral analysis work and the Data management team at CIAT for making the datasets used in this work available with open access.

Supplementary Material

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

Supplementary Figure 1 | Precipitation, maximum, and minimum temperatures during trials at Palmira, Colombia.

Supplementary Figure 2 | Histograms of frequency distribution for Fe and Zn in seven trials of the SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations.

Supplementary Figure 3 | Genome-wide association analysis for Fe accumulation (SdFe) showing Manhattan and Q-Q plot for SMC × SCR populations with 24,785 markers in seven trials. The Bonferroni correction threshold (p = 0.05) is presented as a red horizontal line based on the number of markers for the population.

Supplementary Figure 4 | Genome-wide association analysis for Zn accumulation (SdZn) showing Manhattan and Q-Q plot for SMC × SCR populations with 24,785 markers in seven trials. The Bonferroni correction threshold (p = 0.05) is presented as a red horizontal line based on the number of markers for the population.

Supplementary Figure 5 | Heterozygosity of individuals and heterozygosity of markers in the SMC × SCR populations.

Supplementary Table 1 | Phenotypic variability and broad-sense heritability (H²) for evaluated traits in the trials 2011, 2012, 2013, 2014, and 2015 of the SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations.

Supplementary Table 2 | Component variance, F statistic evaluated for seven trials to SdFe, seed iron; SdZn, seed zinc; PHI, pod harvest index; PdN, pod number; SdN, seed number; Sd/Pd, number of seeds per pod; 100SdW, 100 seed weight in SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations.

Supplementary Table 3 | Pearson’s correlations coefficients between evaluated traits on three populations SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 for each trial and each trait. Significance of correlations indicated as ^***p < 0.0001; ^**p < 0.001; *p < 0.01; ns, not significant. SdFe, seed iron; SdZn, seed zinc; PHI, pod harvest index; PdN, pod number; SdN, seed number; Sd/Pd, number of seeds per pod; 100SdW, 100 seed weight.

Supplementary Table 4 | Details of all QTL identified by inclusive composite interval mapping for additive effects (ICIM-ADD), genetic and physical position, LOD, phenotypic variation explained mapped for six traits in SCR16 × SMC40, SMC33 × SCR16, and SMC44 × SCR9 populations evaluated in 2011, 2012, 2013, 2014, and 2015.

Footnotes

1. ^ http://www.cropontologycurationtool.org/
2. ^ https://hudsonalpha.org/
3. ^ https://lac.harvestplus.org/en/new-biofortified-bean-varieties-bio-101-bio-107/
4. ^ https://inta.gob.ni/wp-content/uploads/2019/11/Brochure-Frijol-Biofortificado-HarvesPlus-2019.pdf
5. ^ https://lac.harvestplus.org/ometepe-biofortificados-nicaragua/

References

Amanullah, S., Gao, P., Agyei, B., Saroj, A., Yang, T., Liu, S., et al. (2021). Genetic linkage mapping and QTLs identification for morphology and fruit quality related traits of melon by SNP based CAPS markers. Sci. Hortic. 278:109849. doi: 10.1016/j.scienta.2020.109849

CrossRef Full Text | Google Scholar

Asfaw, A., Blair, M., and Struik, P. (2012). Multienvironment quantitative trait loci analysis for photosynthate acquisition, accumulation, and remobilization traits in common bean under drought stress. G3 Genes Genomes Genetics 2, 579–595. doi: 10.1534/g3.112.002303

PubMed Abstract | CrossRef Full Text | Google Scholar

Assefa, T., Beebe, S., Rao, I., Cuasquer, J., Duque, M., Rivera, M., et al. (2013). Pod harvest index as a selection criterion to improve drought resistance in white pea bean. Field Crops Res. 148, 24–33. doi: 10.1016/j.fcr.2013.04.008

CrossRef Full Text | Google Scholar

Astudillo, C., Fernandez, A., Blair, M., and Cichy, K. (2013). The phaseolus vulgaris ZIP gene family: identification, characterization, mapping, and gene expression. Front. Plant Sci. 4:286. doi: 10.3389/fpls.2013.00286

PubMed Abstract | CrossRef Full Text | Google Scholar

Bates, D., Mächler, M., Bolker, B., and Walker, S. (2015). Fitting linear mixed-effects models using Lme4. J. Stat. Softw. 67, 1–48. doi: 10.18637/jss.v067.i01

CrossRef Full Text | Google Scholar

Bauer, P., Ling, H., and Guerinot, M. (2007). FIT, the fer-like iron deficiency induced transcription factor in Arabidopsis. Plant Physiol. Biochem. 45, 260–261. doi: 10.1016/j.plaphy.2007.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaver, J., and Osorno, J. (2009). Achievements and limitations of contemporary common bean breeding using conventional and molecular approaches. Euphytica 168, 145–175. doi: 10.1007/s10681-009-9911-x

CrossRef Full Text | Google Scholar

Beebe, S. (2012). Common bean breeding in the tropics. Plant Breed. Rev. 36, 357–426. doi: 10.1002/9781118358566

CrossRef Full Text | Google Scholar

Beebe, S. (2020). Biofortification of common bean for higher iron concentration. Front. Sustain. Food Syst. 4:573449. doi: 10.3389/fsufs.2020.573449

CrossRef Full Text | Google Scholar

Beebe, S., Gonzalez, A., and Rengifo, J. (2000). Research on trace minerals in the common bean. Food Nutr. Bull. 21, 387–391. doi: 10.1177/156482650002100408

CrossRef Full Text | Google Scholar

Beebe, S., Rao, I., Cajiao, C., and Grajales, M. (2008). Selection for drought resistance in common bean also improves yield in phosphorus limited and favorable environments. Crop Sci. 48, 582–592. doi: 10.2135/cropsci2007.07.0404

PubMed Abstract | CrossRef Full Text | Google Scholar

Berny Mier Y Teran, J., Konzen, E., Palkovic, A., Tsai, S., Rao, I., et al. (2019). Effect of drought stress on the genetic architecture of photosynthate allocation and remobilization in pods of common bean (Phaseolus vulgaris L.), a key species for food security. BMC Plant Biol. 19:171. doi: 10.1186/s12870-019-1774-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Blair, M., Astudillo, C., Grusak, M., Graham, R., and Beebe, B. (2009). Inheritance of seed iron and zinc concentrations in common bean (Phaseolus vulgaris L.). Mol. Breed. 23, 197–207. doi: 10.1007/s11032-008-9225-z

CrossRef Full Text | Google Scholar

Blair, M., Astudillo, C., Rengifo, J., Beebe, S., and Graham, R. (2011). QTL analyses for seed iron and zinc concentrations in an intra-genepool population of andean common beans (Phaseolus vulgaris L.). Theor. Appl. Genet. 122, 511–521. doi: 10.1007/s00122-010-1465-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Blair, M., and Izquierdo, P. (2012). Use of the advanced backcross-QTL method to transfer seed mineral accumulation nutrition traits from wild to andean cultivated common beans. Theor. Appl. Genet. 125, 1015–1031. doi: 10.1007/s00122-012-1891-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Blair, M., Medina, J., Astudillo, C., Rengifo, J., Beebe, S., Machado, G., et al. (2010). QTL for seed iron and zinc concentration and content in a mesoamerican common bean (Phaseolus vulgaris L.) population. Theor. Appl. Genet. 121, 1059–1070. doi: 10.1007/s00122-010-1371-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bouis, H. (2003). Micronutrient fortification of plants throught plant breeding: can it improve nutrition in man at low cost? Proc. Nutr. Soc. 62, 403–411. doi: 10.1079/pns2003262

PubMed Abstract | CrossRef Full Text | Google Scholar

Bouis, H., and Welch, R. (2010). Biofortification—a sustainable agricultural strategy for reducing micronutrient malnutrition in the Global South. Crop Sci. 50, S20–S32. doi: 10.2135/cropsci2009.09.0531

PubMed Abstract | CrossRef Full Text | Google Scholar

Broughton, W., Hernandez, G., Blair, M., Beebe, S., Gepts, P., and Vanderleyden, J. (2003). Beans (Phaseolus Spp.) – model food legumes. Plant Soil 252, 55–128. doi: 10.1023/A:1024146710611

CrossRef Full Text | Google Scholar

Campos, M. (2020). Metallic action! The dynamics of a tripartite iron uptake complex in Arabidopsis roots. Plant Physiol. 184, 1212–1213. doi: 10.1104/pp.20.01271

PubMed Abstract | CrossRef Full Text | Google Scholar

Caproni, L., Raggi, L., Talsma, E., Wenzl, P., and Negri, V. (2020). European landrace diversity for common bean biofortification: a genome-wide association Study. Sci. Rep. 10, 1–13. doi: 10.1038/s41598-020-76417-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chasapis, C., Spiliopoulou, C., Loutsidou, A., and Stefanidou, M. (2012). Zinc and human health: an update. Arch. Toxicol. 86, 521–534. doi: 10.1007/s00204-011-0775-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cichy, K., Caldas, G., Snapp, S., and Blair, M. (2009). QTL Analysis of Seed iron, zinc, and phosphorus levels in an andean bean population. Crop Sci. 49, 1742–1750. doi: 10.2135/cropsci2008.10.0605

PubMed Abstract | CrossRef Full Text | Google Scholar

Cichy, K., Fernandez, A., Kilian, A., Kelly, J., Galeano, G., Shaw, S., et al. (2014). QTL analysis of canning quality and color retention in black beans (Phaseolus vulgaris L.). Mol. Breed. 33, 139–154. doi: 10.1007/s11032-013-9940-y

CrossRef Full Text | Google Scholar

Connorton, J., 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

Cullis, B., Smith, A., and Coombes, N. (2006). On the design of the early generation variety trials with correlated data. J. Agric. Biol. Environ. Stat. 11, 381–393. doi: 10.1198/108571106X154443

PubMed Abstract | CrossRef Full Text | Google Scholar

Diapari, M., Sindhu, A., Bett, K., Deokar, A., Warkentin, T., Taran, B., et al. (2014). Genetic diversity and association mapping of iron and zinc concentrations in chickpea (Cicer arietinum L.). Genome 57, 459–468. doi: 10.1139/gen-2014-0108

PubMed Abstract | CrossRef Full Text | Google Scholar

Diaz, L., Ricaurte, J., Cajiao, C., Galeano, C., Rao, I., Beebe, S., et al. (2017). Phenotypic evaluation and QTL analysis of yield and symbiotic nitrogen fixation in a common bean population grown with two levels of phosphorus supply. Mol. Breed. 37:76. doi: 10.1007/s11032-017-0673-1

CrossRef Full Text | Google Scholar

Diaz, L., Ricaurte, J., Tovar, E., Cajiao, C., Terán, H., Grajales, M., et al. (2018). QTL analyses for tolerance to abiotic stresses in a common bean (Phaseolus vulgaris L.) population. PLoS One 13:e0202342. doi: 10.1371/journal.pone.0202342

PubMed Abstract | CrossRef Full Text | Google Scholar

Diaz, S., Ariza-Suarez, D., Izquierdo, P., Lobaton, J., De La Hoz, J., Acevedo, F., et al. (2020). Genetic mapping for agronomic traits in a MAGIC population of common bean (Phaseolus vulgaris L.) under drought conditions. BMC Genomics 21:799. doi: 10.1186/s12864-020-07213-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Diaz, S., Ariza-Suarez, D., Ramdeen, R., Aparicio, J., Arunachalam, N., Hernandez, C., et al. (2021). Genetic architecture and genomic prediction of cooking time in common bean (Phaseolus vulgaris L.). Front. Plant Sci. 11:622213. doi: 10.3389/fpls.2020.622213

PubMed Abstract | CrossRef Full Text | Google Scholar

Elshire, R., Glaubitz, J., Sun, Q., Poland, J., Kawamoto, K., Buckler, E., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One 6:e0019379. doi: 10.1371/journal.pone.0019379

PubMed Abstract | CrossRef Full Text | Google Scholar

Freyre, R., Skroch, P., Geffroy, V., Adam-Blondon, A., Shirmohamadali, A., Johnson, W., et al. (1998). Towards an integrated linkage map of common bean. 4. Development of a core linkage map and alignment of RFLP maps. Theor. Appl. Genet. 97, 847–856. doi: 10.1007/s001220050964

CrossRef Full Text | Google Scholar

Gao, F., Robe, K., Bettembourg, M., Navarro, N., Rofidal, V., Santoni, V., et al. (2020). The transcription factor BHLH121 interacts with BHLH105 (ILR3) and its closest homologs to regulate iron homeostasis in Arabidopsis. Plant Cell 32, 508–524. doi: 10.1105/tpc.19.00541

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez, A., Yuste-Lisbona, F., Fernandez-Lozano, A., Lozano, R., and Santalla, M. (2017). “Genetic mapping and QTL analysis in common bean,” in The Common Bean Genome, eds M. P. de la Vega, M. Santalla, and F. Marsolais (Berlin: Springer), 69–108.

Google Scholar

Green, L., and Rogers, 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

Guild, G., Paltridge, N., Andersson, M., and Stangoulis, J. (2017a). An energy-dispersive X-Ray fluorescence method for analysing Fe and Zn in common bean, maize and cowpea biofortification programs. Plant Soil 419, 457–466. doi: 10.1007/s11104-017-3352-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Guild, G., Parkes, E., Nutti, M., Palacio-Rojas, N., and Stangoulis, J. (2017b). High-throughput measurement methodologies for developing nutrient-dense crops. Afr. J. Food Agric. Nutr. Dev. 17, 11941–11954. doi: 10.18697/ajfand.78.HarvestPlus06

CrossRef Full Text | Google Scholar

Gunjača, J., Carović, K., Lazarević, B., Vidak, M., Petek, M., Liber, Z., et al. (2021). Genome-wide association studies of mineral content in common bean. Front. Plant Sci. 12:636484. doi: 10.3389/fpls.2021.636484

PubMed Abstract | CrossRef Full Text | Google Scholar

Hambidge, M. (2000). Zinc and health: current status and future directions: introduction. J. Nutr. 130(5 Suppl.), 1344–1349.

Google Scholar

Han, K., Lee, H., Ro, N., Hur, O., Lee, J., Kwon, J., et al. (2018). QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum. Plant Biotechnol. J. 16, 1546–1558. doi: 10.1111/pbi.12894

PubMed Abstract | CrossRef Full Text | Google Scholar

Hummel, M., Hallahan, B., Brychkova, G., Ramirez-Villegas, J., Guwela, V., Chataika, B., et al. (2018). Reduction in nutritional quality and growing area suitability of common bean under climate change induced drought stress in Africa. Sci. Rep. 8:16187. doi: 10.1038/s441598-018-33952-4

CrossRef Full Text | Google Scholar

Islam, F., Basford, K., Jara, C., Redden, R., and Beebe, S. (2002). Seed compositional and disease resistance differences among gene pools in cultivated common bean. Genet. Resour. Crop Evol. 49, 285–293. doi: 10.1023/A:1015510428026

CrossRef Full Text | Google Scholar

Islam, M., Thyssen, G., Jenkins, J., Zeng, L., Delhom, C., McCarty, J., et al. (2016). A MAGIC population-based genome-wide association study reveals functional association of GhRBB1_A07 gene with superior fiber quality in cotton. BMC Genomics 17:903. doi: 10.1186/s12864-016-3249-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Izquierdo, P., Astudillo, C., Blair, M., Iqbal, A., Raatz, B., and Cichy, K. (2018). Meta-QTL analysis of seed iron and zinc concentration and content in common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 131, 1645–1658. doi: 10.1007/s00122-018-3104-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Jha, A., and Warkentin, D. (2020). Biofortification of pulses crops: status and future perspectives. Plants 9:73. doi: 10.3390/plants9010073

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, R., Nelson, G., Troyer, J., Lautenberger, J., Kessing, B., Winkler, C., et al. (2010). Accounting for multiple comparisons in a Genome-Wide Association Study (GWAS). BMC Genomics 11:724. doi: 10.1186/1471-2164-11-724

PubMed Abstract | CrossRef Full Text | Google Scholar

Katuuramu, D., Hart, J., Porch, T., Grusak, M., Glahn, R., and Cichy, K. (2018). Genome-wide association analysis of nutritional composition-related traits and iron bioavailability in cooked dry beans (Phaseolus vulgaris L.). Mol. Breed. 38, 1–18. doi: 10.1007/s11032-018-0798-x

CrossRef Full Text | Google Scholar

Kumar, S., and Pandey, G. (2020). Biofortification of pulses and legumes to enhance nutrition. Heliyon 6:e03682. doi: 10.1016/j.heliyon.2020.e03682

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwak, M., Toro, O., Debouck, D., and Gepts, P. (2012). Multiple origins of the determinate growth habit in domesticated common bean (Phaseolus vulgaris). Ann. Bot. 110, 1573–1580. doi: 10.1093/aob/mcs207

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., and Salzberg, S. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Lauer, J., Bijl, C., Grusak, M., Baenziger, P., Boote, K., Lingle, S., et al. (2012). The scientific grand challenges of the 21st century for the crop science society of America. Crop Sci. 52, 1003–1010. doi: 10.2135/cropsci2011.12.0668

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Wang, Y., Zheng, L., Li, Y., Zhou, X., Li, J., et al. (2019). The intracellular transporter AtNRAMP6 is involved in fe homeostasis in Arabidopsis. Front. Plant Sci. 10:1124. doi: 10.3389/fpls.2019.01124

PubMed Abstract | CrossRef Full Text | Google Scholar

Lobaton, J., Miller, T., Gil, J., Ariza, D., de la Hoz, J., Soler, A., et al. (2018). Resequencing of common bean identifies regions of inter-gene pool introgression and provides comprehensive resources for molecular breeding. Plant Genome 11, 1–21. doi: 10.3835/plantgenome2017.08.0068

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, A., Cacoub, P., Macdougall, I., and Peyrin-Biroulet, L. (2016). Iron deficiency anaemia. Lancet 387, 907–916. doi: 10.1016/S0140-6736(15)60865-0

CrossRef Full Text | Google Scholar

Ma, Y., Coyne, C., Grusak, M., Mazourek, M., Cheng, P., Main, D., et al. (2017). Genome-wide SNP identification, linkage map construction and QTL mapping for seed mineral concentrations and contents in pea (Pisum sativum L.). BMC Plant Biol. 17:43. doi: 10.1186/s12870-016-0956-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, L., Li, H., Zhang, L., and Wang, J. (2015). QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 3, 269–283. doi: 10.1016/j.cj.2015.01.001

CrossRef Full Text | Google Scholar

Mukeshimana, G., Butare, L., Cregan, P., Blair, M., and Kelly, J. (2014). Quantitative trait loci associated with drought tolerance in common bean. Crop Sci. 54, 923–938. doi: 10.2135/cropsci2013.06.0427

PubMed Abstract | CrossRef Full Text | Google Scholar

Parker, T., Berny Mier y Teran, J., Palkovic, A., Jernstedt, J., and Gepts, P. (2020). Pod indehiscence is a domestication and aridity resilience trait in common bean. New Phytol. 225, 558–570. doi: 10.1111/nph.16164

PubMed Abstract | CrossRef Full Text | Google Scholar

Perea, C., De La Hoz, J., Cruz, D., Lobaton, J., Izquierdo, P., Quintero, J., et al. (2016). Bioinformatic analysis of genotype by sequencing (GBS) data with NGSEP. BMC Genomics 17(Supple. 5):539. doi: 10.1186/s12864-016-2827-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Petry, N., Boy, E., Wirth, J., and Hurrell, R. (2015). Review: the potential of the common bean (Phaseolus vulgaris) as a vehicle for iron biofortification. Nutrients 7, 1144–1173. doi: 10.3390/nu7021144

PubMed Abstract | CrossRef Full Text | Google Scholar

Polania, J., Poschenrieder, C., Beebe, S., and Rao, I. (2016a). Effective use of water and increased dry matter partitioned to grain contribute to yield of common bean improved for drought resistance. Front. Plant Sci. 7:660. doi: 10.3389/fpls.2016.00660

PubMed Abstract | CrossRef Full Text | Google Scholar

Polania, J., Rao, I., Cajiao, C., Rivera, M., Raatz, B., and Beebe, S. (2016b). Physiological traits associated with drought resistance in andean and mesoamerican genotypes of common bean (Phaseolus vulgaris L.). Euphytica 210, 17–29. doi: 10.1007/s10681-016-1691-5

CrossRef Full Text | Google Scholar

Raatz, B. (2018). Biofortification of grain legumes. Achiev. Sustain. Cultiv. Grain Legumes 1, 177–200. doi: 10.19103/as.2017.0023.06

CrossRef Full Text | Google Scholar

Sab, S., Lokesha, R., Mannur, D. M., Somasekhar, Jadhav, K., Mallikarjuna, B. P., et al. (2020). Genome-wide SNP discovery and mapping QTLs for seed iron and zinc concentrations in chickpea (Cicer arietinum L.). Front. Nutr. 7:559120. doi: 10.3389/fnut.2020.559120

PubMed Abstract | CrossRef Full Text | Google Scholar

Samira, R., Li, B., Kliebenstein, D., Li, C., Davis, E., Gillikin, J., et al. (2018). The BHLH transcription factor ILR3 modulates multiple stress responses in Arabidopsis. Plant Mol. Biol. 97, 297–309. doi: 10.1007/s11103-018-0735-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmutz, J., McClean, P., Mamidi, S., Wu, G., Cannon, S., Grimwood, J., et al. (2014). A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 46, 707–713. doi: 10.1038/ng.3008

PubMed Abstract | CrossRef Full Text | Google Scholar

Sedlar, A., Zupin, M., Maras, M., Razinger, J., Šuštar-Vozlič, J., Pipan, B., et al. (2020). QTL mapping for drought responsive agronomic traits associated with physiology, phenology, and yield in an andean intra-gene pool common bean population. Agronomy 10:225. doi: 10.3390/agronomy10020225

CrossRef Full Text | Google Scholar

Shiva, Y., Singh, U., Prasad, R., and Kaur, R. (2016). Agronomic interventions for micronutrient biofortification of pulses. Indian J. Agron. 61, S161–S172. doi: 10.1002/jsfa.11059

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, M., Veneklaass, E., Polania, J., Rao, I., Beebe, S., and Merchant, A. (2019). Field drought conditions impact yield but not nutritional quality of the seed in common bean (Phaseolus vulgaris L.). PLoS One 14:e0217099. doi: 10.1371/journal.pone.0217099

PubMed Abstract | CrossRef Full Text | Google Scholar

Stangoulis, J., and Sison, C. (2008). Crop sampling protocols for micronutriens analysis. Harvest Plus Tech. Monogr. Ser. 7, 1–20.

Google Scholar

Steinbicker, A., and Muckenthaler, M. (2013). Out of balance-systemic iron homeostasis in iron-related disorders. Nutrients 5, 3034–3061. doi: 10.3390/nu5083034

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Y., Liu, X., Wang, J., Li, M., Wang, Q., Tian, F., et al. (2016). GAPIT Version 2: an enhanced integrated tool for genomic association and prediction. Plant Genome 9, 1–9. doi: 10.3835/plantgenome2015.11.0120

PubMed Abstract | CrossRef Full Text | Google Scholar

Tello, D., Gil, J., Loaiza, C., Riascos, J., Cardozo, J., and Duitama, J. (2019). NGSEP3: accurate variant calling across species and sequencing protocols. Bioinformatics 35, 4716–4723. doi: 10.1093/bioinformatics/btz275

PubMed Abstract | CrossRef Full Text | Google Scholar

Tissot, N., Robe, K., Gao, F., Grant-Grant, S., Boucherez, J., Bellegarde, F., et al. (2019). Transcriptional integration of the responses to iron availability in Arabidopsis by the BHLH factor ILR3. New Phytol. 223, 1433–1446. doi: 10.1111/nph.15753

PubMed Abstract | CrossRef Full Text | Google Scholar

Tryphone, G., and Nchimbi-Msolla, S. (2010). Diversity of common bean (Phaseolus vulgaris l.) genotypes in iron and zinc contents under screenhouse conditions. Afr. J. Agric. Res. 5, 738–747. doi: 10.5897/AJAR10.304

CrossRef Full Text | Google Scholar

Turner, S. (2018). Qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. J. Open Sour. Softw. 3:731. doi: 10.21105/joss.00731

CrossRef Full Text | Google Scholar

Upadhyaya, H., Bajaj, D., Das, S., Kumar, V., Gowda, C. L. L., Sharma, S., et al. (2016). Genetic dissection of seed-iron and zinc concentrations in chickpea. Sci. Rep. 6:24050. doi: 10.1038/srep24050

PubMed Abstract | CrossRef Full Text | Google Scholar

VanRaden, P. (2008). Efficient methods to compute genomic predictions. J. Dairy Sci. 91, 4414–4423. doi: 10.3168/jds.2007-0980

PubMed Abstract | CrossRef Full Text | Google Scholar