Pollen flow in the wildservice tree, Sorbus torminalis (L.) Crantz. IV. Whole interindividual variance of male fecundity estimated jointly with the dispersal kernel

Interindividual variance of male reproductive success (MRS) contributes to genetic drift, which in turn interacts with selection and migration to determine the short-term response of populations to rapid changes in their environment. Individual relative MRS can be estimated through paternity analysis and can be further dissected into fecundity and spatial components. Existing methods to achieve this decomposition either rely on the strong assumption of a random distribution of pollen donors (TwoGener) or estimate only the part of the variance of male fecundity that is explained by few covariates. We developed here a method to estimate jointly the whole variance of male fecundity and the pollen dispersal curve from the genotypic information of sampled seeds and their putative fathers and geographical information of all individuals in the study area. We modelled the relative individual fecundities as a log-normally distributed random effect. We used a Bayesian approach, well suited to the hierarchical nature of the model, to estimate these fecundities. When applied to Sorbus torminalis , the estimated variance of male fecundity corresponded to an effective density of trees 13 times lower than the observed density ( d_obs / d_ep ~  13). This value is between the value ( ~ 2) estimated with a classical mating model including three covariates (neighbourhood density, diameter, flowering intensity) that affect fecundity and the value ( ~ 30) estimated with TwoGener. The estimated dispersal kernel was close to previous results. This approach allows fine monitoring of ongoing genetic drift in natural populations, and quantitative dissection of the processes contributing to drift, including human actions.     

Pollen flow in the wildservice tree, Sorbus torminalis (L.) Crantz. II. Pollen dispersal and heterogeneity in mating success inferred from parent-offspring analysis

Knowing the extent of gene movements from parents to offspring is essential to understand the potential of a species to adapt rapidly to a changing environment, and to design appropriate conservation strategies. In this study, we develop a nonlinear statistical model to jointly estimate the pollen dispersal kernel and the heterogeneity in fecundity among phenotypically or environmentally defined groups of males. This model uses genotype data from a sample of fruiting plants, a sample of seeds harvested on each of these plants, and all males within a circumscribed area. We apply this model to a scattered, entomophilous woody species, Sorbus torminalis (L.) Crantz, within a natural population covering more than 470 ha. We estimate a high heterogeneity in male fecundity among ecological groups, both due to phenotype (size of the trees and flowering intensity) and landscape factors (stand density within the neighbourhood). We also show that fat-tailed kernels are the most appropriate to depict the important abilities of long-distance pollen dispersal for this species. Finally, our results reveal that the spatial position of a male with respect to females affects as much its mating success as ecological determinants of male fecundity. Our study thus stresses the interest to account for the dispersal kernel when estimating heterogeneity in male fecundity, and reciprocally.     

Paternity analysis-based inference of pollen dispersal patterns, male fecundity variation, and influence of flowering tree density and general flowering magnitude in two dipterocarp species

Background and Aims

Knowledge of pollen dispersal patterns and variation of fecundity is essential to understanding plant evolutionary processes and to formulating strategies to conserve forest genetic resources. Nevertheless, the pollen dispersal pattern of dipterocarp, main canopy tree species in palaeo-tropical forest remains unclear, and flowering intensity variation in the field suggests heterogeneity of fecundity.

Methods

Pollen dispersal patterns and male fecundity variation of Shorea leprosula and Shorea parvifolia ssp. parvifolia on Peninsular Malaysian were investigated during two general flowering seasons (2001 and 2002), using a neighbourhood model modified by including terms accounting for variation in male fecundity among individual trees to express heterogeneity in flowering.

Key Results

The pollen dispersal patterns of the two dipterocarp species were affected by differences in conspecific tree flowering density, and reductions in conspecific tree flowering density led to an increased selfing rate. Active pollen dispersal and a larger number of effective paternal parents were observed for both species in the season of greater magnitude of general flowering (2002).

Conclusions

The magnitude of general flowering, male fecundity variation, and distance between pollen donors and mother trees should be taken into account when attempting to predict the effects of management practices on the self-fertilization and genetic structure of key tree species in tropical forest, and also the sustainability of possible management strategies, especially selective logging regimes.     

