Genome-wide SNP discovery and core marker sets for assessment of genetic variations in cultivated pumpkin (Cucurbita spp.)

Abstract

Three pumpkin species Cucurbita maxima, C. moschata, and C. pepo are commonly cultivated worldwide. To identify genome-wide SNPs in these cultivated pumpkin species, we collected 48 F₁ cultivars consisting of 40 intraspecific hybrids (15 C. maxima, 18 C. moschata, and 7 C. pepo) and 8 interspecific hybrids (C. maxima x C. moschata). Genotyping by sequencing identified a total of 37,869 confident SNPs in this collection. These SNPs were filtered to generate a subset of 400 SNPs based on polymorphism and genome distribution. Of the 400 SNPs, 288 were used to genotype an additional 188 accessions (94 F₁ cultivars, 50 breeding lines, and 44 landraces) with a SNP array-based platform. Reliable polymorphisms were observed in 224 SNPs (78.0%) and were used to assess genetic variations between and within the four predefined populations in 223 cultivated pumpkin accessions. Both principal component analysis and UPGMA clustering found four major clusters representing three pumpkin species and interspecific hybrids. This genetic differentiation was supported by pairwise F_st and Nei’s genetic distance. The interspecific hybrids showed a higher level of genetic diversity relative to the other three populations. Of the 224 SNPs, five subsets of 192, 96, 48, 24, and 12 markers were evaluated for variety identification. The 192, 96, and 48 marker sets identified 204 (91.5%), 190 (85.2%), and 141 (63.2%) of the 223 accessions, respectively, while other subsets showed <25% of variety identification rates. These SNP markers provide a molecular tool with many applications for genetics and breeding in cultivated pumpkin.

References

[1] Sanjur, O. I., Piperno, D. R., Andres, T. C. & Wessel-Beaver, L. Phylogenetic relationships among domesticated and wild species of Cucurbita (Cucurbitaceae) inferred from a mitochondrial gene: Implications for crop plant evolution and areas of origin. Proc. Natl. Acad. Sci. USA 99, 535-540 (2002).

[2] Whitaker, T. W. & Bemis, W. P. Evolution in the genus Cucurbita. Evolution 18, 553-559 (1964).

[3] Aamir, H. D., Sofi, S. A. & Rafiq, S. Pumpkin the functional and therapeutic ingredient: a review. Int. J. Food Sci. Nutr. 2, 165-170 (2017).

[4] Zhou, T., Kong, Q., Huang, J., Dai, R. & Li, Q. Characterization of nutritional components and utilization of pumpkin. Food 1, 313-321 (2007).

[5] Food & Agriculture Organization of the United, N. FAOSTAT statistical database. https://search.library.wisc.edu/catalog/999890171702121 (2020).

[6] Jones, H. et al. Evaluation of the use of high-density SNP genotyping to implement UPOV Model 2 for DUS testing in barley. Theor. Appl. Genet. 126, 901-911 (2013).

[7] Jamali, S. H., Cockram, J. & Hickey, L. T. Insights into deployment of DNA markers in plant variety protection and registration. Theor. Appl. Genet. 132, 1911-1929 (2019).

[8] UPOV, International Union for the Protection of New Varieties of Plants . Combining morphological and molecular distance in the management of the reference collection of potato. BMT/13/10, https://www.upov.int/edocs/mdocs/upov/en/bmt_13/bmt_13_10.pdf (2011).

[9] UPOV, International Union for the Protection of New Varieties of Plants . Guidance on the use of biochemical and molecular markers in the examination of distinctness, uniformity and stability (DUS). TGP/15, https://www.upov.int/edocs/tgpdocs/en/tgp_15.pdf (2013).

[10] Jones, H. & Mackay, I. Implications of using genomic prediction within a high-density SNP dataset to predict DUS traits in barley. Theor. Appl. Genet. 128, 2461-2470 (2015).

[11] Bae, K. M. et al. Development of genomic SSR markers and genetic diversity analysis in cultivated radish (Raphanus sativus L.). Hortic. Environ. Biotechnol. 56, 216-224 (2015).

[12] Choi, S. P. et al. Genetic characterisation of commercial Chinese cabbage varieties using SSR markers. Seed Sci. Technol. 44, 595-608 (2016).

[13] Jamali, S., Sadeghi, L. & Sadeghin-Motahhar, S. Identification and distinction of soybean commercial cultivars using morphological and microsatellite markers. Iran. J. Crop. Sci. 13, 131-145 (2011).

[14] Kwon, Y. S. & Choi, K. J. Construction of a DNA profile database for commercial cucumber (Cucumis sativus L.) cultivars using microsatellite marker. Korean J. Hortic. Sci. Technol. 31, 344-351 (2013).

[15] Phan, N. T., Kim, M.-K. & Sim, S.-C. Genetic variations of F1 tomato cultivars revealed by a core set of SSR and InDel markers. Sci. Hortic. 212, 155-161 (2016).

