Genetic diversity and structure of natural fragmented Chamaecyparis obtusa populations as revealed by microsatellite markers

The genetic diversity and population structure of hinoki (Chamaecyparis obtusa) in Japan were investigated by examining the distribution of alleles at 13 microsatellite loci in 25 natural populations from Iwaki in northern Japan to Yakushima Island in southern Japan. On average, 26.9 alleles per locus were identified across all populations and 4.0% of the genetic variation was retained among populations (G _ST = 0.040). According to linkage disequilibrium analysis, estimates of effective population size and detected evidence of bottleneck events, the genetic diversity of some populations may have declined as a result of fragmentation and/or over-exploitation. The central populations located in the Chubu district appear to have relatively large effective population sizes, while marginal populations, such as the Yakushima, Kobayashi and Iwaki populations, have smaller effective population sizes and are isolated from the other populations. Microsatellite analysis revealed the genetic uniqueness of the Yakushima population. Although genetic differentiation between populations was low, we detected a gradual cline in the genetic structure and found that locus Cos2619 may be non-neutral with respect to natural selection.     

Conservation genetics of the endangered endemic Sambucus palmensis Link (Sambucaceae) from the Canary Islands

Five polymorphic microsatellites (simple sequence repeat; SSR) markers were used to estimate the levels of genetic variation within and among natural populations from different islands of the endangered endemic from the Canary Islands Sambucus palmensis Link (Sambucaceae). Genetic data were used to infer potential evolutionary processes that could have led to present genetic differentiation among islands. The levels of genetic variability of S. palmensis were considerably high; proportion of polymorphic loci (P = 100%), mean number of alleles per locus (A = 6.8), average expected heterozygosity (He = 0.499). In spite of its small population size and endemic character, 58 different multilocus genotypes were detected within the 165 individuals analyzed. All samples located in different islands always presented different multilocus genotypes. Principal Coordinates Analysis, genetic differentiation analysis (F _ST and G _ST ′) and Bayesian Cluster Analysis revealed significant genetic differences among populations located in different islands. However, this genetic differentiation was not recorded among Tenerife and La Gomera populations, possibly revealing the uncontrolled transfer of material between both islands. AMOVA analysis attributed 77% of the variance to differences within populations, whereas 8% was distributed between islands. The levels of genetic differentiation observed among populations, and the genetic diversity distribution within populations in S. palmensis, indicate that management should aim to conserve as many of the small populations as possible. Concentrating conservation efforts only on the few large populations would result in the likelihood of loss of genetic variability for the species.     

Molecular diversity of small balsam populations in relation to site characteristics

Climatic shifts within recent decades created favorable conditions for invasive species flourishing in more Northern parts of Europe. Our study was aimed at evaluation of genetic variability and habitat features of Impatiens parviflora populations growing in Lithuania. Twenty one populations were selected and analysed using simple sequence repeat (SSR) and randomly amplified polymorphic DNA (RAPD) assays. Evaluated by SSRs, 315 individuals were all monomorphic and homozygous at 4 loci and heterozygous at 1 locus. RAPD analyses revealed that the percentage of polymorphic DNA loci (% P) per population ranged from 7 to 39% and genetic differentiation between populations was Φ_PT=0.790 (P<0.01). Genetic distances among populations (0.135–0.426) correlated significantly with geographical distances (r=0.183; P<0.008). Populations in overmoistured soil contained higher % P (28.3) when compared to drier soil (18.7; P<0.05). All recorded populations were close to roads; their % P did not depend on proximity to buildings, light intensity or population size. Our RAPD analyses indicate multiple introductions of this species in Lithuania. Analyses of I. parviflora at SSR and RAPD loci show that the invasion process is reflected in genetic structure.     

Demographic history of a common pioneer tree,Zanthoxylum ailanthoides,reconstructed using isolation-with-migration model

