Two-generation analysis of pollen flow across a landscape. III. Impact of adult population structure

Life history and past demography maintain genetic structure, outcrossing rate, contemporary pollen gene flow of an understory herb in a highly fragmented rainforest

INTRODUCTION

Theory predicts that habitat fragmentation, by reducing population size and increasing isolation among remnant populations, can alter their genetic diversity and structure. A cascade of effects is expected: genetic drift and inbreeding after a population bottleneck, changes in biotic interactions that may affect, as in the case of plants, pollen dynamics, mating system, reproductive success. The detection of the effects of contemporary habitat fragmentation on the genetic structure of populations are conditioned by the magnitude of change, given the few number of generations since the onset of fragmentation, especially for long-lived organisms. However, the present-day genetic structure of populations may bear the signature of past demography events. Here, we examine the effects of rainforest fragmentation on the genetic diversity, population structure, mating system (outcrossing rate), indirect gene flow and contemporary pollen dynamics in the understory herb Aphelandra aurantiaca. Also, we assessed its present-day genetic structure under different past demographic scenarios.

METHODS

Twelve populations of A. aurantiaca were sampled in large (4), medium (3), and small (5) forest fragments in the lowland tropical rainforest at Los Tuxtlas region. Variation at 11 microsatellite loci was assessed in 28-30 reproductive plants per population. In two medium- and two large-size fragments we estimated the density of reproductive plants, and the mating system by analyzing the progeny of different mother plants per population.

RESULTS

Despite prevailing habitat fragmentation, populations of A. aurantiaca possess high genetic variation (He  = 0.61), weak genetic structure (Rst  = 0.037), and slight inbreeding in small fragments. Effective population sizes (Ne ) were large, but slightly lower in small fragments. Migrants derive mostly from large and medium size fragments. Gene dispersal is highly restricted but long distance gene dispersal events were detected. Aphelandra aurantiaca shows a mixed mating system (tm  = 0.81) and the outcrossing rate have not been affected by habitat fragmentation. A strong pollen pool structure was detected due to few effective pollen donors (Nep ) and low distance pollen movement, pointing that most plants received pollen from close neighbors. Past demographic fluctuations may have affected the present population genetic structure as Bayesian coalescent analysis revealed the signature of past population expansion, possibly during warmer conditions after the last glacial maximum.

DISCUSSION

Habitat fragmentation has not increased genetic differentiation or reduced genetic diversity of A. aurantiaca despite dozens of generations since the onset of fragmentation in the region of Los Tuxtlas. Instead, past population expansion is compatible with the lack of observed genetic structure. The predicted negative effects of rainforest fragmentation on genetic diversity and population structure of A. aurantiaca seem to have been buffered owing to its large effective populations and long-distance dispersal events. In particular, its mixed-mating system, mostly of outcrossing, suggests high efficiency of pollinators promoting connectivity and reducing inbreeding. However, some results point that the effects of fragmentation are underway, as two small fragments showed higher membership probabilities to their population of origin, suggesting genetic isolation. Our findings underscore the importance of fragment size to maintain genetic connectivity across the landscape.

Limited Pollen Dispersal Contributes to Population Genetic Structure but Not Local Adaptation in Quercus oleoides Forests of Costa Rica

Background

Quercus oleoides Cham. and Schlect., tropical live oak, is a species of conservation importance in its southern range limit of northwestern Costa Rica. It occurs in high-density stands across a fragmented landscape spanning a contrasting elevation and precipitation gradient. We examined genetic diversity and spatial genetic structure in this geographically isolated and genetically distinct population. We characterized population genetic diversity at 11 nuclear microsatellite loci in 260 individuals from 13 sites. We monitored flowering time at 10 sites, and characterized the local environment in order to compare observed spatial genetic structure to hypotheses of isolation-by-distance and isolation-by-environment. Finally, we quantified pollen dispersal distances and tested for local adaptation through a reciprocal transplant experiment in order to experimentally address these hypotheses.

Results

High genetic diversity is maintained in the population and the genetic variation is significantly structured among sampled sites. We identified 5 distinct genetic clusters and average pollen dispersal predominately occurred over short distances. Differences among sites in flowering phenology and environmental factors, however, were not strictly associated with genetic differentiation. Growth and survival of upland and lowland progeny in their native and foreign environments was expected to exhibit evidence of local adaptation due to the more extreme dry season in the lowlands. Seedlings planted in the lowland garden experienced much higher mortality than seedlings in the upland garden, but we did not identify evidence for local adaptation.

