Adaptive potential of maritime pine under contrasting environments.

BACKGROUND: Predicting the adaptability of forest tree populations under future climates requires a better knowledge of both the adaptive significance and evolvability of measurable key traits. Phenotypic plasticity, standing genetic variation and degree of phenotypic integration shape the actual and future population genetic structure, but empirical estimations in forest tree species are still extremely scarce. We analysed 11 maritime pine populations covering the distribution range of the species (119 families and 8 trees/family, ca. 1300 trees) in a common garden experiment planted at two sites with contrasting productivity. We used plant height as a surrogate of fitness and measured five traits (mean and plasticity of carbon isotope discrimination, specific leaf area, needle biomass, Phenology growth index) related to four different strategies (acquisitive economics, photosynthetic organ size, growth allocation and avoidance of water stress). RESULTS: Estimated values of additive genetic variation would allow adaptation of the populations to future environmental conditions. Overall phenotypic integration and selection gradients were higher at the high productivity site, while phenotypic integration within populations was higher at the low productivity site. Response to selection was related mainly to photosynthetic organ size and drought-avoidance mechanisms rather than to water use efficiency. Phenotypic plasticity of water use efficiency could be maladaptive, resulting from selection for height growth. CONCLUSIONS: Contrary to the expectations in a drought tolerant species, our study suggests that variation in traits related to photosynthetic organ size and acquisitive investment of resources drive phenotypic selection across and within maritime pine populations. Both genetic variation and evolvability of key adaptive traits were considerably high, including plasticity of water use efficiency. These characteristics would enable a relatively fast micro-evolution of populations in response to the ongoing climate changes. Moreover, differentiation among populations in the studied traits would increase under the expected more productive future Atlantic conditions.

A claim for a 'next generation' of multisite range-wide forest genetic trials built on the legacy of ecological genetics to anticipate responses to climate.

Realizamos un comentario sobre el trabajo de revisión presentado por Leites & Benito-Garzón (2023) en este número sobre el valor de las redes de ensayos genéticos forestales multi-sitio para anticipar la distribución y las respuestas de las especies forestales a la actual crisis climática. La persistencia de las masas forestales y de los servicios ecosistémicos que proporcionan necesita de un esfuerzo coordinado sin precedentes de la comunidad científica para comprender, predecir y mitigar los impactos de la crisis climática sobre los bosques templados y boreales. Para ello se necesita información precisa sobre la variación genética, plasticidad fenotípica y la variación genética en plasticidad de las especies forestales. Hacemos una llamada de atención sobre la necesidad urgente de establecer una nueva generación de ensayos genéticos forestales de procedencias replicados multisitio a lo largo del rango de distribución de las especies de interés, y enfocado en la resiliencia al cambio climático. Como científicos debemos facilitar las respuestas evolutivas a los tan rápidos que están viviendo las especies forestales. Debemos ser determinantes en proporcionar información de primera clase a los actores políticos y administradores para desarrollar una gestión forestal basada en la evidencia científica.

Local seed sourcing for sustainable forestry.

Seed sourcing strategies are the basis for identifying genetic material meeting the requirements of future climatic conditions and social demands. Specifically, local seed sourcing has been extensively promoted, based on the expected adaptation of the populations to local conditions, but there are some limitations for the application. We analyzed Strict-sense local and Wide-sense local (based on climatic similarity) seed sourcing strategies. We determined species and genetic pools based on these strategies for 40 species and deployment zones in Spain. We also obtained the total number of seed sources and stands for these species in the EU countries. We analyzed the richness of the pools, the relationship with variables related to the use of the species in afforestation, and the availability of seed production areas approved for the production of reproductive material destined to be marketed. This study confirms the existence of extensive species and genetic local pools. Also, that the importance of these pools differs for different species, limitations being derived from the use of forest reproductive material and the existence of approved basic materials. Strategies derived from local seed sourcing approaches are the basis for the use of forest reproductive material because a large number of the species in the area considered in the study are under regulation. However, despite the extensive work done to approve basic materials, limitations based on the availability of seed production areas to provide local material for sustainable forestry are found in those species. Considering a Wide-sense local seed sourcing strategy we provide alternative pools in order to meet social demands under the actual regulations on marketing of reproductive materials.

Polygenic adaptation and negative selection across traits, years and environments in a long-lived plant species (Pinus pinaster Ait., Pinaceae).

