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

Wild snapdragon plant pedigree sheds light on limited connectivity enhanced by higher migrant reproductive success in a fragmented landscape

Background: In contrast with historical knowledge, a recent view posits that a non-negligible proportion of populations thrive in a fragmented landscape. One underlying mechanism is the maintenance of functional connectivity, i.e., the net flow of individuals or their genes moving among suitable habitat patches. Alternatively, functional connectivity might be typically limited but enhanced by a higher reproductive success of migrants. Methods: We tested for this hypothesis in wild snapdragon plants inhabiting six patches separated by seawater in a fragmented Mediterranean scrubland landscape. We reconstructed their pedigree by using a parentage assignment method based on microsatellite genetic markers. We then estimated functional connectivity and the reproductive success of plants resulting from between-patch dispersal events. Results: We found that wild snapdragon plants thrived in this fragmented landscape, although functional connectivity between habitat patches was low (i.e. 2.9%). The progeny resulting from between-patch dispersal events had a higher reproductive success than residents. Conclusion: Our findings imply that low functional connectivity in a fragmented landscapes may have been enhanced by higher reproductive success after migration. This original mechanisms might be partly compensating the negative impact of fragmentation.

Landscape barriers to pollen and seed flow in the dioecious tropical tree Astronium fraxinifolium in Brazilian savannah

Gene flow studies provide information on gene exchange between populations, which is essential for developing genetic conservation strategies. Such analyses enable a better understanding of the life history and seed and pollen dispersal mechanisms of plant species. In this study, we investigate pollen and seed flow in a regenerant population of the pioneer species Astronium fraxinifolium in an area degraded during the construction of a hydroelectric dam. We mapped, sampled, sexed, and genotyped 386 individuals in the regenerant population (RP), as well as 128 adult trees located along two highways adjacent to the degraded area; one in Mato Grosso do Sul State (MS) and other in São Paulo State (SP). Parentage analyses was carried out for 370 individuals of the RP population, using as putative parents 348 individuals from RP and all 128 individuals sampled in MS and SP. Based on parentage analysis and eight microsatellite loci, our analyses revealed that for individuals of the RP with an identified father (pollen donor), 1.1% of the pollen was dispersed up to 532 m, while for those with an identified mother (seed donor), 0.5% of seeds were dispersed up to 4,782 m. However, a large proportion of pollen (76.5%) and seeds (57%) immigrated from trees outside the sampled populations. Pollen and seeds were dispersed through a pattern of isolation by distance. Genetic diversity was significantly similar between adults of both highway populations and individuals from RP, with significant levels of inbreeding detected only in RP. Our results demonstrate that the nearest trees contributed pollen and seeds for the recovery of the degraded area, indicating reproductive spatial isolation among the sampled populations due to the damming of the river. Such results help to understand the process of regeneration for A. fraxinifolium in regenerant populations to inform strategies for conservation and environmental recovery with this species.

Crown defoliation decreases reproduction and wood growth in a marginal European beech population

BACKGROUND AND AIMS

Abiotic and biotic stresses related to climate change have been associated with increased crown defoliation, decreased growth and a higher risk of mortality in many forest tree species, but the impact of stresses on tree reproduction and forest regeneration remains understudied. At the dry, warm margin of species distributions, flowering, pollination and seed maturation are expected to be affected by drought, late frost and other stresses, eventually resulting in reproduction failure. Moreover, inter-individual variation in reproductive performance versus other performance traits (growth, survival) could have important consequences for population dynamics. This study investigated the relationships among individual crown defoliation, growth and reproduction in a drought-prone population of European beech, Fagus sylvatica.

METHODS

We used a spatially explicit mating model and marker-based parentage analyses to estimate effective female and male fecundities of 432 reproductive trees, which were also monitored for basal area increment and crown defoliation over 9 years.

KEY RESULTS

Female and male fecundities varied markedly between individuals, more than did growth. Both female fecundity and growth decreased with increasing crown defoliation and competition, and increased with size. Moreover, the negative effect of defoliation on female fecundity was size-dependent, with a slower decline in female fecundity with increasing defoliation for the large individuals. Finally, a trade-off between growth and female fecundity was observed in response to defoliation: some large trees maintained significant female fecundity at the expense of reduced growth in response to defoliation, while some other defoliated trees maintained high growth at the expense of reduced female fecundity.

CONCLUSIONS

Our results suggest that, while decreasing their growth, some large defoliated trees still contribute to reproduction through seed production and pollination. This non-coordinated decline of growth and fecundity at individual level in response to stress may compromise the evolution of stress-resistance traits at population level, and increase forest tree vulnerability.

At 50, Janzen–Connell Has Come of Age

