Development of new genetic resources for faba bean (Vicia faba L.) breeding through the discovery of gene-based SNP markers and the construction of a high-density consensus map

Transcriptome Analysis of Stigmas of Vicia faba L. Flowers

Pollination in angiosperms depends on complex communication between pollen grains and stigmas, classified as wet or dry, depending on the presence or absence of secretions at the stigma surface, respectively. In species with wet stigma, the cuticle is disrupted and the presence of exudates is indicative of their receptivity. Most stigma studies are focused on a few species and families, many of them with self-incompatibility systems. However, there is scarce knowledge about the stigma composition in Fabaceae, the third angiosperm family, whose stigmas have been classified as semidry. Here we report the first transcriptome profiling and DEGs of Vicia faba L. styles and stigmas from autofertile (flowers able to self-fertilize in the absence of manipulation, whose exudate is released spontaneously) and autosterile (flowers that need to be manipulated to break the cuticle and release the exudates to be receptive) inbred lines. From the 76,269 contigs obtained from the de novo assembly, only 45.1% of the sequences were annotated with at least one GO term. A total of 115,920, 75,489, and 70,801 annotations were assigned to Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) categories, respectively, and 5918 differentially expressed genes (DEGs) were identified between the autofertile and the autosterile lines. Among the most enriched metabolic pathways in the DEGs subset were those related with amino acid biosynthesis, terpenoid metabolism, or signal transduction. Some DEGs have been related with previous QTLs identified for autofertility traits, and their putative functions are discussed. The results derived from this work provide an important transcriptomic reference for style-stigma processes to aid our understanding of the molecular mechanisms involved in faba bean fertilization.

High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-seq data

The genome of faba bean was first published in 2023. To promote future molecular breeding studies, we improved the quality of the faba genome based on high-density genetic maps and the Illumina and Pacbio RNA-seq datasets. Two high-density genetic maps were used to conduct the scaffold ordering and orientation of faba bean, culminating in an increased length (i.e., 14.28 Mbp) of chromosomes and a decrease in the number of scaffolds by 45. In gene model mining and optimisation, the PacBio and Illumina RNA-seq datasets from 37 samples allowed for the identification and correction 121,606 transcripts, and the data facilitated a prediction of 15,640 alternative splicing events, 2,148 lncRNAs, and 1,752 fusion transcripts, thus allowing for a clearer understanding of the gene structures underlying the faba genome. Moreover, a total of 38,850 new genes including 56,188 transcripts were identified compared with the reference genome. Finally, the genetic data of the reference genome was integrated and a comprehensive and complete faba bean transcriptome sequence of 103,267 transcripts derived from 54,753 uni-genes was formed.

Empirical phenotyping and genome-wide association study reveal the association of panicle architecture with yield in Chenopodium quinoa

Chenopodium quinoa manifests adaptability to grow under varying agro-climatic scenarios. Assessing quinoa germplasm's phenotypic and genetic variability is a prerequisite for introducing it as a potential candidate in cropping systems. Adaptability is the basic outcome of ecological genomics of crop plants. Adaptive variation predicted with a genome-wide association study provides a valuable basis for marker-assisted breeding. Hence, a panel of 72 quinoa plants was phenotyped for agro morphological attributes and association-mapping for distinct imperative agronomic traits. Inter simple sequence repeat (ISSR) markers were employed to assess genetic relatedness and population structure. Heatmap analysis showed three genotypes were early maturing, and six genotypes were attributed for highest yield. The SD-121-07 exhibited highest yield per plant possessing green, glomerulate shaped, compact density panicle with less leaves. However, SJrecm-03 yielded less exhibiting pink, intermediate shape, intermediate density panicles with less leaves. The phenotyping revealed strong correlation of panicle architecture with yield in quinoa. A genome-wide association study unraveled the associations between ISSR makers and agro-morphological traits. Mixed linear modes analysis yielded nine markers associated with eight traits at p ≤ 0.01. Moreover, ISSR markers significantly associated with panicle shape and leafiness were also associated with yield per plant. These findings contribute to the provision of authenticity for marker-assisted selection that ultimately would support quinoa breeding programs.

Novel SNP markers for flowering and seed quality traits in faba bean (Vicia faba L.): characterization and GWAS of a diversity panel

