Genomic data and common garden experiments reveal climate-driven selection on ecophysiological traits in two Mediterranean oaks

Improving our knowledge of how past climate-driven selection has acted on present-day trait population divergence is essential to understand local adaptation processes and improve our predictions of evolutionary trajectories in the face of altered selection pressures resulting from climate change. In this study, we investigated signals of selection on traits related to drought tolerance and growth rates in two Mediterranean oak species (Quercus faginea and Q. lusitanica) with contrasting distribution ranges and climatic niches. We genotyped 182 individuals from 24 natural populations of the two species using restriction-site-associated DNA sequencing and conducted a thorough functional characterization in 1602 seedlings from 21 populations cultivated in common garden experiments under contrasting watering treatments. Our genomic data revealed that both Q. faginea and Q. lusitanica have very weak population genetic structure, probably as a result of high rates of pollen-mediated gene flow among populations and large effective population sizes. In contrast, common garden experiments showed evidence of climate-driven divergent selection among populations on traits related to leaf morphology, physiology and growth in both species. Overall, our study suggests that climate is an important selective factor for Mediterranean oaks and that ecophysiological traits have evolved in drought-prone environments even in a context of very high rates of gene flow among populations.

Multi-trait genetic variation in resource-use strategies and phenotypic plasticity correlates with local climate across the range of a Mediterranean oak (Quercus faginea)

Resource-use strategies are hypothesized to evolve along climatic gradients. However, our understanding of the environmental factors driving divergent evolution of resource-use strategies and the relationship between trait genetic variation and phenotypic plasticity is far from complete. Using the Mediterranean tree Quercus faginea as study system, we tested the hypothesis that a conservative resource-use strategy with increased drought tolerance and reduced phenotypic plasticity has evolved in areas with longer and more severe dry seasons. We conducted a glasshouse experiment in which we measured leaf morphological, physiological, growth and allocation traits in seedlings from 10 range-wide climatically contrasting populations, grown under two different watering treatments. Both univariate and multivariate analyses revealed a genetic gradient of resource-use strategies and phenotypic plasticity associated with provenance climate. In particular, populations from harsher (drier and colder) environments had more sclerophyllous leaves, lower growth rates, better physiological performance under dry conditions and reduced multi-trait phenotypic plasticity compared to populations from more mesic and milder environments. Our results suggest that contrasting precipitation and temperature regimes play an important role in the adaptive intraspecific evolution of multivariate phenotypes and their plasticity, resulting in coordinated morphology, physiology, growth and allometry according to alternative resource-use strategies.

Adaptive responses to temperature and precipitation variation at the early-life stages of Pinus sylvestris

Early-stage fitness variation has been seldom evaluated at broad scales in forest tree species, despite the long tradition of studying climate-driven intraspecific genetic variation. In this study, we evaluated the role of climate in driving patterns of population differentiation at early-life stages in Pinus sylvestris and explored the fitness and growth consequences of seed transfer within the species range. We monitored seedling emergence, survival and growth over a 2-yr period in a multi-site common garden experiment which included 18 European populations and spanned 25° in latitude and 1700 m in elevation. Climate-fitness functions showed that populations exhibited higher seedling survival and growth at temperatures similar to their home environment, which is consistent with local adaptation. Northern populations experienced lower survival and growth at warmer sites, contrary to previous studies on later life stages. Seed mass was higher in populations from warmer areas and was positively associated with survival and growth at more southern sites. Finally, we did not detect a survival-growth trade-off; on the contrary, bigger seedlings exhibited higher survival probabilities under most climatic conditions. In conclusion, our results reveal that contrasting temperature regimes have played an important role in driving the divergent evolution of P. sylvestris populations at early-life stages.

Measuring recent effective gene flow among large populations in Pinus sylvestris: Local pollen shedding does not preclude substantial long-distance pollen immigration