Fifty years ago, Janzen (1970) and Connell (1971) independently published a revolutionary idea to explain the hyperdiverse tree communities of the tropics. The essential observations were that seedfall is concentrated in the vicinity of fruiting trees, whereas saplings recruit at a distance from reproductive conspecifics. These observations were encapsulated in a simple focal-tree model constructed of intersecting curves for seedfall and escape from host-specific enemies postulated to attack propagules (seeds and seedlings) in the vicinity of reproductive conspecifics. In conflict with the thinking of the times, the mechanism operates from the top down rather than from the bottom up. A deterrent to broad acceptance has been the giant intuitive leap required to generalize the focal tree model to an entire forest community. Recent theoretical and empirical results have succeeded in bridging the gap between the focal tree model and its community-level implications. With these new findings, Janzen–Connell has come of age.

The total dispersal kernel: a review and future directions

The distribution and abundance of plants across the world depends in part on their ability to move, which is commonly characterized by a dispersal kernel. For seeds, the total dispersal kernel (TDK) describes the combined influence of all primary, secondary and higher-order dispersal vectors on the overall dispersal kernel for a plant individual, population, species or community. Understanding the role of each vector within the TDK, and their combined influence on the TDK, is critically important for being able to predict plant responses to a changing biotic or abiotic environment. In addition, fully characterizing the TDK by including all vectors may affect predictions of population spread. Here, we review existing research on the TDK and discuss advances in empirical, conceptual modelling and statistical approaches that will facilitate broader application. The concept is simple, but few examples of well-characterized TDKs exist. We find that significant empirical challenges exist, as many studies do not account for all dispersal vectors (e.g. gravity, higher-order dispersal vectors), inadequately measure or estimate long-distance dispersal resulting from multiple vectors and/or neglect spatial heterogeneity and context dependence. Existing mathematical and conceptual modelling approaches and statistical methods allow fitting individual dispersal kernels and combining them to form a TDK; these will perform best if robust prior information is available. We recommend a modelling cycle to parameterize TDKs, where empirical data inform models, which in turn inform additional data collection. Finally, we recommend that the TDK concept be extended to account for not only where seeds land, but also how that location affects the likelihood of establishing and producing a reproductive adult, i.e. the total effective dispersal kernel.

Seed limitation in an Amazonian floodplain forest

We monitored a close-spaced grid of 289 seed traps in 1.44 ha for 8.4 yr in an Amazonian floodplain forest. In a tree community containing hundreds of species, a median of just three to four species of tree seeds falls annually into each 0.5-m² establishment site. The number of seed species reaching a given site increased linearly with time for the duration of the monitoring period, indicating a roughly random arrival of seed species in a given site-year. The number of seed species captured each year over the entire grid ranged from one-third to one-half of the total captured over the 8.4 yr of monitoring, revealing a substantial temporal component of variation in the seed rain. Seed rain at the 0.5-m² scale displayed extreme spatial variability when all potentially viable seeds were tallied, whereas the rain of dispersed seeds was scant, more nearly uniform, and better mixed. Dispersal limitation, defined as failure of seeds to reach establishment sites, is ≥99% per year for a majority of species, explaining why seed augmentation experiments are often successful. Dispersal limitation has been evoked as an explanation for distance-dependent species turnover in tropical tree communities, but that interpretation contrasts with the fact that many Amazonian tree species possess large geographical ranges that extend for hundreds or thousands of kilometers. A better understanding of the processes that bridge the gap between the scales of seedling establishment and the regulation of forest composition will require new methodologies for studying dispersal on scales larger than those yet achieved.

The future of parentage analysis: From microsatellites to SNPs and beyond

Parentage analysis is a cornerstone of molecular ecology that has delivered fundamental insights into behaviour, ecology and evolution. Microsatellite markers have long been the king of parentage, their hypervariable nature conferring sufficient power to correctly assign offspring to parents. However, microsatellite markers have seen a sharp decline in use with the rise of next-generation sequencing technologies, especially in the study of population genetics and local adaptation. The time is ripe to review the current state of parentage analysis and see how it stands to be affected by the emergence of next-generation sequencing approaches. We find that single nucleotide polymorphisms (SNPs), the typical next-generation sequencing marker, remain underutilized in parentage analysis but are gaining momentum, with 58 SNP-based parentage analyses published thus far. Many of these papers, particularly the earlier ones, compare the power of SNPs and microsatellites in a parentage context. In virtually every case, SNPs are at least as powerful as microsatellite markers. As few as 100-500 SNPs are sufficient to resolve parentage completely in most situations. We also provide an overview of the analytical programs that are commonly used and compatible with SNP data. As the next-generation parentage enterprise grows, a reliance on likelihood and Bayesian approaches, as opposed to strict exclusion, will become increasingly important. We discuss some of the caveats surrounding the use of next-generation sequencing data for parentage analysis and conclude that the future is bright for this important realm of molecular ecology.

Fruit production is influenced by tree size and size-asymmetric crowding in a wet tropical forest

