A Simple Sequence Repeat (SSR) Marker Comparison of a Large In- and Ex-situ Potato Landrace Cultivar Collection from Peru Reaffirms the Complementary Nature of both Conservation Strategies

An enhanced understanding of the temporal dynamics of intraspecific diversity is anticipated to improve the adequacy of conservation priorities, methods and metrics. We report on the comparative genetic composition of ex- and in-situ landrace cultivar populations from a potato diversity hotspot in the Andes. A total of 989 landrace cultivars belonging to contemporary custodian-farmer in situ collections from central Peru were compared with 173 accessions from a spatially analogous, but temporally differential ex situ composite genotype reference (CGR) set using 15 nuclear microsatellite markers. A total of 173 alleles were detected, with 129 alleles (74.6%) being shared between both populations. Both populations contain exclusive allelic diversity with 32 and 12 unique alleles belonging to the ex- and in-situ population, respectively. The mean unbiased expected heterozygosity values of the ex- and in-situ population are very similar, 0.749 versus 0.727, with a slightly wider range and standard deviation encountered for the in situ population. Analysis of Molecular Variance shows that 98.8% of the total variation is found within both populations, while the fixation index (Fst = 0.01236) corroborates that the populations are not well differentiated. Surprisingly, only 41.0% of the ex situ population encounters a similar landrace cultivar in 23.4% of the in situ population at a non-stringent threshold similarity coefficient of 0.80. While the ex- and in-situ population under comparison show similarities and unique features at the allelic level, their landrace cultivar composition is surprisingly distinct. Results affirm that crop evolution is an ongoing phenomenon and that change in fixed geographies is occurring.


