Targeted bisulfite sequencing of Scots pine adaptation-related genes.

During environmental changes, epigenetic processes can enable adaptive responses faster than natural selection. In plants, very little is known about the role of DNA methylation during long-term adaptation. Scots pine is a widely distributed coniferous species which must adapt to different environmental conditions throughout its long lifespan. Thus, epigenetic modifications may contribute towards this direction. We provide bisulfite next-generation sequencing data from the putative promoters and exons of eight adaptation-related genes (A3IP2, CCA1, COL1, COL2, FTL2, MFT1, PHYO, and ZTL) in three Scots pine populations located in northern and southern parts of Finland. DNA methylation levels were studied in the two seed tissues: the maternal megagametophyte which contributes to embryo viability, and the biparental embryo which represents the next generation. In most genes, differentially methylated cytosines (DMCs) were in line with our previously demonstrated gene expression differences found in the same Scots pine populations. In addition, we found a strong correlation of total methylation levels between the embryo and megagametophyte tissues of a given individual tree, which indicates that DNA methylation can be inherited from the maternal parent. In conclusion, our results imply that DNA methylation differences may contribute to the adaptation of Scots pine populations in different climatic conditions.

Transcriptome profiles reveal response mechanisms and key role of PsNAC1 in Pinus sylvestris var. mongolica to drought stress.

BACKGROUND: Drought stress severely impedes plant growth, and only a limited number of species exhibit long-term resistance to such conditions. Pinus sylvestris var. mongolica, a dominant tree species in arid and semi-arid regions of China, exhibits strong drought resistance and plays a crucial role in the local ecosystem. However, the molecular mechanisms underlying this resistance remain poorly understood. RESULTS: Here, we conducted transcriptome sequence and physiological indicators analysis of needle samples during drought treatment and rehydration stages. De-novo assembly yielded approximately 114,152 unigenes with an N50 length of 1,363 bp. We identified 6,506 differentially expressed genes (DEGs), with the majority being concentrated in the heavy drought stage (4,529 DEGs). Functional annotation revealed enrichment of drought-related GO terms such as response to water (GO:0009415: enriched 108 genes) and response to water deprivation (GO:0009414: enriched 106 genes), as well as KEGG categories including MAPK signaling pathway (K04733: enriched 35 genes) and monoterpenoid biosynthesis (K21374: enriched 27 genes). Multiple transcription factor families and functional protein families were differentially expressed during drought treatment. Co-expression network analysis identified a potential drought regulatory network between cytochrome P450 genes (Unigene4122_c1_g1) and a core regulatory transcription factor Unigene9098_c3_g1 (PsNAC1) with highly significant expression differences. We validated PsNAC1 overexpression in Arabidopsis and demonstrated enhanced drought resistance. CONCLUSIONS: These findings provide insight into the molecular basis of drought resistance in P. sylvestris var. mongolica and lay the foundation for further exploration of its regulatory network.

Stronger genetic differentiation among within-population genetic groups than among populations in Scots pine provides new insights into within-population genetic structuring.

We investigated the presence of spatial genetic groups within forest tree populations and determined if the genetic divergence among these groups is greater than that between populations using Scots pine (Pinus sylvestris) as a model species. We genotyped 890 adult trees of Scots pine in six natural populations in Lithuania at 11 nuclear microsatellite loci. We used a Bayesian clustering approach to identify the within-population genetic groups within each of the six populations. We calculated the differentiation indexes among the genetic groups within each population and among the six populations by ignoring the genetic groups. The Bayesian clustering revealed 2 to 6 distinct genetic groups of varying size as the most likely genetic structures within populations. The genetic differentiation indexes among the genetic groups within populations were nearly tenfold greater (FST = 0.012-0.070) than those between the populations (FST = 0.003). We conclude on the existence of markedly stronger structuring of genetic variation within populations than between populations of Scots pine in large forest tracts of northern Europe. Such genetic structures serve as a contributing factor to large within population genetic diversity in northern conifers. We assume that within population mating in Scots pine is not completely random but rather is stratified into genetic clusters. Our study provides pioneering novel key insights into structuring of genetic variation within populations. Our findings have implications for examining within-population genetic diversity and genetic structure, conservation, and management of genetic resources.

