Heterozygosity and diversity analysis using mapped single nucelotide polymorphisms in a faba bean inbreeding programme

The aim of this study was to convert existing faba bean (Vicia faba L.) single nucleotide polymorphism (SNP) markers from cleaved amplification polymorphic sequences and SNaPshot® formats, which are expensive and time-consuming, to the more convenient KBiosciences competitive allele‐specific PCR (KASP) assay format. Out of 80 assays designed, 75 were validated, though a core set of 67 of the most robust markers is recommended for further use. The 67 best KASP SNP assays were used across two generations of single seed descent to detect unintended outcrossing and to track and quantify loss of heterozygosity, a capability that will significantly increase the efficiency and performance of pure line production and maintenance. This same set of assays was also used to examine genetic relationships between the 67 members of the partly inbred panel, and should prove useful for line identification and diversity studies in the future.

A baseline study of vicine–convicine levels in faba bean (Vicia faba L.) germplasm

Vicine and convicine are anti-nutritional compounds that accumulate in the cotyledons of faba beans. When humans consume beans with high levels of these compounds, it can cause a condition called favism in individuals harbouring a deficiency in the activity of their glucose-6-phosphate dehydrogenase. When faba beans are used in animal feeds, there can be effects on performance. These concerns have resulted in increasing interest within plant breeding in developing low vicine and convicine faba bean germplasm. In order to facilitate this objective, we developed a rapid and robust screening method for vicine and convicine, capable of distinguishing between faba beans that are either high (wild type) or low in vicine and convicine. In the absence of reliable commercial reference materials, we report an adaptation of a previously published method where a biochemical assay and spectral data were used to confirm the identity of our analytes, vicine and convicine. This method could be readily adopted in other facilities and open the way to the efficient exploitation of diverse germplasm in regions where faba beans play a significant role in human nutrition. We screened a collection of germplasm of interest to a collaborative plant breeding programme developing between the National Institute for Agricultural Botany in the UK and L'Institut Nationale d'Agronomie de Tunisie in Tunisia. We report the results obtained and discuss the prospects for developing molecular markers for the low vicine and convicine trait.

Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

Key message We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species. Abstract Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F5 recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene—a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.

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

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 collinear with Mt3. The obvious candidate gene from 77 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene (Mt3g092830, Mt3g092840) controlling anthocyanin production in Medicago and re-sequencing of the Vf orthologue showed a putative causative deletion of the entire 5’ end of the gene.

Flanking SNP markers for vicine–convicine concentration in faba bean (Vicia faba L).

The pyrimidine glycosides, vicine and convicine, limit the use of faba bean (Vicia faba L.) as food and feed. A single recessive gene, vc-, is responsible for a lowered vicine–convicine concentration. The biosynthetic pathway of these closely related compounds is not known, and the nearest available markers are several cM away from vc-. Improved markers would assist breeding and help to identify candidate genes. A segregating population of 210 F5 recombinant inbred lines was developed from the cross of Mélodie/2 (low vicine–convicine) × ILB 938/2 (normal vicine–convicine), and vicine–convicine concentrations were determined twice on each line. The population was genotyped with a set of 188 SNPs. A strong, single QTL for vicine–convicine concentration was identified on chromosome I, flanked by markers 1.0 cM away on one side and 2.6 cM on the other. The interval defined by these markers in the model species Medicago truncatula includes about 340 genes, but no candidate genes were identified. Further fine mapping should lead to the identification of tightly linked markers as well as narrowing down the search for candidate regulatory or biosynthetic genes which could underlie the vc- locus.

Insect pollination reduces yield loss following heat stress in faba bean (Vicia faba L.)

Global food security, particularly crop fertilization and yield production, is threatened by heat waves that are projected to increase in frequency and magnitude with climate change. Effects of heat stress on the fertilization of insect-pollinated plants are not well understood, but experiments conducted primarily in self-pollinated crops, such as wheat, show that transfer of fertile pollen may recover yield following stress. We hypothesized that in the partially pollinator-dependent crop, faba bean (Vicia faba L.), insect pollination would elicit similar yield recovery following heat stress. We exposed potted faba bean plants to heat stress for 5 days during floral development and anthesis. Temperature treatments were representative of heat waves projected in the UK for the period 2021-2050 and onwards. Following temperature treatments, plants were distributed in flight cages and either pollinated by domesticated Bombus terrestris colonies or received no insect pollination. Yield loss due to heat stress at 30°C was greater in plants excluded from pollinators (15%) compared to those with bumblebee pollination (2.5%). Thus, the pollinator dependency of faba bean yield was 16% at control temperatures (18 to 26°C) and extreme stress (34°C), but was 53% following intermediate heat stress at 30°C. These findings provide the first evidence that the pollinator dependency of crops can be modified by heat stress, and suggest that insect pollination may become more important in crop production as the probability of heat waves increases.

Susceptibility of faba bean (Vicia faba L.) to heat stress during floral development and anthesis