Introduction
The potato is a highly heterozygous outcrossing crop with ploidy levels ranging from the diploid to the pentaploid level.In the Andes, four cultivated species occur: (i) Solanum tuberosum (diploid, triploid and tetraploid Andigenum group); (ii) S. ajanhuiri, (iii) S. juzepzukii, (iv) S. curtilobum [1,2].The crop is characterized by high levels of intraspecific diversity with thousands of landrace cultivars occurring from 42° south (Chile) to 5° north (Colombia) of the equator.A landrace cultivar is a farmer-selected variety and a product of the natural and cultural environment in which it has developed [3].It has not been exposed to formal crop breeding schemes and is clearly distinct, uniform and stable in its characteristics and, when propagated by appropriate means, retains those characteristics [4].By far the highest landrace cultivar diversity is found within Solanum tuberosum, particularly the tetraploid Andigenum group, which constitutes 77% of the total accessions of potato landrace cultivars maintained ex situ by the International Potato Center (CIP).This collection has been intensively characterized and rationalized to reduce redundancy.
From the late 1920s till the early 1990s, potato collection expeditions by American [5], Dutch [6], English [7], Peruvian [8][9][10] and Russian [11] scientists, among researchers from many other nationalities, have contributed to the establishment of important ex situ collections that are currently maintained in gene banks around the world.The signature of the Convention of Biological Diversity (CBD) in 1992 and the adoption of national regulatory frameworks for access to genetic resources in Andean countries have led to a significant reduction of international potato collection missions.Nevertheless, these ex situ collections continue to provide valuable genetic diversity to potato breeding programs around the world in order to continue to deliver bred cultivars adjusted to a changing world in its climatic, as well as human dimensions.At the same time, interventions designed to support in situ conservation have increased throughout the Andes [12].In situ conservation of cultivated potato genetic resources in the Andes is a contemporary phenomenon, and diverse landrace cultivar pools continue to be actively managed by farmers [13,14].It is unlikely that present-day landrace cultivar populations and their inherent genetic diversity are an exact copy of what was found on-farm at the time of the collection missions.However, systematic comparisons of the structure of landrace cultivar populations in fixed geographical spaces along a well-defined timeline are scarce, as it is generally difficult to quantify change.A few cases, however, have been documented.For example, Monteros [14] showed that between 2008-2010, more potato landrace cultivars were found in the provinces of Chimborazo and Loja in Ecuador compared to what was reported by collecting missions in the 1970s and 1980s.Vigouroux et al. [15], comparing sorghum cultivars in the same region of Niger between 1976 and 2003, found no major changes in the main cultivars or in their genetic diversity.However, these authors observed a significant shift in adaptive traits.
Temporal variation in genetic structure and landrace cultivar composition in fixed geographical space is a likely consequence of complex biological or social change processes.Initial gaps in baseline data for timeline comparisons, such as species, landrace cultivar, hotspot or infrastructure biases [16], make comparative analysis difficult.The occurrence of potato landrace cultivars in direct proximity to wild relatives, somatic mutations and farmer selection practices are likely to have maintained ongoing diversification and evolution [17][18][19].Also, landrace cultivar loss and possible genetic erosion may have occurred as a result of extreme natural and human-induced events.Rapidly changing social and environmental conditions in the Andes, such as migration, increased farmer mobility and involvement in off-farm employment, diffusion of bred and commercial landrace cultivars and changing seed-tuber procurement practices, are also likely to play their part in the dynamics of landrace cultivar populations in fixed geographical spaces.Indeed, it is now well recognized that the loss of crop genetic diversity is a complex and context-dependent process conditioned by local differences in environment, culture and economy [20].However, potato landrace cultivar loss and genetic erosion are not directly correlated, since the potato is a highly heterozygous crop sharing alleles across geography and ploidy.Such genetic features of the crop make the present study particularly relevant in order to distinguish loss of genotypes from loss of alleles, which are usually confounded for autogamous crops.
Different molecular marker systems have been used to compare the population genetic structure of exand in-situ populations of crops and wild plant species [21][22][23].Simple sequence repeat (SSR) or microsatellite markers have also been used for this purpose, including for maize [24], beans [25] and grasses [26], and, indeed, are a robust tool to research the genetic structure and diversity of crop populations [27][28][29].SSR markers are ideal for paternity determination, identity verification, comparative genetics and population genetic studies, because of their high specificity, high polymorphism, good reproducibility and unambiguous scorability [30,31].For potato, they are widely used for research concerning genetic diversity [32,33], phylogeny and gene pool structure [2,34,35], as well as routine identification of duplicates [36].So far, SSR markers have not been used to compare exand in situ populations of potato landrace cultivars along a well-defined timeline.Here, we provide a benchmark comparison for a spatially restricted area within the potato"s center of origin and diversity.Systematic time-line comparisons of the population genetics and composition of potato landrace cultivars in well-defined hotspots within the potato"s center of origin will enhance our understanding of the temporal dynamics of intraspecific diversity.This, in turn, can help to inform conservation priorities, as well as adequate methods and metrics.Knowing whether the allelic and landrace cultivar in situ diversity has significantly shifted or remained relatively static within a 40-year timeframe may provide lessons for what crop conservation scientists consider to be geographically representative populations and how we conceptualize farmer-driven in situ conservation.Is the long-term conservation of landrace cultivars and inherent allelic diversity static or dynamic?What does temporal (non-)change teach us about spatial dynamics and representativeness of diversity in fixed geographies?When we monitor in situ conservation, should we consider total, new or lost or unique diversity as indicators?Or a combination?These types of research questions are particularly important for the International Potato Center (CIP), as the institution has recently embarked on the establishment of long-term monitoring sites in selected hotspots throughout the potato"s center of origin and diversity with a view to document how agrobiodiversity is evolving in light of socioeconomic and environmental change.

Site Description
Huancavelica is one of Peru"s 24 departments (19 of which produce potato).The department covers 22,131 km 2 and 80% of its territory and is located between 3,000 and 4,500 m above sea level.The total annual potato cropping area covers around 12,000 ha and is comprised of more than 86,000 individual farm units [37].Huancavelica is a well-recognized hotspot of intraspecific diversity of the potato with a total stock of well over 400 distinct landrace cultivars reported to be present in its territory [38].The highlands of Huancavelica used to be very isolated, but have opened-up as peace returned to the countryside in the late 1990s, and newly paved roads were constructed between 2003 and 2010.The earliest and most voluminous CIP collection from Huancavelica was established between 1972 and 1974 (n = 344; 64.8%), with accessions being added to the ex situ collection in 1985 (n = 7; 1.3%) and 1992-1994 (n = 180; 33.9%).Today, it consists of 143 accessions kept as tubers and 108 accessions conserved as in vitro materials for long-term storage and distribution.

