A comparative metabolomics investigation of flavonoid variation in faba bean flowers.

INTRODUCTION: Faba bean (Vicia faba L.) flowers are edible and used as garnishes because of their aroma, sweet flavor and attractive colors. Anthocyanins are the common plant pigments that give flowers their vivid colors, whereas non-anthocyanin flavonoids can serve as co-pigments that can modify the color intensity of flowers. OBJECTIVES: To explore the polyphenol diversity and differences in standard and wing petals of faba bean flowers; and identify glycosylated flavonoids that contribute to flower color. METHODS: Flower standard and wing petals from 30 faba bean genotypes (eight color groups with a total of 60 samples) were used for polyphenol extraction. Samples were analyzed using a targeted method and a semi-untargeted analysis using liquid chromatography-high resolution mass spectrometry (LC-HRMS) combined with photodiode array (PDA) detection. Compound Discoverer software was used for polyphenol identification and multivariate analysis. RESULTS: The semi-untargeted analysis guided by the PDA detected 90 flavonoid metabolites present in faba bean flower petals. Ten anthocyanins largely influenced the flower colors, but other flavonoids (63 flavonols and 12 flavones) found with variable levels in different flower color groups appeared to also influence color, especially in mixed colors. CONCLUSION: Analysis of the different colored faba bean flowers confirmed that the color variation between the flowers was mainly controlled by anthocyanins in brown, red and purple-red flowers. Of the other flavonoids, multiglycosylated kaempferols were abundant in white and brown flowers, monoglycosylated kaempferols were common in red and purple-red flowers, and quercetin and apigenin glycosides were abundant co-pigments in purple-red flowers.

Genetic analysis of global faba bean diversity, agronomic traits and selection signatures.

KEY MESSAGE: We identified marker-trait associations for key faba bean agronomic traits and genomic signatures of selection within a global germplasm collection. Faba bean (Vicia faba L.) is a high-protein grain legume crop with great potential for sustainable protein production. However, little is known about the genetics underlying trait diversity. In this study, we used 21,345 high-quality SNP markers to genetically characterize 2678 faba bean genotypes. We performed genome-wide association studies of key agronomic traits using a seven-parent-MAGIC population and detected 238 significant marker-trait associations linked to 12 traits of agronomic importance. Sixty-five of these were stable across multiple environments. Using a non-redundant diversity panel of 685 accessions from 52 countries, we identified three subpopulations differentiated by geographical origin and 33 genomic regions subjected to strong diversifying selection between subpopulations. We found that SNP markers associated with the differentiation of northern and southern accessions explained a significant proportion of agronomic trait variance in the seven-parent-MAGIC population, suggesting that some of these traits were targets of selection during breeding. Our findings point to genomic regions associated with important agronomic traits and selection, facilitating faba bean genomics-based breeding.

The giant diploid faba genome unlocks variation in a global protein crop.

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) 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.

Identifying physiological and genetic determinants of faba bean transpiration response to evaporative demand.

BACKGROUND AND AIMS: Limiting maximum transpiration rate (TR) under high vapour pressure deficit (VPD) works as a water conservation strategy. While some breeding programmes have incorporated this trait into some crops to boost yields in water-limited environments, its underlying physiological mechanisms and genetic regulation remain unknown for faba bean (Vicia faba). Thus, we aimed to identify genetic variation in the TR response to VPD in a population of faba bean recombinant inbred lines (RILs) derived from two parental lines with contrasting water use (Mélodie/2 and ILB 938/2). METHODS: Plants were grown in well-watered soil in a climate-controlled glasshouse with diurnally fluctuating VPD and light conditions. Whole plant transpiration was measured in a gas exchange chamber that tightly regulated VPD around the shoot under constant light, while whole-plant hydraulic conductance and its components (root and stem hydraulic conductance) were calculated from dividing TR by water potential gradients measured with a pressure chamber. KEY RESULTS: Although TR of Mélodie/2 increased linearly with VPD, ILB 938/2 limited its TR above 2.0 kPa. Nevertheless, Mélodie/2 had a higher leaf water potential than ILB 938/2 at both low (1.0 kPa) and high (3.2 kPa) VPD. Almost 90 % of the RILs limited their TR at high VPD with a break-point (BP) range of 1.5-3.0 kPa and about 10 % had a linear TR response to VPD. Thirteen genomic regions contributing to minimum and maximum transpiration, and whole-plant and root hydraulic conductance, were identified on chromosomes 1 and 3, while one locus associated with BP transpiration was identified on chromosome 5. CONCLUSIONS: This study provides insight into the physiological and genetic control of transpiration in faba bean and opportunities for marker-assisted selection to improve its performance in water-limited environments.