The estimation of recent gene flow rates among vast and often weakly genetically differentiated tree populations remains a great challenge. Yet, empirical information would help understanding the interaction between gene flow and local adaptation in present-day non-equilibrium forests. We investigate here recent gene flow rates between two large native Scots pine (Pinus sylvestris L.) populations in central Iberian Peninsula (Spain), which grow on contrasting edaphic conditions six kilometers apart from each other and show substantial quantitative trait divergence in common garden experiments. Using a sample of 1,200 adult and offspring chloroplast-microsatellite haplotypes and a Bayesian inference model, we estimated substantial male gametic gene flow rates (8 and 21%) between the two natural populations, and even greater estimated immigration rates (42 and 64%) from nearby plantations into the two natural populations. Our results suggest that local pollen shedding within large tree populations does not preclude long-distance pollen immigration from large external sources, supporting the role of gene flow as a homogenizing evolutionary force contributing to low molecular genetic differentiation among populations of widely distributed wind-pollinated species. Our results also indicate the high potential for reproductive connectivity in large fragmented populations of wind-pollinated trees, and draw attention to a potential scenario of adaptive genetic divergence in quantitative traits under high gene flow.

Selection patterns on early-life phenotypic traits in Pinus sylvestris are associated with precipitation and temperature along a climatic gradient in Europe

Understanding the dynamics of selection is key to predicting the response of tree species to new environmental conditions in the current context of climate change. However, selection patterns acting on early recruitment stages and their climatic drivers remain largely unknown in most tree species, despite being a critical period of their life cycle. We measured phenotypic selection on Pinus sylvestris seed mass, emergence time and early growth rate over 2 yr in four common garden experiments established along the latitudinal gradient of the species in Europe. Significant phenotypic plasticity and among-population genetic variation were found for all measured phenotypic traits. Heat and drought negatively affected fitness in the southern sites, but heavy rainfalls also decreased early survival in middle latitudes. Climate-driven directional selection was found for higher seed mass and earlier emergence time, while the form of selection on seedling growth rates differed among sites and populations. Evidence of adaptive and maladaptive phenotypic plasticity was found for emergence time and early growth rate, respectively. Seed mass, emergence time and early growth rate have an adaptive role in the early stages of P. sylvestris and climate strongly influences the patterns of selection on these fitness-related traits.

Adaptive consequences of human-mediated introgression for indigenous tree species: the case of a relict Pinus pinaster population

Human-induced gene movement via afforestation and restoration programs is a widespread phenomenon throughout the world. However, its effects on the genetic composition of native populations have received relatively little attention, particularly in forest trees. Here, we examine to what extent gene flow from allochthonous plantations of Pinus pinaster Aiton impacts offspring performance in a neighboring relict natural population and discuss the potential consequences for the long-term genetic composition of the latter. Specifically, we conducted a greenhouse experiment involving two contrasting watering treatments to test for differences in a set of functional traits and mortality rates between P. pinaster progenies from three different parental origins: (i) local native parents, (ii) exotic parents and (iii) intercrosses between local mothers and exotic fathers (intraspecific hybrids). Our results showed differences among crosses in cumulative mortality over time: seedlings of exotic parents exhibited the lowest mortality rates and seedlings of local origin the highest, while intraspecific hybrids exhibited an intermediate response. Linear regressions showed that seedlings with higher water-use efficiency (WUE, δ(13)C) were more likely to survive under drought stress, consistent with previous findings suggesting that WUE has an important role under dry conditions in this species. However, differences in mortality among crosses were only partially explained by WUE. Other non-measured traits and factors such as inbreeding depression in the relict population are more likely to explain the lower performance of native progenies. Overall, our results indicated that intraspecific hybrids and exotic individuals are more likely to survive under stressful conditions than local native individuals, at least during the first year of development. Since summer drought is the most important demographic and selective filter affecting tree establishment in Mediterranean ecosystems, a potential early selective advantage of exotic and hybrid genotypes would enhance initial steps of introgression of non-native genes into the study relict population of P. pinaster.

Patterns of pollen dispersal in a small population of the Canarian endemic palm (Phoenix canariensis)

