Genetic variability assessment in ‘ Muscat ’ grapevines including ‘ Muscat of Alexandria ’ clones from selection programs

Genetic variability is needed to face environmental changes and pathogen constrains. In addition, the search for intravarietal variability contributes to the avoidance of genetic erosion, preserving clones that are adapted to particular conditions. Variability is also important to diversify grapevine-derived products. In this work, we have analyzed the genetic variability of ‘Muscat germplasm’ including samples from neglected vineyards from Alicante and Valencia provinces, accessions of the germplasm collections of ‘Colección de Vides de El Encín’ (Alcalá de Henares, Madrid) and ‘La Casa de las Vides’ (Agullent, Valencia), accessions supplied by nurseries of Valencia province, and ‘Muscat of Alexandria’ clones selected using differential ampelographic characteristics in selection programs (La Marina, Alicante). Fifteen microsatellites (SSRs) were used to study intervarietal variability. The SSR fingerprinting allowed the identification of some accessions, variants, and synonymies. Amplified Fragment Length Polymorphisms (AFLPs) markers and Microsatellite-AFLPs were used to determine the variability attended in ‘Muscat of Alexandria’ accessions. A CAPs (Cleaved Amplified Polymorphic Sequences) marker, recently developed for the discrimination of ‘Muscat’ flavor genotypes using the SNP1822 G>T, was assessed and showed that all the analyzed accessions were ‘Muscat’ flavored. The variation found among the analyzed germplasm is very interesting because variants within ‘Muscat of Alexandria’, ‘Muscat Italia’, and ‘Muscat d’Istambul’ have been identified. In addition, intravarietal genetic variation was found among the analyzed accessions in ‘Muscat of Alexandria’ from selection programs. Additional keywords: microsatellites; AFLPs; M-AFLPs; ‘Moscatel’, Vitis vinifera. Abbreviations used: AFLPs (Amplified Fragment Length Polymorphisms); CAPS (Cleaved Amplified Polymorphic Sequences); PDO (Protected Denomination of Origin); He (Expected heterozygosity); Ho (Observed heterozygosity); M-AFLPs (MicrosatelliteAmplified Fragment Length Polymorphisms); PCR (Polymerase Chain Reaction); PIC (Polymorphic Information Content); SSR (Simple Sequence Repeat or Microsatellites); UPGMA (Unweighted Pair Group Method with Arithmetic mean); VIVC (Vitis International Variety Catalogue). Authors’ contributions: Conceived the work: JXS and CG. Designed and performed the experiments: RP, AC and CG. Analyzed the data: RP and CG. Contributed with plant materials: FC, JXS and CJ. Wrote the paper: CG. Contributed with critical reading of the manuscript: RP, JXS, CJ, and FC. All authors read and approved the final manuscript. Supplementary material: (Table S1 and Fig. S1) accompanies the paper on SJAR’s website. Citation: Peiró, R.; Soler, J. X.; Crespo, A.; Jiménez, C.; Cabello, F.; Gisbert, C. (2018). Genetic variability assessment in ‘Muscat’ grapevines including ‘Muscat of Alexandria’ clones from selection programs. Spanish Journal of Agricultural Research, Volume 16, Issue 2, e0702. https://doi.org/10.5424/sjar/2018162-12537 Received: 06 Nov 2017. Accepted: 22 May 2018. Copyright © 2018 INIA. This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC-by 4.0) License. Funding: MINECO and FEDER Funds (project CGL2015-70843-R); Conselleria d’Agricultura, Pesca i Alimentació, Generalitat Valenciana, Spain (project 2007TAHALI00073). Competing interests: The authors have declared that no competing interests exist. Correspondence should be addressed to Carmina Gisbert: cgisbert@btc.upv.es


