Contemporary pollen and seed dispersal in a Prunus mahaleb population: patterns in distance and direction

Pollination and seed dispersal determine the spatial pattern of gene flow in plant populations and, for those species relying on pollinators and frugivores as dispersal vectors, animal activity plays a key role in determining this spatial pattern. For these plant species, reported dispersal patterns are dominated by short-distance movements with a significant amount of immigration. However, the contribution of seed and pollen to the overall contemporary gene immigration is still poorly documented for most plant populations. In this study we investigated pollination and seed dispersal at two spatial scales in a local population of Prunus mahaleb (L.), a species pollinated by insects and dispersed by frugivorous vertebrates. First, we dissected the relative contribution of pollen and seed dispersal to gene immigration from other parts of the metapopulation. We found high levels of gene immigration (18.50%), due to frequent long distance seed dispersal events. Second, we assessed the distance and directionality for pollen and seed dispersal events within the local population. Pollen and seed movement patterns were non-random, with skewed distance distributions: pollen tended moved up to 548 m along an axis approaching the N-S direction, and seeds were dispersed up to 990 m, frequently along the SW and SE axes. Animal-mediated dispersal contributed significantly towards gene immigration into the local population and had a markedly nonrandom pattern within the local population. Our data suggest that animals can impose distinct spatial signatures in contemporary gene flow, with the potential to induce significant genetic structure at a local level.     

ECOLOGICAL AND GENETIC CONSEQUENCES OF POLLINATION BY SEXUAL DECEPTION IN THE ORCHID CALADENIA TENTACTULATA

Only orchids affect pollination by the deceptive sexual attraction of male insects, a syndrome particularly well developed in Australia. We examined the ecological and genetic consequences of exclusive pollination by sexually attracted male thynnine wasps in the orchid Caladenia tentaculata. Male wasps respond rapidly to flowers artificially presented in 1 × 1 m² experimental patches. Sixty of 287 wasps approached within centimeters of the flower, but did not land. Of the remaining 79% who made floral contact, only 7.5% attempted copulation, the step critical for pollination. Wasps only rarely moved among patches (19% of flights) and none attempted copulation a second time, resembling observations in natural populations. We confirmed outcrossing and long distance pollen flow by monitoring how colored pollen moved in natural populations. Pollen movements approximated a linear rather than a leptokurtic distribution (mean distance: 17 m; maximum: 58 m). Pollinator visits varied independently of flower density in three of four populations with most solitary flowers being visited. Allozyme analysis revealed within-population fixation indices (F) close to zero and low levels of differentiation (FST) among populations. Despite behavioral evidence for long distance pollen flow, significant local genetic structure exists, perhaps reflecting restricted seed dispersal. Long distance pollen flow in C. tentaculata may therefore promote outbreeding by minimizing pollen transfers among related neighbors. Although this species is self-compatible, outcrossed progeny develop significantly faster than selfed progeny. Effective pollination at low flower densities could accentuate this advantage. The data are consistent with the predictions that deceptive pollination will result in long distance pollen flow, which may be of selective advantage at low density. Comparative studies of how food reward, food deceptive, and sexual deceptive pollination systems vary within a phylogenetic framework could further illuminate the evolution of sexual deception.     

Long-term persistence of species and the SLOSS problem

The single large or several small (SLOSS) problem has been addressed in a large number of empirical and theoretical studies, but no coherent conclusion has yet been reached. Here I study the SLOSS problem in the context of metapopulation dynamics. I assume that there is a fixed total amount A(0) of habitat available, and I derive formulas for the optimal number n and area A of habitat patches, where n=A(0)/A. I consider optimality in two ways. First, I attempt to maximize the time to metapopulation extinction, which is a relevant measure for metapopulation viability for rare and threatened species. Second, I attempt to maximize the metapopulation capacity of the habitat patch network, which corresponds both with maximizing the distance to the deterministic extinction threshold and with maximizing the fraction of occupied patches. I show that in the typical case, a small number of large patches maximizes the metapopulation capacity, while an intermediate number of habitat patches maximizes the time to extinction. The main conclusion stemming from the analysis is that the optimal number of patches is largely affected by the relationship between habitat patch area and rates of immigration, emigration and local extinction. Here this relationship is summarized by a single factor zeta, termed the patch area scaling factor.     

Migration of the clouded Apollo butterfly Parnassius mnemosyne in a network of suitable habitats – effects of patch characteristics