Faba bean (Vicia faba L.) is a legume crop grown in diverse climates worldwide. It has a high potential for increased cultivation to meet the need for more plant-based proteins in human diets, a prerequisite for a more sustainable food production system. Characterization of diversity panels of crops can identify variation in and genetic markers for target traits of interest for plant breeding. In this work, we collected a diversity panel of 220 accessions of faba bean from around the world consisting of gene bank material and commercially available cultivars. The aims of this study were to quantify the phenotypic diversity in target traits to analyze the impact of breeding on these traits, and to identify genetic markers associated with traits through a genome-wide association study (GWAS). Characterization under field conditions at Nordic latitude across two years revealed a large genotypic variation and high broad-sense heritability for eleven agronomic and seed quality traits. Pairwise correlations showed that seed yield was positively correlated to plant height, number of seeds per plant, and days to maturity. Further, susceptibility to bean weevil damage was significantly higher for early flowering accessions and accessions with larger seeds. In this study, no yield penalty was found for higher seed protein content, but protein content was negatively correlated to starch content. Our results showed that while breeding advances in faba bean germplasm have resulted in increased yields and number of seeds per plant, they have also led to a selection pressure towards delayed onset of flowering and maturity. DArTseq genotyping identified 6,606 single nucleotide polymorphisms (SNPs) by alignment to the faba bean reference genome. These SNPs were used in a GWAS, revealing 51 novel SNP markers significantly associated with ten of the assessed traits. Three markers for days to flowering were found in predicted genes encoding proteins for which homologs in other plant species regulate flowering. Altogether, this work enriches the growing pool of phenotypic and genotypic data on faba bean as a valuable resource for developing efficient breeding strategies to expand crop cultivation.

Construction of a high-density genetic map for faba bean (Vicia faba L.) and quantitative trait loci mapping of seed-related traits

Faba bean (Vicia faba L.) is a valuable legume crop and data on its seed-related traits is required for yield and quality improvements. However, basic research on faba bean is lagging compared to that of other major crops. In this study, an F₂ faba bean population, including 121 plants derived from the cross WY7×TCX7, was genotyped using the Faba_bean_130 K targeted next-generation sequencing genotyping platform. The data were used to construct the first ultra-dense faba bean genetic map consisting of 12,023 single nucleotide polymorphisms markers covering 1,182.65 cM with an average distance of 0.098 cM. The map consisted of 6 linkage groups, which is consistent with the 6 faba bean chromosome pairs. A total of 65 quantitative trait loci (QTL) for seed-related traits were identified (3 for 100-seed weight, 28 for seed shape, 12 for seed coat color, and 22 for nutritional quality). Furthermore, 333 candidate genes that are likely to participate in the regulation of seed-related traits were also identified. Our research findings can provide a basis for future faba bean marker-assisted breeding and be helpful to further modify and improve the reference genome.

Genetic Analyses of Mungbean [Vigna radiata (L.) Wilczek] Breeding Traits for Selecting Superior Genotype(s) Using Multivariate and Multi-Traits Indexing Approaches

Mungbean [Vigna radiata (L.) Wilczek] is an important food, feed, and cash crop in rice-based agricultural ecosystems in Southeast Asia and other continents. It has the potential to enhance livelihoods due to its palatability, nutritional content, and digestibility. We evaluated 166 diverse mungbean genotypes in two seasons using multivariate and multi-traits index approaches to identify superior genotypes. The total Shannon diversity index (SDI) for qualitative traits ranged from moderate for terminal leaflet shape (0.592) to high for seed colour (1.279). The analysis of variances (ANOVA) indicated a highly significant difference across the genotypes for most of the studied traits. Descriptive analyses showed high diversity among genotypes for all morphological traits. Six components with eigen values larger than one contributed 76.50% of the variability in the principal component analysis (PCA). The first three PCs accounted for the maximum 29.90%, 15.70%, and 11.20% of the total variances, respectively. Yield per plant, pod weight, hundred seed weight, pod length, days to maturity, pods per plant, harvest index, biological yield per plant, and pod per cluster contributed more to PC1 and PC2 and showed a positive association and positive direct effect on seed yield. The genotypes were grouped into seven clusters with the maximum in cluster II (34) and the minimum in cluster VII (10) along with a range of intra-cluster and inter-cluster distances of 5.15 (cluster II) to 3.60 (cluster VII) and 9.53 (between clusters II and VI) to 4.88 (clusters I and VII), suggesting extreme divergence and the possibility for use in hybridization and selection. Cluster III showed the highest yield and yield-related traits. Yield per plant positively and significantly correlated with pod traits and hundred seed weight. Depending on the multi-trait stability index (MTSI), clusters I, III, and VII might be utilized as parents in the hybridization program to generate high-yielding, disease-resistant, and small-seeded mungbean. Based on all multivariate-approaches, G45, G5, G22, G55, G143, G144, G87, G138, G110, G133, and G120 may be considered as the best parents for further breeding programs.

Homolog of Pea SGR Controls Stay-Green in Faba Bean (Vicia faba L.)