Advances in disease and pest resistance in faba bean.

Faba bean (Vicia faba) is a grain legume crop widely cultivated in temperate areas for food and feed. Its productivity can be constrained by numerous diseases and pests that can be managed by a number of strategies, complemented with the deployment of resistant cultivars in an integrated manner. Few sources of resistance are available to some of them, although their phenotypic expression is usually insufficiently described, and their genetic basis is largely unknown. A few DNA markers have been developed for resistance to rust, ascochyta blight, and broomrape, but not yet for other diseases or pests. Still, germplasm screenings are allowing the identification of resistances that are being accumulated by classical breeding, succeeding in the development of cultivars with moderate levels of resistance. The adoption of novel phenotyping approaches and the unprecedented development of genomic resources along with speed breeding tools are speeding up resistance characterization and effective use in faba bean breeding.

Genomic regions associated with chocolate spot (Botrytis fabae Sard.) resistance in faba bean (Vicia faba L.).

Chocolate spot (CS), caused by Botrytis fabae Sard., is an important threat to global faba bean production. Growing resistant faba bean cultivars is, therefore, paramount to preventing yield loss. To date, there have been no reported quantitative trait loci (QTL) associated with CS resistance in faba bean. The objective of this study was to identify genomic regions associated with CS resistance using a recombinant inbred line (RIL) population derived from resistant accession ILB 938. A total of 165 RILs from the cross Mélodie/2 × ILB 938/2 were genotyped and evaluated for CS reactions under replicated controlled climate conditions. The RIL population showed significant variation in response to CS resistance. QTL analysis identified five loci contributing to CS resistance on faba bean chromosomes 1 and 6, accounting for 28.4% and 12.5%, respectively, of the total phenotypic variance. The results of this study not only provide insight into disease-resistance QTL, but also can be used as potential targets for marker-assisted breeding in faba bean genetic improvement for CS resistance. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01307-7.

Conventional and Molecular Breeding Tools for Accelerating Genetic Gain in Faba Bean (Vicia Faba L.).

Faba bean is a cool-season grain legume crop, which is grown worldwide for food and feed. Despite a decrease in area under faba bean in the past, the interest in growing faba bean is increasing globally due to its high seed protein content and its excellent ecological service. The crop is, however, exposed to diverse biotic and abiotic stresses causing unstable, low grain yield. Although, sources of resistance to main diseases, such as ascochyta blight (Ascochyta fabae Speg.), rust (Uromyces viciae-fabae (Pers.) Schroet.), chocolate spot (Botrytis fabae Sard.) and gall disease (Physioderma viciae), have been identified, their resistance is only partial and cannot prevent grain yield losses without agronomical practices. Tightly associated DNA markers for host plant resistance genes are needed to enhance the level of resistance. Less progress has been made for abiotic stresses. Different breeding methods are proposed, but until now line breeding, based on the pedigree method, is the dominant practice in breeding programs. Nonetheless, the low seed multiplication coefficient and the requirement for growing under insect-proof enclosures to avoid outcrossing hampers breeding, along with the lack of tools such as double haploid system and cytoplasmic male sterility. This reduces breeding population size and speed of breeding hence the chances of capturing rare combinations of favorable alleles. Availability and use of the DNA markers such as vicine-convicine (vc -) and herbicide tolerance in breeding programs have encouraged breeders and given confidence in marker assisted selection. Closely linked QTL for several biotic and abiotic stress tolerance are available and their verification and conversion in breeder friendly platform will enhance the selection process. Recently, genomic selection and speed breeding techniques together with genomics have come within reach to accelerate the genetic gains in faba bean. Advancements in genomic resources with other breeding tools, methods and platforms will enable to accelerate the breeding process for enhancing genetic gain in this species.