A decade of genetic association studies in multiple organisms suggests that most complex traits are polygenic; that is, they have a genetic architecture determined by numerous loci, each with small effect-size. Thus, determining the degree of polygenicity and its variation across traits, environments and time is crucial to understand the genetic basis of phenotypic variation. We applied multilocus approaches to estimate the degree of polygenicity of fitness-related traits in a long-lived plant (Pinus pinaster Ait., maritime pine) and to analyse this variation across environments and years. We evaluated five categories of fitness-related traits (survival, height, phenology, functional, and biotic-stress response) in a clonal common-garden network planted in contrasted environments (over 20,500 trees). Most of the analysed traits showed evidence of local adaptation based on Qst -Fst comparisons. We further observed a remarkably stable degree of polygenicity, averaging 6% (range of 0%-27%), across traits, environments and years. We detected evidence of negative selection, which could explain, at least partially, the high degree of polygenicity. Because polygenic adaptation can occur rapidly, our results suggest that current predictions on the capacity of natural forest tree populations to adapt to new environments should be revised, especially in the current context of climate change.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine.

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

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.

The legacy of climate variability over the last century on populations' phenotypic variation in tree height.

Phenotypic plasticity and local adaptation are the two main processes underlying trait variability. Under rapid environmental change, phenotypic plasticity, if adaptive, could increase the odds for organisms to persist. However, little is known on how environmental variation has shaped plasticity across species ranges over time. Here, we assess whether the portion of phenotypic variation of tree populations linked to the environment is related to the inter-annual climate variability of the last century and how it varies among populations across species ranges and age. To this aim, we used 372,647 individual tree height measurements of three pine species found in low elevation forests in Europe: Pinus nigra Arnold, P. pinaster Aiton and P. pinea L. Measurements were taken in a network of 38 common gardens established in Europe and North Africa with 315 populations covering the distribution range of the species. We fitted linear mixed-effect models of tree height as a function of age, population, climate and competition effects. Models allowed us to estimate tree height response curves at the population level and indexes of populations' phenotypic variation, as a proxy of phenotypic plasticity, at 4, 8 and 16 years old, and relate these indexes to the inter-annual climate variability of the last century. We found that phenotypic variation in tree height was higher in young trees than in older ones. We also found that P. pinea showed the highest phenotypic variation in tree height compared with P. pinaster and P. nigra. Finally, phenotypic variation in tree height may be partly adaptive, and differently across species, as climate variability during the last century at the origin of the populations explained between 51 and 69% of the current phenotypic variation of P. nigra and P. pinea, almost twice of the levels of P. pinaster. MAIN CONCLUSIONS: Populations' phenotypic variation in tree height is largely explained by the climate variability that the populations experienced during the last century, which we attribute to the genetic diversity among populations.

Author Correction: The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe.

The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios.

Incorporating exposure to pitch canker disease to support management decisions of Pinus pinaster Ait. in the face of climate change.