In tropical forest communities, seedling recruitment can be limited by the number of fruit produced by adults. Fruit production tends to be highly unequal among trees of the same species, which may be due to environmental factors. We observed fruit production for ~2,000 trees of 17 species across 3 years in a wet tropical forest in Costa Rica. Fruit production was modeled as a function of tree size, nutrient availability, and neighborhood crowding. Following model selection, tree size and neighborhood crowding predicted both the probability of reproduction and the number of fruit produced. Nutrient availability only predicted only the probability of reproduction. In all species, larger trees were more likely to be reproductive and produce more fruit. In addition, number of fruit was strongly negatively related to presence of larger neighboring trees in 13 species; presence of all neighboring trees had a weak-to-moderate negative influence on reproductive status in 16 species. Among various metrics of soil nutrient availability, only sum of base cations was positively associated with reproductive status, and for only four species. Synthesis Overall, these results suggest that direct influences on fruit production tend to be mediated through tree size and crowding from neighboring trees, rather than soil nutrients. However, we found variation in the effects of neighbors and nutrients among species; mechanistic studies of allocation to fruit production are needed to explain these differences.

Joining and splitting models with Markov melding

Analysing multiple evidence sources is often feasible only via a modular approach, with separate submodels specified for smaller components of the available evidence. Here we introduce a generic framework that enables fully Bayesian analysis in this setting. We propose a generic method for forming a suitable joint model when joining submodels, and a convenient computational algorithm for fitting this joint model in stages, rather than as a single, monolithic model. The approach also enables splitting of large joint models into smaller submodels, allowing inference for the original joint model to be conducted via our multi-stage algorithm. We motivate and demonstrate our approach through two examples: joining components of an evidence synthesis of A/H1N1 influenza, and splitting a large ecology model.

Estimating seed and pollen dispersal kernels from genetic data demonstrates a high pollen dispersal capacity for an endangered palm species

PREMISE OF THE STUDY

Seed and pollen dispersal are key processes shaping plant population dynamics and maintaining genetic diversity. The essence of these processes is the movement of propagules from a parental tree to the site of propagule establishment. The estimation of plant dispersal kernels has remained challenging due to the difficulty of making direct observations. We estimated the dispersal capacity of the endangered palm Euterpe edulis, whose seeds are dispersed by vertebrates and pollen by insects.

METHODS

We used a hierarchical Bayesian model with genetic data from reproductive plants, juveniles, and embryos to estimate dispersal kernels. Our analyses account for genotyping error and uncertainty in parental assignment.

KEY RESULTS

We found that seeds were dispersed at most a few hundred meters, but pollen was dispersed up to several kilometers. We hypothesize that this long-distance pollen dispersal is generated mainly by euglossine bees, whereas the main dispersal vectors for short-distance seed dispersal are thrushes. The long-distance dispersal of pollen suggests a high level of gene flow that should maintain genetic diversity of E. edulis. Despite the relation between long-distance dispersal and genetic diversity, we observed low genetic diversity and inbreeding within the local population, which are probably due to restricted gene flow due to the low density of this population and its aggregated spatial distribution.

CONCLUSIONS

We conclude that if conservation actions are able to restore the population density of E. edulis, the recovery of its genetic diversity will be facilitated because of its high dispersal capacity, especially with regard to pollen.

Seed and pollen gene dispersal in Taxus baccata, a dioecious conifer in the face of strong population fragmentation

Background and Aims

Dispersal is crucial due to its direct impact on dynamics of a species' distribution as well as having a role in shaping adaptive potential through gene flow. In plants forming scarce and small populations, knowledge about the dispersal process is required to assess the potential for colonizing new habitats and connectivity of present and future populations. This study aimed to assess dispersal potential in Taxus baccata, a dioecious gymnosperm tree with a wide but highly fragmented distribution.

Methods

Seed and pollen dispersal kernels were estimated directly in the framework of the spatially explicit mating model, where genealogies of naturally established seedlings were reconstructed with the help of microsatellite markers. In this way, six differently shaped dispersal functions were compared.

Key Results

Seed dispersal followed a leptokurtic distribution, with the Exponential-Power, the Power-law and Weibull being almost equally best-fitting models. The pollen dispersal kernel appeared to be more fat-tailed than the seed dispersal kernel, and the Lognormal and the Exponential-Power function showed the best fit. The rate of seed immigration from the background sources was not significantly different from the rate of pollen immigration (13.1 % vs. 19.7 %) and immigration rates were in agreement with or below maximum predictions based on the estimated dispersal kernels. Based on the multimodel approach, 95 % of seeds travel <109 m, while 95 % of pollen travels <704 m from the source.

Conclusions

The results showed that, at a local spatial scale, yew seeds travel shorter distances than pollen, facilitating a rapid development of a kinship structure. At the landscape level, however, although yew exhibits some potential to colonize new habitats through seed dispersal, genetic connectivity between different yew remnants is strongly limited. Taking into account strong population fragmentation, the study suggests that gene dispersal may be a limiting factor of the adaptability of the species.

Spatial genetic structure within populations and management implications of the South American species Acacia aroma (Fabaceae)