The clouded Apollo Parnassius mnemosyne is a food plant specialist with short but frequent movements between habitat patches. The short average dispersal distances suggest that the probability of colonisation of vacant patches decreases rapidly as the distance between the source and target patches increases, which means that a dense habitat network is needed for the conservation of the species. Both emigration rate and the number of immigrants varied among patches and were not affected only by isolation but also by several other patch characteristics. The model that explained most of the variation in emigration rates among patches included patch area and the number of conspecifics. The area and the population density of the target patch had significant effects on the number of arriving immigrants. Thus, the colonisation of vacant patches is dependent on these patch characteristics. Generally, emigration rates were lower and residence times longer in large patches with many conspecifics. Butterfly density was the most important single factor explaining the variation in the number of immigrants among patches, although the positive effect of the area of the target patch was also significant. As a consequence of the marked positive density dependence caused by conspecific attraction, small patches with higher than average butterfly density, receive more immigrants than could be expected based on the patch area only. Due to conspecific attraction, per capita immigration rates are higher in small than large patches. Thus, immigration may have a more significant effect on the local dynamics of small than large populations.     

Area-dependent migration by ringlet butterflies generates a mixture of patchy population and metapopulation attributes

Interpretation of spatially structured population systems is critically dependent on levels of migration between habitat patches. If there is considerable movement, with each individual visiting several patches, there is one ”patchy population”; if there is intermediate movement, with most individuals staying within their natal patch, there is a metapopulation; and if (virtually) no movement occurs, then the populations are separate (Harrison 1991, 1994). These population types actually represent points along a continuum of much to no mobility in relation to patch structure. Therefore, interpretation of the effects of spatial structure on the dynamics of a population system must be accompanied by information on mobility. We use empirical data on movements by ringlet butterflies, Aphantopus hyperantus, to investigate two key issues that need to be resolved in spatially-structured population systems. First, do local habitat patches contain largely independent local populations (the unit of a metapopulation), or merely aggregations of adult butterflies (as in patchy populations)? Second, what are the effects of patch area on migration in and out of the patches, since patch area varies considerably within most real population systems, and because human landscape modification usually results in changes in habitat patch sizes? Mark-release-recapture (MRR) data from two spatially structured study systems showed that 63% and 79% of recaptures remained in the same patch, and thus it seems reasonable to call both systems metapopulations, with some capacity for separate local dynamics to take place in different local patches. Per capita immigration and emigration rates declined with increasing patch area, while the resident fraction increased. Actual numbers of emigrants either stayed the same or increased with area. The effect of patch area on movement of individuals in the system are exactly what we would have expected if A. hyperantus were responding to habitat geometry. Large patches acted as local populations (metapopulation units) and small patches simply as locations with aggregations (units of patchy populations), all within 0.5 km². Perhaps not unusually, our study system appears to contain a mixture of metapopulation and patchy-population attributes.     

Relative impact of mate versus pollinator availability on pollen limitation and outcrossing rates in a mass‐flowering species

Plant mating systems are driven by several pre‐pollination factors, including pollinator availability, mate availability and reproductive traits. We investigated the relative contributions of these factors to pollination and to realized outcrossing rates in the patchily distributed mass‐flowering shrub Rhododendron ferrugineum. We jointly monitored pollen limitation (comparing seed set from intact and pollen‐supplemented flowers), reproductive traits (herkogamy, flower size and autofertility) and mating patterns (progeny array analysis) in 28 natural patches varying in the level of pollinator availability (flower visitation rates) and of mate availability (patch floral display estimated as the total number of inflorescences per patch). Our results showed that patch floral display was the strongest determinant of pollination and of the realized outcrossing rates in this mass‐flowering species. We found an increase in pollen limitation and in outcrossing rates with increasing patch floral display. Reproductive traits were not significantly related to patch floral display, while autofertility was negatively correlated to outcrossing rates. These findings suggest that mate limitation, arising from high flower visitation rates in small plant patches, resulted in low pollen limitation and high selfing rates, while pollinator limitation, arising from low flower visitation rates in large plant patches, resulted in higher pollen limitation and outcrossing rates. Pollinator‐mediated selfing and geitonogamy likely alleviates pollen limitation in the case of reduced mate availability, while reduced pollinator availability (intraspecific competition for pollinator services) may result in the maintenance of high outcrossing rates despite reduced seed production.     

Dispersal-induced desynchronization: from metapopulations to metacommunities