Fine-scale genetic structure and gene dispersal inferences in 10 neotropical tree species

The extent of gene dispersal is a fundamental factor of the population and evolutionary dynamics of tropical tree species, but directly monitoring seed and pollen movement is a difficult task. However, indirect estimates of historical gene dispersal can be obtained from the fine-scale spatial genetic structure of populations at drift-dispersal equilibrium. Using an approach that is based on the slope of the regression of pairwise kinship coefficients on spatial distance and estimates of the effective population density, we compare indirect gene dispersal estimates of sympatric populations of 10 tropical tree species. We re-analysed 26 data sets consisting of mapped allozyme, SSR (simple sequence repeat), RAPD (random amplified polymorphic DNA) or AFLP (amplified fragment length polymorphism) genotypes from two rainforest sites in French Guiana. Gene dispersal estimates were obtained for at least one marker in each species, although the estimation procedure failed under insufficient marker polymorphism, limited sample size, or inappropriate sampling area. Estimates generally suffered low precision and were affected by assumptions regarding the effective population density. Averaging estimates over data sets, the extent of gene dispersal ranged from 150 m to 1200 m according to species. Smaller gene dispersal estimates were obtained in species with heavy diaspores, which are presumably not well dispersed, and in populations with high local adult density. We suggest that limited seed dispersal could indirectly limit effective pollen dispersal by creating higher local tree densities, thereby increasing the positive correlation between pollen and seed dispersal distances. We discuss the potential and limitations of our indirect estimation procedure and suggest guidelines for future studies.     

Using genetic markers to directly estimate gene flow and reproductive success parameters in plants on the basis of naturally regenerated seedlings

Estimating seed and pollen gene flow in plants on the basis of samples of naturally regenerated seedlings can provide much needed information about "realized gene flow," but seems to be one of the greatest challenges in plant population biology. Traditional parentage methods, because of their inability to discriminate between male and female parentage of seedlings, unless supported by uniparentally inherited markers, are not capable of precisely describing seed and pollen aspects of gene flow realized in seedlings. Here, we describe a maximum-likelihood method for modeling female and male parentage in a local plant population on the basis of genotypic data from naturally established seedlings and when the location and genotypes of all potential parents within the population are known. The method models female and male reproductive success of individuals as a function of factors likely to influence reproductive success (e.g., distance of seed dispersal, distance between mates, and relative fecundity--i.e., female and male selection gradients). The method is designed to account for levels of seed and pollen gene flow into the local population from unsampled adults; therefore, it is well suited to isolated, but also wide-spread natural populations, where extensive seed and pollen dispersal complicates traditional parentage analyses. Computer simulations were performed to evaluate the utility and robustness of the model and estimation procedure and to assess how the exclusion power of genetic markers (isozymes or microsatellites) affects the accuracy of the parameter estimation. In addition, the method was applied to genotypic data collected in Scots pine (isozymes) and oak (microsatellites) populations to obtain preliminary estimates of long-distance seed and pollen gene flow and the patterns of local seed and pollen dispersal in these species.     

Fluorescent dye particles as pollen analogues for measuring pollen dispersal in an insect-pollinated forest herb

In flowering plants, pollen dispersal is often the major contributing component to gene flow, hence a key parameter in conservation genetics and population biology. A cost-effective method to assess pollen dispersal consists of monitoring the dispersal of fluorescent dyes used as pollen analogues. However, few comparisons between dye dispersal and realized pollen dispersal have been performed to validate the method. We investigated pollen dispersal in two small populations of the insect-pollinated herb Primula elatior from urban forest fragments using direct (paternity analyses based on microsatellite DNA markers) and indirect (fluorescent dyes) methods. We compared these methods using two approaches, testing for the difference between the distance distributions of observed dispersal events and estimating parameters of a dispersal model, and related these results to dye dispersal patterns in three large populations. Dye and realized (based on paternity inference) pollen dispersal showed exponential decay distributions, with 74.2?C94.8% of the depositions occurring at <50?m and a few longer distance dispersal events (up to 151?m). No significant difference in curve shape was found between dye and realized pollen dispersal distributions. The best-fitting parameters characterizing the dye dispersal model were consistent with those obtained for realized pollen dispersal. Hence, the fluorescent dye method may be considered as reliable to infer realized pollen dispersal for forest herbs such as P. elatior. However, our simulations reveal that large sample sizes are needed to detect moderate differences between dye and realized pollen dispersal patterns because the estimation of dispersal parameters suffers low precision.     

