# Environment-dependent microevolution in a Mediterranean pine (Pinus pinaster Aiton)

Ricardo Alia, Regina Chambel, Eduardo Notivol, José Climent, Santiago C.

González-Martínez

► To cite this version:

Ricardo Alia, Regina Chambel, Eduardo Notivol, José Climent, Santiago C. González-Martínez. Environment-dependent microevolution in a Mediterranean pine (Pinus pinaster Aiton). BMC Evolutionary Biology, 2014, 14 (1), pp.1-12. 10.1186/s12862-014-0200-5. http://hal-02639262

HAL Id: hal-02639262

https://hal.inrae.fr/hal-02639262

Submitted on 28 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Environment-dependent microevolution in a Mediterranean pine (*Pinus pinaster* Aiton)

Ricardo Alía<fnref n="1" />,<fnref n="2" />, Regina Chambel<fnref n="1" />, Eduardo Notivol³, José Climent<fnref n="1" />,<fnref n="2" /> and Santiago C González-Martínez<fnref n="1" />,<fnref n="2" />,⁴*

## Abstract

**Background:** A central question for understanding the evolutionary responses of plant species to rapidly changing environments is the assessment of their potential for short-term (in one or a few generations) genetic change. In our study, we consider the case of *Pinus pinaster* Aiton (maritime pine), a widespread Mediterranean tree, and (i) test, under different experimental conditions (growth chamber and semi-natural), whether higher recruitment in the wild from the most successful mothers is due to better performance of their offspring; and (ii) evaluate genetic change in quantitative traits across generations at two different life stages (mature trees and seedlings) that are known to be under strong selection pressure in forest trees.

**Results:** Genetic control was high for most traits ($h^2$ = 0.137-0.876) under the milder conditions of the growth chamber, but only for ontogenetic change (0.276), total height (0.415) and survival (0.719) under the more stressful semi-natural conditions. Significant phenotypic selection gradients were found in mature trees for traits related to seed quality (germination rate and number of empty seeds). Moreover, female relative reproductive success was significantly correlated with offspring performance for specific leaf area (SLA) in the growth chamber experiment, and stem mass fraction (SMF) in the experiment under semi-natural conditions, two adaptive traits related to abiotic stress-response in pines. Selection gradients based on genetic covariance of seedling traits and responses to selection at this stage involved traits related to biomass allocation (SMF) and growth (as decomposed by a Gompertz model) or delayed ontogenetic change, depending also on the testing environment.

**Conclusions:** Despite the evidence of microevolutionary change in adaptive traits in maritime pine, directional or disruptive changes are difficult to predict due to variable selection at different life stages and environments. At mature-tree stages, higher female effective reproductive success can be explained by differences in their production of offspring (due to seed quality) and, to a lesser extent, by seemingly better adapted seedlings. Selection gradients and responses to selection for seedlings also differed across experimental conditions. The distinct processes involved at the two life stages (mature trees or seedlings) together with environment-specific responses advice caution when predicting likely evolutionary responses to environmental change in Mediterranean forest trees.

**Keywords:** Adaptive traits, Response to selection, Selection gradients, Phenotypic directional selection, Mediterranean forest

## Background

To confront the rapid environmental change caused by global warming, plants will have either to migrate or to adapt in-situ to new conditions [1]. Mounting evidence on seed dispersal and plant migration rates has revealed that plant migration capability falls short, in various orders

of magnitude, of that needed to track ecological optima in the near future [2], thus highlighting the importance of within-population standing genetic variation, phenotypic plasticity, and adaptive evolution for plant population survival. In the case of foundation species such as forest trees (see [3]), their adaptive potential is of paramount importance, as it allows to predict the evolutionary responses of entire terrestrial ecosystems in the context of global climate change [4]. A central question in the adaptability of forest tree species to changing environments is their potential for short-term genetic change resulting from differences

<footnote>* Correspondence: santiago@inia.es</footnote>

<footnote n="1">¹Department of Forest Ecology and Genetics, INIA-Forest Research Centre (CIFOR), Avda. A Coruña km 7.5, 28040, Madrid, Spain</footnote>

<footnote n="2">²Sustainable Forest Management Research Institute, University of Valladolid-INIA, Madrid, Spain</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
in parental reproductive success and selection at early re-
cruitment stages, when mortality rates are higher.

Adaptive evolutionary change at any given trait requires
an association between the trait and fitness, plus the herit-
able transmission of the trait across generations [5]. The first
component can be estimated by phenotypic and genetic se-
lection gradients, and the second by measuring the response
to selection. Phenotypic selection gradients (defined as the
vector of partial regression coefficients of individual relative
fitness on the traits [6]) have been characterized for different
plant species at population level, especially through parent-
age analysis based on molecular markers (e.g. [7,8]).
They pose, however, some limitations derived from (i)
potential trait multicollinearity, (ii) biases introduced
by failure to include traits that covary with fitness, (iii)
deviations of response variables from multivariate nor-
mality assumed for hypothesis testing, and (iv) failure to
consider environmental factors that may induce spurious
correlations [9]. Furthermore, it is not clear what the im-
plications are of such selection gradients for a population's
adaptive potential, i.e. its ability to respond to selective
pressures [10]. Some of the limitations mentioned can
be overcome by using genetic selection gradients (i.e. the
strength of selection acting on breeding values), as breed-
ing values "are correlated with their corresponding pheno-
typic value but are, by definition, uncorrelated with error
deviations for fitness as long as there is no genetic correl-
ation between the measured traits and the unmeasured,
intervening phenotypic character" [9]. The second compo-
nent of adaptive evolutionary change, the evolutionary
response to selection (represented by *R*), depends on
genetic variation (i.e. heritability) and the selection dif-
ferential (breeders' equation) of the trait considered. This
component is difficult to estimate, especially for long-lived
organisms in natural conditions ([11]; see [12] for some es-
timates in different organisms and the problems involved).

In ecological and evolutionary studies, fruit set has often
been considered a convenient proxy of fitness (e.g. [13]).
Common observation of a large variance of reproductive
success (or relative fertility) in natural populations (e.g.
[14]) supports a relevant role of fecundity in short-term
evolutionary change. Indeed, the mass effect created by
some parents' higher reproductive output, through larger
fruit and seed crops, could lead to higher recruitment of
their offspring [15] and, if heritable genetic differences
exist for adaptive traits among parent trees, then short-
term change in population means is expected for those
traits. However, fruit and seed production may have lim-
ited value as predictors of genetic change across genera-
tions because of differences in seed (e.g. less inbred seed)
or offspring (in terms of performance) quality and typically
strong post-dispersal selection [16-20]. In long-lived
organisms, such as forest trees, the seedling stage is a
particularly critical part of the life cycle because seedlings

are more susceptible to resource limitation, which affects
establishment and early growth [21], and suffer from ex-
tremely high mortality rates (up to 90% during the first
growing season [22,23]). Seedlings are thus expected to
be under strong selective pressure in natural condi-
tions [24]. Success in recruitment, then, results from a
complex combination of abiotic and biotic factors that
are spatially variable, which reduces the utility of fruit
set as fitness proxy and highlights the need to include
seedling performance under contrasted environments
in evolutionary change studies. This is particularly
true in highly heterogeneous landscapes such as the
Mediterranean or tropical forests [20].

In this study, we used *Pinus pinaster* Aiton (maritime or cluster pine) as experimental model species to study short-term evolutionary change under contrasted environmental conditions. Maritime pine is a long-lived organism (up to ~200 years under natural conditions) with a large distribution range and important ecological and economic value in the western Mediterranean basin. Large differences in effective reproductive success and molecular-based (phenotypic) selection gradients for different mother traits (as in [8], but considering successfully established regeneration) have been reported in a typical Mediterranean natural population (Coca, Spain) using neighborhood mating models [25,26]. However, offspring performance of the most successful mothers and, thus, the role of observed differences in female reproductive success on short-term genetic change is unknown.

The objectives of our study are i) to test whether higher recruitment of the most successful mothers is due to better performance of their offspring or, alternatively, simply a result of the mass effect of higher female fertility, and ii) to evaluate genetic changes in quantitative traits across generations at different life stages (mature trees and seedlings) and their expected role for adaptability to future environmental conditions. We addressed these goals by estimating breeding values of mother trees with contrasted effective (i.e. based on successful regeneration) reproductive success via progeny testing under controlled and semi-natural conditions. Our experiment protocols allowed us to hasten ontogenetic change and decompose growth into different components that have a high genetic control [27-29]. As detectable response to selection may depend on experimental conditions, we tested progenies under different levels of water availability and intra-specific competition, two highly-relevant selective drivers in Mediterranean maritime pine forests. Finally, short-term genetic change was evaluated by estimating phenotypic (mature trees) and genetic (seedlings) selection gradients and response to selection for a wide range of adaptive traits.
**Methods**