Introduction
Since the domestication of grapevine (Vitis vinifera L.) between the seventh and the fourth millennia BC, in a geographical area between the Black Sea and Iran (McGovern & Rudolph, 1996;McGovern et al., 1996;Zohary & Hopf, 2000), cultivated forms have been spread by humans in the Near East, Middle East, and Central Europe.From these areas, that are considered as secondary domestication centers (Grassi et al., 2003;Arroyo-García et al., 2006), the culture of grapevine varieties was spread to warm areas worldwide.The accumulation of casual mutations and natural or artificial crossing with other cultivars and wild grapevines have been the drivers of grapevine evolution since its domestication (This et al., 2006;Forni, 2012).Nowadays, grapevine is one of the economically most important crops cultivated around the world.
Among grapevine varieties, the 'Muscat' family includes a wide spread of grapevines having in common a pronounced floral aroma and a typical 'Muscat' flavor.A high concentration of monoterpenes in the grapes is related to the pronounced sweet floral aroma (Mateo & Jimenez, 2000).The most representative and ancient varieties are 'Muscat à Petit Grains blanc' and 'Muscat of Alexandria' (Robinson, 1986).It was confirmed by DNA analysis that 'Muscat of Alexandria' was itself a natural crossing of 'Muscat à Petit Grains blanc' and a white-skinned table grape variety from Greek islands known as 'Axina de Tres Bias' or 'Heptakilo', which is rarely seen outside of Greece, Malta, or Sardinia (Cipriani et al., 2010).'Muscat' varieties have been known and appreciated since ancient times (Scienza et al., 1989) and could be among those named by Plinius (in the 1 st century AD) as 'apianas'; varieties greatly appreciated by bees (Apis) (Columela, cited by Álvarez de Sotomayor, 1979).The exact origins of these varieties are not clear and different theories exist regarding their origins as well as the name 'Muscat' (Dalmasso et al., 1964;Robinson et al., 2012).What is true is that 'Muscat à Petits Grains blanc' and 'Muscat of Alexandria' have been crossed between themselves and with other known or unknown varieties to produce the great family of 'Muscats' (Robinson et al., 2012).In the Vitis International Variety Catalogue (VIVC) database (www.vivc.de),we can find that 'Muscat Hamburg', which resulted from the cross of 'Muscat of Alexandria' and 'Frankenthal' (synonym: 'Schiava Grossa') (Crespan, 2003;Ibáñez et al., 2009;Lacombe et al., 2013), has also actively contributed to the development of new Muscat hybrids.
Grapes from 'Muscat of Alexandria' and 'Muscat à Petit Grains blanc' are consumed as table grapes or for raisin production, as well as being the basis of different appreciated wines.For instance, 'Muscat à Petits Grains blanc', also known as 'Moscato Asti' in Italy, is the main grape variety used in the production of the Italian sparkling wine Asti or the Liqueur Muscat in Australia.'Muscat of Alexandria' is used in the production of the fortified Spanish Moscatel and different Moscato wines, which are also produced in different countries like Australia, Brazil or South Africa.Both varieties are also commonly used in the elaboration of the French vins doux naturels (Robinson, 1986).Robinson et al. (2012) reported 49 and 66 synonymies for 'Muscat of Alexandria' and 'Muscat à Petit Grains blanc', respectively.In the VIVC database, 220 synonyms were found for the first.The high quantity of nouns reflects their great importance and their wide distribution around the world.Even in a single country, several names were found for a particular cultivar.For instance, for 'Muscat of Alexandria', 'Argelino', or 'Muscat de Raf-Raf' are employed in Tunisia; 'Moscato Alexandrias', 'Moscato de Limnou', or 'Anglico' are used in Greece, and 'Muscat à Gros Grains', 'Muscat de Rome', or 'Muscat Grec' are used in France.It is also common to find varieties or products erroneously associated with 'Muscat' grapes, a consequence of naming them for their aromatic character.For instance, several wine grape varieties − such as 'Chardonnay', 'Chasselas', or 'Sauvignon blanc' − are often suffixed with Musqué.The species Vitis rotundifolia is also commonly known as a 'Muscadine' grape, which may induce confusion (Robinson et al., 2012).Recently, Emanuelli et al. (2014) developed a molecular marker related to the VvDXS gene, which confers 'Muscat flavor' in grapevine.This marker can help in the characterization of 'Muscat' and related varieties.
Along the Mediterranean coast of Spain, viticulture has been present from ancient times.At the archeological site of 'L'Alt de Benimaquia' (seventh century BC), an Iberian village located in the mountain 'El Montgó' in Denia (a village of the La Marina Alta, Alicante), significant quantities of vinification residues (high amounts of seeds and tartaric acid) were found (Dies et al., 1993;Hidalgo, 1999).Even though it is not known which variety(ies) correspond(s) to the seeds found at this archeological site, their age indicates the importance of V. vinifera in this Mediterranean area since ancient times.Evidence for the culture of a high number of varieties in La Marina was reported by Cavanilles (1795), Chabás (1972), and Calvo (2003).However, the first report of the name 'Moscatel' in the Iberian Peninsula appeared in a document of Alonso de Herrera (1645) and the first unequivocal reference to 'Moscatel' ('Moscatell') in La Marina Alta dates from 1696 (Cabrera, 1992).From this date, an expansion of 'Moscatell' culture was initiated in this area; mainly for raisin production, which achieved high economic importance in the XIX century when exportations to England, France, Canada and the USA were common (Chabás, 1972;Calvo, 2003).After phylloxera (Daktulosphaira vitifoliae Fitch) infestation, the raisin production decreased and many vineyards were lost, being replaced by buildings.Nowadays, the recuperation of ancient clones and selection of differential clones are among the objectives in the restoration of the importance of this crop in La Marina, which belongs to Protected Denomination of Origin (PDO) Alicante, and is located along the coast in the northern area of this province.'Muscat of Alexandria' is cultured in a higher extent at the Comunitat Valenciana in the area comprising the villages of Catadau, Cheste, Chiva, Godelleta, Llombai, Montroi, Montserrat, Real de Montroi and Turís, belonging to PDO Valencia (around 2,750 ha were under cultivation in 2015 in this area; http://www.dovalencia.info),where different sweet and effervescent wines with international prestige are elaborated (Zegels, 2011).'Muscat of Alexandria' is also an important variety in PDO Málaga, PDO Jeréz, and PDO Condado de Huelva.Vinification of 'Muscat of Alexandria' was also recently authorized in most Spanish Communities.In Spain, the name for 'Muscat of Alexandria' is 'Moscatel de Alejandría' (as synonymies: 'Moscatel', 'Moscatel Romano', or 'Moscatel de Grano Grueso') and in the Comunitat Valenciana it is 'Moscatell' or 'Moscatell romà' (DGAIC, 1891;Favà, 2001).
The aim of this work is to study the genetic variability of 'Muscat' germplasm which include accessions from neglected vineyards from Alicante and Valencia provinces, grapevine samples from 'Colección de Vides de El Encín' and from 'La Casa de las Vides', accessions supplied by nurseries and, 'Muscat of Alexandria' clones selected using differential ampelographic characteristics in selection programs.For this purpose fifteen SSRs (Short Sequences Repeats) will be used to study intervarietal variability and AFLPs (Amplified Fragment Length Polymorphisms) and M-AFLPs (Microsatellite-Amplified Fragment Length Polymorphisms) markers to study intravarietal variability in 'Muscat of Alexandria' clones.Finally, the use of the molecular marker developed for the VvDXS gene, described as useful to discriminate 'Muscat' flavor, will be evaluated.