The identification of factors that structure intraspecific diversity is of particular interest for biological conservation and restoration ecology. All rangelands in Argentina are currently experiencing some form of deterioration or desertification. Acacia aroma is a multipurpose species widely distributed throughout this country. In this study, we used the AFLP technique to study genetic diversity, population genetic structure, and fine-scale spatial genetic structure in 170 individuals belonging to 6 natural Argentinean populations. With 401 loci, the mean heterozygosity (H_E = 0.2) and the mean percentage of polymorphic loci (PPL = 62.1%) coefficients indicated that the genetic variation is relatively high in A. aroma. The analysis with STRUCTURE showed that the number of clusters (K) was 3. With Geneland analysis, the number of clusters was K = 4, sharing the same grouping as STRUCTURE but dividing one population into two groups. When studying SGS, significant structure was detected in 3 of 6 populations. The neighbourhood size in these populations ranged from 15.2 to 64.3 individuals. The estimated gene dispersal distance depended on the effective population density and disturbance level and ranged from 45 to 864 m. The combined results suggest that a sampling strategy, which aims to maintain a considerable part of the variability contained in natural populations sampled here, would include at least 3 units defined by the clusters analyses that exhibit particular genetic properties. Moreover, the current SGS analysis suggests that within the wider management units/provinces, seed collection from A. aroma should target trees separated by a minimum distance of 50 m but preferably 150 m to reduce genetic relatedness among seeds from different trees.

Inverse approach to estimating larval dispersal reveals limited population connectivity along 700 km of wave-swept open coast

Demographic connectivity is fundamental to the persistence and resilience of metapopulations, but our understanding of the link between reproduction and recruitment is notoriously poor in open-coast marine populations. We provide the first evidence of high local retention and limited connectivity among populations spanning 700 km along an open coast in an upwelling system. Using extensive field measurements of fecundity, population size and settlement in concert with a Bayesian inverse modelling approach, we estimated that, on average, Petrolisthes cinctipes larvae disperse only 6.9 km (±25.0 km s.d.) from natal populations, despite spending approximately six weeks in an open-coast system that was once assumed to be broadly dispersive. This estimate differed substantially from our prior dispersal estimate (153.9 km) based on currents and larval duration and behaviour, revealing the importance of employing demographic data in larval dispersal estimates. Based on this estimate, we predict that demographic connectivity occurs predominantly among neighbouring populations less than 30 km apart. Comprehensive studies of larval production, settlement and connectivity are needed to advance an understanding of the ecology and evolution of life in the sea as well as to conserve ecosystems. Our novel approach provides a tractable framework for addressing these questions for species occurring in discrete coastal populations.

Discovery of long-distance gamete dispersal in a lichen-forming ascomycete

Accurate estimates of gamete and offspring dispersal range are required for the understanding and prediction of spatial population dynamics and species persistence. Little is known about gamete dispersal in fungi, especially in lichen-forming ascomycetes. Here, we estimate the dispersal functions of clonal propagules, gametes and ascospores of the epiphytic lichen Lobaria pulmonaria. We use hierarchical Bayesian parentage analysis, which integrates genetic and ecological information from multiannual colonization and dispersal source data collected in a large, old-growth forest landscape. The effective dispersal range of gametes is several hundred metres to kilometres from potential paternal individuals. By contrast, clonal propagules disperse only tens of metres, and ascospores disperse over several thousand metres. Our study reveals the dispersal distances of individual reproductive units; clonal propagules, gametes and ascospores, which is of great importance for a thorough understanding of the spatial dynamics of ascomycetes. Sexual reproduction occurs between distant individuals. However, whereas gametes and ascospores disperse over long distances, the overall rate of colonization of trees is low. Hence, establishment is the limiting factor for the colonization of new host trees by the lichen in old-growth landscapes.

Simulating the Interacting Effects of Intraspecific Variation, Disturbance, and Competition on Climate-Driven Range Shifts in Trees

Climate change is expected to favor shifts in plant distributions; some such shifts are already being observed along elevation gradients. However, the rate of such shifts may be limited by their ability to reach newly suitable areas and by competition from resident species. The degree of local adaptation and genetic variation may also play a role in the interaction between migrants and residents by affecting relative fitness. We used a simulation model to explore the interacting effects of dispersal, fecundity, disturbance, and genetic variation on range-edge dynamics between a pair of demographically similar tree species. Ideal climate for an individual is determined by genotype. The simulated landscape undergoes an 80-year period of climate change in which climate bands shift upslope; subsequently, climate is held constant for 300 years. The presence of a high-elevation competitor caused a significant lag in the range shift of the low-elevation species relative to competition-free scenarios. Increases in fecundity and dispersal distance both helped to speed up the replacement of the high-elevation species by the low-elevation species at their range boundary. While some disturbance scenarios facilitated this transition, frequent canopy disturbance inhibited colonization by removing reproductive adults and led to range contractions in both species. Differences between dispersal scenarios were more pronounced when disturbance was frequent (15 vs. 25 year return interval) and dispersal was limited. When the high-elevation species lacked genetic variation, its range was more-easily invaded by the low-elevation species, while a similar lack of variation in the low-elevation species inhibited colonization-but only when this lack of variation decreased the fitness of the affected species near the range boundary. Our model results support the importance of measuring and including dispersal/fecundity, disturbance type and frequency, and genetic variation when assessing the potential for range shifts and species vulnerability to climate change.

