# A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.)

Anne Webb, Amanda Cottage, Thomas Wood, Khalil Khamassi, Douglas Hobbs, Krystyna Gostkiewicz, Mark White, Hamid Khazaei, Mohamed Ali, Daniel Street, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Anne Webb, Amanda Cottage, Thomas Wood, Khalil Khamassi, Douglas Hobbs, et al. A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.). Plant Biotechnology Journal, 2016, 14 (1), pp.177-185. 10.1111/pbi.12371. hal-02637592

HAL Id: hal-02637592

https://hal.inrae.fr/hal-02637592

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
*Plant Biotechnology Journal* (2016) **14**, pp. 177–185

doi: 10.1111/pbi.12371

# A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (*Vicia faba L.*)

Anne Webb<fnref n="1" />,†, Amanda Cottage<fnref n="1" />,†, Thomas Wood<fnref n="1" />, Khalil Khamassi<fnref n="1" />,‡, Douglas Hobbs<fnref n="1" />, Krystyna Gostkiewicz<fnref n="1" />, Mark White<fnref n="2" />, Hamid Khazaei<fnref n="3" />,‡, Mohamed Ali<fnref n="4" />,§, Daniel Street<fnref n="5" />, Gérard Duc<fnref n="6" />, Fred L. Stoddard<fnref n="3" />, Fouad Maalouf<fnref n="7" />, Francis C. Ogbonnaya<fnref n="7" />,#, Wolfgang Link<fnref n="4" />, Jane Thomas<fnref n="1" /> and Donal M O'Sullivan<fnref n="1" />,<fnref n="8" />,*

Received 8 December 2014;
revised 12 February 2015;
accepted 3 March 2015.

*Correspondence (Tel +44 118 3786729; fax +44 118 926 1244; email d.m.osullivan@reading.ac.uk)

†These authors made an equal contribution.

‡Present address: Plant science Department, University of Saskatchewan, Agriculture Building, 51 Campus Drive, S7N 5A8 Saskatoon, SK, Canada.

§Present address: Agronomy Department, College of Agriculture, Assiut University, Assiut 71526, Egypt

¶Present address: Grains Research and Development Corporation (GRDC) Level 4 East Building, 4 National Circuit, Barton, ACT 2600, Australia.

**Keywords:** legume, faba bean, KASP genotyping, single nucleotide polymorphism, synteny.

## Summary

Faba bean (*Vicia faba L.*) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely related crop species, lentil and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2-cM interval of Vf chromosome 2 which was highly colinear with Mt3. The obvious candidate gene from 78 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene controlling anthocyanin production in Medicago and resequencing of the Vf orthologue showed a putative causative deletion of the entire 5' end of the gene.

## Introduction

Faba bean is a globally significant grain legume, with 4–5 million tonnes annual production worldwide in recent years (FAO, 2014), providing the predominant affordable dietary source of protein to subsistence farmers and urban populations across North Africa, Horn of Africa (including Ethiopia) and the Middle and Near East as well as parts of China. The use of faba bean is deeply embedded in the cuisine and culture of these countries and the great agronomic and morphological diversity of the crop gene-pool has permitted it to become adapted to many agroecological settings and different end uses (Cubero, 1974; Duc et al., 2010). As a nitrogen-fixing legume, it has an important role in maintaining soil fertility in rotational systems (Jensen et al., 2010; Köpke and Nemecek, 2010). However, faba bean production and consumption are becoming increasingly geographically disjointed, with alarming reductions in output in precisely those

countries with rapidly growing populations that depend most heavily on faba beans as a dietary protein source. For example, Egypt's self-sufficiency in faba bean production fell from 96% to 52% between 1998 and 2009 due to widespread infestation of the Nile Delta region by *Orobanche crenata* Forsk. and the near abandonment of faba bean production in the middle Nile Valley region due to the prevalence of faba bean necrotic virus (Mahmoud Zeid, pers. commun.). Moreover, the almost unchecked spread of parasitic weeds of the genus *Orobanche* has led to abandonment for faba bean growing of entire provinces not only in Egypt but across the entire Mediterranean basin (Maalouf et al., 2009). These as well as other production constraints such as drought and extreme temperature sensitivity could potentially be addressed by more effective genomics-enabled genetic analysis/breeding strategies.

Despite the global significance of the crop and the urgency of researchable challenges it faces, faba bean genetic and genomic

<footnote n="1">¹National Institute of Agricultural Botany, Cambridge, UK</footnote>

<footnote n="2">²PGRO, Thornaugh, Peterborough, UK</footnote>

<footnote n="3">³Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland</footnote>

<footnote n="4">⁴Department of Crop Sciences, Georg-August-Universität, Göttingen, Germany</footnote>

<footnote n="5">⁵LGC Ltd (Genomics Division), Hoddesdon, Herts, UK</footnote>

<footnote n="6">⁶INRA, UMR1347 Agroécologie, Dijon, France</footnote>

<footnote n="7">⁷ICARDA, Terbol, Lebanon</footnote>

<footnote n="8">⁸School of Agriculture, Policy and Development, University of Reading, Whiteknights, UK</footnote>
resources have only started to be developed in the last few years, due in part at least to its outcrossing habit and large (approx. 13 Gb) genome size. Initially, mapping studies have been based on random amplified polymorphic DNA (RAPD) marker technology and delivered small numbers of cross-specific trait-linked sequence characterized amplified region (SCAR) markers (Avila et al., 2003, 2004; Diaz-Ruiz et al., 2010; Gutierrez et al., 2006, 2007). Later, the development of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) marker sets potentially more amenable to deployment in genetic studies and breeding have been reported (Ellwood et al., 2008; Gong et al., 2011; Kaur et al., 2012; Ma et al., 2011; Wang et al., 2011; Zeid et al., 2009). The recent development of the gene sequence-based maps of V. faba (Cruz-Izquierdo et al., 2012; Ellwood et al., 2008; El-Rodeny et al., 2014; Kaur et al., 2014a; Satovic et al., 2013) has been a big step forward for faba bean genetics. However, so far, gene-based linkage maps have either utilized multiple low-throughput marker systems to achieve assignment of most major linkage groups to one of the six chromosomes (Satovic et al., 2013) or generated multiple linkage groups per chromosome, mostly unassigned (El-Rodeny et al., 2014; Kaur et al., 2014a,b), and therefore, there is not yet a single genotyping platform for V. faba which permits genomewide alignment to a comprehensive consensus map providing coverage of all six chromosomes and supports systematic targeting of marker development to regions of interest exploiting synteny with sequenced legumes.

We previously demonstrated conversion of set of mapped cleaved amplified polymorphic sequence (CAPS) and SNAPshot assays to Kompetitive Allele Specific PCR (KASP) format in faba bean (Cottage et al., 2012). Here, we significantly increase the SNP coverage of the faba bean genome by discovering, validating and mapping a much larger number of SNP markers, and further pursue the objective of making all useful SNPs readily accessible by opting for the low-cost, easy-to-use, KASP genotyping platform. We find extensive colinearity of the V. faba linkage map with the sequenced model legume M. truncatula and demonstrate the potential for application of this new syntenically anchored SNP map in trait dissection and breeding via the identification of a candidate gene/causative polymorphism for white-coloured flowers.

## Results

### Transcriptome resequencing for SNP discovery

