RESEARCH ARTICLE https://doi.org/10.12972/kjhst.20190011 Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental ( melo L. var. makuwa)

Nam Ngoc Nguyen1†, Yong-Sham Kwon2†, Jae-Ryoung Park2, and Sung-Chur Sim1,3*

1Department of Bioresources Engineering, Sejong University, 05006, Korea 2Department of Molecular Genetics, Dong-A University, Busan 49315, Korea 3Plant Engineering Research Institute, Sejong University, Seoul 05006, Korea

*Corresponding author: [email protected]

†These authors contributed equally to this work.

Abstract

Oriental melon ( L. var. makuwa), also called ‘Korean melon’ or ‘Chamoe’, is widely grown in , including Korea, Japan, and China. We used 521 simple sequence repeats (SSRs) to develop a core set of SSR markers as an effective tool for cultivar identification and seed purity tests in oriental melon. The polymorphism test using a subset of seven cultivars identified 20 expressed sequence tag SSR (EST-SSR) and 23 genomic SSR markers that showed high levels of polymorphism and simple segregation patterns. These SSR markers were used for DNA fingerprinting Received: July 23, 2018 in additional 46 commercial cultivars. The average number of alleles per marker was 2.95 and 4.61 Revised: October 7, 2018 for the EST-SSR and genomic SSR markers in the 53 cultivars, respectively. For polymorphic Accepted:November 10, 2018 information content (PIC), we found 0.411 to 0.816 with an average of 0.536 for the EST-SSR markers and 0.453 to 0.862 with an average of 0.568 for the genomic SSR markers. The unweighted OPEN ACCESS pair group mean algorithm (UPGMA) dendrograms with the three sets of SSR markers (20 EST-SSR, 23 genomic SSR, and all 43 SSR markers) revealed that the 53 cultivars were separated into 4 to 5 HORTICULTURAL SCIENCE and TECHNOLOGY 37(1):119-129, 2019 major clusters. Furthermore, the set of all SSR markers was more effective at identifying the 53 cultivars URL: http://www.kjhst.org relative to the other two marker sets. Therefore, the 43 SSR markers as a core set will be a valuable pISSN : 1226-8763 tool for several applications, including cultivar identification and seed purity tests in oriental melon. eISSN : 2465-8588 Additional key words: breeding, DNA fingerprinting, genetic diversity, molecular markers, polymorphic This is an Open Access article distributed under the terms of the Creative Commons information content Attribution Non-Commercial License which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. Introduction Copyrightⓒ2019 Korean Society for Horticultural Science. Melon (Cucumis melo L.) is an important cultivated species (2n = 2x = 24) in the family that contains about 965 species representing 95 genera (Christenhusz and Byng, 2016). This This work was supported by Korea Institute of Planning and Evaluation for Technology in crop shows significant variations in fruit traits, such as fruit shape (round, ovoid, oblong, slender, and Food, Agriculture, Forestry and Fisheries (IPET) through Agri-Bioindustry Technology Development flat), fruit surface (smooth, warty, wrinkled, and ribbed), flesh color (orange, yellow, green, pink, and Program funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA) (No. 316087-04) white), and fruit weight (50 g to 15 kg) in commercial cultivars (Nunez-Palenius et al., 2008; to S. Sim. Monforte et al., 2014; Leida et al., 2015). With great nutrient values and unique flavors, melon fruit

Horticultural Science and Technology 119 Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

