Genetic structure and ecological niche space of lentil's closest wild relative, Lens orientalis (Boiss.) Schmalh

Crops arose from wild ancestors and to understand their domestication it is essential to compare the cultivated species with their crop wild relatives. These represent an important source of further crop improvement, in particular in relation to climate change. Although there are about 58,000 Lens accessions held in genebanks, only 1% are wild. We examined the geographic distribution and genetic diversity of the lentil's immediate progenitor L. orientalis. We used Genotyping by Sequencing (GBS) to identify and characterize differentiation among accessions held at germplasm collections. We then determined whether genetically distinct clusters of accessions had been collected from climatically distinct locations. Of the 195 genotyped accessions, 124 were genuine L. orientalis with four identified genetic groups. Although an environmental distance matrix was significantly correlated with geographic distance in a Mantel test, the four identified genetic clusters were not found to occupy significantly different environmental space. Maxent modelling gave a distinct predicted distribution pattern centred in the Fertile Crescent, with intermediate probabilities of occurrence in parts of Turkey, Greece, Cyprus, Morocco, and the south of the Iberian Peninsula with NW Africa. Future projections did not show any dramatic alterations in the distribution according to the climate change scenarios tested. We have found considerable diversity in L. orientalis, some of which track climatic variability. The results of the study showed the genetic diversity of wild lentil and indicate the importance of ongoing collections and in situ conservation for our future capacity to harness the genetic variation of the lentil progenitor.

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