**Plant material and seed collection**

We sampled 22 mature trees (see below for details on tree selection) from a natural maritime pine stand located in a continuous forest in Central Spain (Coca, Segovia). No forest management has been applied to this area in the last 40 years except for removal of dead trees. Mature, well-formed cones were collected from each tree for progeny testing. The cones were collected for the same year, to avoid interannual environmental effects, and for each mother tree separately. Thus, each seed lot had an approximate structure of half-sib families (since correlated paternity within maternal progenies is known to be negligible in maritime pine, see [30]). Seeds were extracted from several cones from each mother and stored in dry conditions at 4°C.

**Reproductive success and field measurements**

Twenty-two mother trees were selected to cover the variation in female reproductive success found in the stand. Previously [26], a study on gene dispersal and female reproductive success based on parentage analysis of naturally-regenerated saplings using nuclear microsatellite (nuSSR) molecular markers had been conducted. This study included 380 mature trees sampled in an area of 3.2 ha and 267 naturally-established saplings, considered the basis for next generation's population, collected from a central subplot of 30 m in radius. Female parentage LOD scores were retrieved from this study. Briefly, the most-likely parents and parent pairs were detected using log-likelihood ratios (LOD-scores) [31,32] and population allele frequencies estimated from the whole dataset (n = 581). Only highly reliable parent-offspring matches (LOD-score > 4.4, as shown by simulation) were used (see details in [32,33]). Following González-Martínez et al. [26], the parent closer to the offspring was considered to be the mother in parent pairs, also because pines present a much more restricted seed than pollen dispersal. LOD scores for each mother-offspring relationship were summed across all offspring for each mother and then divided by the sum of LOD scores across all mothers, producing a relative measure of female reproductive success for each mother tree. Estimates of reproductive success based on fractional paternity assignment are normally less biased than those relying on direct parent-offspring identification [34]. To remove biases due to distance effects, individual reproductive success values were corrected by the probability of a propagule dispersed from the focal female landing within the central subplot where saplings were sampled (assumed for mathematical convenience to be a square of equal area than the actual central subplot). This probability was calculated via numerical integration of a bivariate exponential-power dispersal kernel with parameters α = 68.6 and β = 3, which are the best out of three inverse-modeling dispersal kernels for saplings fitted in [26].

For each of the 22 mother-trees with contrasting values of relative reproductive success (see Additional file 1: Table S1 in Supplementary Information for a summary), we measured different fitness-related traits: size (measured by total height –HT– and diameter –DBH–), diameter increment in the last 10 years (DBH10), number of cones, seed weight (based on a 100-seed sample), germination rate (after 35 days), and number of empty seeds (in a 100-seed sample).

**Experimental design for progeny testing**

To estimate breeding values, two different experiments were conducted using progenies from mother trees with contrasting effective reproductive success. The first one (Experiment 1) was a growth-chamber experiment under environment-controlled conditions, the second (Experiment 2), an outdoor experiment under semi-natural conditions. Individual identification of the seeds and seedlings was maintained for the duration of the experiments.

**Experiment 1: Environment-controlled conditions**

A sample of 100 seeds from each of 18 mother trees (see Additional file 1: Table S1 in Supplementary Information) were weighed and germinated under controlled conditions (20°C, 12 h light) using a sand-substrate. After germination, 24 seedlings per family were transplanted to individual 400 cc pots, with a sand-peat (2:1) substrate, and were moved to a growth chamber under controlled conditions. Prior to transplanting, each pot, including the substrate, was weighed (dry weight and full capacity weight) to determine water content at full capacity. The plants were arranged in four blocks (two for each watering regime, see below) and a row-column design (within blocks), with six experimental units (one-seedling plot) for each genetic entry, hence totaling 24 seedlings per family. The conditions of the experiment were set to hasten ontogenetic changes [27], aiming at maximizing developmental differences in a short cultivation time. This, in addition, helped to control root restriction effects due to container to a minimum [35]. The experimental phase was divided in three (acclimation, growing –including the watering treatment period–, and hardening), depending on temperature and photoperiod (see Additional file 1: Table S2 in Supplementary Information). The first watering regime consisted on keeping plants at full-field capacity (W1, 45% in volume), and the second on reducing water availability to 20% of full-field capacity (W2, 10% in volume). The two watering regimes were controlled by weighing two pots per family and block every second day. In total, the experiment lasted for 189 days (84 days under watering treatments).

**Experiment 2: Outdoor semi-natural conditions**