Experiments were conducted over two years to quantify the response of faba bean (Vicia faba L.) to heat stress. Potted winter faba bean plants (cv. Wizard) were exposed to temperature treatments (18/10; 22/14; 26/18; 30/22; 34/26°C day/night) for five days during floral development and anthesis. Developmental stages of all flowers were scored prior to stress, plants were grown in exclusion from insect pollinators to prevent pollen movement between flowers, and yield was harvested at an individual pod scale, enabling effects of heat stress to be investigated at a high resolution. Susceptibility to stress differed between floral stages, flowers were most affected during initial green-bud stages. Yield and pollen germination of flowers present before stress showed threshold relationships to stress, with lethal temperatures (t50) ~28°C and ~32°C, while whole plant yield showed a linear negative relationship to stress with high plasticity in yield allocation, such that yield lost at lower nodes was partially compensated at higher nodal positions. Faba bean has many beneficial attributes for sustainable modern cropping systems but these results suggest that yield will be limited by projected climate change, necessitating the development of heat tolerant cultivars, or improved resilience by other mechanisms such as earlier flowering times.

When are ant-attractant devices a worthwhile investment? Vicia faba extrafloral nectaries and Lasius niger ants

Most studies aiming to determine the beneficial effect of ants on plants simply consider the effects of the presence or exclusion of ants on plant yield. This approach is often inadequate, however, as ants interact with both non-tended herbivores and tended Homoptera. Moreover, the interaction with these groups of organisms is dependent on ant density, and these functional relationships are likely to be non-linear. A model is presented here that segregates plant herbivores into two categories depending on the sign of their numerical response to ants (myrmecophiles increase with ants, non-tended herbivores decline). The changes in these two components of herbivores with increasing ant density and the resulting implications for ant-plant mutualisms are considered. It emerges that a wide range of ant densities needs to be considered as the interaction sign (mutualism or parasitism) and strength is likely to change with ant density. The model is used to interpret the results of an experimental study that varied levels of Aphis fabae infestation and Lasius niger ant attendance on Vicia faba bean plants. Increasing ant density consistently reduced plant fitness and thus, in this location, the interaction between the ants and the plant can be considered parasitic. In the Vicia faba system, these costs of ants are unlikely to be offset by other beneficial agents (e.g., parasitoids), which also visit extrafloral nectaries.

Standardised ileal digestibility of crude protein and amino acids of UK-grown peas and faba beans by broilers

In view of the increasing interest in home-grown legumes as components of diets for non-ruminant livestock and in an attempt to reduce the reliance on imported soya bean meal (SBM), two experiments were conducted to evaluate samples of peas and faba beans for their standardised ileal digestibility (SID) of amino acids determined with young broiler chicks. Experiment 1 evaluated six faba bean and seven pea cultivars and Experiment 2 evaluated two faba bean and three pea cultivars as well as a sample of soya bean meal provided as a reference material. Peas and beans were added at 750g/kg as the only source of protein/amino acids in a semi-synthetic diet containing the inert marker titanium dioxide; SBM was added, in a control diet, at 500g/kg. Each diet was fed to six replicates of a cage containing two Ross-type broilers for 96h at which point birds were culled allowing removal of ileal digesta. Chemical analyses allowed the calculation of the coefficient of SID of amino acids. There were no differences between samples of the same pulse species (P>0.05) but peas had higher values (P<0.05), similar to SBM, than beans. Trypsin inhibitor content (expressed as g trypsin inhibitor units/mg sample) of all pea samples was low and in the range 0.83–1.77mg/kg. There was relatively little variation in bean tannin content and composition amongst the coloured-flowered varieties; however, the white-flowered cultivar had no tannins. There was no correlation between tannin content and coefficient of SID. The content of SID of amino acids (g/kg legume) was higher in SBM when compared with peas and beans by virtue of having higher total concentrations.

Genetic analysis reveals a novel locus in Vicia faba decoupling pigmentation in the flower from that in the extra-floral nectaries

The genetics of the stipule spot pigmentation (SSP) in faba bean (Vicia faba L.) was studied using four inbred lines, of which Disco/2 was zero-tannin (zt2) with colourless stipule spots, ILB938/2 was normal-tannin (ZT2) with colourless stipule spots, and both Aurora/2 and Mélodie/2 were ZT2 with coloured stipule spots. Crosses Mélodie/2 × ILB 938/2, Mélodie/2 × Disco/2, ILB 938/2 × Aurora/2 and ILB 938/2 × Disco/2 (A, B, C and D, respectively) were prepared, along with reciprocals and backcrosses, and advanced through single-seed descent. All F1 hybrid plants had pigmented stipule spots, and in the F2 generation, the segregation ratio fit 3 coloured:1 colourless in crosses A, B and C and 9:7 in cross D. In the F3 generation, the ratio fit 5:3 in crosses A and C and 25:39 in cross D, and in the F4 generation, 9:7 in cross A. SSP was linked to the zero-tannin characteristics (white flower) only in cross B. The results show that coloured stipule spot is dominant to colourless and that colouration is determined by two unlinked complementary recessive genes. We propose the symbols ssp2 for the gene associated with zt2 in Disco/2 and ssp1 for the gene not associated with tannin content in ILB938/2. The novel ssp1 locus was mapped at F5 in cross ‘A’ using Medicago truncatula-derived single-nucleotide polymorphism and was on chromosome 1 of faba bean, in a well-conserved region of M. truncatula chromosome 5 containing some candidate Myb and basic helix–loop–helix transcription factor genes.

