Genome-wide core sets of SNP markers and Fluidigm assays for rapid and effective genotypic identification of Korean cultivars of lettuce (Lactuca sativa L.)

Abstract

Lettuce is one of the economically important leaf vegetables and is cultivated mainly in temperate climate areas. Cultivar identification based on the distinctness, uniformity, and stability (DUS) test is a prerequisite for new cultivar registration. However, DUS testing based on morphological features is time-consuming, labor-intensive, and costly, and can also be influenced by environmental factors. Thus, molecular markers have also been used for the identification of genetic diversity as an effective, accurate, and stable method. Currently, genome-wide single nucleotide polymorphisms (SNPs) using next-generation sequencing technology are commonly applied in genetic research on diverse plant species. This study aimed to establish an effective and high-throughput cultivar identification system for lettuce using core sets of SNP markers developed by genotyping by sequencing (GBS). GBS identified 17 877 high-quality SNPs for 90 commercial lettuce cultivars. Genetic differentiation analyses based on the selected SNPs classified the lettuce cultivars into three main groups. Core sets of 192, 96, 48, and 24 markers were further selected and validated using the Fluidigm platform. Phylogenetic analyses based on all core sets of SNPs successfully discriminated individual cultivars that have been currently recognized. These core sets of SNP markers will support the construction of a DNA database of lettuce that can be useful for cultivar identification and purity testing, as well as DUS testing in the plant variety protection system. Additionally, this work will facilitate genetic research to improve breeding in lettuce.

References

[1] Korean Statistical Information Service (KOSIS) . http://kosis.kr/ (2020).

[2] Reyes-Chin-Wo S, Wang Z, Yang X et al. Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce. Nat Commun. 2017; 8: 14953.

[3] Kwon S, Truco M, Hu J . LS germ OPA, a custom OPA of 384 EST-derived SNPs for high-throughput lettuce (Lactuca sativa L.) germplasm fingerprinting . Mol Breed. 2012; 29: 887-901.

[4] UPOV, International Union for the Protection of New Varieties of Plants . Lettuce Guidelines for the Conduct of Tests for Distinctness Uniformity and Stability. UPOV/TG/13/11 , https://www.upov.int/edocs/tgpdocs/en/tg013.pdf (2017).

[5] Jamali SH, Cockram J, Hickey LT . Insights into deployment of DNA markers in plant variety protection and registration. Theor Appl Genet. 2019; 132: 1911-29.

[6] 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).

[7] UPOV, International Union for the Protection of New Varieties of Plants . Guidelines for DNA -Profiling: Molecular Marker Selection and Database Construction (BMT Guidelines). UPOV/INF/17/1 , http://www.upov.int/edocs/infdocs/en/upov_inf_17_1.pdf (2010).

[8] Reid A, Hof L, Felix G et al. Construction of an integrated microsatellite and key morphological characteristic database of potato varieties on the EU common catalogue. Euphytica. 2011; 182: 239-49.

[9] Vélez MD, Ibáñez J . Assessment of the uniformity and stability of grapevine cultivars using a set of microsatellite markers. Euphytica. 2012; 186: 419-32.

[10] Zhou H, Zhang P, Luo J et al. The establishment of a DNA fingerprinting database for 73 varieties of Lactuca sativa capitate L. using SSR molecular markers. Hortic Environ Biotechnol. 2019; 60: 95-103.

[11] Zhang S, Li B, Chen Y et al. Molecular-assisted distinctness and uniformity testing using SLAF-sequencing approach in soybean. Genes. 2020; 11: 175.

[12] Zhang D, Vega FE, Infante F et al. Accurate differentiation of green beans of Arabica and Robusta coffee using nanofluidic array of single nucleotide polymorphism (SNP) markers. J AOAC Int. 2020; 103: 315-24.

[13] Kishor DS, Song WH, Noh Y et al. Development of SNP markers and validation assays in commercial Korean melon cultivars, using genotyping-by-sequencing and Fluidigm analyses. Sci Hortic. 2020; 263: 1-9.