Twenty-four seeds from each of 21 mother trees (see Additional file 1: Table S1 in Supplementary Information)
were weighed and germinated under controlled conditions (20°C, 12 h light) using a sand substrate and transplanted to outdoor semi-natural conditions. Transplanting was done periodically by blocks, in order to keep the differences due to germination included in the block effect. The experimental site was located in Madrid (40° 27' 23.82" N, 3° 45' 6.91" W, 597 m of altitude) under average environmental conditions similar to those of the seed source stand (rainfall: 452 mm, summer precipitation: 58 mm, temperature: 14.4°C, mean temperature of the coldest month: 5.3°C, mean temperature of the warmest month: 24.8°C). Actual climatic conditions during the outdoor experiment were extremely dry but similar to those in which the previous generation (i.e. the mother trees) of the population studied underwent natural selection (period 1903–1942, data not shown). As in the previous experiment, the plants were arranged in four blocks and a row-column design (within blocks), with six experimental units (one-seedling plot) for each genetic entry (24 seedlings per family). We used two contrasting planting densities; two blocks were planted at a 10x10 cm spacing (D1, 100 plants/m²), and the other two at 5x5 cm (D2, 400 plants/m²), to elicit contrasting inter-family competition effects. Seedlings were protected to avoid predation by animals (rodents, birds, or ants), and no special management was applied during the experiment duration. The experiment lasted 203 days, from April 12th to November 1st.

### Phenotype measurements in seedlings

Phenotypic traits were chosen for their functional and adaptive importance in pines under limiting resource supply: growth and biomass allocation [36-38], specific leaf area [39] and ontogenetic change [40]. These traits present large differences among maritime pine families and populations evaluated in controlled conditions [41-43]. In addition, most of them seemed to have undergone adaptive genetic differentiation (as evaluated by $Q_{ST}$) in response to local climatic conditions, when tested in rangewide experiments [44].

For the two experiments, total height and shoot developmental status were evaluated in all seedlings every second week, the latter following a subjective scale based on the type of axillary shoots, as described in detail in [42]. The height of the plant at ontogenetic stage 3 (HON3) – corresponding to the earliest observation of axillary long shoots– was used for subsequent analyses as a size-independent quantitative index of seedling development. Higher values of HON3 indicate a slower seedling ontogenetic development. Height growth of each seedling was decomposed by adjusting a Gompertz model [29] and analyzing the three main parameters of this curve: HTOT, b, and m. HTOT is the total height at the end of the experiment (asymptote), the b-parameter is related to the maximum growth rate standardized by HTOT, and m, the

inflection point, indicates whether maximum growth is reached at the beginning (low values of m) or the end (high values of m) of the growing period.

At the end of each experiment, survival (SUR) was recorded, all plants were harvested, and dry mass of each fraction (needles, stem, and roots) was determined with 0.001 g precision. The dry mass of each fraction was expressed relative to that of the total plant (leaf mass fraction, LMF; stem mass fraction, SMF; root mass fraction, RMF) [36]. Total dry mass (DW) was obtained as the sum of leaf, stem, and root dry mass. Specific leaf area (SLA) was estimated from 10 needles randomly chosen from each plant. Finally, since at the seedling stage every metamer bears a primary needle, the number of stem units (NSU) was estimated by comparing total primary needle dry mass to dry mass of the 10 needles collected for SLA assessment.

### Quantitative genetic models in progeny tests

For each variable and experiment, a mixed model was adjusted using Restricted Maximum Likelihood (REML):

$$y_{ijkl} = \mu + \phi_i + T_j + \chi_{ij} + B(T)_{k(j)} + c + \varepsilon_{ijkl}$$

where $y_{ijkl}$ is the phenotypic value of the variable for the $l^{th}$ tree from the $i^{th}$ family in the $j^{th}$ treatment located in the $k^{th}$ block, $\mu$ is the overall mean, $\phi_i$ is the effect of the $i^{th}$-family ($i$ varying from 1 to 18–21, depending on the experiment), $T_j$ is the effect of the $j^{th}$ treatment (two watering levels in the growth chamber experiment and two density levels in the outdoor experiment, thus $j$ varying from 1 to 2 in each experiment), $\chi_{ij}$ is the family by treatment interaction, $B_k$ is the effect of the $k^{th}$-block (1–4), $c$ is the covariate effect, and $\varepsilon_{ijkl}$ the residual. For survival, a generalized logistic mixed model (GLMM) was applied with a binomial distribution and a logit link function. For biomass traits (LMF, SMF, RMF), the total dry mass (DW) was used as a covariate to correct for allometric changes [36]. For the rest of the quantitative traits, seed weight (WSEED) was used as a covariate to control for maternal effects. Family and treatment by family interaction were considered random factors. The models were fitted using ASREML [45].

### Genetic control

Narrow-sense heritability ($h^2$) was estimated for each trait and progeny test following standard procedures [46]. For survival (SUR), the residual variance of the GLMM was set to $\pi^2/3$ [47]. The standard error of heritability was estimated using a Taylor’s series approximation [45]. Best Linear Unbiased Predictors (BLUPs) are weighted means commonly used as point estimates of random effects in mixed models. They predict the expected phenotype of a tree’s offspring or an individual by using phenotypic
information collected from relatives, which is useful to evaluate genetic changes in the populations after selection. Mother-tree breeding values (based on offspring performance) for each trait were computed from the heritability of the trait and the BLUPs estimated in the experiments [45,46]. Evolvability, as estimated by the coefficient of additive genetic variation (CV<sub>A</sub>), was also computed [10].

## Phenotypic selection gradients and correlation of female reproductive success with offspring performance

Phenotypic selection gradients were computed based on traits measured directly on the mother-trees in the field, to evaluate the relationship among trait values in adults and their fitness (in terms of relative reproductive success), as obtained from molecular marker-based parentage analysis (see [26] and above). The proxy for fitness of parent trees used in this study (i.e. the actual number of successful naturally-regenerated offspring) is probably more accurate than amount of pollen donated, flower production or fruit set (e.g. [8,13]), as the latter three do not consider post-dispersal mortality and early selection. These selection gradients were computed as the vector of partial regression coefficients of individual relative reproductive success on trait values. Given the reduced number of mother trees in the study (22), we performed a Principal Component analysis to reduce the number of dimensions in the analysis. The following log-linear model was used ([7] and references therein) to define female relative reproductive success ($\lambda_k$):

$$Ln(\lambda_k) = \beta_1 z_{1k} + \beta_2 z_{2k} + \beta_3 z_{3k} + \gamma_1 z_{1k}^2 + \gamma_2 z_{2k}^2 + \gamma_3 z_{3k}^2$$

The subscripts (1 to 3) refer to each of the three first principal components (explaining 69.64% of the variance, see Additional file 1: Table S4 in Supplementary Information). In biological terms, the $\beta_i$ are analogous to linear selection gradients, and the $\gamma_i$ are analogous to quadratic selection gradients as defined by Lande and Arnold [6]. Traits were standardized, and individual fitness was converted to relative fitness by dividing by the mean fitness of all individuals within the subset of data. The statistical significance of each specific selection gradient, $\beta_i$ or $\gamma_i$, was estimated using likelihood-ratio tests, by subtracting the log-likelihood for the model excluding each parameter, one at a time, from the log-likelihood of the full model, which is asymptotically chi-square distributed with one degree of freedom.

Finally, to know to what extent mother reproductive success (and the associated phenotypes) is driving population phenotypic change across generations, we computed the relationship between female reproductive success and breeding values for each trait, as obtained from offspring performance under the two experimental environments,

using log-linear models. Insights into this relationship are especially relevant at early establishment stages, when natural selection is stronger.

## Genetic selection gradients and response to selection at early stages of establishment

Genetic selection gradients at seedling stage were computed in the two experiments using seedling height (HTOT), which is related to competitive ability and crucial for natural regeneration establishment in pines [48], and/or survival (SURV) as fitness proxies. Both fitness proxies were used in the outdoors semi-natural experiment while only seedling height was used to gauge fitness under the optimal growing conditions of the controlled-environment experiment, as no significant genetic differences among families were found for survival in this experiment (see Results). Genetic selection gradients were obtained by bivariate analysis (to reduce the bias in the estimation, see [9,12]), as the genetic covariance between the trait and fitness divided by the genetic variance of the trait, for the traits estimated in the two experiments. In addition, we computed the genetic response to selection (R), following the Price Theorem [5], as the covariance between the trait and fitness. Standard errors were computed using a Taylor's series approximation [45].

## Results

### Genetic control, and interactions and main effects of water regime and density

Seedlings under semi-natural conditions almost doubled the specific leaf area (SLA) of those under controlled conditions (Experiment 1), and also had higher height and biomass (with a higher allocation to the stem fraction), but the development of axillary shoots (i.e. the same onto-genetic stage) was reached at a similar size (HON3). Genetic control varied depending on the trait and experiment (Table 1). In the growth chamber experiment (Experiment 1), narrow-sense heritability ($h^2$) was significantly different from zero for all traits except for b, the standardized maximum growth rate (obtained from decomposition of each seedling height growth using a Gompertz model), and survival (SUR). Under semi-natural conditions (Experiment 2), genetic control was much lower except for maturation (HON3), total height growth (HTOT) and survival (SUR). Strong genetic control of HTOT and SUR was somehow expected, as these two parameters are directly associated with fitness at early stages of establishment in forest trees. A positive correlation ($r = 0.694$) was observed between HTOT and SUR breeding values under growth chamber conditions, but this correlation was strongly negative ($r = -0.637$) under semi-natural ones (see Additional file 1: Table S3 in Supplementary Information). Finally, genetic
<table><caption>Table 1 Genetic control for adaptive traits in maritime pine under two contrasted experimental conditions</caption><thead><tr><th rowspan="2">Trait</th><th colspan="4">Environment-controlled conditions (Experiment 1)</th><th colspan="4">Outdoor semi-natural conditions (Experiment 2)</th></tr><tr><th>h<sup>2</sup></th><th>S.E.(h<sup>2</sup>)</th><th>MEAN</th><th>CV<sub>A</sub></th><th>h<sup>2</sup></th><th>S.E.(h<sup>2</sup>)</th><th>MEAN</th><th>CV<sub>A</sub></th></tr></thead><tbody><tr><td>LMF</td><td>0.350</td><td>0.163</td><td>0.541</td><td>3.4</td><td>0.068</td><td>0.072</td><td>0.469</td><td>5.3</td></tr><tr><td>SMF</td><td>0.357</td><td>0.171</td><td>0.055</td><td>6.3</td><td>0.041</td><td>0.064</td><td>0.163</td><td>12.0</td></tr><tr><td>RMF</td><td>0.332</td><td>0.158</td><td>0.403</td><td>4.8</td><td>0.137</td><td>0.091</td><td>0.367</td><td>5.4</td></tr><tr><td>DW</td><td>0.377</td><td>0.170</td><td>0.892</td><td>5.5</td><td>0.110</td><td>0.083</td><td>0.485</td><td>29.1</td></tr><tr><td>SLA</td><td>0.253</td><td>0.162</td><td>57.035</td><td>2.5</td><td>0.071</td><td>0.073</td><td>82.5</td><td>4.7</td></tr><tr><td>NSU</td><td>0.611</td><td>0.246</td><td>37.287</td><td>8.4</td><td>0.066</td><td>0.071</td><td>85.8</td><td>19.6</td></tr><tr><td>HON3</td><td>0.137</td><td>0.103</td><td>24.951</td><td>1.9</td><td>0.276</td><td>0.126</td><td>27.48</td><td>9.7</td></tr><tr><td>b</td><td>0.000</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>m</td><td>0.297</td><td>0.157</td><td>36.328</td><td>0.9</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>HTOT</td><td>0.876</td><td>0.2785</td><td>19.231</td><td>8.8</td><td>0.415</td><td>0.159</td><td>46.94</td><td>13.1</td></tr><tr><td>SUR</td><td>0.226</td><td>0.2300</td><td>0.767</td><td>NA</td><td>0.719</td><td>0.243</td><td>0.715</td><td>NA</td></tr></tbody></table>

LMF: leaf mass fraction; SMF: stem mass fraction, RMF: root mass fraction; DW: dry weight; SLA: specific leaf area; NSU: number of stem units; HON3: height at first production of axillary shoots; b, m: coefficients of the Gompertz function; HTOT: total height; SUR: Survival; NA: not available.

h²: heritability; S.E. (h²): standard error of h²; MEAN: mean value of the population; CVA: coefficient of additive genetic variation.

Heritability estimates are shown only when family effects were statistically significant.

values were correlated across experiments only for HTOT (positive, r = 0.418) and SLA (negative, r = -0.478).

In both experiments and for all traits, the family per treatment (water regime in *Experiment 1* and density in *Experiment 2*) interaction term had very little importance (less than 5% of the total variance; Table 2). The lack of family per treatment interaction for survival and other traits related to fitness at early life stages, such as total height (HTOT), indicated reduced opportunities for family specialization caused either by intraspecific competition or by watering regime under our experimental conditions.

With respect to main effects of treatments, a highly significant effect of the watering regimes was found in *Experiment 1* for all traits except for maximum height growth rate (parameter *b* of the Gompertz function) (Table 2). The covariates, which accounted for allometric changes (dry weight, DW) and maternal effects (seed weight, WSEED) in this experiment, were also highly significant for most traits (except for stem mass fraction, SMF, specific leaf area, SLA, and ontogenetic change, HON3). Under the semi-natural outdoor conditions (*Experiment 2*), seeding density had little influence on biomass allocation traits,

<table><caption>Table 2 Mixed model analysis of two experiments (under controlled and semi-natural conditions) in maritime pine, including fixed effects (F-value and significance level) and covariates (DW: dry weight, WSEED: seed weight)</caption><thead><tr><th rowspan="3">Trait</th><th colspan="6">Environment-controlled conditions (Experiment 1)</th><th colspan="6">Outdoor semi-natural conditions (Experiment 2)</th></tr><tr><th colspan="3">F-value, significance level<sup>1</sup></th><th rowspan="2">%Var</th><th colspan="2">Mean<sup>2</sup></th><th colspan="3">F-value, significance level<sup>1</sup></th><th colspan="2">Mean<sup>2</sup></th><th rowspan="2">Covariate</th></tr><tr><th>TRT</th><th>BLOCK</th><th>COV</th><th>TRT*FAM</th><th>W1</th><th>W2</th><th>TRT</th><th>BLOCK</th><th>COV</th><th>TRT*FAM</th><th>D1</th><th>D2</th></tr></thead><tbody><tr><td>LMF</td><td>28.70 ***</td><td>7.43 **</td><td>20.19 ***</td><td>0.0 ns</td><td>0.431</td><td>0.486</td><td>0.04 ns</td><td>12.14 ***</td><td>58.01 ***</td><td>1.4 ns</td><td>0.477</td><td>0.461</td><td>DW</td></tr><tr><td>SMF</td><td>9.00 **</td><td>1.81 ns</td><td>0.61 ns</td><td>3.7 ns</td><td>0.054</td><td>0.059</td><td>0.01 ns</td><td>8.71 **</td><td>24.73 ***</td><td>0.0 ns</td><td>0.161</td><td>0.165</td><td>DW</td></tr><tr><td>RMF</td><td>35.24 ***</td><td>7.96 **</td><td>16.16 ***</td><td>0.0 ns</td><td>0.514</td><td>0.454</td><td>0.03 ns</td><td>10.85 ***</td><td>19.60 ***</td><td>1.0 ns</td><td>0.361</td><td>0.379</td><td>DW</td></tr><tr><td>DW</td><td>280.23 ***</td><td>62.17 ***</td><td>29.62 ***</td><td>0.0 ns</td><td>0.7406</td><td>1.0519</td><td>7.13 **</td><td>4.69 **</td><td>9.08 **</td><td>0.0 ns</td><td>0.4257</td><td>0.5725</td><td>WSEED</td></tr><tr><td>SLA</td><td>77.76 ***</td><td>26.80 ***</td><td>1.21 ns</td><td>0.0 ns</td><td>49.59</td><td>56.70</td><td>12.26 ***</td><td>3.67 *</td><td>1.72 ns</td><td>2.8 ns</td><td>84.99</td><td>79.16</td><td>WSEED</td></tr><tr><td>NSU</td><td>40.56 ***</td><td>40.90 ***</td><td>15.82 **</td><td>3.8 ns</td><td>36.53</td><td>45.30</td><td>2.27 ns</td><td>1.99 ns</td><td>9.71 **</td><td>0.0 ns</td><td>82.13</td><td>90.88</td><td>WSEED</td></tr><tr><td>HON3</td><td>7.73 **</td><td>0.28 ns</td><td>3.16 ns</td><td>2.6 ns</td><td>25.91</td><td>27.61</td><td>21.80 ***</td><td>6.60 **</td><td>9.10 **</td><td>0.0 ns</td><td>26.07</td><td>29.45</td><td>WSEED</td></tr><tr><td>b</td><td>0.03 ns</td><td>1.06 ns</td><td>5.22 *</td><td>0.0 ns</td><td>0.0637</td><td>0.0663</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>WSSED</td></tr><tr><td>m</td><td>101.79 ***</td><td>3.57 *</td><td>10.59 **</td><td>2.6 ns</td><td>191.3</td><td>204.2</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>WSEED</td></tr><tr><td>HTOT</td><td>159.48 ***</td><td>4.48 *</td><td>5.01 *</td><td>0.0 ns</td><td>34.47</td><td>42.45</td><td>8.73 **</td><td>2.55 ns</td><td>16.60 ***</td><td>0.0 ns</td><td>44.82</td><td>50.31</td><td>WSEED</td></tr><tr><td>SUR</td><td>14.49 **</td><td>7.11 **</td><td></td><td>4.4 ns</td><td>0.827</td><td>0.594</td><td>0.08 ns</td><td>1.22 ns</td><td></td><td>0.0 ns</td><td>0.722</td><td>0.697</td><td></td></tr></tbody></table>

TRT: treatment (W1/W2, D1/D2; see Methods); COV: covariate; %Var: Percent of the interaction term (TRT*FAM) with respect to the total variance; LMF: leaf mass fraction; SMF: stem mass fraction; RMF: root mass fraction; DW: dry weight; SLA: specific leaf area; NSU: number of stem units; HON3: height at first production of axillary shoots; b, m: coefficients of the Gompertz function; HTOT: total height; SUR: survival; NA: not available.

¹ns: not significant; *,: 0.05 ≤ α < 0.01, **: 0.01 ≤ α < 0.001, ***: α ≤ 0.001.

²Best linear unbiased estimates (BLUE) for watering (W1 and W2) and density (D1 and D2).
survival (SUR) and number of stem units (NSU), but high effects in dry weight (DW), specific leaf area (SLA), onto-genetic change (HON3) and height (HTOT). The mixed effects of density, i.e. intraspecific competition, on adaptive traits highlight the need to evaluate various components of fitness, in this case SUR and HTOT, in microevolutionary studies. As for Experiment 1, the covariates, dry weight (DW) and seed weight (WSEED), were highly significant for most traits (except for specific leaf area, SLA).

## Phenotypic selection gradients and correlation of female reproductive success with offspring performance

Relative effective reproductive success in natural conditions (i.e. based on successfully established offspring) varied one order of magnitude (0.0078-0.0942) among mother trees (as retrieved from a previous study, see [26]). Selection gradients based on mother-tree phenotypic traits rely on direct field measurements (Table 3). Only the linear phenotypic selection gradient for the PC2 (related to germination rate and number of empty seeds, i.e. seed quality) was significant, with the quadratic selection gradient for this principal component being marginally significant (P = 0.0674).

We found significant linear correlations (of different sign) between female reproductive success and offspring performance (as estimated by breeding values) in two traits: negative for specific leaf area (SLA) in *Experiment 1* and positive for stem mass fraction (SMF) in *Experiment 2* (Figure 1). The quadratic term of the log-linear model was also significant for the former. These results, together with contrasted heritability across experiments, illustrate distinct expected microevolutionary changes due to differences in female reproductive success depending on environmental conditions.

## Genetic selection gradients and response to selection at early stages of establishment

Genetic selection gradients were computed at early establishment stages based on genetic covariance of offspring

traits with total height (in both experiments) and survival (only in the outdoor experiment, *Experiment 2*, as this trait was not under genetic control in *Experiment 1*; see above) as fitness proxies (see Methods). When considering total height, selection gradients were significant for stem mass fraction (SMF, positive) and the inflection point of the growth curve (parameter *m*, negative; i.e. plants that reached maximum growth at the beginning of the growing period showed higher fitness) under environment-controlled conditions (*Experiment 1*), and for the ontogeny index (HON3, positive; i.e. plants with slower ontogenetic development showed higher fitness) under outdoor experimental conditions (*Experiment 2*) (Table 4). Response to selection (*R*) for these variables ranged (in absolute values) from 0.08 (*m*) to 0.12 (SMF and HON3). When survival (SUR) was considered as a fitness proxy in the outdoor experimental conditions (*Experiment 2*), the results were similar, with HON3 being the only trait with a significant (and positive) effect on fitness.

## Discussion

This paper explores the potential for micro-evolutionary changes in a natural population of *Pinus pinaster* at two different stages of development: mature trees and seedlings. To our knowledge, there are no studies on the estimation of such changes in long-lived plants, such as forest trees, while some attempts have been made in other organisms [11,12,49]. By using an experimental approach under different environmental conditions, we were able to analyze the phenotypic selection gradients and the correlation between female reproductive success and breeding values, as well as the genetic selection gradients and cross-environmental correlations at the seedling stage, when selection intensity is higher in natural forest tree populations.

Overall, our results suggest that significant genetic change is expected across generations in maritime pine and that those changes may strongly depend on local environmental conditions.

<table><caption>Table 3 Phenotypic selection gradients for maritime pine based on mother-tree traits</caption><thead><tr><th>Trait</th><th>DF</th><th>Selection gradient</th><th>Standard Error</th><th>Wald 95% CIs</th><th>Chi-Square</th><th>Pr > ChiSq</th></tr></thead><tbody><tr><td colspan="7"><em>Linear selection gradient</em></td></tr><tr><td>PC1</td><td>1</td><td>0.2140</td><td>0.1366</td><td>-0.0538</td><td>0.4818</td><td>2.45</td><td>0.1173</td></tr><tr><td>PC2</td><td>1</td><td>0.4313</td><td>0.2045</td><td>0.0306</td><td>0.8321</td><td>4.45</td><td>0.0349</td></tr><tr><td>PC3</td><td>1</td><td>-0.0289</td><td>0.1247</td><td>-0.2733</td><td>0.2156</td><td>0.05</td><td>0.8170</td></tr><tr><td colspan="7"><em>Quadratic selection gradient</em></td></tr><tr><td>PC1<sup>2</sup></td><td>1</td><td>-0.0073</td><td>0.0802</td><td>-0.1645</td><td>0.1499</td><td>0.01</td><td>0.9273</td></tr><tr><td>PC2<sup>2</sup></td><td>1</td><td>0.1090</td><td>0.0596</td><td>-0.0078</td><td>0.2259</td><td>3.35</td><td>0.0674</td></tr><tr><td>PC3<sup>2</sup></td><td>1</td><td>-0.1609</td><td>0.0974</td><td>-0.3518</td><td>0.0301</td><td>2.73</td><td>0.0987</td></tr></tbody></table>

Significant selection gradients are given in italics (Additional file 1: Table S4 for details of the PC analysis).

PC1 main loading factors: age (0.813), height (0.798) and diameter (0.790); PC2 main loading factors: number of empty seeds (-0.768) and germination rate (0.896); PC3 main loading factors: seed weight (0.894).
<figure><img src="image_2.png" /><figcaption>**Figure 1** Relationship between female relative reproductive success and breeding values for a) Specific leaf area (SLA) and b) Stem mass fraction (SMF) measured in environment-controlled conditions (Experiment 1) and outdoor semi-natural conditions (Experiment 2), respectively. Parameters and significance values corresponding to the log-linear model are included; ns: not significant, *: 0.05 ≤ α < 0.01, **: 0.01 ≤ α < 0.001.</figcaption></figure>

