First Glimpse on Spring Starflower Domestication

The cultivation and domestication of plants are human-driven processes that change the biology and attributes of a plant. Ipheion uniflorum is a bulbous geophyte known as Spring Starflower whose cultivation dates back to the first half of the 19th century. At least seven cultivars have been developed from natural stands. However, comparative analyses of wild and cultivated materials are largely missing. In the present study, we provide a morphological evaluation and analyses of the cytological and genetic variability of I. uniflorum that reveal significant levels of differentiation and evidence of artificial selection in the Spring Starflower. Distinctive phenotypic characters in cultivated materials that are rarely found or lacking in wild plants and natural populations, such as pink or violet flowers, together with its reduced heterozygosity and starting genetic differentiation support the view of early mechanisms of domestication acting upon Spring Starflower plants. The probable geographic origin of the cultivated forms is discussed together with perspectives for plant breeding.


Introduction
Domestication is the result of an evolutionary process in which wild species are exposed to new selective environments related to human cultivation and use, partly resulting in the increased fitness of the domesticated species [1]. Thus, domestication is a co-evolutionary, mutualistic process, in which domesticated species are transformed, adapted to human control, and propagated for human benefit, while humans start to become dependent on their domesticates. As domesticated plants became mostly cultivated and established from a relatively small number of founders, there are strong demographic effects (e.g., reduced effective population size, loss of diversity), and effects directed by the conscious or unconscious selection of traits that leads to genetic changes and dramatic impacts on genomic architecture [2]. Such changes drive the acquisition of unique adaptations relative to wild progenitors.
There are approximately 2000 semi-and fully domesticated plant species from 140 taxonomic families [3,4]. Most are relatively young, in some cases existing for just a few centuries. Most domesticated plant species grown for aesthetic enjoyment entered cultivation after 1500 AC [5]. Unlike domesticated crops, such as rice or grapes, in which domestication has been a protracted process driven in large part by unconscious selection of traits, such as seed size, seed retention (non-shattering), loss of seed dormancy or synchronous germination [6,7], and long and gradual periods of effective population size reduction, ornamental domestication was initiated more recently and is largely driven by the conscious selection of aesthetic traits, such as leaf or flower color, size or shape, through human choice and design [1,8]. In the case of Ipheion uniflorum (Graham) Raf. (Amaryllidaceae, Allioideae), ornamental forms and cultivars have been developed that resulted in the acquisition of new phenotypic traits and are widely known as Spring Starflower.

Plant Materials
A total of 56 wild individuals of Ipheion uniflorum were used in the current study, including all individuals available from Sassone et al. [19], plus 5 individuals of different cultivated origin. Following the geographic differentiation within I. uniflorum previously detected [19], we structured our sampling of individuals into three geographic regions (South, Centre and North; see Figure S1) and two conditions (natural and cultivated). In addition, to contrast our sampling against the natural distribution of I. uniflorum and identify possible sampling bias, we downloaded the available information in GBIF [under the name: I. uniflorum (Lindl.) Raf] and kept only specimens bearing geographical coordinates (Pampean region, country codes: AR, BR and UY) [21]. The map was constructed in the R environment with the packages maptools v. 0.9-8 [22] and ggplot2 v. 3.3.3 [23].
Among the cultivated materials, two individuals were sampled in a private garden of Los Andes, Chile (SSF1 and SSF2), which were established from bulbs imported from England. One of these plants was identified as I. uniflorum 'Rolf Fiedler' (SSF2). A third individual was collected in a garden in San Isidro, Buenos Aires, Argentina (SSF3). The last two (SSF4 and SSF5) specimens were cultivated at the Botanical Garden and Botanical Museum in Berlin, Germany. All cultivated samples were compared against all 56 wild individuals from 10 collection sites (corresponding to 13 natural populations) distributed across the natural range of the species ( Figure S1). Both the cultivated and wild individuals were kept alive in a glasshouse at the Darwinion Institute of Botany (Buenos Aires, Argentina) and the Leibniz Institute of Plant Genetics and Crop Plant Research (IPK, Gatersleben, Germany). Voucher specimens from each collection site were stored at SI.