Faba bean is an important legume crop consumed as a vegetable or snack food, and its green cotyledons could present an attractive color for consumers. A mutation in SGR causes stay-green in plants. In this study, vfsgr was identified from a green-cotyledon-mutant faba bean, SNB7, by homologous blast between the SGR of pea and the transcriptome of faba bean. Sequence analysis revealed that a SNP at position 513 of the CDS of VfSGR caused a pre-stop codon, resulting in a shorter protein in the green-cotyledon faba bean SNB7. A dCaps marker was developed according to the SNP that caused the pre-stop, and this marker was completely associated with the color of the cotyledon of faba bean. SNB7 stayed green during dark treatment, while the expression level of VfSGR increased during dark-induced senescence in the yellow-cotyledon faba bean HST. Transient expression of VfSGR in Nicotiana. benthamiana leaves resulted in chlorophyll degradation. These results indicate that vfsgr is the gene responsible for the stay-green of faba bean, and the dCaps marker developed in this study provides a molecular tool for the breeding of green-cotyledon faba beans.

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.

Genome Size Variation across a Cypriot Fabeae Tribe Germplasm Collection

DNA content is an important trait linked to the evolutionary routes of taxa and often connected to speciation. In the present study, we studied C-values variation across the Cypriot Fabeae gene pool. Several hundred plants (Vicia spp., Lens spp., Pisum spp.) were sampled across Cyprus. Accurate estimates were established by flow cytometry and propidium iodine staining for 155 discrete populations/accessions. A ten-fold variation was detected across lineages with 1C DNA content varying from 1.584 pg for V. cretica (ARI02420) to 13.983 pg for V. faba (ARI00187). In general, flow cytometry was precise for the characterization of species, even though there were instances of genome overlapping across taxa. Most analyses in the current work refer to species that have not been characterized before by flow cytometry (or any other DNA content estimation method). Still, a correlation to C-values previously reported in Kew Plant DNA C-values database was attempted. A high degree of correlation except for V. dalmatica was established. The evaluation of genome size trait in relation with the Fabeae phylogeny, revealed that Pisum and Lens genera were rather homogenous, but an astonishing fluctuation was shown for Vicia spp. Moreover, it was established that genome up- or down-scaling was not directly linked to speciation drivers. The genomic size measurements presented here could deliver extra quality control for the identification and characterization of taxa in germplasm collections, particularly in cases where species share morphological characters.

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.

Genetic and Metabolite Variability among Commercial Varieties and Advanced Lines of Vicia faba L

Vicia faba L. (faba bean) is one of the most promising pulse crops due to its nutritional value and high nitrogen fixation capacity. The aim of the present study was to compare the genetic diversity and the seed metabolite profiles of five genetic materials of faba bean. Specifically, three newly developed advanced lines (KK18, KK14 and KK10) and two commercial cultivars (POLIKARPI and TANAGRA), were evaluated for this purpose. Genetic diversity among populations was assessed by SCoT molecular markers. Through UPGMA dendrogram, genetic distances between populations were estimated. Untargeted metabolomics analysis of the seeds was performed employing GC/EI/MS. The cultivar POLYKARPI exhibited the highest polymorphism. All varieties showed a higher within-cultivars and advanced lines variability than between. POLYKARPI and KK14 had the lowest genetic distances, while KK18 and TANAGRA presented the highest ones. The advanced line KK18 displayed the best nutritional profile, the highest concentration of desirable metabolites (lactic acid and trehalose), the lowest concentration of anti-nutritional factors (oxalic acid) and the lowest concentration of saturated fatty acids (palmitic and stearic acid). According to the results of the present study, KK18 line is a very promising material for further exploration and utilization in breeding programs.

RNA-Seq and genetic diversity analysis of faba bean (Vicia faba L.) varieties in China

Background

Faba bean (Vicia faba L) is one of the most important legumes in the world. However, there is relatively little genomic information available for this species owing to its large genome. The lack of data impedes the discovery of molecular markers and subsequent genetic research in faba bean. The objective of this study was to analyze the faba bean transcriptome, and to develop simple sequence repeat (SSR) markers to determine the genetic diversity of 226 faba bean varieties derived from different regions in China.

Methods

Faba bean varieties with different phenotype were used in transcriptome analysis. The functions of the unigenes were analyzed using various database. SSR markers were developed and the polymorphic markers were selected to conduct genetic diversity analysis.

Results

