The GenTree Platform: growth traits and tree-level environmental data in 12 European forest tree species.

BACKGROUND: Progress in the field of evolutionary forest ecology has been hampered by the huge challenge of phenotyping trees across their ranges in their natural environments, and the limitation in high-resolution environmental information. FINDINGS: The GenTree Platform contains phenotypic and environmental data from 4,959 trees from 12 ecologically and economically important European forest tree species: Abies alba Mill. (silver fir), Betula pendula Roth. (silver birch), Fagus sylvatica L. (European beech), Picea abies (L.) H. Karst (Norway spruce), Pinus cembra L. (Swiss stone pine), Pinus halepensis Mill. (Aleppo pine), Pinus nigra Arnold (European black pine), Pinus pinaster Aiton (maritime pine), Pinus sylvestris L. (Scots pine), Populus nigra L. (European black poplar), Taxus baccata L. (English yew), and Quercus petraea (Matt.) Liebl. (sessile oak). Phenotypic (height, diameter at breast height, crown size, bark thickness, biomass, straightness, forking, branch angle, fructification), regeneration, environmental in situ measurements (soil depth, vegetation cover, competition indices), and environmental modeling data extracted by using bilinear interpolation accounting for surrounding conditions of each tree (precipitation, temperature, insolation, drought indices) were obtained from trees in 194 sites covering the species' geographic ranges and reflecting local environmental gradients. CONCLUSION: The GenTree Platform is a new resource for investigating ecological and evolutionary processes in forest trees. The coherent phenotyping and environmental characterization across 12 species in their European ranges allow for a wide range of analyses from forest ecologists, conservationists, and macro-ecologists. Also, the data here presented can be linked to the GenTree Dendroecological collection, the GenTree Leaf Trait collection, and the GenTree Genomic collection presented elsewhere, which together build the largest evolutionary forest ecology data collection available.

How does contemporary selection shape oak phenotypes?

Most existing forests are subjected to natural and human-mediated selection pressures, which have increased due to climate change and the increasing needs of human societies for wood, fibre and fuel resources. It remains largely unknown how these pressures trigger evolutionary changes. We address this issue here for temperate European oaks (Quercus petraea and Q. robur), which grow in mixed stands, under even-aged management regimes. We screened numerous functional traits for univariate selection gradients and for expected and observed genetic changes over two successive generations. In both species, growth, leaf morphology and physiology, and defence-related traits displayed significant selection gradients and predicted shifts, whereas phenology, water metabolism, structure and resilience-related traits did not. However, the direction of the selection response and the potential for adaptive evolution differed between the two species. Quercus petraea had a much larger phenotypic and genetic variance of fitness than Q. robur. This difference raises concerns about the adaptive response of Q. robur to contemporary selection pressures. Our investigations suggest that Q. robur will probably decline steadily, particularly in mixed stands with Q. petraea, consistent with the contrasting demographic dynamics of the two species.

Which oak provenances for the 22nd century in Western Europe? Dendroclimatology in common gardens.

The current distribution area of the two sympatric oaks Quercus petraea and Q. robur covers most of temperate Western Europe. Depending on their geographic location, populations of these trees are exposed to different climate constraints, to which they are adapted. Comparing the performances of trees from contrasting populations provides the insight into their expected resilience to future climate change required for forest management. In this study, the descendants of 24 Q. petraea and two Q. robur provenances selected from sites throughout Europe were grown for 20 years in three common gardens with contrasting climates. The 2420 sampled trees allowed the assessments of the relationship between radial growth and climate. An analysis of 15-year chronologies of ring widths, with different combinations of climate variables, revealed different response patterns between provenances and between common gardens. As expected, provenances originating from sites with wet summers displayed the strongest responses to summer drought, particularly in the driest common garden. All provenances displayed positive significant relationships between the temperature of the previous winter and radial growth when grown in the common garden experiencing the mildest winter temperatures. Only eastern provenances from continental cold climates also clearly expressed this limitation of growth by cold winter temperatures in the other two common gardens. However, ecological distance, calculated on the basis of differences in climate between the site of origin and the common garden, was not clearly related to the radial growth responses of the provenances. This suggests that the gradient of genetic variability among the selected provenances was not strictly structured according to climate gradients. Based on these results, we provide guidelines for forest managers for the assisted migration of Quercus petraea and Q. robur provenances.