VC1 catalyses a key step in the biosynthesis of vicine in faba bean.

Faba bean (Vicia faba L.) is a widely adapted and high-yielding legume cultivated for its protein-rich seeds1. However, the seeds accumulate the pyrimidine glucosides vicine and convicine, which can cause haemolytic anaemia (favism) in 400 million genetically predisposed individuals2. Here, we use gene-to-metabolite correlations, gene mapping and genetic complementation to identify VC1 as a key enzyme in vicine and convicine biosynthesis. We demonstrate that VC1 has GTP cyclohydrolase II activity and that the purine GTP is a precursor of both vicine and convicine. Finally, we show that cultivars with low vicine and convicine levels carry an inactivating insertion in the coding sequence of VC1. Our results reveal an unexpected, purine rather than pyrimidine, biosynthetic origin for vicine and convicine and pave the way for the development of faba bean cultivars that are free of these anti-nutrients.

Recent advances in faba bean genetic and genomic tools for crop improvement.

Faba bean (Vicia faba L.), a member of the Fabaceae family, is one of the important food legumes cultivated in cool temperate regions. It holds great importance for human consumption and livestock feed because of its high protein content, dietary fibre, and nutritional value. Major faba bean breeding challenges include its mixed breeding system, unknown wild progenitor, and genome size of ~13 Gb, which is the largest among diploid field crops. The key breeding objectives in faba bean include improved resistance to biotic and abiotic stress and enhanced seed quality traits. Regarding quality traits, major progress on reduction of vicine-convicine and seed coat tannins, the main anti-nutritional factors limiting faba bean seed usage, have been recently achieved through gene discovery. Genomic resources are relatively less advanced compared with other grain legume species, but significant improvements are underway due to a recent increase in research activities. A number of bi-parental populations have been constructed and mapped for targeted traits in the last decade. Faba bean now benefits from saturated synteny-based genetic maps, along with next-generation sequencing and high-throughput genotyping technologies that are paving the way for marker-assisted selection. Developing a reference genome, and ultimately a pan-genome, will provide a foundational resource for molecular breeding. In this review, we cover the recent development and deployment of genomic tools for faba bean breeding.

Do Faba Bean Genotypes Carrying Different Zero-Tannin Genes (zt1 and zt2) Differ in Phenolic Profiles?

Faba bean is a cool season grain legume that produces seeds with a high protein content. Seed coat tannins limit its use in food and feed. A low-tannin phenotype is controlled by either of two unlinked recessive genes zt1 and zt2. Liquid chromatography-mass spectrometry was used to characterize phenolic profiles of seed coat and flower tissue of three faba bean genotypes: CDC Snowdrop (zt1 gene), Disco/2 (zt2 gene), and ILB 938/2 (tannin-containing). For both tissues, clear differences in phenolic profiles of ILB 938/2 were observed in comparison to both low-tannin lines. Although seed coat phenolic profiles of zt1 and zt2 genotypes were similar, distinct differences were evident in flower tissue, suggesting that the gene action results in some different end products of the phenolic biosynthetic pathway. These distinctive compounds could be used as biochemical markers to distinguish between low-tannin phenotypes.

Seed Protein of Lentils: Current Status, Progress, and Food Applications.

Grain legumes are widely recognized as staple sources of dietary protein worldwide. Lentil seeds are an excellent source of plant-based proteins and represent a viable alternative to animal and soybean proteins for food processing formulations. Lentil proteins provide not only dietary amino acids but are also a source of bioactive peptides that provide health benefits. This review focuses on the current knowledge of seed protein, extraction and isolation methods, bioactive peptides, and food applications of lentil protein. Lentil is the most rapidly expanding crop for direct human consumption, and has potential for greater impact as a protein source for food processing applications. Improvements in lentil protein quality, amino acid composition, and processing fractions will enhance the nutritional quality of this rapidly expanding crop globally.

Genetic Mapping of Milling Quality Traits in Lentil (Lens culinaris Medik.).

