Unusually limited pollen dispersal and connectivity of Pedunculate oak (Quercus robur) refugial populations at the species' southern range margin

Highly clustered mating networks in naturally fragmented riparian tree populations

Understanding how spatial patterns of mating and gene flow respond to habitat loss and geographical isolation is a crucial aspect of forest fragmentation genetics. Naturally fragmented riparian tree populations exhibit unique characteristics that significantly influence these patterns. In this study, we investigate mating patterns, pollen-mediated gene flow, and genetic diversity in relict populations of Frangula alnus in southern Spain by testing specific hypotheses related to the riparian habitat. We employ a novel approach that combines paternity analysis, particularly suited for small and isolated populations, with complex network theory and Bayesian models to predict mating likelihood among tree pairs. Our findings reveal a prevalence of short-distance pollination, resulting in spatially driven local mating clusters with a distinct subset of trees being highly connected in the mating network. Additionally, we observe numerous pollination events over distances of hundreds of metres and considerable pollen immigration. Local neighbourhood density is the primary factor influencing within-population mating patterns and pollen dispersal; moreover, mating network properties reflect the population's size and spatial configuration. Conversely, among-population pollen dispersal is mainly determined by tree size, which influences floral display. Our results do not support a major role of directional pollen dispersal in longitudinal trends of genetic diversity. We provide evidence that long-term fragmented tree populations persist in unique environments that shape mating patterns and impose constraints to pollen-mediated gene flow. Nevertheless, even seemingly strongly isolated populations can maintain functional connectivity over extended periods, especially when animal-mediated mating networks promote genetic diversity, as in this riparian tree species.

Hawkmoth and bee pollinators impact pollen dispersal at the landscape but not local scales in two species of Oenothera

PREMISE

Animal pollinators play an important role in pollen dispersal. Here, we assessed differences in pollen and seed dispersal and the role of pollinator functional groups with different foraging behaviors in generating patterns of genetic diversity over similar geographic ranges for two closely related taxa. We focused on two members of Oenothera section Calylophus (Onagraceae) that co-occur on gypsum outcrops throughout the northern Chihuahuan Desert but differ in floral phenotype and primary pollinator: Oenothera gayleana (bee) and O. hartwegii subsp. filifolia (hawkmoth).

METHODS

We measured breeding system and floral traits and studied gene flow and population differentiation at the local (<13 km; four populations) and landscape (60-440 km; five populations) scales using 10-11 nuclear (pollen dispersal) and three plastid (seed dispersal) microsatellite markers.

RESULTS

Both taxa were self-incompatible and floral traits were consistent with expectations for different pollinators. Seed and pollen dispersal patterns were distinctly different for both species. We found no evidence of genetic structure at the local scale but did at the landscape scale; O. gayleana showed greater differentiation and significant isolation by distance than in O. hartwegii subsp. filifolia. The plastid data were consistent with gravity dispersal of seeds and suggest that pollen dispersal is the principal driver of genetic structure in both species.

CONCLUSIONS

We demonstrated that pollinator functional groups can impact genetic differentiation in different and predictable ways. Hawkmoths, with larger foraging distances, can maintain gene flow across greater spatial scales than bees.

Joint analysis of microsatellites and flanking sequences enlightens complex demographic history of interspecific gene flow and vicariance in rear-edge oak populations

