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

Reference genome sequence and population genomic analysis of peas provide insights into the genetic basis of Mendelian and other agronomic traits

Peas are essential for human nutrition and played a crucial role in the discovery of Mendelian laws of inheritance. In this study, we assembled the genome of the elite vegetable pea cultivar 'Zhewan No. 1' at the chromosome level and analyzed resequencing data from 314 accessions, creating a comprehensive map of genetic variation in peas. We identified 235 candidate loci associated with 57 important agronomic traits through genome-wide association studies. Notably, we pinpointed the causal gene haplotypes responsible for four Mendelian traits: stem length (Le/le), flower color (A/a), cotyledon color (I/i) and seed shape (R/r). Additionally, we discovered the genes controlling pod form (Mendelian P/p) and hilum color. Our study also involved constructing a gene expression atlas across 22 tissues, highlighting key gene modules related to pod and seed development. These findings provide valuable pea genomic information and will facilitate the future genome-informed improvement of pea crops.

afila, the origin and nature of a major innovation in the history of pea breeding

The afila (af) mutation causes the replacement of leaflets by a branched mass of tendrils in the compound leaves of pea - Pisum sativum L. This mutation was first described in 1953, and several reports of spontaneous af mutations and induced mutants with a similar phenotype exist. Despite widespread introgression into breeding material, the nature of af and the origin of the alleles used remain unknown. Here, we combine comparative genomics with reverse genetic approaches to elucidate the genetic determinants of af. We also investigate haplotype diversity using a set of AfAf and afaf cultivars and breeding lines and molecular markers linked to seven consecutive genes. Our results show that deletion of two tandemly arranged genes encoding Q-type Cys(2)His(2) zinc finger transcription factors, PsPALM1a and PsPALM1b, is responsible for the af phenotype in pea. Eight haplotypes were identified in the af-harbouring genomic region on chromosome 2. These haplotypes differ in the size of the deletion, covering more or less genes. Diversity at the af locus is valuable for crop improvement and sheds light on the history of pea breeding for improved standing ability. The results will be used to understand the function of PsPALM1a/b and to transfer the knowledge for innovation in related crops.

Identification of QTLs associated with seed protein concentration in two diverse recombinant inbred line populations of pea

Improving the seed protein concentration (SPC) of pea (Pisum sativum L.) has turned into an important breeding objective because of the consumer demand for plant-based protein and demand from protein fractionation industries. To support the marker-assisted selection (MAS) of SPC towards accelerated breeding of improved cultivars, we have explored two diverse recombinant inbred line (RIL) populations to identify the quantitative trait loci (QTLs) associated with SPC. The two RIL populations, MP 1918 × P0540-91 (PR-30) and Ballet × Cameor (PR-31), were derived from crosses between moderate SPC × high SPC accessions. A total of 166 and 159 RILs of PR-30 and PR-31, respectively, were genotyped using an Axiom® 90K SNP array and 13.2K SNP arrays, respectively. The RILs were phenotyped in replicated trials in two and three locations of Saskatchewan, Canada in 2020 and 2021, respectively, for agronomic assessment and SPC. Using composite interval mapping, we identified three QTLs associated with SPC in PR-30 and five QTLs in PR-31, with the LOD value ranging from 3.0 to 11.0. A majority of these QTLs were unique to these populations compared to the previously known QTLs for SPC. The QTL SPC-Ps-5.1 overlapped with the earlier reported SPC associated QTL PC-QTL-3. Three QTLs, SPC-Ps-4.2, SPC-Ps-5.1, and SPC-Ps-7.2 with LOD scores of 7.2, 7.9, and 11.3, and which explained 14.5%, 11.6%, and 11.3% of the phenotypic variance, respectively, can be used for marker-assisted breeding to increase SPC in peas. Eight QTLs associated with the grain yield were identified with LOD scores ranging from 3.1 to 8.2. Two sets of QTLs, SPC-Ps-2.1 and GY-Ps-2.1, and SPC-Ps-5.1 and GY-Ps-5.3, shared the QTL/peak regions. Each set of QTLs contributed to either SPC or grain yield depending on which parent the QTL region is derived from, thus confirming that breeding for SPC should take into consideration the effects on grain yield.

Physical seed dormancy in pea is genetically separable from seed coat thickness and roughness

Introduction

The seeds of wild pea (Pisum) exhibit marked physical dormancy due to impermeability of the seed coat to water, and the loss of this dormancy is thought to have been critical for domestication. Wild pea seed coats are also notably thick and rough, traits that have also reduced during domestication and are anecdotally linked to increased permeability. However, how these traits specifically interact with permeability is unclear.

Methods

To investigate this, we examined the genetic control of differences in seed coat characteristics between wild P. sativum ssp. humile and a non-dormant domesticated P. s. sativum accession in a recombinant inbred population. QTL effects were confirmed and their locations refined in segregating F4/5 populations.

Results

In this population we found a moderate correlation between testa thickness and permeability, and identified loci that affect them independently, suggesting no close functional association. However, the major loci affecting both testa thickness and permeability collocated closely with Mendel's pigmentation locus A, suggesting flavonoid compounds under its control might contribute significantly to both traits. We also show that seed coat roughness is oligogenic in this population, with the major locus independent of both testa thickness and permeability, suggesting selection for smooth seed was unlikely to be due to effects on either of these traits.

Discussion

Results indicate loss of seed coat dormancy during domestication was not primarily driven by reduced testa thickness or smooth seededness. The close association between major permeability and thickness QTL and Mendel's 'A' warrant further study, particularly regarding the role of flavonoids.

Fine mapping of Ae-Ps4.5, a major locus for resistance to pathotype III of Aphanomyces euteiches in pea

KEY MESSAGE