Material and methods
A total of 40 grapevine samples of 'Muscat' were used in this work.The nomenclature and origin of these materials are indicated in Table 1.In addition, an accession of 'Monastrell' has been used as outgroup in the intervarietal analysis of 'Muscat' germplasm and an accession of 'Pampolat' was included in the analysis of VvDXS gene.'Muscat germplasm' collected in neglected vineyards (V5, V6, V7, V10 and V13) correspond to old isolated vines, some of them putatively planted before phylloxera attack (V13).

DNA extraction and SSR analysis
Fully-expanded leaves from 40 'Muscat' accessions were used for DNA extraction, using the commercial DNeasy Plant Mini Kit (Qiagen) according to the manufacturer's instructions.The DNA quality and quantity were assessed using gel electrophoresis and spectrophotometry.
Fifteen SSR markers (VVS2, VVMD5, VVMD6, VVMD7, VVMD21, VVMD24, VVMD25, VVMD27, VVMD28, VVMD32, VrZAG62, VRzAG64, VrZAG79, VrZAG83, and VMC1b11) were analyzed using two sets of multiplex PCR reactions (Table 2).Each multiplex was carefully assembled according to the compatibility of the SSRs during PCR and the molecular size of their amplicons.The forward primers of the SSR markers were labeled with one of the four fluorescent dyes: carboxy fluorescein (FAM), carboxytetramethylrhodamine (TAMRA), hexachloro-6-carboxyfluorescein (HEX), or 6-carboxytetramethyl rhodamine (ROX) (Table 2).Multiplex PCRs were carried out in a total volume of 11.00 µL, using 1.25 μL of commercial Master Mix PCR Multiplex (Takara Multiplex Hot Short PCR, Takara), 20-40 ng of genomic DNA, 0.5 U of Taq DNA Polymerase (Takara HotStart PCR, Takara) and labeled multiplexed SSR primers (from 5.5 to 35 pmol, Table 2).The amplification was performed in an ABI 9700 thermocycler, and the amplification conditions were 95 °C for 14 min followed by 30 cycles of 95 °C for 30 s, 55 °C for 90 s, and 72 °C for 60 s, and a final extension of 72 °C for 30 min.Previous to PCR fragment size determination, the multiplex PCR products were previsualized using gel electrophoresis.The electrophoresis was carried out on an ABI 3100 platform (Appl.Biosyst., Foster City, CA, USA).For PCR fragment size determinations, 0.13 µL of an internal size standard (GeneScan TM 500 LIZ, Appl.Biosyst.) was mixed with 1 µL of PCR product and 10.87 µL of formamide.The mixture was heated at 94 °C for 3 min and then cooled in icy water.The size of the SSR fragments was determined with the software package GeneScan 3.7 (Appl.Biosyst.).
To estimate the discriminatory power of the microsatellite loci, the Polymorphic Information  The AFLP analysis was performed using as forward primer MseI and as reverse primer EcoRI (Table 3).Each 20-μL PCR reaction contained 0.5 μL of the preamplified DNA, 50 ng of labeled EcoRI+3, 30 ng of unlabeled MseI+3 primer, 2 μL of 10× PCR buffer, 4 mM dNTPs, and 0.4 U of Taq DNA Polymerase.The cycling conditions of labeled-PCR were 1 cycle of 30 s at 94 °C, 30 s at 65 °C, and 1 min at 72 °C and a touchdown profile (11 cycles with -0.7 °C/cycle, annealing temperature) for the annealing step, followed by 23 cycles at 56 °C constant annealing temperature and, finally, an extension cycle of 30 min at 60 °C.
The M-AFLP analysis was performed using a procedure identical to that adopted for the AFLP analysis reported above, except for one primer type used in the second amplification, which was the labeled SSR in combination with an MseI+3 primer.
Similar to the SSR analysis, the electrophoresis was carried out on an ABI 3100 platform (Appl.Content (PIC) value for each locus was estimated by where p i is the frequency of the i th allele and n is the number of alleles (Botstein et al., 1980).Nei (1978) genetic similarities were calculated and an Unweighted Pair Group Method with Arithmetic mean (UPGMA) phenogram was produced using genetic distances, with PowerMaker software (Liu & Muse, 2005), and plotted using TreeView software (Page, 2011).
The SSR profiles obtained for the analyzed materials were compared with different databases or published works (i.e. the VIVC database).Before the comparison, normalization of the SSRs using common varieties was performed.
Amplified fragments derived from the AFLP and M-AFLP analyses were evaluated using the Genographer program v.2.1.4(Banks & Benham, 2008) and only polymorphic, distinct, reproducible, and wellresolved fragments were used.These fragments were scored according to the presence (1) or absence (0) of bands and were then transformed into a binary matrix.The similarity index was estimated using the Dice coefficient of similarity (Nei & Li, 1979): Sij = 2a / (2a + b + c), where Sij is the similarity between two individuals 'i' and 'j', 'a' is the number of bands shared by 'i' and 'j', 'b' is the number of bands amplified exclusively in 'i', and 'c' is the number of bands amplified exclusively in 'j'.Subsequently, cluster analyses were carried out using Unweighted Pair Group Method with Arithmetic mean (UPGMA), with the PHYLIP software package v.3.69 (Felsenstein, 2008).To verify the robustness of the nodes, resampling of the matrix with 1,000 samples and a replacement of 35% of the data was performed.The dendrogram was visualized with the program TreeViewPPC v.1.6.6 (Page, 2011).
A Cleaved Amplified Polymorphic Sequences (CAPS) marker developed by Emanuelli et al. (2014) was used to detect SNP1822 G>T.Briefly, the PCR mixture (20 μL) contained 5-10 ng of genomic DNA, 1.25 μL of commercial Master Mix PCR Multiplex (Takara Multiplex Hot Short PCR, Takara), 40 μM of each dNTP, 0.6 μM of each primer, and 0.5 U of HotStar Taq DNA Polymerase (Takara).Amplification was carried out using an ABI 9700 thermocycler and a touchdown protocol.Thermocycling consisted of an initial denaturation of the template DNA at 95 °C for 15 min, followed by 11 cycles of 95 °C for 45 s, 62 °C (touchdown step from 62 °C to 57 °C) for 45 s, and 72 °C for 1 min, and another 25 cycles of 95°C for 45 s, 57 °C for 45 s, and 72 °C for 1 min, with a final extension of 10 min at 72 °C.An aliquot (5 μL) of each amplicon was examined by agarosegel electrophoresis, stained with ethidium bromide, and photographed using a gel documentation system (UVItec).Then, the VvDXS PCR products were digested using FastDigest Eco147I (Fermentas), according to the manufacturer's protocol.Each digest was examined by agarose-gel electrophoresis: the profiles were detected upon ultraviolet transillumination and photographed.
The PCR amplified products from the 'Moscatel de Alejandría', 'Moscatel de Uzbezquistan', 'Moscatel Dorado' and 'Pampolat' accessions were purified, using a GeneJET Gel Extraction Kit, and sequenced in both directions by standard Sanger sequencing.The nucleotide sequences were aligned and compared with sequences available for the VvDXS gene of grapevine accessions retrieved from GenBank.
The profile obtained for V20 was compared with those found in the database of the VIVC and was confirmed as 'Moscato Giallo', with a similar SSR profile for the comparable SSRs (VVS2, VVMD5, VVMD7, VVMD25, VVMD27, VVMD28, VVMD32, VrZAG62, and VrZAG79).Of the germplasm analyzed, the genotype that differed most was 'Moscatel Dorado', differing from the other germplasm analyzed as much as the outgroup 'Monastrell', followed by M5 ('Muscadelle') and M9 ('Moscatuel').
Table 4 displays the molecular diversity analysis of the 21 'Muscat' accessions showing a specific SSR profile.The average number of genotypes per locus was 9.1, ranging from 6 to 14.The overall observed and expected heterozygosity values were 0.816 and 0.719, respectively.Therefore, the inbreeding coefficient (F) was -0.135.All the accessions showed the VVMD5 SSR in heterozygosity, whereas the VVD24 SSR showed heterozygosity close to 0.50.The estimated frequency of null alleles (NAI) was negative for 13 loci and positive for the other two (VrZAG62 and VVMD21).Almost all the SSRs analyzed were highly informative  1.
since they showed a PIC value ≥ 0.5; on average, 0.682.The Shannon´s index (I) value was close to 1.5.The percentage of loci showing minor allele frequency (MiAF) values > 0.1 was about 0.156 and the percentage of loci showing major allele frequency (MaAF) values was 0.422.

Analysis of genetic variability in 'Muscat of Alexandria'
The intravarietal variability analysis was performed in 10 accessions grouped with the 'Moscatel de Alexandria' reference, using AFLPs and M-AFLPs and including as the outgroup the accession M6 ('Moscatel Ruso').A total of 571 reproducible amplification products were obtained (481 DNA fragments from AFLPs and 90 DNA fragments from M-AFLPs).Out of these, 284 (49.73%) were polymorphic: 220 AFLPs and 64 M-AFLPs.Therefore, the M-AFLPs molecular polymorphisms were more efficient with regard to discriminating the accessions with some plant-specific polymorphisms while the AFLPs showed many monomorphic markers (71.11% vs 45.73%).
A total of 114 polymorphic bands were found in at least six genotypes, while a total of 77 polymorphic bands were specific for one genotype.As expected, the genotype M6 ('Moscatel Ruso') had the highest number of specific bands (46), followed by V15 ('Moscatel'), with 15 bands.The rest of the genotypes had at least one specific band.
The dendrogram results are in accordance with those of the Principal Coordinate Analysis (PCoA), which are shown in Fig. S1 [suppl].The first two principal components account for 41.5 and 19.8% of the total variation, respectively.The results of the AFLP analysis and the M-AFLP analysis are in agreement.

Amplification of the VvDXS gene
We amplified the corresponding segment of the CAPs marker developed by Emanuelli et al. (2014) 1. 'Moscatel Ruso' (M6) was used as outgroup.
All the samples analyzed showed, after cutting with Eco147I, similar banding patterns (the heterozygote G/T), which corresponds to 'Muscat' flavor in Muscat accessions.Sanger's sequencing confirmed the restriction analysis data.
V10 ('Moscatel' ST) and V11 ('Moscatel de Alejandría' clone 11) can be considered a variant of 'Muscat of Alexandria', as in other cases in which two plants showed identical SSR profiles for all the SSR markers studied except for one or two alleles.This In (E) and (F), typology of vines; corresponding to the older clones (V11, V59, or V61) which are less vigorous than the more recent clones,V16 and V30, which have also bigger bunches.Accessions are described in Table 1.

Intravarietal genetic variability in 'Muscat of Alexandria'
In our work SSR markers have identified some intravarietal variability among the analyzed 'Muscat' varieties as occurred in other grapevine studies i.e. in the cultivar 'Nero d'Avola' (Carimi et al., 2011).However, intravarietal studies using AFLP and M-AFLP gave higher resolution and separated 'Muscat of Alexandria' clones and highlighted the existence of genetic differences among them, as reported in other grapevine varieties (Cabezas et al., 2003;Cretazzo et al., 2010;Meneghetti et al., 2011Meneghetti et al., , 2012)).Random amplified polymorphisms (RAPDs) also have been used to analyze genetic diversity in Muscat germplasm (Fanizza et al., 2000).
The subgroups found for the clones selected in the 2007TAHALI00073 project are in agreement with their origin.The 'Moscatel de Alejandría' accessions V11 (clone 11), V59 (clone 59) and V61 (clone 61) corresponded to a group of similar characteristics that may be the most ancient clones of 'Moscatel de Alejandría' in 'La Marina Alta' (that probably arrived from Málaga).Among them, V59 had whiter grapes (Fig. 3C).These three clones are easily differentiated by farmers from the most recent clones V16 ('Moscatel de Alejandría' clone 16) and V30 ('Moscatel de Alejandría' clone 30), both named as 'Malagueños', which correspond to another accession of 'Muscat of Alexandria' from Málaga.The older clones (11, 59 and 61) have lower height and vigor than the clones 16 and 30 introduced most recently (Fig. 3E and 3F).Respect the clone V8 ('Moscatel' JR), with red berries (Fig. 3A), other mutations may be accumulated in addition to a color sport.This accession clustered with the clones of the 2007TAHALI00073 project but is the most different.
The greater polymorphism obtained with M-AFLPs (71.11% vs 45.73% with AFLPs) was observed also by Meneguetti et al. (2012) when analyzing a total of 30 Italian and Croatian Istrian 'Malvasia' genotypes.Both frequencies of polymorphic bands are higher than those obtained by Fanizza et al. (2000) in Muscat germplasm.
The genetic differences found among 'Muscat of Alexandria' accessions and clones selected in the project 2007TAHALI00073 are very promising.Firstly, it confirms the interest of this project, in which ancient clones of 'Muscat of Alexandria' were selected.Secondly, these materials represent autochthonous biotypes, which are already both selected and adapted to their environments, in this case, to the climate and edaphological conditions of 'La Marina Alta' at the Comunitat Valenciana.Different works highlight the preservation of this variability (Cretazzo et al., 2010;Meneghetti et al., 2011Meneghetti et al., , 2012)).