Materials
DNA was extracted from leaf samples from a total of 989 landrace cultivars representing the total diversity conserved by 8 custodian farmers in their family stocks.These custodian farmers are individuals with large family heritage collections of landrace cultivars.Individual farmers contributed a minimum of 57 and maximum of 175 landrace cultivars.All together, they provided the total landrace cultivar sample defined as the in situ population for the purpose of the research reported here.The in situ population covered 3 of Huancavelica"s 7 provinces and 5 of the department"s 565 farmer communities (comunidades campesinas) along a north-south pseudo-transect.The participating custodian farmers were well recognized by their home communities as having particularly rich and representative landrace cultivar collections.Landrace cultivars were established in on-farm trials, with single replicate plots containing 5 to 10 plants per landrace cultivar.A single DNA sample was taken from one representative plant from each farmer-recognized landrace cultivar.The total in situ population has been molecularly characterized with 18 SSR markers [13].CIP routinely uses a well-defined composite genotype reference (CGR) set consisting of 742 accessions for comparative purposes.It concerns a representative subpopulation containing 16% of the landrace accessions of the gene bank collection with adequate geographical, ploidy, species and genetic coverage.The CGR set has been molecularly characterized with 50 different SSR markers [2].We used a spatially restricted subset (the ex situ population) of 173 landrace cultivar accessions from the department of Huancavelica and its 3 neighboring departments (Lima, Juní n and Ayacucho) for the purpose of comparing a spatially well-defined and representative ex situ population with a contemporary in situ population.A wider geospatial coverage for the ex situ sample was chosen to compensate for the fact that the original gene bank collection had been rationalized in the early 1990s.Duplicate accessions from Huancavelica that were eliminated from the CIP gene bank collection as part of this exercise are partially captured by using an expanded geographical range for the representative ex situ population.Landrace accessions from neighboring departments were also included, because past seed-tuber exchange is likely to have resulted in considerable overlap between the landrace cultivar portfolios of these departments (Figure 1).Both the exand in-situ populations are composed of diploid and polyploid landrace cultivars (2n = 2x = 24 to 2n = 5x = 60) belonging to multiple cultivated species and cultivar groups.

Methods
Standard procedures as practiced at CIP were applied, including DNA extraction with DNeasy 96 plant kits, high throughput genotyping with a LI-COR 4300 DNA Analysis System and SSR allele scoring with SAGA Generation 2 software (LI-COR).A genetic identity kit with highly polymorphic SSR markers was used for genetic fingerprinting [36].A primer combination was chosen to cover the whole genome and for high Polymorphism Index content [39].The alleles were scored as present (1) or absent (0), while missing data were scored as ""9"".The total fingerprinted in situ population (n = 989) was compared with the spatially restricted ex situ CGR set (n = 173) from central Peru.The 989 and 173 landrace cultivars from the inand ex-situ populations had good quality data (<0.3%missing data).Data from 15 SSR primer combinations were shared between both populations and used for comparative analysis (Table 1).The ploidy level of most landrace cultivars was determined applying microscopy and/or flow cytometry with a Partec® ploidy analyzer.The complete dataset is available as supplementary data and can be accessed from the website [40]

Statistical Analysis UPGMA
Dissimilarity analysis for the populations was conducted using the Jaccard coefficient and UPGMA clustering method applying NTSYS-pc 2.2 software [41].Principal Component Analysis (PCA) was conducted, and dissimilarity trees (dendrograms) were constructed using an unweighted neighbor joining (NJ) clustering method for a dissimilarity matrix calculated with the Jaccard"s coefficient using DARwin 5.0 software [42].Unbiased expected heterozygosity (He u ) values were calculated for both populations and all SSR markers.Additionally, relative allele frequencies for the inand ex-situ populations were calculated.The population genetic structure of the spatially restricted ex situ and contemporary in situ populations were compared applying Analysis of Molecular Variance (AMOVA) with Arlequin 3.5.1.3software [43].