Pyrosequencing and de novo transcript assembly of RNA from seedling tissues of two contrasting faba bean lines – NV643-4, a partly inbred line derived from the white flowered, white hilum Polish variety 'Albus' and NV648-1, a derivative of the ICARDA Bean Pure Line 'BPL10', a coloured flower, black hilum inbred line – yielded 17 397 and 16 041 contigs >200 bp in length for NV643-4 and NV648-1, respectively. From amongst these 5562 and 5534 contigs of NV643-4 and NV648-1, respectively, aligned to contigs from the other V. faba line as well as to at least one putative gene of M. truncatula. These transcript assemblies are available from the GenBank Transcript Sequence Assembly (TSA)

repository under the accession names GARZ00000000 (NV643-4) and GASA00000000 (NV648-1). The versions described in this paper are the first versions, GARZ01000000 and GASA01 000000, respectively. Tables 1 and 2 show some metrics of the sequencing runs and assembly results.

### Information content of validated KASP SNP assays

Of 1001 selected target sequences (data not shown) containing putative NV643-4/NV648-1 SNPs, assays were successfully designed for 845 SNPs (for details of individual target sequences, Table S1). To provide continuity with previous work (Cottage et al., 2012) and a link to the previously published Vf6 × Vf27 maps (Cruz-Izquierdo et al., 2012; Ellwood et al., 2008), the 48 most informative SNP assays developed from the Vf6 × Vf27 map (Cottage et al., 2012) were included in the subsequent analyses. With the exception of the aforementioned set of 48 previously described SNPs, all SNPs described here were discovered by alignment of sequences from just a single pair of genotypes. Therefore, we wished to assess allele frequencies in a diverse germplasm sample so that assays could be scored for information content and to quantify any ascertainment bias towards rare alleles, peculiar to one or other of the 'discovery' genotypes. A second, overlapping objective was to screen a set of parents of mapping populations that already existed or are under development, so that users of these populations could quickly extend coverage of their linkage maps to encompass polymorphic SNPs from the new set. To fulfil these complementary purposes, a set of 37 inbred lines (the 'validation panel') was assembled and genotyped (Table 3). A total of 891 SNP assays were validated and assigned to one of four quality classes based on call rate and tightness of clustering of the three possible genotypic classes. Figure S1 shows raw genotyping data for two SNPs assigned to quality classes I and IV, respectively. For the remainder of analyses, only the 757 SNPs, henceforth referred to as 'NV648/NV643' SNPs, belonging to quality classes I and II are considered further. Allele calls for only 67.6% 512 of 757 of the high-quality 'NV648/NV643' SNP assays matched the SNP allele predicted from sequencing (Table S1). Of the 245 mismatched allele calls, 223 were associated with NV643, where unfortunately, individuals from distinct sublines of this accession had been used for transcriptome sequencing and crossing, respectively. For 241 of 245 of the class I-II SNPs where discovery and genotyped calls were discordant, both alleles were present in the validation panel despite the lack of the predicted NV643-3 (discovery sequence) allele in NV643-4 (genotyped) DNA, and therefore, we attribute the majority of mismatches in predicted versus observed allele

<table><caption>**Table 1** Statistics for 454 sequencing of V. faba 'Albus' and 'BPL10'</caption><thead><tr><td>Genotype</td><td>Tissue</td><td>Total reads</td><td>Average raw read length</td></tr></thead><tbody><tr><td>NV643-4 (Albus)</td><td>10-day-old seedling</td><td>788 918</td><td>525</td></tr><tr><td>NV648-1 (BPL10)</td><td>10-day-old seedling</td><td>751 023</td><td>527</td></tr></tbody></table>

<table><caption>**Table 2** Statistics for transcriptome assembly of V. faba 'Albus' and 'BPL10' in gsAssembler, implemented in Newbler v. 2.3 (Roche)</caption><thead><tr><td>Genotype</td><td>Total contigs</td><td>N50</td><td>Total coverage %</td><td>Contigs in TSA submission</td><td>GenBank TSA accession</td></tr></thead><tbody><tr><td>NV634-4</td><td>26 709</td><td>1152</td><td>90.10</td><td>17 397</td><td>GARZ00000000</td></tr><tr><td>NV648-1</td><td>21 591</td><td>1119</td><td>89.72</td><td>16 041</td><td>GASA00000000</td></tr></tbody></table>
<table><caption>**Table 3** Inbred lines of *V. faba*, used to validate panel of SNP markers. The original accessions were provided by National Institute of Agricultural Botany (NIAB), United Kingdom, The Georg-August-University, Göttingen (GAUG), Germany; the International Center for Agricultural Research in Dry Areas (ICARDA), Syria; the Instituto de Investigación y Formación Agroalimentaria (IFAPA), Spain; the University of Helsinki (UoH) Finland; and the Institut National de la Récherche Agronomique (INRA), France</caption><thead><tr><th>Accession code</th><th>Other reference</th><th>Donor</th><th>% Homozygosity</th><th>% Heterozygosity</th></tr></thead><tbody><tr><td>NV153-1</td><td>ig12658</td><td>NIAB</td><td>98.73</td><td>1.27</td></tr><tr><td>NV639-1</td><td>Hedin</td><td>GAUG</td><td>98.61</td><td>1.39</td></tr><tr><td>NV643-3</td><td>Albus</td><td>NIAB</td><td>97.85</td><td>2.15</td></tr><tr><td>NV644-1</td><td>Kasztelan</td><td>NIAB</td><td>92.56</td><td>7.44</td></tr><tr><td>NV648-2</td><td>BPL10</td><td>ICARDA</td><td>98.87</td><td>1.13</td></tr><tr><td>NV656-3</td><td>ig101942</td><td>NIAB</td><td>99.37</td><td>0.63</td></tr><tr><td>NV657-2</td><td>INRA-29H</td><td>IFAPA</td><td>99.11</td><td>0.89</td></tr><tr><td>NV658-2</td><td>CGN07715 cf-3</td><td>GAUG</td><td>99.24</td><td>0.76</td></tr><tr><td>NV662-1</td><td>VF136</td><td>IFAPA</td><td>97.22</td><td>2.78</td></tr><tr><td>NV713-1</td><td>Côte d'Or</td><td>GAUG</td><td>98.36</td><td>1.64</td></tr><tr><td>NV714-1</td><td>Hiverna</td><td>GAUG</td><td>99.25</td><td>0.75</td></tr><tr><td>NV715-1</td><td>Webo</td><td>GAUG</td><td>99.10</td><td>0.90</td></tr><tr><td>NV716-1</td><td>Wibo</td><td>GAUG</td><td>99.36</td><td>0.64</td></tr><tr><td>NV717-1</td><td>L79-79</td><td>GAUG</td><td>99.50</td><td>0.50</td></tr><tr><td>NV718-1</td><td>L977-88</td><td>GAUG</td><td>99.75</td><td>0.25</td></tr><tr><td>NV719-1</td><td>L979-S1</td><td>GAUG</td><td>99.62</td><td>0.38</td></tr><tr><td>NV720-1</td><td>Bourdon</td><td>GAUG</td><td>99.24</td><td>0.76</td></tr><tr><td>NV721-1</td><td>Arrisot</td><td>GAUG</td><td>98.99</td><td>1.01</td></tr><tr><td>NV722-1</td><td>Banner</td><td>GAUG</td><td>98.99</td><td>1.01</td></tr><tr><td>NV723-1</td><td>Bulldog</td><td>GAUG</td><td>99.12</td><td>0.88</td></tr><tr><td>NV724-1</td><td>Pietranera</td><td>GAUG</td><td>99.24</td><td>0.76</td></tr><tr><td>NV725-1</td><td>GIZA3-2</td><td>GAUG</td><td>99.36</td><td>0.64</td></tr><tr><td>NV726-1</td><td>GIZA402</td><td>GAUG</td><td>99.11</td><td>0.89</td></tr><tr><td>NV727-1</td><td>BPL4628</td><td>ICARDA</td><td>98.86</td><td>1.14</td></tr><tr><td>NV728-1</td><td>TW</td><td>ICARDA</td><td>99.11</td><td>0.89</td></tr><tr><td>NV729-1</td><td>VF6</td><td>ICARDA</td><td>98.58</td><td>1.42</td></tr><tr><td>NV730-1</td><td>ILB4347-4</td><td>ICARDA</td><td>99.24</td><td>0.76</td></tr><tr><td>NV731-1</td><td>ILB4347-3</td><td>ICARDA</td><td>99.50</td><td>0.50</td></tr><tr><td>NV732-1</td><td>BPL710</td><td>ICARDA</td><td>97.83</td><td>2.17</td></tr><tr><td>NV733-1</td><td>NA112</td><td>ICARDA</td><td>99.12</td><td>0.88</td></tr><tr><td>NV734-1</td><td>ILB938</td><td>UoH</td><td>99.62</td><td>0.38</td></tr><tr><td>NV735-1</td><td>Mélodie</td><td>UoH</td><td>98.99</td><td>1.01</td></tr><tr><td>NV736-1</td><td>Aurora</td><td>UoH</td><td>98.22</td><td>1.78</td></tr><tr><td>NV737-1</td><td>CRB285</td><td>INRA</td><td>97.11</td><td>2.89</td></tr><tr><td>NV738-1</td><td>CRB2516</td><td>INRA</td><td>99.21</td><td>0.79</td></tr><tr><td>NV739-1</td><td>CRB2702</td><td>INRA</td><td>94.63</td><td>5.37</td></tr><tr><td>NV740-1</td><td>CRB100107</td><td>INRA</td><td>99.60</td><td>0.40</td></tr></tbody></table>