Amplification of the VvDXS gene
Finally, we tested the CAPs marker described by Emanuelli et al. (2014) for the discrimination of 'Muscat'-flavored genotypes, since a single nucleotide polymorphism within VvDXS (SNP1822 G>T) causes a dominant gain of function K284N substitution (Emanuelli et al., 2010).This marker discriminates 'Muscat'-flavored genotypes (homozygotes T/T and heterozygotes G/T) from non-aromatic samples (G/G) in Muscat germplasm.Similar banding patterns were found in our work for all the samples analyzed, including 'Moscatel de Uzbequistan'.In the work of Emanuelli et al. (2014), 'Muscat Uzbekistanskii' (acc.1351) was classified as non-aromatic and a different pattern was obtained for this marker.The comparison of SSR profiles for 'Moscatel de Uzbequistan' from 'La Casa de la Vides' and 'Muscat d'Ouzbekistan' 2647 of Lacombe et al. (2013) gave the same profiles for the ten comparable SSRs.In addition, the polymorphism G/T was confirmed using Sanger's sequencing.Probably the accession of 'Muscat Uzbekistanskii' used in Emanuelli et al. (2014) differed from that used in this work.
In conclusion, inter and intravarietal genetic variation was found among the 'Muscat' analyzed accessions.The SSR fingerprinting also allowed the identification of some accessions, variants, synonymies, and mislabeling.The clustering of accessions in relation to their genetic distance allowed the discussion of their relationships.All the accessions analyzed showed the SNP1822 G/T pattern corresponding to 'Muscat' flavor genotypes.The variability found among the germplasm from public or private collections is very interesting: variants were identified in 'Muscat of Alexandria', 'Muscat Italia', and 'Muscat d'Istambul', and variability was found among accessions of 'Muscat of Alexandria'.In the intravarietal analysis, M-AFLPs yielded more polymorphisms than AFLPs.The use of variability is needed to face up to environmental changes and threats from pathogens.The evaluation of variability and the identification of intravarietal variability are of great interest in PDOs because they allow the choice of the most adequate varieties or clones adapted to specific environments, contribute to diversification of wines, and avoid genetic erosion.
magnesium acetate, 100 mM potassium acetate, 5 mM dithiothreitol, 2.5 mg of bovine serum albumin).Then, the pre-amplification was performed using 5 μL of DNA (seven-fold diluted, digested, and ligated) in 20 μL of reaction mixture containing 75 ng of EcoRI+N and MseI+N primers, 1× PCR buffer (50 mM KCl, 1.5 mM MgCl 2 , 10 mM Tris-HCl), 10 mM dNTPs, and 1 U of Taq DNA polymerase (Takara Clontech).Restriction enzymes were provided by Thermo Scientific and primers and adapters were provided by Metabion GmbH.The pre-amplification conditions were an initial denaturation of 5 min at 95 °C, 1 cycle of 45 s at 94 °C, 35 s at 65 °C, and 1 min at 72 °C, and a touchdown profile (13 cycles with -0.7 °C/cycle, annealing temperature) for the annealing step, followed by 18 cycles at 55 °C constant annealing temperature and, finally, an extension cycle of 5 min at 72 °C.

Figure 1 .
Figure 1.UPGMA (Unweighted Pair Group Method with Arithmetic mean) clustering resulting from the analysis of 40 'Muscat' accessions using 15 SSR markers.A 'Monastrell' accession was included as outgroup.Accessions are described in Table1.

Figure 3 .
Figure 3. Detail of grapes from 'Muscat of Alexandria' clones V8 (A), V11 (B), V59 (C), and V61 (D) which differed in the color and compactness.In (E) and (F), typology of vines; corresponding to the older clones (V11, V59, or V61) which are less vigorous than the more recent clones,V16 and V30, which have also bigger bunches.Accessions are described in Table1.

Table 1 .
'Muscat' accessions analyzed in the present study: their berry skin' color, code and origin.

Table 2 .
PCR primers used in the present study for the amplification of SSR markers.

Table 3 .
Primer combinations used for AFLP and M-AFLP amplifications.