A total of 92.43 Gb of sequencing data was obtained in this study, and 133,487 unigene sequences with a total length of 178,152,541 bp were assembled. A total of 5,200 SSR markers were developed on the basis of RNA-Seq analysis. Then, 200 SSR markers were used to evaluate polymorphisms. In total, 103 (51.5%) SSR markers showed significant and repeatable bands between different faba bean varieties. Clustering analysis revealed that 226 faba bean materials were divided into five groups. Genetic diversity analysis revealed that the relationship between different faba beans in China was related, especially in the same region. These results provided a valuable data resource for annotating genes to different categories and developing SSR markers.

Identification of Rice Accessions Having Cold Tolerance at the Seedling Stage and Development of Novel Genotypic Assays for Predicting Cold Tolerance

Rice is susceptible to cold stress at the seedling stage, which can delay growth and decrease yield. We evaluated 187 rice accessions for cold tolerance at the seedling stage and developed genotypic assays for three markers. All japonica (20/20) and 20/140 indica accessions were highly cold tolerant. Two SNP markers specific for COLD1 and LOC_Os10g34840 were practical to use by normal agarose gel. The SNP marker specific for COLD1 was highly specific for predicting cold tolerance. However, the sensitivity of this marker was low as several cold-tolerant indica accessions lacked the cold-tolerant allele. The LOC_Os10g34840 marker was slightly more sensitive than the COLD1 marker for predicting highly cold-tolerant accessions. An insertion/deletion variant in the NAC6 gene was identified as a novel cold tolerance marker. The NAC6 marker predicted more highly cold-tolerant accessions compared with the other two markers. The SNP marker specific for LOC_Os10g34840 and the NAC6 marker were present in several tested subgroups, suggesting their wide effects and distribution. The three markers combined predicted the most highly cold-tolerant accessions, indicating that the marker combination is superior for applications such as marker-assisted breeding. The cold-tolerant accessions and the genotypic marker assays will be useful for future rice breeding.

SeSAM: software for automatic construction of order-robust linkage maps

Background

Genotyping and sequencing technologies produce increasingly large numbers of genetic markers with potentially high rates of missing or erroneous data. Therefore, the construction of linkage maps is more and more complex. Moreover, the size of segregating populations remains constrained by cost issues and is less and less commensurate with the numbers of SNPs available. Thus, guaranteeing a statistically robust marker order requires that maps include only a carefully selected subset of SNPs.

Results

In this context, the SeSAM software allows automatic genetic map construction using seriation and placement approaches, to produce (1) a high-robustness framework map which includes as many markers as possible while keeping the order robustness beyond a given statistical threshold, and (2) a high-density total map including the framework plus almost all polymorphic markers. During this process, care is taken to limit the impact of genotyping errors and of missing data on mapping quality. SeSAM can be used with a wide range of biparental populations including from outcrossing species for which phases are inferred on-the-fly by maximum-likelihood during map elongation. The package also includes functions to simulate data sets, convert data formats, detect putative genotyping errors, visualize data and map quality (including graphical genotypes), and merge several maps into a consensus. SeSAM is also suitable for interactive map construction, by providing lower-level functions for 2-point and multipoint EM analyses. The software is implemented in a R package including functions in C++.

Conclusions

SeSAM is a fully automatic linkage mapping software designed to (1) produce a framework map as robust as desired by optimizing the selection of a subset of markers, and (2) produce a high-density map including almost all polymorphic markers. The software can be used with a wide range of biparental mapping populations including cases from outcrossing. SeSAM is freely available under a GNU GPL v3 license and works on Linux, Windows, and macOS platforms. It can be downloaded together with its user-manual and quick-start tutorial from ForgeMIA (SeSAM project) at https://forgemia.inra.fr/gqe-acep/sesam/-/releases

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-05045-7.

Unlocking the hidden variation from wild repository for accelerating genetic gain in legumes

The fluctuating climates, rising human population, and deteriorating arable lands necessitate sustainable crops to fulfil global food requirements. In the countryside, legumes with intriguing but enigmatic nitrogen-fixing abilities and thriving in harsh climatic conditions promise future food security. However, breaking the yield plateau and achieving higher genetic gain are the unsolved problems of legume improvement. Present study gives emphasis on 15 important legume crops, i.e., chickpea, pigeonpea, soybean, groundnut, lentil, common bean, faba bean, cowpea, lupin, pea, green gram, back gram, horse gram, moth bean, rice bean, and some forage legumes. We have given an overview of the world and India's area, production, and productivity trends for all legume crops from 1961 to 2020. Our review article investigates the importance of gene pools and wild relatives in broadening the genetic base of legumes through pre-breeding and alien gene introgression. We have also discussed the importance of integrating genomics, phenomics, speed breeding, genetic engineering and genome editing tools in legume improvement programmes. Overall, legume breeding may undergo a paradigm shift once genomics and conventional breeding are integrated in the near future.

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.