calls to NV643 stock heterogeneity rather than sequencing errors or flaws in the assay design. Amongst the 512 class I-II SNPs where discovery sequence and SNP genotype were concordant, 438 were nonprivate with respect to the discovery genotypes as both alleles were observed amongst the 35 nondiscovery genotypes, but amongst the remaining 73 SNPs, the clear majority (60) of minor alleles were private to the NV648-1 discovery genotype, indicating that this line is the more 'exotic' of the discovery genotypes with respect to the diversity sampled in the validation panel. Figure S2 shows the genetic distance relationships between lines of the validation panel. As expected, BPL10 is on

the longest branch, reflecting the ascertainment bias of the SNP discovery. However, polymorphism information content (PIC) values were high with an average of 0.316 amongst the 656 nonprivate SNPs. To add SNPs discovered using different discovery genotypes and markers which could be polymorphic both in the previously reported Vf6 × Vf27 gene-based linkage map and maps created using the new SNPs described here, 40 of the markers from the Vf6 × Vf27 map (Cottage et al., 2014), henceforth referred to as 'Vf6/Vf27' SNPs were genotyped on the validation panel. The average PIC of this subset of SNPs was 0.329. Heterozygosity amongst the inbred lines using the combined set of 837 class I-II SNPs averaged 1.35%, consistent with the status of these as relatively pure research lines used for genetic analyses.

## Towards a consensus SNP-based genetic linkage map

The six mapping populations listed in Table 4 were used to construct first individual linkage maps and in a second step, a consensus map for *V. faba*. The framework for the consensus map was provided by the NV648-1 × NV643-4 F₂ linkage map, where 470 of the 757 'NV648/NV643' and 10 of the 'Vf6/Vf27' SNP markers could be placed in just six linkage groups. The other five maps were more fragmented with between 12 and 20 linkage groups and less well populated with between 109 and 231 SNPs. Nonetheless, 207 SNPs not placed in the framework NV648-1 × NV643-4 F₂ map were added to the consensus map through the integration of these additional populations. The final consensus map constructed using Merge Map (Wu et al., 2011), referred to as the Vf 2014 Consensus, consisted of 687 SNPs (653 'NV646/NV643' and 34 'Vf6/Vf27' SNPs) corresponding to 542 unique, that is non-co-segregating, loci on six linkage groups and spanning a total length of 1403.8 cM (Figure 1, Table S2). The average gap size in this map is 2.6 cM, and the largest single gap is 24.7 cM. We postulate that these six linkage groups correspond to the six physical chromosomes. Based on the distribution of 'Vf6/Vf27' SNPs from linkage groups previously assigned to physical chromosomes (Cruz-Izquierdo et al., 2012; Ruiz-Rodríguez et al., 2014), we have numbered these six linkage groups I-VI to reflect the physical assignments already made (Table S3).

## Colinearity of the faba bean genetic map with the sequenced *M. truncatula* genome and the first SNP map of *Lens culinaris*

Synteny between the 'Vf 2014 Consensus' linkage map and the *M. truncatula* Mt3.5 pseudomolecule assembly, and the first comprehensive SNP map for lentil (Sharpe et al., 2013) and *M. truncatula* was investigated by assigning orthologous links to reciprocal best BLAST hits between the Medicago genome, SNP-flanking sequence from *V. faba* and SNP-containing contigs assembled from several lines of lentil. The comparative map view shown in Figure S3 represents the full set of links thus obtained at a BLAST cut-off of $E^{-30}$ and portrays a relatively simple pattern of colinearity between the respective genomes. Mt1 is colinear in a single, virtually uninterrupted block with Vf Chr III and in lentil to the central part of Linkage group 1 (LC_1) which is thought to represent chromosome 1 of *L. culinaris*. Similarly, Mt2 appears colinear to Chr 1 in bean and Lc2 in lentil. Also, Mt3 and Mt6 chromosomes appear colinear to substantial blocks of Vf Chr II and I in faba bean and Lc3 and Lc2 in lentil, respectively. Mt7 is colinear to blocks of Vf Chr V and I in faba bean and Lc 8in lentil, and Mt8 is colinear with large blocks in Vf Chr VI and Chr IV in bean and smaller areas of Lc1, 4 and 7 in lentil. Large parts of Mt4 are syntenous to Vf Chr II,IV and
<table><caption>**Table 4** Mapping populations genotyped with SNP markers developed in this study</caption><thead><tr><th>Population</th><th>Male parent</th><th>Female parent</th><th>progeny analysed</th><th>mapped SNPs</th><th>linkage groups</th><th>map length (cM)</th><th>Segregating traits</th></tr></thead><tbody><tr><td>1</td><td>NV643-1</td><td>NV648-1</td><td>136 F2</td><td>481</td><td>6</td><td>1166.5</td><td>Flower colour (white/wild type)</td></tr><tr><td>2</td><td>NV643-1</td><td>NV657</td><td>165 F2</td><td>231</td><td>13</td><td>983.3</td><td>Flower colour (white/wild type)</td></tr><tr><td>3</td><td>NV639</td><td>NV658</td><td>52 F2</td><td>204</td><td>18</td><td>1198.9</td><td>Flower opening (closed/open)</td></tr><tr><td>4</td><td>NV644-1</td><td>NV153-1</td><td>125 F2</td><td>192</td><td>12</td><td>1164.9</td><td>Plant height (dwarf/tall)</td></tr><tr><td>5</td><td>Mélodie/2</td><td>ILB938/2</td><td>194 F5</td><td>200</td><td>13</td><td>818.2</td><td>Vicine-convicine (high/low), stipule spot (present/absent)</td></tr><tr><td>6</td><td>Côte D'Or/1</td><td>BPL4628/1521</td><td>101 RILs</td><td>108</td><td>20</td><td>158.1</td><td>Frost hardiness</td></tr></tbody></table>