Speeding up ecological and evolutionary computations in R; essentials of high performance computing for biologists

Computation has become a critical component of research in biology. A risk has emerged that computational and programming challenges may limit research scope, depth, and quality. We review various solutions to common computational efficiency problems in ecological and evolutionary research. Our review pulls together material that is currently scattered across many sources and emphasizes those techniques that are especially effective for typical ecological and environmental problems. We demonstrate how straightforward it can be to write efficient code and implement techniques such as profiling or parallel computing. We supply a newly developed R package (aprof) that helps to identify computational bottlenecks in R code and determine whether optimization can be effective. Our review is complemented by a practical set of examples and detailed Supporting Information material (S1–S3 Texts) that demonstrate large improvements in computational speed (ranging from 10.5 times to 14,000 times faster). By improving computational efficiency, biologists can feasibly solve more complex tasks, ask more ambitious questions, and include more sophisticated analyses in their research.

A last stand in the Po valley: genetic structure and gene flow patterns in Ulmus minor and U. pumila

BACKGROUND AND AIMS

Ulmus minor has been severely affected by Dutch elm disease (DED). The introduction into Europe of the exotic Ulmus pumila, highly tolerant to DED, has resulted in it widely replacing native U. minor populations. Morphological and genetic evidence of hybridization has been reported, and thus there is a need for assessment of interspecific gene flow patterns in natural populations. This work therefore aimed at studying pollen gene flow in a remnant U. minor stand surrounded by trees of both species scattered across an agricultural landscape.

METHODS

All trees from a small natural stand (350 in number) and the surrounding agricultural area within a 5-km radius (89) were genotyped at six microsatellite loci. Trees were morphologically characterized as U. minor, U. pumila or intermediate phenotypes, and morphological identification was compared with Bayesian clustering of genotypes. For paternity analysis, seeds were collected in two consecutive years from 20 and 28 mother trees. Maximum likelihood paternity assignment was used to elucidate intra- and interspecific gene flow patterns.

KEY RESULTS

Genetic structure analyses indicated the presence of two genetic clusters only partially matching the morphological identification. The paternity analysis results were consistent between the two consecutive years of sampling and showed high pollen immigration rates (∼0·80) and mean pollination distances (∼3 km), and a skewed distribution of reproductive success. Few intercluster pollinations and putative hybrid individuals were found.

CONCLUSIONS

Pollen gene flow is not impeded in the fragmented agricultural landscape investigated. High pollen immigration and extensive pollen dispersal distances are probably counteracting the potential loss of genetic variation caused by isolation. Some evidence was also found that U. minor and U. pumila can hybridize when in sympatry. Although hybridization might have beneficial effects on both species, remnant U. minor populations represent a valuable source of genetic diversity that needs to be preserved.

The gravity of pollination: integrating at-site features into spatial analysis of contemporary pollen movement

Pollen-mediated gene flow is a major driver of spatial genetic structure in plant populations. Both individual plant characteristics and site-specific features of the landscape can modify the perceived attractiveness of plants to their pollinators and thus play an important role in shaping spatial genetic variation. Most studies of landscape-level genetic connectivity in plants have focused on the effects of interindividual distance using spatial and increasingly ecological separation, yet have not incorporated individual plant characteristics or other at-site ecological variables. Using spatially explicit simulations, we first tested the extent to which the inclusion of at-site variables influencing local pollination success improved the statistical characterization of genetic connectivity based upon examination of pollen pool genetic structure. The addition of at-site characteristics provided better models than those that only considered interindividual spatial distance (e.g. IBD). Models parameterized using conditional genetic covariance (e.g. population graphs) also outperformed those assuming panmixia. In a natural population of Cornus florida L. (Cornaceae), we showed that the addition of at-site characteristics (clumping of primary canopy opening above each maternal tree and maternal tree floral output) provided significantly better models describing gene flow than models including only between-site spatial (IBD) and ecological (isolation by resistance) variables. Overall, our results show that including interindividual and local ecological variation greatly aids in characterizing landscape-level measures of contemporary gene flow.

Assessing intraspecific variation in effective dispersal along an altitudinal gradient: a test in two Mediterranean high-mountain plants

BACKGROUND

Plant recruitment depends among other factors on environmental conditions and their variation at different spatial scales. Characterizing dispersal in contrasting environments may thus be necessary to understand natural intraspecific variation in the processes underlying recruitment. Silene ciliata and Armeria caespitosa are two representative species of cryophilic pastures above the tree line in Mediterranean high mountains. No explicit estimations of dispersal kernels have been made so far for these or other high-mountain plants. Such data could help to predict their dispersal and recruitment patterns in a context of changing environments under ongoing global warming.