The genetic diversity of small populations is greatly influenced by local dispersal patterns and genetic connectivity among populations, with pollen dispersal being the major component of gene flow in many plants species. Patterns of pollen dispersal, mating system parameters and spatial genetic structure were investigated in a small isolated population of the emblematic palm Phoenix canariensis in Gran Canaria island (Canary Islands). All adult palms present in the study population (n=182), as well as 616 seeds collected from 22 female palms, were mapped and genotyped at 8 microsatellite loci. Mating system analysis revealed an average of 5.8 effective pollen donors (Nep) per female. There was strong variation in correlated paternity rates across maternal progenies (ranging from null to 0.9) that could not be explained by the location and density of local males around focal females. Paternity analysis revealed a mean effective pollen dispersal distance of ∼71 m, with ∼70% of effective pollen originating from a distance of <75 m, and 90% from <200 m. A spatially explicit mating model indicated a leptokurtic pollen dispersal kernel, significant pollen immigration (12%) from external palm groves and a directional pollen dispersal pattern that seems consistent with local altitudinal air movement. No evidence of inbreeding or genetic diversity erosion was found, but spatial genetic structure was detected in the small palm population. Overall, the results suggest substantial pollen dispersal over the studied population, genetic connectivity among different palm groves and some resilience to neutral genetic erosion and subsequently to fragmentation.

Effects of forest plantations on the genetic composition of conspecific native Aleppo pine populations

Afforestation is a common and widespread management practice throughout the world, yet its implications for the genetic diversity of native populations are still poorly understood. We examined the effect of Aleppo pine (Pinus halepensis) plantations on the genetic composition of nearby conspecific native populations. We focused on two native populations in Israel with different levels of isolation from the surrounding plantations and compared the genetic diversity of naturally established young trees within the native populations with that of local native adults, using nine nuclear microsatellite markers. We found that the genetic composition of the recruits was significantly different from that of local adults in both populations, with allelic frequency changes between generations that could not be ascribed to random drift, but rather to substantial gene flow from the surrounding planted Aleppo pine populations. The more isolated population experienced a lower gene-flow level (22%) than the less isolated population (49%). The genetic divergence between native populations at the adult-tree stage (F(st) = 0.32) was more than twice as high as that of the young trees naturally established around native adults (F(st) = 0.15). Our findings provide evidence for a rapid genetic homogenization process of native populations following the massive planting efforts in the last decades. These findings have important implications for forest management and nature conservation and constitute a warning sign for the risk of translocation of biota for local biodiversity.

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.

Wind pollination over mesoscale distances: an investigation with Scots pine

There is a gap between the order of magnitude of maximum documented distances of airborne tree pollen transport (up to 10(2)-10(3) km) and effective wind pollination (up to 10(1) km), which may partly derive from greater difficulties in detecting the latter. This study aims to assess wind pollination over scales closer to the maximum observed physical pollen transport distances. The origin of effective pollen immigrants into a strongly isolated Iberian Pinus sylvestris remnant was investigated using paternally inherited microsatellite markers and maximum-likelihood estimation combined with Monte Carlo assessment of parameter uncertainty. The results revealed significant effective pollen flow (up to 4.4%) from a large population located 100 km away, suggesting that the well-known mesoscale airborne transport of viable pine pollen can result in successful pollination over larger scales than previously reported for wind-pollinated tree species. This study supports the view that the gap between documented potential and effective wind pollen dispersal scales might not accurately reflect biological reality. Expanding the expected range of effective wind pollination has an impact on the assessment of a wide range of ecological and evolutionary processes, including reproductive assurance on fragmentation or colonization, metapopulation connectivity and interactions with local adaptation in heterogeneous habitats.

Climatic niche and neutral genetic diversity of the six Iberian pine species: a retrospective and prospective view