[14] Nguyen NN, Kim M, Jung JK et al. Genome-wide SNP discovery and core marker sets for assessment of genetic variations in cultivated pumpkin (Cucurbita spp.) . Hortic Res. 2020; 7: 121.

[15] Kim M, Jung JK, Shim EJ et al. Genome-wide SNP discovery and core marker sets for DNA barcoding and variety identification in commercial tomato cultivars. Sci Hortic. 2021; 276: 1-7.

[16] Park G, Choi Y, Jung JK et al. Genetic diversity assessment and cultivar identification of cucumber (Cucumis sativus L.) using the Fluidigm single nucleotide polymorphism assay . Plan Theory. 2021; 10: 395.

[17] Riar DS, Rustgi S, Burke IC et al. EST-SSR development from 5 Lactuca species and their use in studying genetic diversity among L. serriola biotypes . J Hered. 2011; 102: 17-28.

[18] Simko I . Development of EST-SSR markers for the study of population structure in lettuce (Lactuca sativa L.) . J Hered. 2009; 100: 256-62.

[19] Uwimana B, D’Andrea L, Felber F et al. A Bayesian analysis of gene flow from crops to their wild relatives: cultivated (Lactuca sativa L.) and prickly lettuce (L. serriola L.) and the recent expansion of L. serriola in Europe. Mol Ecol. 2012; 21: 2640-54.

[20] Rauscher G, Simko I . Development of genomic SSR markers for fingerprinting lettuce (Lactuca sativa L.) cultivars and mapping genes . BMC Plant Biol. 2013; 13: 11-1.

[21] Hong JH, Kwon YS, Choi KJ et al. Identification of lettuce germplasms and commercial cultivars using SSR markers developed from EST. Korean J Hort Sci Technol. 2013; 31: 772-81.

[22] Sochor M, Jemelková M, Doležalová I . Phenotyping and SSR markers as a tool for identification of duplicates in lettuce germplasm. Czech J Genet Plant Breed. 2019; 55: 110-9.

[23] Rui S, Qi G, Shuangxi F et al. Analysis of genetic diversity in purple lettuce (Lactuca sativa L.) by SSR markers . Pak J Bot. 2020; 52: 181-96.

[24] Choi SP, Sim SC, Hong JH et al. Genetic characterisation of commercial Chinese cabbage varieties using SSR markers. Seed Sci Technol. 2016; 44: 595-608.

[25] Kong QS, Liu Y, Xie JJ et al. Development of simple sequence repeat markers from de novo assembled transcriptomes of pumpkins. Plant Mol Biol Rep. 2020; 38: 130-6.

[26] Hong JH, Kwon YS, Mishra RK et al. Construction of EST-SSR databases for effective cultivar identification and their applicability to complement for lettuce (Lactuca sativa L.) distinctness test . Am J Plant Sci. 2015; 06: 113-25.

[27] Raulier P, Maudoux O, Notté C et al. Exploration of genetic diversity within Cichorium endivia and Cichorium intybus with focus on the gene pool of industrial chicory . Genet Resour Crop Evol. 2016; 63: 243-59.

[28] Nadeem MA, Nawaz MA, Shahid MQ et al. DNA molecular markers in plant breeding: current status and recent advancements in genomic selection and genome editing. Biotechnol Biotechnol Equip. 2018; 32: 261-85.

[29] Mammadov J, Aggarwal R, Buyyarapu R et al. SNP markers and their impact on plant breeding. Int J Plant Genomics. 2012; 2012: 1-11.

[30] Plomion C, Bartholome J, Lesur I et al. High-density SNP assay development for genetic analysis in maritime pine (Pinus pinaster) . Mol Ecol Resour. 2016; 16: 574-87.

[31] Truong HT, Marcos Ramos A, Yalcin F et al. Sequence-based genotyping for marker discovery and codominant scoring in germplasm and populations. PLoS One. 2012; 7: 1-9.

[32] Kwon S, Simko I, Hellier B et al. Genome-wide association of 10 horticultural traits with expressed sequence tag-derived SNP markers in a collection of lettuce lines. Crop J. 2013; 1: 25-33.

