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.

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, i.e., 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 analyze this variation across environments and years. We evaluated five categories of fitness-related traits (survival, height, phenology, functional, and biotic-stress response traits) in a clonal common-garden network, planted in contrasted environments (over 12,500 trees). Most of the analyzed traits showed evidence of local adaptation based on Q_st -F_st 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.

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.

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.

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.

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.

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

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² = 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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0200-5) contains supplementary material, which is available to authorized users.

Genetic variation of drought tolerance in Pinus pinaster at three hierarchical levels: a comparison of induced osmotic stress and field testing

Understanding the survival capacity of forest trees to periods of severe water stress could improve knowledge of the adaptive potential of different species under future climatic scenarios. In long lived organisms, like forest trees, the combination of induced osmotic stress treatments and field testing can elucidate the role of drought tolerance during the early stages of establishment, the most critical in the life of the species. We performed a Polyethylene glycol-osmotic induced stress experiment and evaluated two common garden experiments (xeric and mesic sites) to test for survival and growth of a wide range clonal collection of Maritime pine. This study demonstrates the importance of additive vs non additive effects for drought tolerance traits in Pinus pinaster, and shows differences in parameters determining the adaptive trajectories of populations and family and clones within populations. The results show that osmotic adjustment plays an important role in population variation, while biomass allocation and hydric content greatly influence survival at population level. Survival in the induced osmotic stress experiment presented significant correlations with survival in the xeric site, and height growth at the mesic site, at population level, indicating constraints of adaptation for those traits, while at the within population level no significant correlation existed. These results demonstrate that population differentiation and within population genetic variation for drought tolerance follow different patterns.

Potential for evolutionary responses to climate change - evidence from tree populations

Evolutionary responses are required for tree populations to be able to track climate change. Results of 250 years of common garden experiments show that most forest trees have evolved local adaptation, as evidenced by the adaptive differentiation of populations in quantitative traits, reflecting environmental conditions of population origins. On the basis of the patterns of quantitative variation for 19 adaptation-related traits studied in 59 tree species (mostly temperate and boreal species from the Northern hemisphere), we found that genetic differentiation between populations and clinal variation along environmental gradients were very common (respectively, 90% and 78% of cases). Thus, responding to climate change will likely require that the quantitative traits of populations again match their environments. We examine what kind of information is needed for evaluating the potential to respond, and what information is already available. We review the genetic models related to selection responses, and what is known currently about the genetic basis of the traits. We address special problems to be found at the range margins, and highlight the need for more modeling to understand specific issues at southern and northern margins. We need new common garden experiments for less known species. For extensively studied species, new experiments are needed outside the current ranges. Improving genomic information will allow better prediction of responses. Competitive and other interactions within species and interactions between species deserve more consideration. Despite the long generation times, the strong background in quantitative genetics and growing genomic resources make forest trees useful species for climate change research. The greatest adaptive response is expected when populations are large, have high genetic variability, selection is strong, and there is ecological opportunity for establishment of better adapted genotypes.

Correlated responses to clonal selection in populations of Daphnia pulicaria: mechanisms of genetic correlation and the creative power of sex

Genetic correlations among traits alter evolutionary trajectories due to indirect selection. Pleiotropy, chance linkage, and selection can all lead to genetic correlations, but have different consequences for phenotypic evolution. We sought to assess the mechanisms contributing to correlations with size at maturity in the cyclic parthenogen Daphnia pulicaria. We selected on size in each of four populations that differ in the frequency of sex, and evaluated correlated responses in a life table. Size at advanced adulthood, reproductive output, and adult growth rate clearly showed greater responses in high-sex populations, with a similar pattern in neonate size and r. This pattern is expected only when trait correlations are favored by selection and the frequency of sex favors the creation and demographic expansion of highly fit clones. Juvenile growth and age at maturity did not diverge consistently. The inter-clutch interval appeared to respond more strongly in low-sex populations, but this was not statistically significant. Our data support the hypothesis that correlated selection is the strongest driver of genetic correlations, and suggest that in organisms with both sexual and asexual reproduction, adaptation can be enhanced by recombination.

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

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

Variations in growth, survival and carbon isotope composition (delta(13)C) among Pinus pinaster populations of different geographic origins

To evaluate differences in growth and adaptability of maritime pine (Pinus pinaster Ait.), we studied growth, polycyclism, needle tissue carbon isotope composition (delta(13)C) as an estimate of water-use efficiency (WUE) and survival of seven populations at 10 years of age growing in a performance trial at a provenance test site in Escaroupim, Portugal. Six populations were from relatively high rainfall sites in Portugal and southwestern France (Atlantic group), and one population was from a more arid Mediterranean site in Spain. There were significant differences between some populations in total height, diameter at breast height, delta(13)C of bulk needle tissue, polycyclism and survival. A population from central Portugal (Leiria, on the Atlantic coast) was the tallest and had the lowest delta(13)C. Overall, the variation in delta(13)C was better explained by the mean minimum temperatures of the coldest month than by annual precipitation at the place of origin. Analyses of the relationships between delta(13)C and growth or survival revealed a distinct pattern for the Mediterranean population, with low delta(13)C (and WUE) associated with the lowest growth potential and reduced survival. There were significant negative correlations between delta(13)C and height or survival in the Atlantic group. Variation in polycyclism was correlated with annual precipitation at the place of origin. Some Atlantic populations maintained a high growth potential while experiencing moderate water stress. A detailed knowledge of the relationships between growth, survival and delta(13)C in contrasting environments will enhance our ability to select populations for forestry or conservation.