<table><caption>**Table 4** Genetic selection gradients for offspring (i.e. based on genetic covariance of traits with seedling height – HTOT- and/or survival –SUR- as proxies of fitness; see Methods) and response to selection (R) in two experimental conditions in maritime pine</caption><thead><tr><th rowspan="3">Trait</th><th colspan="6">Environment controlled conditions (Experiment 1)</th><th colspan="6">Outdoor semi-natural conditions (Experiment 2)</th></tr><tr><th colspan="3">HTOT</th><th rowspan="2">S.E.</th><th rowspan="2">R</th><th rowspan="2">S.E.</th><th colspan="3">HTOT</th><th colspan="3">SUR</th></tr><tr><th>Selection gradient</th><th>S.E.</th><th>R</th><th>Selection gradient</th><th>S.E.</th><th>R</th><th>S.E.</th><th>Selection gradient</th><th>S.E.</th><th>R</th><th>S.E.</th></tr></thead><tbody><tr><td>LMF</td><td>0.2466</td><td>0.4444</td><td>0.020</td><td>0.037</td><td>NE</td><td>NE</td><td>0</td><td>0</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td></tr><tr><td>SMF</td><td>1.4059</td><td>0.3697</td><td>0.121</td><td>0.048</td><td>-1.2667</td><td>5.3920</td><td>-0.008</td><td>0.027</td><td>2.1487</td><td>12.2429</td><td>0.012</td><td>0.062</td></tr><tr><td>RMF</td><td>-0.6386</td><td>0.4441</td><td>-0.048</td><td>0.038</td><td>0.3422</td><td>1.4835</td><td>0.007</td><td>0.029</td><td>-2.0826</td><td>2.9119</td><td>-0.052</td><td>0.069</td></tr><tr><td>DW</td><td>1.0307</td><td>0.5029</td><td>0.048</td><td>0.030</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>3.3621</td><td>1.1801</td><td>0.289</td><td>0.104</td></tr><tr><td>SLA</td><td>0.8369</td><td>0.7142</td><td>0.039</td><td>0.033</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td></tr><tr><td>NSU</td><td>0.7501</td><td>0.3275</td><td>0.078</td><td>0.045</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>4.3425</td><td>2.2636</td><td>0.217</td><td>0.089</td></tr><tr><td>HON3</td><td>1.5480</td><td>0.898</td><td>0.052</td><td>0.033</td><td>1.0725</td><td>0.2516</td><td>0.120</td><td>0.052</td><td>1.8632</td><td>0.6731</td><td>0.229</td><td>0.119</td></tr><tr><td>m</td><td>-1.5200</td><td>0.5959</td><td>-0.082</td><td>0.036</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td><td>NE</td></tr></tbody></table>