The Japanese prickly ash, Zanthoxylum ailanthoides Siebold & Zucc. (Rutaceae), is one of the most common pioneer tree species that grow in disturbed areas, such as canopy gaps, in Japanese warm-temperate evergreen oak forests. The strong genetic structure of its current population might suggest long-term isolation of subpopulations by demographic events, in addition to ecological features of this species. Analyses of genome-wide nucleotide variation using model-based approaches are necessary to achieve a good understanding of the demographic history of a species. In this study, we analyzed the nucleotide variation in natural populations of Z. ailanthoides using a computer program that applied the isolation-with-migration model to quantitatively infer the demographic history. Nucleotide variations at 10 or 26 nuclear loci in six populations from a wide range of the species distribution were analyzed. The maximum likelihood estimate of the divergence time between the current populations in the two hypothetical refugia of Z. ailanthoides was approximately 24,000 years. Unexpectedly, the estimated size of the ancestral population was larger than the sizes of the two current populations. The results suggested a relatively recent divergence of these two populations and rapid formation of strong genetic structure among subpopulations. The large ancestral population may indicate a more complex demographic history during or after the last glacial period than the simple isolation-with-migration model implies.     

Population genetic structure of Bellamya aeruginosa (Mollusca: Gastropoda: Viviparidae) in China: weak divergence across large geographic distances

Bellamya aeruginosa is a widely distributed Chinese freshwater snail that is heavily harvested, and its natural habitats are under severe threat due to fragmentation and loss. We were interested whether the large geographic distances between populations and habitat fragmentation have led to population differentiation and reduced genetic diversity in the species. To estimate the genetic diversity and population structure of B. aeruginosa, 277 individuals from 12 populations throughout its distribution range across China were sampled: two populations were sampled from the Yellow River system, eight populations from the Yangtze River system, and two populations from isolated plateau lakes. We used seven microsatellite loci and mitochondrial cytochrome oxidase I sequences to estimate population genetic parameters and test for demographic fluctuations. Our results showed that (1) the genetic diversity of B. aeruginosa was high for both markers in most of the studied populations and effective population sizes appear to be large, (2) only very low and mostly nonsignificant levels of genetic differentiation existed among the 12 populations, gene flow was generally high, and (3) relatively weak geographic structure was detected despite large geographic distances between populations. Further, no isolation by linear or stream distance was found among populations within the Yangtze River system and no signs of population bottlenecks were detected. Gene flow occurred even between far distant populations, possibly as a result of passive dispersal during flooding events, zoochoric dispersal, and/or anthropogenic translocations explaining the lack of stronger differentiation across large geographic distances. The high genetic diversity of B. aeruginosa and the weak population differentiation are likely the results of strong gene flow facilitated by passive dispersal and large population sizes suggesting that the species currently is not of conservation concern.     

Inter-populational difference in microsatellite–centromere map distances in the loach, <Emphasis Type="Italic">Misgurnus anguillicaudatus</Emphasis>

Microsatellite–centromere recombination rates were estimated at 21 loci in relation to centromere of chromosomes in gynogenetic diploid lines induced from loaches Misgurnus anguillicaudatus of two different populations in Japan. All the microsatellite loci gave allelic segregation according to the Mendelian manner of inheritance in normal diploid families. Since loaches from the Kita population in the southern area of Hokkaido Island and those from the Memanbetsu population in the northern area, Hokkaido, Japan, were reported by previous genetic studies to be genetically diversified, map distances were compared between loaches from the two different populations. Three (Mado7, Mac3 and Mac49) of five loci, which could be compared inter-populationally, gave significantly different recombination rates, i.e., map distances. The results support the presence of genetic difference between the two populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetic diversity and structure of Manila clam (Ruditapes philippinarum) populations from Liaodong peninsula revealed by SSR markers

Manila clam (Ruditapes philippinarum) is an important commercially marine bivalve, and its wild populations have been severely declining in the coast of China during the last decade. In this study, a set of 7 genomic simple sequence repeat (SSR), and 5 expressed sequences tag (EST)-derived SSR markers were analyzed on eight wild R. philippinarum populations to assess the genetic diversity and population differentiation. A total of 114 alleles were detected on 12 loci, and the number of alleles per locus in each population ranged from 2 to 11, and allelic richness per locus varied from 2.00 to 10.88 for each sample. The average of observed and expected heterozygosity ranged from 0.386 to 0.550, and from 0.533 to 0.707, respectively. Pairwise F_ST values indicated that all population pairs had significant genetic differentiation (overall F_ST = 0.242, P < 0.01). Cluster analysis using unweighted pair group method with the arithmetic mean (UPGMA) separated the eight populations into two groups. This study will shed light on the domestication and cultivation on population genetic diversity of R. philippinarum, and also provide the foundation for conservation of R. philippinarum germplasm resources in clam breeding practices.     

