Phenotypic and Genetic Characterization of the Lentil Single Plant-Derived Core Collection for Resistance to Root Rot Caused by Fusarium avenaceum.

Lentil (Lens culinaris) is a pulse crop grown for its amino acid profile, moderate drought tolerance, and ability to fix nitrogen. As the global demand for lentils expands and new production regions emerge so too have the complement of diseases that reduce yield, including the root rot complex. Although the predominant causal pathogen varies based on growing region, Fusarium avenaceum is often found to be an important contributor to disease. This study screened part of the lentil single plant-derived core collection for resistance to F. avenaceum in a greenhouse. Plants were phenotyped for disease severity using three scoring scales and the differences in biomass traits due to pathogen presence were measured. Lentil accessions varied in disease severity and differences in biomass traits were found to be correlated with each visual severity estimate (r = -0.37 to -0.63, P < 0.001), however, heritability estimates were low to moderate among traits (H2 = 0.12 to 0.43). Results of a genome-wide association study (GWAS) using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing revealed 11 quantitative trait loci (QTL) across four chromosomes. Two pairs of QTL colocated for two traits and were found near putative orthologs that have been previously associated with plant disease resistance. The identification of lentil accessions that did not exhibit a difference in biomass traits may serve as parental material in breeding or in the development of biparental mapping populations to further validate and dissect the genetic control of resistance to root rot caused by F. avenaceum.

Limited divergent adaptation despite a substantial environmental cline in wild pea.

Isolation by environment (IBE) is a widespread phenomenon in nature. It is commonly expected that the degree of difference among environments is proportional to the level of divergence between populations in their respective environments. It is therefore assumed that a species' genetic diversity displays a pattern of IBE in the presence of a strong environmental cline if gene flow does not mitigate isolation. We tested this common assumption by analysing the genetic diversity and demographic history of Pisum fulvum, which inhabits contrasting habitats in the southern Levant and is expected to display only minor migration rates between populations, making it an ideal test case. Ecogeographical and subpopulation structure were analysed and compared. The correlation of genetic with environmental distances was calculated to test the effect of isolation by distance and IBE and detect the main drivers of these effects. Historical effective population size was estimated using stairway plot. Limited overlap of ecogeographical and genetic clustering was observed, and correlation between genetic and environmental distances was statistically significant but small. We detected a sharp decline of effective population size during the last glacial period. The low degree of IBE may be the result of genetic drift due to a past bottleneck. Our findings contradict the expectation that strong environmental clines cause IBE in the absence of extensive gene flow.

Dissecting the Genetic Architecture of Aphanomyces Root Rot Resistance in Lentil by QTL Mapping and Genome-Wide Association Study.

Lentil (Lens culinaris Medikus) is an important source of protein for people in developing countries. Aphanomyces root rot (ARR) has emerged as one of the most devastating diseases affecting lentil production. In this study, we applied two complementary quantitative trait loci (QTL) analysis approaches to unravel the genetic architecture underlying this complex trait. A recombinant inbred line (RIL) population and an association mapping population were genotyped using genotyping by sequencing (GBS) to discover novel single nucleotide polymorphisms (SNPs). QTL mapping identified 19 QTL associated with ARR resistance, while association mapping detected 38 QTL and highlighted accumulation of favorable haplotypes in most of the resistant accessions. Seven QTL clusters were discovered on six chromosomes, and 15 putative genes were identified within the QTL clusters. To validate QTL mapping and genome-wide association study (GWAS) results, expression analysis of five selected genes was conducted on partially resistant and susceptible accessions. Three of the genes were differentially expressed at early stages of infection, two of which may be associated with ARR resistance. Our findings provide valuable insight into the genetic control of ARR, and genetic and genomic resources developed here can be used to accelerate development of lentil cultivars with high levels of partial resistance to ARR.

Confirmation of Fusarium root rot resistance QTL Fsp-Ps 2.1 of pea under controlled conditions.