QTL mapping and recombinant screening confirmed the major effect of QTL Ae-Ps4.5 on pea resistance to pathotype III of Aphanomyces euteiches and fine-mapped the QTL to a 3.06-Mb interval. Aphanomyces root rot, caused by Aphanomyces euteiches, is the most important disease of pea (Pisum sativum L.) worldwide. The development of pea-resistant varieties is a major challenge to control the disease. Previous linkage studies identified seven main resistance quantitative trait loci (QTL), including the QTL Ae-Ps4.5 associated with partial resistance in US nurseries infested by the pea pathotype III of A. euteiches. This study aimed to confirm the major effect of Ae-Ps4.5 on A. euteiches pathotype III, refine its interval, and identify candidate genes underlying the QTL. QTL mapping on an updated genetic map from the Puget × 90-2079 pea recombinant inbred line population identified Ae-Ps4.5 in a 0.8-cM confidence interval with a high effect (R2 = 89%) for resistance to the Ae109 reference strain of A. euteiches (pathotype III) under controlled conditions. However, the QTL mapping did not detect Ae-Ps4.5 for resistance to the RB84 reference strain of A. euteiches (pathotype I). Screening 224-pea BC5F2 plant progeny derived from three near-isogenic lines (NILs) carrying the 90-2079 allele at Ae-Ps4.5 in the Puget genetic background with 26 SNP markers identified 15 NILs showing recombination in the QTL interval. Phenotyping of the recombinant lines for resistance to the Ae109 strain of A. euteiches reduced the QTL to a physical interval of 3.06 Mb, containing 50 putative annotated genes on the Caméor pea genome V1a among which three candidate genes highlighted. This study provides closely linked SNP markers and putative candidate genes to accelerate pea breeding for resistant varieties to Aphanomyces root rot.

Recombinant inbred lines derived from wide crosses in Pisum

Genomic resources are becoming available for Pisum but to link these to phenotypic diversity requires well marked populations segregating for relevant traits. Here we describe two such resources. Two recombinant inbred populations, derived from wide crosses in Pisum are described. One high resolution mapping population involves cv Caméor, for which the first pea whole genome assembly was obtained, crossed to JI0281, a basally divergent P. sativum sativum landrace from Ethiopia. The other is an inter sub-specific cross between P. s. sativum and the independently domesticated P. s. abyssinicum. The corresponding genetic maps provide information on chromosome level sequence assemblies and identify structural differences between the genomes of these two Pisum subspecies. In order to visualise chromosomal translocations that distinguish the mapping parents, we created a simplified version of Threadmapper to optimise it for interactive 3-dimensional display of multiple linkage groups. The genetic mapping of traits affecting seed coat roughness and colour, plant height, axil ring pigmentation, leaflet number and leaflet indentation enabled the definition of their corresponding genomic regions. The consequence of structural rearrangement for trait analysis is illustrated by leaf serration. These analyses pave the way for identification of the underlying genes and illustrate the utility of these publicly available resources. Segregating inbred populations derived from wide crosses in Pisum, together with the associated marker data, are made publicly available for trait dissection. Genetic analysis of these populations is informative about chromosome scale assemblies, structural diversity in the pea genome and has been useful for the fine mapping of several discrete and quantitative traits.

Advanced backcross QTL analysis and comparative mapping with RIL QTL studies and GWAS provide an overview of QTL and marker haplotype diversity for resistance to Aphanomyces root rot in pea (Pisum sativum)

Aphanomyces euteiches is the most damaging soilborne pea pathogen in France. Breeding of pea resistant varieties combining a diversity of quantitative trait loci (QTL) is a promising strategy considering previous research achievements in dissecting polygenic resistance to A. euteiches. The objective of this study was to provide an overview of the diversity of QTL and marker haplotypes for resistance to A. euteiches, by integrating a novel QTL mapping study in advanced backcross (AB) populations with previous QTL analyses and genome-wide association study (GWAS) using common markers. QTL analysis was performed in two AB populations derived from the cross between the susceptible spring pea variety "Eden" and the two new sources of partial resistance "E11" and "LISA". The two AB populations were genotyped using 993 and 478 single nucleotide polymorphism (SNP) markers, respectively, and phenotyped for resistance to A. euteiches in controlled conditions and in infested fields at two locations. GWAS and QTL mapping previously reported in the pea-Aphanomyces collection and from four recombinant inbred line (RIL) populations, respectively, were updated using a total of 1,850 additional markers, including the markers used in the Eden x E11 and Eden x LISA populations analysis. A total of 29 resistance-associated SNPs and 171 resistance QTL were identified by GWAS and RIL or AB QTL analyses, respectively, which highlighted 10 consistent genetic regions confirming the previously reported QTL. No new consistent resistance QTL was detected from both Eden x E11 and Eden x LISA AB populations. However, a high diversity of resistance haplotypes was identified at 11 linkage disequilibrium (LD) blocks underlying consistent genetic regions, especially in 14 new sources of resistance from the pea-Aphanomyces collection. An accumulation of favorable haplotypes at these 11 blocks was confirmed in the most resistant pea lines of the collection. This study provides new SNP markers and rare haplotypes associated with the diversity of Aphanomyces root rot resistance QTL investigated, which will be useful for QTL pyramiding strategies to increase resistance levels in future pea varieties.

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 Mechanisms for Hybrid Breeding in Vegetable Crops

To address the complex challenges faced by our planet such as rapidly changing climate patterns, food and nutritional insecurities, and the escalating world population, the development of hybrid vegetable crops is imperative. Vegetable hybrids could effectively mitigate the above-mentioned fundamental challenges in numerous countries. Utilizing genetic mechanisms to create hybrids not only reduces costs but also holds significant practical implications, particularly in streamlining hybrid seed production. These mechanisms encompass self-incompatibility (SI), male sterility, and gynoecism. The present comprehensive review is primarily focused on the elucidation of fundamental processes associated with floral characteristics, the genetic regulation of floral traits, pollen biology, and development. Specific attention is given to the mechanisms for masculinizing and feminizing cucurbits to facilitate hybrid seed production as well as the hybridization approaches used in the biofortification of vegetable crops. Furthermore, this review provides valuable insights into recent biotechnological advancements and their future utilization for developing the genetic systems of major vegetable crops.

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

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

Assembly of the 81.6 Mb centromere of pea chromosome 6 elucidates the structure and evolution of metapolycentric chromosomes