Quaternary climatic fluctuations have left contrasting historical footprints on the neutral genetic diversity patterns of existing populations of different tree species. We should expect the demography, and consequently the neutral genetic structure, of taxa less tolerant to particular climatic extremes to be more sensitive to long-term climate fluctuations. We explore this hypothesis here by sampling all six pine species found in the Iberian Peninsula (2464 individuals, 105 populations), using a common set of chloroplast microsatellite markers, and by looking at the association between neutral genetic diversity and species-specific climatic requirements. We found large variation in neutral genetic diversity and structure among Iberian pines, with cold-enduring mountain species (Pinus uncinata, P. sylvestris and P. nigra) showing substantially greater diversity than thermophilous taxa (P. pinea and P. halepensis). Within species, we observed a significant positive correlation between population genetic diversity and summer precipitation for some of the mountain pines. The observed pattern is consistent with the hypotheses that: (i) more thermophilous species have been subjected to stronger demographic fluctuations in the past, as a consequence of their maladaptation to recurrent glacial cold stages; and (ii) altitudinal migrations have allowed the maintenance of large effective population sizes and genetic variation in cold-tolerant species, especially in more humid regions. In the light of these results and hypotheses, we discuss some potential genetic consequences of impending climate change.

Isolation by distance in a continuous population under stochastic demographic fluctuations

The local density of individuals is seldom uniform in space and time within natural populations. Yet, formal approaches to the process of isolation by distance in continuous populations have encountered analytical difficulties in describing genetic structuring with demographic heterogeneities, usually disregarding local correlations in the movement and reproduction of genes. We formulate exact recursions for probabilities of identity in continuous populations, from which we deduce definitions of effective dispersal () and effective density (D(e)) that generalize results relating spatial genetic structure, dispersal and density in lattice models. The latter claim is checked in simulations where estimates of effective parameters obtained from demographic information are compared with estimates derived from spatial genetic patterns in a plant population evolving in a heterogeneous and dynamic habitat. The simulations further suggest that increasing spatio-temporal correlations in local density reduce and generally decrease the product , with dispersal kurtosis influencing their sensitivity to density fluctuations. As in the lattice model, the expected relationship between the product and the genetic structure statistic a(r) holds under fluctuating density, irrespective of dispersal kurtosis. The product D sigma(2) between observed census density and the observed dispersal rate over one generation will generally be an upwardly biased (up to 400% in simulations) estimator of in populations distributed in spatially aggregated habitats.

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.

Increased selfing and correlated paternity in a small population of a predominantly outcrossing conifer, Pinus sylvestris

Outcrossing rate, the rates of ovule and seed abortion, and levels of correlated paternity were estimated in a small population of Pinus sylvestris, a predominantly outcrossing conifer, and were compared with estimates from two widely dispersed woodlands of the same species, showing a range of densities. On average, seed trees of the small population showed an eight-fold higher selfing rate (25 vs. 3%) and a 100-fold greater incidence of correlated paternity (19.6 vs. 0.2%) than did trees from the large populations. No evidence was found of pollen limitation within the remnant stand, as suggested by ovule abortion rates. Investigation of the mating patterns in the small population, based on the unambiguous genealogy of 778 open-pollinated seeds, showed a large departure from random mating. Only 8% of the possible mating pairs within the stand were observed. Correlated paternity rate within a maternal sibship was negatively associated (rs = -0.398, P < 0.050) with the distance to the nearest neighbour, and shared paternity among maternal sibships was negatively correlated (rs = -0.704, P < 0.001) with the distance between seed trees. Numerical simulations, based on the estimated individual pollen dispersal kernel, suggest that restricted dispersal might have been the key factor affecting mating patterns in the small population and, together with low population density, may account for the observed mating system variation between the small and the large populations. The results of this study show that a severe size reduction may substantially affect the mating system of a wind-pollinated, typically outcrossed plant species.

Patterns of pollen dispersal in a small population of Pinus sylvestris L. revealed by total-exclusion paternity analysis