BACKGROUND: Dry pea production has increased substantially in North America over the last few decades. With this expansion, significant yield losses have been attributed to an escalation in Fusarium root rots in pea fields. Among the most significant rot rotting pathogenic fungal species, Fusarium solani fsp. pisi (Fsp) is one of the main causal agents of root rot of pea. High levels of partial resistance to Fsp has been identified in plant genetic resources. Genetic resistance offers one of the best solutions to control this root rotting fungus. A recombinant inbred population segregating for high levels of partial resistance, previously single nucleotide polymorphism (SNP) genotyped using genotyping-by-sequencing, was phenotyped for disease reaction in replicated and repeated greenhouse trials. Composite interval mapping was deployed to identify resistance-associated quantitative trait loci (QTL). RESULTS: Three QTL were identified using three disease reaction criteria: root disease severity, ratios of diseased vs. healthy shoot heights and dry plant weights under controlled conditions using pure cultures of Fusarium solani fsp. pisi. One QTL Fsp-Ps 2.1 explains 44.4-53.4% of the variance with a narrow confidence interval of 1.2 cM. The second and third QTL Fsp-Ps3.2 and Fsp-Ps3.3 are closely linked and explain only 3.6-4.6% of the variance. All of the alleles are contributed by the resistant parent PI 180693. CONCLUSION: With the confirmation of Fsp-Ps 2.1 now in two RIL populations, SNPs associated with this region make a good target for marker-assisted selection in pea breeding programs to obtain high levels of partial resistance to Fusarium root rot caused by Fusarium solani fsp. pisi.

Genome-wide SNP identification, linkage map construction and QTL mapping for seed mineral concentrations and contents in pea (Pisum sativum L.).

BACKGROUND: Marker-assisted breeding is now routinely used in major crops to facilitate more efficient cultivar improvement. This has been significantly enabled by the use of next-generation sequencing technology to identify loci and markers associated with traits of interest. While rich in a range of nutritional components, such as protein, mineral nutrients, carbohydrates and several vitamins, pea (Pisum sativum L.), one of the oldest domesticated crops in the world, remains behind many other crops in the availability of genomic and genetic resources. To further improve mineral nutrient levels in pea seeds requires the development of genome-wide tools. The objectives of this research were to develop these tools by: identifying genome-wide single nucleotide polymorphisms (SNPs) using genotyping by sequencing (GBS); constructing a high-density linkage map and comparative maps with other legumes, and identifying quantitative trait loci (QTL) for levels of boron, calcium, iron, potassium, magnesium, manganese, molybdenum, phosphorous, sulfur, and zinc in the seed, as well as for seed weight. RESULTS: In this study, 1609 high quality SNPs were found to be polymorphic between 'Kiflica' and 'Aragorn', two parents of an F6-derived recombinant inbred line (RIL) population. Mapping 1683 markers including 75 previously published markers and 1608 SNPs developed from the present study generated a linkage map of size 1310.1 cM. Comparative mapping with other legumes demonstrated that the highest level of synteny was observed between pea and the genome of Medicago truncatula. QTL analysis of the RIL population across two locations revealed at least one QTL for each of the mineral nutrient traits. In total, 46 seed mineral concentration QTLs, 37 seed mineral content QTLs, and 6 seed weight QTLs were discovered. The QTLs explained from 2.4% to 43.3% of the phenotypic variance. CONCLUSION: The genome-wide SNPs and the genetic linkage map developed in this study permitted QTL identification for pea seed mineral nutrients that will serve as important resources to enable marker-assisted selection (MAS) for nutritional quality traits in pea breeding programs.

Association mapping of starch chain length distribution and amylose content in pea (Pisum sativum L.) using carbohydrate metabolism candidate genes.

BACKGROUND: Although starch consists of large macromolecules composed of glucose units linked by α-1,4-glycosidic linkages with α-1,6-glycosidic branchpoints, variation in starch structural and functional properties is found both within and between species. Interest in starch genetics is based on the importance of starch in food and industrial processes, with the potential of genetics to provide novel starches. The starch metabolic pathway is complex but has been characterized in diverse plant species, including pea. RESULTS: To understand how allelic variation in the pea starch metabolic pathway affects starch structure and percent amylose, partial sequences of 25 candidate genes were characterized for polymorphisms using a panel of 92 diverse pea lines. Variation in the percent amylose composition of extracted seed starch and (amylopectin) chain length distribution, one measure of starch structure, were characterized for these lines. Association mapping was undertaken to identify polymorphisms associated with the variation in starch chain length distribution and percent amylose, using a mixed linear model that incorporated population structure and kinship. Associations were found for polymorphisms in seven candidate genes plus Mendel's r locus (which conditions the round versus wrinkled seed phenotype). The genes with associated polymorphisms are involved in the substrate supply, chain elongation and branching stages of the pea carbohydrate and starch metabolic pathways. CONCLUSIONS: The association of polymorphisms in carbohydrate and starch metabolic genes with variation in amylopectin chain length distribution and percent amylose may help to guide manipulation of pea seed starch structural and functional properties through plant breeding.