METHODS

We used an inverse modelling approach to analyse effective seed dispersal patterns in five populations of both Silene ciliata and Armeria caespitosa along an altitudinal gradient in Sierra de Guadarrama (Madrid, Spain). We considered four commonly employed two-dimensional seedling dispersal kernels exponential-power, 2Dt, WALD and log-normal.

KEY RESULTS

No single kernel function provided the best fit across all populations, although estimated mean dispersal distances were short (<1 m) in all cases. S. ciliata did not exhibit significant among-population variation in mean dispersal distance, whereas significant differences in mean dispersal distance were found in A. caespitosa. Both S. ciliata and A. caespitosa exhibited among-population variation in the fecundity parameter and lacked significant variation in kernel shape.

CONCLUSIONS

This study illustrates the complexity of intraspecific variation in the processes underlying recruitment, showing that effective dispersal kernels can remain relatively invariant across populations within particular species, even if there are strong variations in demographic structure and/or physical environment among populations, while the invariant dispersal assumption may not hold for other species in the same environment. Our results call for a case-by-case analysis in a wider range of plant taxa and environments to assess the prevalence and magnitude of intraspecific dispersal variation.

High rates of gene flow by pollen and seed in oak populations across Europe

Gene flow is a key factor in the evolution of species, influencing effective population size, hybridisation and local adaptation. We analysed local gene flow in eight stands of white oak (mostly Quercus petraea and Q. robur, but also Q. pubescens and Q. faginea) distributed across Europe.

Adult trees within a given area in each stand were exhaustively sampled (range [239, 754], mean 423), mapped, and acorns were collected ([17,147], 51) from several mother trees ([3], [47], 23). Seedlings ([65,387], 178) were harvested and geo-referenced in six of the eight stands. Genetic information was obtained from screening distinct molecular markers spread across the genome, genotyping each tree, acorn or seedling. All samples were thus genotyped at 5–8 nuclear microsatellite loci. Fathers/parents were assigned to acorns and seedlings using likelihood methods. Mating success of male and female parents, pollen and seed dispersal curves, and also hybridisation rates were estimated in each stand and compared on a continental scale.

On average, the percentage of the wind-borne pollen from outside the stand was 60%, with large variation among stands (21–88%). Mean seed immigration into the stand was 40%, a high value for oaks that are generally considered to have limited seed dispersal. However, this estimate varied greatly among stands (20–66%). Gene flow was mostly intraspecific, with large variation, as some trees and stands showed particularly high rates of hybridisation.

Our results show that mating success was unevenly distributed among trees. The high levels of gene flow suggest that geographically remote oak stands are unlikely to be genetically isolated, questioning the static definition of gene reserves and seed stands.

Seeing the forest through the trees: comprehensive inference on individual mating patterns in a mixed stand of Quercus robur and Q. petraea

BACKGROUND AND AIMS

Sexual reproduction is one of the most important moments in a life cycle, determining the genetic composition of individual offspring. Controlled pollination experiments often show high variation in the mating system at the individual level, suggesting a persistence of individual variation in natural populations. Individual variation in mating patterns may have significant adaptive implications for a population and for the entire species. Nevertheless, field data rarely address individual differences in mating patterns, focusing rather on averages. This study aimed to quantify individual variation in the different components of mating patterns.

METHODS

Microsatellite data were used from 421 adult trees and 1911 seeds, structured in 72 half-sib families collected in a single mixed stand of Quercus robur and Q. petraea in northern Poland. Using a Bayesian approach, mating patterns were investigated, taking into account pollen dispersal, male fecundity, possible hybridization and heterogeneity in immigrant pollen pools.

KEY RESULTS

Pollen dispersal followed a heavy-tailed distribution (283 m on average). In spite of high pollen mobility, immigrant pollen pools showed strong genetic structuring among mothers. At the individual level, immigrant pollen pools showed highly variable divergence rates, revealing that sources of immigrant pollen can vary greatly among particular trees. Within the stand, the distribution of male fecundity appeared highly skewed, with a small number of dominant males, resulting in a ratio of census to effective density of pollen donors of 5·3. Male fecundity was not correlated with tree diameter but showed strong cline-like spatial variation. This pattern can be attributed to environmental variation. Quercus petraea revealed a greater preference (74 %) towards intraspecific mating than Q. robur (36 %), although mating preferences varied among trees.

CONCLUSIONS

Mating patterns can reveal great variation among individuals, even within a single even-age stand. The results show that trees can mate assortatively, with little respect for spatial proximity. Such selective mating may be a result of variable combining compatibility among trees due to genetic and/or environmental factors.

Evaluating contemporary pollen dispersal in two common grassland species Ranunculus bulbosus L. (Ranunculaceae) and Trifolium montanum L. (Fabaceae) using an experimental approach

