Mining and validation of novel genotyping-by-sequencing (GBS)-based simple sequence repeats (SSRs) and their application for the estimation of the genetic diversity and population structure of coconuts (Cocos nucifera L.) in Thailand

Abstract

Coconut (Cocos nucifera L.) is an important economic crop in tropical countries. However, the lack of a complete reference genome and the limitations of usable DNA markers hinder genomic studies and the molecular breeding of coconut. Here, we present the results of simple sequence repeat (SSR) mining from a high-throughput genotyping-by-sequencing (GBS) study of a collection of 38 coconut accessions. A total of 22,748 SSRs with di-, tri-, tetra-, penta- and hexanucleotide repeats of five or more were identified, 2451 of which were defined as polymorphic loci based on locus clustering in 38 coconut accessions, and 315 loci were suitable for the development of SSR markers. One hundred loci were selected, and primer pairs for each SSR locus were designed and validated in 40 coconut accessions. The analysis of 74 polymorphic markers identified between 2 and 9 alleles per locus, with an average of 3.01 alleles. The assessment of the genetic diversity and genetic relationships among the 40 coconut varieties based on the analysis of population structure, principal coordinate analysis (PCoA), and phylogenetic tree analysis using the 74 polymorphic SSR markers revealed three main groups of coconuts in Thailand. The identified SSR loci and SSR markers developed in this study will be useful for the study of coconut diversity and molecular breeding. The SSR mining approach used in this study could be applied to other plant species with a complex genome regardless of the availability of reference genome.

References

[1] Gunn, B. F., Baudouin, L. & Olsen, K. M. Independent origins of cultivated coconut (Cocos nucifera L.) in the old world tropics . PLoS ONE 6, e21143 (2011).

[2] Xiao, Y. et al. The genome draft of coconut (Cocos nucifera) . Gigascience 6, 1-11 (2017).

[3] Chan, E. & Elevitch, C. R. Cocos nucifera (coconut) . Species profiles Pac. Isl. Agrofor. 2, 1-27 (2006).

[4] Menon K. P. V., Pandalai K. M. The Coconut. A Monograph. (Indian Central Coconut Committee, 1958).

[5] Lebrun, P. Berger, A. Hodgkin, T. Baudouin, L. Biochemical and molecular methods for characterizing coconut diversity. In: Batugal, P. Ramanatha Rao, V. Oliver, J. (eds) Coconut genetic resources. International Plant Genetic Resources Institute, Regional Office for Asia, the Pacific and Oceania (IPGRI-APO), Serdang, Selangor DE, Malaysia, pp 225-245 (2005).

[6] Kamaral, L. C. J., Perera, S. A. C. N., Perera, K. L. N. S. & Dassanayaka, P. N. Genetic diversity of the Sri Lanka yellow dwarf coconut form as revealed by microsatellite markers. Trop. Agric. Res. 26, 131 (2015).

[7] Kamaral, L. C. J., Perera, S., Perera, K. & Dassanayaka, P. N. Characterisation of Sri Lanka Yellow Dwarf Coconut (Cocos nucifera L.) by DNA fingerprinting with SSR markers . J. Natl Sci. Found. 45, 405 (2017).

[8] Xia, W. et al. Development of gene-based simple sequence repeat markers for association analysis in Cocos nucifera . Mol. Breed. 34, 525-535 (2014).

[9] Lebrun, P., N’cho, Y. P., Seguin, M., Grivet, L. & Baudouin, L. Genetic diversity in coconut (Cocos nucifera L.) revealed by restriction fragment length polymorphism (RFLP) markers . Euphytica 101, 103-108 (1998).

[10] Ashburner, G. R., Thompson, W. K. & Halloran, G. M. RAPD analysis of south pacific coconut palm populations. Crop Sci. 37, 992-997 (1997).

[11] Perera, L., Russell, J. R., Provan, J., McNicol, J. W. & Powell, W. Evaluating genetic relationships between indigenous coconut (Cocos nucifera L.) accessions from Sri Lanka by means of AFLP profiling . Theor. Appl. Genet. 96, 545-550 (1998).

[12] Perera, L., Russell, J. R., Provan, J. & Powell, W. Studying genetic relationships among coconut varieties/populations using microsatellite markers. Euphytica 132, 121-128 (2003).

[13] Geethanjali, S., Anitha Rukmani, J., Rajakumar, D., Kadirvel, P. & Viswanathan, P. L. Genetic diversity, population structure and association analysis in coconut (Cocos nucifera L.) germplasm using SSR markers . Plant Genet. Resour.: Charact. Utilization 16, 156-168 (2018).

[14] Grover, A., Aishwarya, V. & Sharma, P. C. Searching microsatellites in DNA sequences: approaches used and tools developed. Physiol. Mol. Biol. Plants 18, 11-19 (2012).

[15] Powell, W., Machray, G. C. & Provan, J. Polymorphism revealed by simple sequence repeats. Trends plant Sci. 1, 215-222 (1996).

[16] Guan, R. et al. Genetic diversity comparison between Chinese and Japanese soybeans (Glycine max (L.) Merr.) revealed by nuclear SSRs . Genet. Resour. Crop Evol. 57, 229-242 (2010).

[17] Ali, M. L. et al. Assessment of genetic diversity and relationship among a collection of US sweet sorghum germplasm by SSR markers. Mol. Breed. 21, 497-509 (2008).

[18] Shoemaker, R. C. et al. Microsatellite discovery from BAC end sequences and genetic mapping to anchor the soybean physical and genetic maps. Genome 51, 294-302 (2008).