[33] Park S, Kumar P, Shi A et al. Population genetics and genome-wide association studies provide insights into the influence of selective breeding on genetic variation in lettuce. Plant Genome. 2021; 14: 1-12.

[34] El-Esawi MA . Molecular genetic markers for assessing the genetic variation and relationships in Lactuca germplasm. Annu Res Rev Biol. 2015; 8: 1-13.

[35] Kandel JS, Peng H, Hayes RJ et al. Genome-wide association mapping reveals loci for shelf life and developmental rate of lettuce. Theor Appl Genet. 2020; 133: 1947-66.

[36] Egan AN, Schlueter J, Spooner DM . Applications of next-generation sequencing in plant biology. Am J Bot. 2012; 99: 175-85.

[37] Thomson MJ . High-throughput SNP genotyping to accelerate crop improvement. Plant Breed Biotechnol. 2014; 2: 195-212.

[38] Kumar S, Banks TW, Cloutier S . SNP discovery through next-generation sequencing and its applications. Int J Plant Genomics. 2012; 2012: 831460.

[39] D’Agostino N, Tripodi P . NGS-based genotyping, high-throughput phenotyping and genome-wide association studies laid the foundations for next-generation breeding in horticultural crops. Diversity. 2017; 9: 38.

[40] Le Nguyen K, Grondin A, Courtois B et al. Next-generation sequencing accelerates crop gene discovery. Trends Plant Sci. 2019; 24: 263-74.

[41] Deschamps S, Llaca V, May GD . Genotyping-by-sequencing in plants. Biology. 2012; 1: 460-83.

[42] Poland JA, Rife TW . Genotyping-by-sequencing for plant breeding and genetics. Plant Genome. 2012; 5: 92-102.

[43] Chung YS, Choi SC, Jun TH et al. Genotyping-by-sequencing: a promising tool for plant genetics research and breeding. Hortic Environ Biotechnol. 2017; 58: 425-31.

[44] Lee KJ, Lee JR, Sebastin R et al. Genetic diversity assessed by genotyping by sequencing (GBS) in watermelon germplasm. Genes. 2019; 10: 822.

[45] Phan NT, Trinh LT, Rho MY et al. Identification of loci associated with fruit traits using genomewide single nucleotide polymorphisms in a core collection of tomato (Solanum lycopersicum L.) . Sci Hortic. 2019; 243: 567-74.

[46] Palumbo F, Qi P, Pinto VB et al. Construction of the first SNP-based linkage map using genotyping-by-sequencing and mapping of the male-sterility gene in leaf chicory. Front Plant Sci. 2019; 10: 276.

[47] Pereira-Dias L, Vilanova SA, Prohens J et al. Genetic diversity, population structure, and relationships in a collection of pepper (Capsicum spp.) landraces from the Spanish Centre of diversity revealed by genotyping-by-sequencing (GBS) . Hortic Res. 2019; 6: 54.

[48] Reuter JA, Spacek DV, Snyder MP . High-throughput sequencing technologies. Mol Cell. 2015; 58: 586-97.

[49] Slatko BE, Gardner AF, Ausubel FM . Overview of next-generation sequencing technologies. Curr Protoc Mol Biol. 2018; 122: e59.

[50] Wang J, Lin M, Crenshaw A et al. High-throughput single nucleotide polymorphism genotyping using nanofluidic dynamic arrays. BMC Genomics. 2009; 10: 561.

[51] Kim H, Yoon JB, Lee J . Development of Fluidigm SNP type genotyping assays for marker-assisted breeding of chili pepper (Capsicum annuum L.) . Hortic Sci Technol. 2017; 35: 465-79.

[52] Choi SR, Heon OS, Dhandapani V et al. Development of SNP markers for marker-assisted breeding in Chinese cabbage using Fluidigm genotyping assays. Hortic Environ Biotechnol. 2020; 61: 327-38.

[53] Seo J, Lee G, Jin Z et al. Development and application of indica- japonica SNP assays using the Fluidigm platform for rice genetic analysis and molecular breeding . Mol Breeding. 2020; 40: 39.