In situ estimation of genetic variation of functional and ecological traits in Quercus petraea and Q.robur.

BACKGROUND: Predicting the evolutionary potential of natural tree populations requires the estimation of heritability and genetic correlations among traits on which selection acts, as differences in evolutionary success between species may rely on differences for these genetic parameters. In situ estimates are expected to be more accurate than measures done under controlled conditions which do not reflect the natural environmental variance. AIMS: The aim of the current study was to estimate three genetic parameters (i.e. heritability, evolvability and genetic correlations) in a natural mixed oak stand composed of Quercus petraea and Quercus robur about 100 years old, for 58 traits of ecological and functional relevance (growth, reproduction, phenology, physiology, resilience, structure, morphology and defence). METHODS: First we estimated genetic parameters directly in situ using realized genomic relatedness of adult trees and parentage relationships over two generations to estimate the traits additive variance. Secondly, we benefited from existing ex situ experiments (progeny tests and conservation collection) installed with the same populations, thus allowing comparisons of in situ heritability estimates with more traditional methods. RESULTS: Heritability and evolvability estimates obtained with different methods varied substantially and showed large confidence intervals, however we found that in situ were less precise than ex situ estimates, and assessments over two generations (with deeper relatedness) improved estimates of heritability while large sampling sizes are needed for accurate estimations. At the biological level, heritability values varied moderately across different ecological and functional categories of traits, and genetic correlations among traits were conserved over the two species. CONCLUSION: We identified limits for using realized genomic relatedness in natural stands to estimate the genetic variance, given the overall low variance of genetic relatedness and the rather low sampling sizes of currently used long term genetic plots in forestry. These limits can be overcome if larger sample sizes are considered, or if the approach is extended over the next generation.

Genetic differentiation in functional traits among European sessile oak populations.

The vulnerability of forest species and tree populations to climate change is related to the exposure of the ecosystem to extreme climatic conditions and to the adaptive capacity of the population to cope with those conditions. Adaptive capacity is a relatively under-researched topic within the forest science community, and there is an urgent need to understand to what extent particular combinations of traits have been shaped by natural selection under climatic gradients, potentially resulting in adaptive multi-trait associations. Thus, our aim was to quantify genetic variation in several leaf and woody traits that may contribute to multi-trait associations in which intra-specific variation could represent a source for species adaptation to climate change. A multi-trait approach was performed using nine Quercus petraea provenances originating from different locations that cover most of the species' distribution range over Europe and that were grown in a common garden. Multiple adaptive differences were observed between oak provenances but also some evolutionary stasis. In addition, our results revealed higher genetic differentiation in traits related to phenology and growth than in those related to xylem anatomy, physiology and hydraulics, for which no genetic differentiation was observed. The multiple associations between those traits and climate variables resulting from multivariate and path analyses suggest a multi-trait association largely involving phenological and growth traits for Q. petraea.

Evolutionary response to coexistence with close relatives: increased resistance against specialist herbivores without cost for climatic-stress resistance.

Why can hosts coexist with conspecifics or phylogenetically proximate neighbours despite sharing specialist enemies? Do the hosts evolve increased enemy resistance? If so, does this have costs in terms of climatic-stress resistance, or in such neighbourhoods, does climatic-stress select for resistances that are multifunctional against climate and enemies? We studied oak (Quercus petraea) descendants from provenances of contrasting phylogenetic neighbourhoods and climates in a 25-year-old common garden. We found that descendants from conspecific or phylogenetically proximate neighbourhoods had the toughest leaves and fewest leaf miners, but no reduction in climatic-stress resistance. Descendants from such neighbourhoods under cold or dry climates had the highest flavonol and anthocyanin levels and the thickest leaves. Overall, populations facing phylogenetically proximate neighbours can rapidly evolve herbivore resistance, without cost to climatic-stress resistance, but possibly facilitating resistance against cold and drought via multifunctional traits. Microevolution might hence facilitate ecological coexistence of close relatives and thereby macroevolutionary conservatism of niches.

