Skip to main content
  • Research article
  • Open access
  • Published:

Marker-based linkage map of Andean common bean (Phaseolus vulgarisL.) and mapping of QTLs underlying popping ability traits



Nuña bean is a type of ancient common bean (Phaseolus vulgaris L.) native to the Andean region of South America, whose seeds possess the unusual property of popping. The nutritional features of popped seeds make them a healthy low fat and high protein snack. However, flowering of nuña bean only takes place under short-day photoperiod conditions, which means a difficulty to extend production to areas where such conditions do not prevail. Therefore, breeding programs of adaptation traits will facilitate the diversification of the bean crops and the development of new varieties with enhanced healthy properties. Although the popping trait has been profusely studied in maize (popcorn), little is known about the biology and genetic basis of the popping ability in common bean. To obtain insights into the genetics of popping ability related traits of nuña bean, a comprehensive quantitative trait loci (QTL) analysis was performed to detect single-locus and epistatic QTLs responsible for the phenotypic variance observed in these traits.


A mapping population of 185 recombinant inbred lines (RILs) derived from a cross between two Andean common bean genotypes was evaluated for three popping related traits, popping dimension index (PDI), expansion coefficient (EC), and percentage of unpopped seeds (PUS), in five different environmental conditions. The genetic map constructed included 193 loci across 12 linkage groups (LGs), covering a genetic distance of 822.1 cM, with an average of 4.3 cM per marker. Individual and multi-environment QTL analyses detected a total of nineteen single-locus QTLs, highlighting among them the co-localized QTLs for the three popping ability traits placed on LGs 3, 5, 6, and 7, which together explained 24.9, 14.5, and 25.3% of the phenotypic variance for PDI, EC, and PUS, respectively. Interestingly, epistatic interactions among QTLs have been detected, which could have a key role in the genetic control of popping.


The QTLs here reported constitute useful tools for marker assisted selection breeding programs aimed at improving nuña bean cultivars, as well as for extending our knowledge of the genetic determinants and genotype x environment interaction involved in the popping ability traits of this bean crop.


Popbean or nuña bean (Phaseolus vulgaris L., Fabaceae) is traditionally grown in the Andean highlands of South America at 2,000-3,000 meters above sea level, where it is commonly sold in local markets or consumed at home, and it is thought to be an ancient pre-ceramic landrace [1]. It seems probable that nuña beans originated in the Andes, where they are sympatric with wild and primitive common bean populations in certain parts of Peru and Bolivia, and they may have been present in the early stages of Andean agriculture [2, 3]. The first selection pressures leading to domestication of common bean could have resulted in the development of popping beans, and it appears that toasting grains was a well-established tradition in the Andes and possibly in Mesoamerica, where early maize races have also been used for popping. However, no evidence of nuña beans has been found in Mesoamerica, most likely due to genetic differences between the Mesoamerican and Andean gene pools [1]. This, among other factors, may explain their contrasting popping ability and photoperiod response [4].

Popbean is tropical in appearance, with aggressive indeterminate type 4 growth habit [5] and day-length sensitivity since it requires nights of at least 11–12 h for flowering induction. Growing nuña bean in temperate areas requires the development of cultivars that are insensitive to photoperiod, like the modern dry bean cultivars. Furthermore, with a view to commercial production, determinate growth habit is desirable, since it reduces growing and harvesting costs. The foremost trait that distinguishes popbean from all other types of bean is the ability to expand the cotyledons after grains explode in response to heating, which is referred to as popping expansion, similar to popcorn, although the popping mechanism is different. In popcorn, the endosperm is liquefied and explosive pressure builds up in the pericarp [6, 7], and both pericarp thickness and endosperm starch type have been attributed to the popping ability. Thus, there are two crucial factors that influence the popping of popcorn: whether a kernel can pop or not, and if so, to what extent. Popping rate and flake size seemingly correspond to these two factors, which together determine high popping volume [8]. In contrast, popping in nuña beans is the result of pressurized steam trapped within and between the mesophyll cells in the cotyledons [9, 10]. The nutritive value analysis of nuña bean revealed that they have a lower mean content of protein, phosphorus, iron, and boron than dry bean varieties, and a higher level of copper and starch, which may be related to their unique texture and taste similar to roasted peanuts. Lectins and other anti-nutritional compounds were higher in raw and boiled nuña samples than in toasted nuñas, while tannin levels did not change from raw to toasted treatments. Overall in-vitro digestibility was slightly lower for toasted nuñas than boiled dry bean [11]. Taken together, the nutritional features of toasted nuña beans make them a healthy snack, although commercial production would require the genetic improvement of other agronomic traits, particularly the day-length sensitivity that has likely restricted production and commercialization of nuña beans in temperate regions [12, 13].

In maize, experimental evidence indicate that popping characteristics are quantitatively inherited [14, 15], controlled by multiple genes [1517], and influenced by environmental effects [18] and popping methods [19]. Both additive and dominant genetic effects play very important roles in the inheritance of popping characteristics [14], and several putative quantitative trait loci (QTLs) have been identified [8, 20, 21]. However, most of these studies on popping ability have focused on the maize crop, and there is little information on the ability to contribute towards popping expansion of the nuña bean. The identification of genomic regions associated with this popping ability would enable breeders to develop improved cultivars using marker assisted selection (MAS). To identify these genomic regions, it is important not only to establish accurate phenotyping methods, but also to develop a saturated molecular marker-based genetic linkage map, and then to detect QTLs associated with these popping traits. Molecular markers have emerged as powerful tools not only for mapping genes/QTLs governing economically important traits in crops [22], but also for unlocking the useful genetic diversity from unadapted/wild/unrelated germplasm [23].

Common bean is a diploid species with a genome size estimated at 450 to 650 Mb [24] that is distributed among 22 chromosomes (n = 11). The first core genetic linkage map of common bean was based on a recombinant inbred line (RIL) population resulting from the cross between representatives of the Mesoamerican (BAT93) and the Andean (Jalo EEP558) gene pools, which included 194 restriction fragment length polymorphic (RFLP) markers [25]. Nowadays the availability of microsatellite markers and their potential for anchoring new genetic maps have allowed a new expanded version of the core linkage map to be created with several hundred of these types of markers, including markers with putative gene functions. Thus, the current core common bean linkage map covers a genetic distance of 1258.8 cM and includes a total of 413 loci placed across 11 linkage groups (LGs) with an average distance between neighboring loci of 3.0 cM [26]. In addition, about fifteen RIL mapping populations and more than twenty five linkage maps have been developed, most of them created from inter-gene pool crosses, which include divergent parents showing high genetic polymorphism [27, 28]. Nevertheless, few populations have been generated from intra-gene pool crosses [2935], likely due to their low polymorphism level [36, 37] although this kind of population does not usually show phenotypic abnormalities or undesirable segregating individuals [38, 39]. Recently, a Mesoamerican saturated intra-gene pool map has been constructed by combining the genetic information from intra- and inter-gene pool segregating populations [35].

Mapping QTLs is of great importance to understand the genetic architecture underlying complex traits. Popping expansion seems to involve several genes, environmental factors, and gene-gene and gene-environment interactions. Therefore, to better understand the genetic control of popping expansion, in addition to detecting single-locus QTL effects, there is a need to identify the interactions among the different loci (epistatic QTLs, E-QTLs), as well as their environment interaction effects (QTLs x Environment, QE; and E-QTLs x Environment, E-QE). Indeed, these interactions have been successfully analysed for other complex traits in several crops species e.g. rice [40], wheat [41], cotton [42] or maize [43].