The conceptualization of fragmented populations in terms of metapopulation theory has become standard over the last three decades. It is well known that increases in between‐patch migration rates cause more synchronous population fluctuations and that this coherence increases the risk of global metapopulation extinction. Because species’ migration rates and the probability of individuals surviving migration events depend on the effective distance between patches, the benefit of improving conservation corridors or the matrix between habitat patches has been questioned. As populations occur in the context of larger communities, moving from a metapopulation to a metacommunity model framework is a natural extension to address the generality of these conclusions. We show how considering a metacommunity can modify the conclusion that decreasing the effective distance between habitat patches (via improving matrix quality or other measures) necessarily increases the degree of metapopulation synchrony. We show that decreases in effective between‐patch distance may deter population synchrony because of the simultaneous effect this change has on the migration patterns of other species. These results indicate that species interactions need to be considered when the effect of conservation measures on population synchrony, and ultimately persistence, is addressed.     

Joint estimation of contemporary seed and pollen dispersal rates among plant populations

There are few statistical methods for estimating contemporary dispersal among plant populations. A maximum-likelihood procedure is introduced here that uses pre- and post-dispersal population samples of biparentally inherited genetic markers to jointly estimate contemporary seed and pollen immigration rates from a set of discrete external sources into a target population. Monte Carlo simulations indicate that accurate estimates and reliable confidence intervals can be obtained using this method for both pollen and seed migration rates at modest sample sizes (100 parents/population and 100 offspring) when population differentiation is moderate (F(ST) ≥ 0.1), or by increasing pre-dispersal samples (to about 500 parents/population) when genetic divergence is weak (F(ST) = 0.01). The method exhibited low sensitivity to the number of source populations and achieved good accuracy at affordable genetic resolution (10 loci with 10 equifrequent alleles each). Unsampled source populations introduced positive biases in migration rate estimates from sampled sources, although they were minor when the proportion of immigration from the latter was comparatively low. A practical application of the method to a metapopulation of the Australian resprouter shrub Banksia attenuata revealed comparable levels of directional seed and pollen migration among dune groups, and the estimate of seed dispersal was higher than a previous estimate based on conservative assignment tests. The method should be of interest to researchers and managers assessing broad-scale nonequilibrium seed and pollen gene flow dynamics in plants.     

Single-species metapopulation dynamics: concepts, models and observations

This paper outlines a conceptual and theoretical framework for single-species metapopulation dynamics based on the Levins model and its variants. The significance of the following factors to metapopulation dynamics are explored: evolutionary changes in colonization ability; habitat patch size and isolation; compensatory effects between colonization and extinction rates; the effect of immigration on local dynamics (the rescue effect); and heterogeneity among habitat patches. The rescue effect may lead to alternative stable equilibria in metapopulation dynamics. Heterogeneity among habitat patches may give rise to a bimodal equilibrium distribution of the fraction of patches occupied in an assemblage of species (the core-satellite distribution). A new model of incidence functions is described, which allows one to estimate species' colonization and extinction rates on islands colonized from mainland. Four distinct kinds of stochasticity affecting metapopulation dynamics are discussed with examples. The concluding section describes four possible scenarios of metapopulation extinction.     

Predicting immigration of two species in contrasting landscapes: effects of scale, patch size and isolation

Migration is a key process for spatially structured populations. We examined how a variety of patch based metrics commonly used to predict the number of immigrants to a habitat patch performed based on data from three different years, in two distinct insect systems. The first system was an herbivorous beetle inhabiting patches of its host plant within a 'typical' patch network. In this system there were numerous patches located relatively close to one another, given the beetle's dispersal ability. The second system consisted of a butterfly inhabiting a series of 17 subalpine meadows. Here, the patches were arranged in a linear fashion and were more distant from each other. Overall, we found that the best models incorporating aspects of patch size and/or isolation explained a large (30–40%) amount of deviance in immigration, but there were considerable differences between the systems. For the first system, we found that metrics including the size of the target patch explained the highest proportion of deviance in immigrant numbers, while metrics based only on interpatch distances explained very little deviance. The situation was reversed for the second system. Metrics including the size of the target patch explained little deviance, while metrics based on the distance between patches explained the bulk of deviance in the number of immigrants. The results of our study show that the effects of patch size and isolation on the number of immigrants are highly important, but dependent on spatial scale, the organism studied, and how it responds to the spatial arrangement of patches. Correspondingly, there will be no single generalized metric to predict immigration for all cases. Given the dependency of the results on the system studied, we recommend that future studies provide explicit data on habitat areas and dispersal distance relative to interpatch distance to allow for meaningful comparison among organisms and systems.     