Milling qualities are key traits for the red lentil ( Medik.) industry as price is largely determined by recovery yield. milling involves removal of the seed coat and splitting of the cotyledon to produce either splits or footballs (cotyledons still attached). The objectives of the study were to determine the heritability of the milling traits dehulling efficiency (DE), milling recovery (MR), and football recovery (FR) and to identify the genomic regions controlling them. We used a lentil recombinant inbred population from the cross 'CDC Robin' × '946a-46', which have contrasting seed characteristics. The mapping population consists of 127 F-derived lentil recombinant inbred lines that were phenotyped for milling quality parameters from four site-years in Saskatchewan, Canada. A total of 534 single nucleotide polymorphism markers, seven simple sequence repeat markers, and four morphological markers were used for quantitative trait locus (QTL) mapping. The broad-sense heritability was moderate for DE and MR and relatively low for FR. Milling quality traits were significantly correlated with seed shape (seed diameter and seed plumpness). Multiple QTLs for milling traits were detected in six of seven linkage groups (LGs). The most stable QTLs governing DE and MR were clustered on LGs 1, 2, 3, and 7, whereas FR QTLs were clustered on LGs 4, 5, 6, and 7. The molecular markers identified for these traits could be used for improving milling quality in lentil breeding programs.

Toward a high-throughput method for determining vicine and convicine levels in faba bean seeds using flow injection analysis combined with tandem mass spectrometry.

Although faba bean provides environmental and health benefits, vicine and convicine (v-c) limit its use as a source of vegetable protein. Crop improvement efforts to minimize v-c concentration require low-cost, rapid screening methods to distinguish between high and low v-c genotypes to accelerate development of new cultivars and to detect out-crossing events. To assist crop breeders, we developed a unique and rapid screening method that uses a 60 s instrumental analysis step to accurately distinguish between high and low v-c genotypes. The method involves flow injection analysis (FIA) coupled with tandem mass spectrometry (i.e., selective reaction monitoring, SRM). Using seeds with known v-c levels as calibrants, measured v-c levels were comparable with liquid chromatography (LC)-SRM results and the method was used to screen 370 faba bean genotypes. Widespread use of FIA-SRM will accelerate breeding of low v-c faba bean, thereby alleviating concerns about anti-nutritional effects of v-c in this crop.

Single Nucleotide Polymorphism Markers Associated with Seed Quality Characteristics of Cultivated Lentil.

The dimensions of lentil ( Medik.) seeds are important quality parameters that are major determinants of market preference, cooking time, and post-harvest milling quality. Knowledge of the genetic control of traits related to seed dimensions would be useful for crop improvement. The principal aim of this study was to identify single nucleotide polymorphism (SNP) markers linked to genes that control seed diameter, seed thickness, and seed plumpness. Association mapping analysis with SNP markers was used to study the seed dimensions of 138 diverse cultivated lentil accessions grown at two locations in Saskatchewan, Canada, in 2011 and 2012. Six marker-trait associations were shown to be significant for the studied seed dimension characteristics. Two SNP markers closely associated with seed diameter across locations and years identified in previous work were validated in this study. Three additional marker-seed thickness associations were identified. Using the association mapping strategy, we confirmed the presence of two genomic regions controlling seed diameter and plumpness. This information can be used worldwide as a resource for lentil seed quality improvement programs.

Quantification of vicine and convicine in faba bean seeds using hydrophilic interaction liquid chromatography.

Faba bean (Vicia faba L.) provides environmental and health benefits; however, the presence of the pyrimidine glycosides vicine and convicine (v-c) in its seeds limits consumption. Low v-c genotypes have been introduced, but the convicine levels in these genotypes have not been quantified. To improve detection, the polar nature of v-c was exploited by implementing hydrophilic interaction liquid chromatography (HILIC). A sample preparation method using a two-step extraction was developed for use with UV and/or tandem mass spectrometry (SRM) detection. The HILIC-UV method was suitable for over three orders of magnitude, covering the range of v-c concentrations in faba bean seeds across all genotypes tested. The linear range of HILIC-SRM was slightly less (∼3 orders of magnitude), but improved sensitivity and selectivity make it more suitable for quantifying low v-c samples. The analysis of 13 genotypes suggests that v-c concentrations in faba bean seeds may be independent quantitative traits.

Bioavailability of coated and uncoated ZnO nanoparticles to cucumber in soil with or without organic matter.

