Microsatellite marker polymorphism and mapping in pea (Pisum sativum L.)

Fine mapping and identification of a Fusarium wilt resistance gene FwS1 in pea

KEY MESSAGE

A Fusarium wilt resistance gene FwS1 on pea chromosome 6 was identified and mapped to a 91.4 kb region by a comprehensive genomic-based approach, and the gene Psat6g003960 harboring NB-ARC domain was identified as the putative candidate gene. Pea Fusarium wilt, incited by Fusarium oxysporum f. sp. pisi (Fop), has always been a devastating disease that causes severe yield losses and economic damage in pea-growing regions worldwide. The utilization of pea cultivars carrying resistance gene is the most efficient approach for managing this disease. In order to finely map resistance gene, F2 populations were established through the cross between Shijiadacaiwan 1 (resistant) and Y4 (susceptible). The resistance genetic analysis indicated that the Fop resistance in Shijiadacaiwan 1 was governed by a single dominant gene, named FwS1. Based on the bulked segregant analysis sequencing analyses, the gene FwS1 was initially detected on chromosome 6 (i.e., linking group II, chr6LG2), and subsequent linkage mapping with 589 F2 individuals fine-mapped the gene FwS1 into a 91.4 kb region. The further functional annotation and haplotype analysis confirmed that the gene Psat6g003960, characterized by a NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) domain, was considered as the most promising candidate gene. The encoding amino acids were altered by a "T/C" single-nucleotide polymorphism (SNP) in the first exon of the Psat6g003960, and based on this SNP locus, the molecular marker A016180 was determined to be a diagnostic marker for FwS1 by validating its specificity in both pea accessions and genetic populations with different genetic backgrounds. The FwS1 with diagnostic KASP marker A016180 could facilitate marker-assisted selection in resistance pea breeding in pea. In addition, a comparison of the candidate gene Psat6g003960 in 74SN3B and SJ1 revealed the same sequences. This finding indicated that 74SN3B carried the candidate gene for FwS1, suggesting that FwS1 and Fwf may be closely linked or an identical resistant gene against Fusarium wilt.

The Fastest and Most Reliable Identification of True Hybrids in the Genus Pisum L

After crosses, the identification of true hybrids is not only the most important step in the initiation of a breeding program but also plays a crucial role in the improvement of hybrid varieties. However, current morphological or molecular-based hybrid identification methods are time-consuming and costly approaches that require knowledge and skill, as well as specific lab equipment. In the current study, xenia, direct or immediate effect of pollen on seeds was used to identify true hybrids in the genus Pisum L. for the first time without growing F1 plants. The current study was therefore aimed to (i) elucidate the xenia effect on seeds in intra- and interspecific crosses between P. sativum L. subsp. sativum var. sativum or var. arvense L. Poir. and its wild relatives, including P. sativum subsp. elatius (M. Bieb.) Aschers & Graebn. and P. fulvum Sibth. & Sm., and (ii) illuminate the beneficialness of the xenia effect in a practical improvement of the genus Pisum L. The pea cultivars, including P. sativum subsp. sativum var. sativum and P. sativum subsp. sativum var. arvense, were therefore crossed with P. sativum subsp. elatius and P. fulvum, and the occurrence of the xenia effect was studied on the seeds of fertilized female plants immediately after the crosses. It was concluded that using the xenia effect for the early detection of true hybrid immediately after crossing was not only the fastest, most reliable, and least expensive option as early selection criteria, but that xenia also provided information about dominant seed and pod traits after double fertilization.

High-density linkage mapping and genetic dissection of resistance to broomrape (Orobanche crenata Forsk.) in pea (Pisum sativum L.)

Pea (Pisum sativum L.) is a widely cultivated legume of major importance for global food security and agricultural sustainability. Crenate broomrape (Orobanche crenata Forsk.) (Oc) is a parasitic weed severely affecting legumes, including pea, in the Mediterranean Basin and the Middle East. Previously, the identification of the pea line "ROR12", displaying resistance to Oc, was reported. Two-year field trials on a segregant population of 148 F₇ recombinant inbred lines (RILs), originating from a cross between "ROR12" and the susceptible cultivar "Sprinter", revealed high heritability (0.84) of the "ROR12" resistance source. Genotyping-by-sequencing (GBS) on the same RIL population allowed the construction of a high-density pea linkage map, which was compared with the pea reference genome and used for quantitative trait locus (QTL) mapping. Three QTLs associated with the response to Oc infection, named PsOcr-1, PsOcr-2, and PsOcr-3, were identified, with PsOcr-1 explaining 69.3% of the genotypic variance. Evaluation of the effects of different genotypic combinations indicated additivity between PsOcr-1 and PsOcr-2, and between PsOcr-1 and PsOcr-3, and epistasis between PsOcr-2 and PsOcr-3. Finally, three Kompetitive Allele Specific PCR (KASP) marker assays were designed on the single-nucleotide polymorphisms (SNPs) associated with the QTL significance peaks. Besides contributing to the development of pea genomic resources, this work lays the foundation for the obtainment of pea cultivars resistant to Oc and the identification of genes involved in resistance to parasitic Orobanchaceae.

Five Regions of the Pea Genome Co-Control Partial Resistance to D. pinodes, Tolerance to Frost, and Some Architectural or Phenological Traits

Evidence for reciprocal links between plant responses to biotic or abiotic stresses and architectural and developmental traits has been raised using approaches based on epidemiology, physiology, or genetics. Winter pea has been selected for years for many agronomic traits contributing to yield, taking into account architectural or phenological traits such as height or flowering date. It remains nevertheless particularly susceptible to biotic and abiotic stresses, among which Didymella pinodes and frost are leading examples. The purpose of this study was to identify and resize QTL localizations that control partial resistance to D. pinodes, tolerance to frost, and architectural or phenological traits on pea dense genetic maps, considering how QTL colocalizations may impact future winter pea breeding. QTL analysis revealed five metaQTLs distributed over three linkage groups contributing to both D. pinodes disease severity and frost tolerance. At these loci, the haplotypes of alleles increasing both partial resistance to D. pinodes and frost tolerance also delayed the flowering date, increased the number of branches, and/or decreased the stipule length. These results question both the underlying mechanisms of the joint control of biotic stress resistance, abiotic stress tolerance, and plant architecture and phenology and the methods of marker-assisted selection optimizing stress control and productivity in winter pea breeding.

Genetic Diversity Assessment of Spanish and Some Endangered Tunisian Pea (Pisum sativum L.) Accessions Based on Microsatellite Markers (SSRs)

In the current investigation, 28 accessions of Spanish and Tunisian peas were characterized by eight SSR polymorphic markers to assess their genetic diversity. Many methods have been applied to evaluate these relationships including diversity indices, analysis of molecular variance, cluster analysis, and population structure. The means of diversity indices, the polymorphism information content (PIC), the allelic richness, and the Shannon information index were 0.51, 3.87, and 0.9, respectively. These results revealed a large polymorphism (84.15 %) which produced a higher degree of genetic distance amongst the accessions. The unweighted pair group approach with arithmetic mean divided the collection of these accessions into three major genetic clusters. Therefore, this article has clearly demonstrated the usefulness of the SSR markers that can significantly contribute to the management and conservation of pea germplasm in these countries, as well as to future reproduction.

In vivo and in vitro validation of powdery mildew resistance in garden pea genotypes

Powdery mildew is one of the serious diseases of garden pea which causes a large number of yield losses. Genetic resistance is quite effective, being cost-effective and environment friendly than fungicide applications. In the present studies an initial attempt has been made to identify resistant genotypes against powdery mildew disease developed from hybridization followed by validation of the disease. The experimental material comprised of 48 genotypes that includes 44 advanced breeding lines was evaluated for powdery mildew incidence in Randomized Complete Block Design with three replications at two locations under field conditions [Palampur (winter 2017-18 and 2018-19) and Kukumseri (summer 2018)] and in vitro at Palampur [detached leaf method and polyhouse conditions]. Ten lines viz., SP7, SN-1, SN-6-1, SN-7-1, SN-2, SN-5-2, SN-6-2, SN-10, SN-21 and SP-281 showed resistant reaction along with check Palam Sumool while 27 lines were identified as moderately resistant in comparison to susceptible check Azad P-1. Besides, six lines namely, SP-2, SP-5, SP-10, SP-24, SA-4 and SP-12-1 gave moderately susceptible reaction along with checks Pb-89 and Palam Priya. Only, SP-19 was categorized as susceptible. The high yielding lines SP-3, SP-6 and SP-22 showed moderately resistant reaction in both natural and artificial conditions. Validation of resistance using molecular markers revealed that neither the parental genotypes nor the progenies possess the er1 gene of JI1559. The er2 linked marker ScOPX-171700 was polymorphic between Palam Sumool and Palam Priya but the marker didn't show polymorphism between er2 harboring line (JI2480). These results suggested that the lines showing resistance under field conditions may have some other genes or alleles for resistance and further confirmation is needed by developing mapping populations with specific gene or gene combinations.