Intercropping with pulses to concentrate nitrogen and sulphur in wheat

The effects of intercropping wheat with faba bean (Denmark, Germany, Italy and UK) and wheat with pea (France), in additive and replacement designs on grain nitrogen and sulphur concentrations were studied in field experiments in the 2002/03, 2003/04 and 2004/05 growing seasons. Intercropping wheat with grain legumes regularly increased the nitrogen concentration of the cereal grain, irrespective of design or location. Sulphur concentration of the cereal was also increased by intercropping, but less regularly and to a lesser extent compared with effects on nitrogen concentration. Nitrogen concentration (g/kg) in wheat additively intercropped with faba bean was increased by 8% across all sites (weighted for inverse of variance), but sulphur concentration was only increased by 4%, so N:S ratio was also increased by 4%. Intercropping wheat with grain legumes increased sodium dodecyl sulphate (SDS)-sedimentation volume. The effect of intercropping on wheat nitrogen concentration was greatest when intercropping had the most deleterious effect on wheat yield and the least deleterious effect on pulse yield. Over all sites and seasons, and irrespective of whether the design was additive or replacement, increases in crude protein concentration in the wheat of 10 g/kg by intercropping with faba bean were associated with 25-30% yield reduction of the wheat, compared with sole-cropped wheat. It was concluded that the increase in protein concentration of wheat grain in intercrops could be of economic benefit when selling wheat for breadmaking, but only if the bean crop was also marketed effectively.

Interactive effect of floral abundance and semi-natural habitats on pollinators in field beans (Vicia faba)

Pollination services are economically important component of agricultural biodiversity which enhance the yield and quality of many crops. An understanding of the suitability of extant habitats for pollinating species is crucial for planning management actions to protect and manage these service providers. In a highly modified agricultural ecosystem, we tested the effect of different pollination treatments (open, autonomous self- and wind-pollination) on pod set, seed set, and seed weight in field beans (Vicia faba). We also investigated the effect of semi-natural habitats and flower abundance on pollinators of field beans. Pollinator sampling was undertaken in ten field bean fields along a gradient of habitat complexity; CORINE land cover classification was used to analyse the land use patterns between 500–3000 m around the sites. Total yield from open-pollination increased by 185% compared to autonomous self-pollination. There was positive interactive effect of local flower abundance and cover of semi-natural habitats on overall abundance of pollinators at 1500 and 2000 m, and abundance of bumblebees (Bombus spp.) at 1000–2000 m. In contrast, species richness of pollinators was only correlated with flower abundance and not with semi-natural habitats. We did not find a link between pod set from open-pollination and pollinator abundance, possibly due to variations in the growing conditions and pollinator communities between sites. We conclude that insect pollination is essential for optimal bean yields and therefore the maintenance of semi-natural habitats in agriculture-dominated landscapes should ensure stable and more efficient pollination services in field beans.

Phylogeographic analysis of the chloroplast DNA variation in wild common bean (Phaseolus vulgaris L.) in the Americas

The wild common bean (Phaseolus vulgaris) is widely but discontinuously distributed from northern Mexico to northern Argentina on both sides of the Isthmus of Panama. Little is known on how the species has reached its current disjunct distribution. In this research, chloroplast DNA polymorphisms in seven non-coding regions were used to study the history of migration of wild P. vulgaris between Mesoamerica and South America. A penalized likelihood analysis was applied to previously published Leguminosae ITS data to estimate divergence times between P. vulgaris and its sister taxa from Mesoamerica, and divergence times of populations within P. vulgaris. Fourteen chloroplast haplotypes were identified by PCR-RFLP and their geographical associations were studied by means of a Nested Clade Analysis and Mantel Tests. The results suggest that the haplotypes are not randomly distributed but occupy discrete parts of the geographic range of the species. The current distribution of haplotypes may be explained by isolation by distance and by at least two migration events between Mesoamerica and South America: one from Mesoamerica to South America and another one from northern South America to Mesoamerica. Age estimates place the divergence of P. vulgaris from its sister taxa from Mesoamerica at or before 1.3 Ma, and divergence of populations from Ecuador-northern Peru at or before 0.6 Ma. As these ages are taken as minimum divergence times, the influence of past events, such as the closure of the Isthmus of Panama and the final uplift of the Andes, on the migration history and population structure of this species cannot be disregarded.