VI and Lc4 and Lc7. This broad picture of high degree of conservation of gene order over few, extensive blocks of colinearity suggests that the gene content and order of *V. faba* genetic intervals currently lacking SNP coverage can be relatively straightforwardly and accurately predicted from aligned stretches of the *M. truncatula* reference genome sequence.

### Mapping of the ZT1 zero-tannin end-use quality trait

Condensed tannins have long been recognized as antinutritional factors when faba bean is used as the main protein source in feed formulations for monogastric livestock and poultry (Bond, 1976). Previous work has shown that the absence of the precursors of condensed tannins, the anthocyanins, in floral tissues, resulting in a highly recognizable white flower phenotype, is pleiotropic to the absence of condensed tannins in the seed coat (Bond, 1976; Picard, 1976). The zero-tannin (zt) trait as it is called, and has been genetically defined as a simple character, governed independently by two complementary recessive genes zt-1 and zt-2 (Cabrera and Martin, 1989; Picard, 1976). SCAR markers linked to the more common zt-1 source of zero tannin/white flower were identified by bulked segregant analysis (Gutierrez et al., 2007) but the resulting linkage group was not associated with any particular chromosome and the sequence of the markers gives no clue as to their location in the *V. faba* genome. The segregation of the white flower phenotype in the NV648 (♀) × NV643 (♂) F₂ population showed a simple 3 : 1 coloured-to-white ratio, indicating a single dominant gene was responsible, with the dominant allele coming from the coloured parent. The trait mapped to Vf chromosome 2, in a region of highly conserved synteny with Medicago chromosome 3, between markers Vf_Mt3g092810_001 and Vf_Mt3g094760_001 (Figure S4A and B).

### Use of synteny to target marker development to regions of interest and identify candidate genes

The Medicago chromosome segment lying between flanking markers Vf_Mt3g092810_001 and Vf_Mt3g094760_001 corresponds to approx. 0.57 Mb of DNA with 78 predicted gene models (Figure S4B). Most of these encode hypothetical proteins associated with various retroelements and just two of these were considered as positional candidates for regulating the expression of flower colour in *V. faba*.

A gene ortholog of the Medicago WD-40 transcription factors, Transparent Testa Glabra 1 (TTG1) (Medtr3g092830, Medtr3g092840), was identified in the interval between Vf_Mt3g92810_001 and Vf_Mt3g094760_001. Roles in master regulation of the anthocyanin biosynthesis pathway for TTG1 orthologs have been suggested in both Arabidopsis thaliana (Walker et al., 1999) and *M. truncatula* (Pang et al., 2009), whereas a deletion in the A2 gene, an ortholog of

Medtr3g092840 in *Pisum sativum*, has been demonstrated to be responsible for a loss of flower colour (Hellens et al., 2010). Their known conserved role in regulation of flower and seed anthocyanin production across Arabidopsis, *Medicago* and pea therefore made these compelling biological candidates for regulating anthocyanin production and thus flower pigmentation in *V. faba*. We therefore proceeded to investigate whether sequence variation at the Vf_Mt3g092840 locus, which we will henceforth refer to as *VfTTG1*, could be responsible for the loss of flower colour in NV643. The coding region of *VfTTG1* was cloned and sequenced from NV648 (Figure S4C). However, the 5' end of the *VfTTG1* sequence could not be amplified in the white-flowered NV643 line, suggesting a sequence difference or deletion affecting the 5' priming site. Genome walking was utilized to obtain the sequence of the 5' region of the *VfTTG1* locus in NV643. This showed that the entire 5' end of *VfTTG1* gene and upstream promoter elements were missing. The 5' flanking sequence instead showing high levels of homology to the *VfENOD11* gene (a 213 bp region) and an adjacent region of 388 bp with homology to the *PsiPT1* locus, encoding adenylate isopentenyltransferase) as illustrated in Figure S4D. Primers designed to amplify across the putative deletion amplified a fragment of the expected size (data not shown). An allelism test conducted by crossing NV643 to known sources of zt-1 and zt-2 allowed us to conclude that *VfTTG1* is in fact, ZT1 (data not shown). Taken together, the allelic variation, positional and functional data from related taxa strongly indicate that a deletion in the coding region of *VfTTG1/ZT1* in *V. faba* is responsible for the loss of floral pigmentation in NV643 and demonstrates how translational functional genomics exploiting well-described synteny between *Medicago* and *V. faba* can be used to quickly and accurately identify causative genes contributing to specific phenotypes in a crop.

## Discussion

### Emergence of *V. faba* from 'orphan legume' status

In the space of just over 1 year, faba bean entries in GenBank have gone from 5000 EST reads (Ray and Georges, 2010) to multiple transcriptomes (Kaur et al., 2014a,b; this study). Like pigeon pea, groundnut, lentil and other legumes, there is now a linkage map available with a density of marker coverage appropriate to most currently available RIL mapping populations.

Legume crop species now enjoying low-medium density SNP coverage of their genomes include cowpea (Muchero et al., 2009), common bean (Cortés et al., 2011), chickpea (Hiremath et al., 2012), lentil (Sharpe et al., 2013), pea (Deulvot et al., 2010) and faba bean (this work).
<figure><img src="image_3.png" /><figcaption>**Figure 1** A 687-locus consensus SNP-based map showing the six linkage groups thought to represent the six physical chromosomes of *Vicia faba*. Markers in red have been developed in previous studies. Graph created using MapChart (Voorrips, 2002).</figcaption></figure>

## Uses and usability of *Vicia faba* SNP markers in marker-assisted breeding

Exploitation of molecular markers depends on many factors: their information content, entry cost, per data point cost,

infrastructural requirements, ability to access information and coherence as a set of markers likely to provide reasonable genome coverage. We consider that we have attempted to address every single one of these bottlenecks. Our data show that PIC values are high, presumably as a consequence of high
outcrossing rate, which suggests that these markers will be broadly useable in each of the main recognized seed size based gene pools ('major', 'equina' and 'minor'). The KASP competitive PCR fluorescence assay chemistry was chosen with two particular criteria in mind: firstly, minimizing the cost to researchers of conducting low- to medium-resolution QTL scans and secondly, to cater for marker-assisted selection in real-world breeding programmes, where the cost per progeny eliminated on the basis of marker genotypes must be just pennies, ruling out any currently available genomewide genotyping technologies. Many mapping populations are available in the community and we hope that the lower cost barriers will make it possible for more groups to be involved in trait mapping than are currently. For example, SNP markers from this study found to be polymorphic in population 5 (Table 4) have been used to map QTL for stomatal traits (Khazaei et al., 2014a), low vicine-convicine (Khamassi et al., 2013; Khazaei et al., 2015) and a novel anthocyanin pigmentation locus, *ssp1* (Khazaei et al., 2014b). Similarly, SNP markers from this study used to generate the linkage map for population 4 (Table 4) have been found linked to a novel dwarfing locus, *Dwf1* (K. Khamassi, K. Gostkiewicz, and D.M. O'Sullivan, Unpublished data) and associated with frost tolerance (Sallam, 2014). Furthermore, to facilitate redesign of individual markers to suit existing laboratory set-ups and locations not easily served by KASP chemistry and commercial genotyping services, an interactive database-driven website (O'Sullivan and Green, 2012) has been set up to allow users to click from a marker name/location through to a choice of viewing or downloading the original sequence alignment on which the SNP assay design was based as well as accessing clustering profiles, genotype scores and links to any associated references.

