Genetic diversity of Poa pratensis L. depending on geographical origin and compared with genetic markers

Introduction

Poa pratensis L. is one of the most common species of meadow grass in Europe. It occurs throughout the country and is highly valued as a grass cultivar that can be used for various purposes (Mirek & Piȩknoś-Mirkowa, 2007; Szenejko, 2014). Due to its high plasticity, high spreading capability and significant expansiveness, P. pratensis is considered a cosmopolitan species capable of occupying a wide range of different habitats. Its wide range of distribution and the ease with which the species adapts to very different conditions has resulted in a great diversity of ecotypes thriving in a miscellany of habitats. In the case of P. pratensis, the selection of particular ecotypes became the basis for breeding work. In Poland, most cultivars of the species were derived from natural grasslands, often from individual plants. Because apomixis is the dominant manner of reproduction in P. pratensis, achieving high genetic variability within individual populations is rather difficult (Müntzing, 1933; Muller, 1964; Felsenstein, 1974; Pamilo, Nei & Li, 1987; Huff & Bara, 1993). In addition to vegetative and sexual reproduction processes, the plant can also reproduce asexually from seeds, which are formed through apospory or diploid parthenogenesis and without any fertilization process (Mazzucato, Den Nijs & Falcinelli, 1996; Albertini et al., 2001; Spillane, Steimer & Grossniklaus, 2001; Albertini et al., 2005; Matzk et al., 2005). As a result, offspring and parent plants may be genetically identical (Albertini et al., 2004; Carneiro, Dusi & Ortiz, 2006). This implies that the initial forms i.e., the different ecotypes and cultivars used for breeding did not produce new combinations of genes. As a result, cultivars derived from such initial material have similar or related genotypes. The results of our earlier research, which concerned the evaluation of genetic variability in cultivars and ecotypes of P. pratensis native to Central Europe using RAPD markers, confirmed the considerable genetic relatedness of those cultivars and ecotypes (Szenejko, Filip & Słominska-Walkowiak, 2009; Szenejko & Rogalski, 2015). The populations evaluated in those studies originated from Lower Silesia and Podlasie, regions of Poland with different habitats and climatic conditions. The ecotypes analyzed in those previous studies differed with respect to their phenotypic traits but showed little genetic variability, whereas the populations in the present study were selected from habitats as diverse as possible, i.e., lowlands, uplands and mountains. The ecotypes in this study originated from the western and southern regions of Poland as well as from the Kujawy and Podlasie regions. It was assumed that long distances between regions would affect habitat conditions, and that natural spatial barriers would cause the greatest genotypic diversity among the research material. In addition, three cultivars (including two of the oldest Polish cultivars of P. pratensis) were included in our research and subjected to the same analyses. All the plant material was evaluated for genetic diversity using two systems of markers based on DNA amplification (RAPD and ISSR). In the present study, the effectiveness of the methods in assessing the genetic diversity of the selected forms of P. pratensis was compared. Furthermore, we investigated whether these methods could be useful in establishing a link between the geographical origin of a given population and their assessed genetic variability. Primers with the greatest differentiating powers correlating with geographical distance were selected for ISSR, the more effective method in that respect. Principal Component Analysis (PCA) was used for this procedure, which was performed on the chosen values of DNA amplification products obtained in the presence of those selected ISSR markers with the highest genetic differentiating power with respect to the studied forms of P. pratensis.

Materials & Methods

Plant material

This assessment of the genetic diversity of 18 selected forms of P. pratensis was carried out in 2014. The study involved 3 cultivars and 15 ecotypes selected from locations as distant as possible from each other. The seeds used in the study came from the national collection of the Botanical Garden of the Institute of Plant Breeding and Acclimatization in Bydgoszcz (IHAR); the caryopses used for establishing this collection were collected during field trips organized by the staff of the Botanical Garden and originated from different geographical regions of Poland, including 8 different voivodships (provinces) (Fig. 1 and Table 1). The main criteria when harvesting the caryopses were the diversity of plant material and variety of habitats. The seed material was collected from wastelands and arable land (meadows and pastures), as well as from environmentally valuable areas. The assessment of genetic diversity was performed for three cultivars used for different purposes, including ‘Eska 46,’ one of the oldest Polish fodder cultivar and two lawn cultivars: ‘Limousine’, an apomictic uniclonal German cultivar, and ‘Alicia,’ the oldest Polish lawn cultivar.