Inference of recent population divergence requires fast evolving markers and necessitates to differentiate shared genetic variation caused by ancestral polymorphism and gene flow. Theoretical research shows that the use of compound marker systems integrating linked polymorphisms with different mutational dynamics, such as a microsatellite and its flanking sequences, can improve estimation of population structure and inference of demographic history, especially in the case of complex population dynamics. However, empirical application in natural populations has so far been limited by lack of suitable methods for data collection. A solution comes from the development of sequence-based microsatellite genotyping which we used to study molecular variation at 36 sequenced nuclear microsatellites in seven Quercus canariensis and four Q. faginea rear-edge populations across Algeria. We aim to decipher their taxonomic relationship, past evolutionary history and recent demographic trajectory. First, we compare the estimation of population genetics parameters and simulation-based inference of demographic history from microsatellite sequence alone, flanking sequence alone or the combination of linked microsatellite and flanking sequence variation. Second, we apply random forest approximate Bayesian computation to identify which of these sequence types is most informative. Whereas analysing microsatellite variation alone indicates recent interspecific gene flow, additional information gained by integrating nucleotide variation in flanking sequences, by reducing homoplasy, suggests ancient interspecific gene flow followed by drift in isolation instead. The weight of each polymorphism in the inference also demonstrates the value of linked variations with contrasted mutation dynamic to improve estimation of both demographic and mutational parameters.

Rules of Plant Species Ranges: Applications for Conservation Strategies

Earth is changing rapidly and so are many plant species’ ranges. Here, we synthesize eco-evolutionary patterns found in plant range studies and how knowledge of species ranges can inform our understanding of species conservation in the face of global change. We discuss whether general biogeographic “rules” are reliable and how they can be used to develop adaptive conservation strategies of native plant species across their ranges. Rules considered include (1) factors that set species range limits and promote range shifts; (2) the impact of biotic interactions on species range limits; (3) patterns of abundance and adaptive properties across species ranges; (4) patterns of gene flow and their implications for genetic rescue, and (5) the relationship between range size and conservation risk. We conclude by summarizing and evaluating potential species range rules to inform future conservation and management decisions. We also outline areas of research to better understand the adaptive capacity of plants under environmental change and the properties that govern species ranges. We advise conservationists to extend their work to specifically consider peripheral and novel populations, with a particular emphasis on small ranges. Finally, we call for a global effort to identify, synthesize, and analyze prevailing patterns or rules in ecology to help speed conservation efforts.

Impact of Gene Flow and Introgression on the Range Wide Genetic Structure of Quercus robur (L.) in Europe

As for most other temperate broadleaved tree species, large-scale genetic inventories of pedunculate oak (Quercus robur L.) have focused on the plastidial genome, which showed the impact of post-glacial recolonization and manmade seed transfer. However, how have pollen mediated gene flow and introgression impacted the large-scale genetic structure? To answer these questions, we did a genetic inventory on 1970 pedunculate oak trees from 197 locations in 13 European countries. All samples were screened with a targeted sequencing approach on a set of 381 polymorphic loci (356 nuclear SNPs, 3 nuclear InDels, 17 chloroplast SNPs, and 5 mitochondrial SNPs). In a former analysis with additional 1763 putative Quercus petraea trees screened for the same gene markers we obtained estimates on the species admixture of all pedunculate oak trees. We identified 13 plastidial haplotypes, which showed a strong spatial pattern with a highly significant autocorrelation up to a range of 1250 km. Significant spatial genetic structure up to 1250 km was also observed at the nuclear loci. However, the differentiation at the nuclear gene markers was much lower compared to the organelle gene markers. The matrix of genetic distances among locations was partially correlated between nuclear and organelle genomes. Bayesian clustering analysis revealed the best fit to the data for a sub-division into two gene pools. One gene pool is dominating the west and the other is the most abundant in the east. The western gene pool was significantly influenced by introgression from Quercus petraea in the past. In Germany, we identified a contact zone of pedunculate oaks with different introgression intensity, likely resulting from different historical levels of introgression in glacial refugia or during postglacial recolonization. The main directions of postglacial recolonization were south to north and south to northwest in West and Central Europe, and for the eastern haplotypes also east to west in Central Europe. By contrast, the pollen mediated gene flow and introgression from Q. petraea modified the large-scale structure at the nuclear gene markers with significant west–east direction.

Small-scale genetic structure and mating patterns in an extensive sessile oak forest (Quercus petraea (Matt.) Liebl.)