Rust (Uromyces viciae-fabae Pers. de-Bary) of Pea (Pisum sativum L.): Present Status and Future Resistance Breeding Opportunities

Uromyces viciae-fabae Pers. de-Bary is an important fungal pathogen causing rust in peas (Pisum sativum L.). It is reported in mild to severe forms from different parts of the world where the pea is grown. Host specificity has been indicated in this pathogen in the field but has not yet been established under controlled conditions. The uredinial states of U. viciae-fabae are infective under temperate and tropical conditions. Aeciospores are infective in the Indian subcontinent. The genetics of rust resistance was reported qualitatively. However, non-hypersensitive resistance responses and more recent studies emphasized the quantitative nature of pea rust resistance. Partial resistance/slow rusting had been described as a durable resistance in peas. Such resistance is of the pre-haustorial type and expressed as longer incubation and latent period, poor infection efficiency, a smaller number of aecial cups/pustules, and lower units of AUDPC (Area Under Disease Progress Curve). Screening techniques dealing with slow rusting should consider growth stages and environment, as both have a significant influence on the disease scores. Our knowledge about the genetics of rust resistance is increasing, and now molecular markers linked with gene/QTLs (Quantitative Trait Loci) of rust resistance have been identified in peas. The mapping efforts conducted in peas came out with some potent markers associated with rust resistance, but they must be validated under multi-location trails before use in the marker-assisted selection of rust resistance in pea breeding programs.

Mapping QTLs for Super-Earliness and Agro-Morphological Traits in RILs Population Derived from Interspecific Crosses between Pisum sativum × P. fulvum

Earliness in crop plants has a vital role in prevention of heat-induced drought stress and in combating global warming, which is predicted to exacerbate in the near future. Furthermore, earliness may expand production into northern areas or higher altitudes, having relatively shorter growing season and may also expand arable lands to meet global food demands. The primary objective of the present study was to investigate quantitative trait loci (QTLs) for super-earliness and important agro-morphological traits in a recombinant inbred line (RIL) population derived from an interspecific cross. A population of 114 RILs developed through single-seed descent from an interspecific cross involving Pisum sativum L. and P. fulvum Sibth. et Sm. was evaluated to identify QTLs for super-earliness and important agro-morphological traits. A genetic map was constructed with 44 SSRs markers representing seven chromosomes with a total length of 262.6 cM. Of the 14 QTLs identified, two were for super-earliness on LG2, one for plant height on LG3, six for number of pods per plant on LG2, LG4, LG5 and LG6, one for number of seeds per pod on LG6, one for pod length on LG4 and three for harvest index on LG3, LG5, and LG6. AA205 and AA372-1 flanking markers for super-earliness QTLs were suggested for marker-assisted selection (MAS) in pea breeding programs due to high heritability of the trait. This is the first study to map QTLs originating from P. sativum and P. fulvum recently identified species with super-earliness character and the markers (AA205 and AA372-1) linked to QTLs were valuable molecular tools for pea breeding.

Molecular Characterizations of the er1 Alleles Conferring Resistance to Erysiphe pisi in Three Chinese Pea (Pisum sativum L.) Landraces

Powdery mildew caused by Erysiphe pisi DC. is a major disease affecting pea worldwide. This study aimed to confirm the resistance genes contained in three powdery mildew-resistant Chinese pea landraces (Suoshadabaiwan, Dabaiwandou, and Guiwan 1) and to develop the functional markers of the novel resistance genes. The resistance genes were identified by genetic mapping and PsMLO1 gene sequence identification. To confirm the inheritance of powdery mildew resistance in the three Landraces, the susceptible cultivars Bawan 6, Longwan 1, and Chengwan 8 were crossed with Suoshadabaiwan, Dabaiwandou, and Guiwan 1 to produce F1, F2, and F2:3 populations, respectively. All F1 plants were susceptible to E. pisi, and phenotypic segregation patterns in all the F2 and F2:3 populations fit the 3:1 (susceptible: resistant) and 1:2:1 (susceptible homozygotes: heterozygotes: resistant homozygotes) ratios, respectively, indicating powdery mildew resistance in the three Landraces were controlled by a single recessive gene, respectively. The analysis of er1-linked markers and genetic mapping in the F2 populations suggested that the recessive resistance genes in three landraces could be er1 alleles. The cDNA sequences of 10 homologous PsMLO1 cDNA clones from the contrasting parents were obtained. A known er1 allele, er1-4, was identified in Suoshadabaiwan. Two novel er1 alleles were identified in Dabaiwandou and Guiwan 1, which were designated as er1-13 and er1-14, respectively. Both novel alleles were characterized with a 1-bp deletion (T) in positions 32 (exon 1) and 277 (exon 3), respectively, which caused a frame-shift mutation to result in premature termination of translation of PsMLO1 protein. The co-dominant functional markers specific for er1-13 and er1-14, KASPar-er1-13, and KASPar-er1-14 were developed and effectively validated in populations and pea germplasms. Here, two novel er1 alleles were characterized and their functional markers were validated. These results provide powerful tools for marker-assisted selection in pea breeding.

SSR-Based Molecular Identification and Population Structure Analysis for Forage Pea (Pisum sativum var. arvense L.) Landraces

Plant genetic diversity has a significant role in providing traits that can help meet future challenges, such as the need to adapt crops to changing climatic conditions or outbreaks of disease. Our aim in this study was to evaluate the diversity of 61 forage pea specimens (P. sativum ssp. arvense L.) collected from the northeastern Anatolia region of Turkey using 28 simple sequence repeat (SSR) markers. These primers generated a total of 82 polymorphic bands. The number of observed alleles (Na) per primer varied from 2 to 4 with a mean of 2.89 alleles/locus. The mean value of expected heterozygosity (Exp-Het = 0.50) was higher than the mean value of observed heterozygosity (Obs-Het = 0.22). The mean of polymorphic information content (PIC) was 0.41 with a range of 0.03-0.70. The mean number of effective alleles (Ne) was found to be 2.15, Nei's expected heterozygosity (H) 0.49, and Shannon's information index (I) 0.81. Cluster analysis through the unweighted pair-group mean average (UPGMA) method revealed that 61 forage pea landraces were divided into three main clusters. Genetic dissimilarity between the genotypes, calculated with the use of NTSYS-pc software, varied between 0.10 (G30 and G34) and 0.66 (G1 and G32). Principal coordinate analysis (PCoA) revealed that three principal coordinates explained 51.54% of the total variation. Moreover, population structure analysis showed that all genotypes formed three sub-populations. Expected heterozygosity values varied between 0.2669 (the first sub-population) and 0.3223 (third sub-population), with an average value of 0.2924. Average population differentiation measurement (Fst) was identified as 0.2351 for the first sub-population, 0.3838 for the second sub-population, and 0.2506 for the third sub-population. In general, current results suggest that SSR markers could be constantly used to illuminate the genetic diversity of forage pea landraces and can potentially be incorporated into future studies that examine the diversity within a larger collection of forage pea genotypes from diverse regions.

Genomics Enabled Breeding Strategies for Major Biotic Stresses in Pea (Pisum sativum L.)