Characterization of ten polymorphic microsatellite markers for the protected Spanish moon moth <Emphasis Type="Italic">Graellsia isabelae</Emphasis> (Lepidoptera: Saturniidae)

Ten polymorphic microsatellite loci were developed for Graellsia isabelae. Polymorphism was assessed for 20 individuals from a Spanish population (Els-Ports-de-Beseit, Catalonia) and 39 more individuals from one population in the French Alps and six other Spanish localities. Overall, the number of alleles per locus ranged from 5 to 24. Els-Ports-de-Beseit showed an average number of alleles per locus of 9.80 (SD = 4.32), observed heterozygosity was 0.71 (SD = 0.226), and expected heterozygosity was 0.788 (SD = 0.146). Genotypic frequencies conformed to Hardy–Weinberg equilibrium at the Catalonian population, and no evidence for linkage disequilibrium was observed. Multilocus genotypes resulting from this set of markers will be useful to determine genetic diversity and differentiation within and among populations of this highly protected moth. Several loci amplified and resulted polymorphic in two related species: two loci in Actias neidhoeferi, and three loci in A. luna.     

Genetic variation of wild litchi (<Emphasis Type="Italic">Litchi chinensis</Emphasis> Sonn. subsp. <Emphasis Type="Italic">chinensis</Emphasis>) revealed by microsatellites

Understanding the amount and distribution of genetic diversity in natural populations can inform the conservation strategy for the species in question. In this study, genetic variation at eight nuclear microsatellite loci was used to investigate genetic diversity and population structure of wild litchi (Litchi chinensis Sonn. subsp. chinensis). Totally 215 individuals were sampled, representing nine populations of wild litchi. All eight loci were polymorphic, with a total of 51 alleles. The expected heterozygosity in the nine populations ranged from 0.367 to 0.638 with an average value of 0.526. Inbreeding within wild litchi populations was indicated by a strong heterozygote defect. Significant bottleneck events were detected in the populations from Yunnan and Vietnam, which could be responsible for lower levels of genetic diversity in these populations. Measures of genetic differentiation (F _ST = 0.269) indicated strong differentiation among wild litchi populations. Significant correlation was found between genetic differentiation and geographical distance (r = 0.655, P = 0.002), indicating a strong isolation by distance in these populations. Bayesian clustering suggested genetic separation among three regional groups, namely, the western group, the central group and the eastern group. Some conservation strategies for wild litchi populations were also proposed based on our results.     

Analysis of genetic diversity in the endangered Hucho taimen from China

Hucho taimen are listed as endangered in China. The population size has declined recently, prompting an increase in the level of listing from grade three in 2002 to grade five in 2006. We analyzed the genetic diversity of wild populations using 17 microsatellite markers to establish a scientific basis for conservation of this species. We collected tissue samples from four populations in the Heilongjiang River basin: Huma River (HM), Hutou (HT), Haiqing (HQ), and Zhuaji (ZJ). A total of 21 loci were amplified, 18 of which were polymorphic. The number of alleles per locus ranged from 2 to 9 (mean: 4.1905). There were 13 highly polymorphic loci and 5 moderately polymorphic loci. Analysis of five genetic diversity parameters (Na, Ne, Ho, He, and PIC) suggested moderate levels of diversity within the populations. The populations were ranked HT > HQ > ZJ > HM, but the differences in diversity were not statistically significant (P > 0.05). A comparison of variation among all four populations suggested Hardy–Weinberg disequilibrium at 20% of the loci. Genetic differentiation (Fst) was 0.0644 and the gene flow among populations was estimated at 3.36 individuals per generation. The majority of diversity (93.88%) occurred among individuals within a population. In contrast, relatively little (6.12%) of the genetic diversity was distributed between the populations. An analysis of genetic differentiation and genetic distance between pairs of populations revealed that both parameters were higher in comparisons of the HM population to the HT, HQ, and ZJ populations than among the three latter populations. This suggests that the HM population has a distinct genetic structure. We hypothesize that habitat degradation and excessive fishing, not low genetic diversity, has caused the decline in H. taimen populations. However, this species should be protected from further declines in genetic diversity.     