To grow or to seed: ecotypic variation in reproductive allocation and cone production by young female Aleppo pine (Pinus halepensis, Pinaceae)

Age and size at the first reproduction and the reproductive allocation of plants are linked to different life history strategies. Aleppo pine only reproduces through seed, and, as such, early female reproduction confers high fitness in its infertile highly fire-prone habitats along the Mediterranean coast because life expectancy is short. We investigated the extent of ecotypic differentiation in female reproductive allocation and examined the relation between early female reproduction and vegetative growth. In a common-garden experiment, the threshold age and size at first female reproduction and female reproductive allocation at age seven differed significantly among Aleppo pine provenances of ecologically distinct origin. Significant correlations among reproductive features of the provenances and the ecological traits of origin were found using different analytical tools. In nonlinear models of cone counts vs. stem volume, medium-sized trees (not the largest trees) produced the highest cone yield, confirming that, at the individual level, early female reproduction is incompatible with fast vegetative growth. The contribution of founder effects and adaptation to contrasting fire regimes may be confounding factors. But considering all traits analyzed, the geographical patterns of resource allocation by Aleppo pine suggest ecotypic specialization for either resource-poor (favoring early reproduction) or resource-rich (favoring vegetative growth) habitats.

Range-wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers

Some 1339 trees from 48 Pinus pinaster stands were characterized by five chloroplast microsatellites, detecting a total of 103 distinct haplotypes. Frequencies for the 16 most abundant haplotypes (p(k) > 0.01) were spatially interpolated over a lattice made by 430 grid points. Fitting of spatially interpolated values on raw haplotype frequencies at the same geographical location was tested by regression analysis. A range-wide 'diversity map' based on interpolated haplotype frequencies allowed the identification of one 'hotspot' of diversity in central and southeastern Spain, and two areas of low haplotypic diversity located in the western Iberian peninsula and Morocco. Principal component analysis (PCA) carried out on haplotypes frequency surfaces allowed the construction of a colour-based 'synthetic' map of the first three PC components, enabling the detection of the main range-scale genetic trends and the identification of three main 'gene pools' for the species: (i) a 'southeastern' gene pool, including southeastern France, Italy, Corsica, Sardinia, Pantelleria and northern Africa; (ii) an 'Atlantic' gene pool, including all the western areas of the Iberian peninsula; and (iii) a 'central' gene pool, located in southeastern Spain. Multivariate and AMOVA analyses carried out on interpolated grid point frequency values revealed the existence of eight major clusters ('gene zones'), whose genetic relationships were related with the history of the species. In addition, demographic models showed more ancient expansions in the eastern and southern ranges of maritime pine probably associated to early postglacial recolonization. The delineation of the gene zones provides a baseline for designing conservation areas in this key Mediterranean pine.

Effective gene dispersal and female reproductive success in Mediterranean maritime pine (Pinus pinaster Aiton)

Understanding population-scale processes that affect allele frequency changes across generations is a long-standing interest in genetic, ecological and evolutionary research. In particular, individual differences in female reproductive success and the spatial scale of gene flow considerably affect evolutionary change and patterns of local selection. In this study, a recently developed maximum-likelihood (ML) method based on established offspring, the Seedling Neighbourhood Model, was applied and exponentially shaped dispersal kernels were fitted to both genetic and ecological data in a widespread Mediterranean pine, Pinus pinaster Aiton. The distribution of female reproductive success in P. pinaster was very skewed (about 10% of trees mothered 50% of offspring) and significant positive female selection gradients for diameter (gamma = 0.7293) and cone crop (gamma = 0.4524) were found. The selective advantage of offspring mothered by bigger trees could be due to better-quality seeds. These seeds may show more resilience to severe summer droughts and microsite variation related to water and nutrient availability. Both approaches, ecological and of parentage, consistently showed a long-distance dispersal component in saplings that was not found in dispersal kernels based on seed shadows, highlighting the importance of Janzen-Connell effects and microenvironmental variation for survival at early stages of establishment in this Mediterranean key forest tree.

Population divergence for heteroblasty in the Canary Island pine (Pinus canariensis, Pinaceae)