[16] Sim, S. C., Hong, J. H. & Kwon, Y. S. DNA profiling of commercial pumpkin cultivars using simple sequence repeat polymorphisms. Hortic. Environ. Biotechnol. 56, 811-820 (2015).

[17] Kong, Q. S., Liu, Y., Xie, J. J. & Bie, Z. L. Development of simple sequence repeat markers from de novo assembled transcriptomes of pumpkins. Plant Mol. Biol. Rep. 38, 130-136 (2020).

[18] Hao, N. et al. Genome-based breeding approaches in major vegetable crops. Theor. Appl. Genet. https://doi.org/10.1007/s00122-019-03477-z (2019).

[19] Blanca, J. et al. Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae). BMC Genomics 12, 104 (2011).

[20] Zhang, G. et al. A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.). BMC Genomics 16, 1101 (2015).

[21] Holdsworth, W. L., LaPlant, K. E., Bell, D. C., Jahn, M. M. & Mazourek, M. Cultivar-based introgression mapping reveals wild species-derived Pm-0, the major powdery mildew resistance locus in squash. PLoS ONE 11, e0167715 (2016).

[22] Montero-Pau, J. et al. An SNP-based saturated genetic map and QTL analysis of fruit-related traits in Zucchini using Genotyping-by-sequencing. BMC Genomics 18, 94 (2017).

[23] Montero-Pau, J. et al. De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol. J. 16, 1161-1171 (2018).

[24] Sun, H. et al. Karyotype stability and unbiased fractionation in the paleo-allotetraploid Cucurbita genomes. Mol. Plant. 10, 1293-1306 (2017).

[25] Xanthopoulou, A. et al. Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits. Hortic. Res. 6, 94 (2019).

[26] Elshire, R. J. et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6, e19379 (2011).

[27] Kishor, D. S. et al. Development of SNP markers and validation assays in commercial Korean melon cultivars, using Genotyping-by-sequencing and Fluidigm analyses. Sci. Hortic. 263, 109113 (2019).

[28] Phan, N. T. et al. Identification of loci associated with fruit traits using genome-wide single nucleotide polymorphisms in a core collection of tomato (Solanum lycopersicum L.). Sci. Hortic. 243, 567-574 (2019).

[29] Spindel, J. et al. Bridging the genotyping gap: using genotyping by sequencing (GBS) to add high-density SNP markers and new value to traditional bi-parental mapping and breeding populations. Theor. Appl. Genet. 126, 2699-2716 (2013).

[30] Gong, L. et al. Genetic relationships and evolution in Cucurbita as viewed with simple sequence repeat polymorphisms: the centrality of C. okeechobeensis. Genet. Resour. Crop Evol. 60, 1531-1546 (2013).

[31] Kong, Q. S. et al. Genetic diversity of Cucurbita rootstock germplasm as assessed using simple sequence repeat markers. Sci. Hortic. 175, 150-155 (2014).

[32] Ferriol, M., Pico, B. & Nuez, F. Genetic diversity of a germplasm collection of Cucurbita pepo using SRAP and AFLP markers. Theor. Appl. Genet. 107, 271-282 (2003).

[33] Paris, H. A proposed subspecific classification for Cucurbita pepo. Phytologia 61, 133-138 (1986).

[34] Kabelka, E., Franchino, B. & Francis, D. M. Two loci from Lycopersicon hirsutum LA407 confer resistance to strains of Clavibacter michiganensis subsp. michiganensis. Phytopathology 92, 504-510 (2002).

[35] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[36] Glaubitz, J. C. et al. TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline. PLoS ONE 9, e90346 (2014).

[37] Botstein, D., White, R. L., Skolnick, M. & Davis, R. W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 32, 314-331 (1980).

[38] Stacklies, W., Redestig, H., Scholz, M., Walther, D. & Selbig, J. pcaMethods-a bioconductor package providing PCA methods for incomplete data. Bioinformatics 23, 1164-1167 (2007).

[39] R: A language and environment for statistical computing (R Foundation for Statistical Computing Vienna, Austria, 2010).

[40] Galili, T. dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering. Bioinformatics 31, 3718-3720 (2015).

[41] Weir, B. S. & Cockerham, C. C. Estimating F-statistics for the analysis of population structure. Evolution 38, 1358-1370 (1984).

[42] Dieringer, D. & Schlotterer, C. MICROSATELLITE ANALYSER (MSA): a platform independent analysis tool for large microsatellite data sets. Mol. Ecol. Notes 3, 167-169 (2003).

[43] Hurlbert, S. H. The nonconcept of species diversity: a critique and alternative parameters. Ecology 52, 577-586 (1971).

[44] El Mousadik, A. & Petit, R. High level of genetic differentiation for allelic richness among populations of the argan tree [Argania spinosa (L.) Skeels] endemic to Morocco. Theor. Appl. Genet. 92, 832-839 (1996).

[45] Nei, M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89, 583-590 (1978).