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.     

Estimating anisotropic pollen dispersal: a case study in Quercus lobata

The pollen dispersal distribution is an important element of the neighbourhood size of plant populations. Most methods aimed at estimating the dispersal curve assume that pollen dispersal is isotropic, but evidence indicates that this assumption does not hold for many plant species, particularly wind-pollinated species subject to prevailing winds during the pollination season. We propose here a method of detecting anisotropy of pollen dispersal and of gauging its intensity, based on the estimation of the differentiation of maternal pollen clouds (TWOGENER extraction), assuming that pollen dispersal is bivariate and normally distributed. We applied the new method to a case study in Quercus lobata, detecting only a modest level of anisotropy in pollen dispersal in a direction roughly similar to the prevailing wind direction. Finally, we conducted a simulation to explore the conditions under which anisotropy can be detected with this method, and we show that while anisotropy is detectable, in principle, it requires a large volume of data.     

Patterns of pollen flow in a dense population of the insect-pollinated canopy tree species Castanopsis sieboldii

Insect pollinations of tree species with high-density populations have rarely been studied. Since the density of adults can affect effective pollen dispersal, short-distance pollination, even by insects, may frequently occur in high-density populations. To test this prediction, we investigated pollination patterns in a high-density population of the insect-pollinated canopy tree species Castanopsis sieboldii by paternity analysis using genotypes at 8 microsatellite loci of 145 adult trees and 439 seeds from 11 seed parents in a 4-ha plot. We then explored their genetic effects on the population by calculating other population genetics parameters. Although C. sieboldii has high potential for long-distance dispersal of pollen (as indicated by a fat-tailed dispersal kernel), the cumulative pollination at the local scale was spatially limited and strongly dependent on the distance between parents due to the high density of adults. Genetic diversity estimates for pollen pools accepted by each seed parent converged on a maximum as the effective number of pollen parents increased. The genetic diversity of pollen pool bulked over all the seed parents from inside the plot did not differ from that of the total pollen pools. Therefore, although pollen flow from distant pollen parents may help to maintain the genetic diversity of offspring, pollen parents neighboring seed parents may be the main contributors to the genetic diversity of the offspring at the seed stage.     

Two-generation analysis of pollen flow across a landscape V: a stepwise approach for extracting factors contributing to pollen structure

Patterns of pollen dispersal are central to both the ecology and evolution of plant populations. However, the mechanisms controlling either the dispersal process itself or our estimation of that process may be influenced by site-specific factors such as local forest structure and nonuniform adult genetic structure. Here, we present an extension of the AMOVA model applied to the recently developed TWOGENER analysis of pollen pool structure. This model, dubbed the Stepwise AMOVA (StAMOVA), focuses on determining to what extent ecological, demographic, and/or environmental factors influence the observed genetic variation in spatially separated pollen pools. The analysis is verified for efficacy, using an extensive battery of simulations, illustrating: (1) how nonuniform adult genetic structure influences the differentiation of spatially separated pollen pools, and (2) how effectively the Stepwise analysis performs in carrying out the appropriate corrections. Finally, the model is applied to a Quercus alba data set, from which we have prior evidence that the adult genetic structure is nonuniformly distributed across the sampling landscape. From this data set, we show how the Stepwise model can be applied to remove the effects of spatial adult genetic structure on pollen pool differentiation and contrast these results with those derived from the original TWOGENER analysis.     

A new method of estimating the pollen dispersal curve independently of effective density

We introduce a novel indirect method of estimating the pollen dispersal curve from mother-offspring genotypic data. Unlike an earlier indirect approach (TwoGener), this method is based on a normalized measure of correlated paternity between female pairs whose expectation does not explicitly depend on the unknown effective male population density (d(e)). We investigate the statistical properties of the new method, by comparison with those of TwoGener, considering the sensitivity to reductions of d(e), relative to census density, resulting from unequal male fecundity and asynchronous flowering. Our main results are: (i) it is possible to obtain reliable estimates of the average distance of pollen dispersal, delta, from indirect methods, even under nonuniform male fecundity and variable flowering phenology; (ii) the new method yields more accurate and more precise delta-estimates than TwoGener under a wide range of sampling and flowering scenarios; and (iii) TwoGener can be used to obtain approximate d(e) estimates, if needed for other purposes. Our results also show that accurately estimating the shape of the tail of the pollen dispersal function by means of indirect methods remains a very difficult challenge.     

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.     

The pollen dispersal kernel and mating system of an insect-pollinated tropical palm,Oenocarpus bataua