In order to identify genetic determinants of popping ability traits in common bean, we have evaluated the popping dimension index (PDI), expansion coefficient (EC) and percentage of unpopped seeds (PUS). The first two traits reflect structural changes of seeds when heating (cotyledon expansion), while the third is of considerable importance for breeding and from a commercial point of view. Besides, we have constructed a genetic map from an Andean RIL population, which includes 193 molecular markers. This has permitted the location of single-locus and epistatic QTLs involved in PDI, EC and PUS, some of which co-localized in four LGs. Markers associated to these QTLs could be used as genetic tools for MAS programs devoted to popping improvement of nuña bean cultivars.


Phenotypic variation of popping traits and correlations in RIL population

Mean values, standard deviations and ranges of variation for popping traits are shown for each environment in Table 1. PDI, EC and PUS were clearly different between both parents and varied significantly among RILs. In addition, these popping traits showed a pattern of continuous distribution, which weakly departed from normality (data not shown). Strong transgressive segregations were observed for these traits, indicating that alleles with positive effects are distributed among the parents. The broad sense heritability value of EC was rather low (0.33 ± 0.04), whereas the heritability values of PDI (0.53 ± 0.04) and PUS (0.48 ± 0.04) were moderate, suggesting that genetic gains could be obtained when selecting for these traits. The phenotypic correlations between PUS and two popping traits showed that PUS was significant and negatively correlated with PDI and EC, while EC was significant and positively correlated with PDI in all environments assayed (Table 2). Highly significant differences were found among genotypes and genotype x environment interactions for all popping traits (Table 3).

Table 1 Mean values of the popping traits analysed in the RIL population PMB0225 x PHA1037
Table 2 Phenotypic correlations among the popping traits evaluated in the RIL population PMB0225 x PHA1037
Table 3 ANOVA for the popping traits measured in the RIL population PMB0225 x PHA1037

Genetic map construction

We screened PMB0225 and PHA1037 parents for DNA polymorphism by using amplified fragment length polymorphic (AFLP), simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers. A total of 92 AFLP primer combinations were assayed, which allowed the amplification over 3,700 AFLP fragments, 279 (7.5%) of which were polymorphic between the parents. According to their inheritance pattern and reliability features, 18 combinations, producing 94 polymorphic loci, were selected for genotyping the RIL population. Furthermore, 1035 SSR and 251 SNP markers were screened, which rendered polymorphism rates of 10.2% and 7.2%, respectively. Thus, 106 SSR and 18 SNP polymorphic loci were also analysed in the RIL population for map construction. Some of the polymorphic markers (9 AFLPs, 11 SSRs and 5 SNPs) were not linked with any other marker in the existing map and so they could not be mapped. Finally, the genetic map developed from the cross PMB0225 x PHA1037 (Figure 1) was constructed with a total of 193 loci (85 AFLP, 95 SSR, and 13 SNP markers), of which 101 were dominant and 92 codominant, resulting in the formation of 12 LGs. These LGs were designated according to Pedrosa-Harand et al. [44] on the basis of 55 previously mapped common SSR markers [26, 3133, 4550], which were used as a guide for the assignment of LG number and orientation. The map spanned a total genetic distance of 822.1 cM, with an average of 68.5 cM per LG, ranging from 16.5 cM (LG 6) to 106.4 cM (LG 3). The density of markers ranged from 1.3 cM (LG 6) to 6.75 cM (LG 7), with an average of 4.3 cM per LG. A detailed description of this map is provided in Figure 1 and Table 4.

Figure 1
figure 1

Genetic linkage map of common bean based on the RIL population PMB0225 x PHA1037. Location of single-locus QTLs and E-QTLs controlling popping traits: popping dimension index (PDI), expansion coefficient (EC), and percentage of unpopped seeds (PUS). Common SSR markers to previously published maps (see text for references) are indicated in bold. Names of QTLs are listed in Tables 5, 6, and 7. Cumulative distances among markers are indicated in cM to the left of the linkage group, names of markers are shown on the right. QTLs are depicted as vertical bars to the right of the linkage groups. QTLs detected by both MapQTL and QTLNetwork software packages are indicated in black, QTLs identified only by MapQTL are shown in red, and QTLs detected only by QTLNetwork are represented in blue. Epistatic interactions between QTLs are represented with numbered stars.

Table 4 Distribution of molecular markers on the linkage map constructed from the RIL population PMB0225 x PHA1037

Single environment QTL analysis of popping traits

QTL analysis based on MQM mapping using MapQTL was undertaken to identify single-locus QTLs in individual environments; thus, a total of sixteen QTLs for popping traits were detected (Figure 1). Five single-locus QTLs were identified for PDI, of which one (PDI3PP) was identified in four environments, two were detected in two environments (PDI5PP and PDI7PP), and the other two were detected in only one environment (PDI1PP and PDI6PP). The percentage of phenotypic variation explained by a single QTL identified for PDI ranged from 7.5 (for PDI3PP in LD09) to 18.7% (for PDI6PP in SD10). For EC, two QTLs were detected in two or more environments (EC5PP and EC7PP) and two QTLs were identified in only one environment (EC3PP and EC6PP); in SD09 no significant QTLs were detected. The percentage of phenotypic variance for EC ranged from 7.3 (for EC7PP in LD09) to 15.9% (for EC6PP in SD10). In the case of PUS, seven single-locus QTLs were identified, four of which were detected in two environment conditions (PUS3PP, PUS5PP, PUS6PP, and PUS7PP) and three in only one environment (PUS1.2PP, PUS9PP, and PUS10PP). The percentage of phenotypic variance for PUS ranged from 6.3 (for PUS1.2PP and PUS10PP, both detected in LD11) to 16.5% (for PUS6PP in SD10). For each popping trait, QTLs with positive (alleles from PHA1037) and negative (alleles from PMB0225) additive values were identified. A complete description of the MQM mapping analysis is provided in Table 5.

Table 5 Single-locus QTLs detected for popping traits using multiple-QTL model mapping for individual environment analysis

Multiple environment QTL analysis of popping traits

In addition to the single-locus QTLs identified for each environment, a single-locus QTL analysis for multi-environment was also carried out by using QTLNetwork. As the result of one-dimensional MCIM analysis, a total of fourteen single-locus QTLs were detected with significant genetic main effects and/or QTLs x Environment (QE) effects (Figure 1). Four single-locus QTLs were detected for PDI, three of which showed only individual additive effects (PDI3PP, PDI5PP, and PDI6PP); and the remaining QTL (PDI7PP) had both individual additive and QE interaction effects. The percentage of phenotypic variation explained for PDI by a single QTL ranged from 2.2 (for PDI6PP) to 7% (for PDI3PP); in the case of QTL PDI7PP, the phenotypic variation explained by QE interaction effects ranged from 0.6 (in SD09) to 1.6% (in SD10). For EC, four single locus QTLs were detected, two of which displayed only individual additive effects (EC3PP and EC5PP) while the other two showed both individual additive and QE interaction effects (EC7PP and EC9PP). The percentage of phenotypic variance for EC ranged from 2.8 (for EC7PP) to 4.6% (for EC5PP); the phenotypic variation explained by QE interaction effects ranged from 0.8 (for EC9PP in SD10) to 1.4% (for EC7PP in SD10). Six QTLs were identified for PUS, two of which had only individual additive effects (PUS3PP and PUS5PP) while the other four displayed both individual additive and QE interaction effects (PUS1.1PP, PUS1.3PP, PUS6PP, and EC7PP). The percentage of phenotypic variance for PUS ranged from 0.6 (for PUS1.3PP) to 7.4% (for PUS3PP); the phenotypic variation explained by QE interaction effects ranged from 0.8 (for PUS1.1PP in SD10) to 1.8% (for PUS1.3PP in SD10). For the three popping traits, single-locus QTLs with positive and negative additive values were detected, indicating that alleles from both parents, PHA1037 and PMB0225, have a positive agronomical effect on the traits. A complete report of the one-dimensional genome scan analysis using QTLNetwork is given in Table 6.