There is a gap of knowledge for the fate, effects and bioavailability of coated and uncoated ZnO nanoparticles (NPs) in soil. Moreover, little is known about the effects of soil properties on effects of NPs on plants. In this study, the availability ZnO NPs in two soils with different organic matter content (one treated with cow manure (CM) and the other as untreated) was compared with their bulk particles. Results showed that coated and uncoated ZnO NPs can be more bioaccessible than their bulk counterpart and despite their more positive effects at low concentration (< 100mgkg-1), they were more phytotoxic for plants compared to the bulk ZnO particles at high concentration (1000mgkg-1) in the soil untreated with CM. The concentration of 1000mgkg-1 of ZnO NPs, decreased shoot dry biomass (52%) in the soil untreated with CM but increased shoot dry biomass (35%) in CM-treated soil compared to their bulk counterpart. In general, plants in the CM-treated soil showed higher Zn concentration in their tissues compared with those in untreated soil. The difference in shoot Zn concentration between CM-treated and untreated soil for NPs treatments was more than bulk particles treatment. This different percentage at 100mgkg-1 of bulk particles was 20.6% and for coated and uncoated NPs were 37% and 32%, respectively. Generally, the distribution of ZnO among Zn fractions in soil (exchangeable, the metal bound to carbonates, Fe-Mn oxides, organic matter and silicate minerals and the residual fraction) changed based on applied Zn concentration, Zn source and soil organic matter content. The root tip deformation under high concentration of NPs (1000mgkg-1 treatment) was observed by light microscopy in plants at the soil untreated with CM. It seems that root tip deformation is one of the specific effects of NPs which in turn inhibits plant growth and nutrients uptake by root. The transmission electron microcopy image showed the aggregation of NPs inside the plant cytoplasm and their accumulation adjacent to the cell membrane.

Marker-Trait Association Analysis of Iron and Zinc Concentration in Lentil (Lens culinaris Medik.) Seeds.

Lentil ( Medik.) seeds are relatively rich in iron (Fe) and zinc (Zn), making lentil a potential crop to aid in the global battle against human micronutrient deficiency. Understanding the genetic basis for uptake of seed Fe and Zn is required to increase sustainable concentrations of these minerals in seeds. The objectives of this study were to characterize genetic variation in seed Fe and Zn concentration and to identify molecular markers associated with these traits across diverse lentil accessions. A set of 138 cultivated lentil accessions from 34 countries were evaluated in four environments (2 sites × 2 yr) in Saskatchewan, Canada. The collection was genotyped using 1150 single-nucleotide polymorphism (SNP) markers that are distributed across the lentil genome. The germplasm tested exhibited a wide range of variation for seed Fe and Zn concentration. The marker-trait association analysis detected two SNP markers tightly linked to seed Fe and one linked to seed Zn concentration (-log10 ≥ 4.36). Additional markers were detected at -log10 ≥ 3.06. A number of putative candidate genes underlying detected loci encode Fe- and Zn-related functions. This study provides insight into the genetics of seed Fe and Zn concentration of lentil and opportunities for marker-assisted selection to improve micronutrient concentration as part of micronutrient biofortification programs.

Genetic Diversity of Cultivated Lentil (Lens culinaris Medik.) and Its Relation to the World's Agro-ecological Zones.

Assessment of genetic diversity and population structure of germplasm collections plays a critical role in supporting conservation and crop genetic enhancement strategies. We used a cultivated lentil (Lens culinaris Medik.) collection consisting of 352 accessions originating from 54 diverse countries to estimate genetic diversity and genetic structure using 1194 polymorphic single nucleotide polymorphism (SNP) markers which span the lentil genome. Using principal coordinate analysis, population structure analysis and UPGMA cluster analysis, the accessions were categorized into three major groups that prominently reflected geographical origin (world's agro-ecological zones). The three clusters complemented the origins, pedigrees, and breeding histories of the germplasm. The three groups were (a) South Asia (sub-tropical savannah), (b) Mediterranean, and (c) northern temperate. Based on the results from this study, it is also clear that breeding programs still have considerable genetic diversity to mine within the cultivated lentil, as surveyed South Asian and Canadian germplasm revealed narrow genetic diversity.

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 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.

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. 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.