Genome-wide association mapping of partial resistance to Aphanomyces euteiches in pea.

BACKGROUND: Genome-wide association (GWA) mapping has recently emerged as a valuable approach for refining the genetic basis of polygenic resistance to plant diseases, which are increasingly used in integrated strategies for durable crop protection. Aphanomyces euteiches is a soil-borne pathogen of pea and other legumes worldwide, which causes yield-damaging root rot. Linkage mapping studies reported quantitative trait loci (QTL) controlling resistance to A. euteiches in pea. However the confidence intervals (CIs) of these QTL remained large and were often linked to undesirable alleles, which limited their application in breeding. The aim of this study was to use a GWA approach to validate and refine CIs of the previously reported Aphanomyces resistance QTL, as well as identify new resistance loci. METHODS: A pea-Aphanomyces collection of 175 pea lines, enriched in germplasm derived from previously studied resistant sources, was evaluated for resistance to A. euteiches in field infested nurseries in nine environments and with two strains in climatic chambers. The collection was genotyped using 13,204 SNPs from the recently developed GenoPea Infinium® BeadChip. RESULTS: GWA analysis detected a total of 52 QTL of small size-intervals associated with resistance to A. euteiches, using the recently developed Multi-Locus Mixed Model. The analysis validated six of the seven previously reported main Aphanomyces resistance QTL and detected novel resistance loci. It also provided marker haplotypes at 14 consistent QTL regions associated with increased resistance and highlighted accumulation of favourable haplotypes in the most resistant lines. Previous linkages between resistance alleles and undesired late-flowering alleles for dry pea breeding were mostly confirmed, but the linkage between loci controlling resistance and coloured flowers was broken due to the high resolution of the analysis. A high proportion of the putative candidate genes underlying resistance loci encoded stress-related proteins and others suggested that the QTL are involved in diverse functions. CONCLUSION: This study provides valuable markers, marker haplotypes and germplasm lines to increase levels of partial resistance to A. euteiches in pea breeding.

Multiple post-domestication origins of kabuli chickpea through allelic variation in a diversification-associated transcription factor.

Chickpea (Cicer arietinum) is among the founder crops domesticated in the Fertile Crescent. One of two major forms of chickpea, the so-called kabuli type, has white flowers and light-colored seed coats, properties not known to exist in the wild progenitor. The origin of the kabuli form has been enigmatic. We genotyped a collection of wild and cultivated chickpea genotypes with 538 single nucleotide polymorphisms (SNPs) and examined patterns of molecular diversity relative to geographical sources and market types. In addition, we examined sequence and expression variation in candidate anthocyanin biosynthetic pathway genes. A reduction in genetic diversity and extensive genetic admixture distinguish cultivated chickpea from its wild progenitor species. Among germplasm, the kabuli form is polyphyletic. We identified a basic helix-loop-helix (bHLH) transcription factor at chickpea's B locus that conditions flower and seed colors, orthologous to Mendel's A gene of garden pea, whose loss of function is associated invariantly with the kabuli type of chickpea. From the polyphyletic distribution of the kabuli form in germplasm, an absence of nested variation within the bHLH gene and invariant association of loss of function of bHLH among the kabuli type, we conclude that the kabuli form arose multiple times during the phase of phenotypic diversification after initial domestication of cultivated chickpea.

From Mendel's discovery on pea to today's plant genetics and breeding : Commemorating the 150th anniversary of the reading of Mendel's discovery.

KEY MESSAGE: This work discusses several selected topics of plant genetics and breeding in relation to the 150th anniversary of the seminal work of Gregor Johann Mendel. In 2015, we celebrated the 150th anniversary of the presentation of the seminal work of Gregor Johann Mendel. While Darwin's theory of evolution was based on differential survival and differential reproductive success, Mendel's theory of heredity relies on equality and stability throughout all stages of the life cycle. Darwin's concepts were continuous variation and "soft" heredity; Mendel espoused discontinuous variation and "hard" heredity. Thus, the combination of Mendelian genetics with Darwin's theory of natural selection was the process that resulted in the modern synthesis of evolutionary biology. Although biology, genetics, and genomics have been revolutionized in recent years, modern genetics will forever rely on simple principles founded on pea breeding using seven single gene characters. Purposeful use of mutants to study gene function is one of the essential tools of modern genetics. Today, over 100 plant species genomes have been sequenced. Mapping populations and their use in segregation of molecular markers and marker-trait association to map and isolate genes, were developed on the basis of Mendel's work. Genome-wide or genomic selection is a recent approach for the development of improved breeding lines. The analysis of complex traits has been enhanced by high-throughput phenotyping and developments in statistical and modeling methods for the analysis of phenotypic data. Introgression of novel alleles from landraces and wild relatives widens genetic diversity and improves traits; transgenic methodologies allow for the introduction of novel genes from diverse sources, and gene editing approaches offer possibilities to manipulate gene in a precise manner.

