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.     

HOW DOES POLLEN VERSUS SEED DISPERSAL AFFECT NICHE EVOLUTION?

In heterogeneous landscapes, the genetic and demographic consequences of dispersal influence the evolution of niche width. Unless pollen is limiting, pollen dispersal does not contribute directly to population growth. However, by disrupting local adaptation, it indirectly affects population dynamics. We compare the effect of pollen versus seed dispersal on the evolution of niche width in heterogeneous habitats, explicitly considering the feedback between maladaptation and demography. We consider two scenarios: the secondary contact of two subpopulations, in distinct, formerly isolated habitats, and the colonization of an empty habitat with dispersal between the new and ancestral habitat. With an analytical model, we identify critical levels of genetic variance leading to niche contraction (secondary contact scenario), or expansion (new habitat scenario). We confront these predictions with simulations where the genetic variance freely evolves. Niche contraction occurs when habitats are very different. It is faster as total gene flow increases or as pollen predominates in overall gene flow. Niche expansion occurs when habitat heterogeneity is not too high. Seed dispersal accelerates it, whereas pollen dispersal tends to retard it. In both scenarios very high seed dispersal leads to extinction. Overall, our results predict a wider niche for species dispersing seeds more than pollen.     

Contrasting patterns of pollen and seed flow influence the spatial genetic structure of sweet vernal grass (Anthoxanthum odoratum) populations

The spatial genetic structure of plant populations is determined by a combination of gene flow, genetic drift, and natural selection. Gene flow in most plants can result from either seed or pollen dispersal, but detailed investigations of pollen and seed flow among populations that have diverged following local adaptation are lacking. In this study, we compared pollen and seed flow among 10 populations of sweet vernal grass (Anthoxanthum odoratum) on the Park Grass Experiment. Overall, estimates of genetic differentiation that were based on chloroplast DNA (cpDNA) and, which therefore resulted primarily from seed flow, were lower (average F(ST) = 0.058) than previously published estimates that were based on nuclear DNA (average F(ST) = 0.095). Unlike nuclear DNA, cpDNA showed no pattern of isolation by adaptation; cpDNA differentiation was, however, inversely correlated with the number of additions (nutrients and lime) that each plot had received. We suggest that natural selection is restricting pollen flow among plots, whereas nutrient additions are increasing seed flow and genetic diversity by facilitating the successful germination and growth of immigrant seeds. This study highlights the importance of considering all potential gene flow mechanisms when investigating determinants of spatial genetic structure, and cautions against the widespread assumption that pollen flow is more important than seed flow for population connectivity in wind-pollinated species.     

Realized gene flow within mixed stands of Quercus robur L. and Q. petraea (Matt.) L. revealed at the stage of naturally established seedling

The estimates of contemporary gene flow assessed based on naturally established seedlings provide information much needed for understanding the abilities of forest tree populations to persist under global changes through migration and/or adaptation facilitated by gene exchange among populations. Here, we investigated pollen‐ and seed‐mediated gene flow in two mixed‐oak forest stands (consisting of Quercus robur L. and Q. petraea [Matt.] Liebl.). The gene flow parameters were estimated based on microsatellite multilocus genotypes of seedlings and adults and their spatial locations within the sample plots using models that attempt to reconstruct the genealogy of the seedling cohorts. Pollen and seed dispersal were modelled using the standard seedling neighbourhood model and a modification—the 2‐component seedling neighbourhood model, with the later allowing separation of the dispersal process into local and long‐distance components. The 2‐component model fitted the data substantially better than the standard model and provided estimates of mean seed and pollen dispersal distances accounting for long‐distance propagule dispersal. The mean distance of effective pollen dispersal was found to be 298 and 463 m, depending on the stand, while the mean distance of effective seed dispersal was only 8.8 and 15.6 m, which is consistent with wind pollination and primarily seed dispersal by gravity in Quercus. Some differences observed between the two stands could be attributed to the differences in the stand structure of the adult populations and the existing understory vegetation. Such a mixture of relatively limited seed dispersal with occasional long distance gene flow seems to be an efficient strategy for colonizing new habitats with subsequent local adaptation, while maintaining genetic diversity within populations.     

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.     

Nonrandom larval dispersal can steepen marine clines