Table 1:

Ecotypes of P. pratensis analyzed in the study and the locations where their caryopses were collected.

Ecotype	Town/Place	Latitude	Longitude	Habitat
pom06 175	Wielkie Stwolno	53°26′34.1″N	18°39′25.8″E	Vistula River oxbow lake
zap00 240	Łukȩcin	54°02′58.4″N	14°52′03.7″E	Beach
lbs05 170	Owczary	52°28′16.3″N	14°39′08.0″E	Ostnicowy nature reserve
lbs05 628	Przemków	51°32′50.6″N	15°46′33.4″E	Meadow
dos01 333	Rybnica Leśna II	50°42′08.5″N	16°18′00.3″E	Pasture
dos01 405	Dziwiszów–Góry Kaczawskie	50°57′39.4″N	15°49′23.6″E	Meadow
dos01 468	Oleszna Podgórska	51°21′02.3″N	15°28′12.1″E	Meadow
bes03 373	Ustroń–Wielka Czantoria	49°40′45.6″N	18°47′42.7″E	Mountain meadow
bes03 379	Ustroń–Polana Goryczkowa	49°41′05.6″N	18°46′58.6″E	Meadow pasture
bes03 354	Jaworzynka	49°31′20.2″N	18°51′31.2″E	Meadow
ora03 178	Lipnica Wielka	49°27′05.3″N	19°37′40.6″E	Meadow
bie04 189	Płonna	49°26′11.7″N	22°07′05.9″E	Wasteland
pod02 215	Lipina	53°19′59.6″N	23°25′28.5″E	Meadow/pasture
pod02 172	Walne	53°58′20.0″N	23°04′04.6″E	Meadow
pod02 317	Pieńki	53°03′32.3″N	23°38′47.2″E	Pasture

DOI: 10.7717/peerj.2489/table-1

Results

RAPD analysis

In the course of RAPD analysis, in total 246 amplification products were obtained for which 100% polymorphism was demonstrated (Table 3). The primers chosen to evaluate polymorphism generated bands in a wide size range: from 149 bp to 2,367 bp in the presence of primer C16. The greatest number of amplicons (43) was obtained as the result of reactions with primers G19 and M14 (Fig. 2). The primers used for the analysis initiated the synthesis of 79 specific amplification products which differentiated the studied forms of P. pratensis. Of these, the largest number of DNA fragments (16) which had a unique composition were obtained in reaction with the M14 primer (Fig. 2). As many as 68 specific products were found for the ecotypes, with the largest number of them in a population from Silesia (bes03 373). In the presence of six primers, we were able to generate 15 products with a unique composition, ranging in size from 149 bp (for C16) to 1478 bp (for M14). Unique DNA fragments of approximately the same size were identified in amplification reactions with different primers. A significant number of specific bands were also identified for two ecotypes, one from Podkarpacie (bie04 189) and the other originating from the borderland region of Lubuskie and Lower Silesian voivodships (provinces) (lbs05 628) 11 and 10 bands respectively.

Table 3:

RAPD and ISSR primers used, total number of scored bands (NSB), number of polymorphic bands (NPB), their unique bands (NUB) and percentage of polymorphism bands (PPB) for each primer.

Name of primer	Size range (bp)	NSB	NPB	NUB	PPB
RAPD
B10	213–1,867	35	35	11	100
B17	209–1,483	32	32	7	100
C16	149–2,367	33	33	15	100
F05	316–2,366	34	34	11	100
G19	201–1,902	43	43	9	100
M14	226–2,243	43	43	16	100
P08	179–1,173	26	26	10	100
Total	149–2,367	246	246	79	–
ISSR
807	341–1,974	23	23	1	100
808	325–3,105	43	43	17	100
810	338–2,523	35	35	8	100
811	303–2,748	40	40	10	100
834	404–2,686	35	34	7	97
840	373–2,967	50	50	8	100
841	373–2,504	37	36	7	97
857	380–2,577	46	45	7	98
888	501–2,414	26	24	5	92
R07	576–2,327	37	37	13	100
M09	353–2,096	47	47	9	100
N01	404–2,256	51	51	15	100
TGA	362–2,209	44	43	17	98
Total	303–3,105	514	508	124	–

DOI: 10.7717/peerj.2489/table-3

ISSR analysis