is rich in soluble sugars, vitamin C, -carotene, minerals, and antioxidant compounds (Lester, 1997; Wolbang et al., 2008; Leida et al., 2015; Zeb, 2016). Melon is divided into six subspecies, including inodorus, reticulatus, cantalupensis, acidulous, saccharinus, and makuwa (Liu et al., 2004). The subspecies C. melo L. var. makuwa (oriental melon) is also called ‘Korean melon’ or ‘Chamoe’. Generally, oriental melon fruit is ellipsoid or oval-oblong in shape, with a yellow to orange-yellow smooth surface and white longitudinal stripes (Lim, 2012). This subspecies originated from East India and has been cultivated widely in East Asia, particularly in the temperate regions of Korea, Japan, and China (Lim, 2012; Kim et al., 2016). Because of its economic value, many efforts have been made to develop elite cultivars in breeding programs. Thus, it has been important to protect the breeders’ rights for new cultivars. The distinctness, uniformity, and stability (DUS) test was established by the International Union for the Protection of New Varieties of Plants (UPOV, 2002). It is a standard method for evaluating new cultivars for variety protection. However, the DUS test requires labor-intensive and time-consuming evaluation based on phenotypes. Molecular markers have been used to assess genetic diversity in crops (Mondini et al., 2009). Of the several different types of molecular markers, simple sequence repeats (SSRs) have been used as a marker for genetic study due to their codominance, multiallelic nature, and high reproducibility (Vieira et al., 2016). The SSR motifs are found on both coding and noncoding sequences. Expressed sequence tag SSRs (EST-SSRs) from coding sequences show lower levels of polymorphism but higher cross-species transferability relative to genomic SSRs from noncoding sequences (Varshney et al., 2005). In the last two decades, a large number of SSR markers have been developed and applied to DNA fingerprinting, genetic diversity analysis, and linkage analysis in crop species (Danin-Poleg et al., 2001; Guo et al., 2007; Fernandez-Silva et al., 2008; Zhang et al., 2014., Phan et al., 2016). In melon, SSR markers have been developed and used to investigate genetic diversity in germplasm collections (Danin-Poleg et al., 2001; Ritschel et al., 2004; Kong et al., 2011; Aragão et al., 2013). The 61 Cucumis SSR markers were developed using genomic and cDNA libraries and database sequences (Danin-Poleg et al., 2001). Of these, the 31 SSR markers showed a total of 111 alleles in a collection of 13 melon accessions and revealed a distinction between the sweet cultivated melon and the ‘exotic’ melon. Ritschel et al. (2004) used an enriched genomic library to develop 144 SSR markers. A set of 25 SSR markers effectively identified 40 melon accessions representing the cultivated subspecies (inodorus, cantalupensis, cantalupensis/reticulatus) and noncultivated subspecies (conomon). Furthermore, the EST-SSR markers developed using the melon unigenes separated thin-skinned from thick-skinned melons in a collection of 27 Chinese melon accessions (Kong et al., 2011). Aragão et al. (2013) used 17 genomic SSR markers for genetic diversity analysis in the 38 melon accessions, which belong to different botanical groups and show characteristics of specific market classes. In this study, we curated and selected a core set of SSR markers to establish an accurate DNA database for a collection of 53 oriental melons. This core set of 23 genomic SSRs and 20 EST-SSRs was effective to distinguish all of these cultivars, suggesting the usefulness of these markers for cultivar identification and seed purity tests in oriental melon. This molecular tool will be valuable to protect the breeders’ intellectual property.

Materials and Methods

Plant Materials and DNA Isolation

A total of 53 oriental melon cultivars were collected from 16 seed companies in Korea (Table 1). Genomic DNA was isolated directly from the 5 to 7 seeds of each cultivar. Seeds were ground in liquid nitrogen, and genomic DNA was

120 Horticultural Science and Technology Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

extracted using the NucleoSpin® Plant II kit (Macherey-Nagel GmbH & Co. Germany). The quality and quantity of isolated DNA was measured by using electrophoresis on agarose gel with a quantitative DNA standard. The DNA solution was diluted to final concentration of 10 ng·µL-1 for SSR analysis.

Table 1. Fifty-three oriental melon cultivars for genetic characterization using SSR markers

No. Cultivar Seed Company No. Cultivar Seed Company 1 Geumssaragieuncheonz Dongbu Farm Hannong 28 Joeundaey Syngenta 2 Geumjatapy Dongbu Farm Hannong 29 Asakeekkulz Syngenta 3 Mallijangseongz Dongbu Farm Hannong 30 Manidaray Syngenta 4 Geumgwanz Dongbu Farm Hannong 31 Daechankkulz Myungsan 5 Jeongpumplusz Dongbu Farm Hannong 32 Hwanggeummulgyeoul Myungsan 6 Goldentopz Dongbu Farm Hannong 33 Palbogkkuly Hyundae 7 Jeongpumz Dongbu Farm Hannong 34 Hwangkkalkkuly Hyundae 8 Tamsrunggulz Dongbu Farm Hannong 35 Seongbakkkulz Sakata 9 Tamsrungoldz Dongbu Farm Hannong 36 Segyeilkkulz Sakata 10 Jewangkkuly Dongbu Farm Hannong 37 Nonghyeopkkulz Nonghyup 11 Bujakkuly Dongbu Farm Hannong 38 Kkulnaraz Nonghyup 12 Chilseongkkuly Dongbu Farm Hannong 39 Jangchoonmuskkkuly Jangchum 13 Superseongjukkulz Dongbu Farm Hannong 40 Myeongmunplusy Jangchum 14 Daewhangy Dongbu Farm Hannong 41 Crunchballz Samsung 15 Osungkkuly Dongbu Farm Hannong 42 Jinpeumkkurz Dongwon 16 Super007gguly Nongwoobio 43 Sungeumdeongeoriz Jeil 17 Arakkulz Nongwoobio 44 Sarangz Daeyeon 18 Geumnodajieuncheonz Nongwoobio 45 Chamhankkuly Hanil 19 Geumboseokkkulz Nongwoobio 46 Geumhongy Kimhoetae 20 Gayagguly Nongwoobio 47 Seongjaiwhangy GARESx 21 007gguly Nongwoobio 48 Seongjaiggany GARES 22 Gohyangz Nongwoobio 49 Seongjaiway GARES 23 Manittabulgguly Asia 50 Seongjaiheukseony GARES 24 Cheonhyangmilgoldy Asia 51 Yeolgoly GARES 25 Magicgoldy Asia 52 Euncheony GARES 26 Goldenbelly Asia 53 Gangseoy GARES 27 Gaeguriz Asia zCommercial variety for seed marketing. yCandidate or registration variety for plant variety protection. xGyeongsangbuk-do Agricultural Research & Extension Services.