Table 6 Single-locus QTLs and QTLs x Environment (QE) effects detected for the popping traits using multi-environment analysis

Epistatic and environmental interactions

Given that popping expansion seem to be a complex polygenic trait, genetic interactions may have significant effects on the phenotypic values. Therefore, two-dimensional genome scan was undertaken for multi-environment analysis using QTLNetwork to identify epistatic and environment interactions among QTLs. A total of ten epistatic QTLs (E-QTLs) involved in six epistatic interactions (Table 7 and Figure 1) were detected for the three evaluated traits. The percentage of phenotypic variance explained by the interaction of these E-QTLs was low, ranging from 0.2 (for EC) to 2.7% (for PDI). Interestingly, the ten E-QTLs identified were previously detected as single-locus QTLs, which indicated that these QTLs not only participated in epistatic interactions, but they also had an individual effect. Regarding PDI, only one epistatic interaction was detected, between QTLs E-PDI3PP and E-PDI7PP, explaining 2.7% of phenotypic variation. For EC, two epistatic interactions were identified involving four E-QTLs; the percentage of phenotypic variance explained by the interactions of E-EC3PP-E-EC9PP and E-EC5PP-E-EC7PP was 0.2% and 2.6%, respectively. E-QTL x Environment (E-QE) interaction effects were not significant for either PDI or EC popping traits. In the case of PUS, three epistatic interactions and four E-QTLs were detected, and two of the interactions showed only genetic effects and explained 0.6% (E-PUS3PP-E-PUS7PP) and 1.8% (E-PUS5PP-E-PUS6PP) of phenotypic variation, respectively. The remaining epistatic interaction (E-PUS3PP-E-PUS6PP) had both genetic and E-QE interaction effects, and the percentage of phenotypic variation obtained was 1.1% and 0.9% for genetic and E-QE interaction effects, respectively. A detailed description of the digenic epistatic interaction analysis is shown in Table 7.

Table 7 Epistatic QTLs (E-QTLs) and E-QTL x Environment (E-QE) effects detected for popping traits using multi-environment analysis


The main goal of the current study was to unravel the genetic architecture of popping ability in nuña bean. Thus, popping traits related to changes in the physical structure of seeds have been analysed on the basis of their similarity to popcorn, whose cotyledons also expand when dry grains are heated. Genetic analysis performed indicated that popping ability traits show a polygenic inheritance, making this the first work to report the genetic control of these traits in common bean. Transgressive segregation was observed for popping traits, suggesting that combinations of alleles from both parents have effects in the same direction; in fact, not only PHA1037 but also PMB0225 bear alleles with a positive effect on popping ability, a finding backed up by QTL analyses. Since transgressive segregation relies on additive genetic variation, the extreme phenotypes can be maintained and fixed through artificial selection, providing the potential for improvement of popping ability. Furthermore, the analysis of variance showed that although the genotype x environment interaction affects popping ability, this effect is fairly uniform across all genotypes, and it does not seriously compromise genotypic main effects, making progress from selection feasible.

A comprehensive QTL analysis was performed to detect single-locus QTLs, epistatic QTLs and their environment interactions on a newly created genetic linkage map. This map was constructed for an Andean intra-gene pool cross involving PMB0225 (dry bean) and PHA1037 (popbean) parents. Despite the morphological diversity observed in the Andean intra-gene pool, the low genetic polymorphism existing in this common bean germplasm hinders the development of genetic linkage maps [31, 33, 51, 52]. Our results confirmed a low polymorphism between the Andean parents, thus the overall polymorphism rate detected was 8.3%. Likewise, Blair et al. [33] screened a total of 700 SSR markers on the Andean parents G21242 and G21078, but only 74 mappable markers were found in that survey resulting in a polymorphism rate of 10.6%, comparable to the level of polymorphism here reported (i.e. 10.2% for SSR markers). However, Cichy et al. [31] found a moderate SSR polymorphism of 30% between the Andean parents used to generate a G19833 x AND696 RIL mapping population.

Interestingly, the genetic map described in this work shares 55 SSR markers with previously published common bean maps [26, 3133, 4550]. In fact, linkage associations have been found in terms of SSR marker mapping in the present map and previous maps, while the collinear order of the commonly mapped SSR loci has been generally observed although some inversions affecting SSR markers located close to one another have been detected. Since these loci were distributed throughout all LGs, they would permit the alignment of homologous LGs between maps and facilitate marker transfer across populations as well as between related species. Hence, these shared markers could be used as anchor points for map merging and syntenic analysis such as Galeano et al. [35] have recently reported for the consensus Mesoamerican intra-gene pool map.

The genetic linkage map developed herein includes 193 loci (85 AFLP, 95 SSR, and 13 SNP markers) across 12 LGs that cover a genetic distance of 822.1 cM, with an average of 4.3 cM per marker. Prior to this work, two Andean maps have been described for QTL analysis; the map depicted by Cichy et al. [31] included 167 loci that spanned a total map length of 1105 cM with an average marker density of 6.6 cM per locus. On the other hand, the map constructed by Blair et al. [33] contained 118 loci with a total map length of 726.0 cM and a mean marker density of 6.2 cM per locus. Therefore, compared to previous Andean maps [31, 33], the genetic map here developed shows a suitable marker density and genome coverage, which has permitted the first identification of popping ability QTLs.

Three closely related popping traits such as PDI, EC, and PUS have been analysed and the reliability of the QTLs associated to these traits has been enhanced by using several software programs, which decreased the risk of detecting false positive and negative QTLs [5355]. Therefore, single and multi-environment QTL analyses were performed to dissect the genetic architecture of popping ability in nuña bean. In summary, a total of nineteen single-locus QTLs were identified by MapQTL and QTLNetwork. Eleven of the fourteen QTLs identified by QTLNetwork were also detected by MapQTL; thus, the two independent approaches converged on the identification of common single-locus QTLs for PDI (PDI3PP, PDI5PP, PDI6PP, and PDI7PP), EC (EC3PP, EC5PP, and EC7PP), and PUS (PUS3PP, PUS5PP, PUS6PP, and PUS7PP). The percentage of phenotypic variation explained by the single-locus QTLs identified by MapQTL for popping traits was comparatively higher than that of the QTLNetwork. The results of multi-environment analyses showed that genetic main effects were sometimes subject to environmental modification; this could explain why we obtained a lower phenotypic variance using QTLNetwork software. Even so, for multi-environment analyses, the percentage of phenotypic variance attributable to genetic effects was as expected for a complex trait, which is governed by several small effect QTLs/genes located in different genomic regions, and where the environment interactions play an important role. Therefore, the four common single-locus QTLs detected for PDI and PUS together explained 22.2 and 21.8% of the phenotypic variance, respectively. Regarding EC, the three common QTLs explained 11.7% of the phenotypic variance in the RIL population. In addition, it was interesting to find that the common QTLs (detected with both software programs) not only were consistent over environments, but they also co-localized with QTLs for the analysed popping traits.

Overall, significant positive correlations between PDI and EC and negative correlations among PUS, and PDI and EC, together with the detection of co-localized QTLs for PDI, EC, and PUS on LGs 3, 5, 6, and 7, suggested that QTLs for popping ability are not evenly dispersed throughout the genome but rather are clustered in several genomic regions. The QTLs sign values of additive effects corresponded to the significant genotypic correlations observed among the three popping traits. Thus, the co-localized QTLs located on LGs 3, 5 and 6 showed positive (alleles from PHA1037) and negative (alleles from PMB0225) values of additive effects for PDI and EC, and PUS, respectively. Meanwhile, the opposite sign values of additive effects were found for the co-localized QTLs located on LG 7, which indicated that PMB0225 also contributes positively to popping ability. To date, research into popping ability has focused on popcorn. Thus, Babu et al. [21] detected four QTLs for popping expansion volume, five for flake volume, and five for percentage of unpopped kernels, and revealed QTLs in overlapping or mostly adjoining regions in the same chromosomes affecting two or three popping traits. Likewise, Li et al. [8] evaluated three important traits for popcorn (i.e. popping volume, flake size, and popping rate), and six chromosome regions were found to control two or three popping traits simultaneously. Hence, as in nuña bean, the detection of co-localized QTLs for popping traits suggested that either pleiotropic QTLs controlled several popping traits, or tightly linked QTLs for different traits are present together in the same genomic regions. The issue of pleiotropy versus tight linkage of QTLs may be resolved in the future through fine mapping of the target genomic regions.