Centromeres in the legume genera Pisum and Lathyrus exhibit unique morphological characteristics, including extended primary constrictions and multiple separate domains of centromeric chromatin. These so-called metapolycentromeres resemble an intermediate form between monocentric and holocentric types, and therefore provide a great opportunity for studying the transitions between different types of centromere organizations. However, because of the exceedingly large and highly repetitive nature of metapolycentromeres, highly contiguous assemblies needed for these studies are lacking. Here, we report on the assembly and analysis of a 177.6 Mb region of pea (Pisum sativum) chromosome 6, including the 81.6 Mb centromere region (CEN6) and adjacent chromosome arms. Genes, DNA methylation profiles, and most of the repeats were uniformly distributed within the centromere, and their densities in CEN6 and chromosome arms were similar. The exception was an accumulation of satellite DNA in CEN6, where it formed multiple arrays up to 2 Mb in length. Centromeric chromatin, characterized by the presence of the CENH3 protein, was predominantly associated with arrays of three different satellite repeats; however, five other satellites present in CEN6 lacked CENH3. The presence of CENH3 chromatin was found to determine the spatial distribution of the respective satellites during the cell cycle. Finally, oligo-FISH painting experiments, performed using probes specifically designed to label the genomic regions corresponding to CEN6 in Pisum, Lathyrus, and Vicia species, revealed that metapolycentromeres evolved via the expansion of centromeric chromatin into neighboring chromosomal regions and the accumulation of novel satellite repeats. However, in some of these species, centromere evolution also involved chromosomal translocations and centromere repositioning.

Mendel: From genes to genome

Two hundred years after the birth of Gregor Mendel, it is an appropriate time to reflect on recent developments in the discipline of genetics, particularly advances relating to the prescient friar's model species, the garden pea (Pisum sativum L.). Mendel's study of seven characteristics established the laws of segregation and independent assortment. The genes underlying four of Mendel's loci (A, LE, I, and R) have been characterized at the molecular level for over a decade. However, the three remaining genes, influencing pod color (GP), pod form (V/P), and the position of flowers (FA/FAS), have remained elusive for a variety of reasons, including a lack of detail regarding the loci with which Mendel worked. Here, we discuss potential candidate genes for these characteristics, in light of recent advances in the genetic resources for pea. These advances, including the pea genome sequence and reverse-genetics techniques, have revitalized pea as an excellent model species for physiological-genetic studies. We also discuss the issues that have been raised with Mendel's results, such as the recent controversy regarding the discrete nature of the characters that Mendel chose and the perceived overly-good fit of his segregations to his hypotheses. We also consider the relevance of these controversies to his lasting contribution. Finally, we discuss the use of Mendel's classical results to teach and enthuse future generations of geneticists, not only regarding the core principles of the discipline, but also its history and the role of hypothesis testing.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Improved pea reference genome and pan-genome highlight genomic features and evolutionary characteristics

Complete and accurate reference genomes and annotations provide fundamental resources for functional genomics and crop breeding. Here we report a de novo assembly and annotation of a pea cultivar ZW6 with contig N50 of 8.98 Mb, which features a 243-fold increase in contig length and evident improvements in the continuity and quality of sequence in complex repeat regions compared with the existing one. Genome diversity of 118 cultivated and wild pea demonstrated that Pisum abyssinicum is a separate species different from P. fulvum and P. sativum within Pisum. Quantitative trait locus analyses uncovered two known Mendel's genes related to stem length (Le/le) and seed shape (R/r) as well as some candidate genes for pod form studied by Mendel. A pan-genome of 116 pea accessions was constructed, and pan-genes preferred in P. abyssinicum and P. fulvum showed distinct functional enrichment, indicating the potential value of them as pea breeding resources in the future.

Recent advancement in OMICS approaches to enhance abiotic stress tolerance in legumes

The world is facing rapid climate change and a fast-growing global population. It is believed that the world population will be 9.7 billion in 2050. However, recent agriculture production is not enough to feed the current population of 7.9 billion people, which is causing a huge hunger problem. Therefore, feeding the 9.7 billion population in 2050 will be a huge target. Climate change is becoming a huge threat to global agricultural production, and it is expected to become the worst threat to it in the upcoming years. Keeping this in view, it is very important to breed climate-resilient plants. Legumes are considered an important pillar of the agriculture production system and a great source of high-quality protein, minerals, and vitamins. During the last two decades, advancements in OMICs technology revolutionized plant breeding and emerged as a crop-saving tool in wake of the climate change. Various OMICs approaches like Next-Generation sequencing (NGS), Transcriptomics, Proteomics, and Metabolomics have been used in legumes under abiotic stresses. The scientific community successfully utilized these platforms and investigated the Quantitative Trait Loci (QTL), linked markers through genome-wide association studies, and developed KASP markers that can be helpful for the marker-assisted breeding of legumes. Gene-editing techniques have been successfully proven for soybean, cowpea, chickpea, and model legumes such as Medicago truncatula and Lotus japonicus. A number of efforts have been made to perform gene editing in legumes. Moreover, the scientific community did a great job of identifying various genes involved in the metabolic pathways and utilizing the resulted information in the development of climate-resilient legume cultivars at a rapid pace. Keeping in view, this review highlights the contribution of OMICs approaches to abiotic stresses in legumes. We envisage that the presented information will be helpful for the scientific community to develop climate-resilient legume cultivars.

Large-Scale Heat-Tolerance Screening and Genetic Diversity of Pea (Pisum sativum L.) Germplasms

Pea (Pisum sativum L.) is an important legume crop. However, the yield of pea is adversely affected by heat stress in China. In this study, heat-tolerant germplasms were screened and evaluated in the field under multi-conditions. The results showed that heat stress could significantly affect pea yield. On the basis of grain weight per plant, 257 heat-tolerant and 175 heat-sensitive accessions were obtained from the first year’s screening, and 26 extremely heat-tolerant and 19 extremely heat-sensitive accessions were finally obtained in this study. Based on SNaPshot technology, two sets of SNP markers, including 46 neutral and 20 heat-tolerance-related markers, were used to evaluate the genetic diversity and population genetic structure of the 432 pea accessions obtained from the first year’s screening. Genetic diversity analysis showed that the average polymorphic information content was lower using heat-tolerance-related markers than neutral markers because of the selective pressure under heat stress. In addition, population genetic structure analysis showed that neutral markers divided the 432 pea accessions into two subpopulations associated with sowing date type and geographical origin, while the heat-tolerance-related markers divided these germplasms into two subpopulations associated with heat tolerance and sowing date type. Overall, we present a comprehensive resource of heat-tolerant and heat-sensitive pea accessions through heat-tolerance screenings in multi-conditions, which could help genetic improvements of pea in the future.