Large-scale microsatellite development in grasspea (Lathyrus sativus L.), an orphan legume of the arid areas.

BACKGROUND: Grasspea (Lathyrus sativus L., 2n = 14), a member of the family Leguminosae, holds great agronomic potential as grain and forage legume crop in the arid areas for its superb resilience to abiotic stresses such as drought, flood and salinity. The crop could not make much progress through conventional breeding in the past, and there are hardly any detailed molecular biology studies due to paucity of reliable molecular markers representative of the entire genome. RESULTS: Using the 454 FLX Titanium pyrosequencing technique, 651,827 simple sequence repeat (SSR) loci were identified and 50,144 nonredundant primer pairs were successfully designed, of which 288 were randomly selected for validation among 23 L. sativus and one L. cicera accessions of diverse provenance. 74 were polymorphic, 70 monomorphic, and 144 with no PCR product. The number of observed alleles ranged from two to five, the observed heterozygosity from 0 to 0.9545, and Shannon's information index ranged from 0.1013 to 1.0980, respectively. The dendrogram constructed by using unweighted pair group method with arithmetic mean (UPGMA) based on Nei's genetic distance, showed obvious distinctions and understandable relationships among the 24 accessions. CONCLUSIONS: The large number of SSR primer pairs developed in this study would make a significant contribution to genomics enabled improvement of grasspea.

Cicer super-pangenome provides insights into species evolution and agronomic trait loci for crop improvement in chickpea.

Chickpea (Cicer arietinum L.)-an important legume crop cultivated in arid and semiarid regions-has limited genetic diversity. Efforts are being undertaken to broaden its diversity by utilizing its wild relatives, which remain largely unexplored. Here, we present the Cicer super-pangenome based on the de novo genome assemblies of eight annual Cicer wild species. We identified 24,827 gene families, including 14,748 core, 2,958 softcore, 6,212 dispensable and 909 species-specific gene families. The dispensable genome was enriched for genes related to key agronomic traits. Structural variations between cultivated and wild genomes were used to construct a graph-based genome, revealing variations in genes affecting traits such as flowering time, vernalization and disease resistance. These variations will facilitate the transfer of valuable traits from wild Cicer species into elite chickpea varieties through marker-assisted selection or gene-editing. This study offers valuable insights into the genetic diversity and potential avenues for crop improvement in chickpea.

QTL meta-analysis provides a comprehensive view of loci controlling partial resistance to Aphanomyces euteiches in four sources of resistance in pea.

BACKGROUND: Development of durable plant genetic resistance to pathogens through strategies of QTL pyramiding and diversification requires in depth knowledge of polygenic resistance within the available germplasm. Polygenic partial resistance to Aphanomyces root rot, caused by Aphanomyces euteiches, one of the most damaging pathogens of pea worldwide, was previously dissected in individual mapping populations. However, there are no data available regarding the diversity of the resistance QTL across a broader collection of pea germplasm. In this study, we performed a meta-analysis of Aphanomyces root rot resistance QTL in the four main sources of resistance in pea and compared their genomic localization with genes/QTL controlling morphological or phenological traits and with putative candidate genes. RESULTS: Meta-analysis, conducted using 244 individual QTL reported previously in three mapping populations (Puget x 90-2079, Baccara x PI180693 and Baccara x 552) and in a fourth mapping population in this study (DSP x 90-2131), resulted in the identification of 27 meta-QTL for resistance to A. euteiches. Confidence intervals of meta-QTL were, on average, reduced four-fold compared to mean confidence intervals of individual QTL. Eleven consistent meta-QTL, which highlight seven highly consistent genomic regions, were identified. Few meta-QTL specificities were observed among mapping populations, suggesting that sources of resistance are not independent. Seven resistance meta-QTL, including six of the highly consistent genomic regions, co-localized with six of the meta-QTL identified in this study for earliness and plant height and with three morphological genes (Af, A, R). Alleles contributing to the resistance were often associated with undesirable alleles for dry pea breeding. Candidate genes underlying six main meta-QTL regions were identified using colinearity between the pea and Medicago truncatula genomes. CONCLUSIONS: QTL meta-analysis provided an overview of the moderately low diversity of loci controlling partial resistance to A. euteiches in four main sources of resistance in pea. Seven highly consistent genomic regions with potential use in marker-assisted-selection were identified. Confidence intervals at several main QTL regions were reduced and co-segregation among resistance and morphological/phenological alleles was identified. Further work will be required to identify the best combinations of QTL for durably increasing partial resistance to A. euteiches.