A QTL approach in faba bean highlights the conservation of genetic control of frost tolerance among legume species

Frost is a major abiotic stress of winter type faba beans (Vica faba L.) and has adverse effects on crop yield. Climate change, far from reducing the incidence of frost events, is making these phenomena more and more common, severe, and prolonged. Despite the important interaction that the environment has in the tolerance of faba bean to frost, this trait seems to have good levels of heritability. Several QTLs for frost tolerance have already been reported, however, a more robust identification is needed to more precisely identify the genomic regions involved in faba bean tolerance to sub-zero temperatures. Several pea (Pisum sativum L.) and barrel medic (Medicago truncatula L.) frost tolerance QTLs appear to be conserved between these two species, furthering the hypothesis that the genetic control of frost tolerance in legume species might be more generally conserved. In this work, the QTL mapping in two faba bean recombinant inbred line (RIL) populations connected by a common winter-type parent has led to the identification of five genomic regions involved in the control of frost tolerance on linkage groups I, III, IV, and V. Among them, a major and robust QTL of great interest for marker-assisted selection was identified on the lower part of the long-arm of LGI. The synteny between the faba bean frost tolerance QTLs and those previously identified in other legume species such as barrel medic, pea or soybean highlighted at least partial conservation of the genetic control of frost tolerance among different faba bean genetic pools and legume species. Four novel RILs showing high and stable levels of tolerance and the ability to recover from freezing temperatures by accumulating frost tolerance QTLs are now available for breeding programs.

A Study on the Phenotypic Variation of 103 Cucumber (Cucumis sativus L.) Landraces for the Development of Desirable Cultivars Suitable for the Changing Climate

The cucumber (Cucumis sativus L.) is one of the most important vegetables in Bangladesh as well as across the globe. However, many of the important cucumber landraces have disappeared in Bangladesh due to climate change, particularly erratic rainfall, extreme temperature, salinity, and drought. Therefore, to protect against the extinction of the cucumber landraces, we collected 103 landraces in different geographical regions of Bangladesh, including drought and saline-prone areas, and studied their divergence for the future breeding programme for the development of cultivars suitable for the climate-changing situations. Data on morphological features, yield, and its components, which include 17 qualitative and quantitative traits, were recorded during the observation. Among the cucumber landraces, the Shannon–Weaver diversity index analysis revealed the presence of genetic diversity in these landraces. The biggest diversity appeared in the fruit-related characteristics, i.e., stem end fruit shape, bottom end fruit shape, fruit shape, and fruit skin colour at the table and harvest maturity. The descriptive statistics and analysis of variance expressed a wide range of variability for quantitative traits. A broad phenotypic variation was also observed for traits such as yield plant⁻¹ [CV (%) 31.88, ranges 0.96 to 3.11 kg] and fruits plant⁻¹ (CV (%), 28.71, ranges, 2.58 to 9.75). High heritability (broad sense) coupled with a high genetic gain was observed for yield and yield-contributing characteristics, indicating that these characteristics are controlled by additive gene effects, and they are more reliable for effective selection. The phenotypic correlation studies showed that fruit yield plant⁻¹ exhibited a positive and significant correlation with fruits plant⁻¹, fruit length, fruit weight, fruit width, branches plant⁻¹, and plant height. All landraces were grouped into six clusters, and the maximum number of landraces were accommodated in cluster VI (30), followed by cluster V (22), cluster III (22), cluster IV (14), cluster I (13), and cluster II (2). Comparing cluster means with studied traits revealed that cluster III with landraces AC-14, AC-97, AC-471, AC-451, and RAI-209 were more divergent for improving average fruit weight, fruit length, and fruit width. On the other hand, cluster IV with landraces AC-201, TT-161, RAI- 217, RAI-215, and TRMR-103 were more divergent for improving average vine length, internode length, and the number of primary branches plant⁻¹, the number of fruits plant⁻¹, and yield plant⁻¹. According to the MGIDI index, AC-14 (G1), AC-201 (G7), AC-471 (G24), AC-97 (G30), RAI-215 (G68) and TT-161 (G 94) may be considered to be the best parents based on their qualitative and quantitative characteristics for the future breeding programme. Moreover, crossing between the landraces, which were collected from saline and drought areas, in clusters I, V, and VI with those in other clusters could produce suitable cucumber varieties for the climatic changing situation.

Radiosensitivity of seedling traits to varying gamma doses, optimum dose determination and variation in determined doses due to different time of sowings after irradiation and methods of irradiation in faba bean genotypes

PURPOSE

Three experiments were conducted to assess the effect of different doses of gamma radiation on various seedling traits; determine the optimum doses of gamma radiation for different faba bean genotypes; find out the variation in optimum doses with respect to the different times of sowings after irradiation and methods of irradiation.