Common garden comparisons confirm inherited differences in sensitivity to climate change between forest tree species.

The natural distribution, habitat, growth and evolutionary history of tree species are strongly dependent on ecological and genetic processes in ecosystems subject to fluctuating climatic conditions, but there have been few experimental comparisons of sensitivity between species. We compared the responses of two broadleaved tree species (Fagus sylvatica and Quercus petraea) and two conifer tree species (Pinus sylvestris and Picea abies) to climatic transfers by fitting models containing the same climatic variables. We used published data from European provenance test networks to model the responses of individual populations nested within species. A mixed model approach was applied to develop a response function for tree height over climatic transfer distance, taking into account the climatic conditions at both the seed source and the test location. The two broadleaved species had flat climatic response curves, indicating high levels of plasticity in populations, facilitating adaptation to a broader range of environments, and conferring a high potential for resilience in the face of climatic change. By contrast, the two conifer species had response curves with more pronounced slopes, indicating a lower resilience to climate change. This finding may reflect stronger genetic clines in P. sylvestris and P. abies, which constrain their climate responses to narrower climatic ranges. The response functions had maxima that deviated from the expected maximum productivity in the climate of provenance towards cooler/moister climate conditions, which we interpreted as an adaptation lag. Unilateral, linear regression analyses following transfer to warmer and drier sites confirmed a decline in productivity, predictive of the likely impact of ongoing climate change on forest populations. The responses to mimicked climate change evaluated here are of considerable interest for forestry and ecology, supporting projections of expected performance based on "real-time" field data.

Back to America: tracking the origin of European introduced populations of Quercus rubra L.

Quercus rubra has been introduced in Europe since the end of the 17th century. It is widely distributed today across this continent and considered invasive in some countries. Here, we investigated the distribution of genetic diversity of both native and introduced populations with the aim of tracing the origin of introduced populations. A large sampling of 883 individuals from 73 native and 38 European locations were genotyped at 69 SNPs. In the natural range, we found a continuous geographic gradient of variation with a predominant latitudinal component. We explored the existence of ancestral populations by performing Bayesian clustering analysis and found support for two or three ancestral genetic clusters. Approximate Bayesian Computations analyses based on these two or three clusters support recent extensive secondary contacts between them, suggesting that present-day continuous genetic variation resulted from recent admixture. In the introduced range, one main genetic cluster was not recovered in Europe, suggesting that source populations were preferentially located in the northern part of the natural distribution. However, our results cannot refute the introduction of populations from the southern states that did not survive in Europe.

Fine-scale species distribution changes in a mixed oak stand over two successive generations.

Large-scale tree distribution changes have received considerable attention but underlying demo-genetic mechanisms are less well documented. We used a diachronic approach to track species shifts in a mixed oak stand (Quercus petraea-Quercus robur) at a fine spatiotemporal scale. Species assignment was made using single nucleotide polymorphism (SNP) fingerprints employing clustering and parentage analysis. Mating patterns and reproductive success were assessed by parentage analysis. Plot-based inventories of soil parameters and sapling densities provided ecological and demographic information, respectively. Sapling density and reproductive success was higher in Q. petraea than in Q. robur, and were correlated with a spatial expansion of Q. petraea (50% to 67% of the area). Admixed trees resulting from hybridization and backcrossing between the two species were more frequent under the Q. robur canopy. We suspect that species' differential responses to ongoing environmental changes and interspecific competition are the predominant factors accounting for the recruitment success of Q. petraea, while human interference, differential reproduction and hybridization (and backcrossings) are probably of more limited importance. We anticipate in mixed Q. petraea-Q. robur stands, under current ongoing environmental change, that these processes will be enhanced, at least in the western part of the distribution of the two species.

Adaptive and plastic responses of Quercus petraea populations to climate across Europe.