Population structure and linkage disequilibrium in Lupinus albus L. germplasm and its implication for association mapping.

White lupin (Lupinus albus L.) has been around since 300 B.C. and is recognized for its ability to grow on poor soils and application as green manure in addition to seed harvest. The seed has very high levels of protein (33-47 %) and oil (6-13 %). It also has many secondary metabolites that are potentially of nutraceutical value to animals and humans. Despite such a great potential, lupins role in modern agriculture began only in the twentieth century. Although a large collection of Lupinus germplasm accessions is available worldwide, rarely have they been genetically characterized. Additionally, scarce genomic resources in terms of recombinant populations and genome information have been generated for L. albus. With the advancement in association mapping methods, the natural populations have the potential to replace the recombinant populations in gene mapping and marker-trait associations. Therefore, we studied the genetic similarity, population structure and marker-trait association in a USDA germplasm collection for their current and future application in this crop improvement. A total of 122 PI (Plant Inventory) lines were screened with 18 AFLP primer pairs that generated 2,277 fragments. A subset of 892 polymorphic markers with MAF >0.05 (minor allele frequency) were used for association mapping. The cluster analysis failed to group accessions on the basis of their passport information, and a weak structure and low linkage disequilibrium (LD) were observed indicating the usefulness of the collection for association mapping. Moreover, we were also able to identify two markers (a p value of 1.53 × 10(-4) and 2.3 × 10(-4)) that explained 22.69 and 20.5 % of seed weight variation determined using R (LR) (2) . The implications of lack of geographic clustering, population structure, low LD and the ability of AFLP to map seed weight trait using association mapping and the usefulness of the PI collections in breeding programs are discussed.

How Could the Use of Crop Wild Relatives in Breeding Increase the Adaptation of Crops to Marginal Environments?

Alongside the use of fertilizer and chemical control of weeds, pests, and diseases modern breeding has been very successful in generating cultivars that have increased agricultural production several fold in favorable environments. These typically homogeneous cultivars (either homozygous inbreds or hybrids derived from inbred parents) are bred under optimal field conditions and perform well when there is sufficient water and nutrients. However, such optimal conditions are rare globally; indeed, a large proportion of arable land could be considered marginal for agricultural production. Marginal agricultural land typically has poor fertility and/or shallow soil depth, is subject to soil erosion, and often occurs in semi-arid or saline environments. Moreover, these marginal environments are expected to expand with ongoing climate change and progressive degradation of soil and water resources globally. Crop wild relatives (CWRs), most often used in breeding as sources of biotic resistance, often also possess traits adapting them to marginal environments. Wild progenitors have been selected over the course of their evolutionary history to maintain their fitness under a diverse range of stresses. Conversely, modern breeding for broad adaptation has reduced genetic diversity and increased genetic vulnerability to biotic and abiotic challenges. There is potential to exploit genetic heterogeneity, as opposed to genetic uniformity, in breeding for the utilization of marginal lands. This review discusses the adaptive traits that could improve the performance of cultivars in marginal environments and breeding strategies to deploy them.

Corrigendum: How could the use of crop wild relatives in breeding increase the adaptation of crops to marginal environments?

[This corrects the article .].

New consistent QTL in pea associated with partial resistance to Aphanomyces euteiches in multiple French and American environments.

Partial resistances, often controlled by quantitative trait loci (QTL), are considered to be more durable than monogenic resistances. Therefore, a precursor to developing efficient breeding programs for polygenic resistance to pathogens should be a greater understanding of genetic diversity and stability of resistance QTL in plants. In this study, we deciphered the diversity and stability of resistance QTL to Aphanomyces euteiches in pea towards pathogen variability, environments and scoring criteria, from two new sources of partial resistance (PI 180693 and 552), effective in French and USA infested fields. Two mapping populations of 178 recombinant inbred lines each, derived from crosses between 552 or PI 180693 (partially resistant) and Baccara (susceptible), were used to identify QTL for Aphanomyces root rot resistance in controlled and in multiple French and USA field conditions using several resistance criteria. We identified a total of 135 additive-effect QTL corresponding to 23 genomic regions and 13 significant epistatic interactions associated with partial resistance to A. euteiches in pea. Among the 23 additive-effect genomic regions identified, five were consistently detected, and showed highly stable effects towards A. euteiches strains, environments, resistance criteria, condition tests and RIL populations studied. These results confirm the complexity of inheritance of partial resistance to A. euteiches in pea and provide good bases for the choice of consistent QTL to use in marker-assisted selection schemes to increase current levels of resistance to A. euteiches in pea breeding programs.