Pea (Pisum sativum L.) is one of the most important and productive cool season pulse crops grown throughout the world. Biotic stresses are the crucial constraints in harnessing the potential productivity of pea and warrant dedicated research and developmental efforts to utilize omics resources and advanced breeding techniques to assist rapid and timely development of high-yielding multiple stress-tolerant-resistant varieties. Recently, the pea researcher's community has made notable achievements in conventional and molecular breeding to accelerate its genetic gain. Several quantitative trait loci (QTLs) or markers associated with genes controlling resistance for fusarium wilt, fusarium root rot, powdery mildew, ascochyta blight, rust, common root rot, broomrape, pea enation, and pea seed borne mosaic virus are available for the marker-assisted breeding. The advanced genomic tools such as the availability of comprehensive genetic maps and linked reliable DNA markers hold great promise toward the introgression of resistance genes from different sources to speed up the genetic gain in pea. This review provides a brief account of the achievements made in the recent past regarding genetic and genomic resources' development, inheritance of genes controlling various biotic stress responses and genes controlling pathogenesis in disease causing organisms, genes/QTLs mapping, and transcriptomic and proteomic advances. Moreover, the emerging new breeding approaches such as transgenics, genome editing, genomic selection, epigenetic breeding, and speed breeding hold great promise to transform pea breeding. Overall, the judicious amalgamation of conventional and modern omics-enabled breeding strategies will augment the genetic gain and could hasten the development of biotic stress-resistant cultivars to sustain pea production under changing climate. The present review encompasses at one platform the research accomplishment made so far in pea improvement with respect to major biotic stresses and the way forward to enhance pea productivity through advanced genomic tools and technologies.

Integrating Phenotypic and Gene Expression Linkage Mapping to Dissect Rust Resistance in Chickling Pea

Rusts are among the most important foliar biotrophic fungal diseases in legumes. Lathyrus cicera crop can be severely damaged by Uromyces pisi, to which partial resistance has been identified. Nevertheless, the underlying genetic basis and molecular mechanisms of this resistance are poorly understood in L. cicera. To prioritise the causative variants controlling partial resistance to rust in L. cicera, a recombinant inbred line (RIL) population, segregating for response to this pathogen, was used to combine the detection of related phenotypic- and expression-quantitative trait loci (pQTLs and eQTLs, respectively). RILs' U. pisi disease severity (DS) was recorded in three independent screenings at seedling (growth chamber) and in one season of exploratory screening at adult plant stage (semi-controlled field conditions). A continuous DS range was observed in both conditions and used for pQTL mapping. Different pQTLs were identified under the growth chamber and semi-controlled field conditions, indicating a distinct genetic basis depending on the plant developmental stage and/or the environment. Additionally, the expression of nine genes related to U. pisi resistance in L. cicera was quantified for each RIL individual and used for eQTL mapping. One cis-eQTL and one trans-eQTL were identified controlling the expression variation of one gene related to rust resistance - a member of glycosyl hydrolase family 17. Integrating phenotyping, gene expression and linkage mapping allowed prioritising four candidate genes relevant for disease-resistance precision breeding involved in adaptation to biotic stress, cellular, and organelle homeostasis, and proteins directly involved in plant defence.

Gene-Based Resistance to Erysiphe Species Causing Powdery Mildew Disease in Peas (Pisum sativum L.)

Globally powdery mildew (PM) is one of the major diseases of the pea caused by Erysiphe pisi. Besides, two other species viz. Erysiphe trifolii and Erysiphe baeumleri have also been identified to infect the pea plant. To date, three resistant genes, namely er1, er2 and Er3 located on linkage groups VI, III and IV respectively were identified. Studies have shown the er1 gene to be a Pisum sativum Mildew resistance Locus 'O' homologue and subsequent analysis has identified eleven alleles namely er1-1 to er1-11. Despite reports mentioning the breakdown of er1 gene-mediated PM resistance by E. pisi and E. trifolii, it is still the most widely deployed gene in PM resistance breeding programmes across the world. Several linked DNA markers have been reported in different mapping populations with varying linkage distances and effectiveness, which were used by breeders to develop PM-resistant pea cultivars through marker assisted selection. This review summarizes the genetics of PM resistance and its mechanism, allelic variations of the er gene, marker linkage and future strategies to exploit this information for targeted PM resistance breeding in Pisum.

Identification of Quantitative Trait Loci Associated With Partial Resistance to Fusarium Root Rot and Wilt Caused by Fusarium graminearum in Field Pea

Fusarium root rot, caused by a complex of Fusarium spp., is a major disease of field pea (Pisum sativum). The development of genetic resistance is the most promising approach to manage the disease, but no pea germplasm has been identified that is completely resistant to root rot. The aim of this study was to detect quantitative trait loci (QTL) conferring partial resistance to root rot and wilting, caused by five fungal isolates representing Fusarium solani, F. avenaceum, F. acuminatum, F. proliferatum, and F. graminearum. Evaluation of the root rot-tolerant cultivar "00-2067" and susceptible cultivar "Reward" was carried out with the five species. There was a significant difference (p < 0.001) between the mean root rot values of the two cultivars inoculated with the F. avenaceum (F4A) and F. graminearum (FG2) isolates. Therefore, in the QTL study, the F8 recombinant inbred line (RIL) population derived from "Reward" × "00-2067" was inoculated in the greenhouse (4 ×) with only F4A and FG2. The parents and F8 population were genotyped using 13.2K single nucleotide polymorphisms (SNPs) and 222 simple sequence repeat (SSR) markers. A significant genotypic effect (p < 0.05) and high heritability (79% to 92.1%) were observed for disease severity, vigor, and plant height following inoculation with F4A and FG2. Significant correlation coefficients were detected among and within all traits. This suggested that a high proportion of the genetic variance was transmitted from the parents to the progeny. However, no significant QTL (LOD > 3) were detected for the RILs inoculated with F4A. In the case of the RILs inoculated with FG2, 5 QTL for root rot severity and 3 QTL each for vigor and plant height were detected. The most stable QTL for plant height (Hgt-Ps3.1) was detected on Chrom5/LGIII. The two most stable QTL for partial resistance to FG2, Fg-Ps4.1, and Fg-Ps4.2 were located in a 15.1-cM and 11.2-cM genomic region, respectively, on Chrom4/LGIV. The most stable QTL for vigor (Vig-Ps4.1) was found in the same region. Twenty-five major and moderate effect digenic epistatic interactions were detected. The identified region on chrom4/LGIV could be important for resistance breeding and marker development.

A Non-Rogue Mutant Line Induced by ENU Mutagenesis in Paramutated Rogue Peas (Pisum sativum L.) Is Still Sensitive to the Rogue Paramutation

The spontaneously emerging rogue phenotype in peas (Pisum sativum L.), characterized by narrow and pointed leaf stipula and leaflets, was the first identified case of the epigenetic phenomenon paramutation. The crosses of homozygous or heterozygous (e.g., F1) rogue plants with non-rogue (wild type) plants, produce exclusively rogue plants in the first and all subsequent generations. The fact that the wild phenotype disappears forever, is in clear contradiction with the Mendelian rules of inheritance, a situation that impedes the positional cloning of genes involved in this epigenetic phenomenon. One way of overcoming this obstacle is the identification of plant genotypes harboring naturally occurring or artificially induced neutral alleles, non-sensitive to paramutation. So far, such alleles have never been described for the pea rogue paramutation. Here, we report the induction via 1-ethyl-1-nitrosourea (ENU) mutagenesis of a non-rogue revertant mutant in the rogue cv. Progreta, and the completely unusual fixation of the induced non-rogue phenotype through several generations. The reversion of the methylation status of two previously identified differentially methylated genomic sequences in the induced non-rogue mutant, confirms that the rogue paramutation is accompanied by alterations in DNA methylation. Nevertheless, unexpectedly, the induced non-rogue mutant showed to be still sensitive to paramutation.

Mapping QTL associated with partial resistance to Aphanomyces root rot in pea (Pisum sativum L.) using a 13.2 K SNP array and SSR markers

KEY MESSAGE

A stable and major QTL, which mapped to an approximately 20.0 cM region on pea chromosome 4, was identified as the most consistent region conferring partial resistance to Aphanomyces euteiches. Aphanomyces root rot (ARR), caused by Aphanomyces euteiches Drechs., is a destructive soilborne disease of field pea (Pisum Sativum L.). No completely resistant pea germplasm is available, and current ARR management strategies rely on partial resistance and fungicidal seed treatments. In this study, an F8 recombinant inbred line population of 135 individuals from the cross 'Reward' (susceptible) × '00-2067' (tolerant) was evaluated for reaction to ARR under greenhouse conditions with the A. euteiches isolate Ae-MDCR1 and over 2 years in a field nursery in Morden, Manitoba. Root rot severity, foliar weight, plant vigor and height were used as estimates of tolerance to ARR. Genotyping was conducted with a 13.2 K single-nucleotide polymorphism (SNP) array and 222 simple sequence repeat (SSR) markers. Statistical analyses of the phenotypic data indicated significant (P < 0.001) genotypic effects and significant G × E interactions (P < 0.05) in all experiments. After filtering, 3050 (23.1%) of the SNP and 30 (13.5%) of the SSR markers were retained for linkage analysis, which distributed 2999 (2978 SNP + 21 SSR) of the markers onto nine linkage groups representing the seven chromosomes of pea. Mapping of quantitative trait loci (QTL) identified 8 major-effect (R2 > 20%), 13 moderate-effect (10% < R2 < 20%) effect and 6 minor-effect (R2 < 10%) QTL. A genomic region on chromosome 4, delimited by the SNP markers PsCam037549_22628_1642 and PsCam026054_14999_2864, was identified as the most consistent region responsible for partial resistance to A. euteiches isolate Ae-MDCR1. Other genomic regions important for resistance were of the order chromosome 5, 6 and 7.