SSR Genotyping

The 521 SSRs used in this study were derived from previous studies (Danin-Poleg et al., 2001; Fazio et al., 2002; Chiba et al., 2003; Ritschel et al., 2004; Gonzalo et al., 2005; Kong et al., 2007; Fernandez-Silva et al., 2008). To identify polymorphic SSRs, these SSRs were tested using a subset of the seven cultivars, including ‘Geumssaragieuncheon’, TM ‘Manidara’, ‘Chilseongkkul’, ‘Gohyang’, ‘Osungkkul’, ‘Superseongjukkul’, and ‘Hwangkkalkkul’ (Table 1). The HAD-GT12

Horticultural Science and Technology 121 Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

Genetic Analyzer (eGEnE, USA) was used for fragment analysis of SSRs. PCR reactions were prepared in a total volume

of 30 µL containing 50 mM KCl, 1.5 mM MgCl2, 10 mM Tris-HCl (pH 8.0), 0.2 µM of each SSR primer, 0.2 mM of each dNTP, 1 unit of Taq DNA polymerase, and 25 ng of DNA template. For amplification, we used a thermal program starting with a cycle at 94°C for 5 min for initial denaturation, followed by 40 cycles at 94°C for 30 s, 50-60°C for 30 s, 72°C for 45 s, and a cycle at 72°C for 10 min for final extension. A set of polymorphic SSR markers with simple segregation patterns was used to genotype additional 46 cultivars with the 3130XL Genetic Analyzer (Applied Biosystems, USA). The 5’ end of each forward primer for this analysis was labelled with FAM, VIC, NED, or PET fluorescent dye (Applied Biosystems, USA), and SSR alleles were determined using GeneMapper software version 3.7 (Applied Biosystems, USA).

Data Analysis

The polymorphism information content (PIC) value for each marker was estimated using the following equation:

         

th in which n and pi are defined as the total number of alleles and the frequency of the i allele, respectively (Anderson et al., 1993). Furthermore, genetic relationships in our collection of oriental melon cultivars were investigated using three sets of SSR markers, including EST-SSRs, genomic SSRs, and all combined SSRs. Genetic distances between cultivars were estimated using the Jaccard’s method (Sneath and Sokal, 1973) as implemented in NTSYSpc software version 2.2 (Rohlf, 2008). The unweighted pair group mean algorithm (UPGMA) was used to construct a dendrogram based on the resulting genetic distances. In addition, the Mantel test (Mantel, 1967) was performed to estimate correlation between the genetic distances from EST-SSRs and genomic SSRs using the MXCOMP module in the NTSYSpc software.

Results and Discussion

A Core Set of SSR Markers

The polymorphism of the 521 SSRs was tested using a subset of seven commercial cultivars. Of these, 185 SSRs were polymorphic in these cultivars and 43 SSRs (20 EST-SSRs and 23 genomic SSRs) showed clear and simple segregation patterns for DNA fingerprinting (Fig. 1). The 20 EST-SSR markers detected 59 alleles with an average of 2.95 alleles per marker in the 53 oriental melon cultivars (Table 2). The seven SSR markers (ECM91, ECM125, ECM197, ECM204, ECM205, ECM217, and TJ10) detected only two alleles, while the 13 SSR markers (CM16, CM38, ECM54, ECM81, ECM123, ECM124, ECM157, ECM172, ECM181, ECM228, ECM130, ECM133, and ECM60) detected 3 to 7 alleles. The PIC values of the EST-SSR markers ranged from 0.411 (ECM157) to 0.816 (ECM60) with an average of 0.536. In a previous study, 19 EST-SSR markers showed an average of 3.0 alleles per marker in 27 melon cultivars (Kong et al., 2011), which is comparable with that of the EST-SSR markers in this study. Interestingly, our EST-SSR markers showed higher average PIC values relative to 0.455 for the 19 EST-SSR markers.

122 Horticultural Science and Technology Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

AB

CD

Fig. 1. Polymorphism of four SSR markers: ECM60 (A), ECM130 (B), ECM181 (C), and CMGAN12 (D). The PCR amplicons were analyzed using the HAD-GT12TM Genetic Analyzer. M, QX DNA size marker (25-450bp); lane 1, ‘Geumssaragieuncheon’; 2, ‘Manidara’; 3, ‘Chilseongkkul’; 4, ‘Gohyang’; 5, ‘Osungkkul’; 6, ‘Superseongjukkul’; and 7, ‘Hwangkkalkkul’.

Table 2. Repeat motif, number of alleles, and PIC value of 20 EST-SSR markers selected for genetic characterization of oriental melon cultivars

