# Fifty years of genetic studies: what to make of the large amounts of variation found within populations?

Ivan Scotti, Santiago C. González-Martínez, Katharina B. Budde, Hadrien Lalagüe

► To cite this version:

Ivan Scotti, Santiago C. González-Martínez, Katharina B. Budde, Hadrien Lalagüe. Fifty years of genetic studies: what to make of the large amounts of variation found within populations?. Annals of Forest Science, 2016, 73 (1), pp.69-75. 10.1007/s13595-015-0471-z. hal-01290976

HAL Id: hal-01290976

https://hal.science/hal-01290976

Submitted on 19 Mar 2016

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>
# Fifty years of genetic studies: what to make of the large amounts of variation found within populations?

Ivan Scotti · Santiago C. González-Martínez · Katharina B. Budde · Hadrien Lalagüe

Received: 21 November 2014 / Accepted: 3 March 2015 / Published online: 20 March 2015
© INRA and Springer-Verlag France 2015

**Keywords** Microgeography · Natural selection · Population genomics · Intra-specific variation

## 1 Introduction

Determining the sources and role of intra-specific genetic variation is a classical focus of evolutionary biology (Mitchell-Olds et al. 2007). Ever since the beginning of population genetic studies in forest trees, the observation of high levels of within-stand phenotypic and, later, molecular diversity has been a commonplace. In these sessile and long-lived organisms, the maintenance of adaptive genetic and phenotypic diversity within populations seems of paramount importance, because the environment is likely to change within their life span (Petit and Hampe 2006).

**Handling Editor:** Jean-Michel Leban

**Contribution of the co-authors** All authors contributed to the writing of the manuscript.

I. Scotti (✉)
INRA, UR629 Ecologie des Forêts Méditerranéennes (URFM), Domaine St Paul, Site Agroparc, 84914 Avignon, France
e-mail: ivan.scotti@paca.inra.fr

S. C. González-Martínez · K. B. Budde
INIA, Forest Research Centre, Carretera de La Coruña km 7.5, 28040 Madrid, Spain

K. B. Budde
Department of Geosciences and Natural Resource Management, University of Copenhagen, Rolighedsvej 23, 1958 Frederiksberg C, Denmark

H. Lalagüe
INRA, UMR745 Ecologie des Forêts de Guyane, Campus Agronomique, 97387 Kourou, French Guiana