Genetic variation and clonal structure of the rare,riparian shrub <Emphasis Type="Italic">Spiraea virginiana</Emphasis> (Rosaceae)

Genetic diversity is often considered important for species that inhabit highly disturbed environments to allow for adaptation. Many variables affect levels of genetic variation; however, the two most influential variables are population size and type of reproduction. When analyzed separately, both small population size and asexual reproduction can lead to reductions in genetic variation, although the exact nature of which can be contrasting. Genetic variables such as allelic richness, heterozygosity, inbreeding coefficient, and population differentiation have opposite predictions depending upon the trait (rarity or clonality) examined. The goal of this study was to quantify genetic variation and population differentiation in a species that resides in a highly stochastic environment and is both rare and highly clonal, Spiraea virginiana, and to determine if one trait is more influential genetically than the other. From populations sampled throughout the natural range of S. virginiana, we used microsatellite loci to estimate overall genetic variation. We also calculated clonal structure within populations, which included genotypic richness, evenness, and diversity. Gene flow was investigated by quantifying the relationship between genetic and geographic distances, and population differentiation (θ) among populations. Observed heterozygosity, genotypic richness, and inbreeding coefficients were found to be representative of high clonal reproduction (averaging 0.505, 0.1, and –0.356, respectively) and the number of alleles within populations was low (range = 2.0–3.6), being more indicative of rarity. Population differentiation (θ) among populations was high (average = 0.302) and there was no relationship between genetic and geographic distances. By examining a species that exhibits two traits that both can lead to reduced genetic variation, we may find an enhanced urgency for conservation. Accurate demographic counts of clonal species are rarely, if ever, possible and genetic exploration for every species is not feasible. Therefore, the conclusions in this study can be potentially extrapolated to other riparian, clonal shrubs that share similar biology as S. virginiana.     

Genetic structure of Amazonian populations of Hevea brasiliensis is shaped by hydrographical network and isolation by distance

Current ex situ collections of the wild rubber tree Hevea brasiliensis were gathered more than 25 years ago mainly from populations in Southwest Amazonia. Information regarding the structure of genetic diversity among these populations is fragmented and needs refinement. Having analyzed polymorphisms at 15 microsatellites loci throughout the Hevea genome, we report on the genetic diversity of 307 clonally propagated individuals from 19 different collection points. All analyzed loci were highly polymorphic, averaging 21.7 alleles per locus, and displayed heterozygote deficiency. Differentiation among populations was analyzed on a subsample of 220 individuals from 14 populations. The overall F _ST is 0.12; although low, this value is significant and suggests moderate differentiation. Mato Grosso (Brazil) populations were genetically more distant from all other populations. A highly significant correlation between genetic and geographical distances over the whole area under study was shown by the Mantel test. The Bayesian algorithm used by STRUCTURE software allowed for the identification of three population clusters that match the boundaries of hydrographical basins of the main Amazon River tributaries. Differentiation among populations could therefore be jointly explained by both geographical location within the hydrographical Amazon network and by isolation by distance among populations belonging to distinct catchments. It is likely that some of these populations no longer exist due to severe deforestation. We provide suggestions to improve genetic resource management and sampling strategy in natural populations of H. brasiliensis.     

Genomic analysis of a cardinalfish with larval homing potential reveals genetic admixture in the Okinawa Islands