y z Tm Product No. of PIC SSR name Primer sequence Repeat motif (°C) size (bp) alleles value F:PET-TGCCTGTTGTGATTGAGGAG CM16 (TTC)8 55 197-213 3 0.529 R:TTCTTCTTACCTCCGCCAAA F:PET-TAGCATCTGATCGGAAAACC CM38 (TCA)15 55 118-129 3 0.564 R:CAACTTCATCCGCCAAGAAT F:VIC-TCGTGATTTTGCGTCCAGTA ECM54 (AG)17 55 78-107 3 0.550 R:TGCAATCAGGAATGAACAACTC F:NED-GGCCCTCAGGCAGTTTTT ECM60 (AC)11 55 139-174 7 0.816 R:TGCTCCTATGGTTGACCTTTTT F:FAM-CAACCATTCCTCCCATTCAT ECM81 (AG)11 55 123-140 3 0.511 R:CACCACCTGTGACATTGTACG F:VIC-GTTGAACATGGAAGAGTCTGCT ECM91 (ATC)16 55 269-291 2 0.482 R:AAAGAACCAGCCCTATCCAAA F:PET-TCTAATGGCGGCTTCAACTTA ECM123 (AAG)13 55 224-242 3 0.501 R:CTCCTTAGTGCATGGCTTCAC F:NED-GCGTCCTAAAAAGGGATAAGG ECM124 (AAG)13 55 268-277 3 0.444 R:ATTTTCACAAAAGGGGGAGAG F:FAM-GGAAAACGCAAAATCAGTGAG ECM125 (AAG)13 55 112-133 2 0.471 R:CTGAACGTGGACGACATTTTT F:PET-GATTGGGAAAAAGGGTATGGA ECM130 (AAAG)9 55 216-243 4 0.689 R:CTGGCTCCTTCACATTGTTGT F:VIC-AAACATCAACACACACCCACA ECM133 (ATC)13 55 94-121 4 0.672 R:TCAGCGACGGTCATCTATTTT F:FAM-TGTCCGTTACACTCAATCTCTC ECM157 (AG)10 55 154-169 3 0.411 R:GGATGACGGTCACAGGAAGT F:PET-CCCCTTCTCCCATTTCTACC ECM172 (AG)9 55 107-115 3 0.521 R:TGAGACGTGGCAGAGAAGAA

Horticultural Science and Technology 123 Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

Table 2. Repeat motif, number of alleles, and PIC value of 20 EST-SSR markers selected for genetic characterization of oriental melon cultivars (Continued)

Tmy Product No. of PIC SSR name Primer sequencez Repeat motif (°C) size (bp) alleles value F:VIC-CCCCATCTCAATTTGTCACC ECM181 (AC)10 55 116-135 3 0.668 R:AGGTGGAGGAAGGGGTTAAA F:VIC-TGTCTTCGCATCAATCCTACC ECM197 (AAG)11 55 80-92 2 0.417 R:TTTTGAATCACACCCTTCTGC F:NED-CTCTCTTCATTTCCCCTCGTT ECM204 (AAG)15 55 308-324 2 0.493 R:TGGCCTGGAAAGTAAGGGTAT F:NED-GATTCCTTGGGCTTGTACCTC ECM205 (AAG)10 55 96-115 2 0.474 R:GCTAAGGAAAGGGTTTGTTCG F:FAM-AGGAGAGAGAGCTCCGAAAAA ECM217 (AAAG)7 55 257-261 2 0.485 R:CAATTCCACCCTTCCAATTCT F:PET-TCTGATCGGAAAACCCACTT ECM228 (CAT)15 55 103-114 3 0.551 R:GCCAAGAATTTTCCCAACAT F:VIC-ACGAGGAAAACGCAAAATCA TJ10 (CTT)5(CT)3 55 115-137 2 0.478 R:TGAACGTGGACGACATTTTT Total 59 10.727 Average 2.950 0.536 zForward primers were labeled with fluorescent dyes. yMelting temperature. SSR marker sources: CM markers from Kong et al. (2007), ECM markers from Fernandez-Silva et al. (2008), and TJ marker from Gonzalo et al. (2005).

For the 23 genomic SSR markers, we detected a total of 106 alleles with an average of 4.61 alleles per marker in the cultivar collection (Table 3). The 11 genomic SSR markers detected over five alleles: CMBR83 and CMBR95 (five alleles); CMTCN6 and CMGAN12 (six alleles); CMBR7, CMBR25, CMBR64, CMBR78, CMBR98, and CMBR109 (seven alleles); and CMTCN18 (eight alleles). These genomic SSR markers showed the PIC values ranging from 0.453 (CMBR83 and CMBR95) to 0.862 (CMGAN12) (Table 3). In the study of Ritschel et al. (2004), 25 genomic SSR markers detected an average of 3.72 alleles per marker in a collection of 40 commercial melon cultivars. Aragão et al. (2013) found an average of 2.41 alleles per marker for 17 genomic SSR markers in 38 melon accessions representing different botanical groups. Park et al. (2013) reported seven genomic SSR markers to identify six oriental melon accessions. These markers only produced a total of 17 alleles with 2 to 3 alleles per marker. In contrast, we detected a considerably higher number of alleles per marker relative to the previous studies (an average of 4.61). The EST-SSR markers are developed using coding sequences that are more conserved than noncoding sequences. Since the genomic SSR markers are often derived from noncoding genomic sequences, higher levels of polymorphism are expected in genomic SSR markers than EST-SSR markers (Guo et al., 2007; Hu et al., 2011; Mujaju et al., 2013; Tabbasam et al., 2013). However, both EST-SSR and genomic SSR markers used in this study showed comparable ranges for the PIC values. This result suggests that the resulting EST-SSR markers in this study are effective as genomic SSR markers to identify genetic variations in the oriental melon cultivars. The core set of 43 SSR markers (20 EST-SSR and 23 genomic SSR markers) in this study detected a total of 165 alleles with an average of 3.84 alleles per marker in our collection of 53 oriental melon cultivars. The average PIC value of these SSR markers was 0.553.