Disentangling the causes of heterogeneity in male fecundity in gynodioecious Beta vulgaris ssp. maritima

Variation among individuals in reproductive success is advocated as a major process driving evolution of sexual polymorphisms in plants, such as gynodioecy where females and hermaphrodites coexist. In gynodioecious Beta vulgaris ssp. maritima, sex determination involves cytoplasmic male sterility (CMS) genes and nuclear restorers of male fertility. Both restored CMS and non-CMS hermaphrodites co-occur. Genotype-specific differences in male fitness are theoretically expected to explain the maintenance of cytonuclear polymorphism. Using genotypic information on seedlings and flowering plants within two metapopulations, we investigated whether male fecundity was influenced by ecological, phenotypic and genetic factors, while taking into account the shape and scale of pollen dispersal. Along with spatially restricted pollen flow, we showed that male fecundity was affected by flowering synchrony, investment in reproduction, pollen production and cytoplasmic identity of potential fathers. Siring success of non-CMS hermaphrodites was higher than that of restored CMS hermaphrodites. However, the magnitude of the difference in fecundity depended on the likelihood of carrying restorer alleles for non-CMS hermaphrodites. Our results suggest the occurrence of a cost of silent restorers, a condition supported by scarce empirical evidence, but theoretically required to maintain a stable sexual polymorphism in gynodioecious species.     

Limited pollen dispersal and biparental inbreeding in Symphonia globulifera in French Guiana

In this paper, we report a study of the mating system and gene flow of Symphonia globulifera, a hermaphroditic, mainly bird-pollinated tree species with a large geographic distribution in the tropical Americas and Africa. Using three microsatellites, we analysed 534 seeds of 28 open pollinated families and 164 adults at the experimental site 'Paracou' in French Guiana. We observed, compared to other tropical tree species, relatively high values for the effective number of alleles. Significant spatial genetic structure was detected, with trees at distances up to 150 m more genetically similar than expected at random. We estimated parameters of the mating system and gene flow by using the mixed mating model and the TwoGener approach. The estimated multilocus outcrossing rate, tm, was 0.920. A significant level of biparental inbreeding and a high proportion of full-sibs were estimated for the 28 seed arrays. We estimated mean pollen dispersal distances between 27 and 53 m according to the dispersal models used. Although the adult population density of S. globulifera in Paracou was relatively high, the joint estimation of pollen dispersal and density of reproductive trees gave effective density estimates of 1.6 and 1.3 trees/ha. The parameters of the mating system and gene flow are discussed in the context of spatial genetic and demographic structures, flowering phenology and pollinator composition and behaviour.     

Temporal variation in pollen dispersal and breeding structure in a bee-pollinated Neotropical tree

Variation among flowering seasons in the time of flowering, synchrony and length of flowering, and fluctuations in the abundance of pollinators may cause a variation in pollen dispersal distance. In this study, we analyzed the temporal variation in pollen dispersal and breeding structure in the Neotropical tree species Tabebuia aurea (Bignoniaceae) and evaluated pollen dispersal between a population inside the reserve and a patch of isolated individuals on the edge of the reserve, and tested the hypothesis that isolated individuals are sinking for pollen. All adult trees (260) within a population of 40 ha and 9 isolated individuals on the edge of the reserve were sampled, and from these adults, 21 open-pollinated progeny arrays were analyzed in 2 flowering seasons (309 seeds in 2004 and 328 in 2005). Genetic analyses were based on the polymorphism at 10 microsatellite loci. A high proportion of self-pollination found in both flowering seasons indicated a mixed-mating system. The mean pollen dispersal distance differed significantly between the two flowering seasons (307.78 m in 2004 and 396.26 m in 2005). Maximum pollen dispersal was 2608 m, but most pollination events (65%) occurred at distances <300 m. Our results also showed that isolated individuals are sinking for pollen, with high pollen flow between the population inside the reserve and individuals on the edge. These results are most likely due to the large pollinator species, which can potentially fly long distances, and also due to temporal variation in individual fecundity and contribution to pollen dispersal.     