Conclusion

Overall, this study indicates that the Costa Rican Q. oleoides population has a rich population genetic history. Despite environmental heterogeneity and habitat fragmentation, isolation-by-distance and isolation-by-environment alone do not explain spatial genetic structure. These results add to studies of genetic structure by examining a common, tropical tree over multiple habitats and provide information for managers of a successional forest in a protected area.

Long-term effects of habitat fragmentation on mating patterns and gene flow of a tropical dry forest tree, Ceiba aesculifolia (Malvaceae: Bombacoideae)

PREMISE OF THE STUDY

Tropical forest loss and fragmentation isolate and reduce the size of remnant populations with negative consequences for mating patterns and genetic structure of plant species. In a 4-yr study, we determined the effect of fragmentation on mating patterns and pollen pool genetic structure of the tropical tree Ceiba aesculifolia in two habitat conditions: isolated trees in disturbed areas (≤3 trees/ha), and trees (≥6 trees/ha) in undisturbed mature forest. •

METHODS

Using six allozyme loci, we estimated the outcrossing rate (tm), the mean relatedness of progeny (rp) within and between fruits, the degree of genetic structure of pollen pools (Φft), and the effective number of pollen donors (Nep). •

KEY RESULTS

The outcrossing rates reflected a strict self-incompatible species. Relatedness of progeny within fruits was similar for all populations, revealing single sires within fruits. However, relatedness of progeny between fruits within trees was consistently greater for trees in fragmented conditions across 4 yr. We found high levels of genetic structure of pollen pools in all populations with more structure in isolated trees. The effective number of pollen donors was greater for trees in undisturbed forest than in disturbed conditions. •

CONCLUSIONS

Our study showed that the progeny produced by isolated trees in disturbed habitats are sired by a fraction of the diversity of pollen donors found in conserved forests. The foraging behavior of bats limits the exchange of pollen between trees, causing higher levels of progeny relatedness in isolated trees.

Extensive contemporary pollen-mediated gene flow in two herb species, Ranunculus bulbosus and Trifolium montanum, along an altitudinal gradient in a meadow landscape

BACKGROUND AND AIMS

Genetic connectivity between plant populations allows for exchange and dispersal of adaptive genes, which can facilitate plant population persistence particularly in rapidly changing environments.

METHODS

Patterns of historic gene flow, flowering phenology and contemporary pollen flow were investigated in two common herbs, Ranunculus bulbosus and Trifolium montanum, along an altitudinal gradient of 1200-1800 m a.s.l. over a distance of 1 km among five alpine meadows in Switzerland.

KEY RESULTS

Historic gene flow was extensive, as revealed by Fst values of 0·01 and 0·007 in R. bulbosus and T. montanum, respectively, by similar levels of allelic richness among meadows and by the grouping of all individuals into one genetic cluster. Our data suggest contemporary pollen flow is not limited across altitudes in either species but is more pronounced in T. montanum, as indicated by the differential decay of among-sibships correlated paternity with increasing spatial distance. Flowering phenology among meadows was not a barrier to pollen flow in T. montanum, as the large overlap between meadow pairs was consistent with the extensive pollen flow. The smaller flowering overlap among R. bulbosus meadows might explain the slightly more limited pollen flow detected.

CONCLUSIONS

High levels of pollen flow among altitudes in both R. bulbosus and T. montanum should facilitate exchange of genes which may enhance adaptive responses to rapid climate change.

Comparison of Bayesian clustering and edge detection methods for inferring boundaries in landscape genetics