MATERIALS AND METHODS

Five faba bean genotypes viz., L-2013-060, L-2013-092, Anandnagar Local, Gazipur Local and Bangla Gangachar were used in these experiments. In Experiment I, seeds of five experimental genotypes were exposed to different doses (100 Gy 200 Gy, 300 Gy, 400 Gy, 500 Gy, 600 Gy, 700 Gy and 800 Gy) of gamma radiation and were sown immediately after irradiation. In Experiment II, seeds of Bangla Gangachar and L-2013-060 were exposed to varying doses (100-800 Gy) of gamma radiation and were sown at seven sowings starting from 0 h to 24 h at 4-h intervals after irradiation. In Experiment III, L-2013-092 genotypes was exposed to different doses (100 -800 Gy) of gamma radiation with two different methods of irradiation.

RESULTS

In Experiment I, the lethal dose 50 (LD₅₀) values have arrived at 140 Gy, 669 Gy, 575 Gy, 386 Gy and 158 Gy for L-2013-060, L-2013-092, Anandnagar Local, Gazipur Local and Bangla Gangachar, respectively. The growth reduction 50 (GR₅₀) doses for different seedling traits ranged from 130 Gy to 320 Gy for L-2013-060, 250 Gy to 480 Gy for L-2013-092, 130 Gy to 370 Gy for Anandnagar Local, 200 Gy to 350 Gy for Gazipur Local and 250 Gy to 400 Gy for Bangla Gangachar. In Experiment II, the values for LD₅₀ of the genotypes Bangla Gangachar and L-2013-060 were significantly singular for different time intervals of sowing. The values of GR₅₀ for most of the seedling traits were found to increase with the delay in sowing after irradiation from 4 to 24 h when compared with the immediately sown seed lots. In Experiment III, LD₅₀ for L-2013-092 was 337 Gy with Method 1 and 669 Gy with Method 2. In Method 1, most of the growth parameters attained GR₅₀ doses lower than Method 2. The first method was found to increase the radiosensitivity of L-2013-092.

CONCLUSION

Every experimental genotype used in these three experiments showed dose-dependent retardation of different seedling traits. These optimized doses may be employed to establish mutant populations for exploiting the novel traits of faba bean. The time of sowing after irradiation and method of irradiation was found to be essential for confirming optimum doses.

RNA-Seq Analysis of Developing Grains of Wheat to Intrigue Into the Complex Molecular Mechanism of the Heat Stress Response

Heat stress is one of the significant constraints affecting wheat production worldwide. To ensure food security for ever-increasing world population, improving wheat for heat stress tolerance is needed in the presently drifting climatic conditions. At the molecular level, heat stress tolerance in wheat is governed by a complex interplay of various heat stress-associated genes. We used a comparative transcriptome sequencing approach to study the effect of heat stress (5°C above ambient threshold temperature of 20°C) during grain filling stages in wheat genotype K7903 (Halna). At 7 DPA (days post-anthesis), heat stress treatment was given at four stages: 0, 24, 48, and 120 h. In total, 115,656 wheat genes were identified, including 309 differentially expressed genes (DEGs) involved in many critical processes, such as signal transduction, starch synthetic pathway, antioxidant pathway, and heat stress-responsive conserved and uncharacterized putative genes that play an essential role in maintaining the grain filling rate at the high temperature. A total of 98,412 Simple Sequences Repeats (SSR) were identified from de novo transcriptome assembly of wheat and validated. The miRNA target prediction from differential expressed genes was performed by psRNATarget server against 119 mature miRNA. Further, 107,107 variants including 80,936 Single nucleotide polymorphism (SNPs) and 26,171 insertion/deletion (Indels) were also identified in de novo transcriptome assembly of wheat and wheat genome Ensembl version 31. The present study enriches our understanding of known heat response mechanisms during the grain filling stage supported by discovery of novel transcripts, microsatellite markers, putative miRNA targets, and genetic variant. This enhances gene functions and regulators, paving the way for improved heat tolerance in wheat varieties, making them more suitable for production in the current climate change scenario.

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.

Genetic Dissection of Heat Stress Tolerance in Faba Bean (Vicia faba L.) Using GWAS

