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.

The adaptability and stability of 37 faba bean (Vicia faba L.) accessions with different levels of tolerance to metribuzin or imazethapyr was assessed across 12 season–location–herbicide experiments. Significant Genotype x environment (GE) interaction was found for the days to flowering (DFLR), plant height (PLHT) and grain yield (GY). Performance and stability of the accessions regarding PLHT and GY were assessed using four different stability parameters: cultivar superiority, static stability, Wricke’s eco-valence and Finlay and Wilkinson’s regression model. The stability parameters ranked these genotypes differently suggesting that PLHT and GY stability should be assessed not only on a single or a few stability parameters but on a combination of them. GGE biplot analysis indicated that the environments representing metribuzin treatment at Marchouch 2014–2015 and the non-treated treatment at Terbol 2018–2019 are the ideal environments for evaluating faba bean genotypes. GGE biplots showed herbicide tolerant accession IG12983 with simultaneous average PLHT, GY and stability across the environments. The performance of other tolerant accessions, namely IG13945, IG13906, IG106453, FB2648, and FB1216 was less stable but superior under specific mega environments. Therefore, utilizing these accessions in faba bean breeding programs would help broaden the adaptability to diverse locations–season–herbicide treatments.

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.

Weeds cause serious constraint to faba bean (Vicia fabaL.) productivity. Broad weed control is hampered by the availability of postemergence herbicides to control them, as the current faba bean cultivars are highly susceptible to them. Therefore, the deployment of postemergence herbicide tolerance is desirable in faba bean. To address this, 130 accessions were screened for their response in mature plants under field conditions to the recommended dosage of two herbicides, metribuzin at 250 g a.i. ha⁻¹and imazethapyr at 75 g a.i. ha⁻¹at Marchouch and Terbol stations. The recorded herbicide damage score (HDS) varied from 1 (no visual damage) to 5 (full damage with death of more than 50% of plants) at both locations. Low but highly significant (p< .01) and positive correlation (+0.26) was obtained between the recorded HDS at both locations. Both herbicides significantly delayed flowering and maturity time occurrence, reduced plant height and grain yield, and increased number of branches. Reduction index (RI) correlated positively with HDS score at Terbol station in different seasons and at Marchouch in 2016–2017 seasons. Eleven tolerant accessions were identified and further evaluated to 1×, 1.5×, and 2× of recommended dose of both herbicides. The results indicated that the harmful effect of herbicides on grain yield reduction intensified from 13.4 to 27.2% and from −7.6 to 1.8% as the dose of metribuzin and imazethapyr increased respectively from 250 to 500 g a.i. ha⁻¹and from 75 to 150 g a.i. ha⁻¹. Tolerance to metribuzin and imazethapyr in eight faba bean accessions was confirmed with no significant reduction in grain yield.

Genomic selection (GS) has potential to accelerate the genetic gain in crop plants. This study was undertaken to assess the accuracy and potential of GS in faba bean [ Vicia faba (L.)] and to enhance its application in breeding programs. A set of 118 diverse faba bean accessions were phenotyped for key agronomic traits under herbicide and heat stress across 16 environments in Morocco, Lebanon, Sudan, and the United States. These accessions were genotyped, revealing 170 single nucleotide polymorphisms (SNPs) strongly associated with target traits. kompetitive allele‐specific PCR (KASP) markers were subsequently designed and validated on 4512 diverse breeding lines. Prediction accuracy (PA) was assessed using the reproducing kernel Hilbert space model, with and without genotype‐by‐environment interactions and taking into consideration two cross‐validation (CV) strategies: CV1 (predicting new lines) and CV2 (predicting complete records from unbalanced data). Additionally, 75 KASP markers with known association with heat tolerance traits were prioritized to estimate the PA of the models. The results showed a comparable PA between the two models, with CV1 outperforming CV2. This highlighted the difficulty in predicting the performance of untested lines in tested environments compared to lines evaluated in some environments but not others. Moreover, SNP subset size and composition significantly impacted PA, especially under heat stress. The highest accuracies were observed for days to flowering and plant height under heat stress and for plant height and grain yield under herbicide environments, indicating that these traits are ideal for training population selection. Optimizing the size and composition of the training population will enhance the effectiveness of GS in faba bean breeding.

Assessing the adaptability and stability of herbicide-tolerant lentil accessions to two broad-spectrum post-emergence herbicides in multi-environment trials has become a must in a breeding program to improve its selection. The adaptability and stability of 42 herbicide-tolerant lentil accessions were investigated using five stability parameters under eight different environments. Significant Genotype–Environment (GE) interaction was found for days to flowering (DFLR), days to maturity (DMAT), and seed yield per plant (SY). The analyzed stability parameters such as Cultivar superiority, Finlay–Wilkinson, Shukla, Static Stability, and Wricke’s Ecovalence ranked the tested accessions differently, confirming the importance of using a combination of stability parameters when evaluating the performance of a group of accessions. GGE biplot of the SY trait accounted for 60.79% of sums of squares of the GE interaction and showed that cool and high rainfall environments are ideal for testing the agronomic performance of tolerant accessions. The GGE biplot of SY showed that IG4605(19), IG195(6), and IG156635(12) were specifically adapted to one mega environment, whereas IG70056(38) was identified as a superior line having a high and stable yield. These lines should be included in lentil crossing programs to develop herbicide-tolerant cultivars adapted to diverse environments.