Omics resources and omics-enabled approaches for achieving high productivity and improved quality in pea (Pisum sativum L.)

Pea (Pisum sativum L.), a cool-season legume crop grown in more than 85 countries, is the second most important grain legume and one of the major green vegetables in the world. While pea was historically studied as the genetic model leading to the discovery of the laws of genetics, pea research has lagged behind that of other major legumes in the genomics era, due to its large and complex genome. The evolving climate change and growing population have posed grand challenges to the objective of feeding the world, making it essential to invest research efforts to develop multi-omics resources and advanced breeding tools to support fast and continuous development of improved pea varieties. Recently, the pea researchers have achieved key milestones in omics and molecular breeding. The present review provides an overview of the recent important progress including the development of genetic resource databases, high-throughput genotyping assays, reference genome, genes/QTLs responsible for important traits, transcriptomic, proteomic, and phenomic atlases of various tissues under different conditions. These multi-faceted resources have enabled the successful implementation of various markers for monitoring early-generation populations as in marker-assisted backcrossing breeding programs. The emerging new breeding approaches such as CRISPR, speed breeding, and genomic selection are starting to change the paradigm of pea breeding. Collectively, the rich omics resources and omics-enable breeding approaches will enhance genetic gain in pea breeding and accelerate the release of novel pea varieties to meet the elevating demands on productivity and quality.

Breeding, Genetics, and Genomics Approaches for Improving Fusarium Wilt Resistance in Major Grain Legumes

Fusarium wilt (FW) disease is the key constraint to grain legume production worldwide. The projected climate change is likely to exacerbate the current scenario. Of the various plant protection measures, genetic improvement of the disease resistance of crop cultivars remains the most economic, straightforward and environmental-friendly option to mitigate the risk. We begin with a brief recap of the classical genetic efforts that provided first insights into the genetic determinants controlling plant response to different races of FW pathogen in grain legumes. Subsequent technological breakthroughs like sequencing technologies have enhanced our understanding of the genetic basis of both plant resistance and pathogenicity. We present noteworthy examples of targeted improvement of plant resistance using genomics-assisted approaches. In parallel, modern functional genomic tools like RNA-seq are playing a greater role in illuminating the various aspects of plant-pathogen interaction. Further, proteomics and metabolomics have also been leveraged in recent years to reveal molecular players and various signaling pathways and complex networks participating in host-pathogen interaction. Finally, we present a perspective on the challenges and limitations of high-throughput phenotyping and emerging breeding approaches to expeditiously develop FW-resistant cultivars under the changing climate.

Identification of potential candidate genes controlling pea aphid tolerance in a Pisum fulvum high-density integrated DArTseq SNP-based genetic map

BACKGROUND

Pea (Pisum sativum) is one of the most important temperate grain legumes in the world, and its production is severely constrained by the pea aphid (Acyrthosiphon pisum). Wild relatives, such as P. fulvum, are valuable sources of allelic diversity to improve the genetic resistance of cultivated pea species against A. pisum attack. To unravel the genetic control underlying resistance to the pea aphid attack, a quantitative trait loci (QTL) analysis was performed using the previously developed high density integrated genetic linkage map originated from an intraspecific recombinant inbred line (RIL) population (P. fulvum: IFPI3260 × IFPI3251).

RESULTS

We accurately evaluated specific resistance responses to pea aphid that allowed the identification, for the first time, of genomic regions that control plant damage and aphid reproduction. Eight QTLs associated with tolerance to pea aphid were identified in LGs I, II, III, IV and V, which individually explained from 17.0% to 51.2% of the phenotypic variation depending on the trait scored, and as a whole from 17.0% to 88.6%. The high density integrated genetic linkage map also allowed the identification of potential candidate genes co-located with the QTLs identified.

CONCLUSIONS

Our work shows how the survival of P. fulvum after the pea aphid attack depends on the triggering of a multi-component protection strategy that implies a quantitative tolerance. The genomic regions associated with the tolerance responses of P. fulvum during A. pisum infestation have provided six potential candidate genes that could be useful in marker-assisted selection (MAS) and genomic assisted breeding (GAB) after functional validation in the future. © 2019 Society of Chemical Industry.

Identification of quantitative trait loci (QTL) controlling resistance to pea weevil (Bruchus pisorum) in a high-density integrated DArTseq SNP-based genetic map of pea

Pea weevil (Bruchus pisorum) is a damaging insect pest affecting pea (Pisum sativum) production worldwide. No resistant cultivars are available, although some levels of incomplete resistance have been identified in Pisum germplasm. To decipher the genetic control underlying the resistance previously identify in P. sativum ssp. syriacum, a recombinant inbred line (RIL F8:9) population was developed. The RIL was genotyped through Diversity Arrays Technology PL's DArTseq platform and screened under field conditions for weevil seed infestation and larval development along 5 environments. A newly integrated genetic linkage map was generated with a subset of 6,540 markers, assembled into seven linkage groups, equivalent to the number of haploid pea chromosomes. An accumulated distance of 2,503 cM was covered with an average density of 2.61 markers cM-1. The linkage map allowed the identification of three QTLs associated to reduced seed infestation along LGs I, II and IV. In addition, a QTL for reduced larval development was also identified in LGIV. Expression of these QTLs varied with the environment, being particularly interesting QTL BpSI.III that was detected in most of the environments studied. This high-saturated pea genetic map has also allowed the identification of seven potential candidate genes co-located with QTLs for marker-assisted selection, providing an opportunity for breeders to generate effective and sustainable strategies for weevil control.

Two Novel er1 Alleles Conferring Powdery Mildew (Erysiphe pisi) Resistance Identified in a Worldwide Collection of Pea (Pisum sativum L.) Germplasms

Powdery mildew caused by Erysiphe pisi DC. severely affects pea crops worldwide. The use of resistant cultivars containing the er1 gene is the most effective way to control this disease. The objectives of this study were to reveal er1 alleles contained in 55 E. pisi-resistant pea germplasms and to develop the functional markers of novel alleles. Sequences of 10 homologous PsMLO1 cDNA clones from each germplasm accession were used to determine their er1 alleles. The frame shift mutations and various alternative splicing patterns were observed during transcription of the er1 gene. Two novel er1 alleles, er1-8 and er1-9, were discovered in the germplasm accessions G0004839 and G0004400, respectively, and four known er1 alleles were identified in 53 other accessions. One mutation in G0004839 was characterized by a 3-bp (GTG) deletion of the wild-type PsMLO1 cDNA, resulting in a missing valine at position 447 of the PsMLO1 protein sequence. Another mutation in G0004400 was caused by a 1-bp (T) deletion of the wild-type PsMLO1 cDNA sequence, resulting in a serine to leucine change of the PsMLO1 protein sequence. The er1-8 and er1-9 alleles were verified using resistance inheritance analysis and genetic mapping with respectively derived F2 and F2:3 populations. Finally, co-dominant functional markers specific to er1-8 and er1-9 were developed and validated in populations and pea germplasms. These results improve our understanding of E. pisi resistance in pea germplasms worldwide and provide powerful tools for marker-assisted selection in pea breeding.

Mapping-by-sequencing using NGS-based 3'-MACE-Seq reveals a new mutant allele of the essential nodulation gene Sym33 (IPD3) in pea (Pisum sativum L.)

