Merging genotyping-by-sequencing data from two ex situ collections provides insights on the pea evolutionary history

Abstract

Pea (Pisum sativum L. subsp. sativum) is one of the oldest domesticated species and a widely cultivated legume. In this study, we combined next generation sequencing (NGS) data referring to two genotyping-by-sequencing (GBS) libraries, each one prepared from a different Pisum germplasm collection. The selection of single nucleotide polymorphism (SNP) loci called in both germplasm collections caused some loss of information; however, this did not prevent the obtainment of one of the largest datasets ever used to explore pea biodiversity, consisting of 652 accessions and 22 127 markers. The analysis of population structure reflected genetic variation based on geographic patterns and allowed the definition of a model for the expansion of pea cultivation from the domestication centre to other regions of the world. In genetically distinct populations, the average decay of linkage disequilibrium (LD) ranged from a few bases to hundreds of kilobases, thus indicating different evolutionary histories leading to their diversification. Genome-wide scans resulted in the identification of putative selective sweeps associated with domestication and breeding, including genes known to regulate shoot branching, cotyledon colour and resistance to lodging, and the correct mapping of two Mendelian genes. In addition to providing information of major interest for fundamental and applied research on pea, our work describes the first successful example of integration of different GBS datasets generated from ex situ collections – a process of potential interest for a variety of purposes, including conservation genetics, genome-wide association studies, and breeding.

References

[1] Sultani MI, Gill MA, Anwar MM, Athar M . Evaluation of soil physical properties as influenced by various green manuring legumes and phosphorus fertilization under rain fed conditions. Int J Environ Sci Technol. 2007; 4: 109-18.

[2] Annicchiarico P. Feed legumes for truly sustainable crop-animal systems. Ital J Agronomy. 2017; 12: 880.

[3] Smýkal P, Kenicer G, Flavell AJ et al. Phylogeny, phylogeography and genetic diversity of the Pisum genus. Plant Genet Resour. 2011; 9: 4-18.

[4] Trněný O, Brus J, Hradilová I et al. Molecular evidence for two domestication events in the pea crop. Genes (Basel). 2018; 9: 535.

[5] Zohary D, Hopf M . Domestication of pulses in the old world. Sci. 1973; 182: 887-94.

[6] Weeden NF . Domestication of pea (Pisum sativum L.): the case of the abyssinian pea. Front Plant Sci. 2018; 9: 515.

[7] Smýkal P et al. The case for common markers and data compatibility. Effort towards a world pea (Pisum sativum L.) germplasm core collection. Pisum Genetics. 2008; 40: 11-4.

[8] Holdsworth WL, Gazave E, Cheng P et al. A community resource for exploring and utilizing genetic diversity in the USDA pea single plant plus collection. Hortic Res. 2017; 4: 17017.

[9] Annicchiarico P, Romani M, Cabassi G et al. Diversity in a pea (Pisum sativum) world collection for key agronomic traits in a rain-fed environment of southern Europe. Euphytica. 2017; 213: 245.

[10] Elshire RJ, Glaubitz JC, Sun Q et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011; 6: e19379.

[11] Taranto F, Nicolia A, Pavan S et al. Biotechnological and digital revolution for climate-smart plant breeding. Agronomy. 2018; 8: 277.

[12] Singh N, Wu S, Raupp WJ et al. Efficient curation of genebanks using next generation sequencing reveals substantial duplication of germplasm accessions. Sci Rep. 2019; 9: 650.

[13] Kreplav J, Madoui MA, Cápal P et al. A reference genome for pea provides insight into legume genome evolution. Nat Genet. 2019; 51: 1411-22.

[14] Pavan S, Delvento C, Ricciardi L et al. Recommendations for choosing the genotyping method and best practices for quality control in crop genome-wide association studies. Front Genet. 2020; 11: 447.

[15] Jing R, Vershinin A, Grzebyta J et al. The genetic diversity and evolution of field pea (Pisum) studied by high throughput retrotransposon based insertion polymorphism (RBIP) marker analysis. BMC Evol Biol. 2010; 10: 44.

[16] Jing R, Ambrose MA, Knox MR et al. Genetic diversity in European Pisum germplasm collections. Theor Appl Genet. 2012; 125: 367-80.

[17] Rodríguez-Ramilo ST, Wang J . The effect of close relatives on unsupervised bayesian clustering algorithms in population genetic structure analysis. Mol Ecol Resour. 2012; 12: 873-84.

[18] Pritchard JK, Stephens M, Donnelly P . Inference of population structure using multilocus genotype data. Genetics. 2000; 155: 945-59.

[19] Evanno G, Regnaut S, Goudet J . Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005; 14: 2611-20.

[20] Beji S, Fontaine V, Devaux R et al. Genome-wide association study identifies favorable SNP alleles and candidate genes for frost tolerance in pea. BMC Genomics. 2020; 21: 536.

[21] Cousin R. Peas (Pisum sativum L.). Field Crops Res. 1997; 53: 111-30.

[22] Kosev V, Pachev I, Angelova S et al. Inheritance of quantitative traits in crosses between two Pisum sativum subspecies with particular reference to their breeding value. Genetika. 2012; 48: 41-6.

[23] Weeden NF . Genetic changes accompanying the domestication of Pisum sativum: is there a common genetic basis to the ’domestication syndrome’ for legumes? Ann Bot. 2007; 100: 1017-25.