Identification of Novel Genes Associated with Partial Resistance to Aphanomyces Root Rot in Field Pea by BSR-Seq Analysis

Aphanomyces root rot, caused by Aphanomyces euteiches, causes severe yield loss in field pea (Pisum sativum). The identification of a pea germplasm resistant to this disease is an important breeding objective. Polygenetic resistance has been reported in the field pea cultivar '00-2067'. To facilitate marker-assisted selection (MAS), bulked segregant RNA-seq (BSR-seq) analysis was conducted using an F8 RIL population derived from the cross of 'Carman' × '00-2067'. Root rot development was assessed under controlled conditions in replicated experiments. Resistant (R) and susceptible (S) bulks were constructed based on the root rot severity in a greenhouse study. The BSR-seq analysis of the R bulks generated 44,595,510~51,658,688 reads, of which the aligned sequences were linked to 44,757 genes in a reference genome. In total, 2356 differentially expressed genes were identified, of which 44 were used for gene annotation, including defense-related pathways (jasmonate, ethylene and salicylate) and the GO biological process. A total of 344.1 K SNPs were identified between the R and S bulks, of which 395 variants were located in 31 candidate genes. The identification of novel genes associated with partial resistance to Aphanomyces root rot in field pea by BSR-seq may facilitate efforts to improve management of this important disease.

Genetic Augmentation of Legume Crops Using Genomic Resources and Genotyping Platforms for Nutritional Food Security

Recent advances in next generation sequencing (NGS) technologies have led the surge of genomic resources for the improvement legume crops. Advances in high throughput genotyping (HTG) and high throughput phenotyping (HTP) enable legume breeders to improve legume crops more precisely and efficiently. Now, the legume breeder can reshuffle the natural gene combinations of their choice to enhance the genetic potential of crops. These genomic resources are efficiently deployed through molecular breeding approaches for genetic augmentation of important legume crops, such as chickpea, cowpea, pigeonpea, groundnut, common bean, lentil, pea, as well as other underutilized legume crops. In the future, advances in NGS, HTG, and HTP technologies will help in the identification and assembly of superior haplotypes to tailor the legume crop varieties through haplotype-based breeding. This review article focuses on the recent development of genomic resource databases and their deployment in legume molecular breeding programmes to secure global food security.

Combining NDVI and Bacterial Blight Score to Predict Grain Yield in Field Pea

Field pea is the most commonly grown temperate pulse crop, with close to 15 million tons produced globally in 2020. Varieties improved through breeding are important to ensure ongoing improvements in yield and disease resistance. Genomic selection (GS) is a modern breeding approach that could substantially improve the rate of genetic gain for grain yield, and its deployment depends on the prediction accuracy (PA) that can be achieved. In our study, four yield trials representing breeding lines' advancement stages of the breeding program (S0, S1, S2, and S3) were assessed with grain yield, aerial high-throughput phenotyping (normalized difference vegetation index, NDVI), and bacterial blight disease scores (BBSC). Low-to-moderate broad-sense heritability (0.31-0.71) and narrow-sense heritability (0.13-0.71) were observed, as the estimated additive and non-additive genetic components for the three traits varied with the different models fitted. The genetic correlations among the three traits were high, particularly in the S0-S2 stages. NDVI and BBSC were combined to investigate the PA for grain yield by univariate and multivariate GS models, and multivariate models showed higher PA than univariate models in both cross-validation and forward prediction methods. A 6-50% improvement in PA was achieved when multivariate models were deployed. The highest PA was indicated in the forward prediction scenario when the training population consisted of early generation breeding stages with the multivariate models. Both NDVI and BBSC are commonly used traits that could be measured in the early growth stage; however, our study suggested that NDVI is a more useful trait to predict grain yield with high accuracy in the field pea breeding program, especially in diseased trials, through its incorporation into multivariate models.

The genetic architecture of flowering time changes in pea from wild to crop

Change in phenology has been an important component in crop evolution, and selection for earlier flowering through a reduction in environmental sensitivity has helped broaden adaptation in many species. Natural variation for flowering in domesticated pea (Pisum sativum L.) has been noted and studied for decades, but there has been no clear account of change relative to its wild progenitor. Here we examined the genetic control of differences in flowering time between wild P. sativum ssp. humile and a typical late-flowering photoperiodic P. s. sativum accession in a recombinant inbred population under long and short photoperiods. Our results confirm the importance of the major photoperiod sensitivity locus Hr/PsELF3a and identify two other loci on chromosomes 1 (DTF1) and 3 (DTF3) that contribute to earlier flowering in the domesticated line under both photoperiods. The domesticated allele at a fourth locus on chromosome 6 (DTF6) delays flowering under long days only. Map positions, inheritance patterns, and expression analyses in near-isogenic comparisons imply that DTF1, DTF3, and DTF6 represent gain-of-function alleles of the florigen/antiflorigen genes FTa3, FTa1, and TFL1c/LF, respectively. This echoes similar variation in chickpea and lentil, and suggests a conserved route to reduced photoperiod sensitivity and early phenology in temperate pulses.

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.

Phenotypic and Genotypic Characterization and Correlation Analysis of Pea (Pisum sativum L.) Diversity Panel

Phenotypic and genotypic characterization were performed to assess heritability, variability, and seed yield stability of pea genotypes used in breeding to increase the pea production area. A European pea diversity panel, including genotypes from North America, Asia, and Australia consisting of varieties, breeding lines, pea, and landraces was examined in 2019 and 2020 in Serbia and Belgium using augmented block design. The highest heritability was for thousand seed weight; the highest coefficient of variation was for seed yield. The highest positive correlation was between number of seeds per plant and number of pods per plant; the highest negative correlation was between seed yield and protein content. Hierarchical clustering separated pea germplasm based on use and type. Different Principal component analysis grouping of landraces, breeding lines, and varieties, as well as forage types and garden and dry peas, confirms that there was an apparent decrease in similarity between the genotypes, which can be explained by their different purposes. Pea breeding should be focused on traits with consistent heritability and a positive effect on seed yield when selecting high-yielding genotypes, and on allowing for more widespread use of pea in various agricultural production systems.

