Genetic diversity and population structure assessed by SSR and SNP markers in a large germplasm collection of grape
© Emanuelli et al.; licensee BioMed Central Ltd. 2013
Received: 26 October 2012
Accepted: 27 February 2013
Published: 7 March 2013
The economic importance of grapevine has driven significant efforts in genomics to accelerate the exploitation of Vitis resources for development of new cultivars. However, although a large number of clonally propagated accessions are maintained in grape germplasm collections worldwide, their use for crop improvement is limited by the scarcity of information on genetic diversity, population structure and proper phenotypic assessment. The identification of representative and manageable subset of accessions would facilitate access to the diversity available in large collections. A genome-wide germplasm characterization using molecular markers can offer reliable tools for adjusting the quality and representativeness of such core samples.
We investigated patterns of molecular diversity at 22 common microsatellite loci and 384 single nucleotide polymorphisms (SNPs) in 2273 accessions of domesticated grapevine V. vinifera ssp. sativa, its wild relative V. vinifera ssp. sylvestris, interspecific hybrid cultivars and rootstocks. Despite the large number of putative duplicates and extensive clonal relationships among the accessions, we observed high level of genetic variation. In the total germplasm collection the average genetic diversity, as quantified by the expected heterozygosity, was higher for SSR loci (0.81) than for SNPs (0.34). The analysis of the genetic structure in the grape germplasm collection revealed several levels of stratification. The primary division was between accessions of V. vinifera and non-vinifera, followed by the distinction between wild and domesticated grapevine. Intra-specific subgroups were detected within cultivated grapevine representing different eco-geographic groups. The comparison of a phenological core collection and genetic core collections showed that the latter retained more genetic diversity, while maintaining a similar phenotypic variability.
The comprehensive molecular characterization of our grape germplasm collection contributes to the knowledge about levels and distribution of genetic diversity in the existing resources of Vitis and provides insights into genetic subdivision within the European germplasm. Genotypic and phenotypic information compared in this study may efficiently guide further exploration of this diversity for facilitating its practical use.
The genus Vitis contains about 60 species, or more strictly ecospecies, since there are no genetic barriers within the genus and the species are inter-fertile. They have a primarily temperate zone distribution, occurring extensively in the Northern Hemisphere. The leading cultivated species by far is V. vinifera L. ssp. sativa (or vinifera), and its wild form V. vinifera L. ssp. sylvestris represents the only Vitis taxon naturally found in Europe. In contrast, numerous species of this genus are indigenous to North America and East Asia. Although these wild species are only peripherally used for human consumption, they are of great economic importance as a source for resistance breeding and as rootstocks for the highly susceptible V. vinifera. Since the beginnings of cultivation, desirable forms of the wild grapevine and spontaneous mutants within cultivated populations have been selected and preserved by vegetative propagation. Additional cultivars have been developed by both deliberate and spontaneous interspecific as well as intraspecific breeding . A significant reduction of genetic diversity in both cultivated and wild grapevines occurred when the phylloxera insect was brought to Europe from America about 150 years ago. European vineyards were saved from extinction by the introduction of several native American, non-vinifera Vitis species, which were used as rootstocks and for breeding disease resistant interspecific hybrids . Over the past few decades the cultivated grapevine has experienced another drastic reduction of diversity, resulting in the disappearance of old local varieties, and the increased focus of global wine companies on only a few major cultivars. Likewise, genetic variation of the wild V. vinifera species has diminished due to loss of natural habitat. On the other hand, in the recent past, many conservation programs of genetic resources have been conducted in grapevine growing countries. As a result, a significant number of minor varieties have been collected and preserved in field collections. However, due to the long time required for field experiments and the lack of information on genetic variation, research efforts that would facilitate the use of existing collections for crop improvement have not been as frequent as the conservation activities.
Molecular characterization is now the favored way to quantify variation within germplasm samples [3–5]. For instance, microsatellites (simple sequence repeats, SSR), because of their polymorphism, reproducibility, and codominant nature, have become the markers of choice for compilation, standardization and exchange of information concerning grapevine genetic resources . Recently SSRs were applied in several surveys of entire germplasm collections [7–9]. These studies provided a broader estimation of genetic diversity in each collection and found a high degree of clonal relationships, synonyms, homonyms, and curation errors. Similar conclusions on the naming accuracy were achieved in the analysis of the USDA grape germplasm collection, using a genome-wide SNP genotyping approach . The authors evaluated haplotype diversity, pattern of population structure and the decay of linkage disequilibrium in V. vinifera accessions, with a set of 5,387 SNPs. Results of the survey suggest that although substantial genetic diversity has been maintained in the grape following domestication, there has been a limited exploration of this diversity. Since it is still unclear to what extent these collections represent an unbiased sample of the worldwide genetic variation, further broad studies of grapevine germplasm are required, as well as the development of a manageable set of materials that will facilitate access to this variation. The present study is part of an effort to characterize and to dissect the genetic structure of one of the largest collections of grape germplasm in Europe, which maintains, amongst the 2700 accessions, hundreds of putative wild V. vinifera individuals and selections of post-phylloxera breeding materials. Our aim is to maximize the potential contribution of the collection dataset to the development of an international database and the creation of composite core collections. We applied the SSR descriptors chosen for the European Vitis Database , the SSRs employed to genotype the largest grape repository in the world , and 384 SNPs spread throughout the genome, which included the set of markers proposed for grapevine cultivar identification . This allowed us 1) to examine the level of genetic diversity, structure and differentiation within the germplasm collection, comparing the usefulness of different marker systems; 2) to sort out genetic core collections from the dataset of V. vinifera cultivars and contrast their genetic variation with that of a sample representative of the collection’s phenological variation, with the intention of justifying a contribution of these samples to association studies.
Results and discussion
Genetic characterization of a Vitisgermplasm collection
Level of redundancy and number of multilocus genotypes identified using 10 SSRs within the entire FEM germplasm collection and its four grape subpopulations
Different SSR genotypes
SSR genotypes represented by one accession
SSR genotypes represented by two or more accessions
Average number of accessions with identical SSR genotype
Summary statistics of genetic variation at 22 SSR loci and 353 SNP loci in the entire FEM germplasm collection and its four grape subpopulations
When considering the four collection subsets (Sativa, Sylvestris, Rootstocks and Hybrids), the diversity parameters were different, compared to those estimated for the total collection. The number of alleles ranged from 234 in Sylvestris to 412 in Rootstocks for SSR loci, and from 669 in Rootstocks to 704 in Sativa for SNP loci. The average MAF of SNPs in Sativa and Hybrids was similar to that calculated for the total collection (0.25), while in Rootstocks and Sylvestris this value was only 0.08 and 0.19, respectively. The average effective number of alleles for SSR loci ranged from 3.61 in Sylvestris to 8.19 in Rootstocks, while these numbers for SNP loci were from 1.15 in Rootstocks to 1.58 in Sativa. The expected heterozygosity estimates of the subsets ranged from 0.69 (Sylvestris) to 0.83 (Rootstocks) for the SSRs and from 0.11 (Rootstocks) to 0.34 (Sativa) for the SNPs. The fixation index F ranged from 0.01 and −0.01 in Hybrids to 0.12 and 0.09 in Rootstocks, for the SSRs and SNPs, respectively. The SSR markers have been used to characterize diversity among cultivars from most of the regions of grape cultivation. The data from those studies are difficult to compare, because the number of accessions and the marker loci used are very different. However, similar to our survey, the other studies showed that the SSR diversity within V. vinifera is very high. In different analyses, the number of alleles per locus in Sativa varied from 8 to 11 for sample sizes ranging from 58 to 366 individuals and HE ranged from 0.62 to 0.85 with an average of 0.76 [7, 8, 14–17]. For instance, Laucou et al. , using 20 SSRs on 2323 cultivated V. vinifera accessions (Sativa), revealed an average number of 16.9 alleles per locus (6–36) and the expected heterozygosity of 0.76. Likewise, although that analysis was intended for cultivated grapevines, the authors included wild grape accessions (Sylvestris), as well as accessions resulted from interspecific crosses with North American Vitis spp. (Hybrids and Rootstocks). The subset of Rootstocks in those studies also revealed the highest number of alleles (405) and the highest heterozygosity (0.86), in spite of their relatively small sample size compared to the subset of Sativa. The Sylvestris sample in this study  presented the lowest number of alleles (203) and heterozygosity (0.62). We have observed a similar trend in our survey: the set of Rootstocks displayed the highest number of different SSR alleles (412). Sativa was less diverse than the set of Rootstocks and more diverse than Sylvestris. Fixation index was low in Sativa and Hybrids collections (0.02 and 0.01, respectively), compared to Sylvestris and Rootstocks subsets (0.10 and 0.12, respectively). Diversity in Sylvestris was lower than in the cultivated grape, because of the small number of unique individuals available in the collection. This supports a previous hypothesis that suggested the scarcity of this endangered subspecies with small populations results in higher inbreeding rates . The lower number of SSR alleles in Hybrids was likely due to a low sample size, however the level of heterozygosity (0.8) is consistent with former observations [9, 18].
Previous studies on SNP variation in grapevine concerned mainly cultivated grapevines and reported a similar level of diversity. In the survey of 1573 SNPs from a group of 11 grape genotypes corresponding to nine ancient unrelated cultivars and two wild grapevines, the expected heterozygosity (HE) ranged from 0 to 0.66 with a mean value of 0.30 . Likewise, in a set of 48 SNPs from a sample containing 151 non-redundant cultivars, HE was 0.404 .
It was shown that level of diversity quantified by heterozygosity based on SNPs is around two times lower than that estimated for SSR markers [20, 21]. This potential disadvantage of SNP can be overcome either by using a large number of markers or by considering haplotypes structure for each locus instead of single SNPs [10, 19]. The differences between SNPs and SSRs in levels of genetic diversity result from mutational properties of these two marker types. Because of the nature of SNP markers, we observed a smaller proportion of rare alleles in the frequency distribution of the SNP data compared to the SSR data. The intermediate frequency alleles in the SNP loci could also be the consequence of ascertainment bias, which is the bias introduced when loci are identified in a small panel of accessions that do not represent the full genetic variation of a genus or species [22, 23]. Moreover, the current high throughput genotyping of SNPs is based on measuring the relative signal strength of two expected alleles. However, in a population different allele types may exist at any locus and this unknown or “null” allele can interfere with exact genotyping of the expected alleles . Since most of the SNPs used in this study were discovered in V. vinifera, this may explain the lower level of diversity in our dataset Rootstocks. Here the expected heterozygosity based on SNP genotypes was around seven times lower, compared to that from microsatellite loci (0.11 and 0.83 respectively). In contrast, several published studies indicate good transferability of SSR markers amongst Vitis species [25–29]. However, the SSR diversity may be an underestimation since sequencing of some microsatellite loci suggested that the polymorphism did not correspond only to a variation in the number of repeats, but also to changes in their architecture and the flanking regions with substitutions and long indels [29–34].
Construction of a genetic core collection of V. viniferasativa
The purpose of developing genetic core collections is to provide a restricted set of accessions, feasible to handle, and representing the genetic variability among individuals in a large source of germplasm. Genetic core collections were constructed to maximize the allelic diversity among Sativa accessions based on microsatellites, as these markers have been shown to provide greater information content compared to SNP markers. Based on the M-method, fifty eight cultivars (core G-58) were sufficient to capture all the 274 alleles occurring in more than 0.5% of the samples analyzed. The core G-58 was then used to design the final genetic core collection retaining 100% of SSR diversity, i.e. 362 alleles. The optimal size of this core was 110 individuals (core G-110), thus 52 accessions were added at this step to retain 88 rare alleles. The M-method sampling strategy showed a superior efficiency compared to random sampling. In particular core G-58 and core G-110 retained 45 (274 vs 229) and 101 (362 vs 261) more alleles compared to random cores of the same size (Additional file 2). These results show that only a small number of accessions are needed to retain the most frequent alleles as well as the whole allelic diversity (8% and 15% of cultivated grapevines in G-58 and G-110, respectively). A previous genetic core collection developed for cultivated grapevine by Le Cunff et al.  using the M-strategy required fewer individuals (92, i.e. 4%) to capture the total allelic diversity of the whole collection (326 alleles). We can assume that the high level of heterozygosity in grapevine is the major factor leading to capture all the genetic diversity with such a small number of individuals. Indeed, similar experiments have required 18% and 31% of individuals to retain the whole genetic diversity for A. thaliana and M. truncatula, respectively [36, 37].
Diversity of genetic and phenological core collections
Descriptive statistics of SSR and SNP diversity within phenological (Core P) and genetic (Core G-58, Core G-110) core collections with the percentage of alleles retained from the entire germplasm collection of cultivated grapevine (Sativa)
Random core G-58R §
Linkage disequilibrium between SNP loci was low (r2 < 0.2), which is consistent with previous findings [10, 19]. Such level of LD is unlikely to affect the analysis of population structure. Both SSR and SNP datasets were independently used for the model-based Bayesian clustering method as implemented in STRUCTURE. The most likely number of clusters (K) was evaluated considering the plateau criterion proposed by Pritchard et al. , using the non-parametric Wilcoxon test  and the ΔK method  (Additional file 3 and Additional file 4). For the SSR dataset the mean log-likelihood curve attained a maximum value around K = 6, beyond which the mean log-likelihood values reached a plateau and the standard deviations associated with the estimates increased. In contrast, for the SNP dataset the mean log-likelihood curve did not reach a plateau and the standard deviations did not increase drastically. The aspect of consistency among different simulations within each preset K can also be visualized through the similarity coefficient between different runs for each preset K according to Nordborg et al. . For both datasets the mean similarity coefficients among different simulations decreased for K larger than 5. The Wilcoxon test determined that best K is 6 for the SSR dataset and 5 for the SNP dataset. This test was also significant at K > 7 for the SNP data, reflecting a continuous increase of the likelihood values with respect to K. When more than five inferred populations are considered, no individual was strongly assigned (Q > 0.8 for SSRs or Q > 0.65 for SNPs) to the additional inferred populations. The ΔK criterion suggested by Evanno et al.  gave the highest value at two groups both for SSR and SNP loci, although peaks of ΔK were found also at K = 3 (for both SSRs and SNPs) and K = 6 (for SSRs only). This method is known to give rise to the first structural level in the data  and in the present study has led to discriminate the Vitis non-vinifera accessions (used as rootstocks) from the V. vinifera accessions (Sativa and Sylvestris), in agreement with the results displayed by PCA (Figure 2).
Hybrids are interspecific selections developed by crossing wild American species resistant to phylloxera, downy mildew or powdery mildew with European grapevine varieties. Several backcrosses with V. vinifera cultivars were required, especially for direct producer hybrids (ungrafted) to obtain superior wine grape cultivars. Accordingly, in the present study, Hybrids showed mixed ancestry with high admixture proportions of Sativa, ranging from 0.6 to 0.8.
Both SNPs and SSRs performed well in detecting the main subpopulations (K = 3): V. vinifera (Sativa and Sylvestris) and Rootstocks (Vitis ssp.). Except for K = 2, where both marker types showed a high percentage of individuals assigned to populations, all assignment percentages for the SNP dataset were lower than for the SSR dataset (Additional file 1: Table S6).
Differences in assignment percentages between SSR and SNP markers were also reported by independent studies in various maize germplasms [38, 48, 49] and were attributed to the greater information content of SSRs . According to Laval et al.  (k-1) times more biallelic markers should theoretically achieve the same genetic distance accuracy as an SSR set of k alleles. In our study an average of about 22 alleles per SSR marker was found, thus [(22–1)*22] = 462 SNP markers should be needed to resolve all the relationships that have been detected using SSR markers; that would mean 109 SNPs in addition to the 353 finally used. Increasing the number of SNP markers will probably improve the inference of population structure, even outperforming the results obtained with SSRs, since using the SNP dataset standard deviations of L(K) were smaller even at high values of K. Nevertheless, the models of either 2 main groups (ΔK method) or 5–6 subpopulations (Wilcoxon test) could be supported by both marker datasets while the distinction between Sativa and Sylvestris was better resolved by using the SNP markers. To understand how comparable are the structure outputs based on SSRs and SNPs, the level of membership correlation was investigated assuming five populations. This was because both inter-specific (S5 and S4, Sylvestris and Rootstocks respectively) and intra-specific (Sativa: S1, S2, S3) subdivisions were detected at K = 5 and, as stated above, no individual was strongly assigned to the additional inferred populations for K > 5. The relationship between membership in the S1, S2, S3, S4 and S5 populations based on SSRs and membership based on SNPs were plotted. Correlations were strong for S5 and S4 (R2 = 0.93 and 0.86, respectively) and moderate for S1, S2, S3 (R2 = 0.74, 0.63 and 0.57 respectively) but there was clearly much more spread along the x-axes (SNPs) than along the y-axes (SSRs). Finally, the classification in five groups was in agreement with NJ analysis since individuals assigned to the same genetic group tended to be close together in the NJ trees (Additional file 5).
Hierarchical population structure
Using the SSR set, a subsequent round (second round) of STRUCTURE separated most of the Sativa accessions from the group of Sylvestris. Following runs of STRUCTURE revealed a further stratification in the Sativa cluster into 3 groups (VV1, VV2 and VV3) and distinguished an additional group of Sativa accessions (VV4) from the Sylvestris genotypes (VS). Clusters VV1, VV2, VV3 and VV4 represent mainly Italian/Balkan wine grapes, Mediterranean table/wine grapes, Muscats (wine/table grapes) and Central European wine grapes, respectively. With the SNP dataset, the second round of STRUCTURE separated well the cultivated from the wild grapevines (VSI). At this point the threshold imposed on the SNP set excluded some Sativa accessions from further clustering (i.e. Pinot Noir, Gewürztraminer, Rhein Riesling, Aromriesling, Sauvignon, Sauvignonasse, Perle and Sacy among the true-to-type individuals). These accessions showed an admixed ancestry ≈ 0.5 of both Sativa and Sylvestris clusters and were grouped in VV4 when the SSR set was used. Following runs of STRUCTURE first separated the Sativa samples in two main groups and those were subsequently subdivided in five clusters (VV1I, VV2aI, VV2bI, VV3I, VV4I). Further clustering did not reveal anything new about the relationships amongst the accessions either because no individuals were strongly assigned (Q > 0.8 for SSRs and Q > 0.65 for SNPs) or a lack of background information on the samples did not allow the divisions to be independently supported. The clusters identified in the hierarchical and canonical STRUCTURE approach were in agreement, despite the different number of individuals assigned to the final populations (623 vs 831 for the SSR set and 443 vs 649 for the SNP set in the hierarchical and canonical STRUCTURE approach respectively). Moreover the hierarchical approach revealed the presence of an additional subgroup (VV2) not detected at K = 5. Thus, VS (VSI), Rs and RsI), VV3 (VV3I), and VV4 (VV4I) correspond to S4, S5, S2 and S3 respectively while VV1 (VV1I) and VV2 (VV2aI, VV2bI) were not separated at K = 5 and were grouped together in S1. Altogether, almost 97% of the individuals clustered in the hierarchical analysis were similarly clustered in the canonical approach at K = 5.
Clusters composition and the origin of cultivated grapevine subpopulations
Clusters detected through the hierarchical approaches were consistent for both marker sets, despite the difference between the group VV2 (SSR) which was split into VV2aI and VV2bI when the SNP set was used. The considerable degree of uncertainty about the variety names of many samples in the collection limits the interpretation of this kind of result. Nevertheless, when only true-to-type accessions are considered it can be seen that SNPs made a distinction between Mediterranean grapevines of table grape cultivars related to ‘Sultanina’ (i.e. ‘Calmeria’, ‘Flame Seedless’) and Spanish wine grape cultivars (i.e. ‘Xarello’, ‘Macabeu’, ‘Parellada’, ‘Beba’). All the clusters identified using the SSR set contained more accessions compared to those defined based on the SNP set, with the greatest differences observed for groups VV1 (129) vs. VV1I (51) and VV4 (112) vs. VV4I (52) (Additional file 1: Table S7). This can be explained again by differences in allelic richness of these marker types, with lower assignment success for SNPs, even when the proportional membership threshold was reduced to 0.65. However, almost half of the accessions grouped in each Sativa cluster using the SNP set were also grouped accordingly when using the SSR set. Additionally, none of the unshared accessions could be found grouped in a different cluster, supporting the robustness of the clustering method.
Relatedness among samples influences the ability of STRUCTURE to correctly detect the genetic stratification of a germplasm [52, 53]. This problem was partially overcome in the present study by investigating different criteria to find the best K in the “standard” STRUCTURE method, and then by applying a hierarchical approach. When using the former method most of the accessions sharing high order of pedigree relationships were not grouped into a specific cluster and were ultimately excluded. Thus, the final clusters contain mainly first degree pedigree relationships and represent the most plausible genetic structure of the germplasm investigated, being the smallest number of populations (K) “that captures the major structure of the data” . STRUCTURE detected additional subpopulations within these groups as a consequence of sample relatedness. For instance, VV3 is further divided into offspring either of ‘Muscat of Alexandrie’ or ‘Moscato Bianco’ that are considered two of the oldest grape varieties still in existence. Moreover, this clustering substantially agrees with the classification of eco-geographic variation proposed by Negrul  and Levadoux  as well as with previous genetic structure analysis performed on a restricted number of cultivated (222) and wild (22) grapevines from a different germplasm . According to Negrul , Italian and Greek wine grapes (VV1) belong to the group pontica and the French and German wine grapes (VV4) belong to the group occidentalis, whereas the Muscat table and wine cultivars (VV3) belong to the group orientalis (sub-proles caspica). The composition of the group VV2 is more heterogeneous since it includes table grape varieties related to ‘Sultanina’ (group orientalis sub-proles antasiatica) and some Spanish wine grapes, whose origin is still unknown. Aradhya et al.  reported similar results of grouping seedless type cultivars along with a number of southern European minor wine varieties into the “Mediterranean table-grape” cluster.
Genetic diversity among clusters
Likewise, the pairwise FST analysis pointed out that these common clusters defined by hierarchical STRUCTURE represent statistically supported subpopulations (Additional file 1: Table S8). Very similar FST values were also found when considering all the accessions, grouped using the SNP and SSR marker sets separately (data not shown). As expected, the highest pairwise genetic differentiation was observed among clusters of Sylvestris and Rootstocks. In the case of the SNP set, the FST values are much higher compared to the SSR set when estimated for relationship with Rootstocks and Sylvestris. This could be due to bias introduced because the panel of SNPs was designed from a small sample of accessions, preventing detection of multiallele polymorphisms or because additional SNPs were located in the vicinity.
Based on SSRs, the group of Muscats (VV3) and the Mediterranean table and wine grapes (VV2) showed slightly higher FST values when compared to Sylvestris (FST = 0.18 and 0.17, respectively) than to accessions of Rootstocks (FST = 0.16 and 0.15, respectively). On the contrary the Italian-Greek wine grapes (VV1) and the French-German wine grapes (VV4) showed a lower level of differentiation (FST = 0.14 and 0.13, respectively) than the cluster of Sylvestris. This result might suggest a moderate genetic exchange between the wild grapevines and only a part of the cultivated grapevines, namely some of those belonging to the groups of pontica and occidentalis.
Despite a very low genetic differentiation (FST = 0.0048) among cultivars of Spanish and France-Central Europe origin being recently reported using SSRs , a moderate genetic differentiation (FST ranging from 0.09 to 0.15) among the cultivated grapevine clusters was detected based on SSRs in the present study. A moderate genetic differentiation was also found by using SNPs (FST ranging from 0.07 to 0.12) which is slightly higher compared to the results obtained among Western, Central and Eastern European cultivars (FST ranging from 0.02 to 0.051; ). Since the low genetic differentiation among cultivars reported up to now has been suggested to be a consequence of their complex pedigree , the higher values of FST estimated in the present study support the consistency of the intra-specific clusters detected by the hierarchical STRUCTURE approach.
Genetic differentiation among clusters of grapevines
Archeological and historical evidences suggest that grape domestication took place in the Near East  and several studies have proposed the existence of secondary domestication events also along the Mediterranean basin [1, 5, 14, 15]. The findings of Myles et al.  support an Eastern origin of most grape cultivars as well as the existence of introgression from wild germplasm in Western region. This was also suggested by de Andrès et al.  based on the genetic relationships between wild and cultivated Spanish grapevines and agrees with Negrul  and Levadoux  who suggest the cultivars of the occidentalis group possess wild morphological characters as evidence of spontaneous introgression from V. sylvestris into cultivated forms of the pontica group. Accordingly, our results put the occidentalis group (VV4), consisting of wine grapes mainly related to Pinot Noir and Traminer, closer to the wild samples compared to pontica (Italian and Greek) wine grapes. The origins of Traminer and Pinot Noir are poorly known and presumably ancient . It has been suggested they could have arisen from hybridization between Roman grapes and local wild populations or from secondary domestication of the latter . The ancient origin of Pinot Noir and Traminer may also be deduced from the evidence that many modern varieties are first-degree relatives of these cultivars [57, 59].
The hierarchical STRUCTURE analysis grouped wild grapevines into a genetically distinct cluster, which, however, included some Sativa accessions (4 using the SSR set and 6 using the SNP set), while no additional subdivision of the cluster was detected. True-to-type Italian cultivars “Lambrusco a foglia frastagliata” (hermaphrodite) and “Lambrusco di Sorbara” (female) were common to both datasets. This apparently agrees with Levadoux et al.  who identified certain forms within the wild "Lambrusque" grapevines displaying traits commonly associated with domestication.
In this study we provided standard marker profiles for the trueness-to-type assessment of grape cultivars and marker-assisted reduction of redundancy in a grape germplasm collection. Using 22 SSR and 384 SNP loci, we showed that both the marker systems are efficient for the evaluation of genetic diversity and population structure, and microsatellites turned out to be well suited for the construction of core collections. This is an important step to sustainable and effective use of available grape genetic resources in basic and applied research. For instance, the core collections may contribute to development of association mapping populations for investigating genotype-phenotype relationships. Our complementary approaches to the analysis of SSR and SNP datasets detected consistent inter- and intraspecific levels of germplasm stratification with four ancestral subpopulations of V. vinifera ssp. sativa. This is in accord with the eco-geographical origin of the cultivars.
Plant material and DNA extraction
A total of 2273 accessions of grape (Vitis spp.) were analyzed in this study. They belong to the FEM grape germplasm collection (ITA362), located in San Michele all'Adige, Italy (46°18’ N, 11°13’ E). All plants were grafted on the rootstock Kober 5BB in five replicates and trained according to the Guyot system. Each accession was classified to a certain category based on the collections record; thus the material consisted of 1659 cultivated grapevines, V. vinifera ssp. sativa (Sativa), 177 wild individuals of V. vinifera ssp. sylvestris (Sylvestris), 127 interspecific hybrids used for fruit production (Hybrids) and 310 accessions of rootstock varieties including wild non-vinifera Vitis species (Rootstocks). Based on phenological data recorded for every Sativa accession through visual inspection in the field, a set of 163 grapevine cultivars was previously sorted from the germplasm collection. This set can be considered representative with regard to the diverse timing of budburst, flowering, véraison and full ripening, and will be referred in this study as the “phenological core” (core P). The range in dates of the major growth stages observed in the core P since 2008 using the modified E-L system  is reported in Additional file 1: Table S2.
Detailed information about each accession is publicly available at the European Vitis Database .
Young leaf tissue of one field grown plant per accession was harvested and stored immediately in 96-well microtube plates. Two controls (Pinot Noir and Sangiovese cultivars) were added to each set of 94 accessions for DNA extraction and successive analyses. Total genomic DNA was isolated from freeze-dried tissue after grinding with the MM 300 Mixer Mill system (Retsch., Germany). DNA extraction was performed using the DNeasy 96 plant mini kit (QIAGEN, Germany). DNA was suspended in TE buffer (pH = 8) and digested with RNase A (QIAGEN) at 37°C for 1 h. Next the DNA samples were diluted to approximately 4 ng/μl before conducting PCR.
SSR selection and genotyping
Twenty two SSR markers previously developed for grape, scattered over the genome, were chosen in this study (Additional file 1: Table S3): twenty SSR markers used by Laucou et al.  and the markers VrZAG62 and VrZAG79 . This set includes the 9 SSR markers proposed by the European Project GrapeGen06 for the characterization of regional cultivars .
Nine multiplex panels of fluorescent-labeled microsatellite loci were used. Simultaneous PCR amplifications were carried out in a final volume of 12.5 μl containing 10 ng of genomic DNA, 0.25 mM of each dNTPs, 2 mM MgCl2, 1.5 U Taq DNA Polymerase (AmpliTaq, Gold™, Applied Biosystems, Foster City, CA). Depending on the locus, primer concentrations ranged from 0.2 to 0.6 μM. Reactions were performed on a GeneAmp PCR System 9700 using the following profile: a hot start of 95°C for 7 min, 30 amplification cycles of 45 sec at 95°C, 1 min at 54°C, 30 sec at 72°C, and a final extension step of 1 hour at 72°C. PCR products (0.5 μl) generated by two or three different fluorescence dye-labeled primers were mixed with 9.3 μl of formamide and 0.2 μl of the GeneScan™ 500 ROX® Size Standard (Applied Biosystems). The DNA fragments were denatured and size fractionated using capillary electrophoresis on an ABI 3130 Genetic Analyzer (Applied Biosystems). Subsequently, GeneMapper v3.5 (Applied Biosystems) was used for the allele size estimation. Rates of missing data (MD) were below 1% for markers included in the multiplex panels 1, 2, 3 and 4, while for the remaining loci MD were below 5% (for details see Additional file 1: Table S3).
SNP selection and genotyping
The set of 384 SNPs used in this study was selected from informative data produced by previous SNP discovery and validation projects [19, 39, 62, 63] and included 35 out of the 48 SNPs proposed by Cabezas et al.  for the identification of grapevine cultivars. The diverse discovery panels included 11 samples of ancient unrelated cultivars and wild V. vinifera for 164 SNPs found by Lijavetzky et al. , 10 cultivated V. vinifera and 7 wild Vitis species for 88 SNPs identified by Myles et al.  and 10 complex interspecific grape hybrids for 9 SNPs discovered by Zyprian et al. (in preparation). Of the remaining SNPs, 122 were identified in the cultivar Pinot Noir  and one was discovered in the Muscat family of grapevine . All these SNPs meet the criteria of having enough upstream and downstream sequence information and of absence of other known SNPs in their vicinity. A designability score calculated for each SNP by Illumina was higher than 0.6, and this predicted high assay conversion rates. Genetic map positions were known for 257 of these SNPs, with 20, 9, 11, 17, 10, 14, 18, 23, 12, 16, 8, 11, 12, 15, 10, 10, 7, 20 and 14 loci placed on linkage groups 1 to 19, respectively. Genomic locations of SNPs on the reference grape genome (inbred Pinot Noir, 8X) are shown in Additional file 6. Chromosomal location and allelic variants for each SNP marker are reported in Additional file 1: Table S4. The genotyping was performed on an Illumina BeadStation 500 G system at Parco Tecnologico Padano (Lodi, Italy), using the protocol supported by Illumina.
Genetic diversity within and among groups of the germplasm collection
The genotypic data were subjected to various within and among groups genetic diversity measures, such as mean number of alleles per locus (a), number of effective alleles (AE, the number of equally frequent alleles required to give the observed level of heterozygosity; ), levels of observed (HO) and expected (HE) heterozygosity , genetic differentiation (FST) and the fixation index (F, inbreeding coefficient; ). All these calculations, together with values of PI (the probability that two individuals in the population share the same genotype, ) and MAF (minor allele frequency), were performed using GenAlex 6.41  and GENETIX .
Construction of genetic core collections
Core collections are subsamples of larger germplasm collections and include accessions chosen to represent the majority of the genetic diversity contained in these larger collections . To construct the genetic core collections we used the M (maximization) method, suggested by Schoen and Brown , implemented in the MSTRAT software . The M strategy selects specific combinations of accessions while maximizing the number of observed alleles at each marker locus and the MSTRAT uses iterative procedures to select samples with the highest allelic diversity. The final number of iterations per MSTRAT run was 150, while the number of repetitions for core sampling was 100. Putative core collections exhibiting the same allelic richness were ranked using Nei’s diversity index . The accessions that were most often present in the 100 replicates were retained as the final core collection. The efficiency of the strategy was assessed by comparing the total number of alleles captured using MSTRAT in samples of increasing size to the number of alleles captured in randomly chosen collections of equal size (ten independent samplings).
Analysis of population structure
The genetic structure of the germplasm collection was analyzed performing Principal Coordinate Analysis (PCoA) and Principal Component Analysis (PCA) implemented in the programs GenAlex 6.41  and GenoDive 2.0b21 , and by using STRUCTURE 2.1 software [47, 52]. PCoA was based on standardized covariance of genetic distances calculated for codominant markers while PCA was calculated by using the variance-covariance matrix of allele frequencies. Missing data were replaced by alleles randomly picked from the allele pool of each population. To avoid bias in the analysis of population structure, pairwise linkage disequilibrium (LD) between SNPs was evaluated using the software TASSEL v2.1  by setting 1000 permutations. The model-based approach implemented in STRUCTURE 2.1 was used at the Bioportal server . This software applies a Bayesian clustering algorithm to identify subpopulations, assign individuals to them, and estimate the population allele frequencies. STRUCTURE sorts individuals into K clusters, according to their genetic similarity. The best K is chosen based on the estimated membership coefficients (Q) for each individual in each cluster. Ten independent runs for K values ranging from 1 to 20 were performed with a burn-in length of 1000000 followed by 1500000 iterations. The admixture model was applied and no prior population information was used. The log-probability of the data, given for each value of K, was calculated and compared across the range of K. The software CLUMPP 1.1  was used to find optimal alignments of independent runs and the output was used directly as input into a program for cluster visualization DISTRUCT 1.1 . The optimal subpopulation model was investigated in several ways: (1) by applying the informal pointers (i.e. the plateau criterion) proposed by Pritchard et al. , (2) by evaluating L(K), the log probability values from ten runs at each K, using non-parametric Wilcoxon test, as implemented in the R software package Rcmdr , (3) by considering 3ΔK, a second order rate change with respect to K, as defined in Evanno et al. , and (4) by plotting the log probability L(K) and ΔK of the data over ten runs, as implemented in STRUCTURE HARVESTER . In addition the following parameters were calculated using the R-script Structure-Sum : the average similarity coefficients for different simulations within each preset K  and the extent of membership in a single cluster measured by the clusteredness statistic . Furthermore, a “hierarchical STRUCTURE analysis”  was applied in this study by running STRUCTURE subsequently on partitioned data, as suggested by Pritchard et al. , using only the individuals suspected of being subdivided into different clusters. For this approach the ΔK method  was used in adjudication for the best K and the individuals with a proportional membership Q > 0.8 (SSRs) and Q > 0.65 (SNPs) in their primary population were considered in the subsequent analysis.
The consistency of the clusters identified through the hierarchical STRUCTURE approach was tested by pairwise FST analysis . In addition, an unweighted neighbour-joining tree was constructed based on dissimilarities between 330 accessions (calculated from 22 SSRs), and ten thousand bootstrap replicates were performed using the Darwin software package v5.0148 .
We are grateful to Maddalena Sordo for technical assistance and to Silvano Clementi, Monica Dallaserra and Alessandra Zatelli for contributing in collection of plant materials, and Diego Lijavetzky for processing part of the SNP sequences. Paul K. Boss and Paula Moreno-Sanz are kindly acknowledged for critical reading of the manuscript. This study was carried out in the framework of the European research project GrapeGen06 "Management & Conservation of Grapevine Genetic Resources" with the financial support of Fondazione Edmund Mach, DGBPF-GAV-AdP 2010–2012. JMMZ was funded by EUI2008-03752 grant from MICINN (Spain).
Joint publication of the COST Action FA1003 “East–West Collaboration for Grapevine Diversity Exploration and Mobilization of Adaptive Traits for Breeding”.
- This P, Lacombe T, Thomas MR: Historical origins and genetic diversity of wine grapes. Trends Genet. 2006, 22: 511-519. 10.1016/j.tig.2006.07.008.PubMedView ArticleGoogle Scholar
- Olmo HP: Grapes. Evolution of Crop Plants. Edited by: Simmonds NW. London: Longman; 1976: 294-298.Google Scholar
- Cervera MT, Cabezas JA, Sancha JC, de Martínez TF, Martínez-Zapater JM: Application of AFLPs to the characterization of grapevine Vitis vinifera L. genetic resources. A case study with accessions from Rioja (Spain). Theor Appl Genet. 1998, 97: 51-59. 10.1007/s001220050866.View ArticleGoogle Scholar
- Sefc KM, Lopes MS, Lefort F, Botta R, Roubelakis-Angelakis KA, Ibáñez J, Pejic I, Wagner HW, Glössl J, Steinkellner H: Microsatellite variability in grapevine cultivars from different European regions and evaluation of assignment testing to assess the geographic origin of cultivars. Theor Appl Genet. 2000, 100: 498-505. 10.1007/s001220050065.View ArticleGoogle Scholar
- Arroyo-García R, Ruiz-Garcia L, Bolling L, Ocete R, Lopez MA, Arnold C, Ergul A, Soylemezoglu G, Uzun HI, Cabello F, Ibáñez J, Aradhya MK, Atanassov A, Atanassov I, Balint S, Cenis JL, Costantini L, Goris-Lavets S, Grando MS, Klein BY, McGovern PE, Merdinoglu D, Pejic I, Pelsy F, Primikirios N, Risovannaya V, Roubelakis-Angelakis KA, Snoussi H, Sotiri P, Tamhankar S: Multiple origins of cultivated grapevine (Vitis vinifera L. ssp sativa) based on chloroplast DNA polymorphisms. Mol Ecol. 2006, 15 (12): 3707-3714. 10.1111/j.1365-294X.2006.03049.x.PubMedView ArticleGoogle Scholar
- This P, Jung A, Boccacci P, Borrego J, Botta R, Costantini L, Crespan M, Dangl GS, Eisenheld C, Ferreira-Monteiro F, Grando S, Ibáñez J, Lacombe T, Laucou V, Magalhaes R, Meredith CP, Milani N, Peterlunger E, Regner F, Zulini L, Maul E: Development of a standard set of microsatellite reference alleles for identification of grape cultivars. Theor Appl Genet. 2004, 109 (7): 1448-1458. 10.1007/s00122-004-1760-3.PubMedView ArticleGoogle Scholar
- Ibanez J, Velez MD, de Andrés MT, Borrego J: Molecular markers for establishing distinctness in vegetatively propagated crops: a case study in grapevine. Theor Appl Genet. 2009, 119: 1213-1222. 10.1007/s00122-009-1122-2.PubMedView ArticleGoogle Scholar
- Cipriani G, Spadotto A, Jurman I, Di Gaspero G, Crespan M, Meneghetti S, Frare E, Vignani R, Cresti M, Morgante M, Pezzotti M, Pe E, Policriti A, Testolin R: The SSR-based molecular profile of 1005 grapevine (Vitis vinifera L.) accessions uncovers new synonymy and parentages, and reveals a large admixture amongst varieties of different geographic origin. Theor Appl Genet. 2010, 121: 1569-1585. 10.1007/s00122-010-1411-9.PubMedView ArticleGoogle Scholar
- Laucou V, Lacombe T, Dechesne F, Siret R, Bruno JP, Dessup M, Dessup T, Ortigosa P, Parra P, Roux C, Santoni S, Varès D, Péros JP, Boursiquot JM, This P: High throughput analysis of grape genetic diversity as a tool for germplasm collection management. Theor Appl Genet. 2011, 122: 1233-1245. 10.1007/s00122-010-1527-y.PubMedView ArticleGoogle Scholar
- Myles S, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK, Prins B, Reynolds A, Chia J-M, Ware D, Bustamante CD, Buckler ES: Genetic structure and domestication history of the grape. PNAS. 2011, 108 (9): 3457-3458. 10.1073/pnas.1019729108.View ArticleGoogle Scholar
- Maul E, Sudharma KN, Kecke S, Marx G, Müller C, Audeguin L, Boselli M, Boursiquot JM, Bucchetti B, Cabello F, Carraro F, Crespan M, De Andrés MT, Dias JE, Ekhvaia J, Gaforio L, Gardiman M, Grando MS, Gyropoulos D, Jandurova O, Kiss E, Kontic J, Kozma P, Lacombe T, Laucou V, Legrand D, Maghradze D, Marinoni D, Maletic E, Moreira F: The European Vitis Database (http://www.eu-vitis.de): a technical innovation through an online uploading and interactive modification system. Vitis. 2012, 51 (2): 79-85.Google Scholar
- Cabezas JA, Ibáñez J, Lijavetzky D, Vélez D, Bravo G, Rodríguez V, Carreño I, Jermakow AM, Carreño J, Ruiz-García L, Thomas MR, Martinez-Zapater JM: A 48 SNP set for grapevine cultivar identification. BMC Plant Biol. 2011, 11: 153. 10.1186/1471-2229-11-153.PubMedPubMed CentralView ArticleGoogle Scholar
- Pelsy F, Hocquigny S, Moncada X, Barbeau G, Forget D, Hinrichsen P, Merdinoglu D: An extensive study of the genetic diversity within seven French wine grape variety collections. Theor Appl Genet. 2010, 120 (6): 1219-1231. 10.1007/s00122-009-1250-8.PubMedView ArticleGoogle Scholar
- Aradhya MK, Dangl GS, Prins BH, Boursiquot JM, Walker MA, Meredith CP, Simon CJ: Genetic structure and differentiation in cultivated grape, Vitis vinifera L. Genet Res Camb. 2003, 81: 179-182. 10.1017/S0016672303006177.View ArticleGoogle Scholar
- Grassi F, Labra M, Imazio S, Spada A, Sgorbati S, Scienza A, Sala F: Evidence of a secondary grapevine domestication centre detected by SSR analysis. Theor Appl Genet. 2003, 107: 1315-1320. 10.1007/s00122-003-1321-1.PubMedView ArticleGoogle Scholar
- Pelsy F: Untranslated leader region polymorphism of Tvv1, a retrotransposon family, is a novel marker useful for analyzing genetic diversity and relatedness in the genus Vitis. Theor Appl Genet. 2007, 116: 15-27. 10.1007/s00122-007-0643-9.PubMedView ArticleGoogle Scholar
- Santana JC, Heuertz M, Arranz C, Rubio JA, Martinez-Zapater JM, Hidalgo H: Genetic structure, origins, and relationships of grapevine cultivars from the Castilian plateau of Spain. Am J Enol Vitic. 2010, 61: 214-224.Google Scholar
- De Andrés MT, Cabezas JA, Cerveza MT, Borrego J, Martinez-Zapater JM, Jouve N: Molecular characterization of grapevine rootstocks maintained in germplasm collections. Am J Enol Vitic. 2007, 58: 75-86.Google Scholar
- Lijavetzky D, Cabezas JA, Ibáñez A, Rodriguez V, Martínez-Zapater JM: High throughput SNP discovery and genotyping in grapevine (Vitis vinifera L.) by combining a re-sequencing approach and SNPlex technology. BMC Genomics. 2007, 8: 424. 10.1186/1471-2164-8-424.PubMedPubMed CentralView ArticleGoogle Scholar
- Jones ES, Sullivan H, Bhattramakki D, Smith JSC: A comparison of simple sequence repeat and single nucleotide polymorphism marker technologies for the genotypic analysis of maize (Zea mays L.). Theor Appl Genet. 2007, 115: 361-371. 10.1007/s00122-007-0570-9.PubMedView ArticleGoogle Scholar
- Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, Lepoittevin C, Malausa T, Revardel E, Salin F, Petit RJ: Current trends in microsatellite genotyping. Mol Ecol Resour. 2011, 11: 591-611. 10.1111/j.1755-0998.2011.03014.x.PubMedView ArticleGoogle Scholar
- Nielsen R: Correcting for ascertainment biases when analyzing SNP data: applications to the estimation of linkage disequilibrium. Theor Popul Biol. 2003, 63: 245-255. 10.1016/S0040-5809(03)00005-4.PubMedView ArticleGoogle Scholar
- Moragues M, Comadran R, Waugh R, Milne I, Flavell AJ, Russell JR: Effects of ascertainment bias and marker number on estimations of barley diversity from high-throughput SNP genotype data. Theor Appl Genet. 2010, 120: 1525-1534. 10.1007/s00122-010-1273-1.PubMedView ArticleGoogle Scholar
- Carlson CS, Smith JD, Stanaway IB, Rieder MJ, Nickerson DA: Direct detection of null alleles in SNP genotyping data. Hum Mol Genet. 2006, 15 (12): 1931-1937. 10.1093/hmg/ddl115.PubMedView ArticleGoogle Scholar
- Thomas MR, Cain P, Scott NS: DNA Typing of Grapevines a Universal Methodology and Database for Describing Cultivars and Evaluating Genetic Relatedness. Plant Mol Biol. 1994, 25: 939-949. 10.1007/BF00014668.PubMedView ArticleGoogle Scholar
- Lamboy WF, Alpha CG: Using simple sequence repeats (SSRs) for DNA fingerprinting germplasm accessions of grape (Vitis L) species. J Am Soc Hort Sci. 1998, 123: 182-188.Google Scholar
- Lin H, Walker MA: Identifying grape rootstocks with simple sequence repeat (SSR) DNA markers. Am J Enol Vitic. 1998, 49: 403-407.Google Scholar
- Sefc KM, Regner F, Glössl J, Steinkellner H: Genotyping of grapevine and rootstock cultivars using microsatellite markers. Vitis. 1998, 37: 15-20.Google Scholar
- Sefc KM, Regner F, Glössl J, Steinkellner H: Identification of microsatellite sequences in Vitis riparia and their applicability for genotyping of different Vitis species. Genome. 1999, 42: 367-373.PubMedView ArticleGoogle Scholar
- Di Gaspero G, Peterluner E, Testolin R, Edwards KJ, Cipriani G: Conservation of microsatellite loci within the genus Vitis. Theor Appl Genet. 2000, 101: 301-308. 10.1007/s001220051483.View ArticleGoogle Scholar
- Arnold C, Rossetto M, McNally J, Henry RJ: The application of SSRs characterized for grape (Vitis vinifera) to conservation studies in Vitaceae. Am J Bot. 2002, 89: 22-28. 10.3732/ajb.89.1.22.PubMedView ArticleGoogle Scholar
- Rossetto M, Mcnally J, Henry RJ: Evaluating the potential of SSR: flanking regions for examining taxonomic relationships in the Vitaceae. Theor Appl Genet. 2002, 104: 61-66. 10.1007/s001220200007.PubMedView ArticleGoogle Scholar
- Fernandez MP, Nunez Y, Ponz F, Hernaiz S, Gallego FJ, Ibanez J: Characterization of sequence polymorphisms from microsatellite flanking regions in Vitis spp. Mol Breed. 2008, 22: 455-465. 10.1007/s11032-008-9189-z.View ArticleGoogle Scholar
- Péros JP, Berger G, Portemont A, Boursiquot JM, Lacombe T: Genetic variation and biogeography of the disjunct Vitis subgenus (Vitaceae). J Biogeogr. 2010, 38 (3): 471-486.View ArticleGoogle Scholar
- Le Cunff L, Fournier-Level A, Laucou V, Vezzulli S, Lacombe T, Adam-Blondon AF, Boursiquot JM, This P: Construction of nested genetic core collections to optimize the exploitation of natural diversity in Vitis vinifera L. subsp. sativa. BMC Plant Biol. 2008, 8: 31. 10.1186/1471-2229-8-31.PubMedPubMed CentralView ArticleGoogle Scholar
- McKhann HI, Camilleri C, Bérard A, Bataillon T, David JL, Reboud X, Le Corre V, Caloustian C, Gut IG, Brunel D: Nested core collections maximizing genetic diversity in Arabidopsis thaliana. Plant J. 2004, 38: 193-202. 10.1111/j.1365-313X.2004.02034.x.PubMedView ArticleGoogle Scholar
- Ellwood SR, D’Souza NK, Kamphuis LG, Burgess TI, Nair RM, Oliver RP: SSR analysis of the Medicago truncatula SARDI core collection reveals substantial diversity and unusual genotype dispersal throughout the Mediterranean basin. Theor Appl Genet. 2006, 112 (5): 977-983. 10.1007/s00122-005-0202-1.PubMedView ArticleGoogle Scholar
- Hamblin MT, Warburton ML, Buckler ES: Empirical Comparison of Simple Sequence Repeats and Single Nucleotide Polymorphisms in Assessment of Maize Diversity and Relatedness. PLoS One. 2007, 2 (12): e1367. 10.1371/journal.pone.0001367.PubMedPubMed CentralView ArticleGoogle Scholar
- Emanuelli F, Battilana J, Costantini L, Le Cunff L, Boursiquot JM, This P, Grando MS: A candidate gene association study on muscat flavor in grapevine (Vitis vinifera L.). BMC Plant Biol. 2010, 10: 241. 10.1186/1471-2229-10-241.PubMedPubMed CentralView ArticleGoogle Scholar
- Fournier-Level A, Lacombe T, Le Cunff L, Boursiquot JM, This P: Evolution of the VvMybA gene family, the major determinant of berry colour in cultivated grapevine (Vitis vinifera L.). Heredity. 2010, 104: 351-362. 10.1038/hdy.2009.148.PubMedView ArticleGoogle Scholar
- de Andrés MT, Benito A, Pérez-Rivera G, Ocete R, Lopez MA, Gaforio L, Muñoz G, Cabello F, Martínez-Zapater JM, Arroyo-García R: Genetic diversity of wild grapevine populations in Spain and their genetic relationships with cultivated grapevines. Mol Ecol. 2012, 21: 800-816. 10.1111/j.1365-294X.2011.05395.x.PubMedView ArticleGoogle Scholar
- Pritchard JK, Falush D: Documentation for Structure Software: Version 2.2. Chicago, USA: University of Chicago; 2007.Google Scholar
- Rosemberg NA, Burke T, Elo K, Feldman MW, Freidlin PJ, Groenen MAM, Hillel J, Mäki-Tanila A, Tixier-Boichard M, Vignal A, Wimmers K, Weigend S: Empirical evaluation of genetic clustering methods using multilocus genotypes from 20 chicken breeds. Genetics. 2001, 159: 699-713.Google Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Nordborg M, Hu TT, Ishino Y, Jhaveri J, Toomajian C, Zheng H, Bakker E, Calabrese P, Gladstone J, Goyal R, Jakobsson M, Kim S, Morozov Y, Padhukasahasram B, Plagnol V, Rosenberg NA, Shah C, Wall JD, Wang J, Zhao K, Kalbfleisch T, Schulz V, Kreitman M, Bergelson J: The Pattern of Polymorphism in Arabidopsis thaliana. PLoS Biol. 2005, 3 (7): e196. 10.1371/journal.pbio.0030196.PubMedPubMed CentralView ArticleGoogle Scholar
- Lia VV, Poggio L, Confalonieri VA: Microsatellite variation in maize landraces from Northwestern Argentina: genetic diversity, population structure and racial affiliations. Theor Appl Genet. 2009, 119: 1053-1067. 10.1007/s00122-009-1108-0.PubMedView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly PJ: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Van Inghelandt D, Melchinger AE, Lebreton C, Stich B: Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers. Theor Appl Genet. 2010, 120 (7): 1289-1299. 10.1007/s00122-009-1256-2.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang X, Gao S, Xu S, Zhang Z, Prasanna BM, Li L, Li J, Yan J: Characterization of a global germplasm collection and its potential utilization for analysis of complex quantitative traits in maize. Mol Breed. 2011, 28 (4): 511-526. 10.1007/s11032-010-9500-7.View ArticleGoogle Scholar
- Lu Y, Yan J, Guimarães CT, Taba S, Hao Z, Gao S, Chen S, Li J, Zhang S, Vivek BS: Molecular characterization of global maize breeding germplasm based on genome-wide single nucleotide polymorphisms. Theor Appl Genet. 2009, 120 (1): 93-115. 10.1007/s00122-009-1162-7.PubMedView ArticleGoogle Scholar
- Laval G, San Cristobal M, Chevalet C: Measuring genetic distances between breeds: use of some distances in various short term evolution models. Genet Sel Evol. 2002, 34: 481-507. 10.1186/1297-9686-34-4-481.PubMedPubMed CentralView ArticleGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of Population Structure Using Multilocus Genotype Data: Linked Loci and Correlated Allele Frequencies. Genetics. 2003, 164: 1567-1587.PubMedPubMed CentralGoogle Scholar
- Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9: 1322-1332. 10.1111/j.1755-0998.2009.02591.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Negrul AM: Evolution of cultivated forms of grapes. CR Acad Sci USSR. 1938, 18: 585-588.Google Scholar
- Levadoux L: Les populations sauvages et cultivées de Vitis vinifera L. Annales de l’amélioration des plantes. 1956, 6: 59-118.Google Scholar
- McGovern PE, Glusker DL, Exner LJ, Voigt MM: Neolithic resinated wine. Nature. 1996, 381: 480-481. 10.1038/381480a0.View ArticleGoogle Scholar
- Regner F, Stadlbauer A, Eisenheld C, Kaserer H: Genetic Relationships Among Pinots and Related Cultivars. Am J Enol Vitic. 2000, 51 (1): 7-14.Google Scholar
- Bouquet A: Origine et évolution de l’encépagement français à travers les siècles. Progrès Agricole et Viticole. 1982, 5: 110-121.Google Scholar
- Boursiquot JM, Lacombe T, Bowers J, Meredith C: Le Gouais, un cépage clé du patrimoine viticole européen. Bulletin de l’OIV. 2004, 77: 799-809.Google Scholar
- Coombe BG: Adoption of a system for identifying grapevine growth stages. Aust J Grape Wine Res. 1995, 1: 100-110. 10.1111/j.1755-0238.1995.tb00085.x.View ArticleGoogle Scholar
- The European Vitis Database: [http://www.eu-vitis.de].
- Troggio M, Malacarne G, Coppola G, Segala C, Cartwright DA, Pindo M, Stefanini M, Mank R, Moroldo M, Morgante M, Grando MS, Velasco R: A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring ‘Pinot Noir’ bacterial artificial chromosome contigs. Genetics. 2007, 176: 2637-2650. 10.1534/genetics.106.067462.PubMedPubMed CentralView ArticleGoogle Scholar
- Myles S, Chia JM, Hurwitz B, Simon C, Zhong GY, Buckler E, Ware D: Rapid Genomic Characterization of the Genus Vitis. PLoS One. 2010, 5 (1): e8219. 10.1371/journal.pone.0008219.PubMedPubMed CentralView ArticleGoogle Scholar
- Kimura M, Crow JF: The number of alleles that can be maintained in a finite population. Genetics. 1964, 49: 725-738.PubMedPubMed CentralGoogle Scholar
- Nei M: Analysis of gene diversity in subdivided populations. PNAS. 1973, 70 (12): 3321-3323. 10.1073/pnas.70.12.3321.PubMedPubMed CentralView ArticleGoogle Scholar
- Wright S: The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution. 1965, 19: 395-420. 10.2307/2406450.View ArticleGoogle Scholar
- Paetkau D, Strobeck C: Microsatellite analysis of genetic variation in black bear populations. Mol Ecol. 1994, 3: 489-495. 10.1111/j.1365-294X.1994.tb00127.x.PubMedView ArticleGoogle Scholar
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Montpellier (France): Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II; 2004.Google Scholar
- Frankel OH: Genetic perspectives of germplasm conservation. Genetic manipulation: impact on man and society. Edited by: Arber W, Llimensee K, Peacock WJ, Starlinger P. Cambridge: Cambridge University Press; 1984: 161-170.Google Scholar
- Schoen DJ, Brown AH: Conservation of allelic richness in wild crop relatives is aided by assessment of genetic markers. PNAS. 1993, 90 (22): 10623-10627. 10.1073/pnas.90.22.10623.PubMedPubMed CentralView ArticleGoogle Scholar
- Gouesnard B, Bataillon TM, Decoux G, Rozale C, Schoen DJ, David JL: MSTRAT: An Algorithm for Building Germ Plasm Core Collections by Maximizing Allelic or Phenotypic Richness. J Hered. 2001, 92 (1): 93-94. 10.1093/jhered/92.1.93.PubMedView ArticleGoogle Scholar
- Nei L: Molecular evolutionary genetics. New York: Columbia University Press; 1987.Google Scholar
- Meirmans PG, Van Tienderen PH: GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. Mol Ecol Notes. 2004, 4: 792-794. 10.1111/j.1471-8286.2004.00770.x.View ArticleGoogle Scholar
- TASSEL v2.1: [http://www.maizegenetics.net/]
- Bioportal University of Oslo: [http://www.bioportal.uio.no/]
- Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23 (14): 1801-1806. 10.1093/bioinformatics/btm233.PubMedView ArticleGoogle Scholar
- Rosenberg NA: DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004, 4: 137-138.View ArticleGoogle Scholar
- Rosenberg NA, Burke T, Elo K, Feldman MW, Freidlin PJ, Groenen MAM, Hillel J, Mäki-Tanila A, Tixier-Boichard M, Vignal A, Wimmers K, Weigend S: Empirical evaluation of genetic clustering methods using multilocus genotypes from 20 chicken breeds. Genetics. 2001, 159: 699-713.PubMedPubMed CentralGoogle Scholar
- Earl DA, von Holdt BM: STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genet Resour. 2012, 4: 359-361. 10.1007/s12686-011-9548-7.View ArticleGoogle Scholar
- Ehrich D: AFLPDAT: a collection of r functions for convenient handling of AFLP data. Mol Ecol Notes. 2006, 6: 603-604. 10.1111/j.1471-8286.2006.01380.x.View ArticleGoogle Scholar
- Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, Feldman MW: Clines, Clusters, and the Effect of Study Design on the Inference of Human Population Structure. PLoS Genet. 2005, 1 (6): e70. 10.1371/journal.pgen.0010070.PubMedPubMed CentralView ArticleGoogle Scholar
- Vähä JP, Erkinaro J, NIemelä E, Primmer CR: Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. Mol Ecol. 2007, 16: 2638-2654. 10.1111/j.1365-294X.2007.03329.x.PubMedView ArticleGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Perrier X, Jacquemoud-Collet JP: DARwin software. 2006. http://darwin.cirad.fr/.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.