124 Horticultural Science and Technology Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

Table 3. Repeat motif, number of alleles, and PIC value of 23 genomic SSR markers selected for genetic characterization of oriental melon cultivars

y SSR z Tm Product No. of PIC Primer sequence Repeat motif name (°C) size (bp) alleles value F:NED-AAAGAAAGAAACGCAACTTCG CMBR5 (TC)18 55 99-109 3 0.501 R:TGGGAAAGAAAGAGAATCAAAA F:FAM-AAAATGAATGGGAGTGCGTG CMBR7 (AG)30 55 87-113 7 0.601 R:GCCTTCCTTTTCACCATCAA F:PET-TGGGGTTGTCAATACAGCAA CMBR25 (TC)24 55 146-173 7 0.615 R:GGAGTGCGTGGAATGTACG F:VIC-CGACAATCACGGGAGAGTTT CMBR40 (CT)15T3CT2(CT)2 55 146-161 3 0.525 R:TTGTTGCATCAAACTAACACAATC F:NED-ATACAGCAGATCCACAGGGG CMBR64 (CT)24 54 134-163 7 0.616 R:ATGGGAGTGTGTGGGATGTA F:NED-AAAATGAATGGGAGTGCGTG CMBR78 (GA)(AG)(GA)21 55 90-115 7 0.616 R:TTGCCTTCCTGTTCACCATC F:VIC-CGGACAAATCCCTCTCTGAA CMBR83 (GA)21(CA)(GA)2 55 127-166 5 0.453 R:GAACAAGCAGCCAAAGACG F:NED-TTGACCTTTTACGGTGGTCC CMBR95 (CT)N(CT)2N2(CT)21 55 104-143 5 0.453 R:CGGACAAATCCCTCTCTGAA F:FAM-CGACAATCACGGGAGAGTTT (CA)2N2(CT)N(GA)N3(CT)N3 CMBR97 55 164-178 3 0.525 R:CATATTAGACCCATATTTGTTGCAT (GA)N5(CAA)2N4(CT) F:PET-ATACAGCAGATCCACAGGGG CMBR98 (CCTT)2T2(CAT)2CA2(CT)21 55 140-168 7 0.604 R:TGAATGGGAGTGTGTGGAAC F:VIC-TGGGTGAGACATCAACATCC (GA)N3(GAAA)2GA(GAAA)N5(GA)7A CMBR99 55 172-181 3 0.501 R:AGTTTATTTGGCCTCTCCCC 2(GA)4G5(GA)2 F:VIC-TGGAATGTACCGTGATGGGT CMBR109 (GA)23 55 125-152 7 0.616 R:ATACAGCAGATCCACAGGGG F:NED-CTGGCCCCCTCCTAAACTAA CMBR120 (CT)2N11(CT)2G2(CT)T(CT)N7(TC)15 55 161-173 3 0.648 R:CAAAAAGCATCAAAATGGTTG F:FAM-CAACCACCGACTACTAAGTC CMACC146 (ACC)9 55 141-153 4 0.631 R:CGACCAAACCCATCCGATAA F:VIC-GCATTGCTCGATCAGTTTTAC CMTCN6 (TC)2TT(TC)16 55 140-162 6 0.619 R:ACTCCGTCAAGATCCCAAAA F:NED-TTTTTGTCGTTATATGAGGG CMGAN12 (GA)21 55 178-206 6 0.862 R:GTTGCATAATGCTAATTTGG F:FAM-ACCAATCCATCACTCTCACT CMTCN18 (TC)11N24(TC)17(TA)12 55 147-215 8 0.702 R:GAAAGAATGGGGGAAAAGAG F:PET-TCTACTTCCATGAATCCATC CMTCN50 (TC)16(TA)12 55 120-134 4 0.526 R:TAGAATGGTTAGGAAACCCT F:VIC-AAAACGGAGATGTGGAGGAC GCM246 (AG)16 55 164-179 3 0.478 R:TTAAGCAAGCAGCCAAAATG F:NED-GCCGGAAAATTCTCTTACCA GCM303 (AG)8 55 159-165 2 0.500 R:TTGACGTTGATATGAGCAAGC F:FAM-TGGTTCAAATCCTCAGTGGT GCM521 (GA)5GG(GA)7 55 146-161 2 0.497 R:CAAGGGATTCTTCCATTTCG F:NED-GATGCGGTGAGAAAGAGTTGAGAGA CMMS12-4 (GA)15(GA)9 55 451-473 2 0.472 R:AGAGGGAGAGAGTTTGTAAAAAAAT F:PET-TTCTGCATACCCTCTCCT CSWCT02B (CT)21(TG)8 55 169-183 2 0.500 R:CACACTTCCAGATGGTTG Total 106 13.061 Average 4.609 0.568 zForward primers were labeled with fluorescent dyes. yMelting temperature. SSR marker sources: CMBR markers from Ritschel et al. (2004), CMACC marker from Danin-Poleg et al. (2001), CMTCN and CMGAN markers from Gonzalo et al. (2005), GCM markers from Fernandez-Silva et al. (2008), CMMS marker from Chiba et al. (2003), and CSWCT marker from Fazio et al. (2002).

