Conservation genetics of Amorpha georgiana (Fabaceae), an endangered legume of the Southeastern United States

Amorpha georgiana (Fabaceae) is an endangered legume species found in longleaf pine savannas in the Southeastern United States. Approximately 900 individuals and 14 populations remain, most of which are concentrated in North Carolina. Eleven microsatellite loci were used to explore genetic diversity, population structure and recent population bottlenecks using genotypic data from 132 individuals collected at ten different localities. Although A. georgiana is quite rare, it exhibited high levels of genetic diversity (17.7 alleles/locus; H _o = 0.65, H _E = 0.75). Most of the genetic variation was found within rather than between populations of this species. The single remaining Georgia population was well differentiated from populations of the Carolinas ( F _ST > 0.1), which had weaker structure among them ( F _ST < 0.1). Only a geographically disjunct population showed strong evidence of a recent population bottleneck, perhaps due to a recent founder event. Hybridization with A. herbacea was also detected. For conservation management plans, A. georgiana populations in each geographic region (North Carolina, South Carolina and Georgia) plus a disjunct population in North Carolina (Holly Shelter) should be treated as separate management units for which in situ conservation, including habitat restoration and use of prescribed burns, should ensure persistence of this species and preservation of its evolutionary potential.     

Extremely low levels of genetic diversity in the terrestrial orchid Epipactis thunbergii (Orchidaceae) in South Korea: implications for conservation

Levels of allozyme variation, population genetic structure, and fine-scale genetic structure (FSGS) of the rare, both sexually and clonally reproducing terrestrial orchid Epipactis thunbergii were examined for eight ( N  = 734) populations in a 20 × 20-km area in South Korea. Twenty-three putative allozyme loci resolved from 15 enzyme systems were used. Extremely low levels of allozyme variation were found within populations: the mean frequency of polymorphic loci was 3.8% [isocitrate dehydrogenase ( Idh-2 ) with two alleles was polymorphic across populations], the mean number of alleles per locus was 1.04, and the mean expected heterozygosity was 0.013. The overall fixation index was not significantly different from zero ( F _IS = 0.069), although the species is self-compatible. However, a significantly high degree of population differentiation was found between populations at Idh-2 ( F _ST = 0.388) in the studied area. Furthermore, spatial autocorrelation analyses revealed a significant FSGS (up to 3 m) within populations. These observations suggest that the main explanatory factors for the extremely low levels of genetic diversity and the shaping of the population genetic structure of E. thunbergii are genetic drift as a result of a small effective population size, a restricted gene flow, and the isolation of populations. Considering the current genetic structure of E. thunbergii , three guidelines are suggested for the development of conservation strategies for the species in South Korea: (1) protection of habitats of standing populations; (2) prohibition by law of any collection of E. thunbergii ; and (3) protection of nearby pollinator populations, given the fact that fruit set in natural habitats is very low.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 161–169.     

Genetic variation and clonal diversity in island Cephalanthera rubra populations from the Biebrza National Park, Poland

Three populations of the orchid species Cephalanthera rubra localized on neighbouring mineral islands in the Biebrza National Park (north-east Poland) were examined using 16 allozyme loci. The percentage of polymorphic loci in these C. rubra populations ranged from 25 to 31.3% and the mean number of alleles per polymorphic locus from 1.25 to 1.31. The levels of observed heterozygosity at polymorphic loci were higher than expected for all populations. The overabundance of heterozygotes was noted in every population. The values of fixation indices were highly negative. The genetic differentiation among the three populations was small ( F _ST = 0.0173), albeit statistically significant ( P  < 0.001). The level of gene flow varied from 13 to 16. Fourteen multilocus genotypes were found among the 401 ramets sampled from the populations. Only four of them were common for all populations. The frequency of genotypes was different in these three populations of C. rubra . Some of the multilocus genotypes dominated in the given populations, others were found sporadically. Genetic variation within C. rubra populations is mainly the result of small population sizes and reflects the influences of breeding system and type of reproduction. On the other hand, these properties are also connected with habitat conditions.  © 2003 The Linnean Society of London, Botanical Journal of the Linnean Society , 2003, 143 , 99–108.     

