Ecological determinants of mating system within and between three Fagus sylvatica populations along an elevational gradient

Propensity for selfing varies within a population of hermaphroditic snails: coexistence of selfers, outcrossers and mixed-mating individuals

To understand mating-system evolution in self-compatible hermaphrodites, variation in selfing rates is highly relevant. Empirical studies are rarely designed to capture variation between individuals, instead often comparing species and populations. Yet, evolution primarily occurs within populations, rendering among-individual variation essential. Observed individual selfing rates depend on the environment (e.g. differences in mate availability) and individuals' propensity for selfing. We quantified individual variation in selfing propensity in the snail Radix balthica by conducting laboratory mating trials that manipulated mate availability (low versus moderate) and estimating selfing rates from progeny arrays. We also measured female lifetime fitness. We found substantial among-individual variation in selfing propensity, including pure selfers (32%), pure outcrossers (31%) and mixed-mating individuals that selfed and outcrossed (37%). Experimental levels of mate availability did not significantly affect selfing rates. Selfers had reduced female liftetime fitness. Our results show that the propensity for selfing can differ considerably among individuals, with similar proportions of selfers, outcrossers and mixed maters. As mate availability did not affect selfing, our 'moderate' experimental level of mate availability might still have been too low to prompt selfers to outcross. This and the observed fitness differences also cautiously suggest that investigating the heritability of selfing propensities might be worthwhile in this population.

A novel synthesis of two decades of microsatellite studies on European beech reveals decreasing genetic diversity from glacial refugia

Genetic diversity influences the evolutionary potential of forest trees under changing environmental conditions, thus indirectly the ecosystem services that forests provide. European beech (Fagus sylvatica L.) is a dominant European forest tree species that increasingly suffers from climate change-related die-back. Here, we conducted a systematic literature review of neutral genetic diversity in European beech and created a meta-data set of expected heterozygosity (He) from all past studies providing nuclear microsatellite data. We propose a novel approach, based on population genetic theory and a min-max scaling to make past studies comparable. Using a new microsatellite data set with unprecedented geographic coverage and various re-sampling schemes to mimic common sampling biases, we show the potential and limitations of the scaling approach. The scaled meta-dataset reveals the expected trend of decreasing genetic diversity from glacial refugia across the species range and also supports the hypothesis that different lineages met and admixed north of the European mountain ranges. As a result, we present a map of genetic diversity across the range of European beech which could help to identify seed source populations harboring greater diversity and guide sampling strategies for future genome-wide and functional investigations of genetic variation. Our approach illustrates how to combine information from several nuclear microsatellite data sets to describe patterns of genetic diversity extending beyond the geographic scale or mean number of loci used in each individual study, and thus is a proof-of-concept for synthesizing knowledge from existing studies also in other species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11295-022-01577-4.

Rapid divergence in vegetative morphology of a wind-pollinated plant between populations at contrasting densities

Plant sexual dimorphism is thought to evolve in response to sex-specific selection associated with competition for access to mates or resources, both of which may be density dependent. In wind-pollinated plants in particular, vegetative traits such as plant size and architecture may influence resource acquisition and both pollen dispersal and receipt, with potential conflict between these two components of fitness. We evaluated the role of plant density in shaping plant traits by measuring evolutionary responses in experimental populations of the sexually dimorphic wind-pollinated plant Mercurialis annua. After three generations of evolution, we observed divergence between high- and low-density populations in several vegetative traits, whereas there was no divergence for reproductive traits. A reversal in the direction of sexually dimorphic traits expressed in young plants evolved in both low- and high-density populations compared to the original population (stored as seeds). Compared to the source population, males at high density evolved to be taller when young, whereas at low density young females tended to become smaller. These results demonstrate that a simple change in plant density can induce age-dependent and sex-specific evolution in the ontogeny of vegetative organs, and illustrates the power of experimental evolution for investigating plant trait evolution.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

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.

Sex-specific spatial variation in fitness in the highly dimorphic Leucadendron rubrum

Sexual dimorphism in plants may emerge as a result of sex-specific selection on traits enhancing access to nutritive resources and/or to sexual partners. Here we investigated sex-specific differences in selection of sexually dimorphic traits and in the spatial distribution of effective fecundity (our fitness proxy) in a highly dimorphic dioecious wind-pollinated shrub, Leucadendron rubrum. In particular, we tested for the effect of density on male and female effective fecundity. We used spatial and genotypic data of parent and offspring cohorts to jointly estimate individual male and female effective fecundity on the one hand and pollen and seed dispersal kernels on the other hand. This methodology was adapted to the case of dioecious species. Explicitly modelling dispersal avoids the confounding effects of heterogeneous spatial distribution of mates and sampled seedlings on the estimation of effective fecundity. We also estimated selection gradients on plant traits while modelling sex-specific spatial autocorrelation in fecundity. Males exhibited spatial autocorrelation in effective fecundity at a smaller scale than females. A higher local density of plants was associated with lower effective fecundity in males but was not related to female effective fecundity. These results suggest sex-specific sensitivities to environmental heterogeneity in L. rubrum. Despite these sexual differences, we found directional selection for wider canopies and smaller leaves in both sexes, and no sexually antagonistic selection on strongly dimorphic traits in L. rubrum. Many empirical studies in animals similarly failed to detect sexually antagonistic selection in species expressing strong sexual dimorphism, and we discuss reasons explaining this common pattern.