Oaks (Quercus) are major components of temperate forest ecosystems in the Northern Hemisphere where they form intermediate or climax communities. Sessile oak (Quercus petraea) forests represent the climax vegetation in eastern Germany and western Poland. Here, sessile oak forms pure stands or occurs intermixed with Scots Pine (Pinus sylvestris). A large body of research is available on gene flow, reproduction dynamics, and genetic structure in fragmented landscapes and mixed populations. At the same time, our knowledge regarding large, contiguous, and monospecific populations is considerably less well developed. Our study is an attempt to further develop our understanding of the reproduction ecology of sessile oak as an ecologically and economically important forest tree by analyzing mating patterns and genetic structure within adult trees and seedlings originating from one or two reproduction events in an extensive, naturally regenerating sessile oak forest. We detected positive spatial genetic structure up to 30 meters between adult trees and up to 40 meters between seedlings. Seed dispersal distances averaged 8.4 meters. Pollen dispersal distances averaged 22.6 meters. In both cases, the largest proportion of the dispersal occurred over short distances. Dispersal over longer distances was more common for pollen but also appeared regularly for seeds. The reproductive success of individual trees was highly skewed. Only 41 percent of all adult trees produced any offspring while the majority did not participate in reproduction. Among those trees that contributed to the analyzed seedling sample, 80 percent contributed 1-3 gametes. Only 20 percent of all parent trees contributed four or more gametes. However, these relatively few most fertile trees contributed 51 percent of all gametes within the seedling sample. Vitality and growth differed significantly between reproducing and nonreproducing adult trees with reproducing trees being more vital and vigorous than nonreproducing individuals. Our study demonstrates that extensive, apparently homogenous oak forests are far from uniform on the genetic level. On the contrary, they form highly complex mosaics of remarkably small local neighborhoods. This counterbalances the levelling effect of long-distance dispersal and may increase the species' adaptive potential. Incorporating these dynamics in the management, conservation, and restoration of oak forests can support the conservation of forest genetic diversity and assist those forests in coping with environmental change.

Extensive pollen-mediated gene flow across intensively managed landscapes in an insect-pollinated shrub native to semiarid habitats

Our knowledge of the impact of landscape fragmentation on gene flow patterns is mainly drawn from tropical and temperate ecosystems, where landscape features, such as the distance of a tree to the forest edge, drive connectivity and mating patterns. Yet, the structure of arid and semiarid plant communities - with open canopies and a scattered distribution of trees - differs greatly from those that are well-characterized in the literature. As a result, we ignore whether the documented consequences of landscape fragmentation on plant mating and gene flow patterns also hold for native plant communities in arid and semiarid regions. We investigated the relative contribution of plant traits, pollinator activity, and individual neighbourhood in explaining variation in mating and gene flow patterns of an insect-pollinated semiarid arborescent shrub, Ziziphus lotus, at three sites embedded in highly altered agriculture landscapes. We used 14 SSRs, seed paternity analyses, and individual mixed effect mating models (MEMMi) to estimate the individual mating variables and the pollen dispersal kernel at each site. Individual spatial location, flower density, and floral visitation rate explained most of the variation of mating variables. Unexpectedly, individual correlated paternity was very low and shrubs surrounded by the most degraded matrix exhibited an increased fraction of pollen immigration and a high effective number of pollen donors per mother shrub. Overall, our results reveal that an active pollinator assemblage ensures highly efficient mating, and maintains pollen-mediated gene flow and notable connectivity levels, even in highly altered landscapes, potentially halting genetic isolation within and between distant sites.

Identification of determinants of pollen donor fecundity using the hierarchical neighborhood model