Microsatellite DNA analysis of northern pike (Esox lucius L.) populations: insights into the genetic structure and demographic history of a genetically depauperate species

The northern pike Esox lucius L. is a freshwater fish exhibiting pronounced population subdivision and low genetic variability. However, there is limited knowledge on phylogeographical patterns within the species, and it is not known whether the low genetic variability reflects primarily current low effective population sizes or historical bottlenecks. We analysed six microsatellite loci in ten populations from Europe and North America. Genetic variation was low, with the average number of alleles within populations ranging from 2.3 to 4.0 per locus. Genetic differentiation among populations was high (overall θ_ST = 0.51; overall ρ_ST = 0.50). Multidimensional scaling analysis of genetic distances between populations and spatial analysis of molecular variance suggested a single phylogeographical race within the sampled populations from northern Europe, whereas North American and southern European populations were highly distinct. A population from Ireland was monomorphic at all loci, presumably reflecting founder events associated with introduction of the species to the island in the sixteenth century. Bayesian analysis of demographic parameters showed differences in θ (a product of effective population size and mutation rate) among populations from large and small water bodies, but the relative differences in θ were smaller than expected, which could reflect population subdivision within the larger water bodies. Finally, the analyses showed drastic population declines on a time scale of several thousand years within European populations, which we ascribe to either glacial bottlenecks or postglacial founder events.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 91–101.     

Two approaches for the selection of Relevant Genetic Units for Conservation in the narrow European endemic steppe plant Boleum asperum (Brassicaceae)

Several methods based on population biology, biogeography, ecology, and genetics have been traditionally used for the identification of units for conservation below the species level. We use a combination of two methods based on population genetic structure estimators and on probabilities of loss of rare alleles to identify the Relevant Genetic Units for Conservation (RGUCs). The aims were to assess the genetic diversity and population structure of the endemic steppe plant Boleum asperum (Brassicaceae), and to determine how many and which populations significantly represent the total genetic diversity and the rarest allelic variation. Despite the high amplified fragment length polymorphism genetic diversity values detected in B. asperum ( h _T = 0.744), caused probably by its hexaploidy and allogamy, moderate spatial genetic differentiation was detected among populations (< 20%) and geographical ranges (> 13%), suggesting the existence of an ancestral continuous distribution range that was fragmented into separate 'islands' in more recent historical times. Five RGUCs, accounting for the 85.10% of the total genetic variation and representative of the entire geographical distribution of the species, were selected for in situ conservation. Ex situ conservation is proposed to complement the preservation of B. asperum . This method of objective selection of populations may be applied to other candidate taxa for conservation with prior adjustment of the threshold values of diversity required for effective protection of each particular species.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 341–354.     

Contrasting levels of genetic diversity between the common, self-compatible Liparis kumokiri and rare, self-incompatible Liparis makinoana (Orchidaceae) in South Korea

Levels of allozyme variation and intrapopulation spatial genetic structure of the two terrestrial clonal orchids Liparis kumokiri , a self-compatible relatively common species, and L. makinoana , a self-incompatible rare species, were examined for 17 ( N  = 1875) and four ( N  = 425) populations, respectively, in South Korea. Populations of L. makinoana harboured high levels of genetic variation ( H _e = 0.319) across 15 loci. In contrast, L. kumokiri exhibited a complete lack of allozyme variation ( H _e = 0.000). Considering the lack of genetic variability, it is suggested that current populations of L. kumokiri in South Korea originated from a genetically depauperate ancestral population. For L. makinoana , a significant deficit of heterozygosity (mean F _IS = 0.198) was found in population samples excluding clonal ramets, suggesting that pollen dispersal is localized, generating biparental inbreeding. The significant fine-scale genetic structuring (≤ 2 m) found in a previous study, in addition to the moderate levels of population differentiation ( F _ST = 0.107) and the significant relationship between genetic and geographical distances ( r  = 0.680) found here, suggests a leptokurtic distribution of seed dispersal for L. makinoana . Although populations of L. makinoana harbour high levels of genetic variation, they are affected by a recent genetic bottleneck. This information suggests that genetic drift and limited gene flow could be the main evolutionary forces for speciation of a species-rich genus such as Liparis .  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 153 , 41–48.     

