The USDA cucumber (Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development

Abstract

Germplasm collections are a crucial resource to conserve natural genetic diversity and provide a source of novel traits essential for sustained crop improvement. Optimal collection, preservation and utilization of these materials depends upon knowledge of the genetic variation present within the collection. Here we use the high-throughput genotyping-by-sequencing (GBS) technology to characterize the United States National Plant Germplasm System (NPGS) collection of cucumber (Cucumis sativus L.). The GBS data, derived from 1234 cucumber accessions, provided more than 23 K high-quality single-nucleotide polymorphisms (SNPs) that are well distributed at high density in the genome (~1 SNP/10.6 kb). The SNP markers were used to characterize genetic diversity, population structure, phylogenetic relationships, linkage disequilibrium, and population differentiation of the NPGS cucumber collection. These results, providing detailed genetic analysis of the U.S. cucumber collection, complement NPGS descriptive information regarding geographic origin and phenotypic characterization. We also identified genome regions significantly associated with 13 horticulturally important traits through genome-wide association studies (GWAS). Finally, we developed a molecularly informed, publicly accessible core collection of 395 accessions that represents at least 96% of the genetic variation present in the NPGS. Collectively, the information obtained from the GBS data enabled deep insight into the diversity present and genetic relationships among accessions within the collection, and will provide a valuable resource for genetic analyses, gene discovery, crop improvement, and germplasm preservation.

References

[1] Tyagi, R. K. & Agrawal, A. Revised genebank standards for management of plant genetic resources. Indian J. Agric. Sci. 85, 157-165 (2015).

[2] Davies, L. R. & Allender, C. J. Who is sowing our seeds? A systematic review of the use of plant genetic resources in research. Genet. Resour. Crop Evol. 64, 1999-2008 (2017).

[3] Fu, Y. B. The vulnerability of plant genetic resources conserved ex situ. Crop Sci. 57, 2314-2328 (2017).

[4] Grumet, R., Garcia-Mas, J. & Katzir, N. Cucurbit genetics and genomics: a look to the future. In Genetics and Genomics of the Cucurbitaceae (eds Grumet, R., Katzir, N., Garcia-Mas, J. ). pp 409-415 (Springer, 2017).

[5] Renner, S. S., Schaefer, H. & Kocyan, A. Phylogenetics of Cucumis (Cucurbitaceae): Cucumber (C. sativus) belongs in an Asian/Australian clade far from melon (C. melo) . BMC Evol. Biol. 7, 58 (2007).

[6] Sebastian, P. M., Schaefer, H., Telford, I. R. H. & Renner, S. S. Cucumber and melon have their wild progenitors in India, and the sister species of Cucumis melo is from Australia . Proc. Nat. Acad. Sci. USA 107, 14269-14273 (2010).

[7] McCreight, J. D., Staub, J. E., Wehner, T. C. & Dhillon, N. P. S. Gone global: familiar and exotic cucurbits have Asian origins. HortScience 48, 1078-1089 (2013).

[8] Naegele, R. P. & Wehner, T. C. Genetic resources of cucumber. In Genetics and genomics of the Cucurbitacae (eds Grumet, R., Katzir, N., Garcia-Mas, J.). pp 61-86 (Springer, 2016).

[9] Lv, J. et al. Genetic diversity and population structure of cucumber (Cucumis sativus L.) . PLoS One 7, e46919 (2012).

[10] Qi, J. et al. A genomic variation map provides insights into the genetic basis of cucumber domestication and diversity. Nat. Genet. 45, 1510-1515 (2013).

[11] Yang, L. M. et al. Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly. Plant J. 71, 895-906 (2012).

[12] Kirkbride, J. H. Biosystematic monograph of the genus Cucumis (Cucurbitaceae) . (Parkway Publishers, Boone, NC, 1993).

[13] Staub, J. E., Serquen, F. C., Horejsi, T. & Chen, J. Genetic diversity in cucumber (Cucumis sativus L.): IV. An evaluation of Chinese germplasm . Genet. Resour. Crop Evol. 46, 297-310 (1999).

[14] Simon, P. W. & Navazio, J. P. Early orange mass 400, early orange mass 402, and late orange mass 404: High-carotene cucumber germplasm. HortScience 32, 144-145 (1997).

[15] Bo, K. L., Ma, Z., Chen, J. F. & Weng, Y. Molecular mapping reveals structural rearrangements and quantitative trait loci underlying traits with local adaptation in semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis Qi et Yuan) . Theor. Appl. Genet. 128, 25-39 (2014).