Large collections of pea symbiotic mutants were accumulated in the 1990s, but the causal genes for a large portion of the mutations are still not identified due to the complexity of the task. We applied a Mapping-by-Sequencing approach including Bulk Segregant Analysis and Massive Analysis of cDNA Ends (MACE-Seq) sequencing technology for genetic mapping the Sym11 gene of pea which controls the formation of symbioses with both nodule bacteria and arbuscular-mycorrhizal fungi. For mapping we developed an F 2-population from the cross between pea line N24 carrying the mutant allele of sym11 and the wild type NGB1238 (=JI0073) line. Sequencing libraries were prepared from bulks of 20 plants with mutant and 12 with wild-type phenotype. MACE-Seq differential gene expression analysis between mutant-phenotype and wild-type-phenotype bulks revealed 2,235 genes, of which 514 (23%) were up-regulated and 1,721 (77%) were down-regulated in plant roots inoculated with rhizobia as a consequence of sym11 mutation. MACE-Seq also detected single nucleotide variants between bulks in 217 pea genes. Using a novel mathematical model we calculated the recombination frequency (RF) between the Sym11 gene and these 217 polymorphic genes. Six genes with the lowest RF were converted into CAPS or dCAPS markers and genetically mapped on the complete mapping population of 108 F 2-plants which confirmed their tight linkage to Sym11 and to each other. The Medicago truncatula Gaertn. (Mt) homologs of these genes are located in a distinct region of Mt chromosome 5, which corresponds to linkage group I of pea. Among 94 candidate genes from this region only one was down-regulated-the pea Sym33 homolog of the Mt IPD3 gene which is essential for nodulation. Sequencing of the Sym33 allele of the N24 (sym11) mutant revealed a single nucleotide deletion (c.C319del) in its third exon resulting in a codon shift in the open reading frame and premature translation termination. Thus, we identified a novel mutant allele sym33-4 most probably responsible for the mutant phenotype of the N24 (sym11) line, thereby demonstrating that mapping by MACE-Seq can be successfully used for genetic mapping of mutations and identification of candidate genes in pea.

Analysis of Genetic Variation and Enhancement of Salt Tolerance in French Pea (Pisum Sativum L.)

Pisum sativum L. (field pea) is a crop of a high nutritional value and seed oil content. The characterization of pea germplasm is important to improve yield and quality. This study aimed at using fatty acid profiling and amplified fragment length polymorphism (AFLP) markers to evaluate the variation and relationships of 25 accessions of French pea. It also aimed to conduct a marker-trait associations analysis using the crude oil content as the target trait for this analysis, and to investigate whether 5-aminolevulinic acid (ALA) could enhance salt tolerance in the pea germplasm. The percentage of crude oil of the 25 pea genotypes varied from 2.6 to 3.5%, with a mean of 3.04%. Major fatty acids in all of the accessions were linoleic acid. Moreover, the 12 AFLP markers used were polymorphic. The cluster analysis based on fatty acids data or AFLP data divided the 25 pea germplasm into two main clusters. The gene diversity of the AFLP markers varied from 0.21 to 0.58, with a mean of 0.41. Polymorphic information content (PIC) of pea germplasm varied from 0.184 to 0.416 with a mean of 0.321, and their expected heterozygosity (He) varied from 0.212 to 0.477 with a mean of 0.362. The AFLP results revealed that the Nain Ordinaire cultivar has the highest level of genetic variability, whereas Elatius 3 has the lowest level. Three AFLP markers (E-AAC/M-CAA, E-AAC/M-CAC, and E-ACA/M-CAG) were significantly associated with the crude oil content trait. The response of the Nain Ordinaire and Elatius 3 cultivars to high salinity stress was studied. High salinity (150 mM NaCl) slightly reduced the photosynthetic pigments contents in Nain Ordinaire leaves at a non-significant level, however, the pigments contents in the Elatius 3 leaves were significantly reduced by high salinity. Antioxidant enzymes (APX-ascorbate peroxidase; CAT-catalase; and POD-peroxidase) activities were significantly induced in the Nain Ordinaire cultivar, but non-significantly induced in Elatius 3 by high salinity. Priming the salt-stressed Nain Ordinaire and Elatius 3 plants with ALA significantly enhanced the pigments biosynthesis, antioxidant enzymes activities, and stress-related genes expression, as compared to the plants stressed with salt alone. In conclusion, this study is amongst the first investigations that conducted marker-trait associations in pea, and revealed a sort of correlation between the diversity level and salt tolerance.

Genotyping-by-sequencing application on diploid rose and a resulting high-density SNP-based consensus map

Roses, which have been cultivated for at least 5000 years, are one of the most important ornamental crops in the world. Because of the interspecific nature and high heterozygosity in commercial roses, the genetic resources available for rose are limited. To effectively identify markers associated with QTL controlling important traits, such as disease resistance, abundant markers along the genome and careful phenotyping are required. Utilizing genotyping by sequencing technology and the strawberry genome (Fragaria vesca v2.0.a1) as a reference, we generated thousands of informative single nucleotide polymorphism (SNP) markers. These SNPs along with known bridge simple sequence repeat (SSR) markers allowed us to create the first high-density integrated consensus map for diploid roses. Individual maps were first created for populations J06-20-14-3×"Little Chief" (J14-3×LC), J06-20-14-3×"Vineyard Song" (J14-3×VS) and "Old Blush"×"Red Fairy" (OB×RF) and these maps were linked with 824 SNPs and 13 SSR bridge markers. The anchor SSR markers were used to determine the numbering of the rose linkage groups. The diploid consensus map has seven linkage groups (LGs), a total length of 892.2 cM, and an average distance of 0.25 cM between 3527 markers. By combining three individual populations, the marker density and the reliability of the marker order in the consensus map was improved over a single population map. Extensive synteny between the strawberry and diploid rose genomes was observed. This consensus map will serve as the tool for the discovery of marker-trait associations in rose breeding using pedigree-based analysis. The high level of conservation observed between the strawberry and rose genomes will help further comparative studies within the Rosaceae family and may aid in the identification of candidate genes within QTL regions.

A High-Density Integrated DArTseq SNP-Based Genetic Map of Pisum fulvum and Identification of QTLs Controlling Rust Resistance

Pisum fulvum, a wild relative of pea is an important source of allelic diversity to improve the genetic resistance of cultivated species against fungal diseases of economic importance like the pea rust caused by Uromyces pisi. To unravel the genetic control underlying resistance to this fungal disease, a recombinant inbred line (RIL) population was generated from a cross between two P. fulvum accessions, IFPI3260 and IFPI3251, and genotyped using Diversity Arrays Technology. A total of 9,569 high-quality DArT-Seq and 8,514 SNPs markers were generated. Finally, a total of 12,058 markers were assembled into seven linkage groups, equivalent to the number of haploid chromosomes of P. fulvum and P. sativum. The newly constructed integrated genetic linkage map of P. fulvum covered an accumulated distance of 1,877.45 cM, an average density of 1.19 markers cM^-1 and an average distance between adjacent markers of 1.85 cM. The composite interval mapping revealed three QTLs distributed over two linkage groups that were associated with the percentage of rust disease severity (DS%). QTLs UpDSII and UpDSIV were located in the LGs II and IV respectively and were consistently identified both in adult plants over 3 years at the field (Córdoba, Spain) and in seedling plants under controlled conditions. Whenever they were detected, their contribution to the total phenotypic variance varied between 19.8 and 29.2. A third QTL (UpDSIV.2) was also located in the LGIVand was environmentally specific as was only detected for DS % in seedlings under controlled conditions. It accounted more than 14% of the phenotypic variation studied. Taking together the data obtained in the study, it could be concluded that the expression of resistance to fungal diseases in P. fulvum originates from the resistant parent IFPI3260.

High-Throughput RNA-Seq Data Analysis of the Single Nucleotide Polymorphisms (SNPs) and Zygomorphic Flower Development in Pea (Pisum sativum L.)

Pea (Pisum sativum L.) is a model plant that has been used in classical genetics and organ development studies. However, its large and complex genome has hindered research investigations in pea. Here, we generated transcriptomes from different tissues or organs of three pea accessions using next-generation sequencing to assess single nucleotide polymorphisms (SNPs), and further investigated petal differentially expressed genes to elucidate the mechanisms regulating floral zygomorphy. Eighteen samples were sequenced, which yielded a total of 617 million clean reads, and de novo assembly resulted in 87,137 unigenes. A total of 9044 high-quality SNPs were obtained among the three accessions, and a consensus map was constructed. We further discovered several dorsoventral asymmetrically expressed genes that were confirmed by qRT-PCR among different petals, including previously reported three CYC-like proliferating cell factor (TCP) genes. One MADS-box gene was highly expressed in dorsal petals, and several MYB factors were predominantly expressed among dorsal, lateral, and/or ventral petals, together with a ventrally expressed TCP gene. In sum, our comprehensive database complements the existing resources for comparative genetic mapping and facilitates future investigations in legume zygomorphic flower development.