Results
A total of 173 alleles were detected within the two populations under comparison with a total of 161 and 141 alleles encountered in the exand in-situ populations, respectively.A total of 129 alleles occurred in both the exand in-situ populations, while 32 and 12 alleles belonged uniquely to the exand in-situ populations, representing 19.9% and 8.5%, respectively, of the total number of alleles accounted for within each population.In all cases, the largest proportion of alleles detected with individual SSR markers were shared between populations (Figure 2).Most markers detected unique alleles for both the exand in-situ populations.Five individual markers detected unique alleles restricted to one of the populations: STI0030, STI0032, STM1052, STG0006 and STI0022.Most of the alleles detected in the exand in-situ population were frequent (Table 2).The in situ population is characterized by a high presence of rare alleles indicating uneven distribution patterns in farmer"s landrace cultivar stocks.All of the in situ population"s unique alleles were rare, with an average frequency of 0.3% (minimum 0.1%, maximum 0.7%).A high proportion of scarce rather than rare alleles within the ex situ population alludes to a comparatively more even allele distribution among different landrace cultivars.Most of the ex situ population"s unique alleles were scarce, with an average frequency of 1.3% (minimum 0.6%, maximum 5.2%).Allelic diversity was slightly higher for the ex situ as compared to the in situ population, with average unbiased expected heterozygosity (He u ) values of 0.749 and 0.727, respectively.The difference of the allelic diversity for these two populations was non-significant (Figure 3).He u values for 13 out of 15 SSR markers were higher for the ex situ as compared to the in situ population, with differences ranging from a minimum of 0.011 for marker STM5127 to a maximum of 0.056 for marker STI0023.He u values for two SSR markers were higher for the in situ population, with a difference of 0.001 and 0.023 for markers STGBSS and STM0031, respectively.
AMOVA provides additional insights into the molecular variation of the exand in-situ populations.The Fst value of 0.01236 indicates that the populations are not very well differentiated.In other words, the two populations share most of the allelic diversity.Most variation is encountered within rather than among populations, reaffirming that most alleles are shared between populations with comparable allele frequencies.Unique alleles of each population were either rare or scarce, partially explaining the modest percentage of variation of 1.24% found among the two populations (Table 3).Beyond allelic structure and diversity, a comparative analysis was conducted for the number and percentage of landrace cultivars with equivalents shared between both populations at different threshold levels as defined by the Jaccard similarity coefficient (Figure 4).Considering a stringent threshold level equivalent to a similarity coefficient of 1.00, a total of 16 accessions from the ex situ population (n = 173) match with 20 putatively identical landrace cultivars from the in situ population (n = 989).Only very few landrace cultivars are shared between the two populations, with 157 out of 173 and 969 out of 989 cultivars relative to the total sample size belonging uniquely to the exand in-situ population.At a non-stringent threshold level, considering a similarity coefficient of 0.80, a total of 71 accessions from the ex situ population match with 231 shared landrace cultivars from the in situ population.In other words, at this level, 41.0% of the ex situ population encounters a relatively similar landrace cultivar in 23.4% of the in situ populations.Even at this non-stringent level, 102 out of 173 and 758 out of 989 cultivars, respectively, belong uniquely to the exand in-situ population.Allelic diversity, which is slightly higher within the ex situ population, is, thus, not automatically matched by a higher percentage and number of unique landrace cultivars when compared with the in situ populations.Although most alleles are shared between the two populations, the same is not true at the level of landrace cultivars.Most landrace cultivars are actually unique to either the exor in-situ populations.In terms of absolute numbers, the in situ population holds most unique landrace cultivars, even though, as expected, dissimilarity analysis within each population indicates that the percentage of similar landrace cultivars within the in situ population is much higher compared to the ex situ population: 7.5 versus 8.2%, 22.5 versus 35.9% and 31.2%versus 64.2% at a similarity coefficient of 1.00, 0.90 and 0.80 of the exand in-situ population, respectively.In other words, at a similarity coefficient of 1.00, 0.90 and 0.80, the exand in-situ populations, respectively, contain 160 and 908, 134 and 634, 119 and 354 unique landrace cultivars.
Figure 5 shows a comparison of the exand in-situ populations based on a Principal Component Analysis (PCA) and dissimilarity analysis with DARwin software.Both graphs visually reaffirm that the in situ population consists of numerous unique landrace cultivars that do not overlap with landrace accessions from the ex situ population.A comparison of each sub-branch of the dissimilarity tree at similarity coefficients of 0.4, 0.5 and 0.6, respectively, evidences that 21 (2.1%), 335 (33.9%) and 570 (57.6%) out of the 989 landrace cultivars from the in situ population are located in isolated subgroups with no accession from the ex situ population present in the sub-branch.This analysis also shows that five out of 39 (12.8%), 55 out of 136 (40.4%) and 126 out of 239 (52.6%) subgroups existing at similarity coefficients of 0.4, 0.5 and 0.6, respectively, are unrepresented by an accession from the CGR set or ex situ population.In other words, selected areas of the PCA graph and branches of the dendrogram are exclusively composed of sub-populations of the overall in situ population, pointing to the presence of unique landrace cultivars that do not have a near-similar match in the ex situ population.