Pollen flow is a key biological process that connects plant populations, preventing genetic impoverishment and inbreeding. Pollen-mediated long-distance dispersal (LDD) events are especially important for plant species in increasingly fragmented landscapes. Patterns of pollen dispersal were directly estimated and dispersal kernels modelled in an experimental population of Ranunculus bulbosus and Trifolium montanum to determine the potential for LDD. Eight and 11 microsatellite markers were used for R. bulbosus and T. montanum, respectively, to run a likelihood-based paternity analysis on randomly chosen offspring (Ntotal  = 180 per species) from five maternal plants. High rates of selfing were found in R. bulbosus (average 45.7%), while no selfing was observed in T. montanum. The majority (60%) of mating events occurred at very short distances: the median of the observed dispersal distances was 0.8 m in both species, and the average distances were 15.9 and 10.3 m in R. bulbosus and T. montanum, respectively. Modelling the pollen dispersal kernel with four different distribution functions (exponential-power, geometric, 2Dt and Weibull) indicated that the best fit for both species was given by a Weibull function. Yet, the tail of the T. montanum pollen dispersal kernel was thinner than in R. bulbosus, suggesting that the probability for LDD is higher in the latter species. Even though the majority of pollen dispersal occurred across short distances, the detection of several mating events up to 362 m (R. bulbosus) and 324 m (T. montanum) suggests that pollen flow may be sufficient to ensure population connectivity in these herb species across fragmented grasslands in Swiss agricultural landscapes.

Assessment of spatial discordance of primary and effective seed dispersal of European beech (Fagus sylvatica L.) by ecological and genetic methods