Low genetic diversity and allozymic evidence for autopolyploidy in the tetraploid Pyrenean endemic larkspur Delphinium montanum (Ranunculaceae)

Allozyme electrophoresis was conducted to survey the levels and distribution of genetic diversity in the tetraploid perennial larkspur Delphinium montanum (Ranunculaceae), which is endemic to the eastern Pyrenees of Spain and France and considered to be endangered. Seven populations were sampled, resolving 14 putative enzymatic loci belonging to eight enzymes. Banding patterns stained in gels revealed several enzymatic duplications attributable to autotetraploidy, such as the presence of both balanced and unbalanced heterozygotes and the lack of fixed heterozygosity. However, variability in D. montanum ( P  = 23.8%, A  = 1.48, and H _e = 0.082) was lower than that expected for an autotetraploid species. This, in addition to the scarcity of loci showing three or four different alleles, could indicate that this species has suffered genetic erosion through population bottlenecks, or, alternatively, that it is undergoing diploidization.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 211–222.     

Genetic evidence for population expansion in Hydrotaea irritans (Fallèn) (Diptera: Muscidae)

Allozyme variation in the sheep headfly Hydrotaea irritans was studied on two spatial scales. Geographic variation among seven Danish and one Dutch population revealed significant but rather low genetic differentiation with F _ST = 0.01 over all loci. The Dutch population was on average not more different from the Danish populations than the Danish populations from each other. Allele frequencies were very skewed with the most common allele always exceeding 0.85 and usually 0.9 in frequency, but with many rare alleles at some loci. Tests for neutrality of the variation at the nine polymorphic loci revealed highly significant deviations from expected homozygosity in this species, which was not found in a comparative analysis of allozyme variation at similar loci of seven other Hydrotaea species. To explain the peculiar observed pattern of allozyme variation in H. irritans , it is suggested that this species has successfully expanded its range and spread through northern and central Europe in the recent past. Alternatively, H. irritans may have recently invaded a new niche, resulting in increased abundance of the species and subsequent dispersal to former areas of the species distribution.     

Low population differentiation and high genetic diversity in the invasive species Carduus acanthoides L. (Asteraceae) within its native range in the Czech Republic

Colonizing species are predicted to suffer from reductions in genetic diversity during founding events. Although there is no unique mode of reproduction that is characteristic of successful plant colonizers, many of them are predominantly self-fertilizing or apomictic species, and almost all outcrossing colonizers are self-compatible. Carduus acanthoides comprises a species of disturbed habitats with wind-dispersed seeds that colonizes open spaces of various sizes. Population genetic diversity was expressed by assessing patterns of variation at nine putatively neutral allozyme loci within and among 20 natural populations in its native distribution range in the Czech Republic. Overall, C. acanthoides displayed high levels of genetic diversity compared to other herbaceous plants. The percentage of polymorphic loci was 84.5, with values of 2.37, 0.330, and 0.364 for the mean number of alleles per polymorphic locus ( A ), observed heterozygosity ( H _o), and expected heterozygosity ( H _e), respectively. There was only weak evidence of inbreeding within populations ( f  = 0.097) and very low genetic differentiation among populations ( θ  = 0.085). Analyses of the data provide strong evidence for isolation-by-distance for the whole study area. Even the colonizing species, C. acanthoides , currently supports a substantial amount of allozyme variation at both the species and population levels.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 596–607.     

Allozyme differentiation and genetic structure of the Lady's slipper (Cypripedium calceolus) island populations in north-east Poland