Heat waves are expected to become more frequent and intense, which will impact faba bean cultivation globally. Conventional breeding methods are effective but take considerable time to achieve breeding goals, and, therefore, the identification of molecular markers associated with key genes controlling heat tolerance can facilitate and accelerate efficient variety development. We phenotyped 134 accessions in six open field experiments during summer seasons at Terbol, Lebanon, at Hudeiba, Sudan, and at Central Ferry, WA, USA from 2015 to 2018. These accessions were genotyped using genotyping by sequencing (GBS), and 10,794 high quality single nucleotide polymorphisms (SNPs) were discovered. These accessions were clustered in one diverse large group, although several discrete groups may exist surrounding it. Fifteen lines belonging to different botanical groups were identified as tolerant to heat. SNPs associated with heat tolerance using single-trait (ST) and multi-trait (MT) genome-wide association studies (GWASs) showed 9 and 11 significant associations, respectively. Through the annotation of the discovered significant SNPs, we found that SNPs from transcription factor helix-loop-helix bHLH143-like S-adenosylmethionine carrier, putative pentatricopeptide repeat-containing protein At5g08310, protein NLP8-like, and photosystem II reaction center PSB28 proteins are associated with heat tolerance.

A Comprehensive Approach to Evaluate Durum Wheat-Faba Bean Mixed Crop Performance

Plant breeding for intercropping is lagging because most varieties currently available in the market are selected for sole cropping systems. The present study analyzed the response of durum wheat (12 varieties) and faba bean (3 varieties) in pure and mixed cropping. Field trials were conducted in 2019 and 2020. The performance of each variety in mixed and pure cropping was evaluated using both univariate and multivariate analyses of the grain yield and land equivalent ratio (LER). For durum wheat, grain protein content was also evaluated. Durum wheat varieties were characterized by good performance in both years, whereas faba bean varieties were more affected by the growing season, suggesting that much breeding effort is warranted to improve the latter as a pure and mixed crop. Moreover, the relative performance of all varieties was affected by their combination in mixed cropping, as evaluated based on the ratio (LERratio) between LER for wheat (LERw) and LER for faba bean (LERfb). To further evaluate the overall performance of wheat and faba bean in mixed cropping, total yield, LERtotal (LERw + LERfb), and ln(LERratio) were subjected to principal component and cluster analyses. The first principal component combined the total yield and LERtotal in a single index of the overall performance of each mixed crop combination. The second principal component, based on ln(LERratio), highlighted the relative performance of varieties in each mixed crop combination. The proposed multivariate approach can be applied in the breeding programs for intercropping to identify variety combinations based on crop performance and the relative importance of the proportion of cereal and legume grains in the total harvest.

Genomic regions associated with herbicide tolerance in a worldwide faba bean (Vicia faba L.) collection

Weeds represent one of the major constraints for faba bean crop. The identification of molecular markers associated with key genes imparting tolerance to herbicides can facilitate and fasten the efficient and effective development of herbicide tolerant cultivars. We phenotyped 140 faba bean genotypes in three open field experiments at two locations in Lebanon and Morocco against three herbicide treatments (T1 metribuzin 250 g ai/ha; T2 imazethapyr 75 g ai/ha; T3 untreated) and one in greenhouse where T1 and T3 were applied. The same set was genotyped using genotyping by sequencing (GBS) which yield 10,794 high quality single nucleotide polymorphisms (SNPs). ADMIXTURE software was used to infer the population structure which revealed two ancestral subpopulations. To identify SNPs associated with phenological and yield related traits under herbicide treatments, Single-trait (ST) and Multi-trait (MT) Genome Wide Association Studies (GWAS) were fitted using GEMMA software, showing 10 and 14 highly significant associations, respectively. Genomic sequences containing herbicide tolerance associated SNPs were aligned against the NCBI database using BLASTX tool using default parameters to annotate candidate genes underlying the causal variants. SNPs from acidic endochitinase, LRR receptor-like serine/threonine-protein kinase RCH1, probable serine/threonine-protein kinase NAK, malate dehydrogenase, photosystem I core protein PsaA and MYB-related protein P-like were significantly associated with herbicide tolerance traits.

Structural conservation of WEE1 and its role in cell cycle regulation in plants

The WEE1 kinase is ubiquitous in plant development and negatively regulates the cell cycle through phosphorylations. However, analogies with the control of the human cell cycle by tyrosine- (Tyr-) phosphorylation of cyclin-dependent kinases (CDKs) are sometimes questioned. In this in silico study, we assessed the structural conservation of the WEE1 protein in the plant kingdom with a particular focus on agronomically valuable plants, the legume crops. We analyzed the phylogenetic distribution of amino-acid sequences among a large number of plants by Bayesian analysis that highlighted the general conservation of WEE1 proteins. A detailed sequence analysis confirmed the catalytic potential of WEE1 proteins in plants. However, some substitutions of an arginine and a glutamate at the entrance of the catalytic pocket, illustrated by 3D structure predictions, challenged the specificity of this protein toward the substrate and Tyr-phosphorylation compared to the human WEE1. The structural differences, which could be responsible for the loss of specificity between human and plants, are highlighted and suggest the involvement of plant WEE1 in more cell regulation processes.