Discrepancies between potential and observed dispersal distances of reef fish indicate the need for a better understanding of the influence of larval behaviour on recruitment and dispersal. Population genetic studies can provide insight on the degree to which populations are connected, and the development of restriction site‐associated sequencing (RAD‐Seq) methods has made such studies of nonmodel organisms more accessible. We applied double‐digest RAD‐Seq methods to test for population differentiation in the coral reef‐dwelling cardinalfish, Siphamia tubifer, which based on behavioural studies, have the potential to use navigational cues to return to natal reefs. Analysis of 11,836 SNPs from fish collected at coral reefs in Okinawa, Japan, from eleven locations over 3 years reveals little genetic differentiation between groups of S. tubifer at spatial scales from 2 to 140 km and between years at one location: pairwise F_ST values were between 0.0116 and 0.0214. These results suggest that the Kuroshio Current largely influences larval dispersal in the region, and in contrast to expectations based on studies of other cardinalfishes, there is no evidence of population structure for S. tubifer at the spatial scales examined. However, analyses of outlier loci putatively under selection reveal patterns of temporal differentiation that indicate high population turnover and variable larval supply from divergent source populations between years. These findings highlight the need for more studies of fishes across various geographic regions that also examine temporal patterns of genetic differentiation to better understand the potential connections between early life‐history traits and connectivity of reef fish populations.     

Extremely low nucleotide polymorphism in Pinus krempfii Lecomte,a unique flat needle pine endemic to Vietnam

Pinus krempfii Lecomte is a morphologically and ecologically unique pine, endemic to Vietnam. It is regarded as vulnerable species with distribution limited to just two provinces: Khanh Hoa and Lam Dong. Although a few phylogenetic studies have included this species, almost nothing is known about its genetic features. In particular, there are no studies addressing the levels and patterns of genetic variation in natural populations of P. krempfii. In this study, we sampled 57 individuals from six natural populations of P. krempfii and analyzed their sequence variation in ten nuclear gene regions (approximately 9 kb) and 14 mitochondrial (mt) DNA regions (approximately 10 kb). We also analyzed variation at seven chloroplast (cp) microsatellite (SSR) loci. We found very low haplotype and nucleotide diversity at nuclear loci compared with other pine species. Furthermore, all investigated populations were monomorphic across all mitochondrial DNA (mtDNA) regions included in our study, which are polymorphic in other pine species. Population differentiation at nuclear loci was low (5.2%) but significant. However, structure analysis of nuclear loci did not detect genetically differentiated groups of populations. Approximate Bayesian computation (ABC) using nuclear sequence data and mismatch distribution analysis for cpSSR loci suggested recent expansion of the species. The implications of these findings for the management and conservation of P. krempfii genetic resources were discussed.     

The phylogeography of Fagus hayatae (Fagaceae): genetic isolation among populations

The beech species Fagus hayatae is an important relict tree species in subtropical China, whose biogeographical patterns may reflect floral responses to climate change in this region during the Quaternary. Previous studies have revealed phylogeography for three of the four Fagus species in China, but study on F. hayatae, the most sparsely distributed of these species, is still lacking. Here, molecular methods based on eight simple sequence repeat (SSR) loci of nuclear DNA (nDNA) and three chloroplast DNA (cpDNA) sequences were applied for analyses of genetic diversity and structure in 375 samples from 14 F. hayatae populations across its whole range. Both nDNA and cpDNA indicated a high level of genetic diversity in this species. Significant fixation indexes and departures from the Hardy–Weinberg equilibrium, with a genetic differentiation parameter of R_st of 0.233, were detected in nDNA SSR loci among populations, especially those on Taiwan Island, indicating strong geographic partitioning. The populations were classified into two clusters, without a prominent signal of isolation‐by‐distance. For the 15 haplotypes detected in the cpDNA sequence fragments, there was a high genetic differentiation parameter (G_st = 0.712) among populations. A high G_st of 0.829 was also detected outside but not within the Sichuan Basin. Consistent with other Fagus species in China, no recent population expansion was detected from tests of neutrality and mismatch distribution analysis. Overall, genetic isolation with limited gene flow was prominent for this species and significant phylogeographic structures existed across its range except for those inside the Sichuan Basin. Our study suggested long‐term geographic isolation in F. hayatae with limited population admixture and the existence of multiple refugia in the mountainous regions of the Sichuan Basin and southeast China during the Quaternary. These results may provide useful information critical for the conservation of F. hayatae and other Chinese beech species.     

Genetic hitchhiking in a subdivided population of <Emphasis Type="Italic">Mytilus edulis</Emphasis>

Background  