The control of natural variation in cytosine methylation in Arabidopsis

The distance of pollen movement is an important determinant of the neighborhood area of plant populations. In earlier studies, we designed a method for estimating the distance of pollen dispersal, on the basis of the analysis of the differentiation among the pollen clouds of a sample of females, spaced across the landscape. The method was based solely on an estimate of the global level of differentiation among the pollen clouds of the total array of sampled females. Here, we develop novel estimators, on the basis of the divergence of pollen clouds for all pairs of females, assuming that an independent estimate of adult population density is available. A simulation study shows that the estimators are all slightly biased, but that most have enough precision to be useful, at least with adequate sample sizes. We show that one of the novel pairwise methods provides estimates that are slightly better than the best global estimate, especially when the markers used have low exclusion probability. The new method can also be generalized to the case where there is no prior information on the density of reproductive adults. In that case, we can jointly estimate the density itself and the pollen dispersal distance, given sufficient sample sizes. The bias of this last estimator is larger and the precision is lower than for those estimates based on independent estimates of density, but the estimate is of some interest, because a meaningful independent estimate of the density of reproducing individuals is difficult to obtain in most cases.     

Using genetic markers to estimate the pollen dispersal curve

Pollen dispersal is a critical process that shapes genetic diversity in natural populations of plants. Estimating the pollen dispersal curve can provide insight into the evolutionary dynamics of populations and is essential background for making predictions about changes induced by perturbations. Specifically, we would like to know whether the dispersal curve is exponential, thin-tailed (decreasing faster than exponential), or fat-tailed (decreasing slower than the exponential). In the latter case, rare events of long-distance dispersal will be much more likely. Here we generalize the previously developed TWOGENER method, assuming that the pollen dispersal curve belongs to particular one- or two-parameter families of dispersal curves and estimating simultaneously the parameters of the dispersal curve and the effective density of reproducing individuals in the population. We tested this method on simulated data, using an exponential power distribution, under thin-tailed, exponential and fat-tailed conditions. We find that even if our estimates show some bias and large mean squared error (MSE), we are able to estimate correctly the general trend of the curve - thin-tailed or fat-tailed - and the effective density. Moreover, the mean distance of dispersal can be correctly estimated with low bias and MSE, even if another family of dispersal curve is used for the estimation. Finally, we consider three case studies based on forest tree species. We find that dispersal is fat-tailed in all cases, and that the effective density estimated by our model is below the measured density in two of the cases. This latter result may reflect the difficulty of estimating two parameters, or it may be a biological consequence of variance in reproductive success of males in the population. Both the simulated and empirical findings demonstrate the strong potential of TWOGENER for evaluating the shape of the dispersal curve and the effective density of the population (d(e)).     

Ecological implications of a flower size/number trade-off in tropical forest trees

BACKGROUND: In angiosperms, flower size commonly scales negatively with number. The ecological consequences of this trade-off for tropical trees remain poorly resolved, despite their potential importance for tropical forest conservation. We investigated the flower size number trade-off and its implications for fecundity in a sample of tree species from the Dipterocarpaceae on Borneo. METHODOLOGY/PRINCIPAL FINDINGS: We combined experimental exclusion of pollinators in 11 species, with direct and indirect estimates of contemporary pollen dispersal in two study species and published estimates of pollen dispersal in a further three species to explore the relationship between flower size, pollinator size and mean pollen dispersal distance. Maximum flower production was two orders of magnitude greater in small-flowered than large-flowered species of Dipterocarpaceae. In contrast, fruit production was unrelated to flower size and did not differ significantly among species. Small-flowered species had both smaller-sized pollinators and lower mean pollination success than large-flowered species. Average pollen dispersal distances were lower and frequency of mating between related individuals was higher in a smaller-flowered species than a larger-flowered confamilial. Our synthesis of pollen dispersal estimates across five species of dipterocarp suggests that pollen dispersal scales positively with flower size. CONCLUSIONS AND THEIR SIGNIFICANCE: Trade-offs embedded in the relationship between flower size and pollination success contribute to a reduction in the variance of fecundity among species. It is therefore plausible that these processes could delay competitive exclusion and contribute to maintenance of species coexistence in this ecologically and economically important family of tropical trees. These results have practical implications for tree species conservation and restoration. Seed collection from small-flowered species may be especially vulnerable to cryptic genetic erosion. Our findings also highlight the potential for differential vulnerability of tropical tree species to the deleterious consequences of forest fragmentation.     

