<title>Background</title> <p>Faba bean (<italic>Vicia faba </italic>L.) is currently the third most important cool-season food legume in the world. Faba bean provides an important source of dietary protein in human diet, edible oil and animal feeds. Like other grain legumes, faba bean contributes to sustainable agriculture in the management of soil fertility and plays an essential role in crop rotation. Faba bean is a diploid with 2n = 2x = 12 chromosomes [##UREF##0##1##,##UREF##1##2##], is partially cross-pollinated ranging from 4 to 84% [##UREF##2##3##], and possesses one of the largest genomes among crop legumes (~13000 Mb). The development of saturated linkage maps as well as the identification and map-based isolation of important qualitative traits or quantitative trait loci is therefore complex and expensive.</p> <p>A saturated genetic linkage map provides an invaluable tool in plant genetic studies and practical breeding. One of the first genetic linkage maps of faba bean was constructed by Van de Ven <italic>et al</italic>. [##UREF##3##4##] with only 17 markers; followed by successively more detailed genetic maps by Torres <italic>et al</italic>. [##UREF##4##5##] with 51 markers, Satovic <italic>et al</italic>. [##UREF##5##6##] with 157 markers, Vaz Patto <italic>et al</italic>. [##UREF##6##7##] with 116 markers, Roman <italic>et al</italic>. [##REF##12502250##8##] with 121 markers, Roman <italic>et al</italic>. [##REF##15067394##9##] with 192 markers and Avila <italic>et al</italic>. [##REF##15067393##10##] with 103 markers. To date, faba bean genetic maps have been restricted to morphological, isozyme, RFLP, RAPD, a few seed protein genes and four SSR markers. These markers have been limited either in number, transferability or in their ability to provide syntenic information with other legume species.</p> <p>The recent generation of abundant genomic and genetic resources focussed around the model species <italic>Medicago truncatula </italic>and <italic>Lotus Japonicus </italic>has opened up abundant opportunities for creating gene-based molecular markers that are ideal for genetic mapping in general and comparative mapping in particular [##REF##15082563##11##, ####REF##15489274##12##, ##REF##16791689##13##, ##REF##17526914##14##, ##REF##17119911##15##, ##REF##15824281##16####15824281##16##]. In this approach, oligonucleotide primers were designed from sequences of conserved regions in gene exons that flank polymorphic regions such as introns or microsatellites. This PCR-based, codominant marker system has remarkably increased the efficiency of transferring genetic information across species. Examples include the comparison of <italic>M. truncatula </italic>with alfalfa, pea, chickpea, soybean, mung bean, lentil and lupins [##REF##15489274##12##, ####REF##16791689##13##, ##REF##17526914##14##, ##REF##17119911##15##, ##REF##15824281##16##, ##REF##16416153##17##, ##REF##16642337##18####16642337##18##].</p> <p>Comparative genomic studies can expose and confirm phylogenetic relationships among species and determine patterns of chromosomal evolution and syntenic relationships. More importantly, comprehensive comparative genomics can facilitate back-and-forth use of genomic resources between different legumes species, and help to reduce cost and increase efficiency in genetic research as well as crop breeding. The use of conserved genome structure to assist in transferring knowledge among related plant species is well established in grasses [##REF##10899971##19##,##REF##10810140##20##] where synteny greatly assists in gene identification among related species.</p> <p>In this paper we report: (1) the application of gene-based markers in faba bean; (2) the development of the first exclusively gene-based genetic and comparative map for the species; (3) analysis of syntenic relationship between faba bean and <italic>M. truncatula</italic>; and (4) the levels of homology existing between faba bean, <italic>M. truncatula </italic>and lentil, a closely related species to faba bean.</p>
<title>Methods</title> <title>Genetic mapping population</title> <p>A gene-based genetic map of faba bean genome was constructed using a population of 94 F<sub>6 </sub>RILs generated from a cross between faba bean line Vf6 (<italic>equina </italic>type) as the pollen recipient and line Vf27 as pollen donor (<italic>paucijuga </italic>type). These accessions have been widely used in previous genetic studies [##UREF##6##7##,##REF##15067394##9##,##UREF##7##21##] and the population was developed at IFAPA, Centro Alameda del Obispo, in Córdoba, Spain, using diploid parental individuals. Total genomic DNA was isolated from each parent and F<sub>6 </sub>individual as previously described [##REF##16402186##22##].</p> <title>Primer design</title> <p>A total of 796 intron-targeted amplified polymorphic markers (ITAPs) were used for this study. These were composed of four sets of ITAPs; 340 ML and Lup primers developed from alignment of <italic>M. truncatula </italic>and <italic>Lupinus </italic>spp. database EST sequences; 160 MLG primers based on alignments between <italic>M. truncatula</italic>, <italic>Lupinus albus</italic>, and <italic>Glycine max </italic>([##REF##17119911##15##], with a subset of these denoted AtMtL-); 143 cross-species makers (MP) developed by the Department of Plant Pathology, University of California, Davis, USA [##REF##15489274##12##]; and 140 GLIP markers created by the European Grain Legumes Integrated Project (GLIP) based on primarily on <italic>M. truncatula </italic>and pea (Andrea Seres, pers. comm.). The majority of the ITAPs markers could be positioned to a physical location in the <italic>M. truncatula </italic>psuedogenome since most primers were designed from genes in characterized chromosomal regions.</p> <title>Polymorphism detection</title> <p>Each primer pair was screened on <italic>V. faba </italic>parental DNA. PCR conditions were optimised to produce clear single amplicons, and single PCR products of the same size were purified and directly sequenced. Different detection methods were used to genotype the F<sub>6 </sub>population dependent on the type of polymorphism [##REF##17119911##15##]. Details of each marker are given in Additional file ##SUPPL##0##1##.</p> <title>Map construction</title> <p>Chi-squared analysis (P &lt; 0.05) was applied to test the segregation of the mapped markers against the expected Mendelian segregation ratio for co-dominant inheritance in the faba bean F<sub>6 </sub>RIL population. Genetic linkage mapping was conducted with MultiPoint v 1.2 software [##UREF##8##23##], with a recombination fraction (rf) of 0.29 (LOD = 9.0) using 5000 bootstraps. Map distances were calculated in cM by applying the "Kosambi" function. Groups of linked markers that were similarly distorted were accepted for linkage mapping. Independent markers showing significant segregation distortion and markers with missing data (&gt; 10%) were rejected for linkage to avoid bias and false linkages. Genetic maps were drawn with the software program MapChart v 2.1 [##UREF##9##24##].</p> <title>Establishment of macrosyntenic relationships between faba bean, lentil and <italic>M. truncatula</italic></title> <p>Markers mapped in faba bean were located on the <italic>M. truncatula </italic>map by aligning the ESTs originally used to design the ITAPs primers with <italic>M. truncatula </italic>BACs in the Medicago pseudogenome Mt2.0 build 8/10/2007. Alignments with a BLAST E values &lt; 1e<sup>-20</sup>, hsp identity ≥ 60%, and hsp length &gt; 50 nt were retained. Precise positions of markers were obtained by aligning the ITAP primer sequences with <italic>M. truncatula </italic>BACs using BLASTN with an expected value &lt; 1e<sup>-4 </sup>(primer length varied from 18bp to 28bp). Approximate positions of markers that had been genetically mapped in <italic>M. truncatula </italic>but not yet positioned on the physical map were obtained from the <italic>M. truncatula </italic>genetic map [##UREF##10##25##].</p> <title>Loci dot plot created via Grid Map</title> <p>Grid Map [##UREF##11##26##] was used to compare the genetic maps of faba bean and <italic>M. truncatula</italic>. Ordered loci from faba bean and <italic>M. truncatula </italic>linkage groups were listed vertically and horizontally, respectively, and dots were positioned on the diagram at the intersection of the locations of the corresponding markers in the two genetic maps.</p>
<title>Results</title> <title>Gene-based marker development</title> <p>Of the 796 markers screened for amplification in genomic faba bean DNA, 19% produced two or more amplicons of different sizes. Five hundred and seven were selected and optimised for single-locus amplification. Fifty percent (254) of these produced clear single band amplicons (Table ##TAB##0##1##). Seven markers that produced two amplicons (AIGP, GLIP172, GLIP429, GLIP621, GLIP651, Lup226 and MMK1) and one that produced three amplicons (LG34) were also mapped as they exhibited convenient length polymorphisms. One hundred and sixty-five polymorphic markers were identified (Table ##TAB##0##1##), of which 151 markers were used to genotype the 94 individuals of the F<sub>6 </sub>RIL population (supplementary Table ##TAB##0##1##). Apart from the GLIP markers which were selected for their ability to amplify faba bean genomic DNA before this study, 'MLG' and 'MP' markers worked equally well in faba bean (83% and 84%, respectively) but 'MP' markers produced fewer polymorphic amplicons. Amplification rate of the 'ML' markers was the poorest (34%), less than half the rate of the 'MLG' and 'MP' groups (Table ##TAB##0##1##). Likewise, the polymorphism level of the 'ML' markers was the lowest. Both 'MLG' and 'GLIP' markers produced very high levels of polymorphism in the mapping parents (70%; Table ##TAB##0##1##). DNA sequences of markers have been submitted to Genbank under accession codes <ext-link ext-link-type="gen" xlink:href="FH893713">FH893713</ext-link> – <ext-link ext-link-type="gen" xlink:href="FH937528">FH937528</ext-link>.</p> <p>Nineteen mapped markers (&gt; 12%) deviated significantly (P &lt; 0.05) from the expected Mendelian inheritance ratio of 1:1. About half of these markers (9) were highly distorted (P &lt; 0.01). Fifty percent of the distorted markers segregated in favour of the Vf6 parent and fifty percent in favour of Vf27. Nine of the 12 faba bean LGs contained one to four of these distorted markers, which tended to be scattered throughout the faba bean genome. However, it is noteworthy that distorted markers grouped together in chromosomes FB-1, 7 and 8 (Figure ##FIG##0##1##).</p> <title>The first gene-based and comparative map of Faba bean</title> <p>A total of 151 genic markers were used to generate the first gene-based and comparative map of faba bean. This map was constructed with a recombination fraction of 0.29 [LOD = 9.0, 27]. The map consists of seven main linkage groups (FB-1 to FB-7) and 5 fragments (FB-8 to FB-12), which varied in length from 23.6 to 324.8 cM, and spans a total of nearly 1686 cM (Figure ##FIG##0##1## and Table ##TAB##1##2##). The number of markers per LG ranged from three to 30 markers. In addition, there was one pair and thirteen unlinked markers. Eight markers co-segregated at four loci, (three loci in FB-1, one in FB-3; Figure ##FIG##0##1## and Table ##TAB##1##2##). The maximum distances between markers ranged from 13.8 cM in FB-12 to 40.2 cM in FB-3 with an overall mean gap distance of 14.6 cM (Table ##TAB##1##2##).</p> <p>Of the 135 genic markers that mapped to the 12 faba bean LGs in Figure ##FIG##0##1##, 127 were assigned to the <italic>M. truncatula </italic>genetic or physical maps [##UREF##10##25##]. One hundred and four of these (82%) were in syntenic regions. Clear evidence of a simple and direct macrosyntenic relationship between the <italic>V. faba </italic>and <italic>M. truncatula </italic>is presented in the dot matrix in Figure ##FIG##1##2##. The formation of clear isoclinic diagonal lines along the linkage groups provides a strong indication of the conservation of gene order in the two legume genomes. The extensive colinearity was particularly prevailing between FB LG 1, 2, 3, 4 and Mt LG 8, 1, 4, 5 where syntenic regions accounted for 90%, 86%, 50% and 47% of the <italic>M. truncatula </italic>pseudogenome, respectively (Table ##TAB##2##3##).</p> <p>However, chromosomal rearrangements were also evident at a moderate level. For example, <italic>M. truncatula </italic>chromosomes 1 and 7 together with 5 and 2 merged to form the faba bean LGs 2 and 4, respectively. Similarly, <italic>M. truncatula </italic>chromosome 2 splits into FB- 4 and 8 and <italic>M. truncatula </italic>chromosome 3 into FB-5 and 9 (Table ##TAB##2##3## and Figure ##FIG##1##2##). Inversions and translocations were also notable among the orthologous markers within each syntenic pair of faba bean and <italic>M. truncatula </italic>LGs (Figure ##FIG##1##2##).</p> <title>Evidence of macrosynteny among faba bean, lentil and Medicago truncatula</title> <p>A high level of co-linearity was found between the faba bean, lentil and <italic>M. truncatula </italic>genomes based on the macro-synteny established between faba bean and <italic>M. truncatula </italic>(this study) or lentil and <italic>M. truncatula </italic>[##REF##17119911##15##] using only common orthologous markers which mapped in all three species (Figures ##FIG##2##3## and ##FIG##3##4##). The pattern of homology between faba bean and <italic>M. truncatula </italic>was similar to that between lentil and <italic>M. truncatula</italic>: for example, two linkage groups FB-1 and Len-II were exclusively syntenic to <italic>M. truncatula </italic>LG_8 and shared seven markers in common. Common markers were evenly distributed in all the three corresponding LGs suggesting that FB-1 and Len-II are essentially co-linear (Figure ##FIG##2##3A##). Likewise, FB-2 and Len-III were both syntenic to Mt-1 and orthologous to each other with nine markers in common (Figure ##FIG##2##3B##). Other examples are FB-3, Len-I and Mt-4; and FB-4, Len-V and Mt-5 (Figure ##FIG##1##2##, ##FIG##3##4A## and Phan et al., 2007 [##REF##17119911##15##]). FB-5 and FB-9 were co-linear with lentil Len-VII and LenVI, respectively, and both pairs of these linkage groups were colinear with Mt-3 (Figure ##FIG##3##4B##). This suggests shared ancestral chromosomal changes in faba bean and lentil compared to <italic>M. truncatula </italic>and confirms their phylogenetically closer relationship.</p> <title>Evidence of fine scale macrosynteny on individual BACs between faba bean, lentil and Medicago truncatula</title> <p>Several ITAP markers used in this study were designed from different genes on the same <italic>M. truncatula </italic>BACs. These markers often mapped in clusters in two or more species and in common syntenic regions of corresponding linkage groups. An example is shown in Figure ##FIG##2##3A## with markers LG89-LG98 designed from <italic>M. truncatula </italic>BAC AC140032 and LG83-LG88 from <italic>M. truncatula </italic>BAC AC138131 in <italic>M. truncatula </italic>chromosome 8. Similarly, with the exception of LG34c, markers LG31-LG43 designed from <italic>M. truncatula </italic>BAC AC152751 in Mt chromosome 1 were located in one area of FB-2 and Len-III as shown in Figure ##FIG##2##3B##; markers AnMtL6-AnMtL8 and LG102-LG112 from <italic>M. truncatula </italic>BAC AC147712 and AC135800 in Mt chromosome 3 were grouped together on corresponding pairs of LGs i.e. FB-5 and Len-VII and FB-9 and Len-VI; respectively (Figure ##FIG##3##4B##).</p>
<title>Discussion and conclusion</title> <p>The first genetic map of faba bean composed exclusively with gene-based co-dominant molecular markers was constructed using a F<sub>6 </sub>RIL population between lines Vf6 and Vf27. The map is also the first to enable the establishment of syntenic relationships between faba bean and the model legume <italic>M. truncatula</italic>, comparison with other legume species, and integration with genetic maps available in faba bean.</p> <p>The map is composed of 12 linkage groups and 151 genetic markers. Although the number of chromosomes in faba bean has been reported as 2n = 12 [##UREF##1##2##], the number of linkage groups in recent genetic maps in the species range from 13 to 18 [##UREF##6##7##, ####REF##12502250##8##, ##REF##15067394##9##, ##REF##15067393##10####15067393##10##] and previously as many as 48 have been reported [##UREF##5##6##]. The high number of linkage groups compared to the number of chromosomes may be due to the fact that faba bean possesses one of the largest genomes among cultivated legumes (~13000 Mb). This compares with other well-characterised species such as <italic>M. truncatula</italic>, chickpea, soybean, lentil and pea which have genomes of ~450 Mb, ~740 Mb, ~1200 Mb, ~4000 Mb and ~4000 Mb respectively [##UREF##13##28##].</p> <p>Of the 24 non-orthologous markers found in this study, eight were from primer pairs where more than one PCR gel band was present and where two or more such amplicons were mapped. In each case, at least one amplicon mapped syntenically. The percentages of markers sequenced in faba bean were lower compared to lentil (63%, 26% and 55% compared to 93%, 69% and 65% for MP, ML and MLG markers respectively, Table ##TAB##0##1##). This was due to higher proportion of markers amplifying multiple bands in faba bean compared to lentil (data not shown), which may imply duplication. Differences in amplification, sequencing and polymorphism rates among different types of markers used for this study reflect the mode of design of the markers. Since 'MP' and 'MLG' markers were often based on the homology of more than two phylogenetically distant species, they are more likely to work in different legume lineages. The same observation was reported for these primer sets in lentil [##REF##17119911##15##].</p> <p>Despite the large differences in genome sizes between <italic>M. truncatula </italic>and <italic>V. faba</italic>, a simple and direct relationship between the two genomes was identified in this study. Given the number of markers used (151), the syntenic regions cover a large proportion of <italic>M. truncatula </italic>pseudogenome with 90%, 87%, 66%, 62% and 47% for <italic>M. truncatula </italic>chromosomes 8, 1, 3, 4 and 5, respectively (Table ##TAB##2##3##). The appearance of clear isoclinic diagonal lines along the linkage groups in Figure ##FIG##1##2## also demonstrates strong evidence for the extensive co-linearity between linkage group pairs of the two species. Similar high levels of conservation have also been reported between <italic>L. culinaris </italic>ssp. <italic>culinaris </italic>and <italic>M. truncatula </italic>[##REF##17119911##15##] and other closely related legumes such as <italic>L. culinaris </italic>ssp. <italic>culinaris </italic>and <italic>P. sativum </italic>[##UREF##14##29##], <italic>M. sativa </italic>and <italic>P. sativum </italic>[##REF##15340836##30##], <italic>M. truncatula </italic>and <italic>P. sativum </italic>[##REF##16416153##17##], <italic>M. truncatula </italic>and <italic>M. sativa </italic>[##REF##15082563##11##]. This study also shows markers originally designed from genes on the same BAC clustered in corresponding syntenic areas in lentil and faba bean. The mapping populations were too small to resolve marker order in lentil and faba bean but extensive conservation of gene order,(and microsynteny) has been shown in previous studies between other legume species at similar or greater phylogenetic distances [##REF##15489274##12##,##REF##12175015##31##, ####REF##14557858##32##, ##REF##16102170##33##, ##REF##17003129##34####17003129##34##], and to some extent between <italic>M. truncatula </italic>and <italic>Arabidopsis </italic>[##REF##16102170##33##, ####REF##17003129##34##, ##REF##12644654##35####12644654##35##].</p> <p>A higher level of homology between <italic>V. faba </italic>and <italic>L. culinaris </italic>ssp. <italic>culinaris </italic>compared to that between <italic>V. faba </italic>and <italic>M. truncatula </italic>could be inferred from this study based on the common markers mapped in the two genomes, common homology with <italic>M. truncatula </italic>and similar pattern of rearrangements (Figures ##FIG##2##3##, ##FIG##3##4## and Phan <italic>et al</italic>., 2007 [##REF##17119911##15##]). This finding agrees with phylogenetic studies that place the genera <italic>Vicia</italic>, <italic>Lens </italic>and <italic>Pisum </italic>within the tribe Viceae while <italic>Medicago </italic>and <italic>Melilotus </italic>form a parallel tribe Trifolieae within the Galegoid or cool season legumes [##REF##12644643##36##], and is consistent with different levels of macrosynteny observed between <italic>M. truncatula</italic>, <italic>P. sativum</italic>, <italic>V. radiata</italic>, <italic>G. max</italic>, and <italic>Phaceolus vulgaris </italic>dependent on phylogenetic distance [##REF##15489274##12##]. However, chromosomal rearrangements were evident (Figures ##FIG##1##2## and ##FIG##2##3##).</p> <p>Rearrangements involving Mt6 and Mt3 in particular may explain the differences in chromosome number between the two species (<italic>M. truncatula</italic>: n = 8; <italic>V. faba</italic>: n = 6). Mt6 might be considered unusual and is largely composed of heterochromatic DNA [##REF##11489182##37##], contains few transcribed genes [##REF##15082563##11##] and a large proportion of resistance gene analogues [##REF##12059101##38##]. In this study no corresponding linkage group was detected in faba bean, as found previously in pea [##REF##15489274##12##] and lupin [##REF##16791689##13##,##REF##17526914##14##], together with less than five percent estimated coverage by the <italic>L. japonicus </italic>genome [##REF##17003129##34##]. In faba bean FB5 and FB9 appear to correspond to Mt3. This configuration is supported by a similar pattern in lentil (Figure ##FIG##3##4B## and Phan <italic>et al</italic>., 2007 [##REF##17119911##15##]), although a larger number of markers are needed to confirm this.</p> <p>The faba bean comparative map constructed here is consistent with the pattern of chromosome conservation previously observed, where different levels of conservation were found to be relatively consistent between <italic>M. truncatula </italic>and other legume species i.e. high conservation of <italic>M. truncatula </italic>chromosomes 1, 5 and 8; moderate conservation in the <italic>M. truncatula </italic>chromosomes 2, 3, 4 and 7 and lowest conservation in the <italic>M. truncatula </italic>chromosome 6 (Figures ##FIG##2##3##, ##FIG##3##4## and [##UREF##15##39##]). As described above, no homology was identified with <italic>M. truncatula </italic>chromosome 6 in this study. The alignment of this faba bean map with lentil and the current <italic>M. truncatula </italic>genome based on <italic>M. truncatula </italic>genome assembly Mt2.0 is slightly different to that based on an earlier assembly [##REF##17119911##15##]. The changes can be observed in Figure ##FIG##2##3B## where orthologous markers which were syntenic to <italic>M. truncatula </italic>chromosome 6 in lentil are now co-linear with <italic>M. truncatula </italic>chromosome 1 in common with faba bean.</p> <p>Genome studies have demonstrated different factors are responsible for genome size variation and speciation. These include ancient polyploidisation events in the case of the Brassicas [##REF##8978073##40##]; segmental or region-specific duplication [##REF##14615182##41##]; and genetic rearrangements, transposable element amplification, or combination of different genome modifications [##REF##17291821##42##]. Large scale rearrangements, duplications, or polyploidisation were not apparent in this study, possibly as a result of the focus on single locus markers, however differences in non-coding repetitive DNA or transposable elements provide a possible explanation for the large differences in genome size. Retroelements are known to account for substantial proportions of these Viceae genomes as shown by extensive studies in pea, for example [##REF##18031571##43##,##UREF##16##44##], and more recently <italic>Vicia </italic>[##REF##16585134##45##,##REF##8602145##46##]. Local genic rearrangements similar to that in found in the grasses (duplications, translocations, and insertions or deletions) may explain multiple PCR amplicons [##REF##10899971##19##,##REF##10810140##20##].</p> <p>The shared macrosynteny among the three species demonstrated here and even higher level of homology between <italic>L. culinaris </italic>and <italic>V. faba </italic>will undoubtedly facilitate the identification of markers closely linked to traits of interest in <italic>V. faba</italic>. Alignment of this map with existing faba bean maps containing important traits with polymorphic SSR markers and/or markers developed in this study, coupled with cross-reference to the abundant genetic information from the Medicago genome sequencing and extensive EST libraries available for the model legume species, will undoubtedly assist this process. As the parental line Vf6 has been used in a number of genetic and QTL mapping projects [##UREF##5##6##,##UREF##6##7##,##REF##15067394##9##], this map can serve as a central reference map. This study has provided a number of significant outcomes for faba bean genomics and legume genomics in general.</p>
<title>Discussion and conclusion</title> <p>The first genetic map of faba bean composed exclusively with gene-based co-dominant molecular markers was constructed using a F<sub>6 </sub>RIL population between lines Vf6 and Vf27. The map is also the first to enable the establishment of syntenic relationships between faba bean and the model legume <italic>M. truncatula</italic>, comparison with other legume species, and integration with genetic maps available in faba bean.</p> <p>The map is composed of 12 linkage groups and 151 genetic markers. Although the number of chromosomes in faba bean has been reported as 2n = 12 [##UREF##1##2##], the number of linkage groups in recent genetic maps in the species range from 13 to 18 [##UREF##6##7##, ####REF##12502250##8##, ##REF##15067394##9##, ##REF##15067393##10####15067393##10##] and previously as many as 48 have been reported [##UREF##5##6##]. The high number of linkage groups compared to the number of chromosomes may be due to the fact that faba bean possesses one of the largest genomes among cultivated legumes (~13000 Mb). This compares with other well-characterised species such as <italic>M. truncatula</italic>, chickpea, soybean, lentil and pea which have genomes of ~450 Mb, ~740 Mb, ~1200 Mb, ~4000 Mb and ~4000 Mb respectively [##UREF##13##28##].</p> <p>Of the 24 non-orthologous markers found in this study, eight were from primer pairs where more than one PCR gel band was present and where two or more such amplicons were mapped. In each case, at least one amplicon mapped syntenically. The percentages of markers sequenced in faba bean were lower compared to lentil (63%, 26% and 55% compared to 93%, 69% and 65% for MP, ML and MLG markers respectively, Table ##TAB##0##1##). This was due to higher proportion of markers amplifying multiple bands in faba bean compared to lentil (data not shown), which may imply duplication. Differences in amplification, sequencing and polymorphism rates among different types of markers used for this study reflect the mode of design of the markers. Since 'MP' and 'MLG' markers were often based on the homology of more than two phylogenetically distant species, they are more likely to work in different legume lineages. The same observation was reported for these primer sets in lentil [##REF##17119911##15##].</p> <p>Despite the large differences in genome sizes between <italic>M. truncatula </italic>and <italic>V. faba</italic>, a simple and direct relationship between the two genomes was identified in this study. Given the number of markers used (151), the syntenic regions cover a large proportion of <italic>M. truncatula </italic>pseudogenome with 90%, 87%, 66%, 62% and 47% for <italic>M. truncatula </italic>chromosomes 8, 1, 3, 4 and 5, respectively (Table ##TAB##2##3##). The appearance of clear isoclinic diagonal lines along the linkage groups in Figure ##FIG##1##2## also demonstrates strong evidence for the extensive co-linearity between linkage group pairs of the two species. Similar high levels of conservation have also been reported between <italic>L. culinaris </italic>ssp. <italic>culinaris </italic>and <italic>M. truncatula </italic>[##REF##17119911##15##] and other closely related legumes such as <italic>L. culinaris </italic>ssp. <italic>culinaris </italic>and <italic>P. sativum </italic>[##UREF##14##29##], <italic>M. sativa </italic>and <italic>P. sativum </italic>[##REF##15340836##30##], <italic>M. truncatula </italic>and <italic>P. sativum </italic>[##REF##16416153##17##], <italic>M. truncatula </italic>and <italic>M. sativa </italic>[##REF##15082563##11##]. This study also shows markers originally designed from genes on the same BAC clustered in corresponding syntenic areas in lentil and faba bean. The mapping populations were too small to resolve marker order in lentil and faba bean but extensive conservation of gene order,(and microsynteny) has been shown in previous studies between other legume species at similar or greater phylogenetic distances [##REF##15489274##12##,##REF##12175015##31##, ####REF##14557858##32##, ##REF##16102170##33##, ##REF##17003129##34####17003129##34##], and to some extent between <italic>M. truncatula </italic>and <italic>Arabidopsis </italic>[##REF##16102170##33##, ####REF##17003129##34##, ##REF##12644654##35####12644654##35##].</p> <p>A higher level of homology between <italic>V. faba </italic>and <italic>L. culinaris </italic>ssp. <italic>culinaris </italic>compared to that between <italic>V. faba </italic>and <italic>M. truncatula </italic>could be inferred from this study based on the common markers mapped in the two genomes, common homology with <italic>M. truncatula </italic>and similar pattern of rearrangements (Figures ##FIG##2##3##, ##FIG##3##4## and Phan <italic>et al</italic>., 2007 [##REF##17119911##15##]). This finding agrees with phylogenetic studies that place the genera <italic>Vicia</italic>, <italic>Lens </italic>and <italic>Pisum </italic>within the tribe Viceae while <italic>Medicago </italic>and <italic>Melilotus </italic>form a parallel tribe Trifolieae within the Galegoid or cool season legumes [##REF##12644643##36##], and is consistent with different levels of macrosynteny observed between <italic>M. truncatula</italic>, <italic>P. sativum</italic>, <italic>V. radiata</italic>, <italic>G. max</italic>, and <italic>Phaceolus vulgaris </italic>dependent on phylogenetic distance [##REF##15489274##12##]. However, chromosomal rearrangements were evident (Figures ##FIG##1##2## and ##FIG##2##3##).</p> <p>Rearrangements involving Mt6 and Mt3 in particular may explain the differences in chromosome number between the two species (<italic>M. truncatula</italic>: n = 8; <italic>V. faba</italic>: n = 6). Mt6 might be considered unusual and is largely composed of heterochromatic DNA [##REF##11489182##37##], contains few transcribed genes [##REF##15082563##11##] and a large proportion of resistance gene analogues [##REF##12059101##38##]. In this study no corresponding linkage group was detected in faba bean, as found previously in pea [##REF##15489274##12##] and lupin [##REF##16791689##13##,##REF##17526914##14##], together with less than five percent estimated coverage by the <italic>L. japonicus </italic>genome [##REF##17003129##34##]. In faba bean FB5 and FB9 appear to correspond to Mt3. This configuration is supported by a similar pattern in lentil (Figure ##FIG##3##4B## and Phan <italic>et al</italic>., 2007 [##REF##17119911##15##]), although a larger number of markers are needed to confirm this.</p> <p>The faba bean comparative map constructed here is consistent with the pattern of chromosome conservation previously observed, where different levels of conservation were found to be relatively consistent between <italic>M. truncatula </italic>and other legume species i.e. high conservation of <italic>M. truncatula </italic>chromosomes 1, 5 and 8; moderate conservation in the <italic>M. truncatula </italic>chromosomes 2, 3, 4 and 7 and lowest conservation in the <italic>M. truncatula </italic>chromosome 6 (Figures ##FIG##2##3##, ##FIG##3##4## and [##UREF##15##39##]). As described above, no homology was identified with <italic>M. truncatula </italic>chromosome 6 in this study. The alignment of this faba bean map with lentil and the current <italic>M. truncatula </italic>genome based on <italic>M. truncatula </italic>genome assembly Mt2.0 is slightly different to that based on an earlier assembly [##REF##17119911##15##]. The changes can be observed in Figure ##FIG##2##3B## where orthologous markers which were syntenic to <italic>M. truncatula </italic>chromosome 6 in lentil are now co-linear with <italic>M. truncatula </italic>chromosome 1 in common with faba bean.</p> <p>Genome studies have demonstrated different factors are responsible for genome size variation and speciation. These include ancient polyploidisation events in the case of the Brassicas [##REF##8978073##40##]; segmental or region-specific duplication [##REF##14615182##41##]; and genetic rearrangements, transposable element amplification, or combination of different genome modifications [##REF##17291821##42##]. Large scale rearrangements, duplications, or polyploidisation were not apparent in this study, possibly as a result of the focus on single locus markers, however differences in non-coding repetitive DNA or transposable elements provide a possible explanation for the large differences in genome size. Retroelements are known to account for substantial proportions of these Viceae genomes as shown by extensive studies in pea, for example [##REF##18031571##43##,##UREF##16##44##], and more recently <italic>Vicia </italic>[##REF##16585134##45##,##REF##8602145##46##]. Local genic rearrangements similar to that in found in the grasses (duplications, translocations, and insertions or deletions) may explain multiple PCR amplicons [##REF##10899971##19##,##REF##10810140##20##].</p> <p>The shared macrosynteny among the three species demonstrated here and even higher level of homology between <italic>L. culinaris </italic>and <italic>V. faba </italic>will undoubtedly facilitate the identification of markers closely linked to traits of interest in <italic>V. faba</italic>. Alignment of this map with existing faba bean maps containing important traits with polymorphic SSR markers and/or markers developed in this study, coupled with cross-reference to the abundant genetic information from the Medicago genome sequencing and extensive EST libraries available for the model legume species, will undoubtedly assist this process. As the parental line Vf6 has been used in a number of genetic and QTL mapping projects [##UREF##5##6##,##UREF##6##7##,##REF##15067394##9##], this map can serve as a central reference map. This study has provided a number of significant outcomes for faba bean genomics and legume genomics in general.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>The development of genetic markers is complex and costly in species with little pre-existing genomic information. Faba bean possesses one of the largest and least studied genomes among cultivated crop plants and no gene-based genetic maps exist. Gene-based orthologous markers allow chromosomal regions and levels of synteny to be characterised between species, reveal phylogenetic relationships and chromosomal evolution, and enable targeted identification of markers for crop breeding. In this study orthologous codominant cross-species markers have been deployed to produce the first exclusively gene-based genetic linkage map of faba bean (<italic>Vicia faba</italic>), using an F<sub>6 </sub>population developed from a cross between the lines Vf6 (<italic>equina </italic>type) and Vf27 (<italic>paucijuga </italic>type).</p> <title>Results</title> <p>Of 796 intron-targeted amplified polymorphic (ITAP) markers screened, 151 markers could be used to construct a comparative genetic map. Linkage analysis revealed seven major and five small linkage groups (LGs), one pair and 12 unlinked markers. Each LG was comprised of three to 30 markers and varied in length from 23.6 cM to 324.8 cM. The map spanned a total length of 1685.8 cM. A simple and direct macrosyntenic relationship between faba bean and <italic>Medicago truncatula </italic>was evident, while faba bean and lentil shared a common rearrangement relative to <italic>M. truncatula</italic>. One hundred and four of the 127 mapped markers in the 12 LGs, which were previously assigned to <italic>M. truncatula </italic>genetic and physical maps, were found in regions syntenic between the faba bean and <italic>M. truncatula </italic>genomes. However chromosomal rearrangements were observed that could explain the difference in chromosome numbers between these three legume species. These rearrangements suggested high conservation of <italic>M. truncatula </italic>chromosomes 1, 5 and 8; moderate conservation of chromosomes 2, 3, 4 and 7 and no conservation with <italic>M. truncatula </italic>chromosome 6. Multiple PCR amplicons and comparative mapping were suggestive of small-scale duplication events in faba bean. This study also provides a preliminary indication for finer scale macrosynteny between <italic>M. truncatula</italic>, lentil and faba bean. Markers originally designed from genes on the same <italic>M. truncatula </italic>BACs were found to be grouped together in corresponding syntenic areas in lentil and faba bean.</p> <title>Conclusion</title> <p>Despite the large size of the faba bean genome, comparative mapping did not reveal evidence for polyploidisation, segmental duplication, or significant rearrangements compared to <italic>M. truncatula</italic>, although a bias in the use of single locus markers may have limited the detection of duplications. Non-coding repetitive DNA or transposable element content provides a possible explanation for the difference in genome sizes. Similar patterns of rearrangements in faba bean and lentil compared to <italic>M. truncatula </italic>support phylogenetic studies dividing these species into the tribes Viceae and Trifoliae. However, substantial macrosynteny was apparent between faba bean and <italic>M. truncatula</italic>, with the exception of chromosome 6 where no orthologous markers were found, confirming previous investigations suggesting chromosome 6 is atypical. The composite map, anchored with orthologous markers mapped in <italic>M. truncatula</italic>, provides a central reference map for future use of genomic and genetic information in faba bean genetic analysis and breeding.</p>
<title>Authors' contributions</title> <p>SE and RPO designed the research. HTTP and SE wrote the manuscript. HTTP, SE and MG performed marker polymorphism discovery and population genotyping in Perth. JH assisted with informatics procedures. AMT, CMA and SC-I were involved in developing recombinant inbred lines and marker polymorphism discovery of GLIP markers in Córdoba. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This research was supported by an ARC Linkage project LP0454871, the New South Wales Departments of Primary Industry and the GLIP project FP6-2002-FOOD-1-506223.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>A gene-based genetic linkage map of faba bean (<italic>Vicia faba </italic>L)</bold>. Marker distance is given in cM. *indicates markers with distorted segregation.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Matrix plot of common gene-based markers mapped in faba bean and <italic>M. truncatula</italic></bold>. The faba bean and <italic>M. truncatula </italic>loci are listed vertically and horizontally, respectively, according to their linkage group order.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Evidence of shared macrosynteny between <italic>V. faba </italic>chromosomes FB-1 and FB-2, <italic>L. culinaris </italic>and <italic>M. truncatula</italic></bold>. Common orthologous markers are depicted by dashed lines and marker distances are provided in centi-Morgans for <italic>M. truncatula </italic>only. The figures exclude markers that could not be positioned in the <italic>M. truncatula </italic>psuedogenome. A ~ indicates markers that map distally in the corresponding <italic>M. truncatula </italic>chromosome relative to faba bean, and markers suffixed a, b or c denotes derivation from primer pairs that produced multiple PCR products.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Evidence of shared macrosynteny between <italic>V. faba </italic>chromosomes FB-3 and FB-5, <italic>L. culinaris </italic>and <italic>M. truncatula</italic></bold>. Common orthologous markers are depicted by dashed lines and marker distances are provided in centi-Morgans for <italic>M. truncatula </italic>only. The figures exclude markers that could not be positioned in the <italic>M. truncatula </italic>psuedogenome. Markers highlighted * were previously mapped [##REF##16642337##18##] but not positioned in the <italic>M. truncatula </italic>psuedogenome. # indicates markers orthologous between lentil and faba bean but unmapped in <italic>M. truncatula </italic>and markers suffixed a, b or c denotes derivation from primer pairs that produced multiple PCR products.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Efficiency of gene-based markers used to construct the comparative genetic linkage map of faba bean</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Marker type</td><td align="left">Screened</td><td align="left">Amplification<sup>a</sup></td><td align="left">Length polymorphism</td><td align="left">Sequenced<sup>b</sup></td><td align="left">CAPS/SNP<sup>c</sup></td><td align="left">Mapped</td></tr></thead><tbody><tr><td align="left">MP</td><td align="left">143</td><td align="left">120 (84%)</td><td align="center">6</td><td align="left">75 (63%)</td><td align="left">45 (60%)</td><td align="left">36</td></tr><tr><td align="left">ML</td><td align="left">340</td><td align="left">115 (34%)</td><td align="center">9</td><td align="left">30 (26%)</td><td align="left">16 (53%)</td><td align="left">25</td></tr><tr><td align="left">MLG</td><td align="left">160</td><td align="left">132 (83%)</td><td align="center">5</td><td align="left">73 (55%)</td><td align="left">50 (70%)</td><td align="left">42</td></tr><tr><td align="left">GLIP</td><td align="left">140</td><td align="left">140</td><td align="center">6</td><td align="left">76 (54%)</td><td align="left">54 (71%)</td><td align="left">48</td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">Total</td><td align="left">796</td><td align="left">507</td><td align="center">26</td><td align="left">254</td><td align="left">165</td><td align="left">151</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Properties of the faba bean comparative map</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Linkage group</td><td align="left">Length of LGs (cM)</td><td align="left">No of markers</td><td align="left">No of loci</td><td align="left">Average marker spacing<sup>a </sup>(cM)</td><td align="left">Largest distance between markers (cM)</td></tr></thead><tbody><tr><td align="center">LG-1</td><td align="center">324.8</td><td align="center">30</td><td align="center">27</td><td align="center">12.5</td><td align="center">32.5</td></tr><tr><td align="center">LG-2</td><td align="center">313.1</td><td align="center">27</td><td align="center">27</td><td align="center">12.0</td><td align="center">31.8</td></tr><tr><td align="center">LG-3</td><td align="center">238.5</td><td align="center">18</td><td align="center">17</td><td align="center">14.9</td><td align="center">40.2</td></tr><tr><td align="center">LG-4</td><td align="center">210.7</td><td align="center">13</td><td align="center">13</td><td align="center">17.6</td><td align="center">29.8</td></tr><tr><td align="center">LG-5</td><td align="center">151.9</td><td align="center">12</td><td align="center">12</td><td align="center">13.8</td><td align="center">26.8</td></tr><tr><td align="center">LG-6</td><td align="center">120.7</td><td align="center">8</td><td align="center">8</td><td align="center">17.2</td><td align="center">36.9</td></tr><tr><td align="center">LG-7</td><td align="center">119.2</td><td align="center">7</td><td align="center">7</td><td align="center">19.9</td><td align="center">35.9</td></tr><tr><td align="center">LG-8</td><td align="center">49.9</td><td align="center">6</td><td align="center">6</td><td align="center">10.0</td><td align="center">29.5</td></tr><tr><td align="center">LG-9</td><td align="center">64.1</td><td align="center">5</td><td align="center">5</td><td align="center">16.0</td><td align="center">22.6</td></tr><tr><td align="center">LG-10</td><td align="center">41.8</td><td align="center">3</td><td align="center">3</td><td align="center">15.4</td><td align="center">30.7</td></tr><tr><td align="center">LG-11</td><td align="center">27.5</td><td align="center">3</td><td align="center">3</td><td align="center">13.8</td><td align="center">19.9</td></tr><tr><td align="center">LG12</td><td align="center">23.6</td><td align="center">3</td><td align="center">3</td><td align="center">11.8</td><td align="center">13.8</td></tr><tr><td colspan="6"><hr/></td></tr><tr><td align="center">Total</td><td align="center">1685.8</td><td align="center">135</td><td align="center">131</td><td/><td/></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Colinearity between the faba bean and <italic>M. truncatula </italic>genomes</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">FB LG</td><td align="center">No. of markers</td><td align="left">No. of non-orthologous markers</td><td align="left">No. of unassigned markers</td><td align="left">Current Mt genetic coverage (LG: cM, No. colinear markers)<sup>1,2</sup></td><td align="left">Current Mt pseudogenome coverage (bp)<sup>1</sup></td><td align="left">Current faba bean genetic coverage (cM)<sup>2</sup></td></tr></thead><tbody><tr><td align="left">FB-1</td><td align="center">30</td><td align="center">6</td><td align="center">2</td><td align="left">8: 0 – 68.2, 11</td><td align="left">8: 736241 – 32433863 (90%)</td><td align="left">0 – 324.8 (100%)</td></tr><tr><td align="left">FB-2</td><td align="center">27</td><td align="center">4</td><td align="center">2</td><td align="left">1: 0 – 58.5, 9 7: 58.4 – 60, 3</td><td align="left">1: 59849 – 26963044 (86%) 7: 25725918 – 26383168 (2%)</td><td align="left">152.6 – 313.1 (51.3%) 157.1 – 210.7 (25.4%)</td></tr><tr><td align="left">FB-3</td><td align="center">18</td><td align="center">1</td><td align="center">3</td><td align="left">4: 48.1 – 59, 11</td><td align="left">4: 14544831 – 34128326 (50%)</td><td align="left">0 – 238.5 (100%)</td></tr><tr><td align="left">FB-4</td><td align="center">13</td><td align="center">2</td><td align="center">0</td><td align="left">5: 0 – 32, 6 2: 12.7 – 37.9, 3</td><td align="left">5: 594897 – 17960204 (47%) 2: 8491236 – 15201671 (24%)</td><td align="left">0 – 142 (67%) 250.5 – 256.3 (1.9%)</td></tr><tr><td align="left">FB-5</td><td align="center">12</td><td align="center">1</td><td align="center">1</td><td align="left">3: 28 – 70.3, 9</td><td align="left">3: 10888141 – 28213807 (46%)</td><td align="left">0 – 140.3 (92%)</td></tr><tr><td align="left">FB6</td><td align="center">8</td><td align="center">4</td><td align="center">0</td><td align="left">1: 2 – 2.2, NA</td><td align="left">1: 4942324 – 5151182 (0.7%)</td><td align="left">-</td></tr><tr><td align="left">FB-7</td><td align="center">7</td><td align="center">0</td><td align="center">0</td><td align="left">7: 22.6 – 52.8, 4</td><td align="left">7: 9999559 – 22150020 (38%)</td><td align="left">0 – 58.9 (49%)</td></tr><tr><td align="left">FB-8</td><td align="center">6</td><td align="center">1</td><td align="center">0</td><td align="left">2: 48 – 57.8, 5</td><td align="left">2: 17523205 – 20646361 (11%)</td><td align="left">0 – 49.9 (100%)</td></tr><tr><td align="left">FB-9</td><td align="center">5</td><td align="center">0</td><td align="center">0</td><td align="left">3: 62 – 72.5, 5</td><td align="left">3: 19000000* – 30712150 (31%)</td><td align="left">0 – 64.1 (100%)</td></tr><tr><td align="left">FB-10</td><td align="center">3</td><td align="center">1</td><td align="center">0</td><td align="left">4: 0 – 7.4, NA</td><td align="left">4: 537530 – 3489843 (8%)</td><td align="left">-</td></tr><tr><td align="left">FB-11</td><td align="center">3</td><td align="center">3</td><td align="center">0</td><td align="left">NA</td><td align="left">NA</td><td align="left">-</td></tr><tr><td align="left">FB-12</td><td align="center">3</td><td align="center">0</td><td align="center">0</td><td align="left">4: 60.4 – 61.1, NA</td><td align="left">4: 36664007 – 38186081 (4%)</td><td align="left">-</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Supplementary table 1. Orthologous PCR markers developed in this study and genetically mapped in a F<sub>6 </sub>faba bean population between lines Vf6 and Vf27.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>a</sup>Figures in parentheses are percentages of amplified markers of the total markers screened.</p><p><sup>b</sup>Figures in parentheses are percentages of sequenced markers of the total amplified markers.</p><p><sup>c</sup>Figures in parentheses are percentages of polymorphic markers of the total sequenced markers.</p></table-wrap-foot> <table-wrap-foot><p><sup>a</sup>Calculated by dividing the length of the chromosome (cM) by the number of space/distance between markers/loci.</p></table-wrap-foot> <table-wrap-foot><p>NA: not applicable</p><p><sup>1</sup>Data was either based on information from <italic>M. truncatula </italic>genome sequencing website [##UREF##10##25##] or from Choi <italic>et al</italic>. [##REF##15489274##12##].</p><p><sup>2</sup>Data is provided for three or more colinear markers and excludes markers on the same chromosome rearranged relative to the order in <italic>M. truncatula</italic>.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-380-1"/> <graphic xlink:href="1471-2164-9-380-2"/> <graphic xlink:href="1471-2164-9-380-3"/> <graphic xlink:href="1471-2164-9-380-4"/>
<media xlink:href="1471-2164-9-380-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Michaelis", "Rieger"], "given-names": ["A", "R"], "article-title": ["Strukturheterozygotie bei Vicia faba"], "source": ["Zuchter"], "year": ["1959"], "volume": ["29"], "fpage": ["354"], "lpage": ["361"], "pub-id": ["10.1007/BF00709943"]}, {"surname": ["Sjodin"], "given-names": ["J"], "article-title": ["Induced morphological variation in Vicia faba L"], "source": ["Hereditas"], "year": ["1971"], "volume": ["67"], "fpage": ["155"], "lpage": ["180"]}, {"surname": ["Bond", "Poulsen", "Hebblethwaite PD"], "given-names": ["DA", "MH"], "article-title": ["Pollination"], "source": ["The Faba Bean"], "year": ["1983"], "publisher-name": [" Butterworths"], "fpage": ["77"], "lpage": ["101"]}, {"surname": ["Van de Ven", "Waugh", "Duncan", "Ramsay", "Dow", "Powell"], "given-names": ["WTG", "R", "N", "G", "N", "W"], "article-title": ["Development of a genetic linkage map in Vicia faba using molecular and biochemical techniques"], "source": ["Aspects Appl Biol"], "year": ["1991"], "volume": ["27"], "fpage": ["49"], "lpage": ["54"]}, {"surname": ["Torres", "Weeden", "Mart\u00edn"], "given-names": ["AM", "NF", "A"], "article-title": ["Linkage among isozyme, RFLP and RAPD markers in Vicia faba"], "source": ["Theor Appl Genet"], "year": ["1993"], "volume": ["85"], "fpage": ["937"], "lpage": ["945"], "pub-id": ["10.1007/BF00215032"]}, {"surname": ["Satovic", "Torres", "Cubero"], "given-names": ["Z", "AM", "JI"], "article-title": ["Genetic mapping of new morphological, isozyme and RAPD markers in Vicia faba L. using trisomics"], "source": ["Theor Appl Genet"], "year": ["1996"], "volume": ["93"], "fpage": ["1130"], "lpage": ["1138"], "pub-id": ["10.1007/BF00230136"]}, {"surname": ["Vaz Patto", "Torres", "Koblizkova", "Macas", "Cubero"], "given-names": ["MC", "AM", "A", "J", "JI"], "article-title": ["Development of a genetic composite map of Vicia faba using F2 populations derived from trisomic plants"], "source": ["Theor Appl Genet"], "year": ["1999"], "volume": ["98"], "fpage": ["736"], "lpage": ["743"], "pub-id": ["10.1007/s001220051129"]}, {"surname": ["Torres", "Satovic", "Canovas", "Cobos", "Cubero"], "given-names": ["AM", "Z", "J", "S", "JI"], "article-title": ["Genetics and mapping of new isozyme loci in Vicia faba L using trisomics"], "source": ["Theor Appl Genet"], "year": ["1995"], "volume": ["91"], "fpage": ["783"], "lpage": ["789"], "pub-id": ["10.1007/BF00220960"]}, {"article-title": ["MultiQTL - the best QTL mapping software"]}, {"article-title": ["MapChart"]}, {"article-title": ["Medicago truncatula sequencing resources"]}, {"article-title": ["Grid Map"]}, {"surname": ["Purcell", "Plomin PR, DeFries PJ, McClearn PG, McGuffin PP"], "given-names": ["S"], "article-title": ["Statistical Methods in Quantitative Genetics"], "source": ["Behavioural Genetics"], "year": ["2001"], "edition": ["4th"], "publisher-name": ["New York , Worth Publishers"]}, {"surname": ["Arumuganathan", "Earle"], "given-names": ["K", "ED"], "article-title": ["Nuclear DNA content of some important plant species"], "source": ["Plant Mol Biol"], "year": ["1991"], "volume": ["9"], "fpage": ["208"], "lpage": ["218"], "pub-id": ["10.1007/BF02672069"]}, {"surname": ["Weeden", "Muehlbauer", "Ladizinsky"], "given-names": ["NF", "FJ", "G"], "article-title": ["Extensive conservation of linkage relationships between pea and lentil genetic maps"], "source": ["J Hered"], "year": ["1992"], "volume": ["83"], "fpage": ["123"], "lpage": ["129"]}, {"surname": ["Phan", "Ellwood", "Ford", "Thomas", "Oliver"], "given-names": ["HTT", "SR", "R", "S", "R"], "article-title": ["Differences in syntenic complexity between Medicago truncatula with Lens culinaris and Lupinus albus"], "source": ["Funct Plant Biol"], "year": ["2006"], "volume": ["33"], "fpage": ["775"], "lpage": ["782"], "pub-id": ["10.1071/FP06102"]}, {"surname": ["Murray", "Peters", "Thompson"], "given-names": ["MG", "DL", "WF"], "article-title": ["Ancient repeated sequences in the pea and mung bean genomes and implications for genome evolution"], "source": ["Journal of Molecular Evolution"], "year": ["1981"], "volume": ["17"], "fpage": ["31"], "lpage": ["42"], "pub-id": ["10.1007/BF01792422"]}]
BMC Genomics. 2008 Aug 9; 9:380