Climate change is gravely affecting forest ecosystems, resulting in large distribution shifts as well as in increasing infection diseases and biological invasions. Accordingly, forest management requires an evaluation of exposure to climate change that should integrate both its abiotic and biotic components. Here we address the implications of climate change in an emerging disease by analysing both the host species (Pinus pinaster, Maritime pine) and the pathogen's (Fusarium circinatum, pitch canker) environmental suitability i.e. estimating the host's risk of habitat loss and the disease`s future environmental range. We constrained our study area to the Spanish Iberian Peninsula, where accurate climate and pitch canker occurrence databases were available. While P. pinaster is widely distributed across the study area, the disease has only been detected in its north-central and north-western edges. We fitted species distribution models for the current distribution of the conifer and the disease. Then, these models were projected into nine Global Climate Models and two different climatic scenarios which totalled to 18 different future climate predictions representative of 2050. Based on the level of agreement among them, we created future suitability maps for the pine and for the disease independently, which were then used to assess exposure of current populations of P. pinaster to abiotic and biotic effects of climate change. Almost the entire distribution of P. pinaster in the Spanish Iberian Peninsula will be subjected to abiotic exposure likely to be driven by the predicted increase in drought events in the future. Furthermore, we detected a reduction in exposure to pitch canker that will be concentrated along the north-western edge of the study area. Setting up breeding programs is recommended in highly exposed and productive populations, while silvicultural methods and monitoring should be applied in those less productive, but still exposed, populations.

Life-history correlations with seasonal cold hardiness in maritime pine.

PREMISE OF THE STUDY: Plants have developed mechanisms to withstand stressful environmental conditions, but the high energetic cost of these mechanisms may involve exchanges with other key functions. While trade-offs between cold hardiness and growth rates are a general assumption, we lack information regarding genetically based trade-offs between cold hardiness and other life-history traits. Such information has strong implications for tree conservation and breeding, especially in the context of ongoing climate change. METHODS: We used a common garden progeny test to examine the relationships between seasonal cold hardiness and life-history traits of growth, reproduction, juvenile ontogeny, and phenology in 75 families of six maritime pine (Pinus pinaster Ait.) populations, three of continental and three of coastal origins. KEY RESULTS: We found a clear differentiation among populations with regard to cold hardiness and life-history traits. Two continental Iberian populations showed high cold tolerance and slower growth, but faster ontogenetic development in relation to both vegetative heteroblastic change in juveniles and the onset of female reproduction. The coastal populations displayed the opposite behavior, while the continental Moroccan population presented a unique combination of traits. We confirmed trade-offs between cold-hardiness and growth at the population level, but not within populations. There were no trade-offs with other life-history traits at either level. CONCLUSIONS: Relevant local adaptation syndromes were identified in the relationship between cold hardiness and life-history traits. These should be considered in developing tree management guidelines aimed at increasing productivity or adaptability under the expected conditions of climate change.

Capturing neutral and adaptive genetic diversity for conservation in a highly structured tree species.

Preserving intraspecific genetic diversity is essential for long-term forest sustainability in a climate change scenario. Despite that, genetic information is largely neglected in conservation planning, and how conservation units should be defined is still heatedly debated. Here, we use maritime pine (Pinus pinaster Ait.), an outcrossing long-lived tree with a highly fragmented distribution in the Mediterranean biodiversity hotspot, to prove the importance of accounting for genetic variation, of both neutral molecular markers and quantitative traits, to define useful conservation units. Six gene pools associated to distinct evolutionary histories were identified within the species using 12 microsatellites and 266 single nucleotide polymorphisms (SNPs). In addition, height and survival standing variation, their genetic control, and plasticity were assessed in a multisite clonal common garden experiment (16 544 trees). We found high levels of quantitative genetic differentiation within previously defined neutral gene pools. Subsequent cluster analysis and post hoc trait distribution comparisons allowed us to define 10 genetically homogeneous population groups with high evolutionary potential. They constitute the minimum number of units to be represented in a maritime pine dynamic conservation program. Our results uphold that the identification of conservation units below the species level should account for key neutral and adaptive components of genetic diversity, especially in species with strong population structure and complex evolutionary histories. The environmental zonation approach currently used by the pan-European genetic conservation strategy for forest trees would be largely improved by gradually integrating molecular and quantitative trait information, as data become available.

Adaptation and plasticity in aboveground allometry variation of four pine species along environmental gradients.

Plant species aboveground allometry can be viewed as a functional trait that reflects the evolutionary trade-off between above- and belowground resources. In forest trees, allometry is related to productivity and resilience in different environments, and it is tightly connected with a compromise between efficiency-safety and competitive ability. A better understanding on how this trait varies within and across species is critical to determine the potential of a species/population to perform along environmental gradients. We followed a hierarchical framework to assess tree height-diameter allometry variation within and across four common European Pinus species. Tree height-diameter allometry variation was a function of solely genetic components -approximated by either population effects or clinal geographic responses of the population's site of origin- and differential genetic plastic responses -approximated by the interaction between populations and two climatic variables of the growing sites (temperature and precipitation)-. Our results suggest that, at the species level, climate of the growing sites set the tree height-diameter allometry of xeric and mesic species (Pinus halepensis, P. pinaster and P. nigra) apart from the boreal species (P. sylvestris), suggesting a weak signal of their phylogenies in the tree height-diameter allometry variation. Moreover, accounting for interpopulation variability within species for the four pine species aided to: (1) detect genetic differences among populations in allometry variation, which in P. nigra and P. pinaster were linked to gene pools -genetic diversity measurements-; (2) reveal the presence of differential genetic variation in plastic responses along two climatic gradients in tree allometry variation. In P. sylvestris and P. nigra, genetic variation was the result of adaptive patterns to climate, while in P. pinaster and P. halepensis, this signal was either weaker or absent, respectively; and (3) detect local adaptation in the exponent of the tree height-diameter allometry relationship in two of the four species (P. sylvestris and P. nigra), as it was a function of populations' latitude and altitude variables. Our findings suggest that the four species have been subjected to different historical and climatic constraints that might have driven their aboveground allometry and promoted different life strategies.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree.

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Nucleotide polymorphisms in a pine ortholog of the Arabidopsis degrading enzyme cellulase KORRIGAN are associated with early growth performance in Pinus pinaster.

We have carried out a candidate-gene-based association genetic study in Pinus pinaster Aiton and evaluated the predictive performance for genetic merit gain of the most significantly associated genes and single nucleotide polymorphisms (SNPs). We used a second generation 384-SNP array enriched with candidate genes for growth and wood properties to genotype mother trees collected in 20 natural populations covering most of the European distribution of the species. Phenotypic data for total height, polycyclism, root-collar diameter and biomass were obtained from a replicated provenance-progeny trial located in two sites with contrasting environments (Atlantic vs Mediterranean climate). General linear models identified strong associations between growth traits (total height and polycyclism) and four SNPs from the korrigan candidate gene, after multiple testing corrections using false discovery rate. The combined genomic breeding value predictions assessed for the four associated korrigan SNPs by ridge regression-best linear unbiased prediction (RR-BLUP) and cross-validation accounted for up to 8 and 15% of the phenotypic variance for height and polycyclic growth, respectively, and did not improve adding SNPs from other growth-related candidate genes. For root-collar diameter and total biomass, they accounted for 1.6 and 1.1% of the phenotypic variance, respectively, but increased to 15 and 4.1% when other SNPs from lp3.1, lp3.3 and cad were included in RR-BLUP models. These results point towards a desirable integration of candidate-gene studies as a means to pre-select relevant markers, and aid genomic selection in maritime pine breeding programs.

Correlated genetic effects on reproduction define a domestication syndrome in a forest tree.

Compared to natural selection, domestication implies a dramatic change in traits linked to fitness. A number of traits conferring fitness in the wild might be detrimental under domestication, and domesticated species typically differ from their ancestors in a set of traits known as the domestication syndrome. Specifically, trade-offs between growth and reproduction are well established across the tree of life. According to allocation theory, selection for growth rate is expected to indirectly alter life-history reproductive traits, diverting resources from reproduction to growth. Here we tested this hypothesis by examining the genetic change and correlated responses of reproductive traits as a result of selection for timber yield in the tree Pinus pinaster. Phenotypic selection was carried out in a natural population, and progenies from selected trees were compared with those of control trees in a common garden experiment. According to expectations, we detected a genetic change in important life-history traits due to selection. Specifically, threshold sizes for reproduction were much higher and reproductive investment relative to size significantly lower in the selected progenies just after a single artificial selection event. Our study helps to define the domestication syndrome in exploited forest trees and shows that changes affecting developmental pathways are relevant in domestication processes of long-lived plants.

Molecular proxies for climate maladaptation in a long-lived tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Adaptive potential of maritime pine (Pinus pinaster) populations to the emerging pitch canker pathogen, Fusarium circinatum.

There is a concern on how emerging pests and diseases will affect the distribution range and adaptability of their host species, especially due to different conditions derived from climate change and growing globalization. Fusarium circinatum, which causes pitch canker disease in Pinus species, is an exotic pathogen of recent introduction in Spain that threatens its maritime pine (P. pinaster) stands. To predict the impact this disease will have on the species, we examine host resistance traits and their genetic architecture. Resistance phenotyping was done in a clonal provenance/progeny trial, using three-year-old cuttings artificially inoculated with the pathogen and maintained under controlled environmental conditions. A total number of 670 ramets were assessed, distributed in 10 populations, with a total of 47 families, 2 to 5 half-sibs per family, and 3-7 ramets per clone. High genetic variation was found at the three hierarchical levels studied: population, family and clone, being both additive and non-additive effects important. Narrow-sense and broad-sense heritability estimates were relatively high, with respective values of 0.43-0.58 and 0.51-0.8, depending on the resistance traits measured (lesion length, lesion length rate, time to wilting, and survival). These values suggest the species' high capacity of evolutionary response to the F. circinatum pathogen. A population originated in Northern Spain was the most resistant, while another from Morocco was the most susceptible. The total number of plants that did not show lesion development or presented a small lesion (length<30 mm) was 224 out of 670, indicating a high proportion of resistant trees in the offspring within the analyzed populations. We found large differences among populations and considerable genetic variation within populations, which should allow, through natural or artificial selection, the successful adaptation of maritime pine to pitch canker disease.

The role of population origin and microenvironment in seedling emergence and early survival in Mediterranean maritime pine (Pinus pinaster Aiton).

Understanding tree recruitment is needed to forecast future forest distribution. Many studies have reported the relevant ecological factors that affect recruitment success in trees, but the potential for genetic-based differences in recruitment has often been neglected. In this study, we established a semi-natural reciprocal sowing experiment to test for local adaptation and microenvironment effects (evaluated here by canopy cover) in the emergence and early survival of maritime pine (Pinus pinaster Aiton), an emblematic Mediterranean forest tree. A novel application of molecular markers was also developed to test for family selection and, thus, for potential genetic change over generations. Overall, we did not find evidence to support local adaptation at the recruitment stage in our semi-natural experiment. Moreover, only weak family selection (if any) was found, suggesting that in stressful environments with low survival, stochastic processes and among-year climate variability may drive recruitment. Nevertheless, our study revealed that, at early stages of recruitment, microenvironments may favor the population with the best adapted life strategy, irrespectively of its (local or non-local) origin. We also found that emergence time is a key factor for seedling survival in stressful Mediterranean environments. Our study highlights the complexity of the factors influencing the early stages of establishment of maritime pine and provides insights into possible management actions aimed at environmental change impact mitigation. In particular, we found that the high stochasticity of the recruitment process in stressful environments and the differences in population-specific adaptive strategies may difficult assisted migration schemes.