Pollen dispersal shapes the local genetic structure of plant populations and determines the opportunity for local selection and genetic drift, but has been well studied in few animal-pollinated plants in tropical rainforests. Here, we characterise pollen movement for an insect-pollinated Neotropical canopy palm, Oenocarpus bataua, and relate these data to adult mating system and population genetic structure. The study covers a 130-ha parcel in which all adult trees (n=185) were mapped and genotyped at 12 microsatellite loci, allowing us to positively identify the source tree for 90% of pollination events (n=287 of 318 events). Mating system analysis showed O. bataua was effectively outcrossed (t_m=1.02) with little biparental inbreeding (t_m−t_s=−0.005) and an average of 5.4 effective pollen donors (N_ep) per female. Dispersal distances were relatively large for an insect-pollinated species (mean=303 m, max=1263 m), and far exceeded nearest-neighbour distances. Dispersal kernel modelling indicated a thin-tailed Weibull distribution offered the best fit to the genetic data, which contrasts with the fat-tailed kernels typically reported for pollen dispersal in trees. Preliminary analyses suggest that our findings may be explained, at least in part, by a relatively diffuse spatial and temporal distribution of flowering trees. Comparison with previously reported estimates of seed movement for O. bataua suggests that pollen and seed dispersal distances may be similar. These findings add to the growing body of information on dispersal in insect-pollinated trees, but underscore the need for continued research on tropical systems in general, and palms in particular.     

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.     

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long-distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like GST are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D2 is linearly inversely related to the power-law exponent, with a slope of approximately -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.     

Paternity analysis of pollen-mediated gene flow for Fraxinus excelsior L. in a chronically fragmented landscape

Paternity analysis based on microsatellite marker genotyping was used to infer contemporary genetic connectivity by pollen of three population remnants of the wind-pollinated, wind-dispersed tree Fraxinus excelsior, in a deforested Scottish landscape. By deterministically accounting for genotyping error and comparing a range of assignment methods, individual-based paternity assignments were used to derive population-level estimates of gene flow. Pollen immigration into a 300 ha landscape represents between 43 and 68% of effective pollination, mostly depending on assignment method. Individual male reproductive success is unequal, with 31 of 48 trees fertilizing one seed or more, but only three trees fertilizing more than ten seeds. Spatial analysis suggests a fat-tailed pollen dispersal curve with 85% of detected pollination occurring within 100 m, and 15% spreading between 300 and 1900 m from the source. Identification of immigrating pollen sourced from two neighbouring remnants indicates further effective dispersal at 2900 m. Pollen exchange among remnants is driven by population size rather than geographic distance, with larger remnants acting predominantly as pollen donors, and smaller remnants as pollen recipients. Enhanced wind dispersal of pollen in a barren landscape ensures that the seed produced within the catchment includes genetic material from a wide geographic area. However, gene flow estimates based on analysis of non-dispersed seeds were shown to underestimate realized gene immigration into the remnants by a factor of two suggesting that predictive landscape conservation requires integrated estimates of post-recruitment gene flow occurring via both pollen and seed.     

Extensive pollen flow in a natural fragmented population of Patagonian cypress Austrocedrus chilensis

Habitat fragmentation might significantly affect mating and pollen dispersal patterns in plant populations, contributing to the decline of remnant populations. However, wind-pollinated species are able to disperse pollen at longer distances after opening of the canopy. Our objectives were to characterize the mating system parameters and to estimate the average distance of effective pollen dispersal in the wind-pollinated conifer Austrocedrus chilensis. We sampled 19 “mother trees,” 200 progeny, and 81 additional adults (both male and female), in a fragmented population at the Argentinean Patagonian steppe. We registered the spatial positions of individuals and genotyped all samples with five microsatellite markers. We found a high genetic diversity, a moderated rate of biparental inbreeding (t _m? ??t _s?=?0.105), and a complete absence of correlated paternity (r _p?=??0.015). The effective number of pollen donors contributing to a single mother (N _ep) was 13.9. Applying TWOGENER, we estimated a low but significant differentiation among the inferred pollen pools (Φ_FT?=?0.036, p?=?0.001) and a very large average pollen dispersal distance (d?=?1,032.3 m). The leptokurtic distribution (b?=?0.18) presumes a potential for even larger dispersal distances. The high genetic diversity, the mating patterns, and the extensive pollen dispersal presume that habitat fragmentation did not have a negative impact on pollen movement in this population of A. chilensis. Genetic connectivity among fragmented populations scattered in the Patagonian region is possible, and we stress the need of management policies at the landscape level.     

Do Spatially-Implicit Estimates of Neutral Migration Comply with Seed Dispersal Data in Tropical Forests?