Horticultural Science and Technology 125 Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

Genetic Diversity Analysis and Cultivar Identification

Three sets of SSR markers (20 EST-SSR, 23 genomic SSR, and all 43 SSR markers) were used to estimate genetic distances in the 53 oriental melon cultivars. The Mantel test showed a significant correlation (r = 0.767) between pairwise genetic distances from the EST-SSR and genomic SSR markers. However, different clustering patterns were observed with these two marker sets (Figs. 2 and 3). In the UPGMA dendrogram, the EST-SSR marker set divided the 53 cultivars into five clusters at Jaccard’s coefficients of 0.42 to 0.54 (Fig. 2). The largest cluster (cluster I) was further separated into three subclusters: I-1 (24 cultivars), I-2 (nine cultivars), and I-3 (two cultivars). The two cultivars ‘Geumnodajieuncheon’ and ‘Nonghyeopkkul’ were not separated in subcluster I-1. Clusters II and III consisted of seven cultivars and one cultivar, respectively. The cultivar in cluster III was ‘Jangchoonmuskkkul’ from the Jangchum seed company. The other cultivar, ‘Myeongmunplus’, developed by the same company, was grouped in cluster II. All seven cultivars collected from Gyeongsangbuk-do Agricultural Research & Extension Services (GARES) were found in a single cluster (cluster V). The 23 genomic SSR markers showed five major clusters at Jaccard’s coefficients of 0.19 to 0.37 (Fig. 3). The 43 cultivars in the cluster I were further divided into three subclusters: I-1 (22 cultivars), I-2 (18 cultivars), and I-3 (three cultivars). In subcluster I, the three cultivars ‘Geumssaragieuncheon’, ‘Jeongpum’, and ‘Geumnodajieuncheon’ were not separated. Additionally, ‘Bujakkul’ was not distinguished from ‘Chilseongkkul’ in subcluster I-2. Cluster II included only two cultivars, ‘Gayaggul’ and ‘Sarang’. These cultivars were clustered with the cultivar ‘Geumhong’ in the dendrogram based on the EST-SSR markers. The genomic SSR markers showed that cluster IV only had cultivar ‘Geumhong’. The seven cultivars from GARES were separated into cluster III (five cultivars) and cluster V (two cultivars).

Fig. 2. Dendrogram depicting the classification of 53 oriental melon cultivars constructed using UPGMA and based on the 20 EST-SSR markers. The numbers (1 to 53) at the right refer to the list of cultivars in Table 1. The scale at the bottom is Jaccard’s coefficient of similarity.

126 Horticultural Science and Technology Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

Fig. 3. Dendrogram depicting the classification of 53 oriental melon cultivars constructed using UPGMA and based on the 23 genomic SSR markers. The numbers (1 to 53) at the right refer to the list of cultivars in Table 1. The scale at the bottom is Jaccard’s coefficient of similarity.

Fig. 4. Dendrogram depicting the classification of 53 oriental melon cultivars constructed using UPGMA and based on the 43 EST-SSR and genomic SSR markers. The numbers (1 to 53) at the right refer to the list of cultivars in Table 1. The scale at the bottom is Jaccard’s coefficient of similarity.

Horticultural Science and Technology 127 Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

In contrast to the EST-SSR and genomic SSR marker sets, clustering analysis using all marker sets revealed four distinct clusters (Fig. 4). Cluster I consisted of 37 cultivars and the other clusters consisted of three to seven accessions. The three sets of SSR markers generated different dendrograms. This discrepancy was previously reported in other crop species and could be due to different genomic regions covered by both the markers and the variation in number of markers (Tabbasam et al., 2013; Parthiban et al., 2018). In addition, the EST-SSR markers are derived from expressed sequences compared to genomic sequences that are spread throughout the genome. The genomic SSR markers are often developed using the noncoding sequences, such as intron and intergenic regions. Using the combined set of 20 EST-SSR and 23 genomic SSR markers, we found that the 53 oriental melon cultivars were clearly differentiated, suggesting effectiveness of these SSR markers for cultivar identification and other applications in oriental melon breeding programs. The use of molecular markers is an efficient strategy for detecting genetic variations in crop species. SSR markers are often used as markers for DNA profiling and cultivar identification in several vegetables, including pumpkin (Sim et al., 2016), tomato (Phan et al., 2016), and (Kwon and Choi, 2013). In this study, we generated the DNA profiles of 53 oriental melon cultivars using a core set of 43 SSR markers. These SSR markers showed fairly high levels of PIC values relative to the previous SSR markers reported in melon. Furthermore, our SSR markers were effective in separating all 53 oriental melon cultivars with narrow genetic backgrounds. Therefore, the core set of 43 SSR markers will be a valuable tool for protecting breeders’ intellectual property. Rapid progress in next-generation sequencing technologies has led to the discover of genome-wide single nucleotide polymorphisms (SNPs) in melon (Natarajan et al., 2016; Sim et al., 2018). Since the SNP markers provide a cost-effective method for high-throughput genotyping, use of SNP markers along with our SSR markers can leverage development of an advanced DNA profiling pipeline for cultivar identification and seed purity tests in oriental melon.