[24] Bordat A, Savois V, Nicolas M et al. Translational genomics in legumes allowed placing in silico 5460 unigenes on the pea functional map and identified candidate genes in Pisum sativum L. G3 (Bethesda). 2011; 1: 93-103.

[25] Bhattacharyya MK, Smith AM, Ellis THN et al. The wrinkled-seed character of pea described by Mendel is caused by a transposon-like insertion in a gene encoding starch-branching enzyme. Cell. 1990; 60: 115-22.

[26] Hellens RP, Moreau C, Lin-Wang K et al. Identification of Mendel’s white flower character. PLoS One. 2010; 5: e13230.

[27] Hradilová I, Duchoslav M, Brus J et al. Variation in wild pea (Pisum sativum subsp. elatius) seed dormancy and its relationship to the environment and seed coat traits. PeerJ. 2019; 2019: e6263.

[28] Siol M, Jacquin F, Chabert-Martinello M et al. Patterns of genetic structure and linkage disequilibrium in a large collection of pea germplasm. G3 (Bethesda). 2017; 7: 2461-71.

[29] Lawson DJ, van Dorp L, Falush DN . A tutorial on how not to overinterpret STRUCTURE and ADMIXTURE bar plots. Nat Commun. 2018; 9: 3258.

[30] Young JPW, Matthews PA . Distinct class of peas (Pisum sativum L.) from Afghanistan that show strain specificity for symbiotic rhizobium. Heredity. 1982; 48: 203-10.

[31] Pavan S, Marcotrigiano AR, Ciani E et al. Genotyping-by-sequencing of a melon (Cucumis melo L.) germplasm collection from a secondary center of diversity highlights patterns of genetic variation and genomic features of different gene pools. BMC Genomics. 2017; 18: 59.

[32] Pavan S, Delvento C, Mazzeo R et al. Almond diversity and homozygosity define structure, kinship, inbreeding, and linkage disequilibrium in cultivated germplasm, and reveal genomic associations with nut and seed weight. Hortic Res. 2021; 8: 15.

[33] Stirnberg P, van de Sande K, Leyser HMO . MAX1 and MAX2 control shoot lateral branching in Arabidopsis. Development. 2002; 129: 1131-41.

[34] Smýkal P, Aubert G, Burstin J et al. Pea (Pisum sativum L.) in the genomic era. Agronomy. 2012; 2: 74-115.

[35] Hofer J, Turner L, Moreau C et al. Tendril-less regulates tendril formation in pea leaves. Plant Cell. 2009; 21: 420-8.

[36] Sato Y, Morita R, Nishimura M et al. Mendel’s green cotyledon gene encodes a positive regulator of the chlorophyll-degrading pathway. Proc Natl Acad Sci U S A. 2007; 104: 14169-74.

[37] Macas J, Neumann P, Navrátilová A . Repetitive DNA in the pea (Pisum sativum L.) genome: comprehensive characterization using 454 sequencing and comparison to soybean and Medicago truncatula. BMC Genomics. 2007; 8: 427.

[38] Ziarsolo P, Hasing T, Hilario R et al. K-seq, an affordable, reliable, and open Klenow NGS-based genotyping technology. Plant Methods. 2021; 17: 30.

[39] Schork AJ, Thompson WK, Pham P et al. All SNPs are not created equal: genome-wide association studies reveal a consistent pattern of enrichment among functionally annotated SNPs. PLoS Genet. 2013; 9: e1003449.

[40] Wallace JG, Bradbury PJ, Zhang N et al. Association mapping across numerous traits reveals patterns of functional variation in maize. PLoS Genet. 2014; 10: e1004845.

[41] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[42] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[43] Puritz JB, Hollenbeck CM, Gold JR . dDocent: a RADseq, variant-calling pipeline designed for population genomics of non-model organisms. PeerJ. 2014; 2: e431.

[44] Quinlan AR, Hall IM . BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26: 841-2.

[45] Garrison E, Marth G . Haplotype-Based Variant Detection from Short-Read Sequencing. Preprint at https://arxiv.org/abs/1207.3907 (2012).

[46] Bradbury PJ, Zhang Z, Kroon DE et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23: 2633-5.

[47] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[48] Untergasser A, Cutcutache I, Koressaar T et al. Primer3-new capabilities and interfaces. Nucleic Acids Res. 2012; 40: e115.

[49] Chhatre VE, Emerson KJ . StrAuto: automation and parallelization of STRUCTURE analysis. Bioinformatics. 2017; 18: 192.

[50] Earl DA, vonHoldt BM . STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012; 4: 359-61.

[51] Kumar S, Stecher G, Li M et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018; 35: 1547-9.

[52] Gao X, Starmer JD . AWclust: point-and-click software for non-parametric population structure analysis. BMC Bioinformatics. 2008; 9: 77.

[53] Kamvar ZN, Tabima JF, Grunwald NJ . Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ. 2014; 2: 1-14.

[54] Paradis E, Claude J, Strimmer K . APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004; 20: 289-90.

[55] Zheng X, Levine D, Shen J et al. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics. 2012; 28: 3326-8.

[56] Guangchuang Y . Scatterpie: Scatter Pie Plot. R package version 0.1.7. https://cran.r-project.org/web/packages/scatterpie/ (2021).

[57] Chen H, Patterson N, Reich D . Population differentiation as a test for selective sweeps. Genome Res. 2010; 20: 393-402.

[58] Kang HM, Sul JH, Service SK et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010; 42: 348-54.