Neutral community models have shown that limited migration can have a pervasive influence on the taxonomic composition of local communities even when all individuals are assumed of equivalent ecological fitness. Notably, the spatially implicit neutral theory yields a single parameter I for the immigration-drift equilibrium in a local community. In the case of plants, seed dispersal is considered as a defining moment of the immigration process and has attracted empirical and theoretical work. In this paper, we consider a version of the immigration parameter I depending on dispersal limitation from the neighbourhood of a community. Seed dispersal distance is alternatively modelled using a distribution that decreases quickly in the tails (thin-tailed Gaussian kernel) and another that enhances the chance of dispersal events over very long distances (heavily fat-tailed Cauchy kernel). Our analysis highlights two contrasting situations, where I is either mainly sensitive to community size (related to ecological drift) under the heavily fat-tailed kernel or mainly sensitive to dispersal distance under the thin-tailed kernel. We review dispersal distances of rainforest trees from field studies and assess the consistency between published estimates of I based on spatially-implicit models and the predictions of the kernel-based model in tropical forest plots. Most estimates of I were derived from large plots (10–50 ha) and were too large to be accounted for by a Cauchy kernel. Conversely, a fraction of the estimates based on multiple smaller plots (1 ha) appeared too small to be consistent with reported ranges of dispersal distances in tropical forests. Very large estimates may reflect within-plot habitat heterogeneity or estimation problems, while the smallest estimates likely imply other factors inhibiting migration beyond dispersal limitation. Our study underscores the need for interpreting I as an integrative index of migration limitation which, besides the limited seed dispersal, possibly includes habitat filtering or fragmentation.     

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.     

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.     

Limited seed dispersal shapes fine-scale spatial genetic structure in a Neotropical dioecious large-seeded palm

Seed and pollen dispersal contribute to gene flow and shape the genetic patterns of plants over fine spatial scales. We inferred fine-scale spatial genetic structure (FSGS) and estimated realized dispersal distances in Phytelephas aequatorialis, a Neotropical dioecious large-seeded palm. We aimed to explore how seed and pollen dispersal shape this genetic pattern in a focal population. For this purpose, we genotyped 138 seedlings and 99 adults with 20 newly developed microsatellite markers. We tested if rodent-mediated seed dispersal has a stronger influence than insect-mediated pollen dispersal in shaping FSGS. We also tested if pollen dispersal was influenced by the density of male palms around mother palms in order to further explore this ecological process in large-seeded plants. Rodent-mediated dispersal of these large seeds occurred mostly over short distances (mean 34.76 ± 34.06 m) while pollen dispersal distances were two times higher (mean 67.91 ± 38.29 m). The spatial extent of FSGS up to 35 m and the fact that seed dispersal did not increase the distance at which male alleles disperse suggest that spatially limited seed dispersal is the main factor shaping FSGS and contributes only marginally to gene flow within the population. Pollen dispersal distances depended on the density of male palms, decreasing when individuals show a clumped distribution and increasing when they are scattered. Our results show that limited seed dispersal mediated by rodents shapes FSGS in P. aequatorialis, while more extensive pollen dispersal accounts for a larger contribution to gene flow and may maintain high genetic diversity. Abstract in Spanish is available with online material.     

Assessing the risk of gene flow from genetically modified trees carrying mitigation transgenes

The increasing use of genetically modified (GM) plants has raised concerns about the escape of transgenes to conventional populations. To counteract possible fitness advantages of GM plants, an advantageous ‘primary transgene’ may be linked to a ‘mitigation’ transgene that is selectively disfavored. A risk related to this technique is the possible break-up of the transgenes. This may lead to the establishment of genotypes that only express the primary transgenic trait and are therefore selectively favored over both conventional and GM genotypes. We study here how the establishment of break-up genotypes in tree populations depends on life history, pollen dispersal kernels, and the linkage between transgenes. To this end, we model the dynamics of an initially homozygous GM population and a conventional population growing next to each other. The simulations are based on a spatially explicit simulation model for population-genotype dynamics. Pollen and seed dispersal are described by fat-tailed dispersal kernels. Break-up of transgenes is due subsequent crossing of GM and conventional genotypes. The spatial pattern of pollen and seed dispersal causes break-up and establishment of new invasive genotypes to occur mostly in the conventional population. The amount of individuals of these genotypes established in the conventional population depends on the interaction between life-history traits, pollen dispersal distances and the linkage of the transgenes. Thus, risk assessment of transgene break-up requires information on the genetics, local dynamics, and dispersal of GM and conventional varieties. The approach presented here provides a tool for combining these different sources of information, and it contributes to the formulation of general rules for the management of GM tree populations.     