[This corrects the article .].

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 (H² = 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.

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.

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.

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.

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.

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.

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.

A reference genome for pea provides insight into legume genome evolution

We report the first annotated chromosome-level reference genome assembly for pea, Gregor Mendel's original genetic model. Phylogenetics and paleogenomics show genomic rearrangements across legumes and suggest a major role for repetitive elements in pea genome evolution. Compared to other sequenced Leguminosae genomes, the pea genome shows intense gene dynamics, most likely associated with genome size expansion when the Fabeae diverged from its sister tribes. During Pisum evolution, translocation and transposition differentially occurred across lineages. This reference sequence will accelerate our understanding of the molecular basis of agronomically important traits and support crop improvement.

Advanced Imaging for Quantitative Evaluation of Aphanomyces Root Rot Resistance in Lentil

Aphanomyces root rot (ARR) is a soil-borne disease that results in severe yield losses in lentil. The development of resistant cultivars is one of the key strategies to control this pathogen. However, the evaluation of disease severity is limited to visual scores that can be subjective. This study utilized image-based phenotyping approaches to evaluate Aphanomyces euteiches resistance in lentil genotypes in greenhouse (351 genotypes from lentil single plant/LSP derived collection and 191 genotypes from recombinant inbred lines/RIL using digital Red-Green-Blue/RGB and hyperspectral imaging) and field (173 RIL genotypes using unmanned aerial system-based multispectral imaging) conditions. Moderate to strong correlations were observed between RGB, multispectral, and hyperspectral derived features extracted from lentil shoots/roots and visual scores. In general, root features extracted from RGB imaging were found to be strongly associated with disease severity. With only three root traits, elastic net regression model was able to predict disease severity across and within multiple datasets (R 2 = 0.45-0.73 and RMSE = 0.66-1.00). The selected features could represent visual disease scores. Moreover, we developed twelve normalized difference spectral indices (NDSIs) that were significantly correlated with disease scores: two NDSIs for lentil shoot section - computed from wavelengths of 1170, 1160, 1270, and 1280 nm (0.12 ≤ |r| ≤ 0.24, P < 0.05) and ten NDSIs for lentil root sections - computed from wavelengths in the range of 630-670, 700-840, and 1320-1530 nm (0.10 ≤ |r| ≤ 0.50, P < 0.05). Root-derived NDSIs were more accurate in predicting disease scores with an R 2 of 0.54 (RMSE = 0.86), especially when the model was trained and tested on LSP accessions, compared to R 2 of 0.25 (RMSE = 1.64) when LSP and RIL genotypes were used as train and test datasets, respectively. Importantly, NDSIs - computed from wavelengths of 700, 710, 730, and 790 nm - had strong positive correlations with disease scores (0.35 ≤r ≤ 0.50, P < 0.0001), which was confirmed in field phenotyping with similar correlations using vegetation index with red edge wavelength (normalized difference red edge, 0.36 ≤ |r| ≤ 0.57, P < 0.0001). The adopted image-based phenotyping approaches can help plant breeders to objectively quantify ARR resistance and reduce the subjectivity in selecting potential genotypes.

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.

Molecular Evidence for Two Domestication Events in the Pea Crop

Pea, one of the founder crops from the Near East, has two wild species: Pisum sativum subsp. elatius, with a wide distribution centered in the Mediterranean, and P. fulvum, which is restricted to Syria, Lebanon, Israel, Palestine and Jordan. Using genome wide analysis of 11,343 polymorphic single nucleotide polymorphisms (SNPs) on a set of wild P. elatius (134) and P. fulvum (20) and 74 domesticated accessions (64 P. sativum landraces and 10 P. abyssinicum), we demonstrated that domesticated P. sativum and the Ethiopian pea (P. abyssinicum) were derived from different P. elatius genepools. Therefore, pea has at least two domestication events. The analysis does not support a hybrid origin of P. abyssinicum, which was likely introduced into Ethiopia and Yemen followed by eco-geographic adaptation. Both P. sativum and P. abyssinicum share traits that are typical of domestication, such as non-dormant seeds. Non-dormant seeds were also found in several wild P. elatius accessions which could be the result of crop to wild introgression or natural variation that may have been present during pea domestication. A sub-group of P. elatius overlaps with P. sativum landraces. This may be a consequence of bidirectional gene-flow or may suggest that this group of P. elatius is the closest extant wild relative of P. sativum.

Genomic diversity and macroecology of the crop wild relatives of domesticated pea

There is growing interest in the conservation and utilization of crop wild relatives (CWR) in international food security policy and research. Legumes play an important role in human health, sustainable food production, global food security, and the resilience of current agricultural systems. Pea belongs to the ancient set of cultivated plants of the Near East domestication center and remains an important crop today. Based on genome-wide analysis, P. fulvum was identified as a well-supported species, while the diversity of wild P. sativum subsp. elatius was structured into 5 partly geographically positioned clusters. We explored the spatial and environmental patterns of two progenitor species of domesticated pea in the Mediterranean Basin and in the Fertile Crescent in relation to the past and current climate. This study revealed that isolation by distance does not explain the genetic structure of P. sativum subsp. elatius in its westward expansion from its center of origin. The genetic diversity of wild pea may be driven by Miocene-Pliocene events, while the phylogenetic diversity centers may reflect Pleisto-Holocene climatic changes. These findings help set research and discussion priorities and provide geographical and ecological information for germplasm-collecting missions, as well as for the preservation of extant diversity in ex-situ collections.

Selecting Lentil Accessions for Global Selenium Biofortification

The biofortification of lentil (Lens culinaris Medikus.) has the potential to provide adequate daily selenium (Se) to human diets. The objectives of this study were to (1) determine how low-dose Se fertilizer application at germination affects seedling biomass, antioxidant activity, and Se uptake of 26 cultivated lentil genotypes; and (2) quantify the seed Se concentration of 191 lentil wild accessions grown in Terbol, Lebanon. A germination study was conducted with two Se treatments [0 (control) and 30 kg of Se/ha] with three replicates. A separate field study was conducted in Lebanon for wild accessions without Se fertilizer. Among cultivated lentil accessions, PI533690 and PI533693 showed >100% biomass increase vs. controls. Se addition significantly increased seedling Se uptake, with the greatest uptake (6.2 µg g⁻¹) by PI320937 and the least uptake (1.1 µg g⁻¹) by W627780. Seed Se concentrations of wild accessions ranged from 0 to 2.5 µg g⁻¹; accessions originating from Syria (0–2.5 µg g⁻¹) and Turkey (0–2.4 µg g⁻¹) had the highest seed Se. Frequency distribution analysis revealed that seed Se for 63% of accessions was between 0.25 and 0.75 µg g⁻¹, and thus a single 50 g serving of lentil has the potential to provide adequate dietary Se (20–60% of daily recommended daily allowance). As such, Se application during plant growth for certain lentil genotypes grown in low Se soils may be a sustainable Se biofortification solution to increase seed Se concentration. Incorporating a diverse panel of lentil wild germplasm into Se biofortification programs will increase genetic diversity for effective genetic mapping for increased lentil seed Se nutrition and plant productivity.

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

A community resource for exploring and utilizing genetic diversity in the USDA pea single plant plus collection

Globally, pea (Pisum sativum L.) is an important temperate legume crop for food, feed and fodder, and many breeding programs develop cultivars adapted to these end-uses. In order to assist pea development efforts, we assembled the USDA Pea Single Plant Plus Collection (PSPPC), which contains 431 P. sativum accessions with morphological, geographic and taxonomic diversity. The collection was characterized genetically in order to maximize its value for trait mapping and genomics-assisted breeding. To that end, we used genotyping-by-sequencing-a cost-effective method for de novo single-nucleotide polymorphism (SNP) marker discovery-to generate 66 591 high-quality SNPs. These data facilitated the identification of accessions divergent from mainstream breeding germplasm that could serve as sources of novel, favorable alleles. In particular, a group of accessions from Central Asia appear nearly as diverse as a sister species, P. fulvum, and subspecies, P. sativum subsp. elatius. PSPPC genotypes can be paired with new and existing phenotype data for trait mapping; as proof-of-concept, we localized Mendel's A gene controlling flower color to its known position. We also used SNP data to define a smaller core collection of 108 accessions with similar levels of genetic diversity as the entire PSPPC, resulting in a smaller germplasm set for research screening and evaluation under limited resources. Taken together, the results presented in this study along with the release of a publicly available SNP data set comprise a valuable resource for supporting worldwide pea genetic improvement efforts.

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.

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.

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.

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

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

Identification of Quantitative Trait Loci Controlling Root and Shoot Traits Associated with Drought Tolerance in a Lentil (Lens culinaris Medik.) Recombinant Inbred Line Population

Drought is one of the major abiotic stresses limiting lentil productivity in rainfed production systems. Specific rooting patterns can be associated with drought avoidance mechanisms that can be used in lentil breeding programs. In all, 252 co-dominant and dominant markers were used for Quantitative Trait Loci (QTL) analysis on 132 lentil recombinant inbred lines based on greenhouse experiments for root and shoot traits during two seasons under progressive drought-stressed conditions. Eighteen QTLs controlling a total of 14 root and shoot traits were identified. A QTL-hotspot genomic region related to a number of root and shoot characteristics associated with drought tolerance such as dry root biomass, root surface area, lateral root number, dry shoot biomass and shoot length was identified. Interestingly, a QTL (QRSratioIX-2.30) related to root-shoot ratio, an important trait for drought avoidance, explaining the highest phenotypic variance of 27.6 and 28.9% for the two consecutive seasons, respectively, was detected. This QTL was closed to the co-dominant SNP marker TP6337 and also flanked by the two SNP TP518 and TP1280. An important QTL (QLRNIII-98.64) related to lateral root number was found close to TP3371 and flanked by TP5093 and TP6072 SNP markers. Also, a QTL (QSRLIV-61.63) associated with specific root length was identified close to TP1873 and flanked by F7XEM6b SRAP marker and TP1035 SNP marker. These two QTLs were detected in both seasons. Our results could be used for marker-assisted selection in lentil breeding programs targeting root and shoot characteristics conferring drought avoidance as an efficient alternative to slow and labor-intensive conventional breeding methods.

Genetic Relationship in Cicer Sp. Expose Evidence for Geneflow between the Cultigen and Its Wild Progenitor

There is a debate concerning mono- or poly-phyletic origins of the Near Eastern crops. In parallel, some authors claim that domestication was not possible within the natural range of the wild progenitors due to wild alleles flow into the nascent crops. Here we address both, the mono- or poly-phyletic origins and the domestications within or without the natural range of the progenitor, debates in order to understand the relationship between domesticated chickpea (Cicer arietinum L.) and its wild progenitor (C. reticulatum Ladizinsky) with special emphasis on its domestication centre in southeastern Turkey. A set of 103 chickpea cultivars and landraces from the major growing regions alongside wild accessions (C. reticulatum, C. echinospermum P.H Davis and C. bijugum K.H. Rech) sampled across the natural distribution range in eastern Turkey were genotyped with 194 SNPs markers. The genetic affinities between and within the studied taxa were assessed. The analysis suggests a mono-phyletic origin of the cultigen, with several wild accession as likely members of the wild stock of the cultigen. Clear separation between the wild and domesticated germplasm was apparent, with negligible level of admixture. A single C. reticulatum accession shows morphological and allelic signatures of admixture, a likely result of introgression. No evidence of geneflow from the wild into domesticated germplasm was found. The traditional farming systems of southeaster Turkey are characterized by occurrence of sympatric wild progenitor-domesticated forms of chickpea (and likewise cereals and other grain legumes). Therefore, both the authentic crop landraces and the wild populations native to the area are a unique genetic resource. Our results grant support to the notion of domestication within the natural distribution range of the wild progenitor, suggesting that the Neolithic domesticators were fully capable of selecting the desired phenotypes even when facing rare wild-domesticated introgression events.

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.

Lentils (Lens culinaris L.), a rich source of folates

The potential for genetic biofortification of U.S.-grown lentils ( Lens culinaris L.) with bioavailable folate has not been widely studied. The objectives of this study were (1) to determine the folate concentration of 10 commercial lentil cultivars grown in Minot and McLean counties, North Dakota, USA, in 2010 and 2011, (2) to determine the genotype (G) × environmental (E) interactions for folate concentration in lentil cultivars, and (3) to compare the folate concentration of other pulses [field peas ( Pisum sativum L.) and chickpea ( Cicer arietinum L.)] grown in the United States. Folate concentration in lentil cultivars ranged from 216 to 290 μg/100 g with a mean of 255 μg/100 g. In addition, lentil showed higher folate concentration compared to chickpea (42-125 μg/100 g), yellow field pea (41-55 μg/100 g), and green field pea (50-202 μg/100 g). A 100 g serving of lentils could provide a significant amount of the recommended daily allowance of dietary folates (54-73%) for adults. A significant year × location interaction on lentil folate concentration was observed; this indicates that possible location sourcing may be required for future lentil folate research.

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.

Mendel, 150 years on

Mendel's paper 'Versuche über Pflanzen-Hybriden' is the best known in a series of studies published in the late 18th and 19th centuries that built our understanding of the mechanism of inheritance. Mendel investigated the segregation of seven gene characters of pea (Pisum sativum), of which four have been identified. Here, we review what is known about the molecular nature of these genes, which encode enzymes (R and Le), a biochemical regulator (I) and a transcription factor (A). The mutations are: a transposon insertion (r), an amino acid insertion (i), a splice variant (a) and a missense mutation (le-1). The nature of the three remaining uncharacterized characters (green versus yellow pods, inflated versus constricted pods, and axial versus terminal flowers) is discussed.

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.

Identification of Mendel's white flower character

Background

The genetic regulation of flower color has been widely studied, notably as a character used by Mendel and his predecessors in the study of inheritance in pea.

Methodology/Principal Findings

We used the genome sequence of model legumes, together with their known synteny to the pea genome to identify candidate genes for the A and A2 loci in pea. We then used a combination of genetic mapping, fast neutron mutant analysis, allelic diversity, transcript quantification and transient expression complementation studies to confirm the identity of the candidates.

Conclusions/Significance

We have identified the pea genes A and A2. A is the factor determining anthocyanin pigmentation in pea that was used by Gregor Mendel 150 years ago in his study of inheritance. The A gene encodes a bHLH transcription factor. The white flowered mutant allele most likely used by Mendel is a simple G to A transition in a splice donor site that leads to a mis-spliced mRNA with a premature stop codon, and we have identified a second rare mutant allele. The A2 gene encodes a WD40 protein that is part of an evolutionarily conserved regulatory complex.

First Report of Bean yellow mosaic virus from Diseased Lupinus luteus in Eastern Washington

Lupine accessions from the Cool Season Food Legume Seed Collection are grown for seed regenerations in Pullman, WA by the Agricultural Research Service, Western Regional Plant Introduction Station. Selected seed was germinated in the greenhouse and assayed by indirect ELISA using antiserum for potyvirus group detection (Agdia, Inc., Elkhart, IN). Healthy transplants were grown for seed collection on outside plots. In July of 2005, more than 90% of 307 Lupinus luteus L. transplants developed severe yellowing, necrosis, and stunting with an estimated 5% plant death. Plants were heavily infested with aphids and leaf sap was serologically positive for potyvirus. Partially purified virus preparations from infected plants contained filamentous particles and a 35-kDa protein that reacted with universal potyvirus antiserum on western blots. Reverse transcription (RT)-PCR using potyvirus universal primers (2) and cDNA derived from virion RNA generated a ~1.7-kbp product that was cloned and sequenced. The sequenced portion of the genomic RNA contained 1,610 nucleotides (nt) on its 3'-terminus (GenBank Accession No. EU144223) that included a partial nuclear inclusion protein, NIb, (1 to 637 nt) with the conserved amino acid (aa) replicase motif GDD (131 to 139 nt), the coat protein (CP) gene of 821 nt (638 to 1,459 nt), and a 171-nt untranslated region (1,460 to 1,630 nt) attached to a poly(A)tail. The CP sequence contained a NAG motif instead of the DAG motif commonly associated with aphid transmission. Searches in the NCBI GenBank database revealed that the CP aa and nt sequences contained conserved domains with isolates of Bean yellow mosaic virus (BYMV). A pairwise alignment (ClustalX) (4) of the CP aa from 20 BYMV isolates with the BYMV-Pullman isolate revealed identities from 96% (BYMV-S, U47033) to 88% (BYMV-MI [X81124)] -MI-NAT [AF434661]). This meets the species demarcation criteria of more than ~80% identity for inclusion with BYMV (1). Virion mechanical inoculations resulted in local lesions on Chenopodium amaranticolor Coste et Reyn and C. quinoa Willd., necrotic blotches on Phaseolus vulgaris L., and yellow spots and systemic movement in L. succulentus Douglas ex. K. Koch, L. texensis 'Bluebonnet', and L. texensis 'Maroon'; BYMV was confirmed by western blots and ELISA. The experimental inoculations represent the first documented report of BYMV in the annual L. succulentus and biennial L. texensis species. Since BYMV is seedborne and transmitted by many aphid species (3), it is possible that several lupine transplants escaped potyvirus detection, and secondary transmission of BYMV to plants occurred by aphids. During the 1950s, BYMV was confirmed in several annual lupines grown as crops in the southeastern United States (3). To our knowledge, this is the first report of BYMV occurring naturally in a lupine species in Washington. BYMV is a destructive virus to lupine species worldwide and has a wide host range in Fabaceae. This research directly contributes toward the maintenance of virus-free lupine seed for distribution to scientists focusing on lupine research. References: (1) P. H. Berger et al. Family Potyviridae. Page 819 in: Virus Taxonomy: Eighth Report of the ICTV. C. M. Fauquet et al. eds., 2005. (2) J. Chen et al. Arch. Virol. 146:757, 2001. (3) R. A. C. Jones and G. D. Mclean, Ann. Appl. Biol. 114:609, 1989. (4) J. D. Thompson et al. Nucleic Acids Res. 24:4878, 1997.

Construction and characterization of two bacterial artificial chromosome libraries of pea (Pisum sativum L.) for the isolation of economically important genes

Pea (Pisum sativum L.) has a genome of about 4 Gb that appears to share conserved synteny with model legumes having genomes of 0.2-0.4 Gb despite extensive intergenic expansion. Pea plant inventory (PI) accession 269818 has been used to introgress genetic diversity into the cultivated germplasm pool. The aim here was to develop pea bacterial artificial chromosome (BAC) libraries that would enable the isolation of genes involved in plant disease resistance or control of economically important traits. The BAC libraries encompassed about 3.2 haploid genome equivalents consisting of partially HindIII-digested DNA fragments with a mean size of 105 kb that were inserted in 1 of 2 vectors. The low-copy oriT-based T-DNA vector (pCLD04541) library contained 55 680 clones. The single-copy oriS-based vector (pIndigoBAC-5) library contained 65 280 clones. Colony hybridization of a universal chloroplast probe indicated that about 1% of clones in the libraries were of chloroplast origin. The presence of about 0.1% empty vectors was inferred by white/blue colony plate counts. The usefulness of the libraries was tested by 2 replicated methods. First, high-density filters were probed with low copy number sequences. Second, BAC plate-pool DNA was used successfully to PCR amplify 7 of 9 published pea resistance gene analogs (RGAs) and several other low copy number pea sequences. Individual BAC clones encoding specific sequences were identified. Therefore, the HindIII BAC libraries of pea, based on germplasm accession PI 269818, will be useful for the isolation of genes underlying disease resistance and other economically important traits.

Consistent Quantitative Trait Loci in Pea for Partial Resistance to Aphanomyces euteiches Isolates from the United States and France

ABSTRACT Development of pea cultivars resistant to Aphanomyces root rot, the most destructive root disease of pea worldwide, is a major disease management objective. In a previous study of a mapping population of 127 recombinant inbred lines (RILs) derived from the cross 'Puget' (susceptible) x '90-2079' (partially resistant), we identified seven genomic regions, including a major quantitative trait locus (QTL), Aph1, associated with partial resistance to Aphanomyces root rot in U.S. fields (21). The objective of the present study was to evaluate, in the same mapping population, the specificity versus consistency of Aphanomyces resistance QTL under two screening conditions (greenhouse and field, by comparison with the previous study) and with two isolates of Aphanomyces euteiches originating from the United States and France. The 127 RILs were evaluated in the greenhouse for resistance to pure culture isolates SP7 (United States) and Ae106 (France). Using the genetic map previously described, a total of 10 QTL were identified for resistance in greenhouse conditions to the two isolates. Among these were Aph1, Aph2, and Aph3, previously detected for partial field resistance in the United States. Aph1 and Aph3 were detected with both isolates and Aph2 with only the French isolate. Seven additional QTL were specifically detected with one of the two isolates and were not identified for partial field resistance in the United States. The consistency of the detected resistance QTL over two screening environments and isolates is discussed with regard to pathogen variability, and disease assessment and QTL detection methods. This study suggests the usefulness of three consistent QTL, Aph1, Aph2, and Aph3, for marker-assisted selection.

Quantitative trait loci for partial resistance to Aphanomyces root rot in pea

Aphanomyces root rot, caused by Aphanomyces euteiches Drechs, is the most-important disease of pea ( Pisum sativum L.) worldwide. No efficient chemicals are available to control the pathogen. To facilitate breeding for Aphanomyces root rot resistance and to better understand the inheritance of partial resistance, our goal was to identify QTLs associated with field partial resistance. A population of 127 RILs from the cross Puget (susceptible) x 90-2079 (partially resistant) was used. The lines were assessed for resistance to A. euteiches under field conditions at two locations in the United States (Pullman, Wash. and LeSueur, Minn.) in 1996 and 1998 for three criteria based on symptom intensity and disease effects on the whole plant. The RILs were genotyped using automated AFLPs, RAPDs, SSRs, ISSRs, STSs, isozymes and morphological markers. The resulting genetic map consisted of 324 linked markers distributed over 13 linkage groups covering 1,094 cM (Kosambi). Twenty seven markers were anchored to other published pea genetic maps. A total of seven genomic regions were associated with Aphanomyces root rot resistance. The first one, located on LG IVb and named Aph1, was considered as "major" since it was highly consistent over the years, locations and resistance criteria studied, and it explained up to 47% of the variation in the 1998 Minnesota trial. Two other year-specific QTLs, namely Aph2 and Aph3, were revealed from different scoring criteria on LG V and Ia, respectively. Aph2 and Aph3 mapped near the r (wrinkled/round seeds) and af (normal/afila leaves) genes, and accounted for up to 32% and 11% of the variation, respectively. Four other "minor" QTLs, identified on LG Ib, VII and B, were specific to one environment and one resistance criterion. The resistance alleles of Aph3 and the two "minor" QTLs on LG Ib were derived from the susceptible parent. Flanking markers for the major Aphanomyces resistance QTL, Aph1, have been identified for use in marker-assisted selection to improve breeding efficiency.