Phylogeny and disparate selection signatures suggest two genetically independent domestication events in pea (Pisum L.)

Domestication is considered a model of adaptation that can be used to draw conclusions about the modus operandi of selection in natural systems. Investigating domestication may give insights into how plants react to different intensities of human manipulation, which has direct implication for the continuing efforts of crop improvement. Therefore, scientists of various disciplines study domestication-related questions to understand the biological and cultural bases of the domestication process. We employed restriction site-associated DNA sequencing (RAD-seq) of 494 Pisum sativum (pea) samples from all wild and domesticated groups to analyze the genetic structure of the collection. Patterns of ancient admixture were investigated by analysis of admixture graphs. We used two complementary approaches, one diversity based and one based on differentiation, to detect the selection signatures putatively associated with domestication. An analysis of the subpopulation structure of wild P. sativum revealed five distinct groups with a notable geographic pattern. Pisum abyssinicum clustered unequivocally within the P. sativum complex, without any indication of hybrid origin. We detected 32 genomic regions putatively subjected to selection: 29 in P. sativum ssp. sativum and three in P. abyssinicum. The two domesticated groups did not share regions under selection and did not display similar haplotype patterns within those regions. Wild P. sativum is structured into well-diverged subgroups. Although Pisum sativum ssp. humile is not supported as a taxonomic entity, the so-called 'southern humile' is a genuine wild group. Introgression did not shape the variation observed within the sampled germplasm. The two domesticated pea groups display distinct genetic bases of domestication, suggesting two genetically independent domestication events.

An Integrated Linkage Map of Three Recombinant Inbred Populations of Pea (Pisum sativum L.)

Biparental recombinant inbred line (RIL) populations are sets of genetically stable lines and have a simple population structure that facilitates the dissection of the genetics of interesting traits. On the other hand, populations derived from multiparent intercrosses combine both greater diversity and higher numbers of recombination events than RILs. Here, we describe a simple population structure: a three-way recombinant inbred population combination. This structure was easy to produce and was a compromise between biparental and multiparent populations. We show that this structure had advantages when analyzing cultivar crosses, and could achieve a mapping resolution of a few genes.

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.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

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

KEY MESSAGE

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

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.

Integration of the SMXL/D53 strigolactone signalling repressors in the model of shoot branching regulation in Pisum sativum

DWARF53 (D53) in rice (Oryza sativa) and its homologs in Arabidopsis (Arabidopsis thaliana), SUPPRESSOR OF MAX2-LIKE 6 (SMXL6), SMXL7 and SMXL8, are well established negative regulators of strigolactone (SL) signalling in shoot branching regulation. Little is known of pea (Pisum sativum) homologs and whether D53 and related SMXLs are specific to SL signalling pathways. Here, we identify two allelic pea mutants, dormant3 (dor3), and demonstrate through gene mapping and sequencing that DOR3 corresponds to a homolog of D53 and SMXL6/SMXL7, designated PsSMXL7. Phenotype analysis, gene expression, protein and hormone quantification assays were performed to determine the role of PsSMXL7 in regulation of bud outgrowth and the role of PsSMXL7 and D53 in integrating SL and cytokinin (CK) responses. Like D53 and related SMXLs, we show that PsSMXL7 can be degraded by SL and induces feedback upregulation of PsSMXL7 transcript. Here we reveal a system conserved in pea and rice, whereby CK also upregulates PsSMXL7/D53 transcripts, providing a clear mechanism for SL and CK cross-talk in the regulation of branching. To further deepen our understanding of the branching network in pea, we provide evidence that SL acts via PsSMXL7 to modulate auxin content via PsAFB5, which itself regulates expression of SL biosynthesis genes. We therefore show that PsSMXL7 is key to a triple hormone network involving an auxin-SL feedback mechanism and SL-CK cross-talk.

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.

Symbiotic Regulatory Genes Controlling Nodule Development in Pisum sativum L

Analyses of natural variation and the use of mutagenesis and molecular-biological approaches have revealed 50 symbiotic regulatory genes in pea (Pisum sativum L.). Studies of genomic synteny using model legumes, such as Medicago truncatula Gaertn. and Lotus japonicus (Regel) K. Larsen, have identified the sequences of 15 symbiotic regulatory genes in pea. These genes encode receptor kinases, an ion channel, a calcium/calmodulin-dependent protein kinase, transcription factors, a metal transporter, and an enzyme. This review summarizes and describes mutant alleles, their phenotypic manifestations, and the functions of all identified symbiotic regulatory genes in pea. Some examples of gene interactions are also given. In the review, all mutant alleles in genes with identified sequences are designated and still-unidentified symbiotic regulatory genes of great interest are considered. The identification of these genes will help elucidate additional components involved in infection thread growth, nodule primordium development, bacteroid differentiation and maintenance, and the autoregulation of nodulation. The significance of symbiotic mutants of pea as extremely fruitful genetic models for studying nodule development and for comparative cell biology studies of legume nodules is clearly demonstrated. Finally, it is noted that many more sequences of symbiotic regulatory genes remain to be identified. Transcriptomics approaches and genome-wide sequencing could help address this challenge.

Meta-analysis of QTL reveals the genetic control of yield-related traits and seed protein content in pea

Pea is one of the most important grain legume crops in temperate regions worldwide. Improving pea yield is a critical breeding target. Nine inter-connected pea recombinant inbred line populations were evaluated in nine environments at INRAE Dijon, France and genotyped using the GenoPea 13.2 K SNP array. Each population has been evaluated in two to four environments. A multi-population Quantitative Trait Loci (QTL) analysis for seed weight per plant (SW), seed number per plant (SN), thousand seed weight (TSW) and seed protein content (SPC) was done. QTL were then projected on the multi-population consensus map and a meta-analysis of QTL was performed. This analysis identified 17 QTL for SW, 16 QTL for SN, 35 QTL for TSW and 21 QTL for SPC, shedding light on trait relationships. These QTL were resolved into 27 metaQTL. Some of them showed small confidence intervals of less than 2 cM encompassing less than one hundred underlying candidate genes. The precision of metaQTL and the potential candidate genes reported in this study enable their use for marker-assisted selection and provide a foundation towards map-based identification of causal polymorphisms.