[54] Lin Y, Yu W, Zhou L et al. Genetic diversity of oolong tea (Camellia sinensis) germplasm based on the nanofluidic array of single-nucleotide polymorphism (SNP) markers . Tree Genet Genomes. 2020; 16: 3.

[55] Yoo K, Jang S . Intraspecific relationships of Lactuca sativa var. capitata cultivars based on RAPD analysis . Korean J Hort Sci Tech. 2003; 21: 273-8.

[56] Sharma S, Kumar P, Gambhir G et al. Assessment of genetic diversity in lettuce (Lactuca sativa L.) germplasm using RAPD markers . 3 Biotech. 2018; 8: 9.

[57] Yang TJ, Jang SW, Kim WB . Genetic relationships of Lactuca spp. revealed by RAPD, inter SSR, AFLP and PCR-RFLP analyses . J Crop Sci Biotechnol. 2007; 10: 27-32.

[58] Al-Redhaimam K, Motawei MI, Shinawy M et al. Assessment of genetic variation and presence of nitrate reductase gene (NR) in different lettuce genotypes using PCR-based markers. J Food Agric Environ. 2005; 3: 134-6.

[59] Lebeda A, Kristkova E, Kitner M et al. Wild Lactuca species, their genetic diversity, resistance to diseases and pests, and exploitation in lettuce breeding. Eur J Plant Pathol. 2014; 138: 597-640.

[60] Zhang L, Su W, Tao R et al. RNA sequencing provides insights into the evolution of lettuce and the regulation of flavonoid biosynthesis. Nat Commun. 2017; 8: 2264.

[61] He J, Zhao X, Laroche A et al. Genotyping-by-sequencing (GBS), an ultimate marker-assisted selection (MAS) tool to accelerate plant breeding. Front Plant Sci. 2014; 5: 484.

[62] Kerbiriou PJ, Maliepaard CA, Stomph TJ et al. Genetic control of water and nitrate capture and their use efficiency in lettuce (Lactuca sativa L.) . Front Plant Sci. 2016; 7: 343.

[63] Pelgrom AJ, Eikelhof J, Elberse J et al. Recognition of lettuce downy mildew effector BLR38 in Lactuca serriola LS102 requires two unlinked loci. Mol Plant Pathol. 2019; 20: 240-53.

[64] Thiel T, Michalek W, Varshney RK et al. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.) . Theor Appl Genet. 2003; 106: 411-22.

[65] Meszaros K, Karsai I, Kuti C et al. Efficiency of different marker systems for genotype fingerprinting and for genetic diversity studies in barley (Hordeum vulgare L.) . S Afr J Bot. 2007; 73: 43-8.

[66] Singh RK, Sharma RK, Singh AK et al. Suitability of mapped sequence tagged microsatellite site markers for establishing distinctness, uniformity and stability in aromatic rice. Euphytica. 2004; 135: 135-43.

[67] Yang X, Ren R, Ray R et al. A genetic diversity and population structure of core watermelon (Citrullus lanatus) genotypes using DArTseq-based SNPs . Plant Genet Resour. 2016; 14: 226-33.

[68] Chen W, Hou L, Zhang Z et al. Genetic diversity, population structure, and linkage disequilibrium of a core collection of Ziziphus jujuba assessed with genome-wide SNPs developed by genotyping-by-sequencing and SSR markers. Front Plant Sci. 2017; 8: 575.

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

[70] Martin M . Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011; 17: 10-2.

[71] Cox MP, Peterson DA, Biggs PJ . SolexaQA: at-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinformatics. 2010; 11: 485.

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

[73] Kim JE, Oh SK, Lee JH et al. Genome-wide SNP calling using next generation sequencing data in tomato. Mol Cells. 2014; 37: 36-42.

[74] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[75] Tamura K, Stecher G, Peterson D et al. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013; 30: 2725-9.

[76] 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.

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

[78] 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.

[79] Rohlf FJ . NTSYSpc (Numerical Taxonomy and Multivariate Analysis System). Version 2.2, Exeter Software, Applied Biostatistics Inc., New York, USA.(2005).