Spatial discordance between primary and effective dispersal in plant populations indicates that postdispersal processes erase the seed rain signal in recruitment patterns. Five different models were used to test the spatial concordance of the primary and effective dispersal patterns in a European beech (Fagus sylvatica) population from central Spain. An ecological method was based on classical inverse modelling (SSS), using the number of seed/seedlings as input data. Genetic models were based on direct kernel fitting of mother-to-offspring distances estimated by a parentage analysis or were spatially explicit models based on the genotype frequencies of offspring (competing sources model and Moran-Clark's Model). A fully integrated mixed model was based on inverse modelling, but used the number of genotypes as input data (gene shadow model). The potential sources of error and limitations of each seed dispersal estimation method are discussed. The mean dispersal distances for seeds and saplings estimated with these five methods were higher than those obtained by previous estimations for European beech forests. All the methods show strong discordance between primary and effective dispersal kernel parameters, and for dispersal directionality. While seed rain was released mostly under the canopy, saplings were established far from mother trees. This discordant pattern may be the result of the action of secondary dispersal by animals or density-dependent effects; that is, the Janzen-Connell effect.

Integrating frugivory and animal movement: a review of the evidence and implications for scaling seed dispersal

General principles about the consequences of seed dispersal by animals for the structure and dynamics of plant populations and communities remain elusive. This is in part because seed deposition patterns emerge from interactions between frugivore behaviour and the distribution of food resources, both of which can vary over space and time. Here we advocate a frugivore-centred, process-based, synthetic approach to seed dispersal research that integrates seed dispersal ecology and animal movement across multiple spatio-temporal scales. To guide this synthesis, we survey existing literature using paradigms from seed dispersal and animal movement. Specifically, studies are discussed with respect to five criteria: selection of focal organisms (animal or plant); measurement of animal movement; characterization of seed shadow; animal, plant and environmental factors included in the study; and scales of the study. Most studies focused on either frugivores or plants and characterized seed shadows directly by combining gut retention time with animal movement data or indirectly by conducting maternity analysis of seeds. Although organismal traits and environmental factors were often measured, they were seldom used to characterize seed shadows. Multi-scale analyses were rare, with seed shadows mostly characterized at fine spatial scales, over single fruiting seasons, and for individual dispersers. Novel animal- and seed-tracking technologies, remote environmental monitoring tools, and advances in analytical methods can enable effective implementation of a hierarchical mechanistic approach to the study of seed dispersal. This kind of mechanistic approach will provide novel insights regarding the complex interplay between the factors that modulate animal behaviour and subsequently influence seed dispersal patterns across spatial and temporal scales.

Causes and consequences of unequal seedling production in forest trees: a case study in red oaks

Inequality in reproductive success has important implications for ecological and evolutionary dynamics, but lifetime reproductive success is challenging to measure in long-lived species such as forest trees. While seed production is often used as a proxy for overall reproductive success, high mortality of seeds and the potential for trade-offs between seed number and quality draw this assumption into question. Parentage analyses of established seedlings can bring us one step closer to understanding the causes and consequences of variation in reproductive success. In this paper we demonstrate a new method for estimating individual seedling production and average percentage germination, using data from two mixed-species populations of red oaks (Quercus rubra, Q. velutina, Q. falcata, and Q. coccinea). We use these estimates to examine the distribution of female reproductive success and to test the relationship between seedling number and individual seed production, age, and growth rate. We show that both seed and seedling production are highly skewed, roughly conforming to zero-inflated lognormal distributions, rather than to the Poisson or negative-binomial distributions often assumed by population genetics analyses. While the number of established offspring is positively associated with mean annual seed production, a lower proportion of seeds from highly fecund individuals become seedlings. Our red oak populations also show evidence of trade-offs between growth rate and reproductive success. The high degree of inequality in seedling production shown here for red oaks, and by previous studies in other species, suggests that many trees may be more vulnerable to genetic drift than previously thought, if immigration in limited by fragmentation or other environmental changes.

Between-site differences in the scale of dispersal and gene flow in red oak

Background

Nut-bearing trees, including oaks (Quercus spp.), are considered to be highly dispersal limited, leading to concerns about their ability to colonize new sites or migrate in response to climate change. However, estimating seed dispersal is challenging in species that are secondarily dispersed by animals, and differences in disperser abundance or behavior could lead to large spatio-temporal variation in dispersal ability. Parentage and dispersal analyses combining genetic and ecological data provide accurate estimates of current dispersal, while spatial genetic structure (SGS) can shed light on past patterns of dispersal and establishment.

Methodology and Principal Findings

In this study, we estimate seed and pollen dispersal and parentage for two mixed-species red oak populations using a hierarchical Bayesian approach. We compare these results to those of a genetic ML parentage model. We also test whether observed patterns of SGS in three size cohorts are consistent with known site history and current dispersal patterns. We find that, while pollen dispersal is extensive at both sites, the scale of seed dispersal differs substantially. Parentage results differ between models due to additional data included in Bayesian model and differing genotyping error assumptions, but both indicate between-site dispersal differences. Patterns of SGS in large adults, small adults, and seedlings are consistent with known site history (farmed vs. selectively harvested), and with long-term differences in seed dispersal. This difference is consistent with predator/disperser satiation due to higher acorn production at the low-dispersal site. While this site-to-site variation results in substantial differences in asymptotic spread rates, dispersal for both sites is substantially lower than required to track latitudinal temperature shifts.

Conclusions

Animal-dispersed trees can exhibit considerable spatial variation in seed dispersal, although patterns may be surprisingly constant over time. However, even under favorable conditions, migration in heavy-seeded species is likely to lag contemporary climate change.

Genetic evidence for hybridization in red oaks (Quercus sect. Lobatae, Fagaceae)

PREMISE OF THE STUDY

Hybridization is pervasive in many plant taxa, with consequences for species taxonomy, local adaptation, and management. Oaks (Quercus spp.) are thought to hybridize readily yet retain distinct traits, drawing into question the biological species concept for such taxa, but the true extent of gene flow is controversial. Genetic data are beginning to shed new light on this issue, but red oaks (section Lobatae), an important component of North American forests, have largely been neglected. Moreover, gene flow estimates may be sensitive to the choice of life stage, marker type, or genetic structure statistic.

METHODS

We coupled genetic structure data with parentage analyses for two mixed-species stands in North Carolina. Genetic structure analyses of adults (including F(ST), R(ST), G'(ST), and structure) reflect long-term patterns of gene flow, while the percentage of seedlings with parents of two different species reflect current levels of gene flow.

KEY RESULTS

Genetic structure analyses revealed low differentiation in microsatellite allele frequencies between co-occurring species, suggesting past gene flow. However, methods differed in their sensitivity to differentiation, indicating a need for caution when drawing conclusions from a single method. Parentage analyses identified >20% of seedlings as potential hybrids. The species examined exhibit distinct morphologies, suggesting selection against intermediate phenotypes.

CONCLUSIONS

Our results suggest that hybridization between co-occurring red oaks occurs, but that selection may limit introgression, especially at functional loci. However, by providing a source of genetic variation, hybridization could influence the response of oaks and other hybridizing taxa to environmental change.

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

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

Pollen and seed dispersal inferred from seedling genotypes: the Bayesian revolution has passed here too

Understanding precisely how plants disperse their seeds and pollen in their neighbourhood is a central question for both ecologists and evolutionary biologists because seed and pollen dispersal governs both the rate of spread of an expanding population and gene flow within and among populations. The concept of a 'dispersal kernel' has become extremely popular in dispersal ecology as a tool that summarizes how dispersal distributes individuals and genes in space and at a given scale. In this issue of Molecular Ecology, the study by Moran & Clark (2011) (M&C in the following) shows how genotypic and spatial data of established seedlings can be analysed in a Bayesian framework to estimate jointly the pollen and seed dispersal kernels and finally derive a parentage analysis from a full-probability approach. This approach applied to red oak shows important dispersal of seeds (138 m on average) and pollen (178 m on average). For seeds, this estimate contrasts with previous results from inverse modelling on seed trap data (9.3 m). This research gathers several methodological advances made in recent years in two research communities and could become a cornerstone for dispersal ecology.