Scots pine - panmixia and the elusive signal of genetic adaptation.

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

Microgeographic variation in early fitness traits of Pinus sylvestris from contrasting soils.

PREMISE: The possibility of fine-scale intraspecific adaptive divergence under gene flow is established by theoretical models and has been confirmed empirically in tree populations distributed along steep altitudinal clines or across extreme edaphic discontinuities. However, the possibility of microgeographic adaptive divergence due to less severe but more frequent kinds of soil variation is unclear. METHODS: In this study, we looked for evidence of local adaptation to calcareous versus siliceous soil types in two nearby Mediterranean Pinus sylvestris populations connected via pollen flow. Using a greenhouse experiment, we tested for variation in early (up to three years of age) seedling performance among open-pollinated maternal families originating from each edaphic provenance when experimentally grown on both types of natural local substrate. RESULTS: Although seedlings were clearly affected by the edaphic environment, exhibiting lower and slower emergence as well as higher mortality on the calcareous than in the siliceous substrate, neither the performance on each substrate nor the plasticity among substrates varied significantly with seedling edaphic provenance. CONCLUSIONS: We found no evidence of local adaptation to a non-extreme edaphic discontinuity over a small spatial scale, at least during early stages of seedling establishment. Future studies on microgeographic soil-driven adaptation should consider long-term experiments to minimize maternal effects and allow a potentially delayed expression of edaphic adaptive divergence.

Manganese Deficiency Suppresses Growth and Photosynthetic Processes but Causes an Increase in the Expression of Photosynthetic Genes in Scots Pine Seedlings.

Manganese deficiency is a serious plant nutritional disorder, resulting in the loss of crop productivity in many parts of the world. Despite the progress made in the study of angiosperms, the demand for Mn in gymnosperms and the physiological responses to Mn deficiency remain unexplored. We studied the influence of Mn deficiency for 24 weeks on Pinus sylvestris L. seedling growth, ion homeostasis, pigment contents, lipid peroxidation, chlorophyll fluorescence indices and the transcript levels of photosynthetic genes and genes involved in chlorophyll biosynthesis. It was shown that Mn-deficient plants demonstrated suppressed growth when the Mn content in the needles decreased below 0.34 µmol/g DW. The contents of photosynthetic pigments decreased when the Mn content in the needles reached 0.10 µmol/g DW. Mn deficiency per se did not lead to a decrease in the nutrient content in the organs of seedlings. Photoinhibition of PSII was observed in Mn-deficient plants, although this was not accompanied by the development of oxidative stress. Mn-deficient plants had an increased transcript abundance of genes (psbO, psbP, psbQ, psbA and psbC), encoding proteins directly associated with the Mn cluster also as other proteins involved in photosynthesis, whose activities do not depend on Mn directly. Furthermore, the transcript levels of the genes encoding the large subunit of Rubisco, light-dependent NADPH-protochlorophyllide oxidoreductase and subunits of light-independent protochlorophyllide reductase were also increased in Mn-deficient plants.

Phenotypic trait variation in a long-term multisite common garden experiment of Scots pine in Scotland.

Multisite common garden experiments, exposing common pools of genetic diversity to a range of environments, allow quantification of plastic and genetic components of trait variation. For tree species, such studies must be long term as they typically only express mature traits after many years. As well as evaluating standing genetic diversity, these experiments provide an ongoing test of genetic variation against changing environmental conditions and form a vital resource for understanding how species respond to abiotic and biotic variation. Finally, quantitative assessments of phenotypic variation are essential to pair with rapidly accumulating genomic data to advance understanding of the genetic basis of trait variation, and its interaction with climatic change. We describe a multisite, population-progeny, common garden experiment of the economically and ecologically important tree species, Scots pine, collected from across its native range in Scotland and grown in three contrasting environments. Phenotypic traits, including height, stem diameter and budburst were measured over 14 growing seasons from nursery to field site. The datasets presented have a wide range of applications.