Comparison of pollen gene flow among four European beech (Fagus sylvatica L.) populations characterized by different management regimes

The study of the dispersal capability of a species can provide essential information for the management and conservation of its genetic variability. Comparison of gene flow rates among populations characterized by different management and evolutionary histories allows one to decipher the role of factors such as isolation and tree density on gene movements. We used two paternity analysis approaches and different strategies to handle the possible presence of genotyping errors to obtain robust estimates of pollen flow in four European beech (Fagus sylvatica L.) populations from Austria and France. In each country one of the two plots is located in an unmanaged forest; the other plots are managed with a shelterwood system and inside a colonization area (in Austria and France, respectively). The two paternity analysis approaches provided almost identical estimates of gene flow. In general, we found high pollen immigration (~75% of pollen from outside), with the exception of the plot from a highly isolated forest remnant (~50%). In the two unmanaged plots, the average within-population pollen dispersal distances (from 80 to 184 m) were higher than previously estimated for beech. From the comparison between the Austrian managed and unmanaged plots, that are only 500 m apart, we found no evidence that either gene flow or reproductive success distributions were significantly altered by forest management. The investigated phenotypic traits (crown area, height, diameter and flowering phenology) were not significantly related with male reproductive success. Shelterwood seems to have an effect on the distribution of within-population pollen dispersal distances. In the managed plot, pollen dispersal distances were shorter, possibly because adult tree density is three-fold (163 versus 57 trees per hectare) with respect to the unmanaged one.     

Mixed Mating System Are Regulated by Fecundity in Shorea curtisii (Dipterocarpaceae) as Revealed by Comparison under Different Pollen Limited Conditions

The maintenance of mixed mating was studied in Shorea curtisii, a dominant and widely distributed dipterocarp species in Southeast Asia. Paternity and hierarchical Bayesian analyses were used to estimate the parameters of pollen dispersal kernel, male fecundity and self-pollen affinity. We hypothesized that partial self incompatibility and/or inbreeding depression reduce the number of selfed seeds if the mother trees receive sufficient pollen, whereas reproductive assurance increases the numbers of selfed seeds under low amounts of pollen. Comparison of estimated parameters of self-pollen affinity between high density undisturbed and low density selectively logged forests indicated that self-pollen was selectively excluded from mating in the former, probably due to partial self incompatibility or inbreeding depression until seed maturation. By estimating the self-pollen affinity of each mother tree in both forests, mother trees with higher amount of self-pollen indicated significance of self-pollen affinity with negative estimated value. The exclusion of self-fertilization and/or inbreeding depression during seed maturation occurred in the mother trees with large female fecundity, whereas reproductive assurance increased self-fertilization in the mother trees with lower female fecundity.     

Fine-scale genetic structure and gene dispersal in Centaurea corymbosa (Asteraceae). II. Correlated paternity within and among sibships

The fine-scale pattern of correlated paternity was characterized within a population of the narrow-endemic model plant species, Centaurea corymbosa, using microsatellites and natural progeny arrays. We used classical approaches to assess correlated mating within sibships and developed a new method based on pairwise kinship coefficients to assess correlated paternity within and among sibships in a spatio-temporal perspective. We also performed numerical simulations to assess the relative significance of different mechanisms promoting correlated paternity and to compare the statistical properties of different estimators of correlated paternity. Our new approach proved very informative to assess which factors contributed most to correlated paternity and presented good statistical properties. Within progeny arrays, we found that about one-fifth of offspring pairs were full-sibs. This level of correlated mating did not result from correlated pollen dispersal events (i.e., pollen codispersion) but rather from limited mate availability, the latter being due to limited pollen dispersal distances, the heterogeneity of pollen production among plants, phenological heterogeneity and, according to simulations, the self-incompatibility system. We point out the close connection between correlated paternity and the "TwoGener" approach recently developed to infer pollen dispersal and discuss the conditions to be met when applying the latter.     