How temperate forests will respond to climate change is uncertain; projections range from severe decline to increased growth. We conducted field tests of sessile oak (Quercus petraea), a widespread keystone European forest tree species, including more than 150 000 trees sourced from 116 geographically diverse populations. The tests were planted on 23 field sites in six European countries, in order to expose them to a wide range of climates, including sites reflecting future warmer and drier climates. By assessing tree height and survival, our objectives were twofold: (i) to identify the source of differential population responses to climate (genetic differentiation due to past divergent climatic selection vs. plastic responses to ongoing climate change) and (ii) to explore which climatic variables (temperature or precipitation) trigger the population responses. Tree growth and survival were modeled for contemporary climate and then projected using data from four regional climate models for years 2071-2100, using two greenhouse gas concentration trajectory scenarios each. Overall, results indicated a moderate response of tree height and survival to climate variation, with changes in dryness (either annual or during the growing season) explaining the major part of the response. While, on average, populations exhibited local adaptation, there was significant clinal population differentiation for height growth with winter temperature at the site of origin. The most moderate climate model (HIRHAM5-EC; rcp4.5) predicted minor decreases in height and survival, while the most extreme model (CCLM4-GEM2-ES; rcp8.5) predicted large decreases in survival and growth for southern and southeastern edge populations (Hungary and Turkey). Other nonmarginal populations with continental climates were predicted to be severely and negatively affected (Bercé, France), while populations at the contemporary northern limit (colder and humid maritime regions; Denmark and Norway) will probably not show large changes in growth and survival in response to climate change.

Impacts of local adaptation of forest trees on associations with herbivorous insects: implications for adaptive forest management.

Disruption of species interactions is a key issue in climate change biology. Interactions involving forest trees may be particularly vulnerable due to evolutionary rate limitations imposed by long generation times. One mitigation strategy for such impacts is Climate matching - the augmentation of local native tree populations by input from nonlocal populations currently experiencing predicted future climates. This strategy is controversial because of potential cascading impacts on locally adapted animal communities. We explored these impacts using abundance data for local native gallwasp herbivores sampled from 20 provenances of sessile oak (Quercus petraea) planted in a common garden trial. We hypothesized that non-native provenances would show (i) declining growth performance with increasing distance between provenance origin and trial site, and (ii) phenological differences to local oaks that increased with latitudinal differences between origin and trial site. Under a local adaptation hypothesis, we predicted declining gallwasp abundance with increasing phenological mismatch between native and climate-matched trees. Both hypotheses for oaks were supported. Provenance explained significant variation in gallwasp abundance, but no gall type showed the relationship between abundance and phenological mismatch predicted by a local adaptation hypothesis. Our results show that climate matching would have complex and variable impacts on oak gall communities.

Cryptic no more: soil macrofossils uncover Pleistocene forest microrefugia within a periglacial desert.

Despite their critical importance for understanding the local effects of global climate change on biodiversity, glacial microrefugia are not well studied because they are difficult to detect by using classical palaeoecological or population genetics approaches. We used soil macrofossil charcoal analysis to uncover the presence of cryptic glacial refugia for European beech (Fagus sylvatica) and other tree species in the Landes de Gascogne (southwestern France). Using botanical identification and direct radiocarbon dating (140 (14) C-dates) of macrofossil charcoal extracted from mineral soils, we reconstructed the glacial and postglacial history of all extant beech stands in the region (n = 11). Soil charcoal macrofossils were found in all sites, allowing the identification of up to at least 14 distinct fire events per site. There was direct evidence of the presence of beech during the last glacial period at three sites. Beech was detected during Heinrich stadial-1, one of the coldest and driest intervals of the last glacial period in Western Europe. Together with previous results on the genetic structure of the species in the region, these findings suggest that beech persisted in situ in several microrefugia through full glacial and interglacial periods up to the present day.

Distinct male reproductive strategies in two closely related oak species.