Estimating seed and pollen movement in a monoecious plant: a hierarchical Bayesian approach integrating genetic and ecological data

The scale of seed and pollen movement in plants has a critical influence on population dynamics and interspecific interactions, as well as on their capacity to respond to environmental change through migration or local adaptation. However, dispersal can be challenging to quantify. Here, we present a Bayesian model that integrates genetic and ecological data to simultaneously estimate effective seed and pollen dispersal parameters and the parentage of sampled seedlings. This model is the first developed for monoecious plants that accounts for genotyping error and treats dispersal from within and beyond a plot in a fully consistent manner. The flexible Bayesian framework allows the incorporation of a variety of ecological variables, including individual variation in seed production, as well as multiple sources of uncertainty. We illustrate the method using data from a mixed population of red oak (Quercus rubra, Q. velutina, Q. falcata) in the NC piedmont. For simulated test data sets, the model successfully recovered the simulated dispersal parameters and pedigrees. Pollen dispersal in the example population was extensive, with an average father-mother distance of 178 m. Estimated seed dispersal distances at the piedmont site were substantially longer than previous estimates based on seed-trap data (average 128 m vs. 9.3 m), suggesting that, under some circumstances, oaks may be less dispersal-limited than is commonly thought, with a greater potential for range shifts in response to climate change.     

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.     

Pollen dispersal of tropical trees (Dinizia excelsa: Fabaceae) by native insects and African honeybees in pristine and fragmented Amazonian rainforest

Tropical rainforest trees typically occur in low population densities and rely on animals for cross-pollination. It is of conservation interest therefore to understand how rainforest fragmentation may alter the pollination and breeding structure of remnant trees. Previous studies of the Amazonian tree Dinizia excelsa (Fabaceae) found African honeybees (Apis mellifera scutellata) as the predominant pollinators of trees in highly disturbed habitats, transporting pollen up to 3.2 km between pasture trees. Here, using microsatellite genotypes of seed arrays, we compare outcrossing rates and pollen dispersal distances of (i) remnant D. excelsa in three large ranches, and (ii) a population in undisturbed forest in which African honeybees were absent. Self-fertilization was more frequent in the disturbed habitats (14%, n = 277 seeds from 12 mothers) than in undisturbed forest (10%, n = 295 seeds from 13 mothers). Pollen dispersal was extensive in all three ranches compared to undisturbed forest, however. Using a twogener analysis, we estimated a mean pollen dispersal distance of 1509 m in Colosso ranch, assuming an exponential dispersal function, and 212 m in undisturbed forest. The low effective density of D. excelsa in undisturbed forest (approximately 0.1 trees/ha) indicates that large areas of rainforest must be preserved to maintain minimum viable populations. Our results also suggest, however, that in highly disturbed habitats Apis mellifera may expand genetic neighbourhood areas, thereby linking fragmented and continuous forest populations.     

Long distance pollen-mediated gene flow at a landscape level: the weed beet as a case study

Gene flow is a crucial parameter that can affect the organization of genetic diversity in plant species. It has important implications in terms of conservation of genetic resources and of gene exchanges between crop to wild relatives and within crop species complex. In the Beta vulgaris complex, hybridization between crop and wild beets in seed production areas is well documented and the role of the ensuing hybrids, weed beets, as bridges towards wild forms in sugar beet production areas have been shown. Indeed, in contrast to cultivated beets that are bi-annual, weed beets can bolt, flower and reproduce in the same crop season. Nonetheless, the extent of pollen gene dispersal through weedy lineages remains unknown. In this study, the focus is directed towards weed-to-weed gene flow, and we report the results of a pollen-dispersal analysis within an agricultural landscape composed of five sugar beet fields with different levels of infestation by weed beets. Our results, based on paternity analysis of 3240 progenies from 135 maternal plants using 10 microsatellite loci, clearly demonstrate that even if weedy plants are mostly pollinated by individuals from the same field, some mating events occur between weed beets situated several kilometres apart (up to 9.6 km), with rates of interfield-detected paternities ranging from 11.3% to 17.5%. Moreover, we show that pollen flow appears to be more restricted when individuals are aggregated as most mating events occurred only for short-distance classes. The best-fit dispersal curves were fat-tailed geometric functions for populations exhibiting low densities of weed beets and thin-tailed Weibull function for fields with weed beet high densities. Thus, weed beet populations characterized by low density with geographically isolated individuals may be difficult to detect but are likely to act as pollen traps for pollen emitted by close and remote fields. Hence, it appears evident that interfield pollen-mediated gene flow between weed beets is almost unavoidable and could contribute to the diffusion of (trans)genes in the agricultural landscape.