Literature Cited

Anderson JA, Churchill GA, Autrique JE, Tanksley SD, Sorrells ME (1993) Optimizing parental selection for genetic linkage maps. Genome 36:181-186. doi:10.1139/g93-024 Aragão FA, Torres Filho J, Nunes GH, Queiroz MA, Bordallo PN, Buso GS, Ferreira MA, Costa ZP, Bezerra Neto F (2013) Genetic divergence among accessions of melon from traditional agriculture of the Brazilian Northeast. Genet Mol Res 12:6356-6371. doi:10.4238/2013.December.6.3 Chiba N, Suwabe K, Nunome T, Hirai M (2003) Development of microsatellite markers in melon (Cucumis melo L.) and their application to major cucurbit crops. Breeding Sci 53:21-27. doi:10.1270/jsbbs.53.21 Christenhusz MJM, Byng JW (2016) The number of known plants species in the world and its annual increase. Phytotaxa 261:201. doi:10.11646/phytotaxa.261.3.1 Danin-Poleg Y, Reis N, Tzuri G, Katzir N (2001) Development and characterization of microsatellite markers in Cucumis. Theor Appl Genet 102:61-72. doi:10.1007/s001220051618 Fazio G, Staub JE, Chung SM (2002) Development and characterization of PCR markers in cucumber. J Amr Soc Hort Sci 127:545-557. Fernandez-Silva I, Eduardo I, Blanca J, Esteras C, Picó B, Nuez F, Arús P, Garcia-Mas J, Monforte AJ (2008) Bin mapping of genomic and EST-derived SSRs in melon (Cucumis melo L.). Theor Appl Genet 118:139-150. doi:10.1007/s00122-008-0883-3 Gonzalo MJ, Oliver M, Garcia-Mas J, Monfort A, Dolcet-Sanjuan R, Katzir N, Arus P, Monforte AJ (2005) Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.). Theor Appl Genet 110:802-811. doi:10.1007/s00122-004-1814-6 Guo W, Cai C, Wang C, Han Z, Song X, Wang K, Niu X, Wang C, Lu K, et al (2007) A microsatellite-based, gene-rich linkage map reveals genome structure, function and evolution in Gossypium. Genetics 176:527-541. doi:10.1534/genetics.107.070375 Hu J, Wang L, Li J (2011) Comparison of genomic SSR and EST-SSR markers for estimating genetic diversity in cucumber. Biol Plant 55:577-580. doi:10.1007/s10535-011-0129-0 Kim H-A, Shin A-Y, Lee M-S, Lee H-J, Lee H-R, Ahn J, Nahm S, Jo S-H, Park J-M, et al (2016) De novo transcriptome analysis of Cucumis melo L. var. makuwa. Mol Cells 39:141-148. doi:10.14348/molcells.2016.2264

128 Horticultural Science and Technology Development of a Core Set of SSR Markers for Cultivar Identification and Seed Purity Tests in Oriental Melon (Cucumis melo L. var. makuwa)