Reproductive strategies of closely related species distributed along successional gradients should differ as a consequence of the trade-off between competition and colonization abilities. We compared male reproductive strategies of Quercus robur and Q. petraea, two partly interfertile European oak species with different successional status. In the studied even-aged stand, trees of the late-successional species (Q. petraea) grew faster and suffered less from intertree competition than trees of the early-successional species (Q. robur). A large-scale paternity study and a spatially explicit individual-based mating model were used to estimate parameters of pollen production and dispersal as well as sexual barriers between species. Male fecundity was found to be dependent both on a tree's circumference and on its environment, particularly so for Q. petraea. Pollen dispersal was greater and more isotropic in Q. robur than in Q. petraea. Premating barriers to hybridization were strong in both species, but more so in Q. petraea than in Q. robur. Hence, predictions based on the competition-colonization trade-off are well supported, whereas the sexual barriers themselves seem to be shaped by colonization dynamics.

Stronger spatial genetic structure in recolonized areas than in refugia in the European beech.

Extant rear-edge populations located in former glacial refugia remain understudied despite their high conservation value. These populations should have experienced strong genetic drift due to their small size and long isolation. Moreover, the prolonged action of isolation by distance in refugial areas should result in stronger regional spatial genetic structure (SGS) than in recolonized areas, but empirical tests of this prediction are scarce. To fill this gap, we first used a set of 16 microsatellite markers to investigate the genetic structure of European beech in France in 65 populations from three refugial areas and one control recolonized (nonrefugial) area. Then, using the same approach, we reanalysed published isozyme data from 375 populations distributed across the entire species range. We found stronger genetic differentiation among populations in refugia than in recolonized areas. However, contrary to expectations, regional SGS was lower within refugia than within recolonized areas. Published studies presenting similar analyses suggest that our results could have generality across different biogeographical settings and types of organisms. Strong and prolonged genetic drift in refugial areas could have erased the signature of range expansions that is still visible in recolonized areas. Our results therefore suggest that Pleistocene population isolation has played a key role in increasing the genetic complexity of extant rear-edge populations.

Genotypic variability and phenotypic plasticity of cavitation resistance in Fagus sylvatica L. across Europe.

Xylem cavitation resistance is a key physiological trait correlated with species tolerance to extreme drought stresses. Little is known about the genetic variability and phenotypic plasticity of this trait in natural tree populations. Here we measured the cavitation resistance of 17 Fagus sylvatica populations representative of the full range of the species in Europe. The trees were grown in three field trials under contrasting climatic conditions. Our findings suggest that the genotypic variability of cavitation resistance is high between genotypes of a given population. By contrast, no significant differences were found for this trait across populations, the mean population cavitation resistance being remarkably constant in each trial. We found a significant site effect and a significant site × population interaction, suggesting that cavitation resistance has a high phenotypic plasticity and that this plasticity is under genetic control. The implications of our findings for beech forest management in a context of climate change are discussed.

MATING SYSTEM AND ASYMMETRIC HYBRIDIZATION IN A MIXED STAND OF EUROPEAN OAKS.

The sessile (Quercus petraea [Matt.] Liebl.) and pedunculate (Quercus robur L.) oaks are two closely related species having a wide sympatric distribution over Europe. Under natural conditions, they frequently form mixed forests, where hybridization is suspected to occur. In this paper, two different approaches have been applied to the study of the mating system and the interspecific gene flow in a mixed stand formed by the two species. The mating systems of both species have been studied separately by means of the mixed-mating model. The relative contribution of the parental species to the progenies have been estimated with two different methods. The first uses the admixture model. The second is an extension of the mixed-mating model and subdivides the outcrossing rate into intra- and interspecific components. The two species were almost completely outcrossing. This high level of outcrossing and interspecific gene flow could play an important role in the maintenance of the genetic diversity in these long-lived forest tree species. The contribution of the sessile oak to the pedunculate oak progenies varied from 17% to 48%. In contrast, ovules of sessile oak trees appear to be preferentially fertilized by other extreme sessile genotypes. We suggest that interspecific and directional gene flow was responsible for such patterns. Pedunculate oak is considered as a pioneer species and is progressively replaced by sessile oak. Our present findings add a further genetic component to this succession scheme, suggesting that unidirectional gene flow reinforces succession between the two species.