Adaptive strategies of Scots pine under shade: Increase in lignin synthesis and ecotypic variation in defense-related gene expression.

Shade is a stressful condition for plants characterized by low Red:Far-Red (R:FR) ratio. The northern latitudes in Sweden daily receive more hours of FR-enriched light (twilight) or shade-like conditions compared to southern forests during the growing season. Scots pine (Pinus sylvestris L.) is a shade-intolerant species. Yet, it is well adapted to this latitudinal variation in light, which is evident by a northward increase in FR requirement to maintain growth. Shade adversely affects plant growth; it makes the plant weak and, therefore, susceptible to pathogen attack. Lignin is involved in plant protection against pathogen invasion mainly by forming a physical barrier. We studied lignin synthesis and expression of defense-related genes (growth-defense trade-offs) under a low R:FR (shade) ratio in Scots pine. A higher number of immunity/defense-related genes were up-regulated in response to shade in northern populations compared to southern ones, which can be viewed as a local adaptation to light quality for optimal growth and survival. Light quality regulates lignin metabolism; light stimulates lignin synthesis, while shade causes a decrease in lignin synthesis in most angiosperms. In contrast, Scots pine shows an increase in lignin synthesis supported by the higher expression of a few key genes in the lignin biosynthetic pathway, a novel finding reported by our study. These findings can be applied to future breeding strategies in forestry to produce disease-resilient trees.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Influence of Light of Different Spectral Compositions on the Growth, Photosynthesis, and Expression of Light-Dependent Genes of Scots Pine Seedlings.

Varying the spectral composition of light is one of the ways to accelerate the growth of conifers under artificial conditions for the development of technologies and to obtain sustainable seedlings required to preserve the existing areas of forests. We studied the influence of light of different quality on the growth, gas exchange, fluorescence indices of Chl a, and expression of key light-dependent genes of Pinus sylvestris L. seedlings. It was shown that in plants growing under red light (RL), the biomass of needles and root system increased by more than two and three times, respectively, compared with those of the white fluorescent light (WFL) control. At the same time, the rates of photosynthesis and respiration in RL and blue light (BL) plants were lower than those of blue red light (BRL) plants, and the difference between the rates of photosynthesis and respiration, which characterizes the carbon balance, was maximum under RL. RL influenced the number of xylem cells, activated the expression of genes involved in the transduction of cytokinin (Histidine-containing phosphotransfer 1, HPT1, Type-A Response Regulators, RR-A) and auxin (Auxin-induced protein 1, Aux/IAA) signals, and reduced the expression of the gene encoding the transcription factor phytochrome-interacting factor 3 (PIF3). It was suggested that RL-induced activation of key genes of cytokinin and auxin signaling might indicate a phytochrome-dependent change in cytokinins and auxins activity.

The evolutionary heritage and ecological uniqueness of Scots pine in the Caucasus ecoregion is at risk of climate changes.

Scots pine is one of the most widely occurring pines, but future projections suggest a large reduction in its range, mostly at the southern European limits. A significant part of its range is located in the Caucasus, a global hot-spot of diversity. Pine forests are an important reservoir of biodiversity and endemism in this region. We explored demographic and biogeographical processes that shaped the genetic diversity of Scots pine in the Caucasus ecoregion and its probable future distribution under different climate scenarios. We found that the high genetic variability of the Caucasian populations mirrors a complex glacial and postglacial history that had a unique evolutionary trajectory compared to the main range in Europe. Scots pine currently grows under a broad spectrum of climatic conditions in the Caucasus, which implies high adaptive potential in the past. However, the current genetic resources of Scots pine are under high pressure from climate change. From our predictions, over 90% of the current distribution of Scots pine may be lost in this century. By threatening the stability of the forest ecosystems, this would dramatically affect the biodiversity of the Caucasus hot-spot.

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.