The population genetic structure in the rare and endangered perennial plant, Cypripedium calceolus , from north-east Poland was investigated. A total of 11 loci were analysed, of which five were polymorphic for the species and in every population. The 11 loci resolved gave a total of 19 alleles with the mean number of 1.73 alleles per locus. The proportion of polymorphic loci (P) was the same for the species and in every population, 45.5%. The average observed and expected heterozygosities for the total Lady's slipper data set were 0.156 and 0.184, respectively. The average expected heterozygosity for the populations under study ranged from 0.176 to 0.188 and did not differ significantly among populations. Genetic diversity within C. calceolus populations from the Biebrza Valley is relatively high compared with rare taxa and taxa with similar life histories. The genetic diversity among populations of C. calceolus studied ( F _ST = 0.016) is much smaller than the genetic diversity within populations. In the species studied, there was no significant correlation between genetic and geographical distances ( r = −0.73, Mantel test). This may indicate the relatively recent origin of the populations studied, or a high level of gene flow among populations. © 2002 The Linnean Society of London, Botanical Journal of the Linnean Society , 2002, 138 , 433–440.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Isozyme variation under different modes of reproduction in two clonal winter annuals, Sedum rosulato-bulbosum and Sedum bulbiferum (Crassulaceae)

We studied isozyme variation in two annual species that produce bulbils, Sedum rosulato-bulbosum , which includes both sexually reproducing plants and obligate clonal plants that result from triploidy (fertile and sterile S. rosulato-bulbosum , respectively), and an obligate clonal plant, Sedum bulbiferum , to examine the relationship between reproductive mode and isozyme variation. The sterile S. rosulato-bulbosum population exhibited no genotypic variation, but showed high genetic variation (gene diversity, H _e  = 0.60) because five of the six loci that we analyzed were heterozygous. Almost all ramets of S. bulbiferum across 20 populations shared an identical isozyme phenotype, although we could not identify the genetic basis of the phenotype. In contrast, fertile S. rosulato-bulbosum exhibited genotypic variation across the species, but comprised genotypically uniform and polymorphic populations whose genotypic variations correlated positively with the genetic variations within the populations ( H _e at the genet level per population ranged from 0.08 to 0.37). Genetic drift and habitat conditions inhibiting seedling recruitment may have caused this among-population variation. The results for sterile and fertile S. rosulato-bulbosum suggest that exclusive clonal reproduction causes low genotypic variation, but maintains genetic variation within individuals. Factors that affect the maintenance of genetic variation in these plants are discussed on the basis of these findings.     

Genetic Variation and Clonal Diversity of the Two Divergent Types of Clonal Populations of Leymus chinensis Tzvel on the Song Liao Steppe in the West of Northeastern China

The genetic variation and clonal diversity of two divergent types (grey-green and yellow-green) of clonal populations of Leymus chinensis Tzvel at 14 loci were compared. Total gene diversity (HT) and the coefficient of genetic differentiation (GST) were all higher for the yellow-green type (HT = 0.270; GST =0.186) than for the grey-green type (HT = 0.250; GST = 0.157) of L. chinensis. Rare alleles usually occurred as heterozygotes rather than homozygotes and significant deviations from Hardy-Weinberg equilibrium were found only at a few loci. This indicated that these two types of populations were mainly out-crossing. Clonal diversity, evenness of clones, and mean clone size were not significantly different between the two types. We found that differences between the clone size and genetic variation of the yellow-green type of populations occurred with different climate and habitat population groups. However, for the grey-green type of populations, these genetic variations decreased under conditions of different climate and habitat population groups.     

Genetics of the peripheral populations of the alpine bullhead, Cottus poecilopus (Scorpaeniformes, Cottidae) in Poland

The alpine bullhead ( Cottus poecilopus ) in Poland is highly polymorphic ( P for examined populations in the range of 0.0–0.45, mean P  = 0.11, P for species 0.64) in seven out of 11 loci studied ( Pgd-1, Mdh-1, Ck-1, Gpd-1, Gpd-2, Gpi-1, Gpi-2 ). However, they reveal very low heterozygosity (mean H _o = 0.007, mean H _e = 0.019), which may result from population subdivisions, small genetically effective size ( N _e) of some populations and frequent inbreeding, the latter being probably mediated through polygyny. The index of genetic differentiation ( G _ST) for the Polish metapopulation equals 0.48, which is relatively high. By far, the strongest are genetic differences between the Sudeten-Carpathian group and the Hańcza population. Gene flow between the Sudeten-Carpathian group and the Hańcza population is very limited and their divergence may have started as early as the last (Eemian) interglacial period.     