Pea Marker Database (PMD) - A new online database combining known pea (Pisum sativum L.) gene-based markers

Pea (Pisum sativum L.) is the oldest model object of plant genetics and one of the most agriculturally important legumes in the world. Since the pea genome has not been sequenced yet, identification of genes responsible for mutant phenotypes or desirable agricultural traits is usually performed via genetic mapping followed by candidate gene search. Such mapping is best carried out using gene-based molecular markers, as it opens the possibility for exploiting genome synteny between pea and its close relative Medicago truncatula Gaertn., possessing sequenced and annotated genome. In the last 5 years, a large number of pea gene-based molecular markers have been designed and mapped owing to the rapid evolution of "next-generation sequencing" technologies. However, the access to the complete set of markers designed worldwide is limited because the data are not uniformed and therefore hard to use. The Pea Marker Database was designed to combine the information about pea markers in a form of user-friendly and practical online tool. Version 1 (PMD1) comprises information about 2484 genic markers, including their locations in linkage groups, the sequences of corresponding pea transcripts and the names of related genes in M. truncatula. Version 2 (PMD2) is an updated version comprising 15944 pea markers in the same format with several advanced features. To test the performance of the PMD, fine mapping of pea symbiotic genes Sym13 and Sym27 in linkage groups VII and V, respectively, was carried out. The results of mapping allowed us to propose the Sen1 gene (a homologue of SEN1 gene of Lotus japonicus (Regel) K. Larsen) as the best candidate gene for Sym13, and to narrow the list of possible candidate genes for Sym27 to ten, thus proving PMD to be useful for pea gene mapping and cloning. All information contained in PMD1 and PMD2 is available at www.peamarker.arriam.ru.

The Pea (Pisum sativum L.) Rogue Paramutation is Accompanied by Alterations in the Methylation Pattern of Specific Genomic Sequences

The spontaneous emergence among common pea (Pisum sativum L.) cultivars of off-type rogue plants exhibiting leaves with narrower and pointed leaflets and stipules and the non-Mendelian inheritance of this new phenotype were first described in the early 20th century. However, so far, no studies at the molecular level of this first identified case of paramutation have been carried out. In this study, we show for the first time that the pea rogue paramutation is accompanied by alterations in the methylation status of specific genomic sequences. Although, no significant differences were observed in the genome-wide DNA methylation in leaves of non-rogue cv. Onward in comparison to its rogue paramutant line JI2723, 22 DNA sequences were identified by methylation-sensitive amplified fragment length polymorphisms (MS-AFLP) analysis as differentially methylated in the two epigenomes. Mitotically inherited through all leaf tissues, the differential methylation patterns were also found to be meiotically inherited and conserved in pollen grains for 12 out of the 22 sequences. Fourteen of the sequences were successfully amplified in cDNA but none of them exhibited significant differential expression in the two contrasting epigenotypes. The further exploitation of the present research results on the way towards the elucidation of the molecular mechanisms behind this interesting epigenetic phenomenon is discussed.

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.

The fa2 gene and molecular markers mapping in the gp segment of the Pisum linkage group V

Review studies on the world Pisum genetic resources have shown that stem fasciation is controlled by three loci, i.e., fa1 (LGIV; Wt 10006 - type line of the Polish Gene Bank), fa2 (LGV, the line Wt 12185), and fas (LGIII, the line Shtambovii). Outstanding advantages of this character (e.g., pods gathered in upper part of a stem) resulted in breeding some cultivars. Preliminary investigations suggested linkages of the newly described fa2 gene within the gp-U segment. Based on the further linkage test crosses, it was stated that the fa2 is localized between the gp and Pis_Gen_9_3_1 markers (in the LGV). Additionally, four molecular markers (AD175, AB146, AC58, and AD280) and the morphological marker lk were also localized in this segment. Moreover, rms5, lum3, and cri were found to map on the other side of gp with tight linkage observed between lum3 and cri.

A novel er1 allele and the development and validation of its functional marker for breeding pea (Pisum sativum L.) resistance to powdery mildew

KEY MESSAGE

A novel er1 allele, er1 -7, conferring pea powdery mildew resistance was characterized by a 10-bp deletion in PsMLO1 cDNA, and its functional marker was developed and validated in pea germplasms. Pea powdery mildew caused by Erysiphe pisi DC is a major disease worldwide. Pea cultivar 'DDR-11' is an elite germplasm resistant to E. pisi. To identify the gene conferring resistance in DDR-11, the susceptible Bawan 6 and resistant DDR-11 cultivars were crossed to produce F1, F2, and F(2:3) populations. The phenotypic segregation patterns in the F2 and F(2:3) populations fit the 3:1 (susceptible:resistant) and 1:2:1 (susceptible homozygotes:heterozygotes:resistant homozygotes) ratios, respectively, indicating that resistance was controlled by a single recessive gene. Analysis of er1-linked markers in the F2 population suggested that the recessive resistance gene in DDR-11 was an er1 allele, which was mapped between markers ScOPE16-1600 and c5DNAmet. To further characterize er1 allele, the cDNA sequences of PsMLO1 from the parents were obtained and a novel er1 allele in DDR-11 was identified and designated as er1-7, which has a 10-bp deletion in position 111-120. The er1-7 allele caused a frame-shift mutation, resulting in a premature termination of translation of PsMLO1 protein. A co-dominant functional marker specific for er1-7 was developed, InDel111-120, which co-segregated with E. pisi resistance in the mapping population. The marker was able to distinguish between pea germplasms with and without the er1-7. Of 161 pea germplasms tested by InDel111-120, seven were detected containing resistance allele er1-7, which was verified by sequencing their PsMLO1 cDNA. Here, a novel er1 allele was characterized and its an ideal functional marker was validated, providing valuable genetic information and a powerful tool for breeding pea resistance to powdery mildew.

Ascochyta blight disease of pea (Pisum sativum L.): defence-related candidate genes associated with QTL regions and identification of epistatic QTL

KEY MESSAGE

Advances have been made in our understanding of Ascochyta blight resistance genetics through mapping candidate genes associated with QTL regions and demonstrating the importance of epistatic interactions in determining resistance. Ascochyta blight disease of pea (Pisum sativum L.) is economically significant with worldwide distribution. The causal pathogens are Didymella pinodes, Phoma medicaginis var pinodella and, in South Australia, P. koolunga. This study aimed to identify candidate genes that map to quantitative trait loci (QTL) for Ascochyta blight field disease resistance and to explore the role of epistatic interactions. Candidate genes associated with QTL were identified beginning with 101 defence-related genes from the published literature. Synteny between pea and Medicago truncatula was used to narrow down the candidates for mapping. Fourteen pea candidate sequences were mapped in two QTL mapping populations, A26 × Rovar and A88 × Rovar. QTL peaks, or the intervals containing QTL peaks, for the Asc2.1, Asc4.2, Asc4.3 and Asc7.1 QTL were defined by four of these candidate genes, while another three candidate genes occurred within 1.0 LOD confidence intervals. Epistasis involving QTL × background marker and background marker × background marker interactions contributed to the disease response phenotypes observed in the two mapping populations. For each population, five pairwise interactions exceeded the 5% false discovery rate threshold. Two candidate genes were involved in significant pairwise interactions. Markers in three genomic regions were involved in two or more epistatic interactions. Therefore, this study has identified pea defence-related sequences that are candidates for resistance determination, and that may be useful for marker-assisted selection. The demonstration of epistasis informs breeders that the architecture of this complex quantitative resistance includes epistatic interactions with non-additive effects.

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.

SNP discovery and genetic mapping using genotyping by sequencing of whole genome genomic DNA from a pea RIL population

BACKGROUND

Progress in genetics and breeding in pea still suffers from the limited availability of molecular resources. SNP markers that can be identified through affordable sequencing processes, without the need for prior genome reduction or a reference genome to assemble sequencing data would allow the discovery and genetic mapping of thousands of molecular markers. Such an approach could significantly speed up genetic studies and marker assisted breeding for non-model species.

RESULTS