Genetic Variability of Tunisian Faba Beans (Vicia faba L.) Based on Seeds’ Morphophysical Properties as Assessed by Statistical Analysis

Faba bean (Vicia faba L.) is a legume crop cultivated for its nutritious seeds that are an important worldwide source of human food and feed. Seeds characterization is a prerequisite step for faba bean quality improvement. The morphophysical characterization of the seeds of twenty-four local faba bean accessions following the UPOV descriptors and the AOAC International standards was carried out and assessed with an approach based on Euclidean statistical model. “205 Bulk” was the unique accession harboring white hilum color seed which is linked to low convicine grain content. Irregular seed shape was the most observed character among the studied accessions except “Badii” displaying an elliptical seed shape; therefore, seed shape did not allow discrimination within our Tunisian germplasm. Interestingly, the physical characters of the seeds showed significant diversity between the accessions for all the measured parameters. A highly significant variability was observed for axial, length, and width dimensions of seeds, with “Memdouh” being the longest and largest seed accession, whereas “01-02” was the shortest and narrowest. Classification of the studied faba bean germplasm accessions based on morphophysical characters using clustering by Euclidean distance revealed three different groups. Moreover, multivariate PCA analysis further classified the faba bean accessions into four main clusters. Correlation study performed by using Spearman’s test established positive correlations within physical parameters of seeds such as between mean length and mean width of seeds. Therefore, using morphophysical parameters screening, valuable phenotypes have been selected for deeper physiological characterization and further breeding programs.

Unraveling the complexity of faba bean (Vicia faba L.) transcriptome to reveal cold-stress-responsive genes using long-read isoform sequencing technology

Faba bean (Vicia faba L.), a globally important grain legume providing a stable source of dietary protein, was one of the earliest plant cytogenetic models. However, the lack of draft genome annotations and unclear structural information on mRNA transcripts have impeded its genetic improvement. To address this, we sequenced faba bean leaf transcriptome using the PacBio single-molecule long-read isoform sequencing platform. We identified 28,569 nonredundant unigenes, ranging from 108 to 9669 bp, with a total length of 94.5 Mb. Many unigenes (3597, 12.5%) had 2-20 isoforms, indicating a highly complex transcriptome. Approximately 96.5% of the unigenes matched sequences in public databases. The predicted proteins and transcription factors included NB-ARC, Myb_domain, C3H, bHLH, and heat shock proteins, implying that this genome has an abundance of stress resistance genes. To validate our results, we selected WCOR413-15785, DHN2-12403, DHN2-14197, DHN2-14797, COR15-14478, and HVA22-15 unigenes from the ICE-CBF-COR pathway to analyze their expression patterns in cold-treated samples via qRT-PCR. The expression of dehydrin-related genes was induced by cold stress. The assembled data provide the first insights into the deep sequencing of full-length RNA from faba bean at the single-molecule level. This study provides an important foundation to improve gene modeling and protein prediction.

Development and application of the Faba_bean_130K targeted next-generation sequencing SNP genotyping platform based on transcriptome sequencing

KEY MESSAGE

Large-scale faba bean transcriptome data are available, and the first genotyping platform based on liquid-phase probe targeted capture technology was developed for genetic and molecular breeding studies. Faba bean (Vicia faba L., 2n = 12) is an important food legume crop that is widely grown for multiple uses worldwide. However, no reference genome is currently available due to its very large genome size (approximately 13 Gb) and limited single nucleotide polymorphism (SNP) markers as well as highly efficient genotyping tools have been reported for faba bean. In this study, 16.7 billion clean reads were obtained from transcriptome libraries of flowers and leaves of 102 global faba bean accessions. A total of 243,120 unigenes were de novo assembled and functionally annotated. Moreover, a total of 1,579,411 SNPs were identified and further filtered according to a selection pipeline to develop a high-throughput, flexible, low-cost Faba_bean_130K targeted next-generation sequencing (TNGS) genotyping platform. A set of 69 Chinese faba bean accessions were genotyped with the TNGS genotyping platform, and the average mapping rate of captured reads to reference transcripts was 93.14%, of which 53.23% were located in the targeted regions. The TNGS genotyping results were validated by Sanger sequencing and the average consistency rate reached 93.6%. Comprehensive population genetic analysis was performed on the 69 Chinese faba bean accessions and identified four genetic subgroups correlated with the geographic distribution. This study provides valuable genomic resources and a reliable genotyping tool that could be implemented in genetic and molecular breeding studies to accelerate new cultivar development and improvement in faba bean.

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.