Polyamine Metabolism in Scots Pine Embryogenic Cells under Potassium Deficiency.

Polyamines (PA) have a protective role in maintaining growth and development in Scots pine during abiotic stresses. In the present study, a controlled liquid Scots pine embryogenic cell culture was used for studying the responses of PA metabolism related to potassium deficiency. The transcription level regulation of PA metabolism led to the accumulation of putrescine (Put). Arginine decarboxylase (ADC) had an increased expression trend under potassium deficiency, whereas spermidine synthase (SPDS) expression decreased. Generally, free spermidine (Spd) and spermine (Spm)/ thermospermine (t-Spm) contents were kept relatively stable, mostly by the downregulation of polyamine oxidase (PAO) expression. The low potassium contents in the culture medium decreased the potassium content of the cells, which inhibited cell mass growth, but did not affect cell viability. The reduced growth was probably caused by repressed metabolic activity and cell division, whereas there were no signs of H2O2-induced oxidative stress or increased cell death. The low intracellular content of K+ decreased the content of Na+. The decrease in the pH of the culture medium indicated that H+ ions were pumped out of the cells. Altogether, our findings emphasize the specific role(s) of Put under potassium deficiency and strict developmental regulation of PA metabolism in Scots pine.

Quantitative analysis of differential dehydrin regulation in pine and spruce seedlings under water deficit.

Dehydrins are well-known components of plant responses to different stresses that cause dehydration, including drought, freezing, salinity, etc. In conifers, the dehydrin gene family is very large, implying that the members of this family have important physiological functions in conifer stress tolerance. However, dehydrin gene expression displays a wide range of responses to stress, from thousand-fold increased expression to decreased expression, and it is generally unknown how regulatory systems are connected at the mRNA and protein levels. Therefore, we studied these aspects of dehydrin regulation in Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) H. Karst) seedlings under polyethylene glycol 6000-induced osmotic stress ranging from relatively low (culture medium water potential of -0.15 MPa) to very high (-1.0 MPa) intensities. In pine, the major dehydrin protein was Dhn1 in both the roots and needles, and in spruce, two isoforms of the Dhn4 protein were the major dehydrins; both of these proteins are AESK-type dehydrins. The genes encoding these major proteins were highly expressed even under control conditions; surprisingly, we also observed several highly expressed dehydrin genes that were not abundantly translated. Under osmotic stress, the most prominent expression changes were observed for the dehydrin genes with low basal expression levels, whereas highly expressed genes generally demonstrated rather modest changes in expression. We report proposed constitutive physiological functions of the AESK-type dehydrins in Pinaceae plants.

Transcript Dynamics in Wounded and Inoculated Scots Pine.

Comparative transcriptome analysis provides a useful tool for the exploration of plant-pathogen interaction by allowing in-depth comparison of gene expression between unaffected, inoculated and wounded organisms. Here we present the results of comparative transcriptome analysis in genetically identical one-year-old Scots pine ramets after wounding and inoculation with Heterobasidion annosum. We identified 230 genes that were more than 2-fold upregulated in inoculated samples (compared to controls) and 116 downregulated genes. Comparison of inoculated samp les with wounded samples identified 32 differentially expressed genes (30 were upregulated after inoculation). Several of the genes upregulated after inoculation are involved in protection from oxidative stress, while genes involved in photosynthesis, water transport and drought stress tolerance were downregulated. An NRT3 family protein was the most upregulated transcript in response to both inoculation and wounding, while a U-box domain-containing protein gene was the most upregulated gene comparing inoculation to wounding. The observed transcriptome dynamics suggest involvement of auxin, ethylene, jasmonate, gibberellin and reactive oxygen species pathways and cell wall modification regulation in response to H. annosum infection. The results are compared to methyl jasmonate induced transcriptome dynamics.