Weak founder effects but significant spatial genetic imprint of recent contraction and expansion of European beech populations

Understanding the ecological and evolutionary processes occurring during species range shifts is important in the current context of global change. Here, we investigate the interplay between recent expansion, gene flow and genetic drift, and their consequences for genetic diversity and structure at landscape and local scales in European beech (Fagus sylvatica L.) On Mont Ventoux, South-Eastern France, we located beech forest refugia at the time of the most recent population minimum, ~150 years ago, and sampled 71 populations (2042 trees) in both refugia and expanding populations over an area of 15,000 ha. We inferred patterns of gene flow and genetic structure using 12 microsatellite markers. We identified six plots as originating from planting, rather than natural establishment, mostly from local genetic material. Comparing genetic diversity and structure in refugia versus recent populations did not support the existence of founder effects: heterozygosity (He = 0.667) and allelic richness (Ar = 4.298) were similar, and FST was low (0.031 overall). Still, significant spatial evidence of colonization was detected, with He increasing along the expansion front, while genetic differentiation from the entire pool (βWT) decreased. Isolation by distance was found in refugia but not in recently expanding populations. Our study indicates that beech capacities for colonization and gene flow were sufficient to preserve genetic diversity despite recent forest contraction and expansion. Because beech has long distance pollen and seed dispersal, these results illustrate a 'best case scenario' for the maintenance of high genetic diversity and adaptive potential under climate-change-related range change.

Microgeographic adaptation and the effect of pollen flow on the adaptive potential of a temperate tree species

In species with long-distance dispersal capacities and inhabiting a large ecological niche, local selection and gene flow are expected to be major evolutionary forces affecting the genetic adaptation of natural populations. Yet, in species such as trees, evidence of microgeographic adaptation and the quantitative assessment of the impact of gene flow on adaptive genetic variation are still limited. Here, we used extensive genetic and phenotypic data from European beech seedlings collected along an elevation gradient, and grown in a common garden, to study the signature of selection on the divergence of eleven potentially adaptive traits, and to assess the role of gene flow in resupplying adaptive genetic variation. We found a significant signal of adaptive differentiation among plots separated by < 1 km, with selection acting on growth and phenological traits. Consistent with theoretical expectations, our results suggest that pollen dispersal contributes to increase genetic diversity for these locally differentiated traits. Our results thus highlight that local selection is an important evolutionary force in natural tree populations and suggest that management interventions to facilitate movement of gametes along short ecological gradients would boost genetic diversity of individual tree populations, and enhance their adaptive potential to rapidly changing environments.

Elevation, Not Deforestation, Promotes Genetic Differentiation in a Pioneer Tropical Tree

The regeneration of disturbed forest is an essential part of tropical forest ecology, both with respect to natural disturbance regimes and large-scale human-mediated logging, grazing, and agriculture. Pioneer tree species are critical for facilitating the transition from deforested land to secondary forest because they stabilize terrain and enhance connectivity between forest fragments by increasing matrix permeability and initiating disperser community assembly. Despite the ecological importance of early successional species, little is known about their ability to maintain gene flow across deforested landscapes. Utilizing highly polymorphic microsatellite markers, we examined patterns of genetic diversity and differentiation for the pioneer understory tree Miconia affinis across the Isthmus of Panama. Furthermore, we investigated the impact of geographic distance, forest cover, and elevation on genetic differentiation among populations using circuit theory and regression modeling within a landscape genetics framework. We report marked differences in historical and contemporary migration rates and moderately high levels of genetic differentiation in M. affinis populations across the Isthmus of Panama. Genetic differentiation increased significantly with elevation and geographic distance among populations; however, we did not find that forest cover enhanced or reduced genetic differentiation in the study region. Overall, our results reveal strong dispersal for M. affinis across human-altered landscapes, highlighting the potential use of this species for reforestation in tropical regions. Additionally, this study demonstrates the importance of considering topography when designing programs aimed at conserving genetic diversity within degraded tropical landscapes.