Few models of genetic hitchhiking in subdivided populations have been developed and the rarity of empirical examples is even more striking. We here provide evidences of genetic hitchhiking in a subdivided population of the marine mussel Mytilus edulis. In the Bay of Biscay (France), a patch of M. edulis populations happens to be separated from its North Sea conspecifics by a wide region occupied only by the sister species M. galloprovincialis. Although genetic differentiation between the two M. edulis regions is largely non-significant at ten marker loci (average F_ST~0.007), a strong genetic differentiation is observed at a single locus (F_ST = 0.25). We validated the outlier status of this locus, and analysed DNA sequence polymorphism in order to identify the nature of the selection responsible for the unusual differentiation.     

Ten new microsatellite loci for analysis of genetic diversity in isolated populations of the Alpine land snail <Emphasis Type="Italic">Cylindrus obtusus</Emphasis>

In this study we characterized 10 polymorphic microsatellite markers for the land snail Cylindrus obtusus, an endemism of the Austrian Alps with a distribution in isolated populations above approximately 1,600 m. The microsatellite loci were analyzed in 44 individuals from two populations. Number of alleles per locus ranged between two and eight. Observed heterozygosity ranged between 0.00 and 1.00, and expected heterozygosity between 0.09 and 0.72. No significant linkage disequilibrium was found between pairs of loci. One of the sampled populations (Dachstein) showed no deviation from Hardy–Weinberg equilibrium and no presence of null alleles, whereas the other one (Schneeberg) did. These diverging results probably reflect differences in population structure rather than characteristics of the microsatellite loci and underline the usefulness of these markers for studying genetic diversity, population structure and differentiation in C. obtusus.     

Genetic divergence in nuclear genomes between populations of <Emphasis Type="Italic">Fagus crenata</Emphasis> along the Japan Sea and Pacific sides of Japan

Genetic diversity and structure in Fagus crenata were studied by analyzing 14 nuclear microsatellite loci in 23 populations distributed throughout the species’ range. Although population differentiation was very low (F _ST = 0.027; R _ST = 0.041), both neighbor-joining tree and Bayesian clustering analyses provided clear evidence of genetic divergence between populations along the Japan Sea (Japan Sea lineage) and Pacific (Pacific lineage) sides of Japan, indicating that physical barriers to migration and gene flow, notably the mountain ranges separating the populations along the Japan Sea and Pacific sides, have promoted genetic divergence between these populations. The two lineages of the nuclear genome are generally consistent with those of the chloroplast genome detected in a previous study, with several discrepancies between the two genomes. Within-population genetic diversity was generally very high (average H _E = 0.839), but decreased in a clinal fashion from southwest to northeast, largely among populations of the Japan Sea lineage. This geographical gradient may have resulted from the late-glacial and postglacial recolonization to the northeast, which led to a loss of within-population genetic diversity due to cumulative founder effects.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Genetic diversity among alfalfa (Medicago sativa) cultivars coming from a breeding program, using SSR markers

Alfalfa (Medicago sativa) is an autotetraploid, allogamous and heterozygous species whose cultivars are synthetic populations. The breeders apply selection pressure for some agronomic traits within a breeding pool to increase the frequency of favorable individuals. The objective of this study was to investigate the differentiation level among seven cultivars originating from one breeding program, and between these cultivars and the breeding pool, with eight SSR markers. These highly polymorphic and codominant markers, together with recent population genetic statistics extended to autotetraploids, offer tools to analyse genetic diversity in alfalfa. The number of alleles per locus varied between 3 and 24. All loci were at a panmictic equilibrium in the cultivars, except one, probably because of null alleles. With seven SSR loci, each cultivar was at panmictic equilibrium. The mean gene diversity was high, ranging from 0.665 to 0.717 in the cultivars. The parameter F _ST indicated a low but significant diversity among cultivars. Among 21 pairs of cultivars, 15 were significantly different. The breeding pool also had a high diversity, and was significantly different from each cultivar except the most recent one. Considering the characteristics of the breeding program and the mode of cultivar elaboration, we found that they were unable to generate a large variety differentiation. Estimation of population genetics parameters at SSR loci can be applied for assessing the differences between cultivars or populations, either for variety distinction or the management of genetic resources.