Estimating pollen flow using SSR markers and paternity exclusion: accounting for mistyping

Highly informative genetic markers, such as simple sequence repeats (SSRs), can be used to directly measure pollen flow by parentage analysis. However, mistyping (i.e. false inference of genotypes caused by the occurrence of null alleles, mutations, and detection errors) can lead to substantial biases in the estimates obtained. Using computer simulations, we evaluated a direct method for estimating pollen immigration using SSR markers and a paternity exclusion approach. This method accounts for mistyping and does not rely on assumptions about the distribution of male reproductive success. If ignored, even minor rates of mistyping (1.5%) resulted in overestimating pollen immigration by up to 150%. When we required at least two mismatching loci before excluding candidate fathers from paternity, the resulting pollen immigration estimates had small biases for rates of mistyping up to 4.5%. Requiring at least three mismatches for exclusion was needed to minimize the upward biases of pollen immigration caused by rates of mistyping up to 10.5%. The minimum number of highly variable SSR loci needed to minimize cryptic gene flow and obtain reliable estimates of pollen immigration varied from five to seven for a sampling scheme applicable to most conifers (i.e. when paternal haplotypes can be unambiguously determined). Between five and nine highly variable SSR loci were needed for a more general sampling scheme that is applicable to all diploid seed plants. With moderately variable SSR markers, consistently accurate estimates of pollen immigration could be obtained only for rates of mistyping up to 4.5%. We developed the POLLEN FLOW (PFL) computer program which can be used to obtain unbiased and precise estimates of pollen immigration under a wide range of conditions, including population sizes as large as 600 parents and mistyping rates as high as 10.5%.     

DIRECT AND INDIRECT ESTIMATES OF NEIGHBORHOOD AND EFFECTIVE POPULATION SIZE IN A TROPICAL PALM,ASTROCARYUM MEXICANUM

To estimate the relative importance of genetic drift, the effective population size ???(N_e) can be used. Here we present estimates of the effective population size and related measures in Astrocaryum mexicanum, a tropical palm from Los Tuxtlas rain forest, Veracruz, Mexico. Seed and pollen dispersal were measured. Seeds are primarily dispersed by gravity and secondarily dispersed by small mammals. Mean primary and secondary dispersal distances for seeds were found to be small (0.78 m and 2.35 m, respectively). A. mexicanum is beetle pollinated and pollen movements were measured by different methods: a) using fluorescent dyes, b) as the minimum distance between active female and male inflorescences, and c) using rare allozyme alleles as genetic markers. All three estimates of pollen dispersal were similar, with a mean of approximately 20 m. Using the seed and pollen dispersal data, the genetic neighborhood area (A) was estimated to be 2,551 m². To obtain the effective population size, three different overlapping generation methods were used to estimate an effective density with demographic data from six permanent plots. The effective density ranged from 0.040 to 0.351 individuals per m². The product of effective density and neighborhood area yields a direct estimate of the neighborhood effective population size (N_b). N_b ranged from 102 to 895 individuals. Indirect estimates of population size and migration rate (Nm) were obtained using F_st for five different allozymic loci for both adults and seeds. We obtained a range of Nm from 1.2 to 19.7 in adults and a range of Nm from 4.0 to 82.6 for seeds. We discuss possible causes of the smaller indirect estimates of Nm relative to the direct and compare our estimates with values from other plant populations. Gene dispersal distances, neighborhood size, and effective population size in A. mexicanum are relatively high, suggesting that natural selection, rather than genetic drift, may play a dominant role in patterning the genetic variation in this tropical palm.     

Parentage versus two-generation analyses for estimating pollen-mediated gene flow in plant populations