A total of 419,024 SNPs were discovered using HiSeq whole genome sequencing of four pea lines, followed by direct identification of SNP markers without assembly using the discoSnp tool. Subsequent filtering led to the identification of 131,850 highly designable SNPs, polymorphic between at least two of the four pea lines. A subset of 64,754 SNPs was called and genotyped by short read sequencing on a subpopulation of 48 RILs from the cross 'Baccara' x 'PI180693'. This data was used to construct a WGGBS-derived pea genetic map comprising 64,263 markers. This map is collinear with previous pea consensus maps and therefore with the Medicago truncatula genome. Sequencing of four additional pea lines showed that 33 % to 64 % of the mapped SNPs, depending on the pairs of lines considered, are polymorphic and can therefore be useful in other crosses. The subsequent genotyping of a subset of 1000 SNPs, chosen for their mapping positions using a KASP™ assay, showed that almost all generated SNPs are highly designable and that most (95 %) deliver highly qualitative genotyping results. Using rather low sequencing coverages in SNP discovery and in SNP inferring did not hinder the identification of hundreds of thousands of high quality SNPs.

CONCLUSIONS

The development and optimization of appropriate tools in SNP discovery and genetic mapping have allowed us to make available a massive new genomic resource in pea. It will be useful for both fine mapping within chosen QTL confidence intervals and marker assisted breeding for important traits in pea improvement.

Discovery of a Novel er1 Allele Conferring Powdery Mildew Resistance in Chinese Pea (Pisum sativum L.) Landraces

Pea powdery mildew, caused by Erysiphe pisi D.C., is an important disease worldwide. Deployment of resistant varieties is the main way to control this disease. This study aimed to screen Chinese pea (Pisum sativum L.) landraces resistant to E. pisi, and to characterize the resistance gene(s) at the er1 locus in the resistant landraces, and to develop functional marker(s) specific to the novel er1 allele. The 322 landraces showed different resistance levels. Among them, 12 (3.73%), 4 (1.24%) and 17 (5.28%) landraces showed immunity, high resistance and resistance to E. pisi, respectively. The other landraces appeared susceptible or highly susceptible to E. pisi. Most of the immune and highly resistant landraces were collected from Yunnan province. To characterize the resistance gene at the er1 locus, cDNA sequences of PsMLO1 gene were determined in 12 immune and four highly resistant accessions. The cDNAs of PsMLO1 from the immune landrace G0005576 produced three distinct transcripts, characterized by a 129-bp deletion, and 155-bp and 220-bp insertions, which were consistent with those of er1-2 allele. The PsMLO1 cDNAs in the other 15 resistant landraces produced identical transcripts, which had a new point mutation (T→C) at position 1121 of PsMLO1, indicating a novel er1 allele, designated as er1-6. This mutation caused a leucine to proline change in the amino acid sequence. Subsequently, the resistance allele er1-6 in landrace G0001778 was confirmed by resistance inheritance analysis and genetic mapping on the region of the er1 locus using populations derived from G0001778 × Bawan 6. Finally, a functional marker specific to er1-6, SNP1121, was developed using the high-resolution melting technique, which could be used in pea breeding via marker-assisted selection. The results described here provide valuable genetic information for Chinese pea landraces and a powerful tool for pea breeders.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

Genomic Tools in Pea Breeding Programs: Status and Perspectives

Pea (Pisum sativum L.) is an annual cool-season legume and one of the oldest domesticated crops. Dry pea seeds contain 22-25% protein, complex starch and fiber constituents, and a rich array of vitamins, minerals, and phytochemicals which make them a valuable source for human consumption and livestock feed. Dry pea ranks third to common bean and chickpea as the most widely grown pulse in the world with more than 11 million tons produced in 2013. Pea breeding has achieved great success since the time of Mendel's experiments in the mid-1800s. However, several traits still require significant improvement for better yield stability in a larger growing area. Key breeding objectives in pea include improving biotic and abiotic stress resistance and enhancing yield components and seed quality. Taking advantage of the diversity present in the pea genepool, many mapping populations have been constructed in the last decades and efforts have been deployed to identify loci involved in the control of target traits and further introgress them into elite breeding materials. Pea now benefits from next-generation sequencing and high-throughput genotyping technologies that are paving the way for genome-wide association studies and genomic selection approaches. This review covers the significant development and deployment of genomic tools for pea breeding in recent years. Future prospects are discussed especially in light of current progress toward deciphering the pea genome.

Validation of QTL for resistance to Aphanomyces euteiches in different pea genetic backgrounds using near-isogenic lines

KEY MESSAGE

Marker-assisted backcrossing was used to generate pea NILs carrying individual or combined resistance alleles at main Aphanomyces resistance QTL. The effects of several QTL were successfully validated depending on genetic backgrounds. Quantitative trait loci (QTL) validation is an important and often overlooked step before subsequent research in QTL cloning or marker-assisted breeding for disease resistance in plants. Validation of QTL controlling partial resistance to Aphanomyces root rot, one of the most damaging diseases of pea worldwide, is of major interest for the future development of resistant varieties. The aim of this study was to validate, in different genetic backgrounds, the effects of various resistance alleles at seven main resistance QTL recently identified. Five backcross-assisted selection programs were developed. In each, resistance alleles at one to three of the seven main Aphanomyces resistance QTL were transferred into three genetic backgrounds, including two agronomically important spring (Eden) and winter (Isard) pea cultivars. The subsequent near-isogenic lines (NILs) were evaluated for resistance to two reference strains of the main A. euteiches pathotypes under controlled conditions. The NILs carrying resistance alleles at the major-effect QTL Ae-Ps4.5 and Ae-Ps7.6, either individually or in combination with resistance alleles at other QTL, showed significantly reduced disease severity compared to NILs without resistance alleles. Resistance alleles at some minor-effect QTL, especially Ae-Ps2.2 and Ae-Ps5.1, were also validated for their individual or combined effects on resistance. QTL × genetic background interactions were observed, mainly for QTL Ae-Ps7.6, the effect of which increased in the winter cultivar Isard. The pea NILs are a novel and valuable resource for further understanding the mechanisms underlying QTL and their integration in breeding programs.

High-Throughput Development of SSR Markers from Pea (Pisum sativum L.) Based on Next Generation Sequencing of a Purified Chinese Commercial Variety

Pea (Pisum sativum L.) is an important food legume globally, and is the plant species that J.G. Mendel used to lay the foundation of modern genetics. However, genomics resources of pea are limited comparing to other crop species. Application of marker assisted selection (MAS) in pea breeding has lagged behind many other crops. Development of a large number of novel and reliable SSR (simple sequence repeat) or microsatellite markers will help both basic and applied genomics research of this crop. The Illumina HiSeq 2500 System was used to uncover 8,899 putative SSR containing sequences, and 3,275 non-redundant primers were designed to amplify these SSRs. Among the 1,644 SSRs that were randomly selected for primer validation, 841 yielded reliable amplifications of detectable polymorphisms among 24 genotypes of cultivated pea (Pisum sativum L.) and wild relatives (P. fulvum Sm.) originated from diverse geographical locations. The dataset indicated that the allele number per locus ranged from 2 to 10, and that the polymorphism information content (PIC) ranged from 0.08 to 0.82 with an average of 0.38. These 1,644 novel SSR markers were also tested for polymorphism between genotypes G0003973 and G0005527. Finally, 33 polymorphic SSR markers were anchored on the genetic linkage map of G0003973 × G0005527 F2 population.

Assessment of genetic diversity in Ethiopian field pea (Pisum sativum L.) accessions with newly developed EST-SSR markers

BACKGROUND

Field pea (Pisum sativum L.) is among the prominent crops in the world as food and feed. There are relatively few simple sequence repeat (SSR) markers developed from expressed sequence tags (ESTs) in P. sativum.

RESULTS

In the present study, 15 new EST-SSR markers were developed from publicly available ESTs. These markers have successfully amplified their target loci across seven Pisum sativum subsp. sativum accessions. Eleven (73%) of these SSRs were trinucleotide repeats, two (13%) dinucleotide and two (13%) were hexanucleotide repeats. Across-taxa transferability of these new markers was also tested on other subspecies of Pisum as well as on P. fulvum, Vicia faba and Lens culinaris. In Pisum sativum subsp. sativum, 13 of the 15 markers were polymorphic and 12 of them subsequently used for genetic diversity analysis. Forty six accessions, of which 43 were from Ethiopia, were subjected to genetic diversity analysis using these newly developed markers. All accessions were represented by 12 individuals except two (NGB103816 and 237508) that were represented by 9 and 11 individuals, respectively. A total of 37 alleles were detected across all accessions. PS10 was the most polymorphic locus with six alleles, and the average number of alleles per locus over the 12 polymorphic loci was 3.1. Several rare and private alleles were also revealed. The most distinct accession (32048) had private alleles at three loci with 100% frequency.

