Skip to main content

A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia fabaL.)



Faba bean (Vicia faba L.) is among the earliest domesticated crops from the Near East. Today this legume is a key protein feed and food worldwide and continues to serve an important role in culinary traditions throughout Middle East, Mediterranean region, China and Ethiopia. Adapted to a wide range of soil types, the main faba bean breeding objectives are to improve yield, resistance to biotic and abiotic stresses, seed quality and other agronomic traits. Genomic approaches aimed at enhancing faba bean breeding programs require high-quality genetic linkage maps to facilitate quantitative trait locus analysis and gene tagging for use in a marker-assisted selection. The objective of this study was to construct a reference consensus map in faba bean by joining the information from the most relevant maps reported so far in this crop.


A combination of two approaches, increasing the number of anchor loci in diverse mapping populations and joining the corresponding genetic maps, was used to develop a reference consensus map in faba bean. The map was constructed from three main recombinant inbreed populations derived from four parental lines, incorporates 729 markers and is based on 69 common loci. It spans 4,602 cM with a range from 323 to 1041 loci in six main linkage groups or chromosomes, and an average marker density of one locus every 6 cM. Locus order is generally well maintained between the consensus map and the individual maps.


We have constructed a reliable and fairly dense consensus genetic linkage map that will serve as a basis for genomic approaches in faba bean research and breeding. The core map contains a larger number of markers than any previous individual map, covers existing gaps and achieves a wider coverage of the large faba bean genome as a whole. This tool can be used as a reference resource for studies in different genetic backgrounds, and provides a framework for transferring genetic information when using different marker technologies. Combined with syntenic approaches, the consensus map will increase marker density in selected genomic regions and will be useful for future faba bean molecular breeding applications.