Population differentiation and gene flow within a metapopulation of a threatened tree, Magnolia stellata (Magnoliaceae)

We examined genetic differentiation among eight local populations of a metapopulation of Magnolia stellata using 10 nuclear and three chloroplast microsatellite (nSSR and cpSSR) markers and evaluated the influence of historical gene flow on population differentiation. The coefficient of genetic differentiation among populations for nSSR (F(ST) = 0.053) was less than half that for cpSSR (0.137). An isolation-by-distance pattern was detected for nSSRs, but not cpSSRs. These results suggest that pollen flow, as well as seed dispersal, has significantly reduced genetic differentiation among populations. We also examined patterns of contemporary pollen flow by paternity analysis of seeds from nine seed parents in one of the populations using the nSSR markers and found it to be greatly restricted by the distance between parents. Although most pollen flow occurred within the population, pollen flow from outside the population accounted for 2.5% of the total. When historical and contemporary pollen flows among populations were compared, the levels of pollen flow seem to have declined recently. We conclude that to conserve M. stellata, it is important to preserve the whole population by maintaining its metapopulation structure and the gene flow among its populations.     

Maintenance of high pollen dispersal in Eucalyptus wandoo, a dominant tree of the fragmented agricultural region in Western Australia

Pollen dispersal was investigated in five remnant populations of Eucalyptus wandoo, a dominant insect-pollinated tree in the fragmented agricultural region of southern Western Australia. Paternity analysis using six microsatellite loci identified a pollen source for 45% of seedlings, and the remainder were assumed to have arisen from pollen sources outside the stands. Outcrossing was variable, ranging from 52 to 89%, and long distance pollen dispersal was observed in all populations with up to 65% of pollen sourced from outside the populations over distances of at least 1 km. Modelling dispersal functions for pollination events within the two larger populations showed little difference between the four two-parameter models tested and indicated a fat-tailed dispersal curve. Similarity of direct and indirect historical estimates of gene flow indicates maintenance of gene flow at levels experienced prior to fragmentation. The study revealed extensive long distance pollen dispersal in remnant patches of trees within a fragmented agricultural landscape in the southern temperate region and highlighted the role of remnant patches in maintaining genetic connectivity at the landscape scale.     

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.     

Transgenic GFP as a molecular marker for approaches to quantify pollination mechanism and gene flow in Arabidopsis thaliana

Arabidopsis thaliana is commonly regarded as a self-pollinated plant. We observed that the stigma in each flower of A. thaliana cannot be pollinated by its own pollen in the early phases of the flowering process, when the anthers had dehisced but the filaments were still too short for the pollen to be deposited on the stigma. In the later stages, after elongation of the filaments, self-pollination can occur. After artificial pollination of the flower of a wild plant with GFP transgenic pollen grains in earlier stages of flowering, GFP expressed within epidermal cells was detected in some of the offspring (26.1-57.1 %). Wind-mediated pollen dispersal was poor but is likely to exist in natural habitats, while insects were observed visiting flowers of A. thaliana in natural and experimental populations. We constructed an experimental population consisting of 28 GFP transgenic plants and 240 wild plants and examined gene flow in the population. The result was that the distance of gene flow was limited to 0.5 m. 22 offspring with expressed GFP were found in 28,299 filial individuals examined, which suggested a relatively low outcrossing rate (0.74%). We conclude that outcrossing in populations of A. thaliana is mainly due to insect pollination. The data on gene flow could be useful to assess the ecological hazards of experimental transgene combinations.     

Metapopulation dynamics of a flightless alpine insect Hemideina maori in a naturally fragmented habitat