Epistatic effects are often involved in complex traits, but they are difficult to confirm because of their usually small effects and environmental interactions. Genetic models for QTL mapping assuming no epistasis could lead to biased estimation of QTL parameters, and subsequently result in considerable loss of response in MAS. In popcorn, Li et al. [8] carried out a preliminary epistatic analysis and detected thirteen pairs of digenic interactions for popping ability. In the present work, several epistatic interactions were found involving all of the evaluated popping traits. A total of ten E-QTLs, involved in six epistatic interactions, were detected, and only one epistatic interaction for PUS showed significant E-QE interaction effect. Although the phenotypic variation explained by each epistatic interaction was found to be small, it is interesting to note that the genomic regions located on LGs 3, 5, 6 and 7 not only harbor QTLs that have individual genetic effects, but are also involved in epistatic interactions. Therefore, QTL analysis revealed that popping ability of nuña bean is controlled by several QTLs, which have only individual additive effects, or may also be involved in epistatic or environmental interactions, indicating that popping is inherited as a polygenic trait, and that epistasis could play a key role.

Nowadays, popping of common bean is considered an interesting agronomic trait, since it allows greater diversification of this crop as well as the commercialization of nuña bean as a new snack product. In popcorn, selection for increased expansion coefficient has been successfully achieved given its high heritability value [20, 21]. Similarly, Vorwald and Nienhuis [56] estimated that the narrow sense heritability values of fully expanded seeds after popping (popping percentage) and expansion coefficient in nuña bean were relatively high, 0.87 ± 0.07 and 0.74 ± 0.09, respectively. The broad sense heritability values calculated in the present work were moderate, suggesting that genetic gain could be obtained for popping ability in this legume species. The introgression of popping and the development of new day-length insensitive popbean cultivars would require genetic tools which facilitate efficient genotyping selection. Conventional phenotype selection methods for popping traits are laborious and time-consuming; consequently, MAS provides an efficient and cost-effective alternative that accelerates the selection of interesting genotypes. However, MAS approaches have been difficult to apply in the case of complex traits such as popping ability, because individual QTLs have small genetic effects which in many cases are also environmentally modulated. Consequently, the identification of potential candidate QTLs for MAS is crucial. Based on the results obtained in our study, the co-localized QTLs located on LGs 3, 5, 6, and 7 are good candidates for MAS, since they showed stability across significantly correlated traits, while also sharing QTLs for more than one trait, and they could be manipulated simultaneously in breeding programs. Breeding of nuña cultivars would require adapting them to temperate regions, and for this purpose it is important to improve their insensitivity to photoperiod. The Ppd gene regulates photoperiod sensitivity and it is located on LG 1 [57], while popping QTLs are located on LGs 3, 5, 6, and 7. Therefore, the use of QTL marker assisted selection would facilitate the introgression of photoperiod insensitivity without loss of popping ability. QTL pyramiding approach would also permit the combination of QTL alleles with positive effects for popping ability on a day-length insensitive genotype through molecular breeding, thus overcoming the main drawback for the production and commercialization of nuña beans in temperate regions. In this research, some RILs showed popping expansion ability and flowered independently of photoperiod conditions. These lines constitute an interesting breeding goal, and they will hopefully allow researchers to isolate the genes and to understand the molecular and physiological mechanisms underlying agronomic traits which are relevant for the genetic improvement of nuña beans.


We have developed a novel Andean genetic linkage map, which has permitted the first identification of popping ability QTLs in common bean. Our results revealed that popping ability of nuña bean is controlled by several QTLs, which have only individual additive effects or may also be involved in epistatic or environmental interactions, indicating that popping is inherited as a polygenic trait, and that epistasis could play a key role in its genetic control. Individual and multi-environment QTL analyses detected a total of nineteen single-locus QTLs, most notably those co-localized for the three popping ability traits placed on LGs 3, 5, 6, and 7. These QTLs showed an individual effect and also participated in epistatic interactions. Consequently, the co-localized QTLs for popping expansion response are useful tools for MAS breeding programs intended to improve production and adaptation of nuña bean cultivars. The results here reported can contribute towards the diversification of the nuña bean crop, which is becoming increasingly relevant as a new food product in the agro-food industry due to its nutritional and healthy properties.


Population development

Parents included two genotypes from the Andean gene pool, PMB0225 and PHA1037, belonging to Nueva Granada and Peru races, respectively. PMB0225 is a Spanish improved cultivar resistant to the bean common mosaic virus, which shows indeterminate erect growth habit type II, white flowers and large seeds. PHA1037 is a photoperiod-sensitive nuña popbean germplasm accession from Bolivia that has purple flowers and large red seeds, and possesses an indeterminate climbing growth habit type IV (Figure 2). Both parents were chosen to facilitate the introgression of the popping trait in common bean cultivar adapted to temperate areas. Individual plants of a F2 segregating population generated from the cross of PMB0225 × PHA1037 were selfed to develop 185 F7 recombinant inbred lines (RILs) by single-seed descent.

Figure 2
figure 2

Unpopped and popped seeds of the parental genotypes. Both PMB0225 (A) and PHA1037 (B) belong to the Andean gene pool of common bean. Scale bar 1 cm.

Experimental design

The 185 RILs and both parents were evaluated for popping seed traits using 15 plants per genotype in five greenhouse environments over three consecutive years (2009, 2010 and 2011). Plants were grown under long-day (LD) and short-day (SD) natural photoperiod conditions with average day and night temperatures of 25°C and 20°C, respectively. Sowing dates of LD experiments were February 20, 2009 (LD09 code), March 15, 2010 (LD10 code), and April 27, 2011 (LD11 code), while sowing dates of SD experiments were August 15, 2009 (SD09 code) and September 21, 2010 (SD10 code). For all environments the experiments were conducted in a randomized complete block design with two or three replicates of single row plots (3.0 × 0.8 m). Each plot was sown with two seeds per hill and adjusted to a crop density of about 30,000 plants/ha. Pods were harvested when they were completely dried. Seeds were removed and cleaned using a mechanical thrasher followed by hand cleaning and winnowing. Then, seeds were stored at 5°C for about one month before initiation of the present research.

Trait measurements

Three popping component traits have been measured: popping dimension index (PDI), popping expansion coefficient (EC), and percentage of unpopped seeds (PUS). Samples of 50 seeds from each treatment, replicate and environment, were popped in a Palson Denver popcorn maker (1200 W, 230 V, 50 Hz) for 150 s. Seed was considered fully popped when cotyledons had expanded sufficiently to shed the seed coat, and unpopped or partially popped when the seed coat failed to crack or no expansion of the cotyledons was observed. Seed dimensions were scored for each individual seed: length (mm) was measured as the longest distance across the seed parallel to the hilum, height (mm) as the longest distance perpendicular to length, and width (mm) as the longest distance across the hilum seed. Seed length, width and thickness before and after popping were determined from a random 10-seed sample, and PDI was recorded as [(∑popped seed dimensions-∑unpopped seed dimensions)/∑unpopped seed dimensions) x 100]. Each 50-seed sample was placed in a graduated cylinder and distilled water was added to a total volume of water and seeds of 100 mL. The total volume of water added was subtracted from the total volume to give the unpopped seed volume (unPV). The seeds were drained and patted dry with paper toweling and immediately popped to minimize absorption of water by the seeds. The volume of seeds after popping (PV) was obtained using a procedure similar to that used for unpopped seeds. The EC was defined as [(PV-unPV)/unPV] x 100. The PUS was calculated as the proportion of 50 seeds that were unpopped.