SSR-Based Population Genetics
A temporal comparison of the populations separated over a near 40-year time period shows relative stability in the sense that three quarters of the alleles are shared between both populations, and 98.8% of variation based on AMOVA is encountered within rather than among populations.The outcome of this time-line comparison is similar to findings from localized spatial analysis in the Andes, which have shown that the majority of allelic variation is contained within geographies and landrace populations [46].It reaffirms that potato landrace cultivars do share a large proportion of alleles across species, ploidy, space and time.This can be attributed to common parentage, interspecific hybridization and seed-tuber exchange [2,[47][48][49].
The fact that most unique alleles are either scarce (ex situ) or rare (in situ) suggests modest genetic change at the margins of the allele frequency distribution pattern.This, in turn, can also be interpreted as evolution with an influx of "new" and efflux of "old" alleles.A total of 32 alleles, or 19.9% of the total number of alleles detected in the ex-situ population, were apparently not detected in situ, indicating that genetic erosion sensu stricto within a geographical confined space might be occurring at the allele level.Strictly speaking, ex situ conservation is not completely static.Seed-tuber and in vitro conservation needs regular regenerations.Therefore, management practices and somaclonal variations could potentially be additional sources of genetic change of the ex situ conserved population.Another factor that may explain some of these differences is the possible mixture of some landrace cultivar accessions during long-term ex-situ conservation.Clearly, the detection of exclusive alleles ex-situ reaffirms the importance of conservation in gene banks.The fact that the comparison of the exand in-situ populations resulted in a non-significant difference of the allelic diversity and representative allelic diversity coverage as shown by the low fixation index (Fst value) suggests that the spatially restricted CGR set is a genetically representative subset of the overall gene bank collection.
Most of the unique alleles that are now restricted to the ex situ population are scarce, and if taken as a mirror of the original base-line situation for the in situ population in the early 1970s, would consequently have been more prone to loss compared to moderate or frequent alleles.Causes of localized allele loss are likely to involve landrace cultivar loss, which in turn, may be a consequence of multiple underlying causes, ranging from pest and disease pressures to socioeconomic change [50].Interestingly, 12 unique alleles, representing 8.5% of the total number of alleles detected in the in situ population, were encountered exclusively within the contemporary in situ population, evidencing the origination of "new" allelic diversity within the pre-defined geographical space of Huancavelica as compared to the base-line situation reflected by the spatially restricted CGR set with landrace cultivars from the exact same region and direct neighboring departments.This finding reaffirms the importance of in situ conservation to allow for seed-tuber flows and ongoing evolution.The dynamic nature of in situ conservation can be expected to have resulted in significant changes of the landrace cultivar composition during several decades of gene flow, mutations and anthropogenic selection [51][52][53][54].This might also have involved the displacement of landrace by bred cultivars, as these have been grown in Peru since the early 1950s.Enhanced integration, road building and temporal migration will surely also have facilitated seed-tuber exchange between formerly isolated regions and, consequently, can also offer a partial explanation for the "new" allelic diversity encountered [55,56].
Conservation measures to avoid genetic erosion should build on an enhanced understanding of socioeconomic and environmental change and its impact on loss or enrichment [57].It may include such diverse areas as youth engagement, seed system innovations, market linkages and farmer cuisine [58].Links with gene banks through regular checks of exagainst in-situ populations (and vice versa) may also help prioritization of appropriate conservation interventions, such as the incorporation of rare landrace cultivars into gene banks or the reintroduction of "lost" landrace cultivars when requested by farmer communities.

SSR-Based Overlap of Landraces
The unexpected finding that most landrace cultivars do not have shared equivalents between the contemporary in situ and spatially restricted CGR set or ex situ population may have several explanations.It may be a direct consequence of the population dynamics of landrace cultivars, which is known to involve such diverse processes as sporadic sexual reproduction, adaptive farmer management, chimerism, rural change and technology options promoted by development interventions [54,59,60].Farmer seed-tuber exchange of landrace cultivars in Huancavelica is typically very limited during normal cropping seasons.However, during periods of seed shortages following extreme weather events, such as hail, frost and drought, seed flows tend to be more voluminous and offer a partial explanation for cultivar change [61].Peru"s central highlands were severely affected by political instability and rural violence between 1981 and 1994 [62].This conflict disrupted agricultural production and is likely to have had a significant impact on the dynamics of cultivar portfolios.
The CGR set from Huancavelica and its direct neighboring departments is likely to have underrepresented the total coverage of landrace cultivars, because by definition, a composite genotype collection is a reduced subset of the complete ex situ gene bank collection.So-called cosmopolitan landrace cultivars that have a wide geographical distribution range generally only have a small number of representatives in the gene bank and may, therefore, have been unrepresented.The rationalization of CIP"s gene bank was conducted without bias for provenance, such that landrace cultivars collected in Huancavelica were eliminated from the collection if they were duplicates with samples from another department.Furthermore, collection gaps may represent an additional cause of the limited overlap of shared landrace cultivars encountered between the exand in-situ populations.This is partially a consequence of the tendency that unique genotypes are microgeographically concentrated [46].Recent in-depth studies at the sub-regional level also tend to sample more extensively compared to earlier potato collectors [13,14].
It is clear that significant allelic coverage between populations is not equivalent to adequate landrace cultivar coverage.The bulk of landrace cultivars are restricted to either the exor in-situ population: 102 out of 173 and 758 out of 989, respectively, at a non-stringent similarity coefficient of 0.80.This has interesting implications for the conservation ecology of the vegetatively propagated potato in its center of origin.The finding that only 41.0% of the ex situ population encounters a similar landrace cultivar in 23.4% of the in situ population at a non-stringent threshold similarity coefficient of 0.80 clearly suggests that the population dynamics of potato landrace cultivars in confined mountain geographies may be much more dynamic than previously thought.How much of this dynamism over a near 40-year timeline can be assigned to a "natural" Darwinistic type of selection versus anthropogenic change processes remains unresolved.Yet, it is well-known that the population dynamics of landrace cultivars as compared to wild species is to a much larger extent explained by human factors [63].A high level of dynamism for landrace cultivars as compared to their inherent allelic diversity also has implications about what we consider as representative for a region based on georeferenced passport data from gene banks.Ex situ landrace cultivar populations from a spatially well-defined hotspot may only represent a partial snapshot of what a situation was at a predefined moment in time.

Conclusions
The SSR marker comparison of a large spatially restricted exand in-situ population along a well-defined timeline shows that 74.6% of the detected alleles are shared between both populations, and only 1.2% of the molecular variation is encountered among rather than within populations.Overall, the population genetic structure in terms of allele presence versus absence and allele frequency variation can be considered modestly stable.However, out of the total number of alleles detected, both the exand in-situ populations also contain 19.9% and 8.5% of unique alleles.These unique alleles are either scarce or rare and seem to represent the tail ends of an evolutionary continuum with a simultaneous influx and efflux of alleles.
The "loss" and origination of "new" alleles within a defined timeframe and confined geographical space do occur in parallel and, thus, provide evidence for the occurrence of ongoing crop evolution within in situ conservation hotspots.The relatively stable but at the same time open nature of the population genetics of the potato in its natural habitat, resulting in both loss and enrichment, reinforces the complementary nature of exand in-situ conservation strategies.Ex situ conservation will ensure the availability and use of unique genetic diversity that may at present be non-adapted in its original collection site, while continued in situ management will result in responsive evolution favoring the presence and frequency of alleles to facilitate adaptation to changing environments.
Low levels of overlap of equivalent landrace cultivars shared between the exand in-situ populations at stringent (1.00, 0.95) to non-stringent (0.90, 0.80) threshold levels of the Jaccard similarity coefficient suggests that the population dynamics of potato landrace cultivars over a near 40-year time period is much more dynamic than anticipated, particularly more so than changes observed in the population genetics at the allelic level.Both the exand in-situ populations predominantly contain unique landrace cultivars, but particularly the contemporary in situ population is characterized by the presence of numerous unique landrace cultivars.A model for more systematic, rigorous and regular monitoring of the landrace cultivar population dynamics in well-defined in situ hotspots is needed to better understand the drivers of change and evolution in different social and environmental settings.Such long-term monitoring would ideally involve interdisciplinary teams of biophysical and social scientists collaborating toward an enhanced understanding of the links between socioeconomic, environmental and crop evolutionary change processes.

Figure 1 .
Figure 1.Geospatial coverage of the samples composing the inand ex-situ population.

Figure 2 .
Figure 2. Absolute number of alleles shared between the exand in-situ populations and absolute number of alleles restricted to each population detected with each of 15 SSR.

Figure 3 .
Figure 3. Unbiased expected heterozygosity mean (He u ; CI 95%) considering all microsatellite loci for the exand in-situ population.

Figure 4 .
Figure 4. Comparison of the percentage (%) of shared equivalent landrace cultivars between the exand in-situ populations based on different threshold values of the Jaccard similarity coefficient.

Table 1 .
Coverage of the simple sequence repeat (SSR) markers applied.

Table 3 .
Comparison of the ex and in-situ populations based on Analysis of Molecular Variance (AMOVA).Fst, fixation index.