Sharp and stable clinal variation is enigmatic when found in species with high gene flow. Classical population genetic models treat gene flow as a random homogenizing force countering local adaptation across habitat discontinuities. Under this view, dispersal over large spatial scales will lower the effectiveness of adaptation by natural selection at finer spatial scales. Thus, random gene flow will create a shallow phenotypic cline across an ecotone in response to a steep selection gradient. In sedentary marine species that disperse primarily as larvae, nonrandom dispersal patterns are expected due to coastal hydrodynamics. Surprisingly sharp phenotypic and genotypic clines have been documented in marine species with high gene flow. We are interested in the extent to which nonrandom dispersal could accentuate such clines. We model a linear species range in which populations have stable and uniform densities along a selection gradient; in contrast to random dispersal, convergent advection of larvae can amplify phenotypic differentiation if coupled with a semipermeable dispersal barrier in the convergence zone. The migration load caused by directional dispersal pushes the phenotypic mean away from the local trait optimum in downstream populations, that is, near the convergence zone. A dispersal barrier is possible as a result of colliding currents if the water and larvae are mostly displaced offshore, away from suitable settlement habitat. Disjunctions in a quantitative trait were enlarged in the convergence zone by faster current flows or a more complete dispersal barrier. With advection of larvae per generation one-third as far as the average dispersal distance by diffusion, convergence on a dispersal barrier with 40% permeability generated a trait disjunction across the convergence zone of two phenotypic standard deviations. Without directional dispersal, similar clines also developed across a habitat gap, where population density was low, or across dispersal barriers with less than 1% permeability. These findings suggest that the types of hydrographic phenomena often associated with marine transition zones can strongly affect the balance between gene flow and selection and generate surprisingly steep clines given the large-scale gene flow expected from larvae.     

Multiscale spatial genetic structure within and between populations of wild cherry trees in nuclear genotypes and chloroplast haplotypes

Spatial genetic structure (SGS) of plants mainly depends on the effective population size and gene dispersal. Maternally inherited loci are expected to have higher genetic differentiation between populations and more intensive SGS within populations than biparentally inherited loci because of smaller effective population sizes and fewer opportunities of gene dispersal in the maternally inherited loci. We investigated biparentally inherited nuclear genotypes and maternally inherited chloroplast haplotypes of microsatellites in 17 tree populations of three wild cherry species under different conditions of tree distribution and seed dispersal. As expected, interpopulation genetic differentiation was 6–9 times higher in chloroplast haplotypes than in nuclear genotypes. This difference indicated that pollen flow 4–7 times exceeded seed flow between populations. However, no difference between nuclear and chloroplast loci was detected in within‐population SGS intensity due to their substantial variation among the populations. The SGS intensity tended to increase as trees became more aggregated, suggesting that tree aggregation biased pollen and seed dispersal distances toward shorter. The loss of effective seed dispersers, Asian black bears, did not affect the SGS intensity probably because of mitigation of the bear loss by other vertebrate dispersers and too few tree generations after the bear loss to alter SGS. The findings suggest that SGS is more variable in smaller spatial scales due to various ecological factors in local populations.     

Selective interactions between short-distance pollen and seed dispersal in self-compatible species

In plants, genes may disperse through both pollen and seeds. Here we provide a first theoretical study of the mechanisms and consequences of the joint evolution of pollen and seed dispersal. We focus on hermaphroditic self-compatible species distributed in structured populations, assuming island dispersal of pollen and seeds among small patches of plants within large populations. Three traits are studied: the rate of among-patch seed dispersal, the rate of among-patch pollen dispersal, and the rate of within-patch pollen movement. We first analytically derive the evolutionary equilibrium state of each trait, dissect the pairwise selective interactions, and describe the joint three-trait evolutionary equilibrium under the cost of dispersal and kin competition. These results are then analytically and numerically extended to the case when selfed seeds suffer from depressed competitiveness (inbreeding depression, no heterosis). Finally individual-based simulations are used to account for a more realistic model of inbreeding load. Pollen movement is shown to generate opposite selection pressures on seed dispersal depending on spatial scale: within-patch pollen movement favors seed dispersal, whereas among-patch pollen dispersal inhibits seed dispersal. Seed dispersal selects for short-distance movements of pollen and it selects against long-distance dispersal. These interactions shape the joint evolution of these traits. Kin competition favors among-patch seed dispersal over among-patch pollen dispersal for low costs of within-patch pollen movement (and vice versa for significant costs of within-patch pollen movement). Inbreeding depression favors allogamy through high rates of within- and among-patch pollen movement. Surprisingly, it may select either for or against seed dispersal depending on the cost of among-patch pollen dispersal. Heterosis favors increased among-patch dispersal through pollen and seeds. But because these two stages inhibit each other, their joint evolution might lead to decreased seed dispersal in the presence of heterosis. Of crucial importance are the costs of dispersal.     

Inbreeding depression and heterosis in a subdivided population: influence of the mating system