Kong Q, Xiang C, Yang J, Yu Z (2011) Genetic variations of Chinese melon landraces investigated with EST-SSR markers. Hortic Environ Biotechnol 52:163-169. doi:10.1007/s13580-011-0087-7 Kong Q, Xiang C, Yu Z, Zhang C, Liu F, Peng C, Peng X (2007) Mining and charactering microsatellites in Cucumis melo expressed sequence tags from sequence database. Mol Ecol Notes 7:281-283. doi:10.1111/j.1471-8286.2006.01580.x Kwon Y-S, Choi K-J (2013) Construction of a DNA profile database for commercial cucumber (Cucumis sativus L.) cultivars using microsatellite marker. J Korean Soc Hortic Sci 31:344-351. doi:10.7235/hort.2013.13008 Leida C, Moser C, Esteras C, Sulpice R, Lunn JE, de Langen F, Monforte AJ, Pico B (2015) Variability of candidate genes, genetic structure and association with sugar accumulation and climacteric behavior in a broad germplasm collection of melon (Cucumis melo L.). BMC Genet 16:28. doi:10.1186/s12863-015-0183-2 Lester G (1997) Melon (Cucumis melo L.) fruit nutritional quality and health functionality. HortTechnology 7:222-227. doi:10.21273/ HORTTECH.7.3.222 Lim TK (2012) Edible medicinal and non-medicinal plants, Vol 2. Fruits. Springer, Dordrecht, Netherlands, pp 219-221 Liu L, Kakihara F, Kato M (2004) Characterization of six varieties of Cucumis melo L. based on morphological and physiological characters, including shelf-life of fruit. Euphytica 135:305. doi:10.1023/B:EUPH.0000013330.66819.6f Mantel N (1967) The detection of disease clustering and a generalized regression approach. Cancer Res 27:209-220 Mondini L, Noorani A, Pagnotta MA (2009) Assessing plant genetic diversity by molecular tools. Diversity 1:19–35. doi:10.3390/d1010019 Monforte AJ, Diaz A, Cano-Delgado A, van der Knaap E (2014) The genetic basis of fruit morphology in horticultural crops: lessons from tomato and melon. J Exp Bot 65:4625-4637. doi:10.1093/jxb/eru017 Mujaju C, Sehic J, Nybom H (2013) Assessment of EST-SSR markers for evaluating genetic diversity in accessions from Zimbabwe. Am J Plant Sci 4:1448-1456. doi:10.4236/ajps.2013.47177 Natarajan S, Kim HT, Thamilarasan SK, Veerappan K, Park JI, Nou IS (2016) Whole genome resequencing and characterization of powdery mildew disease-associated allelic variation in melon. PLoS ONE 11:e0157524. doi:10.1371/journal.pone.0157524 Nunez-Palenius HG, Gomez-Lim M, Ochoa-Alejo N, Grumet R, Lester G, Cantliffe DJ (2008) Melon fruits: genetic diversity, physiology, and biotechnology features. Crit Rev Biotechnol 28:13-55. doi:10.1080/07388550801891111 Park I, Kim J, Lee J, Kim S, Cho O, Yang K, Ahn J, Nahm S, Kim H (2013) Development of SSR markers by next-generation sequencing of Korean landraces of chamoe (Cucumis melo var. makuwa). Mol Biol Rep 40:6855-6862. doi:10.1007/s11033-013-2803-0 Parthiban S, Govindaraj P, Senthilkumar S (2018) Comparision of relative efficiency of genomic SSR and EST-SSR markers in estimating genetic diversity in sugarcane. 3 Biotech 8:144. doi:10.1007/s13205-018-1172-8 Phan NT, Kim M-K, Sim S-C (2016) Genetic variations of F1 tomato cultivars revealed by a core set of SSR and InDel markers. Sci Hortic 212:155-161. doi:10.1016/j.scienta.2016.09.043 Ritschel PS, Lins TCdL, Tristan RL, Buso GSC, Buso JA, Ferreira ME (2004) Development of microsatellite markers from an enriched genomic library for genetic analysis of melon (Cucumis melo L.). BMC Plant Biol 4:9. doi:10.1186/1471-2229-4-9 Rohlf FJ (2008) NTSYSpc: Numerical taxonomy system, version 2.20. Exeter Publishing, Setauket, NY, USA Sim S-C, Hong J-H, Kwon Y-S (2016) DNA profiling of commercial pumpkin cultivars using simple sequence repeat polymorphisms. Hortic Environ Biotechnol 56:811-820. doi:10.1007/s13580-015-0123-0 Sim S-C., Nguyen NN, Kim N, Kim J, Park Y (2018) Whole-genome resequencing reveals genome-wide single nucleotide polymorphisms between orange-fleshed and green-fleshed melons. Hortic Environ Biotechnol 59:275. doi:10.1007/s13580-018-0030-2 Sneath PHA, Sokal RR (1973) Numerical taxonomy: The principles and practice of numerical classification. W.H. Freeman and Company, San Francisco, CA, USA Tabbasam N, Zafar Y, Mehboob-ur-Rahman (2013) Pros and cons of using genomic SSRs and EST-SSRs for resolving phylogeny of the genus Gossypium. Plant Syst Evol 300:559–75. doi:10.1007/s00606-013-0891-x UPOV (2002) General introduction to the examination of distinctness, uniformity and stability and the development of harmonized descriptions of new varieties of plants. UPOV, Geneva, Switzerland Varshney RK, Graner A, Sorrells ME (2005) Genic microsatellite markers in plants: features and applications. Trends Biote chnol 23:48-55. doi:10.1016/j.tibtech.2004.11.005 Vieira ML, Santini L, Diniz AL, Munhoz Cde F (2016) Microsatellite markers: what they mean and why they are so useful. Genet Mol Biol 39:312-328. doi:10.1590/1678-4685-GMB-2016-0027 Wolbang CM, Fitos JL, Treeby MT (2008) The effect of high pressure processing on nutritional value and quality attributes of Cucumis melo L. Innov Food Sci Emerg Technol 9:196-200. doi:10.1016/j.ifset.2007.08.001 Zeb A (2016) Phenolic profile and antioxidant activity of melon (Cucumis melo L.) seeds from Pakistan. 5:67. doi:10.3390/foods5040067 Zhang M, Mao W, Zhang G, Wu F (2014) Development and characterization of polymorphic EST-SSR and genomic SSR markers for tibetan annual wild barley. PLoS ONE 9:e94881. doi:10.1371/journal.pone.0094881

Horticultural Science and Technology 129