LMF: leaf mass fraction; SMF: stem mass fraction, RMF: root mass fraction; DW: dry weight; SLA: specific leaf area; NSU: number of stem units; HON3: height at first production of axillary shoots; m: the inflection point of the growth curve.

S.E. standard error of the estimates. Significant values are given in italics; NE: not estimable.
## Genetic control and environmental effects

Heritable transmission of traits across generations, a key component of evolutionary change [5], differed across environments, with much higher heritability under climatic chamber conditions than in the outdoor experiment –suggesting more potential for genetic change in less stressful environments. The only exception was for ontogenetic heteroblastic change (HON3), a trait that correlates with seedling field survival under drought in other pine species [50], and that was also involved in genetic selection gradients at seedling stage in our study (see below). Evolvability [10] was, however, much higher under outdoor semi-natural conditions, but also very variable across traits (see Table 1). Therefore, further studies on the potential for evolutionary change in forest trees need to include a wide range of traits and environments, especially as climatic models for trees predict different levels of regional maladaptation in the near future [51].

Despite its expected role as microevolutionary driver in heterogeneous environments [52], we did not find genetic variation for phenotypic plasticity (i.e. the treatment per family interaction) for any trait or treatment, either water regime or density. Lack of genetic differences in phenotypic plasticity for key adaptive traits, such as biomass allocation and survival, could be related to the cost of plasticity [52,53] that may limit its evolutionary importance in Mediterranean pine seedlings. However, one more likely explanation is the difficulty of the genotypes to track the optimum in complex and heterogeneous environments, as highlighted by the unexpected phenotypic changes in response to the treatments in our experiments: seedlings under the limited watering regime had lower root biomass allocation and higher height than under full water-capacity, while seedlings under the high density treatment had higher total biomass and height than under low density. Mediterranean environments are highly variable [54], and interactions among drought tolerance, intra- and inter-specific competition, and light use can result in variable life strategies among populations and individuals [55,56].