[16] Pan, Y. et al. QTL mapping of domestication and diversifying selection related traits in round-fruited semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis) . Theor. Appl. Genet. 130, 1531-1548 (2017).

[17] Call, A. D., Criswell, A. D., Wehner, T. C., Ando, K. & Grumet, R. Resistance of current cucumber cultivars to a new strain of cucurbit downy mildew. HortScience 47, 171-178 (2012).

[18] Block, C. & Reitsma, K. R. Powdery mildew resistance in the U.S. national plant germplasm system cucumber collection. HortScience 40, 416-420 (2005).

[19] Colle, M., Straley, E. N., Makela, S. B., Hammar, S. A. & Grumet, R. Screening the cucumber plant introduction collection for young fruit resistance to Phytophthora capsici . HortScience 49, 244-249 (2014).

[20] Gevens, A. J., Ando, K., Lamour, K. H., Grumet, R. & Hausbeck, M. K. Development of a detached cucumber fruit assay to screen for resistance and effect of fruit age on susceptibility to infection by Phytophthora capsici . Plant Dis. 90, 1276-1282 (2006).

[21] Uchneat, M. S. & Wehner, T. C. Resistance to belly rot in cucumber identified through field and detached-fruit evaluations. J. Am. Soc. Hort. Sci. 123, 78-84 (1998).

[22] Wehner, T. C., Shetty, N. V. & Sloane, J. T. Field and detached-fruit screening tests for resistance to belly rot in cucumber. HortScience 38, 149-152 (2004).

[23] Walters, S. A., Wehner, T. C. & Barker, K. R. NC-42 and NC-43: root-knot nematode-resistant cucumber germplasm. HortScience 31, 1246-1247 (1996).

[24] Shetty, N. V. & Wehner, T. C. Screening the cucumber germplasm collection for fruit yield and quality. Crop Sci. 42, 2174-2183 (2001).

[25] Ando, K. & Grumet, R. Evaluation of altered plant architecture as a means to reduce Phytophthora capsici disease incidence on cucumber fruit . J. Am. Soc. Hort. Sci. 131, 491-498 (2006).

[26] Grumet, R., Barczak, M., Tabaka, C. & Duvall, R. Above-ground screening for genotypic differences in cucumber root growth in the greenhouse and field. J. Am. Soc. Hort. Sci. 117, 1006-1011 (1992).

[27] Meglic, V., Serquen, F. & Staub, J. E. Genetic diversity in cucumber (Cucumis sativus L.): I. A reevaluation of the U.S. germplasm collection . Genet. Res. Crop Evol. 43, 533-546 (1996).

[28] Horejsi, T. & Staub, J. E. Genetic variation in cucumber (Cucumis sativus L.) as assessed by random amplified polymorphic DNA . Genet. Resour. Crop Evol. 46, 337-350 (1999).

[29] Jia, J. Z., Li, H. J., Zhang, X. Y., Li, Z. C. & Qiu, L. J. Genomics-based plant germplasm research (GPGR). Crop J. 5, 166-174 (2017).

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

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

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

[33] Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156-2158 (2011).

[34] Lee, T. H., Guo, H., Wang, X., Kim, C. & Paterson, A. H. SNPhylo, a pipeline to construct a phylogenetic tree from huge SNP data. BMC Genom. 15, 162 (2014).

[35] Yu, G. et al. ggtree, an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 8, 28-36 (2017).

[36] Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559-575 (2007).

[37] Falush, D., Stephens, M. & Pritchard, J. K. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164, 1567-1587 (2003).

[38] Stajich, J. E. et al. The bioperl toolkit: perl modules for the life sciences. Genome Res. 12, 1611-1618 (2002).