Genome-environment association study suggests local adaptation to climate at the regional scale in Fagus sylvatica

The evolutionary potential of long-lived species, such as forest trees, is fundamental for their local persistence under climate change (CC). Genome-environment association (GEA) analyses reveal if species in heterogeneous environments at the regional scale are under differential selection resulting in populations with potential preadaptation to CC within this area. In 79 natural Fagus sylvatica populations, neutral genetic patterns were characterized using 12 simple sequence repeat (SSR) markers, and genomic variation (144 single nucleotide polymorphisms (SNPs) out of 52 candidate genes) was related to 87 environmental predictors in the latent factor mixed model, logistic regressions and isolation by distance/environmental (IBD/IBE) tests. SSR diversity revealed relatedness at up to 150 m intertree distance but an absence of large-scale spatial genetic structure and IBE. In the GEA analyses, 16 SNPs in 10 genes responded to one or several environmental predictors and IBE, corrected for IBD, was confirmed. The GEA often reflected the proposed gene functions, including indications for adaptation to water availability and temperature. Genomic divergence and the lack of large-scale neutral genetic patterns suggest that gene flow allows the spread of advantageous alleles in adaptive genes. Thereby, adaptation processes are likely to take place in species occurring in heterogeneous environments, which might reduce their regional extinction risk under CC.

Elevation as a barrier: genetic structure for an Atlantic rain forest tree (Bathysa australis) in the Serra do Mar mountain range, SE Brazil

Distance and discrete geographic barriers play a role in isolating populations, as seed and pollen dispersal become limited. Nearby populations without any geographic barrier between them may also suffer from ecological isolation driven by habitat heterogeneity, which may promote divergence by local adaptation and drift. Likewise, elevation gradients may influence the genetic structure and diversity of populations, particularly those marginally distributed. Bathysa australis (Rubiaceae) is a widespread tree along the elevation gradient of the Serra do Mar, SE Brazil. This self-compatible species is pollinated by bees and wasps and has autochoric seeds, suggesting restricted gene dispersal. We investigated the distribution of genetic diversity in six B. australis populations at two extreme sites along an elevation gradient: a lowland site (80–216 m) and an upland site (1010–1100 m.a.s.l.). Nine microsatellite loci were used to test for genetic structure and to verify differences in genetic diversity between sites. We found a marked genetic structure on a scale as small as 6 km (F_ST = 0.21), and two distinct clusters were identified, each corresponding to a site. Although B. australis is continuously distributed along the elevation gradient, we have not observed a gene flow between the extreme populations. This might be related to B. australis biological features and creates a potential scenario for adaptation to the different conditions imposed by the elevation gradient. We failed to find an isolation-by-distance pattern; although on the fine scale, all populations showed spatial autocorrelation until ∼10-20 m. Elevation difference was a relevant factor though, but we need further sampling effort to check its correlation with genetic distance. The lowland populations had a higher allelic richness and showed higher rare allele counts than the upland ones. The upland site may be more selective, eliminating rare alleles, as we did not find any evidence for bottleneck.

Effects of tree architecture on pollen dispersal and mating patterns in Abies pinsapo Boiss. (Pinaceae)

Plant architecture is crucial to pollination and mating in wind-pollinated species. We investigated the effect of crown architecture on pollen dispersal, mating system and offspring quality, combining phenotypic and genotypic analyses in a low-density population of the endangered species Abies pinsapo. A total of 598 embryos from three relative crown height levels (bottom, middle and top) in five mother plants were genotyped using eleven nuclear microsatellite markers (nSSRs). Paternity analysis and mating system models were used to infer mating and pollen dispersal parameters. In addition, seeds were weighed (N = 16 110) and germinated (N = 736), and seedling vigour was measured to assess inbreeding depression. Overall, A. pinsapo shows a fat-tailed dispersal kernel, with an average pollen dispersal distance of 113-227 m, an immigration rate of 0.84-26.92%, and a number of effective pollen donors (Nep ) ranging between 3.5 and 11.9. We found an effect of tree height and relative crown height levels on mating parameters. A higher proportion of seeds with embryo (about 50%) and a higher rate of self-fertilization (about 60%) were found at the bottom level in comparison with the top level. Seed weight and seedling vigour are positively related. Nevertheless, no differences were found in seed weight or in seedling-related variables such as weight and length of aerial and subterranean parts among the different relative crown height levels, suggesting that seeds from the more strongly inbred bottom level are not affected by inbreeding depression. Our results point to vertical isotropy for outcross-pollen and they suggest that self-pollen may ensure fertilization when outcross-pollen is not available in low-density population.

Detecting short spatial scale local adaptation and epistatic selection in climate-related candidate genes in European beech (Fagus sylvatica) populations