Statistical data analysis

Analysis of variance with the Generalized Linear Model (GLM) procedure was applied to analyse phenotypic data using the SAS Software [58]. Single degree-of-freedom orthogonal contrasts between parents were calculated to show significant differences. For each trait and environment, mean value, standard deviation and range of variation were calculated. The phenotypic correlation coefficients between PUS and EC, and PDI, were estimated by using PROC CORR [58]. Broad sense heritability (H2) was estimated as H2 = σ2G/(σ2e/re + σ2GE/e + σ2G), where σ2G is the estimate of genotypic or RIL variance, σ2e is the estimate of error variance, σ2GE is the estimate of genotype x environment interaction variance, r is the number of replicates per environment, and e is the number of environments. The genetic components of variance were estimated with the MIXED procedure of SAS software [58].

DNA extraction and molecular marker analyses

Total genomic DNA was isolated from young leaves as described by Doyle and Doyle [59]. DNA was kept in sterile water, visualized after electrophoresis in 1% agarose gels in 1X SB buffer (10 mM sodium boric acid), and quantified by comparison with DNA standards (Lambda phage DNA digest with HindIII; Invitrogen Life Technologies). DNA was diluted in sterile water to a stock concentration of 5–10 ng/μL and stored at −20°C for use in PCR analysis.