[39] Huang, X. et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 42, 961-967 (2010).

[40] Lippert, C. et al. FaST linear mixed models for genome-wide association studies. Nat. Methods 8, 833-835 (2011).

[41] Li, M. X., Yeung, J. M., Cherny, S. S. & Sham, P. C. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum. Genet. 131, 747-756 (2012).

[42] Jeong, S. et al. GenoCore: a simple and fast algorithm for core subset selection from large genotype datasets. PLoS One 12, e0181420 (2017).

[43] Whitaker, T. W. & Davis, G. N. Cucurbits: Botany, cultivation, and utilization. (Interscience Publishers, New York, 1962).

[44] Staub, J. E., Serquen, F. C. & McCreight, J. D. Genetic diversity in cucumber (Cucumis sativus L.): III. An evaluation of Indian germplasm . Genet. Resour. Crop Evol. 44, 315-326 (1997).

[45] Staub, J. E., Dane, F., Reitsma, K., Fazio, G. & Lopez-Sese, A. The formation of test arrays and a core collection in cucumber using phenotypic and molecular marker data. J. Am. Soc. Hortic. Sci. 127, 558-567 (2002).

[46] Walters, T. W. Historical overview on domesticated plants in China with special emphasis on the Cucurbitaceae. Econ. Bot. 43, 297-313 (1989).

[47] Paris, H. S., Daunay, M. C. & Janick, J. Occidental diffusion of cucumber (Cucumis sativus) 500-1300 CE: two routes to Europe . Ann. Bot. 109, 117-126 (2012).

[48] Pandey, S., Ansari, W. A., Mishra, V. K., Singh, A. K. & Singh, M. Genetic diversity in Indian cucumber based on microsatellite and morphological markers. Biochem. Syst. Ecol. 51, 19-27 (2013).

[49] Mansfeld, B. & Grumet, R. Metabolomic plasticity of cucumber fruit peel - effects of developmental stage and market class. In Proceedings of Cucurbitaceae 2016 (eds Kozik, E. U., Paris, H. S.). pp 174-179 (2016).

[50] Wang, Y. et al. QTL mapping for downy mildew resistance in cucumber inbred line WI7120 (PI 330628). Theor. Appl. Genet. 129, 1493-1505 (2016).

[51] Wang, Y., VandenLangenberg, K., Wen, C., Wehner, T. C. & Weng, Y. QTL mapping of downy and powdery mildew resistances in PI197088 cucumber with genotyping-by-sequencing in RIL population. Theor. Appl. Genet. 131, 597-611 (2018).

[52] Pan, J. et al. STAYGREEN (CsSGR) is a candidate for the anthracnose (Colletotrichum orbiculare) resistance locus cla in Gy14 cucumber . Theor. Appl. Genet. 131, 1577-1587 (2018).

[53] Zhang, S. P. et al. Inheritance and QTL mapping of resistance to gummy stem blight in cucumber stem. Mol. Breed. 37, 49 (2017).

[54] Wang, Y. H. et al. STAYGREEN, STAY HEALTHY: a loss-of-susceptibility mutation in the STAYGREEN gene provides durable, broad-spectrum disease resistances for over 50 years of US cucumber production. New Phytol. (2018). https://doi.org/10.1111/nph.15353.

[55] Liu, S. L. et al. Genetic analysis and QTL mapping of resistance to gummy stem blight in Cucumis sativus seedling stage . Plant Dis. 101, 1145-1152 (2017).

[56] Fazio, G., Staub, J. E. & Stevens, M. R. Genetic mapping and QTL analysis of horticultural traits in cucumber (Cucumis sativus L.) using recombinant inbred lines . Theor. Appl. Genet. 107, 864-874 (2003).

[57] Hu, Z. et al. A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo . DNA Res. 25, 1-10 (2018).

[58] Yang, L. M. et al. Next-generation sequencing, FISH mapping, and synteny-based modeling reveal mechanisms of dysploid chromosome reduction in Cucumis . Plant J. 77, 16-30 (2004).

[59] Brown, A. H. D. Core collections: a practical approach to genetic resources management. Genome 31, 818-824 (1989).