Faba bean (Vicia faba L.) is an important food and fodder crop worldwide and a staple in Middle East, Central and East Asia and North Africa. In terms of cultivation area, faba bean ranks fourth among the cool-season food legumes (close to 2.5 million hectares per year) after chickpea, pea and lentil ( Its agricultural role is currently increasing as the crop is receiving a renovated interest in European countries, North America and Australia.

Faba bean is a diploid legume crop (2n = 2× = 12) with facultative cross-pollination and has one of the largest described genomes among legumes. The genome size is ~13,000 Mb, more than 25 times larger than that of the model Medicago truncatula, and contains more than 85% of repetitive DNA [1]. The large chromosome size has allowed faba bean to become one of the cytogenetically best characterised plant species. However, it has also difficulted the development of saturated linkage maps and the identification of relevant genes/QTLs (Quantitative Trait Loci). Moreover, it precludes whole-genome shotgun assembly with next generation sequencing technologies.

Genetic linkage maps are essential tools for a wide range of genetic and breeding applications, including the study of inheritance of qualitative and quantitative traits and the identification of markers linked to relevant agronomic traits. The availability of high-density maps enhances the breeding process through the application of association analyses, map-based cloning or marker-assisted approaches. Table 1 presents a summary of the faba bean linkage maps reported so far.

Table 1 Information of relevant faba bean mapping populations highlighting the ones used to construct this composite map

Prior to 1990, only a few morphological and isozyme loci were mapped in the V. faba genome and no extended linkage groups (LGs) had been reported. Primary trisomics and translocation stocks allowed the first assignment of genes and LGs to specific chromosomes [3, 1619]. This approach was further explored to develop physically localized markers and microsatellites (or SSR-Simple Sequence Repeats) from specific chromosomic regions [20]. This led to the integration of the first genetic and physical maps and allowed the unambiguous assignation of LGs to their respective chromosomes.

Preliminary maps constructed with F2 populations were mostly based on dominant markers such as RAPDs together with morphological, isozyme, seed storage protein genes and microsatellites, which saturated different areas of the genome [58, 21, 22]. First attempts to map genes/QTLs for seed weight [6] and resistance to a parasitic plant (Orobanche crenata) and fungal diseases (Ascochyta fabae and Uromyces viciae-fabae) were reported [7, 8, 10, 23]. Using a F2 population from the cross Vf6 × Vf136, a linkage map was developed to locate QTLs controlling crenate broomrape (O. crenata) [7] and A. fabae resistance [8]. Nine of the 16 LGs reported could be assigned to specific chromosomes thanks to markers that were common with those of previous studies. Subsequently, a linkage map of an F2 population from the cross 29H × Vf136, segregating for resistance to the two pathogens, was constructed in which 6 of the 18 LGs were assigned to chromosomes [10].

These faba bean maps did not allow wider mapping comparisons, since they mostly shared dominant and anonymous markers such as RAPDs, with scarce transferrability between genotypes and legume species. Despite this limitation, marker data of 11 F2 populations (Table 1), all sharing the common female parent Vf6, were used to construct a composite linkage map [9]. After joint segregation analysis of 501 markers in 654 individuals, 192 markers were included in 14 major LGs, of which 5 were unambiguously assigned to specific chromosomes (Table 1). This composite map covered 1,559 cM and was one of the most comprehensive faba bean genetic map published to date [9].

These maps with dominant markers in F2 were followed by more precise maps constructed in the corresponding RIL populations, using co-dominant markers. In addition to microsatellites, expressed sequence tags (EST) from other legume species emerged as efficient tools in faba bean. A large number of intron-targeted primer pairs (ITAPs), developed within the Grain Legumes Integrated Project-GLIP (, was tested and mapped in two faba bean inbred populations (Vf6 × Vf136, 29H × Vf136). These were used to validate QTLs underlying broomrape and Ascochyta resistance in different environments and genetic backgrounds [15, 2426]. A third RIL population derived from cross Vf6 × Vf27 was used to construct the first exclusively gene-based genetic map in faba bean. It contained 135 ITAPs joined in 12 unassigned LGs, that spanned 1,685 cM, and allowed for the first time the study of macrosyntenic relationships between V. faba, M. truncatula, Lens culinaris and other legume species [11]. After further saturation, the map was used to identify and validate QTLs controlling flowering time and other yield-related traits [13]. Recently, a new map was reported based on the F2 population from the cross 91825 × K1563, which includes 128 SSRs markers arranged in 15 unassigned LGs [14]. Unfortunately the lack of common markers prevented comparisons with previous mapping studies.

To date 14 major genetic maps have been constructed in faba bean (Table 1). Integrating the information of multiple populations from diverse genetic backgrounds offers several advantages over individual genetic maps: (i) a larger number of loci is mapped than in single crosses, (ii) the relative position of common markers can be determined across the mapping populations, (iii) better genome coverage and opportunities to validate marker order, (iv) better assignment of LGs to chromosomes, (v) it allows comparison of genes/QTLs of interest across maps and, (vi) it provides the basis for comparing genomes between related species [2729]. Consensus genetic maps have been developed in many crops such as wheat [30], maize [31], barley [32] and rice [33], and in the legume crops soybean [34], pea [35], chickpea [36], phaseolus [37], pigeonpea [38], cowpea [39], groundnut [40] and red clover [41].

With the development of genome sequencing projects and expression studies in different model and crop legumes, the construction of a faba bean consensus genetic map has become possible. The objectives of this study were to: (1) saturate the faba bean maps developed in RIL progenies with common gene based markers to facilitate anchoring of linkage groups from different populations, (2) update the position of the most relevant faba bean QTLs controlling resistance and yield related traits using Bulked Segregant Analysis (BSA), and (3) construct a reference map integrating all the genomic information reported so far in this crop. To this aim, we fused information of 11 F2 populations and marker data of three RIL genetic maps to derive a consensus map including 729 markers and covering 4,602 cM. The six main LGs could be unambiguously assigned to their corresponding faba bean chromosomes. The map represents a significant improvement over single-population genetic maps and provides a new tool of reference for faba bean breeding and genomic approaches.


Individual maps and QTL analysis

Cross Vf6 × Vf27 (RIL1)

The first RIL1 map [13] included 258 markers joined in 16 LGs and covering 1,875 cM. The linkage groups were composed of 2–45 loci with an average marker interval of 7.3 cM. The map allowed to identify and validate QTLs controlling 5 flowering and reproductive traits [13]: days to flowering (DF), flowering length (FL), pod length (PL), number of seeds per pod (NSP) and number of ovules per pod (NOP), located mainly in chromosomes (chr.) V and VI (Additional file 1: Table S1).

For the extended RIL1 map constructed herein, 313 polymophic markers were used in the global analysis (Table 2). Of these, 273 were assembled in 19 LGs, 11 of which could be assigned to specific chromosomes. The distance covered by the map was 2,183 cM with an average marker interval of 10 cM. Sixty five of the markers are common with the other two RIL populations, 25 with RIL2 (Vf6 × Vf136) only, 15 with RIL3 (29H × Vf136) only and 25 with both (Table 3).

Table 2 Number and type of markers genotyped in each inbred population
Table 3 Common markers used as anchors for map integration

Cross VF6 × VF136 (RIL2)

The previously published RIL2 map [24, 25] was based on 277 marker loci assembled in 21 LGs (16 consisting of 3 or more markers) that span 2,857 cM with an average marker interval of 12.7 cM. In this population, 2 QTLs controlling ascochyta blight resistance (Af1 and Af2) were identified on chr. III and chr. II [24]. In parallel studies, 2 QTLs (Of1 and Of2) controlling O. foetida resistance and 4r QTLs controlling O. crenata resistance (Oc2-Oc5) were detected [25]. Oc2 and Oc3 were stable in at least two of the three environments, while Oc4 and Oc5 were only detected in one environment and thus appeared to be environment-dependent.

In an attempt to saturate the regions bearing the O. crenata and A. fabae QTLs, a BSA approach based on RAPD markers was applied. Two-hundred and eight of the 748 RAPD primers assayed in the cross revealed promising polymorphisms between at least one pair of bulks and were subsequently used to screen 14 plants individually. Only 41 of the RAPD primers maintained the expected pattern of polymorphism and were used to screen the entire RIL2 population, resolving 39 scorable polymorphic markers. Of these, 31 markers exhibited the expected segregation pattern and were thus included in the RIL data set for mapping and QTL analysis. Thirty of the 31 RAPD markers were mapped, 24 to target regions [11 to chr. VI (Oc2), 8 to chr. II (Oc3 and Af2) and 5 to chr. III (Af1)] and 6 to other LGs (Table 4; Figure 1). To increase the number of common markers between different faba bean crosses, additional ITAP markers were scored and the map was reconstructed using 404 segregating loci.

Table 4 Linkage map of Vf6 × Vf136 (RIL2)
Figure 1

Faba bean consensus map showing: (a) the six main linkage groups or chromosomes (578 loci) and (b) minor linkage groups (30) or small chromosome fragments (7). Additional markers derived from the BSA analysis are underlined. Boxes represent putative locations of QTLs. Black boxes were used for flowering and yield related traits: days to flowering (df), flowering length (fl), pod length (pl), number of ovules per pod (nop) and number of seeds per pod (nsp). Grey boxes: Ascochyta fabae (Af) QTLs. Stripped boxes: Orobanche crenata (Oc) and Orobanche foetida QTLs. Marker distance is given in cM.

The map obtained in this study consists of 364 mapped loci assembled into 21 LGs, of which 13 were assigned to specific chromosomes. The genome distance covered by the map is 3,537 cM, with an average density of one marker locus every 12.6 cM (Table 4). Fifty four markers included in RIL2 are common with the other two populations, 25 with RIL1 only, 4 with RIL3 only and 25 with both (Table 3).

As mentioned above, QTLs Oc2 and Oc3, that confer broomrape resistance, were previously validated in 2 of the 3 field trials assayed [26]. By saturating the target regions we were able to identify significant QTLs in each trial. Thus, the conservation of QTLs both across generations and environments was confirmed. In case of Oc2 (Additional file 1: Table S1B) a newly added marker, OPAG11956, was the closest to LOD peak value in both Córdoba 2003 and Córdoba 2004 datasets. Three additional QTLs were identified in chr. II (Oc14_C3 and Oc15_M4) and V (Oc15_C4). None of these QTLs was stable in the different field assays or years, pointing towards an environment-dependent expression. In the case of ascochyta blight resistance, the analysis which was based on a more saturated map yielded higher LOD scores and narrower confidence intervals for both QTLs (Af1 and Af2). The new marker OPZ08530 was the closest to LOD peak value of both Af1 QTLs identified in leaves and stems (Additional file 1: Table S1B). Addition of new markers in target areas of the map was an efficient method to increase the genome coverage and to obtain more defined QTLs.

Cross 29H × VF136 (RIL3)

The third map used in the present study was reported recently and includes 172 markers spanning 1402 cM [15]. The linkage groups were composed of 2 to 25 loci with a marker interval of 9.87 cM. Seven QTLs for O. crenata (Oc7 to Oc13) and 3 QTLs for O. foetida (Of3 to Of5) were identified in this map. Oc7 was detected along three years, explaining between 22% and 33% of the phenotypic variation. It has been suggested that Oc2 (previously reported in RIL2) and Oc7, which are both located in chr. VI and validated in different environments and genetic backgrounds, might correspond to the same QTL region (Figure 1; Additional file 1: Table S1). The new analysis considered 205 marker loci segregating in this population, of which 25 were common to both other RIL populations, 25 to RIL1 only and 4 to RIL3 only (Table 3).

Consensus linkage map

Three sets of faba bean mapping data were used in the construction of a consensus map connecting information of 11 F2 populations, marker data of 3 RILs, as well as new markers genotyped in the present study (Tables 1 and 2). The number of individual marker loci ranged from 313 in RIL1 to 404 in RIL2 and 205 in RIL3 (Table 2). Chi-square test was performed on new marker genotyping data for individual mapping population,s to test the null hypothesis of segregation ratios of 1:1. A variable percentage of distorted markers (P < 0.01) was observed, ranging from 2.93% in RIL3 to 7.35% in RIL1. A list of the marker loci is provided in Additional file 2: Table S2. QTL regions characterized in previous studies were also covered by the present consensus map.

The number of anchor markers for pairwise comparisons was initially small, with less than 20 markers in common. However, after the new marker analysis the number triplicated to 69 (2 isozymes, 4 SSRs, 2 legumines, 4 RGAs, 1 DR gene and 56 ITAPs). Sixty five of these, present in RIL1, were primarily used as bridges to integrate the individual maps into a single consensus map (Table 3). Twenty five markers were common in the 3 maps, while the remaining 44 were shared by at least 2 mapping populations. As a result, marker segregation data were assembled for a total of 828 marker loci, 759 of which were unique i.e. mapped only in one population (Table 5; Additional file 2: Table S2).

Table 5 Composite map of faba bean ( Vicia faba L.) genome

A total of 729 marker loci were assembled into the 43 LGs constituting the consensus map (Figure 1; Table 5), while 99 markers remained unlinked. The 6 major LGs contained between 53 (chr. V) and 165 marker loci (chr. I), and were assigned to the corresponding chromosomes. Seven additional LGs (Ia to IVa) could also be assigned thanks to the presence of loci previously located in individual chromosomes. Thirty one LGs consisted of 2–5 markers, and the remaining 6 LGs contained between 6 and 21 loci (Table 5; Additional file 2: Table S2). The total length of the consensus genetic linkage map was 4,613 cM, of which 3,442 cM were covered by the 6 main LGs/chromosomes. The length of these major LGs ranged from 323 cM (chr. V) to 1041 cM (the large metacentric chr. I). The entire consensus map had an average marker density of one marker per 10.7 cM, which was reduced to 6 cM when considering only the 6 main LGs. The marker order of the integrated map was largely collinear with the three individual maps, although a few local inversions and marker rearrangements over short intervals were observed.

Integration of QTL information

The number of QTL studies in faba bean is relatively low compared to other major legume species. Most traits have been genetically characterized in only one or two different mapping experiments, which limits the meta-analysis of QTLs in this species. Moreover, QTL intervals did not always include the minimum of two anchor markers, which is required for their projection onto the consensus map. Nevertheless, by comparing the maps published to date we provide a synthetic view of the most relevant loci controlling polygenic traits in faba bean. Further mapping of common markers between maps will be crucial to enhance the comparison of QTL positions from different mapping studies and to refine the localization of hot-spot genomic regions.

The 5 faba bean mapping experiments in the 3 RIL populations reported so far, identified 37 QTLs for 9 traits (Additional file 1: Table S1). Most of the QTL analyses focused on biotic stresses (e.g. broomrape and ascochyta blight resistance). The number of QTLs for broomrape resistance is 15 for O. crenata and 5 for O. foetida. Meanwhile the A. fabae resistance QTLs were reduced to 2, Af1 and Af2, that were conserved among populations (F2 and RIL) and environments. Regarding flowering and yield related traits, the number of stable QTLs reported were 1 (NOP), 2 (FL), 4 (PL), 3 (NSP) and 5 (DF) (Additional file 1: Table S1). The distribution of these 37 QTLs varied from 9 in chr. I, 8 in chr. II, 4 in chr. III and VI, 3 in chr. IV, 6 in chr. V. The last three QTLs remained unassigned.


During the last decade, significant progress was made in the development of genotyping tools. This allowed the addition of a large number of robust and transferrable marker loci in the genetic maps of relevant crop species. In faba bean, mapping studies were initiated in the 90's with the development of the first maps in F2 populations using mostly RAPDs together with SSRs, isozymes and morphological markers. Previous studies comparing these linkage maps have been reported [9]. The use of a recurrent parent (Vf6) in all the F2 populations then allowed to join data from different progenies by means of common markers in the female parent. Moreover, the use of trisomic families for chromosomes III, IV, V and VI allowed allocation of LGs to chromosomes [9]. Ever since, attempts have been made to increase marker density using new SSRs and gene-based markers in RIL populations. The main objective of the present study was the development of a high density consensus genetic map that integrates all the relevant maps reported so far and serves as a reference map for the international faba bean community.

Building a consensus map is not possible without common or bridge loci on each LG or chromosome. For this reason, a number of additional markers was genotyped in each mapping population to increase the number of common markers among them. A bridge marker was considered as such when its name and position were the same in the different mapping populations. The genetic map was created combining two approaches: (a) increasing the number of anchor marker loci in the different populations, (b) merging the resulting genetic maps through markers common to three RIL populations with MergeMap [42] as reported in many other crop species [39, 4347].

Using this approach, segregation data for 729 marker were assembled on 43 LGs. In terms of marker order, the consensus map contains few changes compared to the individual maps. Small discrepancies in the marker order or position in some LGs might be due to (i) different population sizes used, (ii) weak linkages existing in the different maps, or (iii) missing or poor quality data, rather than to real chromosome rearrangements. As reported in previous studies in Vitis vinifera L. [47], phaseolus [37] or Brassica napus [44], the faba bean consensus marker order is significantly more reliable than that of the individual maps because of the higher number of individuals and recombination events occuring across three or more populations.

Based on previous LG/marker allocations, 13 of the larger LGs could be assigned to specific chromosomes while 30 LGs remained unassigned. Considering the enormous size of the faba bean genome, unassigned LGs may be due to recombination gaps at the distal ends of the main LGs because of a lack of marker loci. None of the main LGs differed considerably in marker density. The length of our core map was 3,431 cM, which is higher than the single RIL maps. In many other species the increased size of the composite map was attributed to an improved coverage of the chromosome ends [37, 4850].

All the individual maps reported to date allocate LGs to 5 of the 6 faba bean chromosomes, excluding chr. IV. After acknowledging the erroneous assignment of LG I.B, which actually corresponds to chr. IV [51], the new integrated map anchors for the first time the main LGs to the whole chromosome complement of the species. This information was used here to update the reported large-scale synteny between LGs and/or chromosomes of M. truncatula and cool season grain legumes such as pea, chickpea, lens and faba bean [52]. Figure 2 shows the main syntenic blocks and rearrangements among these species and their correspondence to the six faba bean chromosomes.

Figure 2

Schematic representation of large-scale synteny blocks between chromosomes and chromosome segments of M. truncatula (Mt) and main cool season grain legumes (Source [52], with modifications). Chickpea (Cicer arietinum; Ca), faba bean (Vicia faba; Vf), lens (Lens culinaris; Lc), and pea (Pisum sativum; Ps). Bars representing Medicago and pea homologous chromosomal regions are shown with the same gray intensity or pattern. Arrows in the boxes indicate the orientation of the chromosomes (short arm - long arm) in the case of Medicago. The corresponding synteny blocks of faba bean, chickpea and lentil are represented by blank bars. The bars do not reflect the relative sizes of chromosome or chromosome segments and the break points of chromosomes are indicated approximately. The figure integrates data from [11, 13, 5358] and this study.

The integration of a high number of coding regions in the consensus map provides an excellent framework for downstream analyses, including comparisons between the locations of major genes for important traits or QTL positions between populations from different crosses. Moreover, combining syntenic studies with a consensus map will contribute to increase marker density in genomic regions of interest for indirect selection or for map based cloning [37]. Successful application of consensus maps for synteny based candidate gene identification and/or definition of QTL location has been extensively used both in cereals [5962] and in legumes [35, 37, 63, 64].

The RIL mapping populations considered in the consensus map were used previously for detecting QTLs of agronomically important traits. These are displayed in Figure 1, together with the QTLs identified in the improved Vf6 × Vf136 map and detected in the present study. In order to increase the density of loci around the OTLs for broomrape and ascochyta blight resistance, we used BSA in contrasted DNA pools. The BSA approach has been applied in numerous studies and provides a platform for high-resolution genetic analysis [65]. In the present study, 24 of 31 RAPD markers were mapped to the major linkage groups and allowed more accurate determination of QTL locations and effects. These results highlight the usefulness of BSA based on markers flanking QTLs, as an efficient tool for saturation of targeted regions, opening the possibility of future marker-assisted selection for these traits.

Faba bean has been considered a “genomic orphan” crop with a huge and complex genome and limited availability of genetic and genomic resources. At present, the situation has greatly improved thanks to the technological advances in high-throughput sequencing and genotyping, together with the access to genomic and transcriptomic tools. Genome-wide transcription profiling by deepSuperSAGE was recently used for quantifying the transcriptional changes elicited by A. fabae and to identify candidate resistance genes governing faba bean responses to this fungal pathogen [66]. Several genome libraries have been constructed and characterized for putative SSR sequences using the Roche 454 GS FLX Titanium Sequencing Platform [67, 68]. These transcriptomic studies provide a foundation for the identification of novel regulators associated with faba bean-pathogen interactions and also a valuable source of markers for molecular breeding applications in this crop.

Translation of genomic resources from the model species M. truncatula or other sequenced related legume species such as chickpea, should be further exploited to raise the prospects in molecular faba bean breeding programs. The availability of large sets of conserved ESTs from model or related species constitutes a valuable source of markers that are physically associated with coding regions. These are good candidates for gene cloning or faba bean marker assisted selection. This is the approach used in this study to integrate all the faba bean genomic information so far reported, and to generate a new tool of reference for faba bean breeding and genomics approaches.


We have constructed the first marker consensus genetic linkage map for faba bean by integrating segregation data from three recombinant inbred line populations, together with new common markers genotyped in this study. The final integrated map has allowed to (i) join a larger number of markers than in any previous individual map, (ii) obtain a more complete coverage of the faba bean genome, (iii) fill a number of gaps in previous independent maps, and (iv) improve the resolution of key QTLs. The colinearity of the consensus map was well maintained and will serve as reference for future faba bean multiple-line cross QTL mapping studies. Since 60% of the markers in the most developed map (RIL1) corresponded to coding regions, this consensus enhanced-density faba bean map provides a functional framework for candidate gene studies, expression analysis, comparative genomics, evolution studies and anchoring of the future faba bean genome sequences.


Mapping populations

The most recent maps of three RIL mapping populations,Vf6 × Vf27 (RIL1), Vf6 × Vf136 (RIL 2) and 29H × Vf136 (RIL3), were used to develop an integrated map of faba bean (Table 1). Vf6 was a common female parent in two progenies, Vf6 × VF136, which segregates for broomrape and Ascochyta resistance [7, 8, 2426] and Vf6 × VF27, first reported by [6] and further used to construct the first exclusively gene-based genetic map in the species [11] and to identify and validate QTLs controlling flowering time and other yield-related traits [13]. Vf136 was the common male parent with the third population, 29H × Vf136, segregating for resistance to broomrape and A. fabae [10, 15]. The populations consisted of 124 RILs for RIL1, 165 individual lines for RIL2, and 119 for RIL3.

The female parents, Vf6 and 29H are equina medium-seeded field beans with beige seed coat and resistance to A. fabae, Vf136 is also an equina type with reported resistance to broomrape and Vf27 is a black and small-seeded paucijuga form, supposedly close to a putative wild faba bean progenitor [69].

Marker analyses

Genomic DNA was extracted from young leaves using liquid nitrogen and the procedure was as described by [3]. To increase marker density and to provide common markers to anchor the LGs from different populations, new markers from different sources were assayed. A set of SSRs, resistant gene analogs (RGAs), defence delated (DR) genes and ITAPs designed from different legume ESTs (M. truncatula, pea, lentil, lupin and soybean), were tested in the parental lines and the polymorphic ones genotyped in the corresponding RIL population.

SSR markers

Fifty four faba bean SSRs [70] and 41 pea SSRs [71], were assayed for polymorphism among the parental lines using their respective protocols. In case of pea SSRs the reaction mixture was modified slightly by using 2,5 mM of MgCl2 and 1U Taq polymerase instead, to facilitate the orthologous amplification. SSRs revealing consistent and easily scorable bands were genotyped in the whole populations after electrophoresis in 2.5% - 3% agarose gels.

RGAs and DR genes

Ten RGA classes were tested using PCR conditions described by [72]. To reveal polymorphism, amplification products for each RGA class were digested with a set of restriction enzymes according to the manufacturer’s instructions to obtain CAPS (Cleaved Amplified Polymorphic Sequences). Twelve additional RGAs [73, 74] along with 12 DR genes, cloned and mapped in different legume species and mapped in pea [74], were also assayed. Amplifications and PCR conditions were as described by [72].

Intron-targeted amplified polymorphic markers (ITAPs)

A total of 635 EST derived markers developed within the Grain Legumes Integrated Project (GLIP-Food-CT-2004-506223), were tested for polymorphism among the parental lines using the amplification protocols reported by [13, 26]. Special efforts were focused on genotyping the ITAPs previously mapped in the most advanced Vf6 × Vf27 map [11, 13]. As mentioned above, when no polymorphism was detected on agarose gels, PCR products amplified from both parents were digested with a range of restriction endonucleases in order to detect a SNP as a CAPS (Cleaved Amplified Polymorphic Sequence) marker that was further genotyped in the corresponding population.

Saturation mapping

In order to saturate targeted regions conferring broomrape or ascochyta blight resistance we applied the BSA [75] based on previous QTL mapping information in cross Vf6 × Vf136 [24, 26]. BSA has been widely adopted as a method to rapidly identify molecular makers in specific genome regions. The BSA principle consists in pooling DNAs from individuals from a segregating population according to two phenotypic classes. The resulting DNA bulks are equivalent to those from two Near Isogenic Lines (NILs) for which is assumed to generate a random genetic background at all other unlinked loci. In this study segregating individuals were grouped according to the genotype of markers flanking already localized QTLs. The contrasting pools were then screened with new markers in order to identify recombinants within each QTL interval.

QTLs underlying resistance to O. crenata and A. fabae were named Oc and Af, respectively [24, 26]. Accordingly, bulks of plants fixed for alleles of the two markers flanking four QTL regions were selected from the RIL population: OPN071409 and OPAI131018 (flanking Oc2 on linkage group VI.B), OPC191059 and OPD12425 (Oc3 on LG II.A), OPF08710 and OPW15576 (Af1 on chr. III), OPAG5737 and MER021469 (Af2 on chr. II). A total of 748 RAPD primers was used in search for polymorphisms between the two bulks. For a given target region, markers showing expected differences between the pair of bulks were subsequently used to screen 14 plants individually. Markers that maintained the expected pattern of polymorphisms, were then used to screen the entire RIL population.

Quantitative traits

Traits and QTL information were selected from seven published works [7, 8, 10, 13, 15, 24, 26] and supplemented by the bulked segregant analysis (in cross VF6 × VF136) and the saturation process described above. Trait descriptions, evaluation methods and abbreviations were assigned according to the previous references. Thus, the nine traits considered (Figure 1; Additional file 2), were the following: (1) Orobanche crenata resistance (trait abbreviation OC); (2) Orobanche foetida resistance (OF); (3) Ascochyta fabae resistance: disease severity on leaves (DSL); (4) Ascochyta faba resistance: disease severity on stems (DSS); (5) Days to flowering (DF); (6) Flowering length (FL); (7) Pod length (PL); (8) Number of ovules per pod (NOP); (9) Number of seeds per pod (NSP). Most of these QTLs showed to be stable as were identified and validated in different environments or genetic background. Each QTL was treated independently, making it possible to notice the number of times that a QTL is reported in a similar genomic location across independent experiments (Figure 1; Additional file 2).

Data analysis

Map construction and QTL analysis in Vf6 × VF136 (RIL 2)

MAPMAKER 3.0 [76] was used to identify linkage groups using an LOD score of four as the threshold for considering significant linkage. MSTMap software [77] was used to determine maker orders by finding the minimum spanning tree of a graph for each linkage group. MAPMAKER was used to confirm marker orders determined by MSTMap and to convert the recombination fractions to centiMorgans (cM) using the mapping function of Kosambi [78].

QTL analysis was conducted using composite interval mapping (CIM) and multiple interval mapping (MIM) in Windows QTL Cartographer V2.5 [79]. Markers to be used as cofactors for CIM were selected by forward-backward stepwise regression. The number of markers controlling the genetic background in CIM was set to five. The thresholds for the detection of QTLs were estimated by permutations analysis [80] using 1,000 permutations.

Consensus map construction

Three individual genetic maps were used to generate a consensus map using MergeMap [42] by converting the individual maps into directed acyclic graphs (DAGs) that are then merged in consensus graph on the basis of their shared vertices [81]. As MergeMap tends to inflate genetic distances in the consensus genetic map [42, 82], marker data from different mapping populations were pooled together and the order of each consensus linkage group as established by MergeMap was set, in order to calculate genetic distances using MAPMAKER. The consensus map for each linkage group was visualized by MapChart [83].


  1. 1.

    Flavell RB, Bennett MD, Smith JB, Smith DB: Genome size and proportion of repeated nucleotide sequence DNA in plants. Biochem Genet. 1974, 12: 257-269. 10.1007/BF00485947.

    Article  CAS  PubMed  Google Scholar 

  2. 2.

    van de Ven WTG, Waugh R, Duncan N, Ramsay G, Dow N, Powell W: Development of a genetic linkage map in Vicia faba using molecular and biochemical techniques. Aspects Appl Biol. 1991, 27: 49-54.

    Google Scholar 

  3. 3.

    Torres AM, Weeden NF, Martín A: Linkage among isozyme, RFLP and RAPD markers in Vicia faba. Theor Appl Genet. 1993, 85: 937-945.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Ramsay G, van de Ven W, Waugh R, Griffiths DW, Powel W: Mapping quantitative trait loci in faba beans. Improving production and utilisation of grain legumes. Edited by: AEP. 1995, Copenhagen, Denmark: 2nd European Conference on Grain Legumes, 444-445.

    Google Scholar 

  5. 5.

    Satovic Z, Torres AM, Cubero JI: Genetic mapping of new morphological, isozyme and RAPD markers in Vicia faba L. using trisomics. Theor Appl Genet. 1996, 93: 1130-1138. 10.1007/BF00230136.

    Article  CAS  PubMed  Google Scholar 

  6. 6.

    Vaz Patto MC, Torres AM, Koblizkova A, Macas J, Cubero JI: Development of a genetic composite map of Vicia faba using F2 populations derived from trisomic plants. Theor Appl Genet. 1999, 98: 736-743. 10.1007/s001220051129.

    Article  CAS  Google Scholar 

  7. 7.

    Román B, Torres AM, Rubiales D, Cubero JI, Satovic Z: Mapping of quantitative trait loci controlling broomrape (Orobanche crenata Forsk.) resistance in faba bean (Vicia faba L.). Genome. 2002, 45: 1057-1063. 10.1139/g02-082.

    Article  PubMed  Google Scholar 

  8. 8.

    Román B, Satovic Z, Avila CM, Rubiales D, Moreno MT, Torres AM: Locating genes associated with Ascochyta fabae resistance in Vicia faba L. Aust J Agric Res. 2003, 54: 85-90. 10.1071/AR02034.

    Article  Google Scholar 

  9. 9.

    Román B, Satovic Z, Pozarkova D, Macas J, Dolezel J, Cubero JI, Torres AM: Development of a composite map in Vicia faba L. breeding applications and future prospects. Theor Appl Genet. 2004, 108: 1079-1088. 10.1007/s00122-003-1515-6.

    Article  PubMed  Google Scholar 

  10. 10.

    Avila CM, Satovic Z, Sillero JC, Rubiales D, Moreno MT, Torres AM: Isolate and organ-specific QTLs for ascochyta blight resistance in faba bean. Theor Appl Genet. 2004, 108: 1071-1078. 10.1007/s00122-003-1514-7.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Ellwood SR, Phan HTT, Jordan M, Torres AM, Avila CM, Cruz-Izquierdo S, Oliver RP: Construction of a comparative genetic map in faba bean (Vicia faba L.); conservation of genome structure with Lens culinaris. BMC Genomics. 2008, 9: 380-10.1186/1471-2164-9-380.

    PubMed Central  Article  PubMed  Google Scholar 

  12. 12.

    Arbaoui M, Link W, Satovic Z, Torres AM: Quantitative trait loci of frost tolerance and physiologically related traits in faba bean (Vicia faba L.). Euphytica. 2008, 164: 93-104. 10.1007/s10681-008-9654-0.

    Article  CAS  Google Scholar 

  13. 13.

    Cruz-Izquierdo S, Avila CM, Satovic Z, Palomino C, Gutierrez N, Ellwood SR, Phan HTT, Cubero JI, Torres AM: Comparative genomics to bridge Vicia faba with model and closely-related legume species: stability of QTLs for flowering and yield-related traits. Theor Appl Genet. 2012, 125: 1767-1782. 10.1007/s00122-012-1952-1.

    Article  CAS  PubMed  Google Scholar 

  14. 14.

    Ma Y, Bao S, Yang T, Hu J, Guan J, He Y, Wang X, Wan Y, Sun X, Jiang J, Gong C, Zong X: Genetic linkage map of Chinese native variety faba bean (Vicia faba L.) based on simple sequence repeat markers. Plant Breeding. 2013, doi:10.1111/pbr.12074

    Google Scholar 

  15. 15.

    Gutiérrez N, Palomino C, Satovic Z, Ruiz-Rodríguez MD, Vitale S, Gutiérrez MV, Rubiales D, Kharrat M, Amri M, Emeran A, Cubero JI, Atienza SG, Torres AM, Avila CM: QTLs for Orobanche spp. resistance in faba bean: Identification and validation across different environments. Molecular Breeding. 2013, in press

    Google Scholar 

  16. 16.

    Cabrera A, Cubero JI, Martín A: Genetic mapping using trisomics in Vicia faba L. FABIS Newslett. 1989, 23: 5-7.

    Google Scholar 

  17. 17.

    Macas J, Dolezel J, Lucretti S, Pich U, Meister A, Fuchs J, Schubert I: Localization of seed genes on flow-sorted field bean chromosomes. Chromosome Res. 1993, 1: 107-115. 10.1007/BF00710033.

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Macas J, Weschke W, Baumlein H, Pich U, Houben A, Wobus U, Schubert I: Localization of vicilin genes via polymerase chain reaction on microisolated field bean chromosomes. Plant J. 1993, 3: 883-886. 10.1111/j.1365-313X.1993.00883.x.

    Article  CAS  PubMed  Google Scholar 

  19. 19.

    Fuchs J, Pich U, Meister A, Schubert I: Differentiation of field bean heterochromatin by in situ hybridization with a repeated FokI sequence. Chromosome Res. 1994, 2: 25-28. 10.1007/BF01539450.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Pozarkova D, Koblizkova A, Román B, Torres AM, Lucretti S, Lysak M, Dolezel J, Macas J: Development and characterization of microsatellite markers from chromosome 1-specific DNA libraries of Vicia faba. Biol Plantarum. 2002, 45: 337-345. 10.1023/A:1016253214182.

    Article  CAS  Google Scholar 

  21. 21.

    Torres AM, Satovic Z, Cánovas J, Cobos S, Cubero JI: Genetics and mapping of new isozyme loci in Vicia faba L. using trisomics. Theor Appl Genet. 1995, 91: 783-789.

    Article  CAS  PubMed  Google Scholar 

  22. 22.

    Torres AM, Vaz Patto MC, Satovic Z, Cubero JI: New isozyme loci in Faba bean (Vicia faba L.): Genetic analysis and mapping using trisomics. J Hered. 1998, 89: 271-274. 10.1093/jhered/89.3.271.

    Article  CAS  Google Scholar 

  23. 23.

    Avila CM, Sillero JC, Rubiales D, Moreno MT, Torres AM: Identification of RAPD markers linked to Uvf-1 gene conferring hypersensitive resistance against rust (Uromyces viciae-fabae) in Vicia faba L. Theor Appl Genet. 2003, 107: 353-358. 10.1007/s00122-003-1254-8.

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Díaz-Ruiz R, Satovic Z, Avila CM, Alfaro CM, Gutierrez MV, Torres AM, Román B: Confirmation of QTLs controlling Ascochyta fabae resistance in different generations of faba bean (Vicia faba L.). Crop Pasture Sci. 2009, 60: 353-361. 10.1071/CP08190.

    Article  Google Scholar 

  25. 25.

    Díaz-Ruiz R, Torres A, Gutierrez MV, Rubiales D, Cubero JI, Kharrat M, Satovic Z, Román B: Mapping of quantitative trait loci controlling Orobanche foetida Poir. resistance in faba bean (Vicia faba L.). African J Biotechnol. 2009, 8: 2718-2724.

    Google Scholar 

  26. 26.

    Díaz R, Torres AM, Satovic Z, Gutierrez MV, Cubero JI, Román B: Validation of QTLs for Orobanche crenata resistance in faba bean (Vicia faba L.) across environments and generations. Theor Appl Genet. 2010, 120: 909-919. 10.1007/s00122-009-1220-1.

    Article  Google Scholar 

  27. 27.

    Truco MJ, Antonise R, Lavelle D, Ochoa O, Kozik A, Witsenboer H, Fort S, Jeuken MJW, Kesseli RV, Lindhout P, Michelmore RW: A high-density, integrated genetic linkage map of lettuce (Lactuca spp.). Theor Appl Genet. 2007, 115: 735-746. 10.1007/s00122-007-0599-9.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Varshney RK, Marcel TC, Ramsay L, Russell J, Röder MS, Stein N, Waugh R, Langridge P, Niks RE, Graner A: A high density barley microsatellite consensus map with 775 SSR loci. Theor Appl Genet. 2007, 114 (6): 1091-103. 10.1007/s00122-007-0503-7.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Shirasawa K, Bertioli DJ, Varshney RK, Moretzsohn MC, Leal-Bertioli SC, Thudi M, Pandey MK, Rami JF, Foncéka D, Gowda MV, Qin H, Guo B, Hong Y, Liang X, Hirakawa H, Tabata S, Isobe S: Integrated consensus map of cultivated peanut and wild relatives reveals structures of the A and B genomes of Arachis and divergence of the legume genomes. DNA Res. 2013, 20 (2): 173-184. 10.1093/dnares/dss042. doi:10.1093/dnares/dss042

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  30. 30.

    Somers JD, Isaac P, Edwards K: A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theor Appl Genet. 2004, 109: 1105-1114. 10.1007/s00122-004-1740-7.

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Davis GL, McMullen MD, Baysdorfer C, Musket T, Grant D, Staebell M, Xu G, Polacco M, Koster L, Melia-Hancock S, Houchins K, Chao S, Coe EH: A maize map standard with sequenced core markers, grass genome reference points and 932 expressed sequence tagged sites (ESTs) in a 1736-locus map. Genetics. 1999, 152: 1137-1172.

    PubMed Central  CAS  PubMed  Google Scholar 

  32. 32.

    Close TJ, Bhat PR, Lonardi S, Wu Y, Rostoks N, Ramsay L, Druka A, Stein N, Svensson JT, Wanamaker S, Bozdag S, Roose ML, Moscou MJ, Chao S, Varshney RK, Szucs P, Sato K, Hayes PM, Matthews DE, Kleinhofs A, Muehlbauer GJ, DeYoung J, Marshall DF, Madishetty K, Fenton RD, Condamine P, Graner A, Waugh R: Development and implementation of high-throughput SNP genotyping in barley. BMC Genomics. 2009, 10: 582-10.1186/1471-2164-10-582.

    PubMed Central  Article  PubMed  Google Scholar 

  33. 33.

    Antonio BA, Inoue T, Kajiya H, Nagamura Y, Kurata N, Minobe Y, Yano M, Nakagahra M, Sasaki T: Comparison of genetic distance and order of DNA markers in five populations of rice. Genome. 1996, 39: 946-956. 10.1139/g96-119.

    Article  CAS  PubMed  Google Scholar 

  34. 34.

    Hwang T-Y, Sayama T, Takahashi M, Takada Y, Nakamoto Y, Funatsuki H, Hisano H, Sasamoto S, Sato S, Tabata S, Kono I, Hoshi M, Hanawa M, Yano C, Xia Z, Harada K, Kitamura K, Ishimoto M: High-density integrated linkage map based on SSR markers in soybean. DNA Res. 2009, 16: 213-225. 10.1093/dnares/dsp010. doi:10.1093/dnares/dsp010

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  35. 35.

    Bordat A, Savois V, Nicolas M, Salse J, Chauveau A, Bourgeois M, Potier J, Houtin H, Rond C, Murat F: Translational genomics in legumes allowed placing in silico 5460 unigenes on the Pea functional Map and identified candidate genes in Pisum sativum L. G3-Genes genomes. Genetics. 2011, 13 (2): 93-103.

    Google Scholar 

  36. 36.

    Millan T, Winter P, Jüngling R, Gil J, Rubio J, Cho S, Cobos MJ, Iruela M, Rajesh PN, Tekeoglu M, Kahl G, Muehlbauer FJ: A consensus genetic map of chickpea (Cicer arietinum L.) based on 10 mapping populations. Euphytica. 2010, 175: 175-189. 10.1007/s10681-010-0157-4.

    Article  CAS  Google Scholar 

  37. 37.

    Galeano CH, Fernandez AC, Franco N, Cichy K, McClean P, 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): e28135-10.1371/journal.pone.0028135. doi:10.1371/journal.pone.0028135

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  38. 38.

    Bohra A, Saxena RK, Gnanesh BN, Kulbhushan Saxena K, Byregowda M, Rathore A, KaviKishor PB, Cook DR, Varshney RK: An intra-specific consensus genetic map of pigeonpea [Cajanus cajan (L.) Millspaugh] derived from six mapping populations. Theor Appl Genet. 2012, 125 (6): 1325-1338. 10.1007/s00122-012-1916-5.

    PubMed Central  Article  PubMed  Google Scholar 

  39. 39.

    Muchero W, Diop NN, Bhat PR, Fenton RD, Wanamaker S, Pottorff M, Hearne S, Cisse N, Fatokun C, Ehlers JD, Roberts PA, Close TJ: A consensus genetic map of cowpea [Vigna unguiculata (L) Walp.] and synteny based on EST-derived SNPs. Proc Natl Acad Sci U S A. 2009, 106: 18159-18164. 10.1073/pnas.0905886106.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  40. 40.

    Gautami B, Foncéka D, Pandey MK, Moretzsohn MC, Sujay V, Qin H, Hong Y, Faye I, Chen X, Prakash AB, Shah TM, Gowda MVC, Nigam SN, Liang X, Hoisington DA, Guo B, Bertioli DJ, Ram JF, Varshney RK: An international reference consensus genetic map with 897 marker loci based on 11 mapping populations for tetraploid groundnut (Arachis hypogaea L.). PLoS One. 2012, 7 (7): e41213-10.1371/journal.pone.0041213. doi:10.1371/journal.pone.0041213. 2012

    PubMed Central  Article  PubMed  Google Scholar 

  41. 41.

    Isobe S, Kölliker R, Hisano H, Sasamoto S, Wada T, Klimenko I, Okumura K, Tabata S: Construction of a consensus linkage map for red clover (Trifolium pratense L). BMC Plant Biol. 2009, 9: 57-10.1186/1471-2229-9-57. doi:10.1186/1471-2229-9-57

    PubMed Central  Article  PubMed  Google Scholar 

  42. 42.

    Wu Y, Close TJ, Lonardi S: Accurate construction of consensus genetic maps via integer linear programming. IEEE/ACM Trans. Comput Biol Bioinform. 2011, 8: 381-394. doi:10.1109/TCBB.2010.35

    Google Scholar 

  43. 43.

    Rajaram E, Nepolean T, Senthilvel S, Varshney RK, Vadez V, Srivastava RK, Shah TM, Supriya A, Kumar S, Kumari BR, Bhanuprakash A, Narasu ML, Riera-Lizarazu O, Hash CH: Pearl millet [Pennisetum glaucum (L.) R. Br.] consensus linkage map constructed using four RIL mapping populations and newly developed EST-SSRs. BMC Genomics. 2013, 14: 159-10.1186/1471-2164-14-159. doi:10.1186/1471-2164-14-15

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  44. 44.

    Wang J, Lydiate D, Parkin I, Falentin C, Delourme R, Carion P, King GJ: Integration of linkage maps for the Amphidiploid Brassica napus and comparative mapping with Arabidopsis and Brassica rapa. BMC Genomics. 2011, 12: 101-10.1186/1471-2164-12-101.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  45. 45.

    Khan MA, Han Y, Zhao YF, Troggio M, Korban SS: A multi-population consensus genetic map reveals inconsistent marker order among maps likely attributed to structural variations in the apple genome. PLoS One. 2012, 7 (11): e47864-10.1371/journal.pone.0047864. doi:10.1371/journal.pone.0047864

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  46. 46.

    Han Y, Zheng D, Vimolmangkang S, Khan MA, Beever JE, Korban SS: Integration of physical and genetic maps in apple confirms whole-genome and segmental duplications in the apple genome. J Exp Bot. 2011, 62: 5117-5130. 10.1093/jxb/err215. doi:10.1093/jxb/err215

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  47. 47.

    Vezzulli S, Troggio M, Coppola G, Jermakow A, Cartwright D, Zharkikh A, Stefanini M, Grando MS, Viola R, Adam-Blondon AF, Thomas M, This P, Velasco R: A reference integrated map for cultivated grapevine (Vitis vinifera L.) from three crosses, based on 283 SSR and 501 SNP-based markers. Theor Appl Genet. 2008, 117: 499-511. 10.1007/s00122-008-0794-3.

    Article  CAS  PubMed  Google Scholar 

  48. 48.

    Mace ES, Rami J-F, Bouchet S, Klein P, Klein RR, Kilian A, Wenzl P, Xia L, Halloran K, Jordan DR: A consensus genetic map of sorghum that integrates multiple component maps and high-throughput diversity array technology (DArT) markers. BMC Plant Biol. 2009, 9: 13-10.1186/1471-2229-9-13. doi:10.1186/1471-2229-9-13

    PubMed Central  Article  PubMed  Google Scholar 

  49. 49.

    Spiller M, Linde M, Hibrand-Saint Oyant L, Tsai C-J, Byrne DH, Smulders MJ, Foucher F, Debener T: Towards a unified genetic map for diploid roses. Theor Appl Genet. 2011, 122: 489-500. 10.1007/s00122-010-1463-x.

    Article  PubMed  Google Scholar 

  50. 50.

    Marone D, Laido’ G, Gadaleta A, Colasuonno P, Ficco DBM, Giancaspro A, Giove S, Panio G, Russo MA, De Vita P, Cattivelli L, Papa R, Blanco A, Mastrangelo AM: A high-density consensus map of A and B wheat genomes. Theor Appl Genet. 2012, 125: 1619-1638. 10.1007/s00122-012-1939-y. doi:10.1007/s00122-012-1939-y

    PubMed Central  Article  PubMed  Google Scholar 

  51. 51.

    Ruiz-Rodriguez MD, Avila CM, Torres AM, Fuchs J, Schubert I: Anchoring of genetic linkage maps to the chromosome complement of Vicia faba L. Mole breed. in press doi:10.1007/s11032-013-9969-y

  52. 52.

    Kaló P, Choi HK, Ellis N, Kiss GB: Synteny and Comparative Genomics between Model and Cool Season Grain Legumes. Genetics, Genomics and Breeding in Crop Plants: Cool Season Food Legumes. Edited by: de la Vega P, Torres AM, Cubero JI, Kole C. 2011, New Hampshire, Jersey, Plymouth: Science Pubs Inc, 285-302.

    Google Scholar 

  53. 53.

    Choi HK, Mun JH, Kim DJ, Zhu H, Baek JM, Mudge J, Roe B, Ellis N, Doyle J, Kiss GB, Young ND, Cook DR: Estimating genome conservation between crop and model legume species. Proc Natl Acad Sci U S A. 2004, 101: 15289-15294. 10.1073/pnas.0402251101.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  54. 54.

    Kaló P, Seres A, Taylor SA, Jakab J, Kevei Z, Kereszt A, Endre G, Ellis THN, Kiss GB: Comparative mapping between Medicago sativa and Pisum sativum. Mol Genet Genom. 2004, 272: 235-246. 10.1007/s00438-004-1055-z.

    Article  Google Scholar 

  55. 55.

    Zhu H, Choi HK, Cook DR, Shoemaker RC: Bridging model and crop legumes through comparative genomics. Plant Physiol. 2005, 137: 1189-1196. 10.1104/pp.104.058891.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  56. 56.

    Nelson MN, Phan HTT, Ellwood SR, Moolhuijzen PM, Hane J, Williams A, O‘Lone CE, Fosu-Nyarko J, Scobie M, Cakir M, Jones MGK, Bellgard M, Ksiarkiewicz M, Wolko B, Barker SJ, Oliver RP, Cowling WA: The first gene-based map of Lupinus angustifolius L. -location of domestication genes and conserved synteny with Medicago truncatula. Theor Appl Genet. 2006, 113: 225-238. 10.1007/s00122-006-0288-0.

    Article  CAS  PubMed  Google Scholar 

  57. 57.

    Phan HTT, Ellwood SR, Ford R, Thomas S, Oliver R: Differences in syntenic complexity between Medicago truncatula with Lens culinaris and Lupinus albus. Funct Plant Biol. 2006, 33: 775-782. 10.1071/FP06102.

    Article  CAS  Google Scholar 

  58. 58.

    Nayak SN, Zhu HY, Varghese N: Integration of novel SSR and gene-based SNP marker loci in the chickpea genetic map and establishment of new anchor points with Medicago truncatula genome. Theor Appl Genet. 2010, 120: 1415-1441. 10.1007/s00122-010-1265-1.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  59. 59.

    Marcel TC, Varshney RK, Barbieri M, Jafary H, Kock MJD, Graner A, Niks RE: A high-density consensus map of barley to compare the distribution of QTLs for partial resistance to Puccinia hordei and of defence gene homologues. Theor Appl Genet. 2007, 114: 487-500. 10.1007/s00122-006-0448-2. doi:10.1007

    Article  CAS  PubMed  Google Scholar 

  60. 60.

    Cockram J, White J, Zuluaga DL, Smith D, Comadran J, Macaulay M, Luo Z, Kearsey MJ, Werner P, Harrap D, Tapsell C, Liu H, Hedley PE, Stein N, Schulte D, Steuernagel B, Marshall DF, Thomas WTB, Ramsay L, Mackay I, Balding DJ, Waugh R, O’Sullivan DM, The AGOUEB Consortium: Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome. Proc Natl Acad Sci U S A. 2010, 107 (50): 21611-21616. 10.1073/pnas.1010179107. doi:10.1073/pnas.1010179107

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  61. 61.

    Sehgal D, Rajaram V, Vadez V, Hash CT, Yadav RS: Integration of gene-based markers in pearl millet genetic map for identification of candidate genes underlying drought tolerance quantitative trait loci. BMC Plant Biol. 2012, 12: 9-10.1186/1471-2229-12-9.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  62. 62.

    Swamy BP, Vikram P, Dixit S, Ahmed HU, Kumar A: Meta-analysis of grain yield QTL identified during agricultural drought in grasses showed consensus. BMC Genomics. 2011, 2: 319-doi:10.1186/1471-2164-12-319

    Article  Google Scholar 

  63. 63.

    Muchero W, Ehlers JD, Close TJ, Roberts PA: Genic SNP markers and legume synteny reveal candidate genes underlying QTL for Macrophomina phaseolina resistance and maturity in cowpea [Vigna unguiculata (L) Walp.]. BMC Genomics. 2011, 12: 8-10.1186/1471-2164-12-8.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  64. 64.

    Pottorff M, Ehlers JD, Christian Fatokun, JD, Philip AR, Timothy JC: Leaf morphology in Cowpea [Vigna unguiculata (L.) Walp]: QTL analysis, physical mapping and identifying a candidate gene using synteny with model legume species. BMC Genomics. 2012, 13: 234-10.1186/1471-2164-13-234.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  65. 65.

    Chen X, Hedley PE, Morris J, Liu H, Niks RE, Waugh R: Combining genetical genomics and bulked segregant analysis-based differential expression: an approach to gene localization. Theor App Gen. 2011, 122: 1375-1383. 10.1007/s00122-011-1538-3.

    Article  Google Scholar 

  66. 66.

    Madrid E, Palomino C, Plötner A, Horres R, Jüngling R, Rotter B, Winter P, Kahl G, Torres AM: DeepSuperSage analysis of the Vicia faba transcriptome in response to Ascochyta fabae infection. Phytopathologia Mediterranea. 2013, 52 (1): 166-182.

    CAS  Google Scholar 

  67. 67.

    Yang T, Bao S, Ford R, Jia T, Guan J, He Y, Sun X, Jiang J, Hao J, Zhang X, Zong X: High-throughput novel microsatellite marker of faba bean via next generation sequencing. BMC Genomics. 2012, 13: 602-10.1186/1471-2164-13-602.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  68. 68.

    Kaur S, Pembleton L, Cogan N, Savin K, Leonforte T, Paull J, Materne M, Forster J: Transcriptome sequencing of field pea and faba bean for discovery and validation of SSR genetic markers. BMC Genomics. 2012, 13: 104-10.1186/1471-2164-13-104.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  69. 69.

    Cubero JI: On the evolution of Vicia faba L. Theor Appl Genet. 1974, 45: 47-51.

    Article  CAS  PubMed  Google Scholar 

  70. 70.

    Zeid M, Mitchell S, Link W, Carter M, Nawar A, Fulton T, Kresovich S: Simple sequence repeats (SSRs) in faba bean: new loci from Orobanche-resistant cultivar ‘Giza 402’. Plant Breeding. 2009, 128: 149-155. 10.1111/j.1439-0523.2008.01584.x.

    Article  CAS  Google Scholar 

  71. 71.

    Loridon K, McPhee K, Morin J, Dubreuil P, Pilet-Nayel ML, Aubert G, Rameau C, Baranger A, Coyne C, Lejeune-Henaut I, Burstin J: Microsatellite marker polymorphism and mapping in pea (Pisum sativum L.). Theor Appl Genet. 2005, 111: 1022-103. 10.1007/s00122-005-0014-3.

    Article  CAS  PubMed  Google Scholar 

  72. 72.

    Palomino C, Fernández-Romero MD, Rubio J, Torres A, Moreno MT, Millan T: Integration of new CAPS and dCAPS-RGA markers into a composite chickpea genetic map and their association with disease resistance. Theor Appl Genet. 2009, 118: 671-682. 10.1007/s00122-008-0928-7.

    Article  CAS  PubMed  Google Scholar 

  73. 73.

    Timmerman-Vaughan GM, Frew TJ, Weerden NF: Characterization and linkage mapping of R-gene analogous DNA sequences in pea (Pisum sativum L). Theor Appl Genet. 2000, 101: 241-247. 10.1007/s001220051475.

    Article  CAS  Google Scholar 

  74. 74.

    Prioul-Gervais S, Deniot G, Receveur EM, Frankewitz A, Fourmann M, Rameau C, Pilet-Nayel ML, Baranger A: Candidate genes for quantitative resistance to Mycosphaerella pinodes in pea (Pisum sativum L.). Theor App Gen. 2007, 114: 971-984. 10.1007/s00122-006-0492-y.

    Article  CAS  Google Scholar 

  75. 75.

    Michelmore RW, Paran I, Kesseli RV: Identification of markers linked to disease-resistance genes by bulked segregate analysis: A rapid method to detect markers in specific genomic regions by using segregating populations. Proc Natl Acad Sci U S A. 1991, 88: 9828-9832. 10.1073/pnas.88.21.9828.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  76. 76.

    Lander ES, Green P, Abramson J, Barlow A, Dali MJ, Lincoln DE, Newburg L: MAPMAKER: an interactive computer program for constructing genetic linkage maps of experimental and natural populations. Genomics. 1987, 1: 174-181. 10.1016/0888-7543(87)90010-3.

    Article  CAS  PubMed  Google Scholar 

  77. 77.

    Wu Y, Bhat PR, Close TJ, Lonardi S: Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 2008, 4: e1000212-10.1371/journal.pgen.1000212.

    PubMed Central  Article  PubMed  Google Scholar 

  78. 78.

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

    Article  Google Scholar 

  79. 79.

    Wang S, Basten CJ, Zeng ZB: Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC. 2010,,

    Google Scholar 

  80. 80.

    Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.

    PubMed Central  CAS  PubMed  Google Scholar 

  81. 81.

    Yap IV, Schneider D, Kleinberg J, Matthews D, Cartinhour S, McCouch SR: A graph-theoretic approach to comparing and integrating genetic, physical and sequence-based maps. Genetics. 2003, 165: 2235-2247.

    PubMed Central  CAS  PubMed  Google Scholar 

  82. 82.

    Muñoz-Amatriaín M, Moscou MJ, Bhat PR, Svensson JT, Bartoš J, Suchánková P, Šimková P, Endo TR, Fenton RD, Lonardi S, Castillo AM, Chao S, Cistué L, Cuesta-Marcos A, Forrest KL, Hayden MJ, Hayes PM, Horsley RD, Makoto K, Moody D, Sato D, Vallés MP, Wulff BBH, Muehlbauer GJ, Doležel J, Close TJ: An improved consensus linkage map of barley based on flow sorted chromosomes and SNP markers. The Plant Genome. 2011, 4: 238-249. 10.3835/plantgenome2011.08.0023.

    Article  Google Scholar 

  83. 83.

    Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 7-78.

    Article  Google Scholar 

Download references


This research has been funded by the Spanish projects AGL2008-02305/AGR, RTA2010-00059 and IPT-2011-1259-010000, co-financed by FEDER, and by the EU project GLIP (FP6- FOOD-CT-2004-506223). The authors thank Juan Prieto, Gloria de Lara and Ana Pozo for excellent technical assistance in the laboratory, greenhouse and field. We are grateful to A. Di Pietro for carefully reading the manuscript.

Author information



Corresponding author

Correspondence to Ana M Torres.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ZS performed the statistical analyses for map merging, integrated QTL information, prepared tables and graphic representations and helped to draft the manuscript. CMA generated marker and QTL data for the mapping populations, participated in the design of the study and assisted with manuscript preparation and editing. SCI, RDR performed additional marker genotyping and provided new QTL data. GGR, CP, NG, SV, SOM and MVG provided new marker data for mapping saturation and anchoring points for map integration. JIC assisted with manuscript editing. AMT coordinated the map integration study, contribute to the analysis and interpretation of data, provided the marker and QTL data and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1: Putative QTLs for flowering time and yield related traits detected in the faba bean RIL population Vf6 × Vf27 (from Cruz-Izquierdo et al., 2012 with modifications). (B) Putative QTLs for Ascochyta fabae, Orobanche crenata and Orobanche foetida resistance detected in the faba bean RIL population Vf6 × Vf136 (Díaz-Ruíz et al., 2009a; 2009b; 2010 and this study). (C) Putative QTLs for Orobanche crenata and Orobanche foetida resistance detected in the faba bean RIL population 29H × Vf136 (from Gutierrez et al., 2013 with modifications). (DOC 103 KB)

Additional file 2: Table S2: Information on the markers used in this study and mapped in the three faba bean RIL populations. (XLS 201 KB)

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

Authors’ original file for figure 3

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 credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Cite this article

Satovic, Z., Avila, C.M., Cruz-Izquierdo, S. et al. A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia fabaL.). BMC Genomics 14, 932 (2013).

Download citation


  • Faba bean
  • Consensus map
  • Gene based markers
  • Quantitative trait loci
  • Comparative mapping
  • Molecular breeding
  • Marker-assisted selection
  • Genomics