## Determinants of female effective reproductive success

Higher recruitment from the most successful mothers is mainly due to their seed quality. Significant and positive linear phenotypic selection gradients were detected at mature stage (parent trees) for the second principal component (PC2), mostly related to germination rate and number of empty seeds. Seed quality and parental phenotype are often correlated in plants, and a clear relation between mother-tree size, and seed output and quality exists [57], in agreement with sexual-selection theory predicting that female reproduction is limited by physical resources [58]. The phenotypic selection gradient for seed quality traits points also to a potential role of inbreeding. However, maritime pine is generally less affected by inbreeding than

other pines in terms of general vigor [59], and often exhibits low inbreeding coefficients in natural populations, as estimated by molecular markers [60]. Despite phenotypic traits being reduced to independent PCAs, phenotypic selection gradients detected using Lande & Arnold's [6] multivariate fitness analysis could be false-positives resulting from correlation with traits under selection that were not measured. However, this is unlikely in our study, as we included most of the adaptive traits that have shown differences for Mediterranean pines at the individual and population levels in several previous quantitative genetic studies [39,42-44].

The mass effect of seed quality on female effective reproductive success is modified by a distinct recruitment across families, seemingly derived from better adapted seedlings. Significant correlations among female relative reproductive success and breeding values (based on offspring performance) were detected for specific leaf area (SLA) under environment-controlled conditions and stem mass fraction (SMF) under semi-natural outdoor conditions, two adaptive traits related to abiotic stress-response in pines. The functional significance of these traits and their performance in the two environments provided insights into some potential mechanisms of evolutionary change in maritime pine. SLA was negatively correlated across the two experiments, indicating a different genetic basis. The lower SLA under controlled conditions for offspring from the mothers with higher reproductive success, could be related to the production of shorter needles with less transpirational surface area, while this same trait under the more stressful conditions of the semi-natural experiment may represent a largely plastic response to drought stress (as indicated by the lower heritability). SLA is related to the modulation of photosynthetic exchange (at equal levels of N in the needles [61]) and thus is a main component of relative growth rates. The adaptive importance of SLA under different stressful conditions is also well-known [62,63] and negative selection gradients for this trait have been found in other plants (e.g. [64]). With respect to SMF, changes in biomass allocation have been described as a strategy for drought tolerance in pine [38,44] and is correlated with fitness proxies (HTOT, see Methods) in genetic selection gradients under controlled conditions (see below). But under semi-natural conditions, despite significant genetic differences among families (see Table 1), the expected evolutionary potential is reduced due to the low genetic control of the trait ($h^2$ of 0.041). Such small genetic effects are expected in natural environments [65]. Variability in the correlation among reproductive success and functional or adaptive traits due to offspring's environment can foster intra-population genetic diversity, even in the absence of spatial heterogeneity, as recruitment in forest trees depends largely on highly-variable annual environmental conditions [23].
## Genetic change and responses to selection

Genetic selection gradients (i.e. those based on genetic co- variance between the trait and fitness) using seedling height as a fitness proxy indicated a different microevolu- tionary pattern between the two experimental conditions. Higher fitness is expected from plants with higher biomass allocation to stems (SMF, positive selection gradient) and a greater growth at the beginning of the growing season ( *m* negative selection gradient) under (mild) environment- controlled conditions, but under more stressful outdoor conditions, seedlings with delayed ontogenetic develop- ment (HON3, positive selection gradient) will perform better. A strong negative correlation between height and survival in the outdoors experiment also suggests that a slow-growing drought-efficient phenotype is favored under stressful conditions. A bias in the estimation of selection gradients, which could have affected our esti- mates across the two environmental conditions [21], is unlikely, as we used a bivariate analysis [9], so these environment-dependent shifts in correlation among traits and fitness are rather the likely outcome of fluc- tuating selection in variable environments [66-68]. Low cross-environmental correlations in heterogeneous envi- ronments (see Additional file 1: Table S3) favor a more rapid approach to the optima and a more independent (less constrained) trait evolution [59].

The response to selection, *R*, associated with traits in- volved in significant genetic selection gradients (see above) was moderate when height was used as fitness proxy (~8- 12%), but much higher when survival under semi-natural conditions was considered (~23%). In particular, ontogen- etic change (HON3), a trait under relatively high genetic control ($h^2$ of 0.276), appeared to be involved in microevo- lutionary change under these close-to-nature conditions. In Mediterranean environments, drought stress and inter- specific competition play overarching roles [23,55] and, therefore, individuals with a delayed developmental status (i.e. that reach a bigger size before ontogenetic changes take place) could be more efficiently established, as sug- gested by our results.

## Conclusions

Microevolutionary change due to differences in maternal reproductive success and early selection is expected in Mediterranean maritime pine populations, and could be an important mechanism mitigating the negative conse- quences of climate change. Substantial genetic variation was detected for adaptive traits related to growth, onto- genetic change and biomass allocation in maritime pine, but no clear differences in phenotypic plasticity among seedlings when two different factors (water availability and density) were considered. At mature-tree stages, higher fe- male effective reproductive success can be mainly explained by a mass effect (seed quality), but was also modulated by

the production of better-adapted seedlings. Selection gradi- ents and responses to selection for seedlings differed across experimental conditions. Thus, although we have collected evidence of short-term genetic change for some relevant adaptive traits, directional changes are difficult to predict due to variable selection at different life stages and environ- ments. In fact, contrasting experimental results highlight the relevance of testing genetic value under very different conditions to assess the possible trajectories of the re- sponse to selection. The distinct processes involved at the two life stages (mature trees or seedlings) together with environment-specific responses advise caution when predicting likely evolutionary responses to environmental change in Mediterranean forest trees.

## Availability of supporting data

The data sets supporting the results of this article are available in the LabArchives repository, doi:10.6070/H4JM27KW.

## Additional file 1

### Additional file 1: Supporting Information.

### Abbreviations

a: Significance level; b: Coefficient of the Gompertz function related to the maximum growth rate; CV<sub>A</sub>: Coefficient of additive genetic variation;

COV: Covariance; D1: Outdoor experiment density (400 plants/m²);

D2: Outdoor experiment density (10,000 plants/m²); DBH: Diameter (mother trees); DBH10: Diameter increment in the last 10 years (mother trees);

DW: Total dry mass; $h^2$: Narrow-sense heritability; HON3: Height at first