As the result of amplification reactions involving 13 microsatellite primers, a total of 514 products ranging in size from 303 bp (for 811) up to 3,105 bp (for 808) were obtained (Table 3 and Fig. 3). In reactions with primers N01 and 840 a large number of bands were generated (51 and 50 respectively), whereas in the presence of primer 807 only 23 products (the smallest number) were obtained. A high level of polymorphism was observed in the amplified DNA fragments and in the presence of eight primers only polymorphic products were obtained. We identified a total of 124 specific amplification products in the ecotypes and cultivars of P. pratensis analyzed in this study. The largest number (17) was obtained in reaction with primer 808, while in the presence of primer 807 only one product with a unique composition (649 bp) for ecotype lbs05 170 was obtained (Fig. 3). In total, 95 such products were identified for the ecotypes. The largest numbers were recorded for populations lbs05 170 (15 bands) and lbs05 628 (13 bands), and for ecotype pom06 175. Specific products of similar molecular size were identified in the presence of different primers. A significant number of unique bands were also identified for the lawn cultivars, including the German ‘Limousine’ (11).

Evaluation of RAPD and ISSR polymorphisms

Using RAPD and ISSR markers, a significant number of specific amplification products were obtained. 98.6% of the bands obtained using ISSR were polymorphic, whereas those obtained in RAPD reactions were 100% polymorphic (Table 3). A large number of products were obtained using both methods in the presence of single primers (on average 30 products for each method) (Table 4). For both RAPD and ISSR methods the polymorphism information content (PIC) was calculated; the average PIC value for each method was almost identical and was 0.264 for RAPD and 0.270 for ISSR. The value of the assay efficiency index, which indicates the average number of polymorphic products identified in the presence of a single primer, was also calculated for both RAPD and ISSR marker systems. The obtained AEI value was higher for ISSR markers. Specific DNA fragments were obtained for almost all the studied populations (Tables 3 and 4); RAPD primers were more effective in identifying them. In a reaction with a single RAPD primer, we were able to identify as many as 11 unique bands. Dice’s genetic similarity matrix was determined for all the studied forms of P. pratensis. Independently of each other, the RAPD and ISSR methods revealed little genetic relationship between the studied populations and cultivars (Table 4). On average Dice’s coefficient for RAPD was 0.371, while it was 0.367 for ISSR markers. The greatest genetic relationship was demonstrated between two populations from Podlaskie voivodship, pod02 215 and pod02 317 (0.733 for RAPD, 0.650 for ISSR and 0.678 for RAPD + ISSR). The lowest values of Dice’s coefficient for RAPD (0.154) were obtained for populations lbs05 170 and bie04 189, while in the cases of ISSR and RAPD + ISSR the smallest degree of genetic relatedness was found between ecotype lbs05 628 and the cultivar ‘Alicia’. The following ecotypes were found to be most genetically distant: bie04 189, lbs05 628, and lbs05 170. To highlight the compatibility of the results of the RAPD and ISSR methods, the correlation between them was calculated. Using Spearman’s coefficient (0.461), the statistical correlation between the genetic similarity matrices obtained on the basis of electrophoretic images for RAPD and ISSR was calculated (Fig. 4).

Table 4:

Evaluation of RAPD and ISSR polymorphisms.

The average Number of Scored Bands (NSB) including Number of Unique Bands (NUB) obtained in the presence of a single primer, the value of Polymorphism Information Content (PIC), Assay Efficiency Index (AEI), and Genetic Similarity (Si).

Parameter/Index	RAPD	ISSR	RAPD + ISSR
NSB	35.1	38.0	36.8
NUB	11.3	7.7	9.3
AEI	35.1	37.3	36.4
PIC	0.264	0.270	0.268
Si	0.371	0.367	0.369
(Range)	(0.154–0.733)	(0.218–0.650)	(0.222–0.678)

DOI: 10.7717/peerj.2489/table-4

Cluster analysis

