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

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

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

[This corrects the article .].

Multi-trait genomic prediction improves selection accuracy for enhancing seed mineral concentrations in pea.

Multi-trait genomic selection (MT-GS) has the potential to improve predictive ability by maximizing the use of information across related genotypes and genetically correlated traits. In this study, we extended the use of sparse phenotyping method into the MT-GS framework by split testing of entries to maximize borrowing of information across genotypes and predict missing phenotypes for targeted traits without additional phenotyping expenditure. Using 300 advanced breeding lines from North Dakota State University (NDSU) pulse breeding program and ∼200 USDA accessions that were evaluated for 10 nutritional traits, our results show that the proposed sparse phenotyping aided MT-GS can further improve predictive ability by >12% across traits compared with univariate (UNI) genomic selection. The proposed strategy departed from the previous reports that weak genetic correlation is a limitation to the advantage of MT-GS over UNI genomic selection, which was evident in the partially balanced phenotyping-enabled MT-GS. Our results point to heritability and genetic correlation between traits as possible metrics to optimize and further improve the estimation of model parameters, and ultimately, prediction performance. Overall, our study offers a new approach to optimize the prediction performance using the MT-GS and further highlight strategy to maximize the efficiency of GS in a plant breeding program. The sparse-testing-aided MT-GS proposed in this study can be further extended to multi-environment, multi-trait GS to improve prediction performance and further reduce the cost of phenotyping and time-consuming data collection process.

We extended the use of sparse phenotyping into the multi-trait genomic selection (MT-GS) framework by split testing of entries. The sparse-phenotyping-aided MT-GS can increase predictive ability by >12% across traits. Heritability and genetic correlation are possible metrics to optimize and further improve prediction performance of MT-GS. The sparse-testing-aided MT-GS can be further extended to multi-environment, multi-trait GS framework.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

A reference genome for pea provides insight into legume genome evolution.

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

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

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

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

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

Molecular Evidence for Two Domestication Events in the Pea Crop.

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

Author Correction: Genomic diversity and macroecology of the crop wild relatives of domesticated pea.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

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

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

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

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

Selecting Lentil Accessions for Global Selenium Biofortification.

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