Cryptic speciation, genetic diversity and gene flow in the California turret spider Atypoides riversi (Araneae: Antrodiaetidae)

Populations of the turret spider Atypoides riversi from eight central Californian sites were compared based on variation at ten allozyme loci. Multidimensional scaling of interpopulation genetic distances defined four population units (Coast Range, Sierran, Valley, Jenness Camp), corroborating the distinctness of Coast Range and Sierran populations indicated by a prior study. While the species status of these units has yet to be determined, Jenness Camp is the most likely to represent a new species, given its clear genetic uniqueness (two fixed allelic differences). Populations in all units were generally in Hardy–Weinberg equilibrium with no evidence of inbreeding, though variability was minimal (mean H _o = 2.8%, mean P  = 15.4%). Reduced variability in these populations may be the result of repeated bottlenecks, environmental homogeneity, and/or directional selection. Interpopulation differentiation within units was significant in the absence of intervening forest habitat and was substantially less in its presence, indicating that gene flow is likely only when forest corridors exist. To foster preservation of the existing gene pool and enhance participation in it, management of the units of At. riversi should focus on maintaining as many populations in situ as possible and facilitating connections between them, while also creating or restoring habitat for potential colonization.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 27–37.     

Genetic diversity in the endangered dysploid larkspur Delphinium bolosii and its close diploid relatives in the series Fissa of the Western Mediterranean area

Allozyme diversity was evaluated in four closely-related taxa of the Delphinium series Fissa distributed throughout the Western Mediterranean area. All are considered threatened plants. Delphinium bolosii and Delphinium mansanetianum are narrowly endemic to the Eastern Iberian Peninsula, whereas Delphinium fissum ssp. sordidum is found in a few populations across the Peninsula. Delphinium fissum ssp. fissum is more widely distributed but often in small and isolated populations. In this group, Delphinium bolosii is dysploid (2 n  = 18) whereas the other taxa are diploid (2 n  = 16). A total of 12 populations were surveyed, including all known locations for D. bolosii , D. mansanetianum , and D. fissum ssp. sordidum . Eleven enzyme systems were assayed and 15 loci were resolved. Markedly depauperate values for genetic diversity were obtained for D. mansanetianum ( H _e = 0.013) and D. fissum ssp. sordidum ( H _e = 0.044). The estimates for D. fissum ssp. fissum ( H _e = 0.071) were below the values expected for widespread species. Small population size and marginal distribution have probably contributed to the low variability observed in this group. By contrast, D. bolosii exhibited comparatively larger populations and greater genetic diversity ( H _e = 0.138). We suggest that, apart from population size and local adaptation, genetic diversity during speciation may have been promoted by dysploidy through genomic recombination.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 773–784.     

Genetic structure and differentiation of Gentiana atuntsiensis W. W. Smith and G. striolata T. N. Ho (Gentianaceae) as revealed by ISSR markers

The Hengduan Mountains are the core region of the Himalaya hotspot, and are renowned for their high levels of endemism. Gentiana atuntsiensis and G. striolata are two closely related and morphologically similar species endemic to this region. In this study, inter-simple sequence repeat (ISSR) markers were used to investigate the genetic structure within these two species, as well as the differentiation between them. An analysis of molecular variance-derived estimate demonstrated only 13.5% of genetic differentiation between the two species. Considering their adjacent distribution patterns, low genetic divergence, and clear clustering into two groups, it is hypothesized that the two species arose from a rapid and recent speciation event induced mainly by geographical isolation. A relatively high level of genetic diversity was revealed in each species (Shannon's index of diversity: H _sp = 0.324 and 0.391; H _pop = 0.225 and 0.274; for G. atuntsiensis and G. striolata , respectively). Most of the genetic variation was partitioned within populations ( Ø _ST = 0.232 and 0.226 in G. atuntsiensis and G. striolata , respectively). The large population sizes, outcrossing breeding system, and small, light seeds that disperse widely may explain the genetic structure in both species.  © 2007 The Linnean Society of London. Botanical Journal of the Linnean Society , 2007, 154 , 225–232.     

Chloroplast DNA variation in the shrub Justicia areysiana (Acanthaceae) endemic to the monsoon affected coastal mountains of the southern Arabian Peninsula

The monsoon affected mountains of the southern Arabian Peninsula harbour in climatically favoured refugia vegetation elements of palaeo-African origin. To understand better the temporal and spatial differentiation of these refugia, chloroplast variation in Justicia areysiana Deflers (Acanthaceae), a shrub species endemic to the Yemeni and Omani mountains close to the Arabian Sea, was studied using PCR-RFLP and chloroplast microsatellite diversity. Eleven haplotypes were characterized and show a distinct geographical distribution pattern with a deep split between populations from south Yemeni fog oases and those from Hawf Mountains/Dhofar region in east Yemen and south Oman. Very limited haplotype diversity within populations (h_S = 0.15) and a high level of population differentiation (G_ST = 0.81) demonstrate the strong genetic isolation of populations from each other. Past oscillations between humid and arid periods connected with glacial and interglacial episodes in the Pleistocene and Holocene are considered responsible for the observed patterns of genetic variation.  © 2005 The Linnean Society of London, Botanical Journal of the Linnean Society , 2005, 148 , 437–444.     

Genetic diversity and structure of the cycad Zamia loddigesii Miq. (Zamiaceae): implications for evolution and conservation

The genetic diversity and structure of four populations of the cycad Zamia loddigesii were studied throughout its range in Mexico. Allozyme electrophoresis of 15 loci was conducted. The mean number of alleles per locus was 1.80 ± 0.09, the percentage of polymorphic loci was 66.6 ± 5.4, and the expected heterozygosity was 0.266 ± 0.02. The results indicated that the genetic diversity was relatively higher, with respect to tropical tree species and other cycads. The genetic variation explained by differences among populations was 18%. On average, gene flow between paired populations was similar ( Nm  = 1.6) to other tropical forest trees and cycad species. Our results indicated that the geographical isolation among populations of Z. loddigesii generated allele loss, as well as a clinal variation in the frequencies of two loci ( MDH and MNR2 ), in relation to the latitudinal distribution of populations. The populations have become fragmented due to increasingly higher pressure of habitat conversion and disturbance. The importance of the establishment of sanctuaries and protected areas and a reduction in deforestation is highlighted in this research as a way of preserving the high genetic diversity of this and other endemic species.  © 2006 The Linnean Society of London, Botanical Journal of the Linnean Society , 2006, 152 , 533–544.     

Population genetics of Dioon edule Lindl. (Zamiaceae, Cycadales): biogeographical and evolutionary implications

Dioon edule Lindl. (Zamiaceae) is a cycad endemic to Mexico, that occurs as one species D. edule and the geographical variety D. edule var. angustifolium (Miq.) Miq. Dioon edule has a north to south distribution in eastern Mexico. In this study, we analysed 14 allozymic loci in eight populations of D. edule from its total distribution range by sampling all known populations. Patterns of diversity and genetic variability, within and among populations, were obtained. The mean number of alleles per locus ( A ) was 1.44 and the percentage of polymorphic loci was relatively high ( P  = 54.78). The mean observed ( H _O ) and expected heterozygosity ( H _E ) were 0.27 and 0.24, respectively. F -statistics revealed an excess of heterozygous genotypes, locally and globally ( F  = −0.17 and f  = −0.27, respectively). The genetic variation explained by differences among populations was only 7.5%. We also detected a negative relationship between genetic diversity and latitude. On average, the gene flow between population pairs was relatively high ( Nm  = 2.98); furthermore, gene flow between population pairs was significantly correlated with geographical distances ( r  = −0.38, P  = 0.025). Therefore, patterns of genetic diversity in D. edule appear to be associated with the post-Pleistocene spread of the species, from its southerly (origin) to its northerly range (derived populations, including its central distribution). The biogeographical and evolutionary aspects of the results of this study are discussed. We recognize Dioon angustifolium Miq. for the northernmost disjunct populations.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80 , 457–467