A heteroblastic (or vegetative phase) change is an abrupt manifestation in the general heteroblastic development during the ontogeny of plants. The Canary Island pine undergoes an especially marked and delayed heteroblastic change, including both the formation of secondary needles on dwarf shoots and the onset of preformed growth. To assess genetic and environmental effects on the heteroblastic change in this species, we followed plants from 19 populations at a dry site and a wetter site. Comparing juvenile and adult needles from the same individuals, the adult had a significantly lower rate of water loss and higher leaf mass per area. Pooling data from all seed sources, the heteroblastic change took place when plants reached a critical height, on average, at 4 years of age at the dry site and 1 year earlier at the wet site. Within a subsample of individuals of equal size, mortality was significantly higher in juvenile plants than in mature plants. However, the juvenile phase was longer in plants from dry regions when compared to plants from highly productive, wet regions. This apparent contradiction might be explained through differential resource allocation and the cost of sclerophylly and resprouting ability. Considering the life strategy of the Canary Island pine, we interpret the prolonged juvenile phase as an unavoidable trade-off for the high tolerance of adults to harsh environments.

Low temperature during winter elicits differential responses among populations of the Mediterranean evergreen cork oak (Quercus suber)

Populations of cork oak (Quercus suber L.) were assessed for seasonal and inter-population variability in, and temperature responses of, the ratio between light-induced variable and maximum fluorescence of chlorophyll, Fv/Fm, considered a surrogate for the maximum photochemical efficiency of photosystem II (PSII). Seedlings from 10 populations throughout the distribution range of Q. suber in the Mediterranean basin were grown in a common garden in central Spain. The Fv/Fm ratio of dark-adapted leaves was measured at dawn every month for 2 years. Air temperature was recorded at a nearby climatic station. During the summer, when maximum air temperatures reached 40 degrees C, there were no significant differences in Fv/Fm among populations, but significant differences were seen during the winter. In colder months, Fv/Fm ranged in all populations between 0.5-0.6 and 0.2-0.3 in 2001 and 2002, respectively. The variance explained by the population effect was greatest during winter months, especially in 2002, reaching a peak value of 10% when minimum air temperature was below -10 degrees C. Populations originating from warmer sites showed the largest decline in Fv/Fm between the end of 2001 and the beginning of 2002. Thus, a negative linear relationship was established between mean annual temperature at the population source and population mean Fv/Fm recorded in the coldest month in 2002 and normalized by the Fv/Fm spring measurement.

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

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

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

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

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

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

Population genetic structure in a Mediterranean pine (Pinus pinaster Ait.): a comparison of allozyme markers and quantitative traits

F-statistics were employed to analyse quantitative and allozyme variation among 19 native populations of maritime pine (Pinus pinaster Ait.). Fourteen polymorphic allozyme loci were used to provide an empirical basis for constructing a null hypothesis to test natural selection as a determinant of quantitative evolution in stem form, total height growth and survival at 30 years old. Hidden biases, that may result in a difference between quantitative (Q(ST)) and allozyme (F(ST)) differentiation which are not because of the action of natural selection, were avoided by comparing pairs of populations using linear models. All quantitative traits showed higher differentiation than allozymes. The highest divergence was found in stem form, whereas divergences in total height and survival were significantly lower. Differential adaptation to regional and local patterns of precipitation, temperature and soil type seem to be the best explanation of the different structure found in quantitative traits and allozyme loci. Possible bias in the estimation of Q(ST) due to the level of quantitative within-population diversity and the role of adaptation of maritime pine after the last glaciation to highly diverse ecological conditions are discussed with special reference to the actual geographical structure of gene diversity in the species' native range.

Seed gene flow and fine-scale structure in a Mediterranean pine ( Pinus pinaster Ait.) using nuclear microsatellite markers

The Mediterranean populations of maritime pine ( Pinus pinaster Ait.) are typically small and have a scattered distribution, being threatened by human activities and forest fires. In the framework of the genetic-resources conservation program of this species, a native multi-age stand located in a Mediterranean area (central Spain) was studied using three highly polymorphic nuclear microsatellites (SSRs). Spatial autocorrelation analysis was conducted using Moran's index in order to detect fine-scale structure in both natural regeneration and mature trees. The spatial pattern of seed flow based on dispersed progeny was studied using a highly reliable subset of parent-offspring matches obtained by means of parentage analysis and simulation-based calculation of statistical confidence. Maritime pine showed a fine-scale structure at the seedling stage. In natural regeneration, the autocorrelograms indicated a patch size of approximately 10 m. The fine-scale structure seems to be produced by a restricted seed gene flow. In fact, there was an excess of parent-offspring matches in a radius of 15 m from the parent trees. Pines with a heavy seed, such as P. pinaster, are expected to have a short dispersal distance, thus producing a fine-scale structure. However, the fine-scale structure did not persist in the mature trees. Within-population genetic structure in Mediterranean pines may be affected by a number of post-dispersal events (e.g. mortality due to the severity of the Mediterranean climate and animal-mediated secondary dispersal during the summer period). Thus, great alteration in the pattern produced by the initial seed rain and differences in genetic structure between tree cohorts are expected.