Morphological Assessment
Morphological exploration focused on identifying traits likely selected in cultivated plants and of potential interest for horticulturalists. Flower color, shape of the tepal apices and the size of the flower were studied in individuals from the natural populations. We also collected extra information about the analyzed traits from an exhaustive review of the literature [9,11,12,16] and commercial websites marketing the cultivars. In addition, flower color variability and tepal apices were also evaluated in ca. 1400 iNaturalist records (https: //www.inaturalist.org/, accessed on 27 December 2021), 60 records from wild materials and ca. 1340 from cultivated materials (>1000 from the U.S.A., ca. 300 from Europe, and the rest from Australia, Japan, New Zealand and South Africa). The collected data are summarized in Tables S1 and S2, including the codification of the morphological characters.

Karyological Analyses and Ploidy Estimations
Fresh root tips were collected and treated for 6 hr at room temperature with a solution of 0.05% colchicine, fixed for 12-48 hr in ethanol/acetic acid (3:1, v/v) and stored in 70% ethanol at 4 • C. Root tips were then hydrolysed for 10 min at 60 • C in a solution of 1 m HCl, and stained for 20 min or more using Schiff's reagents (Feulgen reaction; Sigma-Aldrich, St. Louis, MO, USA). Root meristems were macerated in a drop of 2% acetic orcein and analyzed using a Leica DM5500B microscope (Leica Microsystems GmbH, Wetzlar, Germany). Chromosome measurements were performed using three to six selected mitotic metaphases with the Leica Application Suite software (LAS version 4.1.0, Leica Microsystems AG, Heerbrugg, Switzerland). Chromosomes were classified as metacentric (r = 1.00-1.49), submetacentric (r = 1.50-2.99), acrocentric (r ≥ 3.00) or telocentric (r = ∞), where r = length of the long arm/length of the short arm.
Karyological differentiation among individuals was inferred by the mean length of chromosome complements and karyotype symmetries. For possible interchromosomal asymmetries, the A2 index [24] was used. For intrachromosomal asymmetries, the total form percentage (TF%; [25]), the karyotype asymmetry index percentage (Ask%; [26]), the symmetric index (Syi; [27]) and the intrachromosomal asymmetry index (A1; [24]) were estimated. In addition, Stebbins' categories [28] were added for the characterization of intrachromosomal asymmetries. For each asymmetry index, the variation between mean values of natural and cultivated individuals was contrasted statistically by ANOVA followed by Tukey's test.
Ploidy estimations were based on cytological analysis (above) and by measures of genome-wide heterozygosity [19]. Briefly, the R package "gbs2ploidy" [29] was used to infer ploidy of individuals based on the allelic ratios of heterozygous single nucleotide polymorphisms (SNPs) during variant calling. Cases of inferred ploidy were corroborated using samples with previous cytological data and known ploidy.

Genotyping by Sequencing (GBS) Analyses and Population Genomics of Natural and Cultivated Materials
DNA extraction, library preparation and sequencing were carried out as described in Sassone et al. [19]. GBS loci were assembled using the ipyrad version 0.9.77 pipeline [30]. The different parameters were defined after exploring alternative settings following the recommendations of Eaton et al. [31,32]. Sequence reads for the GBS Illumina runs are stored in the European Nucleotide Archive (ENA, http://www.ebi.ac.uk/ena/data/view/ PRJEB43919, accessed on 29 December 2021).
A discriminant analysis of principal components (DAPC, [33]) was carried out using the detected number of clusters of natural populations (three groups [19]). The DAPC first performs a principal component analysis (PCA) to summarize genotypic variation among individuals, and then a discriminant analysis to categorize the PCA results, thus maximizing the variation among a predefined set of groups and minimizing variation within them. The DAPC was also used to identify the genomic similarity of the five cultivated specimens to each other and to the samples collected in nature. The analyses were performed using the R package 'adegenet' [34] and confirmed through a cross-validation test with 1000 permutations. The first three principal components (48.46% of variance conserved) of the PCA and three discriminant eigenvalues were retained. These values were confirmed by cross-validation analysis. The genetic diversity among accessions was also assessed with the phylogenetic network using the neighbor-net method [35]. The Bayesian clustering method STRUCTURE version 2.2.4 [36] was used to determine the number of distinct genetic clusters (K) and with a burn-in period of 500,000 followed by 2,000,000 repetitions, as implemented in the ipyrad API. Ten replicate analyses were performed with values of K = 1-10. Population genetic summary statistics (H O , H E , H T , F ST , and F IS ) were calculated to describe and compare overall and population-specific genetic diversity using the R package 'hierfstat' [37], and file conversions were facilitated with the R package 'dartR' [38]. We also used the R package 'hierfstat' to calculate the pairwise genetic distance (F ST ) for each of the subpopulations, according to Weir and Cockerham [39].

Morphological Characters within Ipheion uniflorum
When assessing qualitative characters among natural versus cultivated materials including the revised literature, an interesting bias in the distribution of characters was found. While natural Ipheion uniflorum plants generally have white flowers (in 86% of records, n = 43) and lilac or blueish flowers (these colors appear in some individuals within the populations, Tables S1 and S2), sky-blue-like flowers with rounded tepals have only been reported in two populations from Uruguay (Punta del Este and Piriapolis [11]) and they resemble the deep periwinkle blue flowers that are a distinctive trait of the I. uniflorum 'Rolf Fiedler' (Figure 1). This color is only found in Spring Starflowers. Likewise, pink and pinkish colors exhibited in the I. uniflorum 'Jessie' and 'Tessa' as well as the violet colors of I. uniflorum 'Tessa' were not recorded in our dataset (Tables S1 and S2). Among the studied materials, variation in tepal apices morphology (acute or apiculate, obtuse or round apices) is also found across natural populations, the obtuse tepal apices being the most widespread ( Figure 1). variation within the species from 2 to 3.5 cm length with a tepal fusion from 0.9 to (Tables S1 and S2). In addition, our assessment also indicates that the variatio morphological characters found among wild individuals is not significantly different that observed among cultivated individuals and show no extreme values except for plant size measurements (Table S2). Data from different studies and growers (Tab indicate that cultivated plants are taller, often reaching up to 30 or 40 cm height, wh wild plants sizes range from 6 to 20 cm height.  When assessing quantitative characters, the size of the flower exhibited a broad variation within the species from 2 to 3.5 cm length with a tepal fusion from 0.9 to 2 cm (Tables  S1 and S2). In addition, our assessment also indicates that the variation in morphological characters found among wild individuals is not significantly different from that observed among cultivated individuals and show no extreme values except for a few plant size measurements (Table S2). Data from different studies and growers (Table S2) indicate that cultivated plants are taller, often reaching up to 30 or 40 cm height, whereas wild plants sizes range from 6 to 20 cm height.

Cytological Variability within I. uniflorum
Cytological information of I. uniflorum was gathered from this study and previous analysis ( [19] and references therein) for six natural and five cultivated specimens ( Figure 2; Table 1). Overall comparative analyses of ploidy, genome size and karyotype formula were highly stable among individuals, except for the finding of a polyploid individual among the cultivated materials (sample AS183; Table 1). Detailed karyological data were extracted for seven individuals from the South and Centre geographic areas, while the very few individuals available from the northern region and from cultivated materials did not provide enough good quality metaphases. The different chromosomal parameters support the cytological data of Table 1 and add sources of karyological variation between individuals from the South and Central regions.
The quantification of the heterogeneity in chromosome size through intrachromosomal asymmetry indexes showed little variation between individuals and regions, all cases displaying A1 = 0.99, falling in the 1C Stebbins category, and non-significant differences for TF%, Ask% and Syi indexes ( Table 2). The interchromosomal asymmetry index A2 also displayed non-significant differences with values close to zero for individuals in the southern (A2 = 0.077) as in the central (A2 = 0.093) areas, reflecting conservation among chromosome size in the karyotype. Yet, variation in some chromosomal parameters from individuals between regions was observed. The largest/smallest chromosome ratio (R) were more variable among individuals from the South (R = 4.0; range 2.3-4.9) in comparison to individuals from the central region (R = 3.4; range 3.1-3.6). Likewise, a wider variation was found in the mean chromosome ratio (long/short arms; r) between individuals from the southern (4.2 ± 0.91, min 3.1, max 4.7) compared to those from the central (4.0 ± 0.27, min 3.6 max 4.2) regions (Table 2), as well as the number of chromosome pairs with arm ratio > 2 (r > 2), ranging 2-4 and 2-3 between individuals from the South and the Centre, respectively.  ; c = mean chromosome length (μm); c max, c min = maximum and minimum chromosome length (μm); i = mean centromeric index; TF% = ratio between the total sum of short arms (p) and the total length of a chromosome set ×100; Ask% = ratio

Genomic Diversity within I. uniflorum
In total, 20,641 loci were retained after filtering, and the concatenated GBS matrix was 1,864,491 bp length. Bayesian clustering analysis in STRUCTURE based on 11,355 SNPs (2302 unlinked) suggested that the best grouping number following ∆K was K = 3 ( Figure S2a), and this is mostly in concordance with the 3 geographical areas identified in a previous study [19] ( Figure S2b).
The DAPC was performed with four a priori groupings (three corresponding to the geographically differentiated natural populations plus one more for cultivated specimens (Figure 1).
The results show that the genetic similarity among wild individuals from natural populations and cultivated specimens is congruent among methods (Figures 3 and S1). The cultivated materials SSF1and SSF2 show a genomic composition matching specimens from populations in the Centre of I. uniflorum distribution (LP, AS271, AS284; TA541, AS202) and reveals high levels of admixture with Northern populations (UR; AS545) (Figures 3a and S2b). Regarding SSF3, this genotype shows a genomic composition and level of admixture comparable to those found in Southern I. uniflorum populations (BB, AS135, AS138; MBu, AS260) (Figures 3a and S2b). The other two genotypes, the SSF4 and SSF5, show a distinct pattern with slight differences between methods. While in the STRUCTURE analysis these cultivated materials are grouped with wild individuals from Southern/Centre populations ( Figure S2b), the DAPC shows them as clearly differentiated from wild individuals in all populations (Figure 3a).
The neighbour-net analysis depicts high genetic diversity and a phylogenetic network concordant with the above results ( Figure 3b). The analysis also shows the genotype SSF1 segregated from individuals from the Centre and more connected to individuals from the South (SCa; AS236, AS298, AS297) (Figure 3b). The genotype SSF3 is connected to Southern populations, but also it appears close to one individual from the Centre (SA; AS188) (  Table 3), while with northern populations, the observed level of genetic differentiation was the lowest (F ST = 0.139; Table 3). ; c = mean chromosome length (µm); c max, c min = maximum and minimum chromosome length (µm); i = mean centromeric index; TF% = ratio between the total sum of short arms (p) and the total length of a chromosome set ×100; Ask% = ratio between the total sum of long arms (q) and the total length of a chromosome set ×100; Syi = ratio between the mean length of the short arms (p) and the mean length of the long arms (q) ×100. A1, A2 = intrachromosomal and interchromosomal asymmetry indices, respectively; R = largest/smallest chromosome ratio; r = mean chromosome ratio (long/short arms).  The neighbour-net analysis depicts high genetic diversity and a phylogenetic network concordant with the above results (Figure 3b). The analysis also shows the genotype SSF1 segregated from individuals from the Centre and more connected to individuals from the South (SCa; AS236, AS298, AS297) (Figure 3b). The genotype SSF3 is connected to Southern populations, but also it appears close to one individual from the Centre (SA; AS188) ( Figure 3b). Despite SSF4 and SSF5 remain segregated, they are far away connected to an individual from a population of the North of I. uniflorum distribution (UR; AS545) ( Figure 3b).

Discussion
Domestication is a process involving the accumulation of trait differences over time, and, therefore, domestication mechanisms are most evident and widely known among the earliest and intensively cultivated species [2]. Domestication has dramatic impacts on genomic architecture reflecting historical demographic transitions (e.g., small number of founders and decreased effective population sizes) and changes in gene flow (e.g., directed selection) [2]. In more recently or less intensively cultivated plants, such as ornamentals, trait differences to wild relatives are expectedly less evident, and the genomic changes driven by incipient domestication processes are less known. In this paper, we present an initial analysis aimed at understanding the effects of human cultivation and artificial selection on Ipheion uniflorum that suggest the action of early domestication mechanisms in the Spring Starflower.

Natural Variation within Ipheion uniflorum
Ipheion uniflorum is the most relevant species within the genus Ipheion. From an evolutionary viewpoint, the species shows the broadest distribution, largest genetic and morphological diversity, and has played a pivotal role in the origin of its sister species [19]. From an economic viewpoint, the species is a resource of the different cultivated types of Spring Starflower.
The species grows in the Pampas of South America, a region subjected to intensive agricultural activities. Pampa's grasslands have only been preserved in areas with unfavorable edaphic or climatic conditions for agriculture, such as hills (Ventania and Tandilia systems), floodplains (e.g., the Salado basin) or at rural roadsides [40].
The observed within (and among) population morphological variability (e.g., Figure 1, Table S1) is more relevant in traits such as flower color and size, plant robustness and the presence of bulbils [11,12], and has support in extant cytological and molecular diversity. Even though the species is karyologically and genetically differentiated from the sister species I. tweedieanum (Griseb.) Traub and I. recurvifolium (C.H.Wright) Traub, at the intraspecific level, no chromosomal structural changes between the studied I. uniflorum samples and minor differences in cytological and karyological parameters were found, particularly among individuals from the South and Centre of the distribution (Figure 2; Tables 1 and 2). In the genetic variation analysis, the use of Evanno's ∆K method and best-fit population models at K = 2 and 3 revealed three genetically differentiated regional groups that follow a latitudinal geographical distribution [19]. Individuals in populations from the South (BB, MBu, SCa), Centre (MP, LaPer, SA, AZ, TA), and North (UR, LP) of the species distribution show substantial differences in private loci and relatively low admixture among them [19]. Wild I. uniflorum individuals across populations showed values of heterozygosity (H O = 0.084) lower than expected (H E = 0.178), with overall levels of genetic diversity H T = 0.267. The fixation indices and inbreeding coefficients showed values of F ST = 0.334 and F IS = 0.527, respectively, with the strongest differentiation decreasing from south towards north [19]. Overall, the data support the view that the natural variability of I. uniflorum in morphological, karyological, and genetic attributes being geographically structured with a gradual latitudinal differentiation trend.

Evidence for Early Domestication Mechanisms in I. uniflorum
As mentioned above, the domestication of ornamental plants was initiated more recently than other intensively cultivated crops and is largely dependent on aesthetically driven human selection [1,8]. This implies that morphological changes in aesthetic traits by human choice and design and accompanying changes in the genomic architecture will show a variable incidence of the domestication syndrome depending on the cultivation and domestication history of each ornamental.
Ornamental forms and cultivars of I. uniflorum have a recent origin dating back to the first half of the 19th century. Despite the relatively recent cultivation start, I. uniflorum's cultivars have acquired new phenotypic traits that differentiate to a certain extent the cultivated materials from the wild ones. The most striking one relates to flower color and the shape of tepal apices. The cultivar I. uniflorum 'Rolf Fiedler' displays a deep periwinkle blue flower color in combination with rounded tepals apices that only bear a resemblance to individuals from two Uruguayan populations (but see comments in the next section). Similarly, the distinctive pink flower color of the cultivar 'Jessie', or the pink-to-violet flower color of the recently released cultivar 'Tessa' have not been observed among natural populations.
While in nature, plant sizes range from 6-20 cm height, individuals in cultivation often reach up to 25-40 cm. Phenotypic variation for quantitative genetic traits, such as plant height, is continuously distributed in natural populations and often responds slowly to artificial selection [41,42]. In our study, the tendency to bigger plant sizes observed in Spring Starflowers might suggest a visible effect of artificial selection on genes influencing vegetative growth. Of course, an effect of nutrient availability in soils used for cultivation cannot be discarded, but we can point out that soils in the geographic range where I. uniflorum grows naturally are among the richest worldwide [43]. So, if there is an effect of nutrient availability, we might expect it to be low. This assumption could be easily tested in future studies by growing and evaluating different accessions under the same environmental conditions.
Regarding the genetic architecture of the Spring Starflower materials, our population structure analysis shows the cultivated materials are integrated among the three regional genetic clusters ( Figure S1). Whereas the wild I. uniflorum populations show best-fit population models at K = 2 and 3, when adding the cultivated specimens, Evanno's ∆K shifts the best-fit population model to K = 3, enhancing the differentiation between North and Centre populations ( Figure S2a).
This agrees with the DAPC and neighbor-net analyses that place three of the cultivated genotypes in-between individuals from the different regions by holding variable levels of admixture and segregates the other two (Figure 3a,b). For these two cultivated individuals (SSF4 and SSF5), the DAPC analysis identifies a new genomic composition of sequences that differs from the explored wild populations (Figure 3a,b). This indicates substantial genetic differences between these two cultivated genotypes and the natural genotypes studied in this paper. Whether SSF4 and SSF5 are derived from an otherwise extinct population/genotype no longer found today remains an open question.
Genetic diversity parameters indicate that, by including five cultivated specimens  [19] to F ST = 0.213, and the inbreeding coefficient changes from F IS = 0.527 [19] to F IS = 0.636. On the one hand, the data suggest a slight reduction in the genetic diversity of the dataset after including the cultivated samples, likely due to inbreeding. On the other hand, the gain in overall genetic diversity and decay in fixation indices indicate some level of genetic differentiation, which agrees with the observed segregation of two Spring Starflower specimens from the natural I. uniflorum populations mentioned above.
The overall data on morphological traits and genome sequence diversity point to the occurrence of phenotypic responses to artificial selection in the cultivated genotypes and suggest the action of early domestication mechanisms as a probable side effect of human cultivation.
Among cultivated species, gradual declines in effective population sizes and genetic bottlenecks often occur during initial domestication, and a further reduction in genetic diversity is caused by artificial selection [1,2,44]. In our analysis, a substantial reduction in levels of heterozygosity in cultivated Spring Starflower relative to wild materials was found. While in domesticated crops it is possible to distinguish between a decrease in genetic diversity caused either by a genetic bottleneck (genetic drift results in genomewide reduction in diversity) or by artificial selection (expected to reduce diversity in domestication-associated genomic regions) (e.g., [45,46]), in our case it was not possible to draw a precise conclusion. Considering that we are not aware of any active genetic improvement program using the Spring Starflower, the reduced genetic diversity observed in cultivated samples is likely caused by genetic drifts initiated through pre-domestication human management of wild stands. A wider sampling, including cultivated materials from diverse collections and markets, is needed to properly evaluate the extent of such genetic bottlenecks due to reduction in effective population sizes and the contribution of artificial selection to the observed decrease in heterozygosity.

A Glance at the Origin and Future of the Spring Starflower
Our genomic analysis provides little resolution on the original populations from which Spring Starflower cultivars used in the present study were first harvested. Yet, at least two cultivated individuals (SSF1 and SSF3) can be assigned to natural populations of a geographic group. One showing a genomic composition grouping with populations in the Centre and Northern regions (SSF1) and the other one with populations in the Southern region (SSF3). Most importantly, the proposed Uruguayan origin of the I. uniflorum 'Rolf Fiedler' [11] is corroborated, which here shows high levels of admixture to a single Northern population of I. uniflorum (UR; AS545).
A point that deserves consideration and is recurrent in color profiling for plant biology studies [47] is the lack of a systematic evaluation for the color of flowers found within I. uniflorum individuals. The occurrence of natural variation of floral color and the natural changes in flower pigmentation during flower development, maturation and climatic conditions makes it difficult to obtain standardized measurements of flower colors, which would help scientists to manage biodiversity stands. Therefore, without a numerical classification of flower colors obtained, for example, by a colorimeter (e.g., [48]), the current categorization of flower color in natural and cultivated plants remains subjected to human interpretation bias, especially when stating the occurrence of colors in the lilac-blueish range.
From a plant breeding perspective, the summary statistics of Spring Starflower specimens (p < 0.001) show low levels of heterozygosity (H O = 0.051) compared to the expected value (H E = 0.090), and relatively high inbreeding (F IS = 0.340). However, these parameters might show bias due to the small sampling size. The pairwise fixation indices showed the strongest differentiation of cultivated specimens to populations in the South (F ST = 0.485), and the lowest genetic differentiation with populations in the North (F ST = 0.139). This pattern parallels the genetic differentiation trend observed among I. uniflorum natural populations [19]. Thus, wild populations of I. uniflorum might well work as reservoirs of morphological and genetic variability for the further development of Spring Starflower cultivars. The continued development of new cultivars (e.g., I. uniflorum 'Tessa') depict both the interest of horticulturalists and the appreciation of the public (four cultivars had won prizes of the Royal Horticultural Society, https://www.rhs.org.uk/plants/pdfs/agm-lists/agm-ornamentals.pdf, accessed on 27 December 2021). In an era of accelerated human-induced species losses, many natural populations are today threatened and may become extinct in the next few years if no conservation measures are taken. The adoption of measures to protect and conserve the remaining populations of I. uniflorum throughout the distribution of the species (perhaps considering geographic regions as separate genetic reservoirs) is of great importance for future breeding plans and the release of new cultivars.

Supplementary Materials:
The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/genes13020243/s1, Figure S1: Distribution map of the naturally occurring populations; Figure S2: Population structure analysis; Table S1: Subsampling of the morphological data; Table S2: Summary of quantitative and qualitative morphological data.