[19] Stàgel, A., Portis, E., Toppino, L., Rotino, G. L. & Lanteri, S. Gene-based microsatellite development for mapping and phylogeny studies in eggplant. BMC Genomics 9, 357 (2008).

[20] Lebrun, P. et al. Construction of a linkage map of the Rennell Island Tall coconut type (Cocos nucifera L.) and QTL analysis for yield characters . Genome 44, 962-970 (2001).

[21] Baudouin, L. et al. QTL analysis of fruit components in the progeny of a Rennell Island Tall coconut (Cocos nucifera L.) individual . Theor. Appl. Genet. 112, 258-268 (2006).

[22] Neeraja, C. N. et al. A marker-assisted backcross approach for developing submergence-tolerant rice cultivars. Theor. Appl. Genet. 115, 767-776 (2007).

[23] Danson, J., Lagat, M., Kimani, M. & Kuria, A. Quantitative trait loci (QTLs) for resistance to gray leaf spot and common rust diseases of maize. Afr J Biotechnol. 7, 3247-3254 (2008).

[24] Vieira, M. L. C., Santini, L., Diniz, A. L., Munhoz, C. & de, F. Microsatellite markers: what they mean and why they are so useful. Genet Mol. Biol. 39, 312-328 (2016).

[25] Squirrell, J. et al. How much effort is required to isolate nuclear microsatellites from plants? Mol. Ecol. 12, 1339-1348 (2003).

[26] Poland, J. A. & Rife, T. W. Genotyping-by-sequencing for plant breeding and genetics. The plant genome. Journal 5, 92 (2012).

[27] Chung, Y. S., Choi, S. C., Jun, T.-H. & Kim, C. Genotyping-by-sequencing: a promising tool for plant genetics research and breeding. Hortic. Environ. Biotechnol. 58, 425-431 (2017).

[28] Catchen, J. M., Amores, A., Hohenlohe, P., Cresko, W. & Postlethwait, J. H. Stacks: building and genotyping Loci de novo from short-read sequences. G3 (Bethesda) 1, 171-182 (2011).

[29] Beier, S., Thiel, T., Münch, T., Scholz, U. & Mascher, M. MISA-web: a web server for microsatellite prediction. Bioinformatics 33, 2583-2585 (2017).

[30] Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584 (2016).

[31] You, F. M. et al. BatchPrimer3: a high throughput web application for PCR and sequencing primer design. BMC Bioinforma. 9, 253 (2008).

[32] Yang, X. et al. Characterization of a global germplasm collection and its potential utilization for analysis of complex quantitative traits in maize. Mol. Breed. 28, 511-526 (2011).

[33] Perrier X., Jacquemoud-Collet J. P. DARwin software. 2006. 2015.

[34] Jombart, T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403-1405 (2008).

[35] Peakall, R. O. D. & Smouse, P. E. GenAlEx 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 6, 288-295 (2006).

[36] Liu, X., Tang, H., Li, D. & Hou, L. Genetic diversity of coconut cultivars in China by microsatellite (SSR) markers. Mol. Plant Breeding. 2, 83-91 (2011).

[37] Loiola, C. M. et al. Genetic relationships among tall coconut palm (Cocos nucifera L.) accessions of the international coconut genebank for Latin America and the Caribbean (ICG-LAC), Evaluated using microsatellite markers (SSRs) . PLoS ONE 11, e0151309 (2016).

[38] Martinez, R. T., Baudouin, L., Berger, A. & Dollet, M. Characterization of the genetic diversity of the Tall coconut (Cocos nucifera L.) in the Dominican Republic using microsatellite (SSR) markers . Tree Genet. Genomes 6, 73-81 (2010).

[39] Meerow, A. W. et al. Analysis of genetic diversity and population structure within Florida coconut (Cocos nucifera L.) germplasm using microsatellite DNA, with special emphasis on the Fiji Dwarf cultivar . Theor. Appl. Genet. 106, 715-726 (2003).

[40] Perera, L., Russell, J. R., Provan, J. & Powell, W. Use of microsatellite DNA markers to investigate the level of genetic diversity and population genetic structure of coconut (Cocos nucifera L.) . Genome 43, 15-21 (2000).

[41] Teulat, B. et al. An analysis of genetic diversity in coconut (Cocos nucifera) populations from across the geographic range using sequence-tagged microsatellites (SSRs) and AFLPs . TAG Theor. Appl. Genet. 100, 764-771 (2000).

[42] Lantican, D. V. et al. De novo genome sequence assembly of dwarf coconut (Cocos nucifera L.‘Catigan Green Dwarf’) provides insights into genomic variation between coconut types and related palm species . G3: Genes, Genomes. Genetics 9, 2377-2393 (2019).

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

[44] Lu, F. et al. Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol. PLoS Genet 9, e1003215 (2013).

[45] Liu, K. & Muse, S. V. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21, 2128-2129 (2005).

[46] Nei, M. Genetic distance between populations. Am. Nat. 106, 283-292 (1972).

[47] Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945-959 (2000).

[48] Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611-2620 (2005).

[49] Earl, D. A. & vonHoldt, B. M. Structure Harvester: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359-361 (2012).

[50] Jombart, T., Pavoine, S., Devillard, S. & Pontier, D. Putting phylogeny into the analysis of biological traits: a methodological approach. J. Theor. Biol. 264, 693-701 (2010).

[51] Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 27, 209-220 (1967).