Genome-wide association study identifies favorable SNP alleles and candidate genes for frost tolerance in pea

BACKGROUND

Frost is a limiting abiotic stress for the winter pea crop (Pisum sativum L.) and identifying the genetic determinants of frost tolerance is a major issue to breed varieties for cold northern areas. Quantitative trait loci (QTLs) have previously been detected from bi-parental mapping populations, giving an overview of the genome regions governing this trait. The recent development of high-throughput genotyping tools for pea brings the opportunity to undertake genetic association studies in order to capture a higher allelic diversity within large collections of genetic resources as well as to refine the localization of the causal polymorphisms thanks to the high marker density. In this study, a genome-wide association study (GWAS) was performed using a set of 365 pea accessions. Phenotyping was carried out by scoring frost damages in the field and in controlled conditions. The association mapping collection was also genotyped using an Illumina Infinium® BeadChip, which allowed to collect data for 11,366 single nucleotide polymorphism (SNP) markers.

RESULTS

GWAS identified 62 SNPs significantly associated with frost tolerance and distributed over six of the seven pea linkage groups (LGs). These results confirmed 3 QTLs that were already mapped in multiple environments on LG III, V and VI with bi-parental populations. They also allowed to identify one locus, on LG II, which has not been detected yet and two loci, on LGs I and VII, which have formerly been detected in only one environment. Fifty candidate genes corresponding to annotated significant SNPs, or SNPs in strong linkage disequilibrium with the formers, were found to underlie the frost damage (FD)-related loci detected by GWAS. Additionally, the analyses allowed to define favorable haplotypes of markers for the FD-related loci and their corresponding accessions within the association mapping collection.

CONCLUSIONS

This study led to identify FD-related loci as well as corresponding favorable haplotypes of markers and representative pea accessions that might to be used in winter pea breeding programs. Among the candidate genes highlighted at the identified FD-related loci, the results also encourage further attention to the presence of C-repeat Binding Factors (CBF) as potential genetic determinants of the frost tolerance locus on LG VI.

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.

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

Faba bean (Vicia faba L.) is a pulse crop of high nutritional value and high importance for sustainable agriculture and soil protection. With the objective of identifying gene-based SNPs, transcriptome sequencing was performed in order to reduce faba bean genome complexity. A set of 1,819 gene-based SNP markers polymorphic in three recombinant line populations was selected to enable the construction of a high-density consensus genetic map encompassing 1,728 markers well distributed in six linkage groups and spanning 1,547.71 cM with an average inter-marker distance of 0.89 cM. Orthology-based comparison of the faba bean consensus map with legume genome assemblies highlighted synteny patterns that partly reflected the phylogenetic relationships among species. Solid blocks of macrosynteny were observed between faba bean and the most closely-related sequenced legume species such as pea, barrel medic or chickpea. Numerous blocks could also be identified in more divergent species such as common bean or cowpea. The genetic tools developed in this work can be used in association mapping, genetic diversity, linkage disequilibrium or comparative genomics and provide a backbone for map-based cloning. This will make the identification of candidate genes of interest more efficient and will accelerate marker-assisted selection (MAS) and genomic-assisted breeding (GAB) in faba bean.

Development and Proof-of-Concept Application of Genome-Enabled Selection for Pea Grain Yield under Severe Terminal Drought

Terminal drought is the main stress limiting pea (Pisum sativum L.) grain yield in Mediterranean environments. This study aimed to investigate genotype × environment (GE) interaction patterns, define a genomic selection (GS) model for yield under severe drought based on single nucleotide polymorphism (SNP) markers from genotyping-by-sequencing, and compare GS with phenotypic selection (PS) and marker-assisted selection (MAS). Some 288 lines belonging to three connected RIL populations were evaluated in a managed-stress (MS) environment of Northern Italy, Marchouch (Morocco), and Alger (Algeria). Intra-environment, cross-environment, and cross-population predictive ability were assessed by Ridge Regression best linear unbiased prediction (rrBLUP) and Bayesian Lasso models. GE interaction was particularly large across moderate-stress and severe-stress environments. In proof-of-concept experiments performed in a MS environment, GS models constructed from MS environment and Marchouch data applied to independent material separated top-performing lines from mid- and bottom-performing ones, and produced actual yield gains similar to PS. The latter result would imply somewhat greater GS efficiency when considering same selection costs, in partial agreement with predicted efficiency results. GS, which exploited drought escape and intrinsic drought tolerance, exhibited 18% greater selection efficiency than MAS (albeit with non-significant difference between selections) and moderate to high cross-population predictive ability. GS can be cost-efficient to raise yields under severe drought.

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.

Genome-Wide Association Mapping for Agronomic and Seed Quality Traits of Field Pea (Pisum sativum L.)

Genome-wide association study (GWAS) was conducted to identify loci associated with agronomic (days to flowering, days to maturity, plant height, seed yield and seed weight), seed morphology (shape and dimpling), and seed quality (protein, starch, and fiber concentrations) traits of field pea (Pisum sativum L.). A collection of 135 pea accessions from 23 different breeding programs in Africa (Ethiopia), Asia (India), Australia, Europe (Belarus, Czech Republic, Denmark, France, Lithuania, Netherlands, Russia, Sweden, Ukraine and United Kingdom), and North America (Canada and USA), was used for the GWAS. The accessions were genotyped using genotyping-by-sequencing (GBS). After filtering for a minimum read depth of five, and minor allele frequency of 0.05, 16,877 high quality SNPs were selected to determine marker-trait associations (MTA). The LD decay (LD_1/2max,90) across the chromosomes varied from 20 to 80 kb. Population structure analysis grouped the accessions into nine subpopulations. The accessions were evaluated in multi-year, multi-location trials in Olomouc (Czech Republic), Fargo, North Dakota (USA), and Rosthern and Sutherland, Saskatchewan (Canada) from 2013 to 2017. Each trait was phenotyped in at least five location-years. MTAs that were consistent across multiple trials were identified. Chr5LG3_566189651 and Chr5LG3_572899434 for plant height, Chr2LG1_409403647 for lodging resistance, Chr1LG6_57305683 and Chr1LG6_366513463 for grain yield, Chr1LG6_176606388, Chr2LG1_457185, Chr3LG5_234519042 and Chr7LG7_8229439 for seed starch concentration, and Chr3LG5_194530376 for seed protein concentration were identified from different locations and years. This research identified SNP markers associated with important traits in pea that have potential for marker-assisted selection towards rapid cultivar improvement.

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.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization

Background

Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies.

Results

By applying a two-step design process, which consisted of the construction of a first ‘draft’ array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping.

Conclusions

The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5712-3) contains supplementary material, which is available to authorized users.

Triarabinosylation is required for nodulation-suppressive CLE peptides to systemically inhibit nodulation in Pisum sativum

Legumes form root nodules to house beneficial nitrogen-fixing rhizobia bacteria. However, nodulation is resource demanding; hence, legumes evolved a systemic signalling mechanism called autoregulation of nodulation (AON) to control nodule numbers. AON begins with the production of CLE peptides in the root, which are predicted to be glycosylated, transported to the shoot, and perceived. We synthesized variants of nodulation-suppressing CLE peptides to test their activity using petiole feeding to introduce CLE peptides into the shoot. Hydroxylated, monoarabinosylated, and triarabinosylated variants of soybean GmRIC1a and GmRIC2a were chemically synthesized and fed into recipient Pisum sativum (pea) plants, which were used due to the availability of key AON pathway mutants unavailable in soybean. Triarabinosylated GmRIC1a and GmRIC2a suppressed nodulation of wild-type pea, whereas no other peptide variant tested had this ability. Suppression also occurred in the supernodulating hydroxyproline O-arabinosyltransferase mutant, Psnod3, but not in the supernodulating receptor mutants, Pssym29, and to some extent, Pssym28. During our study, bioinformatic resources for pea became available and our analyses identified 40 CLE peptide-encoding genes, including orthologues of nodulation-suppressive CLE peptides. Collectively, we demonstrated that soybean nodulation-suppressive CLE peptides can function interspecifically in the AON pathway of pea and require arabinosylation for their activity.

Unleashing meiotic crossovers in crops

Improved plant varieties are important in our attempts to face the challenges of a growing human population and limited planet resources. Plant breeding relies on meiotic crossovers to combine favourable alleles into elite varieties¹. However, meiotic crossovers are relatively rare, typically one to three per chromosome², limiting the efficiency of the breeding process and related activities such as genetic mapping. Several genes that limit meiotic recombination were identified in the model species Arabidopsis thaliana². Mutation of these genes in Arabidopsis induces a large increase in crossover frequency. However, it remained to be demonstrated whether crossovers could also be increased in crop species hybrids. We explored the effects of mutating the orthologues of FANCM³, RECQ4⁴ or FIGL1⁵ on recombination in three distant crop species, rice (Oryza sativa), pea (Pisum sativum) and tomato (Solanum lycopersicum). We found that the single recq4 mutation increases crossovers about three-fold in these crops, suggesting that manipulating RECQ4 may be a universal tool for increasing recombination in plants. Enhanced recombination could be used with other state-of-the-art technologies such as genomic selection, genome editing or speed breeding⁶ to enhance the pace and efficiency of plant improvement.

Unleashing meiotic crossovers in crops

Improved plant varieties are important in our attempts to face the challenges of a growing human population and limited planet resources. Plant breeding relies on meiotic crossovers to combine favourable alleles into elite varieties1. However, meiotic crossovers are relatively rare, typically one to three per chromosome2, limiting the efficiency of the breeding process and related activities such as genetic mapping. Several genes that limit meiotic recombination were identified in the model species Arabidopsis thaliana2. Mutation of these genes in Arabidopsis induces a large increase in crossover frequency. However, it remained to be demonstrated whether crossovers could also be increased in crop species hybrids. We explored the effects of mutating the orthologues of FANCM3, RECQ44 or FIGL15 on recombination in three distant crop species, rice (Oryza sativa), pea (Pisum sativum) and tomato (Solanum lycopersicum). We found that the single recq4 mutation increases crossovers about three-fold in these crops, suggesting that manipulating RECQ4 may be a universal tool for increasing recombination in plants. Enhanced recombination could be used with other state-of-the-art technologies such as genomic selection, genome editing or speed breeding6 to enhance the pace and efficiency of plant improvement.

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.

Construction of high-density linkage maps for mapping quantitative trait loci for multiple traits in field pea (Pisum sativum L.)

BACKGROUND

The objective of this research was to map quantitative trait loci (QTLs) of multiple traits of breeding importance in pea (Pisum sativum L.). Three recombinant inbred line (RIL) populations, PR-02 (Orb x CDC Striker), PR-07 (Carerra x CDC Striker) and PR-15 (1-2347-144 x CDC Meadow) were phenotyped for agronomic and seed quality traits under field conditions over multiple environments in Saskatchewan, Canada. The mapping populations were genotyped using genotyping-by-sequencing (GBS) method for simultaneous single nucleotide polymorphism (SNP) discovery and construction of high-density linkage maps.

RESULTS

After filtering for read depth, segregation distortion, and missing values, 2234, 3389 and 3541 single nucleotide polymorphism (SNP) markers identified by GBS in PR-02, PR-07 and PR-15, respectively, were used for construction of genetic linkage maps. Genetic linkage groups were assigned by anchoring to SNP markers previously positioned on these linkage maps. PR-02, PR-07 and PR-15 genetic maps represented 527, 675 and 609 non-redundant loci, and cover map distances of 951.9, 1008.8 and 914.2 cM, respectively. Based on phenotyping of the three mapping populations in multiple environments, 375 QTLs were identified for important traits including days to flowering, days to maturity, lodging resistance, Mycosphaerella blight resistance, seed weight, grain yield, acid and neutral detergent fiber concentration, seed starch concentration, seed shape, seed dimpling, and concentration of seed iron, selenium and zinc. Of all the QTLs identified, the most significant in terms of explained percentage of maximum phenotypic variance (PVmax) and occurrence in multiple environments were the QTLs for days to flowering (PVmax = 47.9%), plant height (PVmax = 65.1%), lodging resistance (PVmax = 35.3%), grain yield (PVmax = 54.2%), seed iron concentration (PVmax = 27.4%), and seed zinc concentration (PVmax = 43.2%).