Assessment of contemporary pollen-mediated gene flow in plants is important for various aspects of plant population biology, genetic conservation and breeding. Here, through simulations we compare the two alternative approaches for measuring pollen-mediated gene flow: (i) the NEIGHBORHOOD model--a representative of parentage analyses, and (ii) the recently developed TWOGENER analysis of pollen pool structure. We investigate their properties in estimating the effective number of pollen parents (N(ep)) and the mean pollen dispersal distance (delta). We demonstrate that both methods provide very congruent estimates of N(ep) and delta, when the methods' assumptions considering the shape of pollen dispersal curve and the mating system follow those used in data simulations, although the NEIGHBORHOOD model exhibits generally lower variances of the estimates. The violations of the assumptions, especially increased selfing or long-distance pollen dispersal, affect the two methods to a different degree; however, they are still capable to provide comparable estimates of N(ep). The NEIGHBORHOOD model inherently allows to estimate both self-fertilization and outcrossing due to the long-distance pollen dispersal; however, the TWOGENER method is particularly sensitive to inflated selfing levels, which in turn may confound and suppress the effects of distant pollen movement. As a solution we demonstrate that in case of TWOGENER it is possible to extract the fraction of intraclass correlation that results from outcrossing only, which seems to be very relevant for measuring pollen-mediated gene flow. The two approaches differ in estimation precision and experimental efforts but they seem to be complementary depending on the main research focus and type of a population studied.     

Pollen-mediated gene dispersal within continuous and fragmented populations of a forest understorey species, Trillium cuneatum

Pollen movement plays a critical role in the distribution of genetic variation within and among plant populations. Direct measures of pollen movement in the large, continuous populations that characterize many herbaceous plant species are often technically difficult and biologically unreliable. Here, we studied contemporary pollen movement in four large populations of Trillium cuneatum. Three populations, located in the Georgia Piedmont, are exposed to strong anthropogenic disturbances, while the fourth population, located in the Southern Appalachian Mountains, is relatively undisturbed. Using the recently developed TwoGener analysis, we extracted estimates of the effective number of pollen donors (N(ep)), effective mating neighbourhood size (A(ep)) and the average distance of pollen movement (delta) for each population. We extended the TwoGener method by developing inference on the paternal gametic contribution to the embryo in situations where offspring genotypes are inferred from seeds and elaiosomes of species with bisporic megagametogenesis. Our estimates indicate that maternal plants do not sample pollen randomly from a global pool; rather, pollen movement in all four populations is highly restricted. Although the effective number of pollen donors per maternal plant is low (1.22-1.66) and pollen movement is highly localized in all populations, N(ep) in the disturbed Piedmont populations is higher and there is more pollen movement than in the mountains. The distance pollen moves is greater in disturbed sites and fragmented populations, possibly due to edge effects in Trillium habitats.     

Two-generation analysis of pollen flow across a landscape. I. Male gamete heterogeneity among females

Gene flow is a key factor in the spatial genetic structure in spatially distributed species. Evolutionary biologists interested in microevolutionary processess and conservation biologists interested in the impact of landscape change require a method that measures the real time process of gene movement. We present a novel two-generation (parent-offspring) approach to the study of genetic structure (TwoGener) that allows us to quantify heterogeneity among the male gamete pools sampled by maternal trees scattered across the landscape and to estimate mean pollination distance and effective neighborhood size. First, we describe the model's elements: genetic distance matrices to estimate intergametic distances, molecular analysis of variance to determine whether pollen profiles differ among mothers, and optimal sampling considerations. Second, we evaluate the model's effectiveness by simulating spatially distributed populations. Spatial heterogeneity in male gametes can be estimated by phiFT, a male gametic analogue of Wright's F(ST) and an inverse function of mean pollination distance. We illustrate TwoGener in cases where the male gamete can be categorically or ambiguously determined. This approach does not require the high level of genetic resolution needed by parentage analysis, but the ambiguous case is vulnerable to bias in the absence of adequate genetic resolution. Finally, we apply TwoGener to an empirical study of Quercus alba in Missouri Ozark forests. We find that phiFT = 0.06, translating into about eight effective pollen donors per female and an effective pollination neighborhood as a circle of radius about 17 m. Effective pollen movement in Q. alba is more restricted than previously realized, even though pollen is capable of moving large distances. This case study illustrates that, with a modest investment in field survey and laboratory analysis, the TwoGener approach permits inferences about landscape-level gene movements.