Detecting signatures of selection in tree populations threatened by climate change is currently a major research priority. Here, we investigated the signature of local adaptation over a short spatial scale using 96 European beech (Fagus sylvatica L.) individuals originating from two pairs of populations on the northern and southern slopes of Mont Ventoux (south-eastern France). We performed both single and multilocus analysis of selection based on 53 climate-related candidate genes containing 546 SNPs. FST outlier methods at the SNP level revealed a weak signal of selection, with three marginally significant outliers in the northern populations. At the gene level, considering haplotypes as alleles, two additional marginally significant outliers were detected, one on each slope. To account for the uncertainty of haplotype inference, we averaged the Bayes factors over many possible phase reconstructions. Epistatic selection offers a realistic multilocus model of selection in natural populations. Here, we used a test suggested by Ohta based on the decomposition of the variance of linkage disequilibrium. Overall populations, 0.23% of the SNP pairs (haplotypes) showed evidence of epistatic selection, with nearly 80% of them being within genes. One of the between gene epistatic selection signals arose between an FST outlier and a nonsynonymous mutation in a drought response gene. Additionally, we identified haplotypes containing selectively advantageous allele combinations which were unique to high or low elevations and northern or southern populations. Several haplotypes contained nonsynonymous mutations situated in genes with known functional importance for adaptation to climatic factors.

Pollen, wind and fire: how to investigate genetic effects of disturbance-induced change in forest trees

Understanding the consequences of habitat disturbance on mating patterns although pollen and seed dispersal in forest trees has been a long-standing theme of forest and conservation genetics. Forest ecosystems face global environmental pressures from timber exploitation to genetic pollution and climate change, and it is therefore essential to comprehend how disturbances may alter the dispersal of genes and their establishment in tree populations in order to formulate relevant recommendations for sustainable resource management practices and realistic predictions of potential adaptation to climate change by means of range shift or expansion (Kremer et al. 2012). However, obtaining reliable evidence of disturbance-induced effects on gene dispersal processes from empirical evaluation of forest tree populations is difficult. Indeed, tree species share characteristics such as high longevity, long generation time and large reproductive population size, which may impede the experimenter's ability to assess parameters at the spatial and time scales at which any change may occur (Petit and Hampe 2006). It has been suggested that appropriate study designs should encompass comparison of populations before and after disturbance as well as account for demonstrated variation in conspecific density, that is, the spatial distribution of mates, and forest density, including all species and relating to alteration in landscape openness (Bacles & Jump 2011). However, more often than not, empirical studies aiming to assess the consequences of habitat disturbance on genetic processes in tree populations assume rather than quantify a change in tree densities in forests under disturbance and generally fail to account for population history, which may lead to inappropriate interpretation of a causal relationship between population genetic structure and habitat disturbance due to effects of unmonitored confounding variables (Gauzere et al. 2013). In this issue, Shohami and Nathan (2014) take advantage of the distinctive features of the fire-adapted wind-pollinated Aleppo pine Pinus halepensis (Fig. 1) to provide an elegant example of best practice. Thanks to long-term monitoring of the study site, a natural stand in Israel, Shohami and Nathan witnessed the direct impact of habitat disturbance, here taking the shape of fire, on conspecific and forest densities and compared pre- and postdisturbance mating patterns estimated from cones of different ages sampled on the same surviving maternal individuals (Fig. 2). This excellent study design is all the more strong that Shohami and Nathan took further analytical steps to account for confounding variables, such as historical population genetic structure and possible interannual variation in wind conditions, thus giving high credibility to their findings of unequivocal fire-induced alteration of mating patterns in P. halepensis. Most notably, the authors found, at the pollen pool level, a disruption of local genetic structure which, furthermore, they were able to attribute explicitly to enhanced pollen-mediated gene immigration into the low-density fire-disturbed stand. This cleverly designed research provides a model approach to be followed if we are to advance our understanding of disturbance-induced dispersal and genetic change in forest trees.

Space, time and complexity in plant dispersal ecology

Dispersal of pollen and seeds are essential functions of plant species, with far-reaching demographic, ecological and evolutionary consequences. Interest in plant dispersal has increased with concerns about the persistence of populations and species under global change. We argue here that advances in plant dispersal ecology research will be determined by our ability to surmount challenges of spatiotemporal scales and heterogeneities and ecosystem complexity. Based on this framework, we propose a selected set of research questions, for which we suggest some specific objectives and methodological approaches. Reviewed topics include multiple vector contributions to plant dispersal, landscape-dependent dispersal patterns, long-distance dispersal events, spatiotemporal variation in dispersal, and the consequences of dispersal for plant communities, populations under climate change, and anthropogenic landscapes.