CONCLUSION

These newly developed EST-SSR primer-pairs successfully amplified expected loci in P. sativum subsp. sativum as well as in other subspecies of the genus Pisum and related genera. High levels of genetic variation were detected in field pea accessions from Ethiopia using these markers. This result implies the potential of the Ethiopian field pea gene pool for improvement of field peas in various desirable traits. In addition, these markers could be a valuable asset in resolving the inconsistency in the taxonomic status of the different subspecies of genus Pisum as well as for characterization of field pea accessions in different gene banks around the world for breeding and conservation purposes.

Validation of SSR markers associated with rust (Uromyces fabae) resistance in pea (Pisum sativum L.)

Pea rust is a devastating disease of peas especially in the sub-tropical regions of the world and greatly influenced by the environmental conditions during disease development. Molecular markers associated with pea rust resistance would be useful in marker assisted selection (MAS). Utility of molecular markers associated with the pea rust resistance were evaluated in 30 diverse pea genotypes using four SSR markers (AA446 and AA505 flanking the major QTL Qruf; AD146 and AA416 flanking the minor QTL, Qruf1). QTL, Qruf flanking markers were able to identify all the resistant genotypes when used together, except Pant P 31. While, SSR markers AD146 and AA416 flanking the minor QTL, Qruf1 were able to identify all the pea resistant genotypes used for validation, except for HUDP-11 by AD146 and Pant P 31 by AA416. Similarly, SSR markers AA446 and AA505 were able to identify all the susceptible pea genotypes, except IPFD 99-13, HFP 9415 and S- 143. SSR markers AD146 and AA416 were together able to identify all the pea susceptible genotypes used for validation, except KPMR 526, KPMR 632 and IPFD 99-13. On the basis of marker allele analysis it may be concluded that SSR markers (AA446, AA505, AD146 and AA416) can be used in MAS of pea rust resistance.

Genetic diversity and trait genomic prediction in a pea diversity panel

BACKGROUND

Pea (Pisum sativum L.), a major pulse crop grown for its protein-rich seeds, is an important component of agroecological cropping systems in diverse regions of the world. New breeding challenges imposed by global climate change and new regulations urge pea breeders to undertake more efficient methods of selection and better take advantage of the large genetic diversity present in the Pisum sativum genepool. Diversity studies conducted so far in pea used Simple Sequence Repeat (SSR) and Retrotransposon Based Insertion Polymorphism (RBIP) markers. Recently, SNP marker panels have been developed that will be useful for genetic diversity assessment and marker-assisted selection.

RESULTS

A collection of diverse pea accessions, including landraces and cultivars of garden, field or fodder peas as well as wild peas was characterised at the molecular level using newly developed SNP markers, as well as SSR markers and RBIP markers. The three types of markers were used to describe the structure of the collection and revealed different pictures of the genetic diversity among the collection. SSR showed the fastest rate of evolution and RBIP the slowest rate of evolution, pointing to their contrasted mode of evolution. SNP markers were then used to predict phenotypes -the date of flowering (BegFlo), the number of seeds per plant (Nseed) and thousand seed weight (TSW)- that were recorded for the collection. Different statistical methods were tested including the LASSO (Least Absolute Shrinkage ans Selection Operator), PLS (Partial Least Squares), SPLS (Sparse Partial Least Squares), Bayes A, Bayes B and GBLUP (Genomic Best Linear Unbiased Prediction) methods and the structure of the collection was taken into account in the prediction. Despite a limited number of 331 markers used for prediction, TSW was reliably predicted.

CONCLUSION

The development of marker assisted selection has not reached its full potential in pea until now. This paper shows that the high-throughput SNP arrays that are being developed will most probably allow for a more efficient selection in this species.

Genetic diversity and population structure among pea (Pisum sativum L.) cultivars as revealed by simple sequence repeat and novel genic markers

Field pea (Pisum sativum L.) is an important cool season legume crop widely grown around the world. This research provides a basis for selection of pea germplasm across geographical regions in current and future breeding and genetic mapping efforts for pea improvement. Eleven novel genic markers were developed from pea expressed sequence tag (EST) sequences having significant similarity with gene calls from Medicago truncatula spanning at least one intron. In this study, 96 cultivars widely grown or used in breeding programs in the USA and Canada were analyzed for genetic diversity using 31 microsatellite or simple sequence repeat (SSR) and 11 novel EST-derived genic markers. The polymorphic information content varied from 0.01-0.56 among SSR markers and 0.04-0.43 among genic markers. The results showed that SSR and EST-derived genic markers displayed one or more highly reproducible, multi-allelic, and easy to score loci ranging from 200 to 700 bp in size. Genetic diversity was assessed through unweighted neighbor-joining method, and 96 varieties were grouped into three main clusters based on the dissimilarity matrix. Four subpopulations were determined through STRUCTURE analysis with no significant geographic separation of the subpopulations. The findings of the present study can be used to select diverse genotypes to be used as parents of crosses aimed for breeding improved pea cultivars.

Gene-based SNP discovery and genetic mapping in pea

KEY MESSAGE

Gene-based SNPs were identified and mapped in pea using five recombinant inbred line populations segregating for traits of agronomic importance. Pea (Pisum sativum L.) is one of the world's oldest domesticated crops and has been a model system in plant biology and genetics since the work of Gregor Mendel. Pea is the second most widely grown pulse crop in the world following common bean. The importance of pea as a food crop is growing due to its combination of moderate protein concentration, slowly digestible starch, high dietary fiber concentration, and its richness in micronutrients; however, pea has lagged behind other major crops in harnessing recent advances in molecular biology, genomics and bioinformatics, partly due to its large genome size with a large proportion of repetitive sequence, and to the relatively limited investment in research in this crop globally. The objective of this research was the development of a genome-wide transcriptome-based pea single-nucleotide polymorphism (SNP) marker platform using next-generation sequencing technology. A total of 1,536 polymorphic SNP loci selected from over 20,000 non-redundant SNPs identified using deep transcriptome sequencing of eight diverse Pisum accessions were used for genotyping in five RIL populations using an Illumina GoldenGate assay. The first high-density pea SNP map defining all seven linkage groups was generated by integrating with previously published anchor markers. Syntenic relationships of this map with the model legume Medicago truncatula and lentil (Lens culinaris Medik.) maps were established. The genic SNP map establishes a foundation for future molecular breeding efforts by enabling both the identification and tracking of introgression of genomic regions harbouring QTLs related to agronomic and seed quality traits.

QTL analysis of frost damage in pea suggests different mechanisms involved in frost tolerance

KEY MESSAGE

Avoidance mechanisms and intrinsic resistance are complementary strategies to improve winter frost tolerance and yield potential in field pea. The development of the winter pea crop represents a major challenge to expand plant protein production in temperate areas. Breeding winter cultivars requires the combination of freezing tolerance as well as high seed productivity and quality. In this context, we investigated the genetic determinism of winter frost tolerance and assessed its genetic relationship with yield and developmental traits. Using a newly identified source of frost resistance, we developed a population of recombinant inbred lines and evaluated it in six environments in Dijon and Clermont-Ferrand between 2005 and 2010. We developed a genetic map comprising 679 markers distributed over seven linkage groups and covering 947.1 cM. One hundred sixty-one quantitative trait loci (QTL) explaining 9-71 % of the phenotypic variation were detected across the six environments for all traits measured. Two clusters of QTL mapped on the linkage groups III and one cluster on LGVI reveal the genetic links between phenology, morphology, yield-related traits and frost tolerance in winter pea. QTL clusters on LGIII highlighted major developmental gene loci (Hr and Le) and the QTL cluster on LGVI explained up to 71 % of the winter frost damage variation. This suggests that a specific architecture and flowering ideotype defines frost tolerance in winter pea. However, two consistent frost tolerance QTL on LGV were independent of phenology and morphology traits, showing that different protective mechanisms are involved in frost tolerance. Finally, these results suggest that frost tolerance can be bred independently to seed productivity and quality.