We investigate the joint effects of gene flow and selfing on the level of inbreeding depression, heterosis and genetic load in a subdivided population at equilibrium. Low gene flow reduces inbreeding depression and substantially increases heterosis. However, in highly self-fertilizing populations, inbreeding depression is independent of the amount of gene flow. When migration occurs via pollen, consanguinity of the reproductive system could have a negative influence on subpopulation persistence, in contrast to the case of isolated populations. However, with only seed migration, genetic load and heterosis depend mildly on the mating system. From an evolutionary point of view, we reach two main conclusions: first, outcrossing is selected for if gene flow is low; second, intermediate levels of gene flow could promote mixed mating systems, especially when migration occurs through pollen.     

Adaptive gradients and isolation-by-distance with postglacial migration in Picea sitchensis

Fossil pollen records suggest rapid migration of tree species in response to Quaternary climate warming. Long-distance dispersal and high gene flow would facilitate rapid migration, but would initially homogenize variation among populations. However, contemporary clinal variation in adaptive traits along environmental gradients shown in many tree species suggests that local adaptation can occur during rapid migration over just a few generations in interglacial periods. In this study, we compared growth performance and pollen genetic structure among populations to investigate how populations of Sitka spruce (Picea sitchensis) have responded to local selection along the historical migration route. The results suggest strong adaptive divergence among populations (average Q(ST)=0.61), corresponding to climatic gradients. The population genetic structure, determined by microsatellite markers (R(ST)=0.09; F(ST)=0.11), was higher than previous estimates from less polymorphic genetic markers. The significant correlation between geographic and pollen haplotype genetic (R(ST)) distances (r=0.73, P<0.01) indicates that the current genetic structure has been shaped by isolation-by-distance, and has developed in relatively few generations. This suggests relatively limited gene flow among populations on a recent timescale. Gene flow from neighboring populations may have provided genetic diversity to founder populations during rapid migration in the early stages of range expansion. Increased genetic diversity subsequently enhanced the efficiency of local selection, limiting gene flow primarily to among similar environments and facilitating the evolution of adaptive clinal variation along environmental gradients.     

Effect of gene flow on spatial genetic structure in the riparian canopy tree Cercidiphyllum japonicum revealed by microsatellite analysis

Few studies have analyzed pollen and seed movements at local scale, and genetic differentiation among populations covering the geographic distribution range of a species. We carried out such a study on Cercidiphyllum japonicum; a dioecious broad-leaved tree of cool-temperate riparian forest in Japan. We made direct measurement of pollen and seed movements in a site, genetic structure at the local scale, and genetic differentiation between populations covering the Japanese Archipelago. Parentage analysis of seedlings within a 20-ha study site indicated that at least 28.8% of seedlings were fertilized by pollen from trees outside the study site. The average pollination distance within the study site was 129 m, with a maximum of 666 m. The genotypes of 30% of seedlings were incompatible with those of the nearest female tree, and the maximum seed dispersal distance within the study site was over 300 m. Thus, long-distance gene dispersal is common in this species. The correlation between genetic relatedness and spatial distance among adult trees within the population was not significant, indicating an absence of fine-scale genetic structure perhaps caused by high levels of pollen flow and overlapping seed shadows. Six populations sampled throughout the distribution of C. japonicum in Japan showed significant isolation-by-distance but low levels of genetic differentiation (F(ST) = 0.043), also indicating long-distance gene flow in C. japonicum. Long-distance gene flow had a strong influence on the genetic structure at different spatial scales, and contributes to the maintenance of genetic diversity in C. japonicum.     

Conceptual issues in local adaptation

Studies of local adaptation provide important insights into the power of natural selection relative to gene flow and other evolutionary forces. They are a paradigm for testing evolutionary hypotheses about traits favoured by particular environmental factors. This paper is an attempt to summarize the conceptual framework for local adaptation studies. We first review theoretical work relevant for local adaptation. Then we discuss reciprocal transplant and common garden experiments designed to detect local adaptation in the pattern of deme × habitat interaction for fitness. Finally, we review research questions and approaches to studying the processes of local adaptation – divergent natural selection, dispersal and gene flow, and other processes affecting adaptive differentiation of local demes. We advocate multifaceted approaches to the study of local adaptation, and stress the need for experiments explicitly addressing hypotheses about the role of particular ecological and genetic factors that promote or hinder local adaptation. Experimental evolution of replicated populations in controlled spatially heterogeneous environments allow direct tests of such hypotheses, and thus would be a valuable way to complement research on natural populations.     

Population structure of an endangered species living in contrasted habitats: Parnassia palustris (Saxifragaceae)