production of axillary shoots; HT: Total height of the mother tree; HTOT: Total height of seedlings in experiments; LMF: Leaf mass fraction;

m: Coefficient of the Gompertz function related to the inflection point;

MOR: Mortality; NA: Not available; NE: Not estimable; ns: Not significant;

NSU: Number of stem units; R: Response to selection; RMF: Root mass

fraction; S.E: Standard error; SLA: Specific leaf area; SMF: Stem mass fraction; TRT: Treatment; W1: Watering regime at full-field capacity (45% in volume);

W2: Watering regime at 20% of full-field capacity (10% in volume);

WSEED: Seed weight (mg).

### Competing interests

The authors declare that they have no competing interests.

### Authors' contributions

RA was responsible for study design, data analysis, manuscript writing, and co-funding of the study. RC was responsible for experimental design, data measurements, and manuscript review. EN contributed in data analysis and manuscript review. JC contributed to experimental design and manuscript writing. SCG-M participated in the study design, data analysis, manuscript writing, and co-funding of the study. All authors read and approved the final manuscript.

### Acknowledgments

We are grateful to INIA's field team (F. del Caño, N. Godoy and D. Barba) for experiment establishment and trait measurement, and to J.J. Robledo-Arnuncio for revising an earlier version of the manuscript and for producing a script to correct bias due to distance effects in estimates of female reproductive success. Thanks are extended to an anonymous reviewer whose positive and constructive comments helped us to improve early versions of the manuscript. The text has been revised by a professional scientific editor, P. C. Grant. This research was financially supported by grants from the Spanish Ministry of Education and Science to Santiago C. González-Martínez (project REPROFOR, AGL2005-07440-C02-01), and to R. Alía (Project RTA2010-00120-C02-01).
Author details

Received: 6 June 2014 Accepted: 8 September 2014

Published online: 23 September 2014

References

1. Aitken SN, Yeaman S, Holliday JA, Wang T, Curtis-McLane S: Adaptation, migration or extirpation: climate change outcomes for tree populations. *Evol Appl* 2008, 1:95-111.

2. Petit RJ, Hu FS, Dick CW: Forests of the past: a window to future changes. *Science* 2008, **320**:1450-1452.

3. Whitham TG, Bailey JK, Schweitzer JA, Shuster SM, Bangert RK, Leroy CJ, Lonsdorf EV, Allan GJ, DiFazio SP, Potts BM, Fischer DG, Gehring CA, Lindroth RL, Marks JC, Hart SC, Wimp GM, Wooley SC: A framework for community and ecosystem genetics: from genes to ecosystems. *Nat Rev Genet* 2006, **7**:510-523.

4. Hamrick JL: Response of forest trees to global environmental changes. *For Ecol Manage* 2004, **197**:323-335.

5. Price GR: Selection and covariance. *Nature* 1970, **227**:520-521.

6. Lande R, Arnold SJ: The measurement of selection on correlated characters. *Evolution* 1983, **37**:1210-1226.

7. Wright JW, Meagher TR: Selection on floral characters of natural Spanish populations of *Silene latifolia*. *J Evol Biol* 2004, **17**:382-395.

8. Morgan MT, Conner JK: Using genetic markers to directly estimate male selection gradients. *Evolution* 2001, **55**:272-281.

9. Rausher MD: The measurement of selection of quantitative traits: biases due to environmental covariances between traits and fitness. *Evolution* 1992, **46**:616-626.

10. Houle D: Comparing evolvability and variability of quantitative traits. *Genetics* 1992, **130**:195-204.

11. Kruuk LEB: Estimating genetic parameters in natural populations using the "animal model". *Philos Trans R Soc Lond B* 2004, **359**:873-890.

12. Hadfield JD, Wilson AJ, Garant D, Sheldon BC, Kruuk LEB: The misuse of BLUP in ecology and evolution. *Am Nat* 2010, **175**:116-125.

13. Campbell DR: Effects of floral traits on sequential components of fitness in *Ipomopsis aggregata*. *Am Nat* 1991, **137**:713-737.

14. Oddou-Muratorio S, Klein EK, Austerlitz F: Pollen flow in the wildservice tree, *Sorbus torminalis* (L.) Crantz. II. Pollen dispersal and heterogeneity in mating success inferred from parent-offspring analysis. *Mol Ecol* 2005, **14**:4441-4452.

15. Castro J: Seed mass versus seedling performance in Scots pine: a maternally dependent trait. *New Phytol* 1999, **144**:153-161.

16. Naito Y, Konuma A, Iwata H, Suyama Y, Seiwa K, Okuda T, Lee SL, Muhammad N, Tsumura Y: Selfing and inbreeding depression in seeds and seedlings of *Neobalanocarpus heimii* (Dipterocarpaceae). *J Plant Res* 2005, **118**:423-430.

17. Kärkkäinen K, Savolainen O, Koski V: Why do plants abort so many developing seeds: bad offspring or bad maternal genotypes? *Evol Ecol* 1999, **13**:305-317.

18. Steven DDE, Wright SJ: Consequences of variable reproduction for seedling recruitment in three neotropical tree species. *Ecology* 2011, **83**:2315-2327.

19. Kolb A: Reduced reproductive success and offspring survival in fragmented populations of the forest herb *Phyteuma spicatum*. *J Ecol* 2005, **93**:1226-1237.

20. Nathan R, Casagrandi R: A simple mechanistic model of seed dispersal, predation and plant establishment: Janzen-Connell and beyond. *J Ecol* 2004, **92**:733-746.

21. Leck MA, Parker VT, Simpson RL: *Seedling Ecology and Evolution*. Cambridge: Cambridge University Press; 2008:534.

22. Miguel I, González-Martínez SC, Alia R, Gil LA: Growth phenology and mating system of maritime pine (*Pinus pinaster* Aiton) in central Spain. *Inv Agrar Sist y Rec For* 2002, **11**:193-204.

23. Castro J, Zamora R, Hodar JA, Gomez JM: Seedling establishment of a boreal tree species (*Pinus sylvestris*) at its southernmost distribution limit:

consequences of being in a marginal Mediterranean habitat. *J Ecol* 2004, **92**:266-277.

24. Petit RJ, Hampe A: Some evolutionary consequences of being a tree. *Ann Rev Evol Ecol Syst* 2006, **37**:187-214.

25. González-Martínez SC, Gerber S, Cervera MT, Martínez-Zapater JM, Alia R, Gil LA: Selfing and sibship structure in a two-cohort stand of maritime pine (*Pinus pinaster* Ait.) using nuclear SSR markers. *Ann For Sci* 2003, **60**:115-121.

26. González-Martínez SC, Burczyk J, Nathan R, Nanos N, Gil LA, Alia R: Effective gene dispersal and female reproductive success in Mediterranean maritime pine (*Pinus pinaster* Aiton). *Mol Ecol* 2006, **15**:4577-4588.

27. Lascoux DM, Kremer A, Dormling I: Growth and phenology of 1-year-old maritime pine (*Pinus pinaster*) seedlings under continuous light - implications for early selection. *Can J For Res* 1993, **23**:1325-1336.

28. Jonsson A, Eriksson G, Ye ZH, Yeh FC: A retrospective early test of *Pinus sylvestris* seedlings grown at wide and dense spacing. *Can J For Res* 2000, **30**:1443-1452.

29. Notivol E, Garcia-Gil MR, Alia R, Savolainen O: Genetic variation of growth rhythm traits in the limits of a latitudinal cline in Scots pine. *Can J For Res* 2007, **37**:540-551.

30. De-Lucas AI, Robledo-Arnuncio JJ, Hidalgo E, González-Martínez SC: Mating system and pollen gene flow in Mediterranean maritime pine. *Heredity* 2008, **100**:390-399.

31. Meagher TR, Thompson EA: The relationship between single parent and parent pair genetic likelihoods in genealogy reconstruction. *Theor Popul Biol* 1986, **29**:87-106.

32. Gerber S, Mariette S, Streiff R, Bodénès C, Kremer A: Comparison of microsatellites and amplified fragment length polymorphism markers for parentage analysis. *Mol Ecol* 2000, **9**:1037-1048.