The applied measure of genetic similarity was used to construct UPGMA dendrograms (Fig. 5). There was greater compatibility with respect to the clustered forms for the ISSR and RAPD + ISSR dendrograms. Regardless of which method was used, two asymmetrical and distinct major clusters were clearly distinguishable in each dendrogram. The majority of the studied forms were assigned to the second cluster, within which two subclusters of similarity were distinguishable 11 forms for RAPD and RAPD + ISSR and 12 forms for ISSR. In the ISSR dendrogram, the first main cluster included two ecotypes from Kujawsko-Pomorskie (pom06 175) and from West Pomerania voivodships (zap00 240). The combined RAPD + ISSR analysis also included ecotype dos01 405 from Lower Silesia. In the case of RAPD, the first main cluster contained two other ecotypes (bie04 189 and ora03 178), while pom06 175 and zap00 240 were assigned to the second main cluster and to both common subclusters (Fig. 5). Two ecotypes from Lubuskie voivodship (lbs05 170 and lbs05 628) turned out to be genetically most distant. These were not assigned to any pair or subcluster; this is visible in all the dendrograms (Fig. 5). These populations also differed among themselves; this was confirmed by the value of Dice’s similarity index (0.224 for RAPD + ISSR). In the ISSR and RAPD + ISSR dendrograms, other genetically different forms included an ecotype from Podkarpacie region (bie04 189), while in the RAPD dendrogram, a form from Silesia (bes03 373) was genetically different. Among the cultivars, ‘Eska 46’ was the most distinctly different: in the RAPD and ISSR dendrograms, it was connected loosely with the second cluster, while the image obtained for RAPD + ISSR depicted it as an unpaired form. In each of the dendrograms the ecotypes from Podlasie and the lawn cultivars ‘Alicia’ and ‘Limousine’ formed a common subcluster of forms. In the RAPD dendrogram, the Silesian population bes03 354 joined them, whereas in the ISSR and RAPD + ISSR dendrograms this form was included in the second subcluster together with other ecotypes from Silesia and populations from Lesser Poland and Lower Silesia. In the ISSR diagram, dos01 405 was placed with ecotypes from Lower Silesia, while in the RAPD diagram it was only loosely connected with the second cluster and in the RAPD + ISSR diagram it was included in a different cluster.

Discussion