Analysis of AFLP markers was carried out according to the procedure described by Vos et al. [60], with some modifications [61]. A total of 500 ng of genomic DNA was digested with 5 U of MseI and EcoRI enzymes for 2 h at 37°C in a final volume of 40 μL. The DNA fragments were ligated to appropriate adapters via addition of 1 U of T4-DNA Ligase (Roche) and incubated for 1 h at 37°C. The pre-amplification reactions were performed in a volume of 20 μL using A as the selective nucleotide for the Eco primer (Eco + A) and four different Mse primers (Mse + A, Mse + C, Mse + G, and Mse + T). The PCR cycling parameters were 20 cycles at 94°C for 30 s, 56°C for 60 s, and 72°C for 60 s. Subsequent PCR amplifications were performed with primers that included three selective bases in their sequences. To detect AFLP fragments, Eco primers were labelled using a fluorescent dye (FAM, NED, PET or VIC), and the following PCR cycling parameters were used for selective amplifications: an initial cycle at 94°C for 30 s, 65°C for 30 s, and 72°C for 60 s. During the next 12 cycles, the annealing temperature was lowered by 0.7°C per cycle. The temperature conditions for the next 23 cycles were 94°C for 30 s, 56°C for 30 s, and 72°C for 60 s. The PCR products of selective amplifications were separated by capillary electrophoresis using a DNA sequencer (ABI PRISM® 3130 XL Genetic Analyzer, Applied Biosystems, USA). An internal size marker, GeneScan 500 LIZ (35–500 bp; Applied Biosystems) was added, allowing the co-loading of different labelled reactions. Data regarding selectively amplified DNA fragments were analysed with GeneMapper Software 3.7 (Applied Biosystems). Each AFLP marker name included a Keygene primer code ( followed by the fragment size in base pairs.

Different sets of SSR and SNP markers previously reported (references are shown below) were tested for polymorphism in the parental genotypes, and polymorphic loci were used for the construction of the genetic linkage map. SSR markers were named according to the respective authors (IAC-, [49, 62, 63]; ATA-, [45]; BM-, GATS-, [64, 65]; BMb-, [47]; BMc-, [46, 52]; BMd-, [66]; PV-, [67]; PVBR-, [68, 69]; PVEST-, [50]; PvM-, FJ-, [26, 70]). PCR reactions were carried out following the protocols described in the publications mentioned above, although PCR conditions were changed for some SSR markers. Data analysis of the SSR markers was performed by using either gel electrophoresis or capillary electrophoresis in an ABI PRISM® 3130 XL Genetic Analyzer (Applied Biosystems, USA). SNP markers tested in our mapping population were designated as BSNP- [71], Leg- [72], and SNP- [73]. High resolution melting technology (HRM) was employed to analyse the SNP markers using a LightScanner instrument (Idaho Technology), following the protocols described by Montgomery et al. [74].

Linkage map construction and QTL analyses

JoinMap® 4.0 software [75] was used to generate the linkage maps. Marker data were assigned to LG using a minimum logarithm of odds ratio (LOD) score of 6.0, and a recombination frequency value of 0.3. The Kosambi map function [76] was used to calculate the genetic distance between markers. The LGs were designated according to Pedrosa-Harand et al. [44].

Candidate QTL regions for popping traits were identified by using two different mapping software packages, MapQTL® 5.0 [77] and QTLNetwork 2.0 [78]. Interval mapping and multiple QTL model (MQM) mapping were used to detect single-locus QTLs for each environment separately by MapQTL. Thus, once potential QTLs were detected by interval mapping analysis, markers with higher LOD scores were selected as cofactors and tested using the automatic cofactor selection procedure (default P value cut off for elimination of a cofactor set of 0.02). Using the set of selected cofactors, MQM mapping analyses were carried out. A permutation test (1000 cycles) was used to determine the LOD threshold score at which the QTL was deemed to be present in a particular genomic region with a confidence interval of 95%. In addition, QTLNetwork software was used to identify single locus QTLs, epistatic QTLs (E-QTL) and their environment interaction effects (QTLs x Environment, QE; and E-QTLs x Environment, E-QE) across environments. The mixed-model based composite interval mapping method (MCIM) was carried out for one-dimensional genome scan to detect putative QTLs and their environment interactions, and for two-dimensional genome scan to identify epistatic effects. An experimental-wise significant level of 0.05 was designated for candidate interval selection, putative QTL detection, and QTL effect. Both testing and filtration window size were set at 10 cM, with a walk speed of 1 cM. The critical F-value to declare putative QTLs was determined by 1000 permutation test. The effects of QTLs and environment interactions were estimated by Markov Chain Monte Carlo method [79]. QTLs with only genetic effects indicated that these were expressed in the same way across environments. In addition, QTLs with environment interaction effects suggested that their expressions were environmentally dependent. The detected QTLs were designated as recommended by Miklas and Porch [80]. The genetic map and the QTLs detected were drawn using the MapChart 2.2 software [81].



Amplified fragment length polymorphic


Expansion coefficient


E-QTL x Environment effect


Epistatic QTL




Linkage groups


Logarithm of odds ratio


Marker assisted selection


Mixed-model based composite interval mapping


Multiple QTL model mapping


Popping dimension index


Percentage of unpopped seeds


QTL x Environment effect


Quantitative trait loci


Restriction fragment length polymorphic


Recombinant inbred lines




Single nucleotide polymorphism


Simple sequence repeats.


  1. Tohme J, Toro O, Vargas J, Debouck DG: Variability in Andean Nuña common beans (Phaseolus vulgaris, Fabaceae). Econ Bot. 1995, 49: 78-95. 10.1007/BF02862280.

    Article  Google Scholar 

  2. Gepts P, Osborn T, Rashka K, Bliss F: Phaseolin-protein variability in wild forms and landraces of the common bean (Phaseolus vulgaris L.): evidence for multiple centers of domestication. Econ Bot. 1986, 40: 451-468. 10.1007/BF02859659.

    Article  CAS  Google Scholar 

  3. Voysest O: Variedades de fríjol en América Latina y su origen. Cali, Colombia: Spanish Centro Internacional de Agricultura Tropical: 1983,

    Google Scholar 

  4. Kornegay J, White JW, Dominguez JR, Tejada G, Cajiao C: Inheritance of photoperiod response in Andean and Mesoamerican common bean. Crop Sci. 1993, 33: 977-984. 10.2135/cropsci1993.0011183X003300050021x.

    Article  Google Scholar 

  5. CIAT: Morfología de la planta de fríjol común (Phaseolus vulgaris L.). Guía de estudio. Cali, Colombia: Spanish. Centro Internacional de Agricultura Tropical: 1984,

    Google Scholar 

  6. Hoseney RC, Zeleznak K, Abdelrahman A: Mechanism of popcorn popping. J Cereal Sci. 1983, 1: 43-52. 10.1016/S0733-5210(83)80007-1.

    Article  Google Scholar 

  7. da Silva WJ, Vidal BC, Martins MEQ, Vargas H, Pereira C, Zerbetto M, Miranda LCM: What makes popcorn pop. Nature. 1993, 362: 417-

    Article  Google Scholar 

  8. Li YL, Dong YB, Niu SZ, Cui DQ: QTL for popping characteristics in popcorn. Plant Breed. 2007, 126: 509-514. 10.1111/j.1439-0523.2007.01372.x.

    Article  Google Scholar 

  9. Spaeth SC, Debouck DG, Tohme J, van Beem J: Microstructure of nuñas: Andean popping beans (Phaseolus vulgaris L.). Food Microstruct. 1989, 8: 263-269.

    Google Scholar 

  10. Vorwald J, Nienhuis J: Effects of seed moisture content, cooking time, and chamber temperature on nuña bean (Phaseolus vulgaris L.) popping. HortScience. 2009, 44: 135-137.

    Google Scholar 

  11. van Beem J, Kornegay J, Lareo L: Nutritive value of the nuña popping bean. Econ Bot. 1992, 46: 164-170. 10.1007/BF02930631.

    Article  CAS  Google Scholar 

  12. Kmiecik K, Nienhuis J: Development of nuña beans with temperate zone adaptation. Annu Rpt Bean Improv Coop. 1997, 40: 36-37.

    Google Scholar 

  13. Ogg B, Brick MA, Pearson C: Introgression of popping ability into dry beans adapted to Colorado. Annu Rpt Bean Improv Coop. 1998, 41: 35-36.

    Google Scholar 

  14. Dofing S, D'Croz-Mason MN, Thomas-Compton MA: Inheritance of expansion volume and yield in two dent corn × popcorn crosses. Crop Sci. 1991, 31: 715-718. 10.2135/cropsci1991.0011183X003100030035x.

    Article  Google Scholar 

  15. Ashman RB: Popcorn. Cooperative extension service, Plant disease control, Bulletin BP-4. West Lafayette: Purdue University: 1984.

    Google Scholar 

  16. Clary G: A study of the inheritance of expansion in popcorn. West Lafayette, IN: PhD thesis. Purdue University :1954.

    Google Scholar 

  17. Zeigler KE: Popcorn. Specialty Corns 2nd edition. Edited by: Hallauer AR. Boca Raton, FL, USA: CRC Press: 2001,199-234.

    Google Scholar 

  18. Li YL, Wu SW, Niu SZ, Li ZQ: Analysis of combining ability of popcorn with normal corn inbreds belonging to different heterotic groups. J Henan Agric Univ. 2003, 37: 1-5. In Chinese

    Google Scholar 

  19. Dofing SM, Thomas-Compton MA, Buck JS: Genotype × popping method interaction for expansion volume in popcorn. Crop Sci. 1990, 30 (1): 62-65. 10.2135/cropsci1990.0011183X003000010014x.

    Article  Google Scholar 

  20. Lu HJ, Bernardo R, Ohm HW: Mapping QTL for popping expansion volume in popcorn with simple sequence repeat markers. Theor Appl Genet. 2003, 106: 423-427.

    PubMed  CAS  Google Scholar 

  21. Babu R, Nair SK, Kumar A, Rao HS, Verma P, Gahalain A, Singh IS, Gupta HS: Mapping QTLs for popping ability in a popcorn x flint corn cross. Theor Appl Genet. 2006, 112: 1392-1399. 10.1007/s00122-006-0242-1.

    Article  PubMed  CAS  Google Scholar 

  22. Dudley JW: Molecular markers in plant improvement: manipulation of genes affecting quantitative traits. Crop Sci. 1993, 33: 660-668. 10.2135/cropsci1993.0011183X003300040003x.

    Article  CAS  Google Scholar 

  23. Tanksley SD, Nelson JC: Advanced backcross QTL analysis: a method for the simultaneous discovery and transfer of valuable QTLs from unadapted germplasm into elite breeding lines. Theor Appl Genet. 1996, 92: 191-203. 10.1007/BF00223376.

    Article  PubMed  CAS  Google Scholar 

  24. Bennett MD, Leitch IJ: Nuclear DNA amounts in angiosperms. Ann Bot. 1995, 76: 113-176. 10.1006/anbo.1995.1085.

    Article  CAS  Google Scholar 

  25. Freyre R, Skroch PW, Geffroy V, Adam-Blondon AF, Shirmohamadali A, Johnson W, Llaca V, Nodari R, Pereira P, Tsai SM, Tohme J, Dron M, Nienhuis J, Vallejos C, Gepts P: Towards an integrated linkage map of common bean. 4. Development of a core linkage map and alignment of RFLP. Theor Appl Genet. 1998, 97: 847-856. 10.1007/s001220050964.

    Article  CAS  Google Scholar 

  26. Hanai LR, Santini L, Aranha LEC, Pelegrinelli MHF, Gepts P, Tsai SM, Carneiro ML: Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers. Mol Breeding. 2010, 25: 25-45. 10.1007/s11032-009-9306-7.

    Article  CAS  Google Scholar 

  27. Broughton WJ, Hernandez G, Blair M, Beebe S, Gepts P, Vanderleyden J: Beans (Phaseolus spp.) - model food legumes. Plant Soil. 2003, 252: 55-128.

    Article  CAS  Google Scholar 

  28. Kelly JD, Gepts P, Miklas PN, Coyne DP: Tagging and mapping of genes and QTL and molecular marker-assisted selection for traits of economic importance in bean and cowpea. Field Crops Res. 2003, 82: 135-154. 10.1016/S0378-4290(03)00034-0.

    Article  Google Scholar 

  29. Kolkman JM, Kelly JD: QTL conferring resistance and avoidance to white mold in common bean. Crop Sci. 2003, 43: 539-548. 10.2135/cropsci2003.0539.

    Article  CAS  Google Scholar 

  30. Frei A, Blair MW, Cardona C, Beebe SE, Gu H, Dorn S: QTL mapping of resistance to Thrips palmi Karny in common bean. Crop Sci. 2005, 45: 379-387. 10.2135/cropsci2005.0379.

    Article  CAS  Google Scholar 

  31. Cichy KA, Blair MW, Galeno C, Snapp SS, Kelly JD: QTL analysis of root architecture traits and low phosphorus tolerance in an Andean bean population. Crop Sci. 2009, 49: 59-68. 10.2135/cropsci2008.03.0142.

    Article  CAS  Google Scholar 

  32. Blair MW, Medina JI, Astudillo C, Rengifo J, Beebe SE, Machado G, Graham R: QTL for seed iron and zinc concentration and content in a Mesoamerican common bean (Phaseolus vulgaris L.) population. Theor Appl Genet. 2010, 121: 1059-1070. 10.1007/s00122-010-1371-0.

    Article  PubMed  CAS  Google Scholar 

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

    Article  PubMed  CAS  Google Scholar 

  34. Blair MW, Galeano CH, Tovar E, Muñoz-Torres MC, Castrillón AV, Beebe SE, Rao IM: Development of a Mesoamerican intra-genepool genetic map for quantitative trait loci detection in a drought tolerant x susceptible common bean (Phaseolus vulgaris L.) cross. Mol Breeding. 2012, 29: 71-88. 10.1007/s11032-010-9527-9.

    Article  Google Scholar 

  35. Galeano CH, Fernandez AC, Franco-Herrera N, Cichy KA, McClean PE, Vanderleyden J, Blair MW: Saturation of an intra-gene pool linkage map: towards a unified consensus linkage map for fine mapping and synteny analysis in common bean. PLoS One. 2011, 6: 12-

    Article  Google Scholar 

  36. Welsh W, Bushuk W, Roca W, Singh SP: Characterization of agronomic traits and markers of recombinant inbred lines from intra- and interracial populations of Phaseolus vulgaris L. Theor Appl Genet. 1995, 91: 169-177.

    Article  PubMed  CAS  Google Scholar 

  37. González AM, Rodiño AP, Santalla M, De Ron AM: Genetics of intra-gene pool and inter-gene pool hybridization for seed traits in common bean (Phaseolus vulgaris L.) germplasm from Europe. Field Crops Res. 2009, 112: 66-76. 10.1016/j.fcr.2009.02.003.

    Article  Google Scholar 

  38. Koinange EMK, Gepts P: Hybrid weakness in wild Phaseolus vulgaris L. J Hered. 1992, 83: 135-139.

    Google Scholar 

  39. Singh SP, Gutierrez JA: Geographical distribution of the DL1 and DL2 genes causing hybrid dwarfism in Phaseolus vulgaris L., their association with seed size, and their significance to breeding. Euphytica. 1984, 33: 337-345. 10.1007/BF00021130.

    Article  Google Scholar 

  40. Xing YZ, Tan YF, Hua JP, Sun XL, Xu CG, Zhang QF: Characterization of the main effects, epistatic effects and their environmental interactions of QTLs on the genetic basis of yield traits in rice. Theor Appl Genet. 2002, 105: 248-257. 10.1007/s00122-002-0952-y.

    Article  PubMed  CAS  Google Scholar 

  41. Kulwal PL, Singh R, Balyan HS, Gupta PK: Genetic basis of pre-harvest sprouting tolerance using single-locus and two-locus QTL analyses in bread wheat. Funct Integr Genom. 2004, 4: 94-101. 10.1007/s10142-004-0105-2.

    Article  CAS  Google Scholar 

  42. Shen X, Zhang T, Guo W, Zhu X, Zhang X: Mapping fiber and yield QTLs with main, epistatic and QTL x environment interaction effects in recombinant inbred lines of upland cotton. Crop Sci. 2006, 46: 61-66. 10.2135/cropsci2005.0056.

    Article  CAS  Google Scholar 

  43. Yang X, Guo Y, Yan J, Zhang J, Song T, Rocheford T, Li JS: Major and minor QTL and epistasis contribute to fatty acid compositions and oil concentration in high-oil maize. Theor Appl Genet. 2009, 120: 665-678.

    Article  PubMed  Google Scholar 

  44. Pedrosa-Harand A, Porch T, Gepts T: Standard nomenclature for common bean chromosomes and linkage groups. Annu Rep Bean Improv Coop. 2008, 51: 106-107.

    Google Scholar 

  45. Blair MW, Buendía HF, Giraldo MC, Métais I, Peltier D: Characterization of AT-rich microsatellites in common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2008, 118: 91-103. 10.1007/s00122-008-0879-z.

    Article  PubMed  CAS  Google Scholar 

  46. Blair MW, Hurtado N, Chavarro CM, Muñoz-Torres M, Giraldo MC, Pedraza F, Tomkins J, Wing R: Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: an integration of the BMc series. BMC Plant Bio. 2011, 11: 50-10.1186/1471-2229-11-50.

    Article  CAS  Google Scholar 

  47. Córdoba JM, Chavarro C, Schlueter JA, Jackson SA, Blair MW: Integration of physical and genetic maps of common bean through BAC-derived microsatellite markers. BMC Genomics. 2010, 11: 436-10.1186/1471-2164-11-436.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Pérez-Vega E, Pañeda A, Rodríguez-Suárez C, Campa A, Giraldez R, Ferreira JJ: Mapping of QTLs for morpho-agronomic and seed quality traits in a RIL population of common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2010, 120: 1367-1380. 10.1007/s00122-010-1261-5.

    Article  PubMed  Google Scholar 

  49. Campos T, Oblessuc PR, Sforça DA, Cardoso JMK, Baroni RM, Sousa ACB, Carbonell SAM, Chioratto AF, Garcia AAF, Rubiano LB, Souza AP: Inheritance of growth habit detected by genetic linkage analysis using microsatellites in the common bean (Phaseolus vulgaris L.). Mol Breeding. 2011, 27: 549-560. 10.1007/s11032-010-9453-x.

    Article  Google Scholar 

  50. García RAV, Rangel PN, Brondani C, Martins WS, Melo LC, Carneiro MS, Borba TCO, Brondani RPV: The characterization of a new set of EST-derived simple sequence repeat (SSR) markers as a resource for the genetic analysis of Phaseolus vulgaris. BMC Plant Bio. 2011, 12: 41-

    Google Scholar 

  51. Blair MW, Díaz JM, Hidalgo R, Díaz LM, Duque MC: Microsatellite characterization of Andean races of common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2007, 116: 29-43. 10.1007/s00122-007-0644-8.

    Article  PubMed  CAS  Google Scholar 

  52. Blair MW, Muñoz-Torres M, Giraldo MC, Pedraza F: Development and diversity assessment of Andean-derived, gene-based microsatellites for common bean (Phaseolus vulgaris L.). BMC Plant Bio. 2009, 9: 100-10.1186/1471-2229-9-100.

    Article  Google Scholar 

  53. McElroy JP, Zhang W, Koehler KJ, Lamont SJ, Dekkers JC: Comparison of methods for analysis of selective genotyping survival data. Genet Sel Evol. 2006, 38: 637-655.

    PubMed  CAS  PubMed Central  Google Scholar 

  54. Mackay I, Powell W: Methods for linkage disequilibrium mapping in crops. Trends Plant Sci. 2007, 12: 57-63. 10.1016/j.tplants.2006.12.001.

    Article  PubMed  CAS  Google Scholar 

  55. Ravi K, Vadez V, Isobe S, Mir RR, Guo Y, Nigam SN, Gowda MVC, Radhakrishnan T, Bertioli DJ, Knapp SJ, Varshney RK: Identification of several small main-effect QTLs and a large number of epistatic QTLs for drought tolerance related traits in groundnut (Arachis hypogaea L.). Theor Appl Genet. 2011, 122: 1119-1132. 10.1007/s00122-010-1517-0.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  56. Vorwald J, Nienhuis J: Heritability and correlations among popping characteristics and seed weight in a broad-based temperate-adapted nuña bean (Phaseolus Vulgaris L.) population. Hortscience. 2009, 44: 214-216.

    Google Scholar 

  57. Gu W, Zhu J, Wallace DH, Singh SP, Weeden NF: Analysis of genes controlling photoperiod sensitivity in common bean using DNA markers. Euphytica. 1998, 102: 125-132. 10.1023/A:1018340514388.

    Article  CAS  Google Scholar 

  58. SAS Institute Inc: SAS version 9.02. Cary, NC, USA: SAS Institute Inc: 2010.

    Google Scholar 

  59. Doyle JJ, Doyle JL: Isolation of plant DNA from fresh tissue. Focus. 1990, 12: 13-15.

    Google Scholar 

  60. Vos P, Hogers R, Bleeker M, Reijans M, Van de Lee T, Hormes M, Frijters A, Pot J, Peleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA-fingerprinting. Nucleic Acids Res. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  61. Noguera FJ, Capel J, Alvarez JI, Lozano R: Development and mapping of a codominant SCAR marker linked to the andromonoecious gene of melon. Theor Appl Genet. 2005, 110: 714-720. 10.1007/s00122-004-1897-0.

    Article  PubMed  CAS  Google Scholar 

  62. Benchimol LL, Campos T, Carbonell SAM, Colombo CA, Chioratto AF, Formighieri EF, Gouvea LRL, Souza AP: Structure of genetic diversity among common bean (Phaseolus vulgaris L.) varieties of Mesoamerican and Andean origins using new developed microsatellite markers. Genet Resour Crop Evol. 2007, 54: 1747-1762. 10.1007/s10722-006-9184-3.

    Article  CAS  Google Scholar 

  63. Cardoso JMK, Oblessuc PR, Campos T, Sforc¸a DA, Carbonell SAM, Chioratto AF, Formighieri EE, Souza AP, Benchimol LL: New microsatellite markers developed from an enriched microsatellite common bean library. Pesq Agropec Bras. 2008, 43: 929-936. 10.1590/S0100-204X2008000700019.

    Article  Google Scholar 

  64. Gaitán-Solís E, Duque MC, Edwards KJ, Tohme J: Microsatellite repeats in common bean (Phaseolus vulgaris): isolation, characterization, and cross-species amplification in Phaseolus ssp. Crop Sci. 2002, 42: 2128-2136. 10.2135/cropsci2002.2128.

    Article  Google Scholar 

  65. Blair MW, Muñoz-Torres M, Pedraza F, Giraldo MC, Buendía HF, Hurtado N: Development of microsatellite markers for common bean (Phaseolus vulgaris L.) based on screening of non-enriched small insert genomic libraries. Genome. 2009, 52: 772-787. 10.1139/G09-053.

    Article  PubMed  CAS  Google Scholar 

  66. Blair MW, Pedraza F, Buendía HF, Gaitán-Solís E, Beebe SE, Gepts P, Tohme J: Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2003, 107: 1362-1374. 10.1007/s00122-003-1398-6.

    Article  PubMed  CAS  Google Scholar 

  67. Yu K, Park SJ, Poysa V, Gepts P: Integration of simple sequence repeat (SSR) markers into a molecular linkage map of common bean (Phaseolus vulgaris L.). J Hered. 2000, 91: 429-434. 10.1093/jhered/91.6.429.

    Article  PubMed  CAS  Google Scholar 

  68. Buso GSC, Amaral ZPS, Brondani RPV, Ferreira ME: Microsatellite markers for the common bean Phaseolus vulgaris. Mol Ecol Notes. 2006, 6: 252-254. 10.1111/j.1471-8286.2006.01210.x.

    Article  CAS  Google Scholar 

  69. Grisi MCM, Blair MW, Gepts P, Brondani C, Pereira PAA, Brondani RPV: Genetic mapping of a new set of microsatellite markers in a reference common bean (Phaseolus vulgaris) population BAT93 x Jalo EEP558. Genet Mol Res. 2007, 6: 691-706.

    PubMed  CAS  Google Scholar 

  70. Hanai LR, Campos T, Camargo LEA, Benchimol LL, Souza AP, Melotto M, Carbonell SAM, Chioratto AF, Consoli L, Formighieri EF, Siqueira MVBM, Tsai SM, Vieira MLC: Development, characterization, and comparative analysis of polymorphism at common bean SSR loci isolated from genic and genomic sources. Genome. 2007, 50: 266-277. 10.1139/G07-007.

    Article  PubMed  CAS  Google Scholar 

  71. Galeano CH, Fernández AC, Gómez M, Blair MW: Single strand conformation polymorphism based SNP and Indel markers for genetic mapping and synteny analysis of common bean (Phaseolus vulgaris L.). BMC Genomics. 2009, 10: 629-10.1186/1471-2164-10-629.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Hougaard BK, Madsen LH, Sandal N, Moretzsohm MC, Fredslund J, Schauser L, Nielsen AM, Rohde T, Sato S, Tabata S, Bertioli DJ, Stougaard J: Legume anchor markers link syntenic regions between Phaseolus vulgaris, Lotus japonicus, Medicago truncatula and Arachis. Genetics. 2008, 179: 2299-2312. 10.1534/genetics.108.090084.

    Article  PubMed  PubMed Central  Google Scholar 

  73. McConnell M, Mamidi S, Lee R, Chikara S, Rossi M, Papa R, McClean P: Syntenic relationships among legumes revealed using a gene-based genetic linkage map of common bean (Phaseolus vulgaris L.). Theor Appl Genet. 2010, 121: 1103-1116. 10.1007/s00122-010-1375-9.

    Article  PubMed  Google Scholar 

  74. Montgomery J, Wittwer CT, Palais R, Luming Z: Simultaneous mutation scanning and genotyping by high-resolution DNA melting analysis. Nature. 2007, 2: 59-66.

    CAS  Google Scholar 

  75. van Ooijen JH: JoinMap®4, software for the calculation of genetic linkage maps in experimental populations. Wageningen, Netherlands: Kyazma BV: 2006,

    Google Scholar 

  76. Kosambi DD: The estimation of map distances from recombination values. Ann Eugen. 1944, 12: 172-175.

    Article  Google Scholar 

  77. van Ooijen JH: MapQTL®5, software for the mapping of quantitative trait loci in experimental populations. Wageningen, Netherlands: Kyazma BV: 2004.

    Google Scholar 

  78. Yang J, Hu CC, Hu H, Yu RD, Xia Z, Ye XZ, Zhu J: QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations. Bioinformatics. 2008, 24: 721-723. 10.1093/bioinformatics/btm494.

    Article  PubMed  Google Scholar 

  79. Wang CS, Rutledge JJ, Gianola D: Bayesian analysis of mixed linear models via Gibbs sampling with an application to litter size in Iberian pigs. Genet Sel Evol. 1994, 26: 91-115. 10.1186/1297-9686-26-2-91.

    Article  PubMed Central  Google Scholar 

  80. Miklas PN, Porch T: Guidelines for common bean QTL nomenclature. Annu Rep Bean Improv Coop. 2010, 53: 202-204.

    Google Scholar 

  81. Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 77-78. 10.1093/jhered/93.1.77.

    Article  PubMed  CAS  Google Scholar 

Download references


The authors thank Quival-Frutos Secos El Nogal (Pontevedra, Spain) for technical support and Diputación de Pontevedra for farm facilities. We also thank Rosana Pereira Vianello Brondani from Embrapa Arroz e Feijão, CNPq (Brasil) for supplying some microsatellite primers. MDLF was supported by a research contract of the Xunta de Galicia. This work has been funded by grants PET2008_0167, EUI2009-04052 and AGL2011-25562 of the Ministerio de Ciencia e Innovación and PGIDI03RAG16E of the Xunta de Galicia.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Rafael Lozano.

Additional information

Authors’ contributions

FY-L carried out genetic and molecular marker analyses, developed the genetic map and drafted the manuscript. MS developed the recombinant inbred line population, performed the field trials and carried out the phenotypic evaluation of plant material, and collaborating in writing the manuscript. CC and MG-A collaborated in marker genetic analyses and contributed to data analysis. MDLF collaborated in the phenotypic analysis of segregating population and the selection of microsatellite markers. JC participated in the design of genotyping experiments, supported mapping methodologies and contributed to a critical review of the manuscript. AMDR provided significant information about the plant material, collaborated in the analysis of the phenotypic and genetic data and reviewed the manuscript. RL conceived the project, planned the research work, assisted in analysis and interpretation of the data, and edited the manuscript. All authors have read and approved the final manuscript.

Fernando J Yuste-Lisbona, Marta Santalla contributed equally to this work.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Yuste-Lisbona, F.J., Santalla, M., Capel, C. et al. Marker-based linkage map of Andean common bean (Phaseolus vulgarisL.) and mapping of QTLs underlying popping ability traits. BMC Plant Biol 12, 136 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: