Genetic spatial autocorrelation can readily detect sex-biased dispersal

Sex-biased dispersal is expected to generate differences in the fine-scale genetic structure of males and females. Therefore, spatial analyses of multilocus genotypes may offer a powerful approach for detecting sex-biased dispersal in natural populations. However, the effects of sex-biased dispersal on fine-scale genetic structure have not been explored. We used simulations and multilocus spatial autocorrelation analysis to investigate how sex-biased dispersal influences fine-scale genetic structure. We evaluated three statistical tests for detecting sex-biased dispersal: bootstrap confidence intervals about autocorrelation r values and recently developed heterogeneity tests at the distance class and whole correlogram levels. Even modest sex bias in dispersal resulted in significantly different fine-scale spatial autocorrelation patterns between the sexes. This was particularly evident when dispersal was strongly restricted in the less-dispersing sex (mean distance <200 m), when differences between the sexes were readily detected over short distances. All tests had high power to detect sex-biased dispersal with large sample sizes (n ≥ 250). However, there was variation in type I error rates among the tests, for which we offer specific recommendations. We found congruence between simulation predictions and empirical data from the agile antechinus, a species that exhibits male-biased dispersal, confirming the power of individual-based genetic analysis to provide insights into asymmetries in male and female dispersal. Our key recommendations for using multilocus spatial autocorrelation analyses to test for sex-biased dispersal are: (i) maximize sample size, not locus number; (ii) concentrate sampling within the scale of positive structure; (iii) evaluate several distance class sizes; (iv) use appropriate methods when combining data from multiple populations; (v) compare the appropriate groups of individuals.     

The influence of dispersal on macroecological patterns of Lesser Antillean birds

Aim Dispersal is often assumed to be a major force in shaping macroecological patterns, but this is rarely tested. Here I describe macroecological patterns for two groups of Lesser Antillean birds and then use population genetic data to assess if differences in dispersal ability could be responsible for the groups’ contrasting patterns. Importantly, the population genetic data are derived independently from any data used to generate the macroecological patterns. Location The Lesser Antilles, Caribbean. Methods I used data from the literature to construct species–area curves and evaluate the decline in species compositional similarity with geographic distance (hereafter distance–decay) for two sets of bird communities in the Lesser Antilles, those found in rain forest and those in dry forest. I then used mitochondrial DNA sequences from island populations to assess the dispersal ability of rain forest and dry forest species. Results Rain forest species show steeper species–area curves and greater distance–decay in community similarity than dry forest species, patterns that could be explained by rain forest species having more limited dispersal ability. Both conventional analyses of M, the number of migrants per generation between populations, and alternative analyses of D_A, the genetic distance between populations, suggest that rain forest species disperse between islands less frequently than dry forest species. Main conclusions Differences in dispersal ability are a plausible explanation for the contrasting macroecological patterns of rain forest and dry forest species. Additionally, historical factors, such as the taxon cycle and Pleistocene climate fluctuations, may have played a role in shaping the distribution patterns of Lesser Antillean birds.     

Microgeographic population genetic structure in the northern water snake, Nerodia sipedon sipedon detected using microsatellite DNA loci

We describe the isolation and genetic characterization of eight microsatellite DNA loci from the northern water snake, Nerodia sipedon sipedon and use these loci to analyse levels of genetic differentiation between local (< 2 km apart) populations of these snakes in Ontario. These loci are variable, with expected heterozygosities ranging from 0.28 to 0.91, and can correctly exclude nonsires in parentage analyses with a high probability (0.998). Population analyses reveal significant deviation from expected heterozygosity levels for one population, probably a result of a null allele(s) at a single locus and small but significant levels of genetic differentiation among all three populations. This demonstrates that microgeographic genetic structure exists in this species, possibly due to limited dispersal.     

Slight differences among individuals and the unified neutral theory of biodiversity

The unified neutral theory of biodiversity provides a very simple and counterintuitive explanation of species diversity patterns. By specifying speciation, community size and dispersal, and completely ignoring differences among individual organisms and species, it generates biodiversity patterns that remarkably resemble natural ones. Here I show that adding even slight differences among organisms generates very different patterns and predictions. In large communities with widespread dispersal, heritable differences in viability among individual organisms lead to biodiversity patterns characterised by the overdominance of a single species comprising organisms with relatively high fitness. In communities with local dispersal, the same differences produce rapid community extinction. I conclude that the unified neutral theory is not robust to slight deviations from its most controversial assumption.     

Is the population turnover of patchy‐distributed annuals determined by dormancy dynamics or dispersal processes?

In species with fragmented distribution, regional turnover dynamics is given by the processes of local population extinction and patch (re)colonization by migrants spreading from neighboring occupied patches. In plants with dormant stages (e.g. seeds) and limited dispersal capacity, regional dynamics based on dispersal processes can be overridden by pseudo-turnover determined by signals inducing or breaking dormancy (e.g. due to changes in habitat quality) resulting in a low importance of habitat configuration and size.
In this study, I investigated the turnover dynamics of 5 annual plant species growing on ant mounds of Lasius flavus over three years. I analyzed whether the grassland-scale dynamics of these annuals is influenced by dispersal processes, or alternatively, by pseudo-turnover of soil seed populations. For that purpose I 1) searched for populations formed from soil seeds only, 2) compared the relative contribution of the soil seed bank and seed rain for population restoration after disappearance from the vegetation and 3) investigated whether colonization and extinction events are affected by patch isolation. I assumed if population turnover was rather a result of the soil seed bank dynamics then spatial effects would be hard to detect.
In spite of the presence of populations formed from soil seed and the relatively more important soil seed bank for potential population reestablishment, turnover dynamics followed the predictions of metapopulation theory. Population appearance was more probable in larger and less isolated patches. Probability of disappearance increased with decrease of population size that was negatively influenced by the patch size and its isolation. These findings indicate dispersal processes to be important in the turnover dynamics and only limited contribution of soil seed populations. Their small effectiveness is probably related to the low chance of recurrent disturbance on the mound surface.     

Seascape genetics: a coupled oceanographic-genetic model predicts population structure of Caribbean corals

Population genetics is a powerful tool for measuring important larval connections between marine populations [1-4]. Similarly, oceanographic models based on environmental data can simulate particle movements in ocean currents and make quantitative estimates of larval connections between populations possible [5-9]. However, these two powerful approaches have remained disconnected because no general models currently provide a means of directly comparing dispersal predictions with empirical genetic data (except, see [10]). In addition, previous genetic models have considered relatively simple dispersal scenarios that are often unrealistic for marine larvae [11-15], and recent landscape genetic models have yet to be applied in a marine context [16-20]. We have developed a genetic model that uses connectivity estimates from oceanographic models to predict genetic patterns resulting from larval dispersal in a Caribbean coral. We then compare the predictions to empirical data for threatened staghorn corals. Our coupled oceanographic-genetic model predicts many of the patterns observed in this and other empirical datasets; such patterns include the isolation of the Bahamas and an east-west divergence near Puerto Rico [3, 21-23]. This new approach provides both a valuable tool for predicting genetic structure in marine populations and a means of explicitly testing these predictions with empirical data.     

Larval transport modeling of deep-sea invertebrates can aid the search for undiscovered populations

Background

Many deep-sea benthic animals occur in patchy distributions separated by thousands of kilometres, yet because deep-sea habitats are remote, little is known about their larval dispersal. Our novel method simulates dispersal by combining data from the Argo array of autonomous oceanographic probes, deep-sea ecological surveys, and comparative invertebrate physiology. The predicted particle tracks allow quantitative, testable predictions about the dispersal of benthic invertebrate larvae in the south-west Pacific.

Principal Findings

In a test case presented here, using non-feeding, non-swimming (lecithotrophic trochophore) larvae of polyplacophoran molluscs (chitons), we show that the likely dispersal pathways in a single generation are significantly shorter than the distances between the three known population centres in our study region. The large-scale density of chiton populations throughout our study region is potentially much greater than present survey data suggest, with intermediate ‘stepping stone’ populations yet to be discovered.

Conclusions/Significance

We present a new method that is broadly applicable to studies of the dispersal of deep-sea organisms. This test case demonstrates the power and potential applications of our new method, in generating quantitative, testable hypotheses at multiple levels to solve the mismatch between observed and expected distributions: probabilistic predictions of locations of intermediate populations, potential alternative dispersal mechanisms, and expected population genetic structure. The global Argo data have never previously been used to address benthic biology, and our method can be applied to any non-swimming larvae of the deep-sea, giving information upon dispersal corridors and population densities in habitats that remain intrinsically difficult to assess.     

Female-biased dispersal, low female recruitment, unpaired males, and the extinction of small and isolated bird populations

Small and isolated populations are usually assumed to be at a high risk of extinction due to environmental or demographic stochasticity, genetic problems, or too little immigration. In birds, natal dispersal is usually female-biased, but the consequences of such a pattern on vulnerability to extinction of isolated populations has not received much attention before. In this paper I derive predictions as to how female-biased natal dispersal may differentially affect the extinction risk of populations and species with contrasting distributions, migratory behaviours, life histories and mating systems. Female-biased dispersal will lead to male-biased sex ratios in small, isolated or fragmented populations, in particular because recent research has shown that females often have a limited ability to search for mates and may therefore effectively be lost from the breeding population if they disperse into areas empty of conspecifics. I reviewed published studies on birds and found that a high proportion of unpaired males is common in isolated populations or populations in small habitat fragments. Dispersal of females may therefore increase the vulnerability to extinction of small or isolated populations, or populations at the periphery of a species' distribution range. I also predict that vulnerability to extinction should be greater for migratory than for resident species and greater for short-lived than for long-lived species because of differences in the time available for females to locate unpaired males. Further, extinction risk may also be greater for birds than for mammals due to differences in which sex disperses and patterns of parental care. Finally, mating system will also affect vulnerability to extinction when natal dispersal leads to biased sex ratios. I review available evidence for these predictions (e.g. songbird declines in North America) and discuss implications for conservation.     

Dispersal,kin competition,and the ideal free distribution in a spatially heterogeneous population

In this paper, we reanalyze simple models of the evolution of dispersal in a heterogeneous landscape. Previous analyses concluded that without temporal variability, dispersal can evolve only if it is not costly and if it is conditional on the habitat. If both conditions hold, these models predict that selection on dispersal should lead to balanced dispersal between habitats (the number of immigrants equals the number of emigrants in each habitat). To evaluate the generality of these conclusions, we extended the analysis of these models to finite populations. This requires us to establish fitness measures for finite class-structured populations. These fitness measures allow us to take kin competition into account. Our analysis shows that even without temporal variability, conditional dispersal and the absence of a dispersal cost are not necessary conditions for dispersal to evolve. In the absence of a dispersal cost, we predict that selection on conditional dispersal will always lead to panmixia and not simply to balanced dispersal. When dispersal is costly, we show that the ideal free distribution (IFD) and balanced dispersal do not occur. Our results show that the deviations from IFD are of the order of the dispersal cost. We propose an approach to test our predictions.     

The evolution of growth rates on an expanding range edge

Individuals in the vanguard of a species invasion face altered selective conditions when compared with conspecifics behind the invasion front. Assortment by dispersal ability on the expanding front, for example, drives the evolution of increased dispersal, which, in turn, leads to accelerated rates of invasion. Here I propose an additional evolutionary mechanism to explain accelerating invasions: shifts in population growth rate (r). Because individuals in the vanguard face lower population density than those in established populations, they should (relative to individuals in established populations) experience greater r-selection. To test this possibility, I used the ongoing invasion of cane toads (Bufo marinus) across northern Australia. Life-history theory shows that the most efficient way to increase the rate of population growth is to reproduce earlier. Thus, I predict that toads on the invasion front will exhibit faster individual growth rates (and thus will reach breeding size earlier) than those from older populations. Using a common garden design, I show that this is indeed the case: both tadpoles and juvenile toads from frontal populations grow around 30 per cent faster than those from older, long established populations. These results support theoretical predictions that r increases during range advance and highlight the importance of understanding the evolution of life history during range advance.     

Metapopulation dynamics in an aphid-parasitoid system

Metapopulation theory makes a number of predictions concerning the effects of dispersal on the persistence of predator-prey or host-parasitoid systems. While the stabilising effects of dispersal have been shown in a number of laboratory studies, evidence from field studies remains scarce due to a lack of suitable model systems. I describe a host-parasitoid system that shows a classical metapopulation structure with frequent extinctions and colonisations consisting of the aphidiid Lysiphlebus hirticornisand the aphid Metopeurum fuscoviride. Both the parasitoid and the aphid are specialists on their respective hosts. I followed the dynamics of host and parasitoid on individually marked tansy (Tanacetum vulgare) plants, the host of M. fuscoviride. Dynamics of host and parasitoid populations were characterized by frequent extinctions and colonisations. Mean longevity of aphid colonies was only 3.1 weeks. Parasitism by L. hirticorniswas a main cause of extinction for the aphid as rates of parasitism often reached 100%, in particular towards the end of the field season. Patchiness in this system occurs at two spatial scales. Aphid colonies form on single tansy ramets = shoots but movements of aphid individuals among ramets within a particular tansy genet are frequent. Because aphids can persist on a genet for a large numer of generations, it is argued that local populations form on genets rather than ramets. The number of host and parasitoid extinctions described in this study exceeds the number of extinctions usually observed in field studies of host-parasitoid metapopulations. It is suggested that aphid-parasitoid systems such as the one studied in this paper may be good models to test the predictions of metapopulation theory.     

Panbiogeography from tracks to ocean basins: evolving perspectives

Misconceptions arising from efforts to translate panbiogeography into terms used in other biogeographic and evolutionary theories are discussed with respect to 5 ) critique of panbiogeography. Croizat’s rejection of ‘Darwinian dispersal’ applies only to efforts to utilize this concept as a general explanation for biogeographic patterns. The conceptual difference between distribution and panbiogeographic dispersal maps is illustrated to show that Croizat did not synonymize distribution and dispersal. Croizat’s position on continental drift and plate tectonics does not support 5 claim that Croizat ‘for a long time’ refused to accept the theory of plate tectonics. The methodological relationship between panbiogeographic analysis and geology suggests an independence of methodology that prevents geological theory from falsifying panbiogeographic predictions. Panbiogeographic predictions for the eastern Pacific are shown to be in agreement with current historical geological models. Claims by 5 that the panbiogeographic method is variable and questionable are evaluated with respect to the biogeographic homology of primitive frogs, ratite birds, and southern beeches to demonstrate the consistent application of minimal distance, main massing, phylogenetic affinity and baseline criteria. Panbiogeographic classification concepts are contrasted with the Darwinian system (supported by Cox) utilizing a concept of unitary geographical area based on the language of Roman military rule. Inconsistent positions expressed in recent critiques of panbiogeography may indicate an underlying and implicit acceptance of the empirical and theoretical progress generated by panbiogeography within modern biogeography. ‘The formation of groups has an invigorating effect in all spheres of human striving, perhaps mostly due to the struggle between the convictions and aims represented by the different groups’ (Einstein, 1938. Collier’s, 26 November).     

Integrating the effects of area, isolation, and habitat heterogeneity on species diversity: a unification of island biogeography and niche theory

We present an analytical model that unifies two of the most influential theories in community ecology, namely, island biogeography and niche theory. Our model captures the main elements of both theories by incorporating the combined effects of area, isolation, stochastic colonization and extinction processes, habitat heterogeneity, and niche partitioning in a unified, demographically based framework. While classical niche theory predicts a positive relationship between species richness and habitat heterogeneity, our unified model demonstrates that area limitation and dispersal limitation (the main elements of island biogeography) may create unimodal and even negative relationships between species richness and habitat heterogeneity. We attribute this finding to the fact that increasing heterogeneity increases the potential number of species that may exist in a given area (as predicted by niche theory) but simultaneously reduces the amount of suitable area available for each species and, thus, increases the likelihood of stochastic extinction. Area limitation, dispersal limitation, and low reproduction rates intensify the latter effect by increasing the likelihood of stochastic extinction. These analytical results demonstrate that the integration of island biogeography and niche theory provides new insights about the mechanisms that regulate the diversity of ecological communities and generates unexpected predictions that could not be attained from any single theory.     

Satyrinae butterflies from Sardinia and Corsica show a kaleidoscopic intraspecific biogeography (Lepidoptera,Nymphlidae)

The Mediterranean islands of Sardinia and Corsica are known for their multitude of endemics. Butterflies in particular have received much attention. However, no comprehensive studies aiming to compare populations of butterflies from Sardinia and Corsica with those from the neighbouring mainland and Sicily have been carried out. In the present study, the eleven Satyrinae species inhabiting Sardinia and Corsica islands were examined and compared with continental and Sicilian populations by means of geometric morphometrics of male genitalia. Relative warp computation, discriminant analyses, hierarchical clustering, and cross‐validation tests were used to identify coherent distributional patterns including both islands and mainland populations. The eleven species showed multifaceted distributional patterns, although three main conclusions can be drawn: (1) populations from North Africa and Spain are generally different from those belonging to the Italian Peninsula; (2) populations from Sardinia and Sicily often resemble the North Africa/Spain ones; Corsica shows transitional populations similar to those from France; and (3) sea barriers represent filters to dispersal, although their efficacy appears to be unrelated to their extension. Indeed, the short sea straits between Sardinia and Corsica and between Sicily and the Italian Peninsula revealed a strong effectiveness with respect to preventing faunal exchanges; populations giving onto sea channels between Corsica and Northern Italy and between Sicily and Tunisia showed a higher similarity. A comparison of island and mainland distributions of the eleven taxa have helped to unravel the complex co‐occurrence of historical factors, refugial dynamics, and recent (post‐glacial) dispersal with respect to shaping the populations of Mediterranean island butterflies. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 195–212.     

Characteristics of polymorphism at a VNTR locus 3' to the apolipoprotein B gene in five human populations.

We have analyzed the allele frequency distribution at the hypervariable locus 3' to the apolipoprotein B gene (ApoB 3' VNTR) in five well-defined human populations (Kacharis of northeast India, New Guinea Highlanders of Papua New Guinea, Dogrib Indians of Canada, Pehuenche Indians of Chile, and a relatively homogeneous Caucasian population of northern German extraction) by using the PCR technique. A total of 12 segregating alleles were detected in the pooled sample of 319 individuals. A fairly consistent bimodal pattern of allele frequency distribution, apparent in most of these geographically and genetically diverse populations, suggests that the ApoB 3' VNTR polymorphism predates the geographic dispersal of ancestral human populations. In spite of the observed high degree of polymorphism at this locus (expected heterozygosity levels 55%-78%), the genotype distributions in all populations (irrespective of their tribal or cosmopolitan nature) conform to their respective Hardy-Weinberg predictions. Furthermore, analysis of the congruence between expected heterozygosity and the observed number of alleles reveals that, in general, the allele frequency distributions at this locus are in agreement with the predictions of the classical mutation-drift models. The data also show that alleles that are shared by all populations have the highest average frequency within populations. These findings demonstrate the potential utility of highly informative hypervariable loci such as the ApoB 3' VNTR locus in population genetic research, as well as in forensic medicine and determination of biological relatedness of individuals.     

Population structure of the yellow-footed rock-wallaby Petrogale xanthopus (Gray, 1854) inferred from mtDNA sequences and microsatellite loci

The yellow-footed rock-wallaby Petrogale xanthopus is considered to be potentially vulnerable to extinction. This wallaby inhabits naturally disjunct rocky outcrops which could restrict dispersal between populations, but the extent to which that occurs is unknown. Genetic differences between populations were assessed using mitochondrial DNA (control region) sequencing and analysis of variation at four microsatellite loci among three geographically close sites in south-west Queensland (P. x. celeris) and, for mtDNA only, samples from South Australia (P. x. xanthopus) as well. Populations from South Australia and Queensland had phylogenetically distinct mtDNA, supporting the present classification of these two groups as evolutionarily distinct entities. Within Queensland, populations separated by 70 km of unsuitable habitat differed significantly for mtDNA and at microsatellite loci. Populations separated by 10 km of apparently suitable habitat had statistically homogeneous mtDNA, but a significant difference in allele frequency at one microsatellite locus. Tests for Hardy-Weinberg equilibrium and micro-geographical variation at microsatellite loci did not detect any substructuring between two wallaby aggregations within a colony encircling a single rock outcrop. Although the present study was limited by small sample sizes at two of the three Queensland locations examined, the genetic results suggest that dispersal between colonies is limited, consistent with an ecological study of dispersal at one of the sites. Considering both the genetic and ecological data, we suggest that management of yellow-footed rock-wallabies should treat each colony as an independent unit and that conservation of the Queensland and South Australian populations as separate entities is warranted.     

Scale‐dependent role of demography and dispersal on the distribution of populations in heterogeneous landscapes

Both dispersal and local demographic processes determine a population's distribution among habitats of varying quality, yet most theory, experiments, and field studies have focused on the former. We use a generic model to show how both processes contribute to a population's distribution, and how the relative importance of each mechanism depends on scale. In contrast to studies only considering habitat‐dependent dispersal, we show that predictions of ideal free distribution (IFD) theory are relevant even at landscape scales, where the assumptions of IFD theory are violated. This is because scales that inhibit one process, promote the other's ability to drive populations to the IFD. Furthermore, because multiple processes can generate IFDs, the pattern alone does not specify a causal mechanism. This is important because populations with IFDs generated by dispersal or demography respond much differently to shifts in resource distributions.     

Population size and fragmentation thresholds for the maintenance of genetic diversity in the herbaceous endemic Scutellaria montana (Lamiaceae)

Abstract.-The level and distribution of genetic variation is thought to be affected primarily by the size of individual populations and by gene flow among populations. Although the effects of population size have frequently been examined, the contributions of regional gene flow to levels of genetic variation are less well known. Here I examine the effects of population size and the number of neighboring populations (metapopulation density) on the distribution and maintenance of genetic diversity in an endemic herbaceous perennial. Reductions in the proportion of polymorphic loci and the effective number of alleles per locus were apparent for many populations with a census size of less than 100 individuals, but no effects of population size on levels of inbreeding were detected. I assess the effects of regional population density on levels of diversity and inbreeding using stepwise regression analysis of metapopulation diameter (i.e., the size of a circle within which population density is estimated). This procedure provides a spatially explicit evaluation of the effects of metapopulation size on population genetic parameters and indicates the critical number of neighboring populations (fragmentation threshold) for the regional maintenance of genetic diversity. Stepwise regression analyses revealed fragmentation thresholds at two levels; at a scale of 2 km, where small metapopulations resulted in greater levels of selfing or sibling mating, and at a scale of 8 km, where metapopulation size was positively associated with higher levels of genetic diversity. I hypothesize that the smaller fragmentation threshold may reflect higher levels of selfing in isolated populations because of the absence of pollinators. The larger threshold probably indicates the maximum distance over which pollen dispersal rates are high enough to counteract genetic drift. This study demonstrates that the regional distribution of populations can be an important factor for the long-term maintenance of genetic variation.     

Deeply divergent lineages of the widespread New Zealand amphipod Paracalliope fluviatilis revealed using allozyme and mitochondrial DNA analyses

1. We evaluated the population genetic structure of the common New Zealand amphipod Paracalliope fluviatilis using eight allozyme loci, and the mitochondrial cytochrome oxidase c subunit I (COI) gene locus. Morphological analyses were also conducted to evaluate any phenotypic differences. Individuals belonging to P. fluviatilis were collected from a total of 14 freshwater fluvial habitats on the North and South Islands, New Zealand. 2. We found evidence for strong genetic differentiation among locations (Wright's F_ST > 0.25), and fixed differences (non‐shared alleles) at two of the eight allozyme loci indicating the possibility of previously unknown species. Analysis of a 545‐bp fragment of the COI locus was mostly congruent with the allozyme data and revealed the same deeply divergent lineages (sequence divergences up to 26%). 3. Clear genetic breaks were identified between North Island and South Island populations. North Island populations separated by <100 km also showed genetic differences between east and west draining watersheds (sequence divergence >12%). Accordingly, present‐day dispersal among hydrologically isolated habitats appears minimal for this taxon. 4. Although population differences were clearly shown by allozyme and mtDNA analyses, individuals were morphologically indistinguishable. This suggests that, as in North American and European taxa (e.g. Hyalella and Gammarus), morphological conservatism may be prevalent among New Zealand's freshwater amphipods. We conclude that molecular techniques, particularly the COI gene locus, may be powerful tools for resolving species that show no distinctive morphological differences.     

Microevolution of S-allele frequencies in wild cherry populations: respective impacts of negative frequency dependent selection and genetic drift

Negative frequency dependent selection (NFDS) is supposed to be the main force controlling allele evolution at the gametophytic self-incompatibility locus (S-locus) in strictly outcrossing species. Genetic drift also influences S-allele evolution. In perennial sessile organisms, evolution of allelic frequencies over two generations is mainly shaped by individual fecundities and spatial processes. Using wild cherry populations between two successive generations, we tested whether S-alleles evolved following NFDS qualitative and quantitative predictions. We showed that allelic variation was negatively correlated with parental allelic frequency as expected under NFDS. However, NFDS predictions in finite population failed to predict more than half S-allele quantitative evolution. We developed a spatially explicit mating model that included the S-locus. We studied the effects of self-incompatibility and local drift within populations due to pollen dispersal in spatially distributed individuals, and variation in male fecundity on male mating success and allelic frequency evolution. Male mating success was negatively related to male allelic frequency as expected under NFDS. Spatial genetic structure combined with self-incompatibility resulted in higher effective pollen dispersal. Limited pollen dispersal in structured distributions of individuals and genotypes and unequal pollen production significantly contributed to S-allele frequency evolution by creating local drift effects strong enough to counteract the NFDS effect on some alleles.