Patterns of pollen dispersal were investigated in a small, isolated, relict population of Pinus sylvestris L., consisting of 36 trees. A total-exclusion battery comprising four chloroplast and two nuclear microsatellites (theoretical paternity exclusion probability EP=0.996) was used to assign paternity to 813 seeds, collected from 34 trees in the stand. Long-distance pollen immigration accounted for 4.3% of observed matings. Self-fertilization rate was very high (0.25), compared with typical values in more widespread populations of the species. The average effective pollen dispersal distance within the stand was 48 m (or 83 m excluding selfs). Half of effective pollen was dispersed within 11 m, and 7% beyond 200 m. A strong correlation was found between the distance to the closest tree and the mean mating-distance calculated for single-tree progenies. The effective pollen dispersal distribution showed a leptokurtic shape, with a large and significant departure from that expected under uniform dispersal. A maximum-likelihood procedure was used to fit an individual pollen dispersal distance probability density function (dispersal kernel). The estimated kernel indicated fairly leptokurtic dispersal (shape parameter b=0.67), with an average pollen dispersal distance of 135 m, and 50% of pollen dispersed beyond 30 m. A marked directionality pattern of pollen dispersal was found, mainly caused by the uneven distribution of trees, coupled with restricted dispersal and unequal male success. Overall, results show that the number and distribution of potential pollen donors in small populations may strongly influence the patterns of effective pollen dispersal.

Cross-amplification and sequence variation of microsatellite loci in Eurasian hard pines

Microsatellite transfer across coniferous species is a valued methodology because de novo development for each species is costly and there are many species with only a limited commodity value. Cross-species amplification of orthologous microsatellite regions provides valuable information on mutational and evolutionary processes affecting these loci. We tested 19 nuclear microsatellite markers from Pinus taeda L. (subsection Australes) and three from P. sylvestris L. (subsection Pinus) on seven Eurasian hard pine species ( P. uncinata Ram., P. sylvestris L., P. nigra Arn., P. pinaster Ait., P. halepensis Mill., P. pinea L. and P. canariensis Sm.). Transfer rates to species in subsection Pinus (36-59%) were slightly higher than those to subsections Pineae and Pinaster (32-45%). Half of the trans-specific microsatellites were found to be polymorphic over evolutionary times of approximately 100 million years (ten million generations). Sequencing of three trans-specific microsatellites showed conserved repeat and flanking regions. Both a decrease in the number of perfect repeats in the non-focal species and a polarity for mutation, the latter defined as a higher substitution rate in the flanking sequence regions close to the repeat motifs, were observed in the trans-specific microsatellites. The transfer of microsatellites among hard pine species proved to be useful for obtaining highly polymorphic markers in a wide range of species, thereby providing new tools for population and quantitative genetic studies.

Cross-amplification and sequence variation of microsatellite loci in Eurasian hard pines

Microsatellite transfer across coniferous species is a valued methodology because de novo development for each species is costly and there are many species with only a limited commodity value. Cross-species amplification of orthologous microsatellite regions provides valuable information on mutational and evolutionary processes affecting these loci. We tested 19 nuclear microsatellite markers from Pinus taeda L. (subsection Australes) and three from P. sylvestris L. (subsection Pinus) on seven Eurasian hard pine species ( P. uncinata Ram., P. sylvestris L., P. nigra Arn., P. pinaster Ait., P. halepensis Mill., P. pinea L. and P. canariensis Sm.). Transfer rates to species in subsection Pinus (36-59%) were slightly higher than those to subsections Pineae and Pinaster (32-45%). Half of the trans-specific microsatellites were found to be polymorphic over evolutionary times of approximately 100 million years (ten million generations). Sequencing of three trans-specific microsatellites showed conserved repeat and flanking regions. Both a decrease in the number of perfect repeats in the non-focal species and a polarity for mutation, the latter defined as a higher substitution rate in the flanking sequence regions close to the repeat motifs, were observed in the trans-specific microsatellites. The transfer of microsatellites among hard pine species proved to be useful for obtaining highly polymorphic markers in a wide range of species, thereby providing new tools for population and quantitative genetic studies.