Intra- and interspecific diversity analyses in the genus Eremurus in Iran using genotyping-by-sequencing reveal geographic population structure

Abstract

Eremurus species, better known as ‘Foxtail Lily’ or ‘Desert Candle’, are important worldwide in landscaping and the cut-flower industry. One of the centers of highest diversity of the genus Eremurus is Iran, which has seven species. However, little is known about the genetic diversity within the genus Eremurus. With the advent of genotyping-by-sequencing (GBS), it is possible to develop and employ single nucleotide polymorphism (SNP) markers in a cost-efficient manner in any species, regardless of its ploidy level, genome size or availability of a reference genome. Population structure and phylogeographic analyses of the genus Eremurus in Iran using a minimum of 3002 SNP markers identified either at the genus level or at the species level from GBS data showed longitudinal geographic structuring at the country scale for the genus and for the species E. spectabilis and E. luteus, and at the regional scale for E. olgae. Our analyses furthermore showed a close genetic relatedness between E. olgae and E. stenophyllus to the extent that they should be considered subspecies within an E. olgae/stenophyllus species complex. Their close genetic relatedness may explain why crosses between these two (sub)species have been found in the wild and are exploited extensively as ornamentals. Last, current species identification, while robust, relies on flower morphology. A subset of seven SNPs with species-specific (private) alleles were selected that differentiate the seven Eremurus species. The markers will be especially useful for cultivar protection and in hybrid production, where true hybrids could be identified at the seedling stage.

References

[1] Wendelbo, P. Further notes on Eremurus (Liliaceae) in Afghanistan. Acta Hortic. Gotobg. 28, 57-63 (1966).

[2] Kamenetsky, R. & Rabinowitch, E. Flowering response of Eremurus to postharvest temperatures. Sci. Hortic. 79, 75-86 (1999).

[3] Gaggeri, R. et al. Towards elucidating Eremurus root remedy: chemical profiling and preliminary biological investigations of Eremurus persicus and Eremurus spectabilis root ethanolic extracts. J. Med. Plants Res. 9, 1038-1048 (2015).

[4] Tuzcu, Z. et al. Antioxidant, antimicrobial and anticancer effects of different extracts from wild edible plant Eremurus spectabilis leaves and roots. Int J. Clin. Exp. Med. 10, 4787-4797 (2017).

[5] Rossi, D. et al. (R)-(-)-Aloesaponol III 8-methyl ether from Eremurus persicus: a novel compound against Leishmaniosis. Molecules 22, 519 (2017).

[6] Aysu, T., Demirbaş, A., Bengü, A. Ş. & Küçük, M. M. Evaluation of Eremurus spectabilis for production of bio-oils with supercritical solvents. Process Saf. Environ. 94, 339-349 (2015).

[7] Eghtedarnejad, N. Building wooden panels glued with a combination of natural adhesive of tannin/Eremurus root (syrysh). Eur. J. Wood Wood Prod. 74, 269-272 (2016).

[8] Naderi Safar, K., Kazempour Osaloo, S., Assadi, M., Zarrei, M. & Khoshsokhan Mozaffar, M. Phylogenetic analysis of Eremurus, Asphodelus, and Asphodeline (Xanthorrhoeaceae-Asphodeloideae) inferred from plastid trnL-F and nrDNA ITS sequences. Biochem Syst. Ecol. 56, 32-39 (2014).

[9] Baird, N. A. et al. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE 3, e3376 (2008).

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

[11] Qi, P. et al. UGbS-Flex, a novel bioinformatics pipeline for imputation-free SNP discovery in polyploids without a reference genome: finger millet as a case study. BMC Plant Biol. 18, 117 (2018).

[12] Cheng, L. & Zhang, Y. J. Chromosome number and karyotype of Eremurus chinensis Fedtsch. J. Wuhan. Bot. Res. 11, 281-282 (1993).

[13] Ling, W., Lingji, W., Miao, M. & Aijun, Q. Karyotype analysis of Eremurus anisopterus. J. Shihezi Univ. 22, 417-418 (2004).

[14] Smarda, P. et al. Ecological and evolutionary significance of genomic GC content diversity in monocots. Proc. Natl Acad. Sci. USA 11, E4096-E4102 (2014).

[15] Wendelbo, P. Asphodeloideae: Asphodelus, Asphodeline & Eremerus. Flora Iran. 151, 3-31 (1982).

[16] Jalili, A. & Jamzad, Z. Red Data Book of Iran: A Preliminary Survey of Endemic, Rare and Endangered Plant Sepcies in Iran. (Research Inst of Forests and Rangelands, Tehran, Iran, 1999).

[17] Fan, J. & Ma, M. A study on the pollination ecology of Eremurus inderiensis. J. Sustain Dev. 1, 133-137 (2008).

[18] Ma, M., Fan, J. & Li, J. Pollination characteristics of the ephemeroid plant Eremurus anisopterus. Front Biol. 3, 315-319 (2008).

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

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

[21] Berthouly-Salazar, C. et al. Genotyping-by-sequencing SNP identification for crops without a reference genome: using transcriptome based mapping as an alternative strategy. Front Plant Sci. 7, 777 (2016).

[22] Melo, A. T. O., Guthrie, R. S. & Hale, I. GBS-based deconvolution of the surviving North American collection of cold-hardy kiwifruit (Actinidia spp.) germplasm. PLoS ONE 12, e0170580 (2017).

[23] Poland, J. A., Brown, P. J., Sorrells, M. E. & Jannink, J.-L. Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS ONE 7, e32253 (2012).

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

[25] Baker, J. G. Revision of the genera and species of Anthericeae and Eriospermeae. Bot. J. Linn. Soc. 15, 253-363 (1876).

[26] Doyle, J. J. & Doyle, J. L. Isolation of plant DNA from fresh tissue. Focus 12, 13-15 (1990).

[27] Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

[28] DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491-498 (2011).

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

[30] Peakall, R. & Smouse, P. E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28, 2537-2539 (2012).

[31] Perrier, X. & Jacquemoud-Collet, J. P. DARwin software. http://darwin.cirad.fr. (2006).

[32] R Development Core Team . R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2012).