Phenotype diversity has been recurrently reported in forest trees (Borghetti et al. 1988; Cornelius 1994; Aitken et al. 1996; Howe et al. 2003; Wightman et al. 2008; O'Reilly-Wapstra et al. 2013); variability in life history traits (in particular, phenology) and adaptation to stress have consistently been a major subject in forestry. Such studies have typically focussed on variation at the regional or species range level, in relation to large-scale environmental (mostly climatic) gradients (Savolainen et al. 2007; Alberto et al. 2013). Examples abound in several tree species: to quote only a few, we may mention studies on oak bud burst (Ducousso et al. 1996), poplar leaf phenology (Hall et al. 2007), spruce bud burst and bud set (Mimura and Aitken 2010), Scots pine bud burst and frost resistance (Hurme et al. 2000) and growth performance in Spanish cedar (Navarro et al. 2004). The question of the adaptive relevance of clinal genetic variation was explicitly addressed in most such studies or underlay their experimental setup. The amount of genetically based trait variance due to differentiation between populations, as measured by $Q_{ST}$ (Spitze 1993; Prout and Barker 1993), varies from moderate to high depending on the trait; nevertheless, a substantial portion of variance is explained by family differences within populations, as $Q_{ST}$ is generally far smaller than one (see review for forest trees in Latta 2004). At the molecular level, intra-population variance for nuclear loci generally accounts for around 95 % of total variance, as observed by Hamrick and Godt (1990) and largely confirmed by all subsequent studies (see review in Petit et al. 2005). Early molecular marker-based studies reported, in addition to large amounts of intra-population diversity, significant fine-scale spatial genetic structure (FSGS) in several species; the first such study was performed by Brunel and Rodolphe (1985) in *Picea abies*, followed by a steady stream of publications ever since. Processes leading to the establishment of FSGS have been studied in detail in trees, and modelling of such processes has often been carried out by tree population geneticists
(Hardy and Vekemans 1999; Smouse et al. 2001; Sagnard et al. 2007). Such detailed knowledge of FSGS has mainly insisted on neutral processes: it was implicitly assumed that genetic variation was a *marker* for tracing individual relatedness, not for adaptive processes, and the question of whether within-stand genetic variation may have any adaptive meaning was only seldom asked (but see, e.g. Epperson 1992). This was reasonable, given the nature of the markers used (almost always putatively neutral), but ignored growing evidence of quantitative trait aggregation within populations (e.g. Nanos et al. 2004).

## 2 The adaptive meaning of local variability

FSGS studies take explicitly into account the spatial distribution of genotypes. They are the (neutral) foundation of the study of “microgeographic” processes, which focuses on evolutionary mechanism occurring within spatially structured populations (Richardson et al. 2014). The “microgeography” keyword recurrently appeared in the plant population literature, starting in the 1980s. In its broader meaning, it refers to studies carried out at geographical scales compatible with substantial amounts of gene dispersal; this spatial scale generally corresponds to the occurrence of continuous populations.

On a landscape scale, two strategies have recently been developed to relate genetic and genomic patterns of diversity to landscape features. On the one hand, environmental features and the related neutral and adaptive processes that shape genetic divergence patterns were screened using multivariate methods (e.g. reviewed in Balkenhol et al. 2009; Wang and Bradburd 2014), while on the other hand, researchers aimed to identify phenotypes and loci under selection and environmental factors driving local adaptation (Sork et al. 2013). Recent landscape genetic studies came closer to analysing the causes of microgeographic genetic variation; nevertheless, the conceptual frame of reference of landscape genomics is to compare distinct populations sampled across a gradient or in contrasting environments, rather than looking at within-stand variation, and thus cannot help explain intra-population patterns. Yet, plant population biologists have known for long that the distribution of plant species, phenotypes and genotypes can be shaped by stark environmental contrasts even at very short spatial scales: Janis Antonovics et al. have published a long series of articles both from the theoretical (Dickinson and Antonovics 1973) and the empirical point of view on heavy metal tolerance in *Anthoxanthum odoratum* at local scales (Antonovics 2006 and all the previous papers of the “Evolution in closely adjacent plant populations” series quoted therein); oat (*Avena barbata*) was also the target of similar studies, based on isozyme data (Hamrick and Allard 1972); the long-term Park Grass Experiment has shown microgeographic adaptation in-the-

making (Gould et al. 2014); the “Evolutionary Canyon” was the subject of several studies on microsite adaptation to moisture gradients at the molecular level, mostly on cereal plants (Li et al. 2000 and citations therein); more recently, Delph and Kelly (2014) provided an interpretation of plant intra-population genetic diversity as clearly adaptive; and Pannell and Fields (2014) explicitly highlighted local scale processes in their recent review on local adaptation.

## 3 Microgeographic variation: a neglected field in forest science

Contrary to annuals and herbs, the study of microgeographic adaptation has drawn only little attention from forest geneticists. The reasons are unclear: an early focus on provenance variation may have contributed to the generalised thought that the genetic variation among populations is the one that matters (see, for example, Morgenstern 1996). Furthermore, the long life span of trees, during which not only broad climate patterns but also local environmental conditions can change, makes it difficult to identify the specific critical moments and conditions which shape within-stand genetic structure. Finally, perhaps estimates of dispersal distances convinced the forest genetics community that selection would be systematically overridden by gene flow at spatial scales smaller than the landscape, and the advent of neutral theories (Kimura 1985; Hubbell 2001) pushed forest scientists to favour a view in which selection is rarely observed. Even though microgeographic adaptation studies in trees have been historically rare, forest science has consistently (and sensibly) supported the idea that maintaining high within-population genetic diversity (or large population sizes, which is a prerequisite for intra-population diversity) is a major target in the conservation and management of tree populations (e.g. Lefèvre et al. 2013). Populations must maintain sufficient levels of diversity to ensure survival, but the understanding of the origin and mechanisms underlying the maintenance of such large amounts of diversity within forest stands is limited. Moreover, the “option value” of genetic diversity (i.e. the idea that given that we do not know which combinations of genetic variants will be favoured by future selection, the best strategy is to preserve as much diversity as possible; Koskela et al. 2013), while intuitive and strongly supported by quantitative genetics theory, has not been rigorously demonstrated. Proving experimentally that some level of genetic diversity increases fitness in forest trees is elusive, because one should be able to tear apart population-level and individual effects; to prove such effects, one should manipulate diversity and perform long-term monitoring of the evolution of fitness in multiple, distinct environments (Hughes et al. 2008), which can be extremely difficult with long-lived forest trees. At the heart of the study of population-level genetic diversity in trees, there
seems to lie a conundrum: we often explain its existence with purely neutral considerations, yet the “genetic diversity option value” tenet links long-term maintenance of forest stands (that is, their fitness) to their levels of genetic diversity. How can this be? Add to this that the heritable basis of several fitness-related traits has been known for decades, and the knowledge gap becomes self-evident: we know that there is ample within-stand diversity, we know that some substantial part of phenotypic diversity is heritable and we assume that diversity is the fuel of adaptation, but we do not know how to relate such variation to current or past adaptive processes in natural stands. So, how should this large genetic diversity be interpreted and managed? Has it got any adaptive meaning? How to prevent losing the portion that “matters”, if any? Keeping overall high levels of genetic diversity within populations, in accordance with option value arguments (see Jump et al. 2009), is sensible conservation-wise but suffers from two possible flaws: it does not provide any guidance in case we need to make “no regret” choices (a highly probable scenario given limited resources), and it assumes that all we know about past and current selection is useless to predict future evolution. To overcome these limitations, we need to understand to what extent standing genetic diversity is maintained by adaptive processes and how to develop targeted measures to preserve, or engineer, adaptation in a precise manner and over short time spans. Besides the fundamental endeavour of acknowledging the role of selection in shaping local variability patterns, such questions have also an applied side: managing diversity by “keeping as much as possible” is a blunt tool; if we have to act here and now to preserve forests in the context of rapid global change, we better hold a scalpel, to complement general genetic diversity preservation strategies with specific, more refined actions.

## 4 Microgeographic variation in forest trees: state of the art

What do we know about the adaptive significance of within-stand genetic diversity? How much of it is shaped by environmental variation and adaptation? First, let us look at the larger picture: are tree species distributed according to any environmental parameter varying at the geographical scale of individual trees? There are some reasons to think this is the case, at least in landscapes displaying microgeographic variation in parameters affecting tree fitness. John et al. (2007) have shown that soil properties can explain a non-negligible portion of the spatial distribution of tree species in tropical forests in Panama, Colombia and Ecuador. At higher latitudes, soil depth and clay content in typically sandy soils determines the presence of either *Pinus pinaster* or *Pinus pinea* in continuous mixed conifer forests in central Spain (Allué and Allué 1994). These examples illustrate that neutral processes

are not sufficient to explain species distributions, even at small spatial scales. It is not a big conceptual leap to imagine that the same holds true for genotypes within populations of a single species, given that within-population phenotypic variance can be as large as, or larger than, interspecific variance (Bangert et al. 2005; Shuster et al. 2006; Brousseau et al. 2013). Proofs that microgeographic environmental gradients and contrasts influence genotype and phenotype distribution are rare but consistent: examples are known from a variety of species and habitats. Table 1 summarises the microgeographic studies in forest trees described below. In Engelmann spruce (*Picea engelmannii*), Mitton et al. (1989) found association between genotypes at the phosphoglucomutase (PGM) isozyme locus and drought stress in sub-populations separated by few tens of metres; the same research group found associations between isozyme genotypes at the glycerate dehydrogenase (GLY) locus and stomatal size, on the one side, and soil moisture, on the other side, in pinyon pine (*Pinus edulis*) over distances in the order of hundreds of metres (Cobb et al. 1994; Mitton et al. 1998); Kelly et al. (2003) found that a single birch (*Betula pendula*) stand could be sorted into genetic groups associated to temperature at seedling establishment, according to 358 amplified fragment length polymorphism (AFLP) markers. In a similar study, in European beech (*Fagus sylvatica*), Jump et al. (2006) found one AFLP marker (out of 254 screened) associated with spatial and temporal variations in temperature at seedling establishment and differentiation between sub-populations less than 3 km apart; in the same species, Pluess and Weber (2012) found three AFLPs (among 517) associated with moisture across populations separated by only 0.5 km, and Csilléry et al. (2014) identified contrasted signatures of selection for multi-locus single nucleotide polymorphism (SNP) combinations in a set of 546 SNPs from 53 climate response genes, in sub-populations sampled over an elevation gradient spanning less than 5 km. A similar study in Norway spruce (Scalfi et al. 2014) found two SNP loci (out of 384) associated with elevation; Brousseau et al. (2013) found extensive divergence for functional traits in contiguous sub-populations of tropical *Eperua falcata* and *Eperua grandiflora*, interspersed and located on seasonally flooded or vertically drained soils; in the same populations of *E. falcata*, Audigeos et al. (2013) found that genotypes and alleles were associated to soil type at eight out of 88 SNPs from seven candidate genes and for three of the same genes at the haplotype level. In a further study, six AFLP markers out of about 1,200 were found to be divergent between sub-populations of the same species, submitted to the same environmental contrasts (Brousseau et al. 2015). In *P. pinaster*, Vizcaíno-Palomar et al. (2014) found higher fitness (as evaluated by survival in the field) for seedlings sampled from a foreign population growing in harsher climate than for those from the local population but only in particular microenvironments (Vizcaíno-Palomar et al. 2014); in this same species,

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

<figure><img src="image_3.png" /></figure>
<table><caption>**Table 1** Summary of the studies about microgeographic genetic variation in forest trees</caption><thead><tr><th>Species</th><th>Genetic marker</th><th>Phenotype</th><th>Environmental factor</th><th>Scale (km)</th><th>Citation</th></tr></thead><tbody><tr><td>Eperua falcata</td><td>—</td><td>Functional traits</td><td>Soil moisture</td><td>&lt;0.5</td><td>Brousseau et al. (2013)</td></tr><tr><td>Eperua grandiflora</td><td>—</td><td>Functional traits</td><td>Soil moisture</td><td>&lt;0.5</td><td>Brousseau et al. (2013)</td></tr><tr><td>Pinus edulis</td><td>Isozymes</td><td>Functional traits</td><td>Temporal soil moisture variation</td><td>&lt;0.1</td><td>Cobb et al. (1994)</td></tr><tr><td>Pinus edulis</td><td>Isozymes</td><td>Stomatal size</td><td>Soil moisture</td><td>&lt;0.6</td><td>Mitton et al. (1998)</td></tr><tr><td>Picea engelmannii</td><td>Isozymes</td><td>—</td><td>Drought stress</td><td>&lt;0.1</td><td>Mitton et al. (1989)</td></tr><tr><td>Betula pendula</td><td>AFLPs</td><td>—</td><td>Temperature at seedling establishment</td><td>&lt;0.2</td><td>Kelly et al. (2003)</td></tr><tr><td>Fagus sylvatica</td><td>AFLPs</td><td>—</td><td>Elevational gradient and temperature at seedling establishment</td><td>&lt;3</td><td>Jump et al. (2006)</td></tr><tr><td>Fagus sylvatica</td><td>AFLPs</td><td>—</td><td>Moisture</td><td>&lt;0.5</td><td>Pluess and Weber (2012)</td></tr><tr><td>Eperua falcata</td><td>SNPs</td><td>—</td><td>Soil moisture</td><td>&lt;0.4</td><td>Audigeos et al. (2013)</td></tr><tr><td>Eperua falcata</td><td>AFLPs</td><td>—</td><td>Soil moisture</td><td>&lt;0.4</td><td>Brousseau et al. (2015)</td></tr><tr><td>Pinus pinaster</td><td>SNPs</td><td>—</td><td>Elevational gradient</td><td>&lt;1.8</td><td>Budde (2014)</td></tr><tr><td>Picea abies</td><td>SNPs</td><td>—</td><td>Elevational gradient</td><td>&lt;2.6</td><td>Scalfi et al. (2014)</td></tr><tr><td>Fagus sylvatica</td><td>SNPs</td><td>—</td><td>Elevational gradient</td><td>&lt;5</td><td>Csilléry et al. (2014)</td></tr></tbody></table>

Budde (2014) found association between individual genetic distance and small elevation differences in a continuous forest, independent of geographic distance. These studies indicate that microgeographic variation contributes to shape within-stand genetic structure. Most of them explicitly showed that local adaptation occurs and can be detected, provided that the necessary trait, genetic and environmental information are gathered. The approach shared by all these studies is straightforward: identify an environmental gradient or contrast and look for genotype- (or trait-) environment associations. Note that these gradients, although identified by a single variable, probably represent bundles of correlated variables (e.g. soil moisture is associated with soil fertility and concentration of particular compounds, and elevation is associated with water availability and soil structure).

The analysis of microgeographic adaptation can also be undertaken with trait heritability as a starting point, focussing on intra-population genetic variation for fitness-related traits. In this case, the rationale is that variance in such traits may indicate different coexisting optima—a sign of spatially or temporally variable selection and adaptation. Some pioneering studies have shown that natural forest populations studied in situ display heritability for fitness-related traits, such as chemical traits associated to herbivore resistance in eucalypts (O'Reilly-Wapstra et al. 2013) and white spruce (Mageroy et al. 2014), serotiny and growth traits in Mediterranean pines (Budde et al. 2014; Alía et al. 2014) and leaf functional traits in cork oak (Ramírez-Valiente et al. 2014). More results of the same kind are available in beech (Bontemps et al., in preparation). These studies show that there is a sizeable amount of additive variance in adaptive traits expressed in natural populations, which can be modelled by natural selection.

## 5 How much selection?

Is natural selection on adaptive traits enough to explain the divergence patterns described above? Or is it reinforced by life history or neutral processes? And what portion of within-population genetic diversity is typically maintained by natural selection? Phenological divergence by differential response to environmental cues can promote restriction of gene flow (Soularue and Kremer 2012), which in turn can promote adaptive divergence, which could reinforce isolation (isolation-by-adaptation; Nosil et al. 2008), and limitations to dispersal of both pollen and seeds may contribute to reinforce directional selection processes that would be otherwise too weak, such mechanisms would not only induce divergence at several (combinations of) adaptive loci but also promote divergence at the genome level, mimicking widespread selection. It is therefore very likely that most polymorphisms found in a population are effectively neutral and that most within-population divergence patterns are maintained by neutral processes (Ohsawa and Ide 2008). One of the biggest empirical (and methodological) challenges here is to distinguish genuine adaptive processes from neutral processes and from neutral processes enhancing or generalising (to the genome or to other traits) the effect of adaptive processes. Another limitation is that the intensity and direction of selection are likely to vary in time, both because different life stages can have different fitness optima and because environmental conditions change over the life span of an individual tree. These factors can seriously limit our ability to detect selection; however, if we are able to setup studies taking into account temporal data (e.g. Jump et al. 2006) and multiple life stages (e.g. Alía et al. 2014), the breakdown of effects of selection into

<figure><img src="image_4.png" /></figure>
temporal and ontogenetic components should be possible. Note that forest trees, with their clear-cut life stages and long life cycle, are rather well-suited for such an approach.

The studies cited above prove that genetic diversity at the allelic, genotypic, multilocus or trait level is associated to environmental heterogeneity (Yeaman and Jarvis 2006). Detailed knowledge on how the environment varies in space and time is often lacking at fine spatial scales, and we may miss major relationships between genetic and environmental heterogeneity just because we have characterised the former but not the latter. The current and historical focus on large-scale patterns may actually be the result of constraints in the availability (and accessibility) of environmental data at the local scale. The amount of genetic variants found to be linked to microgeographic environmental variation is expected to increase as we acquire better data on environmental variables, i.e. microgeographic environmental data currently act as a “limiting factor”. A similar argument holds for genetic studies typically based on relatively few traits (<5) or genetic markers (<1,000). New phenotyping capacities (but see below) and the astonishing increase of genomic resources made recently available for forest trees will contribute to the identification of additional links between genetic divergence and environment. This will nevertheless fulfil only part of the task: as already stated above, part of the environment–genotype and environment–trait relationships that will be detected will occur by chance or due to neutral processes. We will need to finely assess false discovery rates and appropriately take demographic (e.g. dispersal) parameters into account. It is likely that the last word on adaptive processes will come from a combination of methods, necessarily including experimental approaches and/or functional validation steps which are often neglected in forestry field.

We do not know how common intra-population divergence
may be, and it is likely that underpowered studies having
sought environment-related genetic patterns, and having
found none, have never been published. Yet compared
to the common-sense intuition that adaptive patterns
should be really hard to detect at such short scales,
the few studies quoted above have found such patterns
with relatively simple off-the-shelf methods and in a
multiplicity of forest environments, from temperate to
tropical. This may indicate that, if studied with suffi-
cient power, more case studies may show habitat-
related genotype and phenotype distributions, suggesting
widespread microgeographic adaptation. Nevertheless,
one should keep in mind that studies with very large
sample numbers, genotypes, traits and environmental data
are unlikely to arise soon, and therefore statistical power will
be limited. Experimental design will have to take this into
account, with special attention to the problem of multiple test-
ing of potentially non-independent effects.

**6 Perspectives**

Accurately characterising microgeographic adaptive process-
es in such large, long-lived organisms as trees will not be an
easy task, but if we are to explore the adaptive potential to
future environmental change that resides in standing genetic
variation, this goal is inescapable, because 95 % of a tree's
genetic variability resides within populations. The first goal
will be to assess how common and how strong intra-
population selection is. Next, we will need to turn this knowl-
edge into policies, if we are to manage adaptive diversity; this
requires the intermediate step of evaluating how easily we can
drive adaptive processes. Several advances may help us ac-
quiring a better understanding of microgeographic adaptive
variation. One is obviously the widespread implementation
of high-throughput genotyping methods (currently restricted
to tree species with large commercial value), which will per-
mit the fine screening of genomic variants. As it is likely that
selection may target multilocus combinations rather than indi-
vidual polymorphisms (Le Corre and Kremer 2012; Csilléry
et al. 2014), methods to detect selection taking multilocus
combinations into account should be sought. More efficient
methods to characterise functional traits in situ will also be
required to identify the major environmental axes along which
genetic differentiation occurs, and there is a great need for the
development of phenotyping methods that are able to handle
large number of individuals and traits; some options were
reviewed by Houle et al. (2010), but such methods are still
restricted to model species. To make this task more feasible,
broadly defined groups of “candidate” traits may be identified.
An educated guess is obviously to start from those traits that
can respond strongly and rapidly to environmental variables
changing over microgeographic distances (e.g. soil
hydromorphy, composition and fertility and biotic interac-
tions, particularly with microorganisms). To quote only a cou-
ple of traits, leaf nitrogen content and specific area covary
with soil composition and fertility (Ordoñez et al. 2009;
Orwin et al. 2010). Furthermore, the identification and char-
acterisation of the environmental factors underlying genotypic
and phenotypic variation is crucial. Multivariate methods
(such as redundancy analysis (RDA), partial Mantel tests or
multiple regression of distance matrices (MRDM), reviewed
in Balkenhol et al. 2009) have been mostly used to detect
relationships between genetic structure and environmental
factors at the landscape scale (Gram and Sork 2001; Eckert
et al. 2010; Manel et al. 2012; Orsini et al. 2013), but they
have already proven useful at a microgeographic scale (e.g.
Andrew et al. 2012) and could be promising to test within tree
populations (Budde 2014). Finally, population dynamics
models integrating all aspects of dispersal and selection, in-
cluding the genetic and environmental basis of traits and their
impact on fitness (Oddou-Muratorio and Davi 2014), will be
very useful to estimate the strength of selection processes and

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

<figure><img src="image_6.png" /></figure>
to evaluate the efficiency of management policies over few generations. If we succeed in getting access to these pieces of information, and to merge them into a unified view of population processes, then we will hold a very powerful analytical and applied tool for forest evolutionary biology and ecology and for helping to manage forests in a sustainable way.

**Acknowledgments** The authors wish to thank Caroline Scotti-Saintagne for fruitful discussions and for critically reading the manuscript.

**Funding** This work was funded by the ANR-12-ADAP-0007-01 "FLAG" grant.

## References

Aitken SN, Adams WT, Schermann N, Fuchigami LH (1996) Family variation for fall cold hardiness in two Washington populations of coastal Douglas-fir (*Pseudotsuga menziesii* var. menziesii (Mirb.) Franco). For Ecol Manag 80:187-195

Alberto FJ, Aitken SN, Alía R et al (2013) Potential for evolutionary responses to climate change—evidence from tree populations. Glob Chang Biol 19:1645-1661

Alía R, Chambel R, Notivol E et al (2014) Environment-dependent microevolution in a Mediterranean pine (*Pinus pinaster* Aiton). BMC Evol Biol 14:200

Allué C, Allué M (1994) Notes on the managed development of the «Pinar Viejo» pine forest (Coca, Segovia) from 1901 to 1990. Invest Agric Sist Recur For 3:259-281

Andrew RL, Ostevik KL, Ebert DP, Rieseberg LH (2012) Adaptation with gene flow across the landscape in a dune sunflower. Mol Ecol 21:2078-2091

Antonovics J (2006) Evolution in closely adjacent plant populations: long-term persistence of prereproductive isolation at a mine boundary. Heredity 97:33-37

Audigeos D, Brousseau L, Traissac S et al (2013) Molecular divergence in tropical tree populations occupying environmental mosaics. J Evol Biol 26:529-544

Balkenhol N, Waits LP, Dezzani RJ (2009) Statistical approaches in landscape genetics: an evaluation of methods for linking landscape and genetic data. Ecography 32:818-830

Bangert RK, Turek RJ, Martinsen GD et al (2005) Benefits of conservation of plant genetic diversity to arthropod diversity. Conserv Biol 19:379-390

Borghetti M, Giannini R, Menozzi P (1988) Geographic variation in cones of Norway spruce (*Picea abies* (L.) Karst.). Silvae Genet 37:178-184

Brousseau L, Bonal D, Cigna J, Scotti I (2013) Highly local environmental variability promotes intra-population divergence of quantitative traits: an example from tropical rainforest trees. Ann Bot 112: 1169-1179

Brousseau L, Foll M, Scotti-Saintagne C, Scotti I (2015) Neutral and adaptive drivers of microgeographic genetic divergence within continuous populations: the case of the Neotropical tree *Eperua falcata* Aubl. PLoS One, (in press)

Brunel D, Rodolphe F (1985) Genetic neighbourhood structure in a population of *Picea abies* L. TAG Theor Appl Genet 71:101-110

Budde KB (2014) Genetic structure of forest trees in biodiversity hotspots at different spatial scales. Complutense University, Madrid

Budde KB, Heuertz M, Hernández-Serrano A et al (2014) In situ genetic association for serotiny, a fire-related trait, in Mediterranean maritime pine (*Pinus pinaster*). New Phytol 201:230-241

Cobb N, Mitton J, Whitham TG (1994) Genetic variation associated with chronic water and nutrient stress in Pinyon pine. Am J Bot 81: 936-940

Cornelius J (1994) Heritabilities and additive genetic coefficients of variation in forest trees. Can J For Res 24:372-379

Csilléry K, Lalagüe H, Vendramin GG et al (2014) Detecting short spatial scale local adaptation and epistatic selection in climate-related candidate genes in European beech (*Fagus sylvatica*) populations. Mol Ecol 23:4696-4708

Delph LF, Kelly JK (2014) On the importance of balancing selection in plants. New Phytol 201:45-56

Dickinson H, Antonovics J (1973) The effects of environmental heterogeneity on the genetics of finite populations. Genetics 73:713-735

Ducousso A, Guyon JP, Kremer A (1996) Latitudinal and altitudinal variation of bud burst in western populations of sessile oak (*Quercus petraea* (Matt) Liebl). Ann Sci For 53:775-782

Eckert AJ, van Heerwaarden J, Wegrzyn JL et al (2010) Patterns of population structure and environmental associations to aridity across the range of loblolly pine (*Pinus taeda* L., Pinaceae). Genetics 185: 969-982

Epperson BK (1992) Spatial structure of genetic variation within populations of forest trees. New For 6:257-278

Gould B a., McCouch S, Geber M a. (2014) Variation in soil aluminum tolerance genes is associated with rapid evolution at the Park Grass Experiment. Mol Ecol 23:5931-5933

Gram WK, Sork VL (2001) Association between environmental and genetic heterogeneity in forest tree populations. Ecology 82: 2012-2021

Hall D, Luquez V, Garcia VM et al (2007) Adaptive population differentiation in phenology across a latitudinal gradient in European aspen (*Populus tremula*, L.): a comparison of neutral markers, candidate genes and phenotypic traits. Evolution 61:2849-2860

Hamrick JL, Allard RW (1972) Microgeographical variation in allozyme frequencies in *Avena barbata*. Proc Natl Acad Sci U S A 69:2100-2104

Hamrick JL, Godt MJW (1990) Allozyme diversity in plant species. In: Brown AHD, Clegg MT, Kahler A, Weir B (eds) Plant population genetics, breeding, and genetic resources. Sinauer Associates, Sunderland, pp 43-63

Hardy OJ, Vekemans X (1999) Isolation by distance in a continuous population: reconciliation between spatial autocorrelation analysis and population genetics models. Heredity 83:145-154

Houle D, Govindaraju DR, Omholt S (2010) Phenomics: the next challenge. Nat Rev Genet 11:855-866

Howe GT, Aitken SN, Neale DB et al (2003) From genotype to phenotype: unraveling the complexities of cold adaptation in forest trees. Can J Bot 81:1247-1266

Hubbell SP (2001) The unified neutral theory of biodiversity and biogeography. Princeton University Press, Princeton

Hughes AR, Inouye BD, Johnson MTJ et al (2008) Ecological consequences of genetic diversity. Ecol Lett 11:609-623

Hurme P, Sillanpaa MJ, Arjas E et al (2000) Genetic basis of climatic adaptation in Scots pine by Bayesian quantitative trait locus analysis. Genetics 156:1309-1322

John R, Dalling JW, Harms KE et al (2007) Soil nutrients influence spatial distributions of tropical tree species. Proc Natl Acad Sci U S A 104:864-869

Jump AS, Hunt JM, Martínez-Izquierdo JA, Peñuelas J (2006) Natural selection and climate change: temperature-linked spatial and temporal trends in gene frequency in *Fagus sylvatica*. Mol Ecol 15:3469-3480

Jump AS, Marchant R, Peñuelas J (2009) Environmental change and the option value of genetic diversity. Trends Plant Sci 14:51-58

Kelly CK, Chase MW, de Bruijn A et al (2003) Temperature-based population segregation in birch. Ecol Lett 6:87-89

<figure><img src="image_7.png" /></figure>
Kimura M (1985) The neutral theory of molecular evolution. Cambridge University Press, Cambridge

Koskela J, Lefèvre F, Schueler S et al (2013) Translating conservation genetics into management: Pan-European minimum requirements for dynamic conservation units of forest tree genetic diversity. Biol Conserv 157:39–49

Latta RG (2004) Relating processes to patterns of genetic variation across landscapes. For Ecol Manag 197:91–102

Le Corre V, Kremer A (2012) The genetic differentiation at quantitative trait loci under local adaptation. Mol Ecol 21:1548–1566

Lefèvre F, Koskela J, Hubert J et al (2013) Dynamic conservation of forest genetic resources in 33 European countries. Conserv Biol 27:373–384

Li YC, Röder MS, Fahima T et al (2000) Natural selection causing microsatellite divergence in wild emmer wheat at the ecologically variable microsite at Ammiad, Israel. TAG Theor Appl Genet 100:985–999

Mageroy MH, Parent G, Germanos G et al (2014) Expression of the β-glucosidase gene Pgβglu-1 underpins natural resistance of white spruce against spruce budworm. Plant J 81:68–80

Manel S, Gugerli F, Thuiller W et al (2012) Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation. Mol Ecol 21:3729–3738

Mimura M, Aitken SN (2010) Local adaptation at the range peripheries of Sitka spruce. J Evol Biol 23:249–258

Mitchell-Olds T, Willis JH, Goldstein DB (2007) Which evolutionary processes influence natural genetic variation for phenotypic traits? Nat Rev Genet 8:845–856

Mitton J, Stutz H, Schuster W, Shea K (1989) Genotypic differentiation at PGM in Engelmann spruce from wet and dry sites. Silvae Genet 38:5–6

Mitton JB, Grant MC, Yoshino AM (1998) Variation in allozymes and stomatal size in pinyon (Pinus edulis, Pinaceae), associated with soil moisture. Am J Bot 85:1262–1265

Morgenstern EK (1996) Geographical variation in forest trees: genetic basis and application of knowledge in silviculture. UBC Press, Vancouver

Nanos N, González-Martínez SC, Bravo F (2004) Studying within-stand structure and dynamics with geostatistical and molecular marker tools. For Ecol Manag 189:223–240

Navarro C, Montagnini F, Hernández G (2004) Genetic variability of *Cedrela odorata* Linnaeus: results of early performance of provenances and families from Mesoamerica grown in association with coffee. For Ecol Manag 192:217–227

Nosil P, Egan SP, Funk DJ (2008) Heterogeneous genomic differentiation between walking-stick ecotypes: “isolation by adaptation” and multiple roles for divergent selection. Evolution 62:316–336

O'Reilly-Wapstra JM, Miller AM, Hamilton MG et al (2013) Chemical variation in a dominant tree species: population divergence, selection and genetic stability across environments. PLoS One 8:e58416

Oddou-Muratorio S, Davi H (2014) Simulating local adaptation to climate of forest trees with a physio-demo-genetics model. Evol Appl 7:453–467

Ohsawa T, Ide Y (2008) Global patterns of genetic variation in plant species along vertical and horizontal gradients on mountains. Glob Ecol Biogeogr 17:152–163

Ordoñez JC, Van Bodegom PM, Witte J-PM et al (2009) A global study of relationships between leaf traits, climate and soil measures of nutrient fertility. Glob Ecol Biogeogr 18:137–149

Orsini L, Mergeay J, Vanoverbeke J, De Meester L (2013) The role of selection in driving landscape genomic structure of the waterflea *Daphnia magna*. Mol Ecol 22:583–601

Orwin KH, Buckland SM, Johnson D et al (2010) Linkages of plant traits to soil properties and the functioning of temperate grassland. J Ecol 98:1074–1083

Pannell JR, Fields PD (2014) Evolution in subdivided plant populations: concepts, recent advances and future directions. New Phytol 201: 417–432

Petit RJ, Hampe A (2006) Some Evolutionary Consequences of Being a Tree. Annu Rev Ecol Evol Syst 37:187–214

Petit RJ, Duminil J, Fineschi S et al (2005) Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Mol Ecol 14:689–701

Pluess AR, Weber P (2012) Drought-adaptation potential in *Fagus sylvatica*: linking moisture availability with genetic diversity and dendrochronology. PLoS One 7:e33636

Prout T, Barker J (1993) F statistics in *Drosophila buzzatii*: selection, population size and inbreeding. Genetics 134:369–375

Ramírez-Valiente JA, Valladares F, Aranda I (2014) Exploring the impact of neutral evolution on intrapopulation genetic differentiation in functional traits in a long-lived plant. Tree Genet Genomes 10: 1181–1190

Richardson JL, Urban MC, Bolnick DI, Skelly DK (2014) Microgeographic adaptation and the spatial scale of evolution. Trends Ecol Evol 29: 165–176

Sagnard F, Pichot C, Dreyfus P et al (2007) Modelling seed dispersal to predict seedling recruitment: recolonization dynamics in a plantation forest. Ecol Modell 203:464–474

Savolainen O, Pyhäläri T, Knürr T (2007) Gene flow and local adaptation in trees. Annu Rev Ecol Evol Syst 38:595–619

Scalfi M, Mosca E, Di Pierro EA et al (2014) Micro- and macrogeographic scale effect on the molecular imprint of selection and adaptation in Norway spruce. PLoS One 9:e115499

Shuster SM, Lonsdorf EV, Wimp GM et al (2006) Community heritability measures the evolutionary consequences of indirect genetic effects on community structure. Evolution 60:991–1003

Smouse PE, Dyer RJ, Westfall RD, Sork VL (2001) Two-generation analysis of pollen flow across a landscape. I. Male gamete heterogeneity among females. Evolution 55:260–271

Sork VL, Aitken SN, Dyer RJ et al (2013) Putting the landscape into the genomics of trees: approaches for understanding local adaptation and population responses to changing climate. Tree Genet Genomes 9:901–911

Soularue J-P, Kremer A (2012) Assortative mating and gene flow generate clinal phenological variation in trees. BMC Evol Biol 12:79

Spitze K (1993) Population structure in *Daphnia obtusa*: quantitative genetic and allozymic variation. Genetics 135:367–374

Vizcaíno-Palomar N, Revuelta-Eugercios B, Zavala MA et al (2014) The role of population origin and microenvironment in seedling emergence and early survival in Mediterranean maritime pine (*Pinus pinaster* Aiton). PLoS One 9:e109132

Wang IJ, Bradburd GS (2014) Isolation by environment. Mol Ecol 23: 5649–5662

Wightman KE, Ward SE, Haggar JP et al (2008) Performance and genetic variation of big-leaf mahogany (*Swietenia macrophylla* King) in provenance and progeny trials in the Yucatan Peninsula of Mexico. For Ecol Manag 255:346–355

Yeaman S, Jarvis A (2006) Regional heterogeneity and gene flow maintain variance in a quantitative trait within populations of lodgepole pine. Proc Biol Sci 273:1587–1593

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

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