In endangered species, it is critical to analyse the level at which populations interact (i.e. dispersal) as well as the levels of inbreeding and local adaptation to set up conservation policies. These parameters were investigated in the endangered species Parnassia palustris living in contrasted habitats. We analysed population structure in 14 populations of northern France for isozymes, cpDNA markers and phenotypic traits related to fitness. Within population genetic diversity and inbreeding coefficients were not correlated to population size. Populations seem not to have undergone severe recent bottleneck. Conversely to pollen migration, seed migration seems limited at a regional scale, which could prevent colonization of new sites even if suitable habitats appear. Finally, the habitat type affects neither within-population genetic diversity nor genetic and phenotypic differentiation among populations. Thus, even if unnoticed local adaptation to habitats exists, it does not influence gene flow between populations.     

LOCAL ADAPTATION WHEN COMPETITION DEPENDS ON PHENOTYPIC SIMILARITY

Recent work incorporating demographic–genetic interactions indicates the importance of population size, gene flow, and selection in influencing local adaptation. This work typically assumes that density‐dependent survival affects individuals equally, but individuals in natural population rarely compete equally. Among‐individual differences in resource use generate stronger competition between more similar phenotypes (frequency‐dependent competition) but it remains unclear how this additional form of selection changes the interactions between population size, gene flow, and local stabilizing selection. Here, we integrate migration–selection dynamics with frequency‐dependent competition. We developed a coupled demographic‐quantitative genetic model consisting of two patches connected by dispersal and subject to local stabilizing selection and competition. Our model shows that frequency‐dependent competition slightly increases local adaptation, greatly increases genetic variance within patches, and reduces the amount that migration depresses population size, despite the increased genetic variance load. The effects of frequency‐dependence depend on the strength of divergent selection, trait heritability, and when mortality occurs in the life cycle in relation to migration and reproduction. Essentially, frequency‐dependent competition reduces the density‐dependent interactions between migrants and residents, the extent to which depends on how different and common immigrants are compared to residents. Our results add new dynamics that illustrate how competition can alter the effects of gene flow and divergent selection on local adaptation and population carrying capacities.     

The role of gene flow in shaping genetic structures of the subtropical conifer species Araucaria angustifolia

The morphological features of pollen and seed of Araucaria angustifolia have led to the proposal of limited gene dispersal for this species. We used nuclear microsatellite and AFLP markers to assess patterns of genetic variation in six natural populations at the intra- and inter-population level, and related our findings to gene dispersal in this species. Estimates of both fine-scale spatial genetic structure (SGS) and migration rate suggest relatively short-distance gene dispersal. However, gene dispersal differed among populations, and effects of more efficient dispersal within population were observed in at least one stand. In addition, even though some seed dispersal may be aggregated in this principally barochorous species, reasonable secondary seed dispersal, presumably facilitated by animals, and overlap of seed shadows within populations is suggested. Overall, no correlation was observed between levels of SGS and inbreeding, density or age structure, except that a higher level of SGS was revealed for the population with a higher number of juvenile individuals. A low estimate for the number of migrants per generation between two neighbouring populations implies limited gene flow. We expect that stepping-stone pollen flow may have contributed to low genetic differentiation among populations observed in a previous survey. Thus, strategies for maintenance of gene flow among remnant populations should be considered in order to avoid degrading effects of population fragmentation on the evolution of A. angustifolia.     

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.     

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.     

Fragmentation reduces regional‐scale spatial genetic structure in a wind‐pollinated tree because genetic barriers are removed

Gene flow strongly influences the regional genetic structuring of plant populations. Seed and pollen dispersal patterns can respond differently to the increased isolation resulting from habitat fragmentation, with unpredictable consequences for gene flow and population structuring. In a recently fragmented landscape we compared the pre‐ and post‐fragmentation genetic structure of populations of a tree species where pollen and seed dispersal respond differentially to forest fragmentation generated by flooding. Castanopsis sclerophylla is wind‐pollinated, with seeds that are dispersed by gravity and rodents. Using microsatellites, we found no significant difference in genetic diversity between pre‐ and post‐fragmentation cohorts. Significant genetic structure was observed in pre‐fragmentation cohorts, due to an unknown genetic barrier that had isolated one small population. Among post‐fragmentation cohorts this genetic barrier had disappeared and genetic structure was significantly weakened. The strengths of genetic structuring were at a similar level in both cohorts, suggesting that overall gene flow of C. sclerophylla has been unchanged by fragmentation at the regional scale. Fragmentation has blocked seed dispersal among habitats, but this appears to have been compensated for by enhanced pollen dispersal, as indicated by the disappearance of a genetic barrier, probably as a result of increased wind speeds and easier pollen movement over water. Extensive pollen flow can counteract some negative effects of fragmentation and assist the long‐term persistence of small remnant populations.     

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.