Abstract 1. Despite widespread acceptance of metapopulation theory, the effects that inter-patch dispersal and variability in patch size have on metapopulation dynamics in insects are two issues that require further study. In addition, previous studies of metapopulations have tended to focus on organisms with high dispersal capabilities such as some species of butterfly and bird.
2. Mountain stone weta Hemideina maori are a long-lived, flightless orthopteran that live on island rock outcrops or tors in the alpine region of southern New Zealand. A total of 480 adults and 789 juveniles was marked over three seasons on four large and 14 small tors to assess the effects of habitat fragmentation on the population dynamics of H. maori .
3. Only 12 adults (2.5% of marked adults and 4.0% of recaptured adults) and two juveniles (0.3% of marked juveniles and 0.7% of recaptured juveniles) dispersed between tors. The mean dispersal distance was 361 m (range = 36–672 m). Larger tors supported larger populations and had a higher number of emigrants and immigrants while smaller tors had proportionally higher emigration and immigration rates. Although adults on large and small tors had similar mean lifespans, five extinction events and three recolonisation events occurred during the study period, all on small tors.
4. Hemideina maori conform to many of the predictions of metapopulation theory even though they are flightless, show relatively low dispersal rates, and occur at low densities. Extinction and colonisation events are more common on small tors but may be relatively unimportant for the long-term survival of the metapopulation because they occur on the smallest habitat patches, which support the smallest proportion of the overall population.     

Evolutionary suicide and evolution of dispersal in structured metapopulations

 We study the evolution of dispersal in a structured metapopulation model. The metapopulation consists of a large (infinite) number of local populations living in patches of habitable environment. Dispersal between patches is modelled by a disperser pool and individuals in transit between patches are exposed to a risk of mortality. Occasionally, local catastrophes eradicate a local population: all individuals in the affected patch die, yet the patch remains habitable. We prove that, in the absence of catastrophes, the strategy not to migrate is evolutionarily stable. Under a given set of environmental conditions, a metapopulation may be viable and yet selection may favor dispersal rates that drive the metapopulation to extinction. This phenomenon is known as evolutionary suicide. We show that in our model evolutionary suicide can occur for catastrophe rates that increase with decreasing local population size. Evolutionary suicide can also happen for constant catastrophe rates, if local growth within patches shows an Allee effect. We study the evolutionary bifurcation towards evolutionary suicide and show that a discontinuous transition to extinction is a necessary condition for evolutionary suicide to occur. In other words, if population size smoothly approaches zero at a boundary of viability in parameter space, this boundary is evolutionarily repelling and no suicide can occur. Received: 10 November 2000 / Revised version: 13 February 2002 / Published online: 17 July 2002     

Population biology of Avena

Summary Pollen and seed dispersal patterns were analyzed in both natural and experimental populations of Avena barbata. Localized estimates of gene flow rates and plant densities gave estimates of neighborhood size in the range of 40 to 400 plants; the estimates of mean rate and distance of gene flow seemed to vary widely due to variable wind direction, rodent activity, microsite heterogeneity, etc. The relative sizes of neighborhoods in several populations were correlated with the patchy distribution of different genotypes (scored for lemma color and leaf sheath hairiness) within short distances, but patch sizes had a wide range among different sites. Highly localized gene flow patterns seemed to account for the observed pattern of highly patchy variation even when the dispersal curves for both pollen and seed were platykurtic in many cases. Measures of the stability of patches in terms of their size, dispersion in space and genetic structure in time are needed in order to sort out the relative roles of founder effects, random drift (due to small neighborhood size), and highly localized selection. However, our observations suggest that many variables and stochastic processes are involved in such studies so as to allow only weak inference about the underlying role of natural selection, drift and factors of population regulatien.     

Spatial Heterogeneity in Habitat Quality and Cross-Scale Interactions in Metapopulations

Abstract Integration of habitat heterogeneity into spatially realistic metapopulation approaches reveals the potential for key cross-scale interactions. Broad-scale environmental gradients and land-use practices can create autocorrelation of habitat quality of suitable patches at intermediate spatial scales. Patch occupancy then depends not only on habitat quality at the patch scale but also on feedbacks from surrounding neighborhoods of autocorrelated patches. Metapopulation dynamics emerge from how demographic and dispersal processes interact with relevant habitat heterogeneity. We provide an empirical example from a metapopulation of round-tailed muskrats (Neofiber alleni) in which habitat quality of suitable patches was spatially autocorrelated most strongly within 1,000 m, which was within the expected dispersal range of the species. After controlling for factors typically considered in metapopulation studies—patch size, local patch quality, patch connectivity—we use a cross-variogram analysis to demonstrate that patch occupancy by muskrats was correlated with habitat quality across scales ≤1,171 m. We also discuss general consequences of spatial heterogeneity of habitat quality for metapopulations related to potential cross-scale interactions. We focus on spatially correlated extinctions and metapopulation persistence, hierarchical scaling of source–sink dynamics, and dispersal decisions by individuals in relation to information constraints.     

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F_ST=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5–30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.