These data can be used to pick informative sets of SNPs discovered from four diverse genetic backgrounds and used to efficiently capture genomewide genetic relationships, measure diversity, map traits and once mapped, use marker-assisted selection schemes to breed for those traits.

### Synteny-based mapping

Our results broadly support, but refine and extend the conclusions of Cruz-Izquierdo *et al.* (2012) and Kaur *et al.* (2014a) regarding synteny between the eight *M. truncatula* chromosomes and the six *V. faba* chromosomes. The simple pattern of few, large blocks of colinearity between *V. faba* and *M. truncatula* reported here therefore greatly simplifies the task of translational functional genomics.

Synteny-based candidate gene discovery as practised for example in cowpea (Muchero *et al.*, 2009) is now becoming a reality in faba bean. The identification of a candidate for the zero-tannin trait illustrates this very well as the SNP markers most closely linked to the zero-tannin locus occur within a large block of colinearity between Mt3 and Vf LG2 in which a plausible candidate gene is found.

### Towards isolation of causative genes for important trait variation

Tannins (condensed polyphenols) are considered antinutritional compounds limiting the inclusion rates of faba bean in some animal feed formulations. It has been established for a long time that tannins are absent from the testae of white-flowered bean lines, both phenotypes being manifestations of the absence of anthocyanin production. Recessive alleles at two separate com-

plementary loci (*zt1* and *zt2*) act independently to abolish anthocyanin production, testa tannins and flower colour markings. RAPD-based SCAR markers linked to *zt1* were previously isolated using a bulked segregant analysis approach although these SCARs were not diagnostic for white flower colour in a panel of 37 European inbred lines and the linkage groups produced in this study did not contain any markers appearing on subsequent *Vicia faba* linkage maps (Cruz-Izquierdo *et al.*, 2012). Our results show that when the zero-tannin trait was placed on this syntenically anchored map, a rapid inspection of annotated gene content of the colinear region of the Mt3.5 genome revealed tandem copies of the WD40-repeat containing transcription factor, whose *Arabidopsis* orthologue had been functionally characterized as the *TRANSPARENT TESTA GLABRA1* gene (Pang *et al.*, 2009; Qi *et al.*, 2011; Sagasser *et al.*, 2002; Walker *et al.*, 1999; Zhao *et al.*, 2008). Although a causative polymorphism controlling determinate flowering trait in faba bean was previously isolated by amplification of the presumed functional ortholog of the *TFL1* gene (Avila *et al.*, 2007), this is, to our knowledge, the first time in faba bean that it has been possible to pinpoint a plausible causal polymorphism by combining a syntenically anchored SNP consensus map and sequence and gene function information from Medicago. The syntenic framework presented here provides a basis for the reconstruction of predicted Vf gene order encompassing the large majority of genes that fall within genetic intervals of nearly uninterrupted synteny. A 'genome zipper' of this nature is an essential stepping stone to exploit the power of population-based and bulk segregant mapping-by-sequencing approaches (Gardiner *et al.*, 2014).

## Experimental procedures

### Plant material

The faba bean inbred lines used in this study either for transcriptome resequencing or SNP validation and their sources are listed in Table 3, as well as all mapping population parents. Genetic linkage mapping was conducted using six biparental populations, further details of which are presented in Table 4.

### RNA library preparation and sequencing

RNA was extracted from complete 7-day-old seedlings of part inbred line NV643-4 (Albus) and inbred line NV648-1 (BPL10) using the Qiagen RNeasy Plant Mini Kit, and sequenced using Roche 454 GS-FLX Titanium chemistry. Raw reads were assembled in gsAssembler, based on Newbler v2.3, available free from Roche Life Sciences at http://454.com/products/analysis-software/index.asp (Anonymous, 2010).

CDS, whole transcript and protein sequences from *M. truncatula* were obtained online from the *M. truncatula* Genome Sequencing Project, led by the International Medicago Genome Annotation Group (IMGAG). Version Mt3.5 (Data freeze: 12/31/2009, Assembly release date: 02/10/2010, annotation release date: 05/27/2010) of the annotated genome of *M. truncatula*, consisting of 8 chromosomes and some unanchored BACs was downloaded from http://www.medicagohapmap.org/downloads.php on 01 December 2010.

### SNP discovery and validation

Assembled contigs from NV643-4 and NV648-1 were aligned to each other and to transcript sequences of *M. truncatula* using
BLASTn. All contigs which gave unique reciprocal best hits over at least 100 bases at a BLASTn cut-off of $E^{-30}$ in the respective faba bean transcriptome assemblies and to a single predicted *M. truncatula* gene over more than 100 bases were selected for further analysis.