DNA markers have become a useful tool for breeding work and for evaluating the genetic diversity of many different plant species, including the family Poaceae. Molecular analysis, in conjunction with the characteristics of morphological traits, may be used to classify plants and select the initial forms for creating new cultivars or breeding lines. In the case of P. pratensis, this was confirmed by Shortell, Meyer & Bonos (2009), Tamkoc & Arslan (2010), Raggi et al. (2015) and Yuan et al. (2015). In this study, two markers (RAPD and ISSR) based on DNA amplification reactions were used to assess the genetic diversity of ecotypes and cultivars of P. pratensis. These methods do not require any prior information on the analyzed DNA sequence. They do not take a lot of time and are also relatively simple and therefore quite universal. Furthermore they require only a small amount of DNA; this is why they are often employed in the genetic evaluation of different plant species. They can be used either in combination (Nkongolo, Michael & Demers, 2005; Alam, Naik & Mishra, 2009; De Lima et al., 2011; Singh et al., 2011; Debajit et al., 2015) or independently of each other (Al-Humaid & Motawei, 2004; Wang, 2010; Guo et al., 2009; Bednarskaya et al., 2014). In grasses, RAPD and ISSR markers are most commonly used to evaluate genetic variability within and between populations, and for molecular characterization and identification of different species and their hybrids, cultivars and genotypes (Posselt et al., 2006; Pivorienė et al., 2008; Al-Humaid, Ibrahim & Motawei, 2011; Motawei & Al-Ghumaiz, 2012; Madesis et al., 2014; Yuan et al., 2014). In P. pratensis, genetic analyses employing RAPD and ISSR focus mainly on the identification of sexual and apomictic genotypes (Huff & Bara, 1993; Mazzucato et al., 1995; Barcaccia et al., 1997; Barcaccia, Veronesi & Falcinelli, 1998; Stephens et al., 2006); the assessment of genetic variability and identification of genotypes and cultivars (Lickfeldt, Voigt & Hamblin, 2002; Ning et al., 2005; Liang et al., 2009; Fard et al., 2012; Wang et al., 2012; Yuan et al., 2015); determining the genetic relationships between different genotypes and species belonging to the genus Poa and their hybrids (Johnson et al., 2002; Curley & Jung, 2004; Patterson, Larson & Johnson, 2005; Goldman, 2008; Goldman, 2013); and identification of pathogens responsible for fungal diseases in grasses (Hsiang et al., 2000). Both marker systems can generate high resolution band patterns and high levels of polymorphism as high as 90%. The presence of these markers resulted in a significant number of amplification products in our study: 246 for RAPD and 514 for ISSR. However, ISSR gives better reproducibility of results, which can be explained by its use of longer primers (microsatellites) and higher annealing temperatures (Gostimsky, Kokaiva & Konovalov, 2005; Reddy, Sarla & Siddiq, 2002; Fernández, Figueiras & Benito, 2002; Vaillancourt et al., 2008). Similarly, in previous studies concerning the assessment of the morphological and genetic variability of domestic ecotypes of P. pratensis originating from the areas of Lower Silesia and Podlasie (Szenejko & Rogalski, 2015), RAPD primers F05, G19 and M14 proved to be most useful in the evaluation of polymorphisms. The highest number of DNA fragments were obtained in the presence of these primers. In the case of ISSR markers, as in the studies of Bednarskaya et al. (2014), better results were obtained in the presence of primers with dinucleotide repeat sequences (AG)₂,(GA)₂ and (AC)₂, anchored at the end of 3’. As suggested by Al-Humaid et al. (2004) and Al-Humaid, Ibrahim & Motawei (2011), dinucleotide and trinucleotide repeats are abundant in genomes of fodder grasses and thus might constitute potential markers for assessing their genetic diversity and identifying their fungal diseases using the ISSR method. In grasses, (CAC)_n, (CTC)_n and (GA)_n and (CA)_n are particularly abundant. In the genus Poa, ISSR primers with dinucleotide (AC)_n (CA)_n (GA)_n and (AG)_n repeats turned out to be particularly useful for the identification and evaluation of genetic relationships between various species, including P. angustifolia, P. trivialis, P. arachnifera, P. pratensis, P. compressa, P. ligularis, and P. secunda and their hybrids (Goldman, 2008; Arslan & Tamkoç, 2011). A high level of polymorphism was observed in numerous amplification products generated in the course of our analyses. The nature of ISSR, which targets regions especially rich in microsatellites, may also explain the high level of polymorphism, since those regions are known to accumulate a large number of mutations due to DNA polymerase slippage during replication and unequal crossing-over (Kashi & King, 2006; Da Costa et al., 2011). In the case of ISSR markers, this can lead to longer dinucleotide repeat sequences and insertion into at least one of them. In this way new variants of microsatellite sequences of different lengths appear, which enriches the collection of alleles identified as polymorphic in a population. However, according to many researchers, in the course of RAPD and ISSR analyses significant levels of polymorphism can be generated, and for various forms and species belonging to the genus Poa this can exceed 90% (Liang et al., 2009; Arslan & Tamkoç, 2011; Fard et al., 2012; Wang et al., 2012; Cichorz, Gośka & Litwiniec, 2014). In our analyses, a slightly greater percentage of polymorphic DNA fragments was obtained in reaction with RAPD primers. The data available in the literature indicates that the RAPD method is somewhat superior due to its higher levels of generated polymorphism in the presence of various molecular marker systems which is explained, inter alia, by their polyallelic nature (Alam, Naik & Mishra, 2009; Bhattacharya, Bandopadhyay & Ghosh, 2010; Naik et al., 2010). A possible explanation for the difference in resolution of RAPD and ISSR is that the two-marker techniques target different portions of the genome (Debajit et al., 2015). Although a greater percentage of polymorphic DNA products was identified in the presence of RAPD primers, it was microsatellites that we found to be more effective in detecting them. This finding was confirmed by calculating the average efficiency index (AEI), which was slightly higher for ISSR markers. The slight differences in the assessment of RAPD and ISSR polymorphisms did not rule out the possibility of establishing a correlation between them. Even though Spearman’s coefficient was not very high (0.461), the result was statistically very significant. Hence, it can be concluded that the results presented in this study obtained using RAPD and ISSR markers are comparable. Budak et al. (2004), in a study examining the phylogenetic relationships between seed and vegetative biotypes of Buchloe dactyloides, used various systems of molecular markers (ISSR, SSR, RAPD and SRAP), for which they determined correlations. In the case of RAPD and ISSR, the Spearman’s coefficient (0.461) calculated by Budak was similar to ours, but their results turned out to be statistically insignificant. Both RAPD and ISSR markers revealed a small genetic similarity among the studied ecotypes and cultivars of P. pratensis. The average value of Dice’s coefficient for both methods combined (RAPD + ISSR) was low (0.369), which together with the high number of polymorphic products may indicate substantial genetic differences between the studied populations and their low level of relatedness. It can be speculated that natural spatial barriers might have inhibited and minimized the transfer of genes among these populations; this would certainly have contributed to their diversity. At the same time, we can surmise that there were populations among the studied ecotypes in which apomixis was the main mode of reproduction. The ecotypes originating from three voivodships (Podkarpackie, Lower Silesian and Lubuskie (bie04 189, lbs05 628 and lbs05 170)), were determined to be those genetically most distant from the other forms. These populations had the lowest similarity index and included the most uncommon genotypes, for which the largest number of specific amplification products were identified. In comparison to the other cultivars, ‘Eska 46’ was clearly different in terms of its genetic structure, and, overall, few amplification products were generated for it. It should be stressed that in previous studies aimed at evaluating RAPD polymorphism in different populations of P. pratensis, such a high genetic diversity among ecotypes had never been found, while the average Dice’s coefficient was about 0.7 (Szenejko & Rogalski, 2015). Other authors also reported significant genetic relatedness among the different genotypes and cultivars of P. pratensis (Johnson et al., 2002; Ning et al., 2005; Tamkoc & Arslan, 2010; Fard et al., 2012). UPGMA dendrograms based on genetic similarity matrices revealed discrepancies in the mode of clustering of the studied forms of P. pratensis. Greater similarity in that respect was demonstrated for the ISSR and RAPD + ISSR methods. This difference in the clustering process between the RAPD and ISSR dendrograms was also demonstrated in research by Bhattacharya, Bandopadhyay & Ghosh (2010), who were assessing the genetic diversity of Cymbopogon winterianus and by Souframanien & Gopalakrishna (2004), who were studying genotypes of Vigna mungo. Similarly to our results, these authors found greater clustering conformity for ISSR and RAPD + ISSR. Looking at the distribution in individual UPGMA dendrograms, it appears that the ISSR dendrogram depicts the relationship between the geographical origin of the ecotypes and their genetic similarity in the clearest way. With the exception of populations bie04 189, lbs05 628 and lbs05 170, which differed from the others genetically, the remaining ecotypes were clustered according to their geographical distribution. In the course of PCA, and depending on their geographical origin, a similar tendency in the distribution of the studied populations in two-dimensional space was observed. It was noted that ecotypes from the same voivodships (provinces) or from regions located not far from one another were depicted similarly in two-dimensional space as in the UPGMA dendrogram constructed for the ISSR process. Those populations from eastern Poland and from West Pomerania and Kujawy clearly differed from the ecotypes originating from the western and southern parts of the country and were located on the opposite side of the horizontal axis of coordinates. They also differed from one another, which can be seen in their arrangement on both sides of the vertical axis. Some discrepancies were found between the UPGMA and PCA clusters of the studied forms, which were related, inter alia, to the distribution of the cultivars. Only a modest relationship between genetic divergence and geographical origin in P. pratensis was reported by several authors (Johnson et al., 2002; Fard et al., 2012; Raggi et al., 2015). For comparison, Abbaszade et al. (2013), who were studying the genetic diversity of Lolium multiflorum, and Rahmati et al. (2013), who were studying Festuca arundinacea, demonstrated a greater result consistency between UPGMA and PCA. However, neither of these research teams was able to demonstrate a correlation between the genetic variability of the grass genotypes they analyzed and their geographical origin.

Conclusions

Our results suggest that both the systems of dominant RAPD and ISSR markers that were used in our research can be used to assess the genetic variability of selected ecotypes and cultivars of P. pratensis. However, given their better reproducibility of results, greater efficiency in detecting polymorphic amplification products in the presence of a single primer and also their differentiating properties (which make them capable of distinguishing studied forms with respect to their genetic structure), ISSR markers seem to be a more useful and reliable tool for this kind of analysis. The plant material evaluated in this study was characterized by a high genetic variability, which seems to confirm the hypothesis that ecotypes of P. pratensis originating from different parts of Central Europe, each with different terrain structures and habitat conditions, can provide great genetic variability and a large collection of alleles. Microsatellites with dinucleotide repeat sequences (AG)_n, (GA)_n anchored at the end of 3’ proved particularly useful for assessing the genetic diversity of these ecotypes and determining their main directions of variability. PCA helped to establish a correlation between the geographical origin of the studied ecotypes and their genetic variation.

Supplemental Information