CONCLUSION

We have identified highly significant and reproducible QTLs for several agronomic and seed quality traits of breeding importance in pea. The QTLs identified will be the basis for fine mapping candidate genes, while some of the markers linked to the highly significant QTLs are useful for immediate breeding applications.

Accelerating genetic gains in legumes for the development of prosperous smallholder agriculture: integrating genomics, phenotyping, systems modelling and agronomy

Grain legumes form an important component of the human diet, provide feed for livestock, and replenish soil fertility through biological nitrogen fixation. Globally, the demand for food legumes is increasing as they complement cereals in protein requirements and possess a high percentage of digestible protein. Climate change has enhanced the frequency and intensity of drought stress, posing serious production constraints, especially in rainfed regions where most legumes are produced. Genetic improvement of legumes, like other crops, is mostly based on pedigree and performance-based selection over the past half century. To achieve faster genetic gains in legumes in rainfed conditions, this review proposes the integration of modern genomics approaches, high throughput phenomics, and simulation modelling in support of crop improvement that leads to improved varieties that perform with appropriate agronomy. Selection intensity, generation interval, and improved operational efficiencies in breeding are expected to further enhance the genetic gain in experimental plots. Improved seed access to farmers, combined with appropriate agronomic packages in farmers' fields, will deliver higher genetic gains. Enhanced genetic gains, including not only productivity but also nutritional and market traits, will increase the profitability of farming and the availability of affordable nutritious food especially in developing countries.

Co-inoculation of a Pea Core-Collection with Diverse Rhizobial Strains Shows Competitiveness for Nodulation and Efficiency of Nitrogen Fixation Are Distinct traits in the Interaction

Pea forms symbiotic nodules with Rhizobium leguminosarum sv. viciae (Rlv). In the field, pea roots can be exposed to multiple compatible Rlv strains. Little is known about the mechanisms underlying the competitiveness for nodulation of Rlv strains and the ability of pea to choose between diverse compatible Rlv strains. The variability of pea-Rlv partner choice was investigated by co-inoculation with a mixture of five diverse Rlv strains of a 104-pea collection representative of the variability encountered in the genus Pisum. The nitrogen fixation efficiency conferred by each strain was determined in additional mono-inoculation experiments on a subset of 18 pea lines displaying contrasted Rlv choice. Differences in Rlv choice were observed within the pea collection according to their genetic or geographical diversities. The competitiveness for nodulation of a given pea-Rlv association evaluated in the multi-inoculated experiment was poorly correlated with its nitrogen fixation efficiency determined in mono-inoculation. Both plant and bacterial genetic determinants contribute to pea-Rlv partner choice. No evidence was found for co-selection of competitiveness for nodulation and nitrogen fixation efficiency. Plant and inoculant for an improved symbiotic association in the field must be selected not only on nitrogen fixation efficiency but also for competitiveness for nodulation.

Comparative Genome-Wide-Association Mapping Identifies Common Loci Controlling Root System Architecture and Resistance to Aphanomyces euteiches in Pea

Combining plant genetic resistance with architectural traits that are unfavorable to disease development is a promising strategy for reducing epidemics. However, few studies have identified root system architecture (RSA) traits with the potential to limit root disease development. Pea is a major cultivated legume worldwide and has a wide level of natural genetic variability for plant architecture. The root pathogen Aphanomyces euteiches is a major limiting factor of pea crop yield. This study aimed to increase the knowledge on the diversity of loci and candidate genes controlling RSA traits in pea and identify RSA genetic loci associated with resistance to A. euteiches which could be combined with resistance QTL in breeding. A comparative genome wide association (GWA) study of plant architecture and resistance to A. euteiches was conducted at the young plant stage in a collection of 266 pea lines contrasted for both traits. The collection was genotyped using 14,157 SNP markers from recent pea genomic resources. It was phenotyped for ten root, shoot and overall plant architecture traits, as well as three disease resistance traits in controlled conditions, using image analysis. We identified a total of 75 short-size genomic intervals significantly associated with plant architecture and overlapping with 46 previously detected QTL. The major consistent intervals included plant shoot architecture or flowering genes (PsLE, PsTFL1) with putative pleiotropic effects on root architecture. A total of 11 genomic intervals were significantly associated with resistance to A. euteiches confirming several consistent previously identified major QTL. One significant SNP, mapped to the major QTL Ae-Ps7.6, was associated with both resistance and RSA traits. At this marker, the resistance-enhancing allele was associated with an increased total root projected area, in accordance with the correlation observed between resistance and larger root systems in the collection. Seven additional intervals associated with plant architecture overlapped with GWA intervals previously identified for resistance to A. euteiches. This study provides innovative results about genetic interdependency of root disease resistance and RSA inheritance. It identifies pea lines, QTL, closely-linked markers and candidate genes for marker-assisted-selection of RSA loci to reduce Aphanomyces root rot severity in future pea varieties.

Revolution in Genotyping Platforms for Crop Improvement

In the past decade, the application of high-throughput sequencing to crop genotyping has given rise to novel platforms capable of genotyping tens of thousands of genome-wide DNA markers. Coupled with the decreasing costs of sequencing, this rapid increase in markers allows accelerated and highly accurate genotyping of entire crop populations and diversity sets using single nucleotide polymorphisms (SNPs). These revolutionary advances accelerate crop improvement by facilitating a more precise connection of phenotype to genotype through association studies, linkage mapping and diversity analysis. The platforms driving the advances in genotyping are array technologies and genotyping by sequencing (GBS) methods, which include both low-coverage whole genome resequencing (skim sequencing) and reduced representation sequencing (RRS) approaches. Here, we outline and compare these genotyping platforms and provide a perspective on the promising future of crop genotyping. While SNP arrays provide high quality, simple handling, and unchallenging analysis, the lower cost of RRS and the greater data volume produced by skim sequencing suggest that use of GBS will become more prevalent in crop genomics as sequencing costs decrease and data analysis becomes more streamlined. Graphical Abstract.

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.