Imprints of selection in peripheral and ecologically marginal central-eastern European Scots pine populations.

Knowledge of the molecular mechanisms underlying the stress response in plants is essential to understand evolutionary processes that result in long-term persistence of populations. Populations inhabiting marginal ecological conditions at the distribution range periphery may have preserved imprints of natural selection that have shaped functional genetic variation of the species. Our aim was to evaluate the extent of selection processes in the extremely fragmented, peripheral and isolated populations of Scots pine in central-eastern Europe. Autochthonous populations of the Carpathian Mts. and the Pannonian Basin were sampled and drought stress-related candidate genes were re-sequenced. Neutrality tests and outlier detection approaches were applied to infer the effect and direction of selection. Populations retained high genetic diversity by preserving a high number of alleles and haplotypes, many of them being population specific. Neutrality tests and outlier detection highlighted nucleotide positions that are under divergent selection and may be involved in local adaptation. The detected genetic pattern confirms that natural selection has played an important role in shaping modern-day genetic variation in marginal Scots pine populations, allowing for the long-term persistence of populations. Selection detected at functional regions possibly acts to maintain diversity and counteract the effect of genetic erosion.

Divergent patterns between phenotypic and genetic variation in Scots pine.

In boreal forests, autumn frost tolerance in seedlings is a critical fitness component because it determines survival rates during regeneration. To understand the forces that drive local adaptation in this trait, we conducted freezing tests in a common garden setting for 54 Pinus sylvestris (Scots pine) populations (>5000 seedlings) collected across Scandinavia into western Russia, and genotyped 24 of these populations (>900 seedlings) at >10 000 SNPs. Variation in cold hardiness among populations, as measured by QST , was above 80% and followed a distinct cline along latitude and longitude, demonstrating significant adaptation to climate at origin. In contrast, the genetic differentiation was very weak (mean FST 0.37%). Despite even allele frequency distribution in the vast majority of SNPs among all populations, a few rare alleles appeared at very high or at fixation in marginal populations restricted to northwestern Fennoscandia. Genotype-environment associations showed that climate variables explained 2.9% of the genetic differentiation, while genotype-phenotype associations revealed a high marker-estimated heritability of frost hardiness of 0.56, but identified no major loci. Very extensive gene flow, strong local adaptation, and signals of complex demographic history across markers are interesting topics of forthcoming studies on this species to better clarify signatures of selection and demography.

Phosphate availability and ectomycorrhizal symbiosis with Pinus sylvestris have independent effects on the Paxillus involutus transcriptome.

Many plant species form symbioses with ectomycorrhizal fungi, which help them forage for limiting nutrients in the soil such as inorganic phosphate (Pi). The transcriptional responses to symbiosis and nutrient-limiting conditions in ectomycorrhizal fungal hyphae, however, are largely unknown. An artificial system was developed to study ectomycorrhizal basidiomycete Paxillus involutus growth in symbiosis with its host tree Pinus sylvestris at different Pi concentrations. RNA-seq analysis was performed on P. involutus hyphae growing under Pi-limiting conditions, either in symbiosis or alone. We show that Pi starvation and ectomycorrhizal symbiosis have an independent effect on the P. involutus transcriptome. Notably, low Pi availability induces expression of newly identified putative high-affinity Pi transporter genes, while reducing the expression of putative organic acid transporters. Additionally, low Pi availability induces a close transcriptional interplay between P and N metabolism. GTP-related signalling was found to have a positive effect in the maintenance of ectomycorrhizal symbiosis, whereas multiple putative cytochrome P450 genes were found to be downregulated, unlike arbuscular mycorrhizal fungi. We provide the first evidence of global transcriptional changes induced by low Pi availability and ectomycorrhizal symbiosis in the hyphae of P. involutus, revealing both similarities and differences with better-characterized arbuscular mycorrhizal fungi.

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.