The complete contigs of each NV643-4:NV648-1 pair were aligned using ClustalW (Chenna et al., 2003) and alignments containing SNPs more than 50 bp from either end of a contig were flagged for marker development. For each SNP, a consensus sequence with the target SNP and at least 50 bp either side of it was then searched against the database of *M. truncatula* pseudomolecules v.3.5, using the BLASTn function (http://www.medicagohapmap.org/advanced_search_page.php?seq). Any consensus which matched multiple chromosomes of *M. truncatula* or in multiple copies on the same chromosome or whose *M. truncatula* alignments contained introns located closer than 50 bp of either side of the SNP were excluded. Sequences with more than 50 bp between the target SNP and the nearest intron site on the corresponding *M. truncatula* contig were truncated at the nearest intron-exon junctions to avoid primers being designed across intron-exon splice sites which would not anneal to genomic (intron-containing) sequence. The remaining SNP-containing sequences were then searched in BLASTn against all contigs of NV643-4 and NV648-1. Any SNP-containing alignments with hits to more than one contig from each line were excluded. A total of 887 KASP probes were designed by LGC Genomics Ltd. and each assay given an internal LGC Genomics ID as well as an assay name following the format *Vf_MtXgYYYYYY_ZZZ* where X indicates the Medicago chromosome number, YYYYYYY the predicted gene identifier of the reciprocal best BLAST hit *M. truncatula* orthologue as given in version Mt3.5 of the Medicago genome assembly (Young et al., 2011), and ZZZ the number of the assay (Table S1).

A panel of 37 lines of *V. faba* (Table 3) chosen to represent the parents/founders of key mapping population resources was assembled and DNA extracted from a seedling leaf of a single plant per line using a modified Tanksley method (Fulton et al., 1995). Designed KASP assays were genotyped across this panel according to the manufacturer's instructions.

## SNP-marker quality filtering

Clustering of alternate homozygous SNP alleles and heterozygous genotypes was manually inspected using LGC Genomics SNP viewer. Each SNP assay was scored on a 1–IV scale according to tightness of clustering and the number of missing data calls. As an additional criterion on which to rank assays, gene diversity, often referred to as expected heterozygosity, a measure of how informative individual assays are, was calculated as 2pq where *p* is the frequency of the first allele and *q* (the frequency of the alternate allele). Lines or markers with more than 6% no-calls were excluded in an iterative process. Markers were regarded as very skewed and removed, if either one of the homozygous parental alleles was absent or the homozygous to heterozygous ratio was >2.1 or <0.4 for F2 populations 1–4 and >0.25 or <0.045 for populations 5 and 6, respectively (Table 4).

## Linkage mapping

MapDisto 1.7.5.1 for Windows (Excel 2007 version) was used to find and refine linkage groups (Lorieux, 2012). The LOD was set to five and the map was generated using the procedure

Automap. For calculation of all individual linkage maps, the maximum permitted recombination fraction, the parameter *r* in MapDisto, was set to 0.3, except for population four where the use of a more stringent setting of 0.2 was required to avoid spurious linkage. The consensus map was generated in MergeMap (Wu et al., 2008) from the linkage groups generated in MapDisto. Graphs were produced using MapChart (Voorrips, 2002).

## Synteny analysis

The flanking sequences of all mapped SNP markers and all lentil contigs were searched in BLASTn against *M. truncatula* CDS version Mt3.5 with an E-value cut-off $E^{-37}$. The results were displayed in Strudel (Bayer et al., 2011). Discontinuous megablast was used to detect sequence homologies in the flanking region of the recessive VfWD40-1 allele from NV643.

## Trait targeting by synteny

The coding sequence of Vf_Mt3g092840 (*VfTTG1*) was amplified using standard PCR from NV648 using oligonucleotides: VfTTG1.F1: 5'-ATGAGATCTAAAACTACGCCTGTGG/VfTTG1.R1: 5'-TCAAACTCGAACCCTCAAAAGC (0.5 µM) using Phusion taq polymerase (N.E.B., Hitchin, UK). Thermal cycling conditions: 96 °C for 5 min, (30 × 96 °C/15 s, 58 °C/15 s, 72 °C/30 s) and 72 °C/10 min final extension. Products were cloned into a pGEM-T Easy vector (Promega, Southampton, UK) and sequenced using BigDye V3.0 (Life Technologies, Paisley, UK).

A genome walking kit (Clontech, Mountain View (CA), USA) was used in conjunction with the VfTTG1.R1 primer to generate primary amplification products for walking. A nested primer (VfTTG1.N) 5'-TCAAACCCTCAAAAGCTGCATCTTGTTAG was then used for walking upstream generating a 982-bp product. This was cloned and sequenced as described previously.

## Acknowledgements

We would like to acknowledge funding from the Defra Pulse Crop Genetic Improvement Network and Technology Strategy Board grant TS/001301 to DOS and JET. We are grateful to Peter Smith of Wherry & Sons Ltd. for access to mapping populations and to Melanie Febrer and Jane Rogers of The Genome Analysis Centre (TGAC), Norwich for generating the NV643 and NV648 transcriptome sequence data sets.

## References

* Anonymous (2010) 454 Sequencing System Software. 454 Life Sciences Corp., A Roche Company, Branford, CT 06405, USA. http://454.com/products/analysis-software/index.asp.

* Avila, C.M., Sillero, J.C., Rubiales, D., Moreno, M.T. and Torres, A.M. (2003) Identification of RAPD markers linked to the *Uvf-1* gene conferring hypersensitive resistance against rust (*Uromyces viciae-fabae*) in *Vicia faba* L. *Theor. Appl. Genet.* **107**, 353–358.

* Avila, C.M., Satovic, Z., Sillero, J.C., Rubiales, D., Moreno, M.T. and Torres, A.M. (2004) Isolate and organ-specific QTLs for ascochyta blight resistance in faba bean (*Vicia faba* L). *Theor. Appl. Genet.* **108**, 1071–1078.

* Avila, C., Atienza, S., Moreno, M. and Torres, A. (2007) Development of a new diagnostic marker for growth habit selection in faba bean (*Vicia faba* L.) breeding. *Theor. Appl. Genet.* **115**, 1075–1082.

* Bayer, M., Milne, I., Stephen, G., Shaw, P., Cardle, L., Wright, F. and Marshall, D. (2011) Comparative visualization of genetic and physical maps with Strudel. *Bioinformatics*, **27**, 1307–1308.
Bond, D.A. (1976) In vitro digestibility of the testa in tannin-free field beans (*Vicia faba* L). *J. Agric. Sci.* **86**, 561-566.

Cabrera, A. and Martin, A. (1989) Genetics of tannin content and its relationship with flower and testa colours in *Vicia faba* L. *J. Agric. Sci.* **113**, 93-98.

Chenna, R., Sugawara, H., Koike, T., Lopez, R., Gibson, T.J., Higgins, D.G. and Thompson, J.D. (2003) Multiple sequence alignment with the Clustal series of programs. *Nucleic Acids Res.* **31**, 3497–3500.

Cortés, A., Chavarro, M. and Blair, M. (2011) SNP marker diversity in common bean (*Phaseolus vulgaris* L.). *Theor. Appl. Genet.* **123**, 827–845.

Cottage, A., Gostkiewicz, K., Thomas, J., Borrows, R., Torres, A.M. and O'Sullivan, D. (2012) Heterozygosity and diversity analysis using mapped single nucleotide polymorphisms in a faba bean inbreeding programme. *Mol. Breeding*, **30**, 1799–1809.

Cruz-Izquierdo, S., Avila, C.M., Satovic, Z., Palomino, C., Gutierrez, N., Ellwood, S.R., Phan, H.T., Cubero, J.I. and Torres, A.M. (2012) 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.* **125**, 1767–1782.

Cubero, J.I. (1974) On the evolution of *Vicia faba* L. *Theor. Appl. Genet.* **45**, 47–51.

Deulvot, C., Charrel, H., Marty, A., Jacquin, F., Donnadieu, C., Lejeune-Henaut, I., Burstin, J. and Aubert, G. (2010) Highly-multiplexed SNP genotyping for genetic mapping and germplasm diversity studies in pea. *BMC Genom.* **11**, 468.

Diaz-Ruiz, R., Torres, A.M., Satovic, Z., Gutierrez, M.V., Cubero, J.I. and Roman, B. (2010) Validation of QTLs for *Orobanche crenata* resistance in faba bean (*Vicia faba* L.) across environments and generations. *Theor. Appl. Genet.* **120**, 909–919.

Duc, G., Bao, S., Baum, M., Redden, B., Sadiki, M., Suso, M.J., Vishniakova, M. and Zong, X. (2010) Diversity maintenance and use of *Vicia faba* L. genetic resources. *Field. Crop. Res.* **115**, 270–278.

Ellwood, S., Phan, H., Jordan, M., Hane, J., Torres, A., Avila, C., Cubero-Izquierdo, S. and Oliver, R. (2008) Construction of a comparative genetic map in faba bean (*Vicia faba* L.); conservation of genome structure with *Lens culinaris*. *BMC Genet.* **9**, 380.

El-Rodeny, W., Kimura, M., Hirakawa, H., Sabah, A., Shirasawa, K., Sato, S., Tabata, S., Sasamoto, S., Watanabe, A., Kawashima, K., Kato, M., Wada, T., Tsuruoka, H., Takahashi, C., Minami, C., Nanri, K., Nakayama, S., Kohara, M., Yamada, M., Kishida, Y., Fujishiro, T. and Isobe, S. (2014) Development of EST-SSR markers and construction of a linkage map in faba bean (*Vicia faba*). *Breed. Sci.* **64**, 252–263.

FAO (2014) FAOSTAT. Food and agriculture organisation of the United Nations http://faostat.fao.org/faostat.

Fulton, T.M., Chunwongse, J. and Tanksley, S.D. (1995) Microprop protocol for extraction of DNA from tomato and other herbaceous plants. *Plant Mol. Biol. Rep.* **13**, 207–209.

Gardiner, L.-J., Gawroński, P., Olohan, L., Schnurbusch, T., Hall, N. and Hall, A. (2014) Using genic sequence capture in combination with a syntenic pseudo genome to map a deletion mutant in a wheat species. *Plant J.* **80**, 895–904.

Gong, Y.-M., Xu, S.-C., Mao, W.-H., Hu, Q.-Z., Zhang, G.-W., Ding, J. and Li, Z.-Y. (2011) Generation and characterization of 11 novel EST derived microsatellites from *Vicia faba* (Fabaceae). *Am. J. Bot.* **98**, e22–e24.

Gutierrez, N., Avila, C.M., Duc, G., Marget, P., Suso, M.J., Moreno, M.T. and Torres, A.M. (2006) CAPs markers to assist selection for low vicine and convicine contents in faba bean (*Vicia faba* L.). *Theor. Appl. Genet.* **114**, 59–66.

Gutierrez, N., Avila, C., Rodriguez-Suarez, C., Moreno, M. and Torres, A. (2007) Development of SCAR markers linked to a gene controlling absence of tannins in faba bean. *Mol. Breeding*, **19**, 305–314.

Hellens, R.P., Moreau, C., Lin-Wang, K., Schwinn, K.E., Thomson, S.J., Fiers, M.W.E.J., Frew, T.J., Murray, S.R., Hofer, J.M.I., Jacobs, J.M.E., Davies, K.M., Allan, A.C., Bendahmane, A., Coyne, C.J., Timmerman-Vaughan, G.M. and Ellis, T.H.N. (2010) Identification of Mendel's White Flower Character. *PLoS ONE*, **5**, e13230.

Hiremath, P.J., Kumar, A., Penmetsa, R.V., Farmer, A., Schlueter, J.A., Chamarthi, S.K., Whaley, A.M., Carrasquilla-Garcia, N., Gaur, P.M., Upadhyaya, H.D., Kavi Kishor, P.B., Shah, T.M., Cook, D.R. and Varshney, R.K. (2012) Large-scale development of cost-effective SNP marker assays for diversity assessment and genetic mapping in chickpea and comparative mapping in legumes. *Plant Biotechnol. J.* **10**, 716–732.

Jensen, E.S., Peoples, M.B. and Hauggaard-Nielsen, H. (2010) Faba bean in cropping systems. *Field. Crop. Res.* **115**, 203–216.

Kaur, S., Pembleton, L.W., Cogan, N.O., Savin, K.W., Leonforte, T., Paull, J., Materne, M. and Forster, J.W. (2012) Transcriptome sequencing of field pea and faba bean for discovery and validation of SSR genetic markers. *BMC Genom.* **13**, 104.

Kaur, S., Kimber, R.B.E., Cogan, N.O.I., Materne, M., Forster, J.W. and Paull, J.G. (2014a) SNP discovery and high-density genetic mapping in faba bean (*Vicia faba* L.) permits identification of QTLs for ascochyta blight resistance. *Plant Sci.* **217–218**, 47–55.

Kaur, S., Cogan, N., Forster, J. and Paull, J. (2014b) Assessment of genetic diversity in faba bean based on single nucleotide polymorphism. *Diversity*, **6**, 88–101.

Khamassi, K., Ben Jeddi, F., Hobbs, D., Irigoyen, J., Stoddard, F.L., O'Sullivan, D.M. and Jones, H. (2013) A baseline study of vicine-convicine levels in faba bean (*Vicia faba* L.) germplasm. *Plant Genet. Resour.* **11**, 250–257.

Khazaei, H., O'Sullivan D, M., Sillanpää, M.J. and Stoddard, F. (2014) Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (*Vicia faba* L.). *Theor. Appl. Genet.* **127**, 2371–2385.

Khazaei, H., O'Sullivan, D.M., Jones, H., Pitts, N., Sillanpää, M.J., Pärssinen, P., Manninen, O. and Stoddard, F.L. (2015) Flanking SNP markers for vicine-convicine concentration in faba bean (*Vicia faba* L.). *Mol. Breeding*, **35**, 38.

Khazaei, H., O'Sullivan, D.M., Sillanpää, M. and Stoddard, F.L. (2014b) Genetic analysis reveals a novel locus in *Vicia faba* decoupling pigmentation in the flower from that in the extra-floral nectaries. *Mol. Breeding*, **34**, 1507–1513.

Köpke, U. and Nemecek, T. (2010) Ecological services of faba bean. *Field. Crop. Res.* **115**, 217–233.

Lorieux, M. (2012) MapDisto: fast and efficient computation of genetic linkage maps. *Mol. Breeding*, **30**, 1231–1235.

Ma, Y., Yang, T., Guan, J., Wang, S., Wang, H., Sun, X. and Zong, X. (2011) Development and characterization of 21 EST-derived microsatellite markers in *Vicia faba* (fava bean). *Am. J. Bot.* **98**, e22–e24.

Maalouf, F., Hamdi, A., Cubero, J.I., Khalifa, G.E., Jarsso, M., Kemal, S. and Karajeh, F. (2009) *Summary review of faba bean research in the Nile Valley, Red Sea and sub-Saharan region, and future recommendations: ICARDA*.

Muchero, W., Diop, N.N., Bhat, P.R., Fenton, R.D., Wanamaker, S., Pottorff, M., Hearne, S., Cisse, N., Fatokun, C., Ehlers, J.D., Roberts, P.A. and Close, T.J. (2009) A consensus genetic map of cowpea [Vigna unguiculata (L) Walp.] and synteny based on EST-derived SNPs. *Proc. Natl Acad. Sci. USA*, **106**, 18159–18164.

O'Sullivan, D.M. and Green, S. (2012) *The Vicia Toolbox*. www.viciatoolbox.org.

Pang, Y., Wenger, J.P., Saathoff, K., Peel, G.J., Wen, J., Huhman, D., Allen, S.N., Tang, Y., Cheng, X., Tadege, M., Ratet, P., Mysore, K.S., Sumner, L.W., Marks, M.D. and Dixon, R.A. (2009) A WD40 repeat protein from *Medicago truncatula* is necessary for tissue-specific anthocyanin and proanthocyanidin biosynthesis but not for trichome development. *Plant Physiol.* **151**, 1114–1129.

Picard, J. (1976) Aperu sur l'herédite du caractère absence de tannins dans les graines de feverole (*Vicia faba* L.). *Annales de l'Amelioration des Plantes*, **26**, 101–106.

Qi, T., Song, S., Ren, Q., Wu, D., Huang, H., Chen, Y., Fan, M., Peng, W., Ren, C. and Xie, D. (2011) The Jasmonate-ZIM-domain proteins interact with the WD-repeat/bHLH/MYB complexes to regulate jasmonate-mediated anthocyanin accumulation and trichome initiation in *Arabidopsis thaliana*. *Plant Cell*, **23**, 1795–1814.

Ray, H. and Georges, F. (2010) A genomic approach to nutritional, pharmacological and genetic issues of faba bean (*Vicia faba*): prospects for genetic modifications. *GM Crops*, **2010**, 99–106.

Ruiz-Rodriguez, M.D., Avila, C.M., Torres, A.M., Fuchs, J. and Schubert, I. (2014) Anchoring of genetic linkage maps to the chromosome complement of *Vicia faba* L. *Mol. Breeding*, **33**, 743–748.

Sagasser, M., Lu, G.H., Hahlbrock, K. and Weisshaar, B. (2002) A-thaliana TRANSPARENT TESTA 1 is involved in seed coat development and defines the WIP subfamily of plant zinc finger proteins. *Genes Dev.* **16**, 138–149.

Sallam, A. (2014) *Detailed genetic approach to improve frost tolerance of German winter faba beans*. PhD thesis, Göttingen: Georg-August Universität.
Satovic, Z., Avila, C.M., Cruz-Izquierdo, S., Diaz-Ruiz, R., Garcia-Ruiz, G.M., Palomino, C., Gutierrez, N., Vitale, S., Ocana-Moral, S., Gutierrez, M.V., Cubero, J.I. and Torres, A.M. (2013) A reference consensus genetic map for molecular markers and economically important traits in faba bean (*Vicia faba* L.). *BMC Genom.* **14**, 932.

Sharpe, A.G., Ramsay, L., Sanderson, L.A., Fedoruk, M.J., Clarke, W.E., Li, R., Kagale, S., Vijayan, P., Vandenberg, A. and Bett, K.E. (2013) Ancient orphan crop joins modern era: gene-based SNP discovery and mapping in lentil. *BMC Genom.* **14**, 192.

Voorrips, R.E. (2002) MapChart: Software for the graphical presentation of linkage maps and QTLs. *J. Hered.* **93**, 77–78.

Walker, A.R., Davison, P.A., Bolognesi-Winfield, A.C., James, C.M., Srinivasan, N., Blundell, T.L., Esch, J.J., Marks, M.D. and Gray, J.C. (1999) The TRANSPARENT TESTA GLABRA1 locus, which regulates trichome differentiation and anthocyanin biosynthesis in arabidopsis, encodes a WD40 repeat protein. *Plant Cell*, **11**, 1337–1349.

Wang, H., Yang, T., Guan, J., Ma, Y., Sun, X. and Zong, X. (2011) Development and characterization of 20 novel polymorphic STS markers in *Vicia faba* (fava bean). *Am. J. Bot.* **98**, e189–e191.

Wu, Y., Close, T.J. and Lonardi, S. (2008) On the accurate construction of consensus genetic maps. *Comput. Syst. Bioinformatics Conf.* **7**, 285–296.

Wu, Y., Close, T.J. and Lonardi, S. (2011) Accurate construction of consensus genetic maps via integer linear programming. *IEEE/ACM Trans. Comput. Biol. Bioinform.* **8**, 381–394.

Young, N.D., Debelle, F., Oldroyd, G.E.D., Geurts, R., Cannon, S.B., Udvardi, M.K., Benedito, V.A., Mayer, K.F.X., Gouzy, J., Schoof, H., Van de Peer, Y., Proost, S., Cook, D.R., Meyers, B.C., Spannagl, M., Cheung, F., De Mita, S., Krishnakumar, V., Gundlach, H., Zhou, S., Mudge, J., Bharti, A.K., Murray, J.D., Naoumkina, M.A., Rosen, B., Silverstein, K.A.T., Tang, H., Rombauts, S., Zhao, P.X., Zhou, P., Barbe, V., Bardou, P., Bechner, M., Bellec, A., Berger, A., Berges, H., Bidwell, S., Bisseling, T., Choisne, N., Couloux, A., Denny, R., Deshpande, S., Dai, X., Doyle, J.J., Dudez, A.-M., Farmer, A.D., Fouteau, S., Franken, C., Gibelin, C., Gish, J., Goldstein, S., Gonzalez, A.J., Green, P.J., Hallab, A., Hartog, M., Hua, A., Humphray, S.J., Jeong, D.-H., Jing, Y., Jocker, A., Kenton, S.M., Kim, D.-J., Klee, K., Lai, H., Lang, C., Lin, S., Macmil, S.L., Magdelenat, G., Matthews, L., McCorrison, J., Monaghan, E.L., Mun, J.-H., Najar, F.Z., Nicholson, C., Noirot, C., O'Bleness, M., Paule, C.R., Poulain, J., Prion, F., Qin, B., Qu, C., Retzel, E.F., Riddle, C., Sallet, E., Samain, S., Samson, N., Sanders, I., Saurat, O., Scarpelli, C., Schiex, T., Segurens, B., Severin, A.J., Sherrier, D.J., Shi, R., Sims, S., Singer, S.R., Sinharoy, S., Sterck, L., Viollet, A., Wang, B.-B., Wang, K., Wang, M., Wang, X., Warfsmann, J., Weissenbach, J., White, D.D., White, J.D., Wiley, G.B., Wincker,

P., Xing, Y., Yang, L., Yao, Z., Ying, F., Zhai, J., Zhou, L., Zuber, A., Denarie, J., Dixon, R.A., May, G.D., Schwartz, D.C., Rogers, J., Quetier, F., Town, C.D. and Roe, B.A. (2011) The Medicago genome provides insight into the evolution of rhizobial symbioses. *Nature*, **480**, 520–524.

Zeid, M., Mitchell, S., Link, W., Carter, M., Nawar, A., Fulton, T. and Kresovich, S. (2009) Simple sequence repeats (SSRs) in faba bean: new loci from *Orobanche*-resistant cultivar 'Giza 402'. *Plant Breeding*, **128**, 149–155.

Zhao, M., Morohashi, K., Hatlestad, G., Grotewold, E. and Lloyd, A. (2008) The TTG1-bHLH-MYB complex controls trichome cell fate and patterning through direct targeting of regulatory loci. *Development*, **135**, 1991–1999.

## Supporting information

Additional Supporting information may be found in the online version of this article:

**Figure S1** Examples of SNP calls for the SNP validation panel of 37 inbred *V. faba* lines (Table 3) displayed in SNPViewer (LGC Ltd., UK) for a marker assigned to quality class I (top) and IV (bottom).

**Figure S2** UPGMA dendrogram of the SNP-validation panel of 37 inbred *V. faba* lines, based on the panel of SNP markers developed in this study.

**Figure S3** Display of synteny between chromosomes of *M. truncatula* (centre) and linkage groups likely to represent chromosomes of *V. faba* (left) and *L. culinaris* (Sharpe et al., 2013) (right).

**Figure S4** A portion of the linkage map from the $F_2$ population 1 (NV643-1 × NV648-1) showing the map location of ZT1 between flanking markers Mt3g094760_001 and Vf_Mt3g092810_001.

**Table S1** SNPs identified in *V. faba* lines Albus and BPL10, their flanking sequences and alleles, quality score and PIC values.

**Table S2** Consensus linkage map created in MergeMap from separate individual maps listed in Table S4.

**Table S3** Markers mapped here and in previous studies and their assignment to linkage groups of three maps created for *V. faba*.

**Table S4** Linkage groups from individual maps used to produce the consensus map in this study. Individual maps from populations listed in Table 2 were produced in MapDisto and given a weighting factor.