Integrating Genomic and Phenomic Approaches to Support Plant Genetic Resources Conservation and Use.

Plant genebanks provide genetic resources for breeding and research programs worldwide. These programs benefit from having access to high-quality, standardized phenotypic and genotypic data. Technological advances have made it possible to collect phenomic and genomic data for genebank collections, which, with the appropriate analytical tools, can directly inform breeding programs. We discuss the importance of considering genebank accession homogeneity and heterogeneity in data collection and documentation. Citing specific examples, we describe how well-documented genomic and phenomic data have met or could meet the needs of plant genetic resource managers and users. We explore future opportunities that may emerge from improved documentation and data integration among plant genetic resource information systems.

Harnessing Genetic Diversity in the USDA Pea Germplasm Collection Through Genomic Prediction.

Phenotypic evaluation and efficient utilization of germplasm collections can be time-intensive, laborious, and expensive. However, with the plummeting costs of next-generation sequencing and the addition of genomic selection to the plant breeder's toolbox, we now can more efficiently tap the genetic diversity within large germplasm collections. In this study, we applied and evaluated genomic prediction's potential to a set of 482 pea (Pisum sativum L.) accessions-genotyped with 30,600 single nucleotide polymorphic (SNP) markers and phenotyped for seed yield and yield-related components-for enhancing selection of accessions from the USDA Pea Germplasm Collection. Genomic prediction models and several factors affecting predictive ability were evaluated in a series of cross-validation schemes across complex traits. Different genomic prediction models gave similar results, with predictive ability across traits ranging from 0.23 to 0.60, with no model working best across all traits. Increasing the training population size improved the predictive ability of most traits, including seed yield. Predictive abilities increased and reached a plateau with increasing number of markers presumably due to extensive linkage disequilibrium in the pea genome. Accounting for population structure effects did not significantly boost predictive ability, but we observed a slight improvement in seed yield. By applying the best genomic prediction model (e.g., RR-BLUP), we then examined the distribution of genotyped but nonphenotyped accessions and the reliability of genomic estimated breeding values (GEBV). The distribution of GEBV suggested that none of the nonphenotyped accessions were expected to perform outside the range of the phenotyped accessions. Desirable breeding values with higher reliability can be used to identify and screen favorable germplasm accessions. Expanding the training set and incorporating additional orthogonal information (e.g., transcriptomics, metabolomics, physiological traits, etc.) into the genomic prediction framework can enhance prediction accuracy.

Environmental and genetic determinants of amphicarpy in Pisum fulvum, a wild relative of domesticated pea.

Pisum fulvum is an annual legume native to Syria, Lebanon, Israel and Jordan. In certain locations, P. fulvum individuals were documented to display a reproductive dimorphism - amphicarpy, with both above and below ground flowers and pods. Herein we aimed to study the possible role of soil texture on amphicarpy in P. fulvum, to investigate the possible bio-climatic associations of P. fulvum amphicarpy and to identify genetic markers associated with this phenotype. A set of 127 germplasm accessions sampled across the Israeli distribution range of the species was phenotyped in two common garden nurseries. Land use and bioclimatic data were used to delineate the eco-geographic clustering of accession's sampling sites. Single nucleotide polymorphism (SNP) markers were employed in genome-wide association study to identify associated loci. Amphicarpy was subject to strong experimental site x genotype interaction with higher phenotypic expression in fine textured soil relative to sandy loam. Amphicarpy was more prevalent among accessions sampled in eastern Judea and Samaria and was weakly associated with early phenology and relatively modest above ground biomass production. Twelve SNP markers were significantly associated with amphicarpy, each explaining between 8 and 12 % of the phenotypic variation. In P. fulvum amphicarpy seems to be a polygenetic trait controlled by an array of genes that is likely to be affected by environmental stimuli. The probable selective advantage of the association between amphicarpy and early flowering is in line with its relative prevalence in drought prone territories subject to heavy grazing.