Recently, techniques available for identifying clusters of individuals or boundaries between clusters using genetic data from natural populations have expanded rapidly. Consequently, there is a need to evaluate these different techniques. We used spatially-explicit simulation models to compare three spatial Bayesian clustering programs and two edge detection methods. Spatially-structured populations were simulated where a continuous population was subdivided by barriers. We evaluated the ability of each method to correctly identify boundary locations while varying: (i) time after divergence, (ii) strength of isolation by distance, (iii) level of genetic diversity, and (iv) amount of gene flow across barriers. To further evaluate the methods' effectiveness to detect genetic clusters in natural populations, we used previously published data on North American pumas and a European shrub. Our results show that with simulated and empirical data, the Bayesian spatial clustering algorithms outperformed direct edge detection methods. All methods incorrectly detected boundaries in the presence of strong patterns of isolation by distance. Based on this finding, we support the application of Bayesian spatial clustering algorithms for boundary detection in empirical datasets, with necessary tests for the influence of isolation by distance.

Influence of cryptic population structure on observed mating patterns in the wild progenitor of maize (Zea mays ssp. parviglumis)

Indirect two-generation analysis of pollen flow has proven to be an effective alternative to exhaustive paternity analysis in numerous plant populations. In this investigation, the method is extended to an annual wild maize species, Zea mays ssp. parviglumis (Poaceae). Our analysis of mating system in parviglumis revealed high levels of outcrossing and higher biparental inbreeding than typically observed in grasses. Pollen dispersal analysis suggested low levels of long-distance dispersal. Given previous evidence for intrapopulation genetic structure in parviglumis populations, we explored the impact of cryptic population structure on estimates of mating system and pollen flow. Subpopulations inferred through spatially explicit Bayesian assignment showed markedly different values for both mating system parameters and pollen flow than the entire population. Finally, a novel method of pollen haplotype assignment revealed nonrandom mating consistent with intrapopulation structure. These results indicate parviglumis could be particularly susceptible to habitat fragmentation currently occurring throughout Mexico due to recent changes in land use.

Restricted pollen flow of Dieffenbachia seguine populations in fragmented and continuous tropical forest

Habitat fragmentation can change the ecological context of populations, rupturing genetic connectivity among them, changing genetic structure, and increasing the loss of genetic diversity. We analyzed mating system and pollen structure in two population fragments and two continuous forest populations of Dieffenbachia seguine (Araceae), an insect-pollinated understory herb in the tropical rain forest of Los Tuxtlas, México, using nine allozyme loci. Mating system analysis indicated almost complete outcrossing but some inbreeding among the adults. Pollen structure analysis indicated highly restricted pollen flow, both within and among populations. We showed that the effective pollination neighborhood was small in all populations, and slightly (though not significantly) smaller in fragments, partially as a consequence of an increase in density of reproductive individuals in those fragments. Using assignment analysis, we showed that all populations were strongly structured, suggesting that pollen and seed flow across the Los Tuxtlas landscape has been spatially restricted, though sufficient to maintain connectedness. Forest fragmentation at Los Tuxtlas has (so far) had limited impact on pollen dynamics, despite the changing ecological context, with reduced pollinator abundance being partially offset by increased flowering density in fragments. Continued outcrossing and limited pollen immigration, coupled with more extensive seed migration, should maintain genetic connectedness in D. seguine, if fragmentation is not further exacerbated by additional deforestation.

Powers of discerning: challenges to understanding dispersal processes in natural populations

In this issue of Molecular Ecology, authors Robledo-Arnuncio & Garcia present a compelling approach for quantifying seed dispersal in plant populations. Building upon methods previously used for quantification of pollen dispersal, the authors not only examine the behaviour of the model with respect to sample sizes, dispersal distance, and the kurtosis of the dispersal function but also provide an empirical example using Prunus mahaleb.

Extensive pollen dispersal in a bird-pollinated shrub,Calothamnus quadrifidus, in a fragmented landscape

Pollen dispersal was investigated in six populations of Calothamnus quadrifidus, a bird-pollinated shrub in the fragmented agricultural region of southern Western Australia. Paternity analysis using six microsatellite loci identified a pollen source within populations for 67% of seedlings, and the remainder were assumed to have arisen from pollen sources outside the populations. Outcrossing was variable, ranging from 5% to 82%, and long-distance pollen dispersal was observed in all populations with up to 43% of pollen sourced from outside the populations over distances of up to 5 km. This extensive pollen immigration was positively associated with population size but not isolation. Comparison of two populations of similar size but different density showed greater internal pollination and less selfing in the denser population, suggesting an influence of density on pollinator behaviour. The study revealed extensive long-distance pollen dispersal for C. quadrifidus within this fragmented agricultural landscape and highlighted the interaction between reserve populations and isolated road verge remnants in maintaining genetic connectivity at the landscape scale.

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‐mediated gene dispersal within continuous and fragmented populations of a forest understorey species,Trillium cuneatum

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