Individual differences in male reproductive success drive genetic drift and natural selection, altering genetic variation and phenotypic trait distributions in future generations. Therefore, identifying the determinants of reproductive success is important for understanding the ecology and evolution of plants. Here, based on the spatially explicit mating model (the neighborhood model), we develop a hierarchical probability model that links co-dominant genotypes of offspring and candidate parents with phenotypic determinants of male reproductive success. The model accounts for pollen dispersal, genotyping errors as well as individual variation in selfing, pollen immigration, and differentiation of immigrant pollen pools. Unlike the classic neighborhood model approach, our approach is specially designed to account for excessive variation (overdispersion) in male fecundity. We implemented a Bayesian estimation method (the Windows computer program available at: https://www.ukw.edu.pl/pracownicy/plik/igor_chybicki/1806/) that, among others, allows for selecting phenotypic variables important for male fecundity and assessing the fraction of variance in fecundity (R² ) explained by selected variables. Simulations showed that our method outperforms both the classic neighborhood model and the two-step approach, where fecundities and the effects of phenotypic variables are estimated separately. The analysis of two data examples showed that in wind-pollinated trees, male fecundity depends on both the amount of produced pollen and the ability to pollen spread. However, despite that the tree size was positively correlated with male fecundity, it explained only a fraction of the total variance in fecundity, indicating the presence of additional factors. Finally, case studies highlighted the importance of accounting for pollen dispersal in the estimation of fecundity determinants.

Candidate gene SNP variation in floodplain populations of pedunculate oak (Quercus robur L.) near the species' southern range margin: Weak differentiation yet distinct associations with water availability

Populations residing near species' low-latitude range margins (LLMs) often occur in warmer and drier environments than those in the core range. Thus, their genetic composition could be shaped by climatic drivers that differ from those occurring at higher latitudes, resulting in potentially adaptive variants of conservation value. Such variants could facilitate the adaptation of populations from other portions of the geographical range to similar future conditions anticipated under ongoing climate change. However, very few studies have assessed standing genetic variation at potentially adaptive loci in natural LLM populations. We investigated standing genetic variation at single nucleotide polymorphisms (SNPs) located within 117 candidate genes and its links to putative climatic selection pressures across 19 pedunculate oak (Quercus robur L.) populations distributed along a regional climatic gradient near the species' southern range margin in southeastern Europe. These populations are restricted to floodplain forests along large lowland rivers, whose hydric regime is undergoing significant shifts under modern rapid climate change. The populations showed very weak geographical structure, suggesting extensive genetic connectivity and gene flow or shared ancestry. We identified eight (6.2%) positive F_ST -outlier loci, and genotype-environment association analyses revealed consistent associations between SNP allele frequencies and several climatic variables linked to water availability. A total of 61 associations involving 37 SNPs (28.5%) from 35 annotated genes provided important insights into putative functional mechanisms in our system. Our findings provide empirical support for the role of LLM populations as sources of potentially adaptive variation that could enhance species' resilience to climate change-related pressures.

Great Genetic Diversity but High Selfing Rates and Short-Distance Gene Flow Characterize Populations of a Tree (Foetidia; Lecythidaceae) in the Fragmented Tropical Dry Forest of the Mascarene Islands

Following the global trend of deforestation and degradation, tropical dry forests in the Mascarenes archipelago on Reunion has undergone harsh reduction and fragmentation within 3 centuries of human occupation. We investigated the genetic diversity, mating system, and gene flow in fragmented populations of the native tree Foetidia mauritiana (Lecythidaceae) on Reunion, using microsatellite genotyping of adults (in- and ex situ) and seed progenies (in situ only). To test genetic isolation between the Mascarene islands, we also genotyped conspecific adults on Mauritius, and trees of Foetidia rodriguesiana on Rodrigues. We found a high genetic diversity among the trees on Reunion, but no population structure (G'ST: 0.039-0.090), and an increase of the fixation index (FIS) from adults to progenies. A subsequent analysis of mating systems from progeny arrays revealed selfing rates >50% in fragmented populations and close to 100% in lone trees. A paternity analysis revealed pollen flow ranging from 15.6 to 296.1 m within fragments. At broader scale, the populations of F. mauritiana on Reunion and Mauritius are genetically differentiated. The morphologically allied taxa F. rodriguesiana and F. mauritiana are clearly isolated. Therefore, this case study shows that genetic diversity may persist after deforestation, especially in long-lived tree species, but the reproductive features may be deeply altered during this process. This would explain the low seed production and the absence of recruitment in F. mauritiana. Restoration programs should take into account these features, as well as the importance that trees ex situ represent in restoring and conserving diversity.

Combining rare alleles and grouped pollen donors to assign paternity in pollen dispersal studies

PREMISE

Pollen dispersal plays a critical role in gene flow of seed plants. Most often, pollen dispersal is measured using paternity assignment. However, this approach can be time-consuming because it typically entails genotyping all pollen donors, receptors, and offspring at several molecular markers.

METHODS

We developed a faster, simpler protocol to track paternity, using pollen receptors and grouped pollen donors that possess rare alleles. We tested this approach using wind-pollinated Amaranthus tuberculatus and insect-pollinated Solanum lycopersicum. After screening potential markers for rare alleles, we grew both species in experimental arrays under field conditions.

RESULTS

All tested A. tuberculatus seeds and 97% of S. lycopersicum fruits could be assigned to the grouped pollen donors using each of two markers. From these results, we could infer paternity of untested offspring and assess pollen dispersal patterns in each array.

DISCUSSION

By combining rare alleles and grouped pollen donors, we could assess pollen dispersal for both species and across all arrays after genotyping a small number of pollen donors and a representative subset of offspring. While directly applicable to A. tuberculatus and S. lycopersicum, this approach could be used in other species to assess pollen dispersal under field conditions.

An analysis of mating biases in trees

Assortative mating is a deviation from random mating based on phenotypic similarity. As it is much better studied in animals than in plants, we investigate for trees whether kinship of realized mating pairs deviates from what is expected from the set of potential mates and use this information to infer mating biases that may result from kin recognition and/or assortative mating. Our analysis covers 20 species of trees for which microsatellite data is available for adult populations (potential mates) as well as seed arrays. We test whether mean relatedness of observed mating pairs deviates from null expectations that only take pollen dispersal distances into account (estimated from the same data set). This allows the identification of elevated as well as reduced kinship among realized mating pairs, indicative of positive and negative assortative mating, respectively. The test is also able to distinguish elevated biparental inbreeding that occurs solely as a result of related pairs growing closer to each other from further assortativeness. Assortative mating in trees appears potentially common but not ubiquitous: nine data sets show mating bias with elevated inbreeding, nine do not deviate significantly from the null expectation, and two show mating bias with reduced inbreeding. While our data sets lack direct information on phenology, our investigation of the phenological literature for each species identifies flowering phenology as a potential driver of positive assortative mating (leading to elevated inbreeding) in trees. Since active kin recognition provides an alternative hypothesis for these patterns, we encourage further investigations on the processes and traits that influence mating patterns in trees.

Refining predictions of population decline at species' rear edges

According to broad-scale application of biogeographical theory, widespread retractions of species' rear edges should be seen in response to ongoing climate change. This prediction rests on the assumption that rear edge populations are "marginal" since they occur at the limit of the species' ecological tolerance and are expected to decline in performance as climate warming pushes them to extirpation. However, conflicts between observations and predictions are increasingly accumulating and little progress has been made in explaining this disparity. We argue that a revision of the concept of marginality is necessary, together with explicit testing of population decline, which is increasingly possible as data availability improves. Such action should be based on taking the population perspective across a species' rear edge, encompassing the ecological, geographical and genetic dimensions of marginality. Refining our understanding of rear edge populations is essential to advance our ability to monitor, predict and plan for the impacts of environmental change on species range dynamics.

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.

Drivers of tree fecundity in pedunculate oak (Quercus robur) refugial populations at the species' southwestern range margin

The current low latitudinal range margins of many extra-tropical plant species consist of small and scattered populations that persist locally in microrefugia. It remains poorly understood how their refugial distribution affects mating patterns and reproductive success. Here we examine flower and acorn production and their determinants in refugial populations of the widespread European forest tree pedunculate oak (Quercus robur). We monitored male flower, female flower and acorn production in 159 adult trees from 12 oak stands over 2 years. We related these and derived parameters to a series of ecological and genetic predictor variables extrinsic (stand size, density and isolation as well as elevation, topography and forest cover) or intrinsic (size, phenology and several genotypic measures) to the target tree. Tree fertility was unrelated to extrinsic factors but determined by tree size, although we detected size-independent variation in reproductive investment. Female flower number accurately predicted acorn crop size. Fruit set differed between years, evidencing the existence of pollen limitation at the landscape but not at the local scale. Fruit set also tended to increase with the number of mates of the target tree. We detected no other evidence for genetic constraints on mating. Reproduction was triggered by a combination of small-scale and landscape-scale drivers. Although short-distance mating prevailed, limited pollen flow did not appear to significantly constrain reproductive success. The high intrinsic ability of populations to maintain their reproductive capacity may help explain their successful long-term persistence in an adverse broader environment.

European Forest Cover During the Past 12,000 Years: A Palynological Reconstruction Based on Modern Analogs and Remote Sensing

Characterization of land cover change in the past is fundamental to understand the evolution and present state of the Earth system, the amount of carbon and nutrient stocks in terrestrial ecosystems, and the role played by land-atmosphere interactions in influencing climate. The estimation of land cover changes using palynology is a mature field, as thousands of sites in Europe have been investigated over the last century. Nonetheless, a quantitative land cover reconstruction at a continental scale has been largely missing. Here, we present a series of maps detailing the evolution of European forest cover during last 12,000 years. Our reconstructions are based on the Modern Analog Technique (MAT): a calibration dataset is built by coupling modern pollen samples with the corresponding satellite-based forest-cover data. Fossil reconstructions are then performed by assigning to every fossil sample the average forest cover of its closest modern analogs. The occurrence of fossil pollen assemblages with no counterparts in modern vegetation represents a known limit of analog-based methods. To lessen the influence of no-analog situations, pollen taxa were converted into plant functional types prior to running the MAT algorithm. We then interpolate site-specific reconstructions for each timeslice using a four-dimensional gridding procedure to create continuous gridded maps at a continental scale. The performance of the MAT is compared against methodologically independent forest-cover reconstructions produced using the REVEALS method. MAT and REVEALS estimates are most of the time in good agreement at a trend level, yet MAT regularly underestimates the occurrence of densely forested situations, requiring the application of a bias correction procedure. The calibrated MAT-based maps draw a coherent picture of the establishment of forests in Europe in the Early Holocene with the greatest forest-cover fractions reconstructed between ∼8,500 and 6,000 calibrated years BP. This forest maximum is followed by a general decline in all parts of the continent, likely as a result of anthropogenic deforestation. The continuous spatial and temporal nature of our reconstruction, its continental coverage, and gridded format make it suitable for climate, hydrological, and biogeochemical modeling, among other uses.

Extensive recent secondary contacts between four European white oak species

Historical trajectories of tree species during the late Quaternary have been well reconstructed through genetic and palaeobotanical studies. However, many congeneric tree species are interfertile, and the timing and contribution of introgression to species divergence during their evolutionary history remains largely unknown. We quantified past and current gene flow events between four morphologically divergent oak species (Quercus petraea, Q. robur, Q. pyrenaica, Q. pubescens), by two independent inference methods: diffusion approximation to the joint frequency spectrum (∂a∂i) and approximate Bayesian computation (ABC). For each pair of species, alternative scenarios of speciation allowing gene flow over different timescales were evaluated. Analyses of 3524 single nucleotide polymorphisms (SNPs) randomly distributed in the genome, showed that these species evolved in complete isolation for most of their history, but recently came into secondary contact, probably facilitated by the most recent period of postglacial warming. We demonstrated that: there was sufficient genetic differentiation before secondary contact for the accumulation of barriers to gene flow; and current European white oak genomes are a mosaic of genes that have crossed species boundaries and genes impermeable to gene flow.