33. González-Martínez SC, Gerber S, Cervera MT, Martínez-Zapater J, Gil LA, Alia R, Martínez-Zapater M: Seed gene flow and fine-scale structure in a Mediterranean pine (*Pinus pinaster* Ait.) using nuclear microsatellite markers. *Theor Appl Genet* 2002, **104**:1290-1297.

34. Devlin B, Roeder K, Ellstrand NC: Fractional paternity assignment: theoretical development and comparison to other methods. *Theor Appl Genet* 1988, **76**:369-380.

35. Climent J, Chambel MR, Pardos M, Lario F, Villar-Salvador P: Biomass allocation and foliage heteroblasty in hard pine species respond differentially to reduction in rooting volume. *Eur J For Res* 2011, **130**:841-850.

36. Poorter H, Nagel O: The role of biomass allocation in the growth response of plants to different levels of light, CO₂, nutrients and water: a quantitative review. *Aust J Plant Physiol* 2000, **27**:595-607.

37. Cregg BM, Zhang JW: Physiology and morphology of *Pinus sylvestris* seedlings from diverse sources under cyclic drought stress. *For Ecol Manage* 2001, **154**:131-139.

38. Nguyen-Queyrens A, Ferhi A, Loustau D, Guehl J-M: Within-ring δ₁₃C spatial variability and interannual variations in wood cellulose of two contrasting provenances of *Pinus pinaster*. *Can J For Res* 1998, **28**:766-773.

39. Fernández M, Gil LA, Pardos JA: Response of *Pinus pinaster* Ait. provenances at early age to water supply. I. Water relation parameters. *Ann For Sci* 1999, **56**:179-187.

40. Climent J, Silva FCE, Chambel MR, Pardos M, Almeida MH: Freezing injury in primary and secondary needles of Mediterranean pine species of contrasting ecological niches. *Ann For Sci* 2009, **66**:407-414.

41. Harfouche A: Retrospective early test for adult vigor of *Pinus pinaster* families grown under two water regimes. Implications for early selection. *Ann For Sci* 2003, **60**:539-547.

42. Chambel MR, Climent J, Alia R: Divergence among species and populations of Mediterranean pines in biomass allocation of seedlings grown under two watering regimes. *Ann For Sci* 2007, **64**:87-97.

43. Aranda I, Alia R, Ortega U, Dantas ÁK, Majada J: Intra-specific variability in biomass partitioning and carbon isotopic discrimination under moderate drought stress in seedlings from four *Pinus pinaster* populations. *Tree Genet Genomes* 2010, **6**:169-178.

44. Sanchez-Gomez D, Majada J, Alia R, Feito I, Aranda I: Intraspecific variation in growth and allocation patterns in seedlings of *Pinus pinaster* Ait. submitted to contrasting watering regimes: can water availability explain regional variation? *Ann For Sci* 2010, **67**:505-513.

45. Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release 4.0. Hemel Hempstead, HP1. 1ES. UK: VSN International Ltd; 2012:374. January.

46. Lynch M, Walsh B: *Genetics and Analysis of Quantitative Traits*. Massachusetts: Sinauer Assoc; 1998:980.

<footnote n="1">¹Department of Forest Ecology and Genetics, INIA-Forest Research Centre (CIFOR), Avda. A Coruña km 7.5, 28040, Madrid, Spain. ²Sustainable Forest Management Research Institute, University of Valladolid-INIA, Madrid, Spain.</footnote>

<footnote n="3">³Unit of Forest Resources, CITA. Carr. Montañana s/n, 50059 Zaragoza, Spain.</footnote>

<footnote n="4">⁴Department of Ecology and Evolution, University of Lausanne, 1015, Lausanne, Switzerland.</footnote>
47. Fahrmeir L, Tutz G: Multivariate Statistical Modelling Based on Generalized Linear Models. New York: Springer-Verlag; 1994:548.

48. Costa e Silva J, Kerr RJ: Accounting for competition in genetic analysis, with particular emphasis on forest genetic trials. *Tree Genet Genomes* 2013, **9**:1–17.

49. Grant PR, Grant BR: Predicting microevolutionary responses to directional selection on heritable variation. *Evolution* 1995, **49**:241–251.

50. Climent J, Chambel MR, Lopez R, Mutke S, Alia R, Gil LA: Population divergence for heteroblasty in the Canary Island pine (*Pinus canariensis*, *Pinaceae*). *Am J Bot* 2006, **93**:840–848.

51. Benito-Garzón M, Alia R, Robson TM, Zavala MA: Intra-specific variability and plasticity influence potential tree species. *Glob Ecol Biogeogr* 2011, **20**:766–788.

52. Chevin L-M, Lande R, Mace GM: Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory. *PLoS Biol* 2010, **8**:e1000357.

53. Valladares F, Gianoli E, Gomez JM: Ecological limits to plant phenotypic plasticity. *New Phytol* 2007, **176**:749–763.

54. Chevin L-M, Collins S, Lefèvre F: Phenotypic plasticity and evolutionary demographic responses to climate change: taking theory out to the field. *Funct Ecol* 2013, **27**:967–979.

55. Sánchez-Gómez D, Valladares F, Zavala MA: Performance of seedlings of Mediterranean woody species under experimental gradients of irradiance and water availability: trade-offs and evidence for niche differentiation. *New Phytol* 2006, **170**:795–806.

56. Gaspar MJ, Velasco T, Feito I, Alia R, Majada J: Genetic variation of drought tolerance in *Pinus pinaster* at three hierarchical levels: a comparison of induced osmotic stress and field testing. *PLoS One* 2013, **8**:e79094.

57. Obeso JR: A hierarchical perspective in allocation to reproduction from whole plant to fruit and seed level. *Perspect Plant Ecol Evol Syst* 2004, **6**:217–225.

58. Bateman AJ: Intrasexual selection in Drosophila. *Heredity* 1948, **2**:349–368.

59. Durel CE, Bertin P, Kremer A: Relationship between inbreeding depression and inbreeding coefficient in maritime pine (*Pinus pinaster*). *Theor Appl Genet* 1996, **92**:347–356.

60. De-Lucas AI, González-Martínez SC, Vendramin GG, Hidalgo E, Heuertz M: Spatial genetic structure in continuous and fragmented populations of *Pinus pinaster* Aiton. *Mol Ecol* 2009, **18**:4564–4576.

61. Reich PB, Ellsworth DS, Walters MB: Leaf structure (specific leaf area) modulates photosynthesis-nitrogen relations: evidence from within and across species and functional groups. *Funct Ecol* 1998, **12**:948–958.

62. Steinger T, Roy BA, Stanton ML: Evolution in stressful environments II: adaptive value and costs of plasticity in response to low light in *Sinapis arvensis*. *J Evol Biol* 2003, **16**:313–323.

63. Galmes J, Cifre J, Medrano H, Flexas J, Galmés J: Modulation of relative growth rate and its components by water stress in Mediterranean species with different growth forms. *Oecologia* 2005, **145**:21–31.

64. Agrawal AA, Erwin AC, Cook SC: Natural selection on and predicted responses of ecophysiological traits of swamp milkweed (*Asclepias incarnata*). *J Ecol* 2008, **96**:536–542.

65. Schmid B, Dolt C: Effects of maternal and paternal environment and genotype on offspring phenotype in *Solidago altissima* L. *Evolution* 1994, **48**:1525–1549.

66. Sgrò CM, Hoffmann AA: Genetic correlations, tradeoffs and environmental variation. *Heredity* 2004, **93**:241–248.

67. Bell G: Fluctuating selection: the perpetual renewal of adaptation in variable environments. *Phil Trans R Soc Lond B* 2010, **365**:87–97.

68. Nathan R, Ne'eman G: Spatiotemporal dynamics of recruitment in Aleppo pine (*Pinus halepensis* Miller). *Plant Ecol* 2004, **171**:123–137.

doi:10.1186/s12862-014-0200-5

Cite this article as: Alía et al.: Environment-dependent microevolution in a Mediterranean pine (*Pinus pinaster* Aiton). BMC Evolutionary Biology 2014 **14**:200.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central