Genetic Structure Inside a Declining Red Oak Community in Old-Growth Forest

Problems with oak regeneration have been documented in the last 50 years at numerous sites in the Midwestern United States. We applied nuclear microsatellites to examine the demographic and fine-scale spatial genetic structure of red oaks in two old-growth stands in Indiana. Oaks in one stand have declined in numbers over the past several decades whereas oaks in the other, smaller stand have increased. Large amounts of genetic variation were maintained within stands, and there was slight but significant differentiation among stands. There was significant but weak isolation by distance genetic structure within the large stand, likely reflecting family structure. No significant differences exist in allele frequencies or in levels of genetic diversity between cohorts that remain well represented within each stand, even between medium-sized adults and those antedating European settlement of the area. However, a virtual absence of smaller size classes in the forest interior of the large stand represents the early stages of a genetic bottleneck in what had been the core habitat of this stand. Whether future generations of this old-growth stand will retain the present genetic character depends on the oaks regenerating at the forest margins, absent any major changes in disturbance regimes. Similar demographic and genetic dynamics are likely occurring in a large number of remnant oak forests across the Midwest.

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.

Measuring the genetic structure of the pollen pool as the probability of paternal identity

Contemporary pollen flow in forest plant species is measured by the probability of paternal identity (PPI) for two randomly sampled offspring, drawn from a single female, and contrasting that with PPI for two random offspring, drawn from different females. Two different estimation strategies have emerged: (a) an indirect approach, using the 'genetic structure' of the pollen received by different mothers and (b) a direct approach, based on parentage analysis. The indirect strategy is somewhat limited by the assumptions, but is widely useful. The direct approach is most appropriate where a large majority of the true fathers can be identified exactly, which is sometimes possible with high-resolution SSR markers. Using the parentage approach, we develop estimates of PPI, showing that the obvious estimates are severely biased, and providing an unbiased alternative. We then illustrate the methods with SSR data from a 36-tree isolated population of Pinus sylvestris from the Meseta region of Spain, for which categorical paternity assignment was available for over 95% of offspring. For all the females combined, we estimate that PPI=0.0425, indicating uneven male reproductive contributions. Different (but overlapping) arrays of males pollinate different females, and for the average female, PPI=0.317, indicating substantial 'pollen structure' for the population. We also relate the direct measures of PPI to those available from indirect approaches, and show that they are generally comparable.

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.

Pollen movement in declining populations of California Valley oak, Quercus lobata: where have all the fathers gone?

The fragmented populations and reduced population densities that result from human disturbance are issues of growing importance in evolutionary and conservation biology. A key issue is whether remnant individuals become reproductively isolated. California Valley oak (Quercus lobata) is a widely distributed, endemic species in California, increasingly jeopardized by anthropogenic changes in biota and land use. We studied pollen movement in a savannah population of Valley oak at Sedgwick Reserve, Santa Barbara County, to estimate effective number of pollen donors (Nep) and average distance of effective pollen movement (delta). Using twogener, our recently developed hybrid model of paternity and genetic structure treatments that analyses maternal and progeny multilocus genotypes, we found that current Nep = 3.68 individuals. Based on an average adult density of d= 1.19 stems/ha, we assumed a bivariate normal distribution to model current average pollen dispersal distance (delta) and estimated delta= 64.8 m. We then deployed our parameter estimates in spatially explicit models of the Sedgwick population to evaluate the extent to which Nep may have changed, as a consequence of progressive stand thinning between 1944 and 1999. Assuming that pollen dispersal distance has not changed, we estimate Nep was 4.57 individuals in 1944, when stand density was 1.48. Both estimates indicate fewer effective fathers than one might expect for wind-pollinated species and fewer than observed elsewhere. The results presented here provide a basis for further refinements on modelling pollen movement. If the trends continue, then ongoing demographic attrition could further reduce neighbourhood size in Valley oak resulting in increased risk of reproductive failure and genetic isolation.

Two-generation analysis of pollen flow across a landscape. IV. Estimating the dispersal parameter

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.