Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity¹. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value². Faba bean (Vicia fabaL.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

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 (Vicia faba), using an F6 population developed from a cross between the lines Vf6 (equina type) and Vf27 (paucijuga type).

Orobanche crenata Forsk. is a root parasite that produces devastating effects on many crop legumes and has become a limiting factor for faba bean production in the Mediterranean region. The efficacy of available control methods is minimal and breeding for broomrape resistance remains the most promising method of control. Resistance seems to be scarce and complex in nature, being a quantitative characteristic difficult to manage in breeding programmes. To identify and map the QTLs (quantitative trait loci) controlling the trait, 196 F₂plants derived from the cross between a susceptible and a resistant parent were analysed using isozymes, RAPD, seed protein genes, and microsatellites. F₂- derived F₃lines were studied for broomrape resistance under field conditions. Of the 130 marker loci segregating in the F₂population, 121 could be mapped into 16 linkage groups. Simple interval mapping (SIM) and composite interval mapping (CIM) were performed using QTL Cartographer. Composite interval mapping using the maximum number of markers as cofactors was clearly the most efficient way to locate putative QTLs. Three QTLs for broomrape resistance were detected. One of the three QTLs explained more than 35% of the phenotypic variance, whereas the others accounted for 11.2 and 25.5%, respectively. This result suggests that broomrape resistance in faba bean can be considered a polygenic trait with major effects of a few single genes.Key words: Orobanche crenata, Vicia faba, QTL, broomrape resistance.