Adaptive gradients and isolation-by-distance with postglacial migration in Picea sitchensis

The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms

In tree species native to temperate and boreal regions, the activity-dormancy cycle is an important adaptive trait both for survival and growth. We discuss recent research on mechanisms controlling the overlapping developmental processes that define the activity-dormancy cycle, including cessation of apical growth, bud development, induction, maintenance and release of dormancy, and bud burst. The cycle involves an extensive reconfiguration of metabolism. Environmental control of the activity-dormancy cycle is based on perception of photoperiodic and temperature signals, reflecting adaptation to prevailing climatic conditions. Several molecular actors for control of growth cessation have been identified, with the CO/FT regulatory network and circadian clock having important coordinating roles in control of growth and dormancy. Other candidate regulators of bud set, dormancy and bud burst have been identified, such as dormancy-associated MADS-box factors, but their exact roles remain to be discovered. Epigenetic mechanisms also appear to factor in control of the activity-dormancy cycle. Despite evidence for gibberellins as negative regulators in growth cessation, and ABA and ethylene in bud formation, understanding of the roles that plant growth regulators play in controlling the activity-dormancy cycle is still very fragmentary. Finally, some of the challenges for further research in bud dormancy are discussed.

Predicting adaptive phenotypes from multilocus genotypes in Sitka spruce (Picea sitchensis) using random forest

Climate is the primary driver of the distribution of tree species worldwide, and the potential for adaptive evolution will be an important factor determining the response of forests to anthropogenic climate change. Although association mapping has the potential to improve our understanding of the genomic underpinnings of climatically relevant traits, the utility of adaptive polymorphisms uncovered by such studies would be greatly enhanced by the development of integrated models that account for the phenotypic effects of multiple single-nucleotide polymorphisms (SNPs) and their interactions simultaneously. We previously reported the results of association mapping in the widespread conifer Sitka spruce (Picea sitchensis). In the current study we used the recursive partitioning algorithm ‘Random Forest’ to identify optimized combinations of SNPs to predict adaptive phenotypes. After adjusting for population structure, we were able to explain 37% and 30% of the phenotypic variation, respectively, in two locally adaptive traits—autumn budset timing and cold hardiness. For each trait, the leading five SNPs captured much of the phenotypic variation. To determine the role of epistasis in shaping these phenotypes, we also used a novel approach to quantify the strength and direction of pairwise interactions between SNPs and found such interactions to be common. Our results demonstrate the power of Random Forest to identify subsets of markers that are most important to climatic adaptation, and suggest that interactions among these loci may be widespread.

Divergent selection and heterogeneous migration rates across the range of Sitka spruce (Picea sitchensis)

Gene flow and effective population size (N(e)) should depend on a population's position within its range: those near the edges are expected to have smaller N(e) and lower relative emigration rates, whereas those nearer the centre should have larger N(e) and higher relative emigration rates. In species with continuous ranges, this phenomenon may limit the ability of peripheral populations to respond to divergent selection. Here, we employ Sitka spruce as a model to test these predictions. We previously genotyped 339 single nucleotide polymorphisms (SNPs) in 410 individuals from 13 populations, and used these data to identify putative targets of divergent selection, as well as to explore the extent to which central-peripheral structure may impede adaptation. Fourteen SNPs had outlier F(ST) estimates suggestive of divergent selection, of which nine were previously associated with phenotypic variation in adaptive traits (timing of autumn budset and cold hardiness). Using coalescent simulations, we show that populations from near the centre of the range have higher effective populations sizes than those from the edges, and that central populations contribute more migrants to marginal populations than the reverse. Our results suggest that while divergent selection appears to have shaped allele frequencies among populations, asymmetrical movement of alleles from the centre to the edges of the species range may affect the adaptive capacity of peripheral populations. In southern peripheral populations, the movement of cold-adapted alleles from the north represents a significant impediment to adaptation under climate change, while in the north, movement of warm-adapted alleles from the south may enhance adaptation.

Metabolic dynamics during autumn cold acclimation within and among populations of Sitka spruce (Picea sitchensis)

• Autumnal cold acclimation in conifers is a complex process, the timing and extent of which vary widely along latitudinal gradients for many tree species and reflect local adaptation to climate. Although previous studies have detailed some aspects of the metabolic remodelling that accompanies cold acclimation in conifers, little is known about global metabolic dynamics, or how these changes vary among phenotypically divergent populations. • Using untargeted GC-MS metabolite profiling, we monitored metabolic dynamics during autumnal cold acclimation in three populations of Sitka spruce from the southern, central, and northern portions of the species range, which differ in both the timing and extent of cold acclimation. • Latitudinal variation was evident in the nature, intensity, and timing of metabolic events. Early development of strong freezing tolerance in the northern population was associated with a transient accumulation of amino acids. By late autumn, metabolic profiles were highly similar between the northern and central populations, whereas profiles for the southern population were relatively distinct. • Our results provide insight into the metabolic architecture of latitudinal adaptive variation in autumn acclimation and show that different mechanisms are the basis of early October cold hardiness and autumn-acclimated cold hardiness.

Long-distance gene flow and adaptation of forest trees to rapid climate change

Forest trees are the dominant species in many parts of the world and predicting how they might respond to climate change is a vital global concern. Trees are capable of long-distance gene flow, which can promote adaptive evolution in novel environments by increasing genetic variation for fitness. It is unclear, however, if this can compensate for maladaptive effects of gene flow and for the long-generation times of trees. We critically review data on the extent of long-distance gene flow and summarise theory that allows us to predict evolutionary responses of trees to climate change. Estimates of long-distance gene flow based both on direct observations and on genetic methods provide evidence that genes can move over spatial scales larger than habitat shifts predicted under climate change within one generation. Both theoretical and empirical data suggest that the positive effects of gene flow on adaptation may dominate in many instances. The balance of positive to negative consequences of gene flow may, however, differ for leading edge, core and rear sections of forest distributions. We propose future experimental and theoretical research that would better integrate dispersal biology with evolutionary quantitative genetics and improve predictions of tree responses to climate change.

Low genetic diversity, moderate local adaptation, and phylogeographic insights in Cornus nuttallii (Cornaceae)

PREMISE OF THE STUDY

Genetic knowledge is completely lacking for Pacific dogwood (Cornus nuttallii), a western North American temperate tree that is pollinated and dispersed by biological vectors. We investigated how history, geography, and climate have affected population genetic structure, local adaptation, and the phylogeography of this species.

METHODS

We examined patterns and levels of diversity in nuclear microsatellites (SSRs) and cpDNA haplotypes in populations from across the species range. We compared these results to population differentiation and genetic clines in phenotypic traits in a common garden.

KEY RESULTS

Genetic diversity was low for both nuclear SSRs and cpDNA. There was a lack of population structure (F(ST) = 0.090) in the coastal portion of the species range, with estimates of population genetic diversity in microsatellite markers decreasing with latitude from California to British Columbia. A disjunct interior population in Idaho 450 km from the coastal range had the lowest diversity but the highest divergence of all populations studied. Only a single nucleotide polymorphism was discovered after sequencing 5547 base pairs in seven noncoding regions of cpDNA. Both cpDNA haplotypes were widely distributed throughout the species range. Quantitative variation among populations was moderate (0.11 ≤ Q(ST) ≤ 0.63), and weak but significant adaptive clines were found between quantitative traits and population climatic variables (0.09 ≤ R(2) ≤ 0.34).

CONCLUSIONS

Cornus nuttallii likely faced a population bottleneck in a single southern refugium during the Last Glacial Maximum. Despite low genetic diversity, it is weakly to moderately locally adapted.

Environmental and climatic determinants of molecular diversity and genetic population structure in a coenagrionid damselfly

Identifying environmental factors that structure intraspecific genetic diversity is of interest for both habitat preservation and biodiversity conservation. Recent advances in statistical and geographical genetics make it possible to investigate how environmental factors affect geographic organisation and population structure of molecular genetic diversity within species. Here we present a study on a common and wide ranging insect, the blue tailed damselfly Ischnuraelegans, which has been the target of many ecological and evolutionary studies. We addressed the following questions: (i) Is the population structure affected by longitudinal or latitudinal gradients?; (ii) Do geographic boundaries limit gene flow?; (iii) Does geographic distance affect connectivity and is there a signature of past bottlenecks?; (iv) Is there evidence of a recent range expansion and (vi) what is the effect of geography and climatic factors on population structure? We found low to moderate genetic sub-structuring between populations (mean F(ST) = 0.06, D(est) = 0.12), and an effect of longitude, but not latitude, on genetic diversity. No significant effects of geographic boundaries (e.g. water bodies) were found. F(ST)-and D(est)-values increased with geographic distance; however, there was no evidence for recent bottlenecks. Finally, we did not detect any molecular signatures of range expansions or an effect of geographic suitability, although local precipitation had a strong effect on genetic differentiation. The population structure of this small insect has probably been shaped by ecological factors that are correlated with longitudinal gradients, geographic distances, and local precipitation. The relatively weak global population structure and high degree of genetic variation within populations suggest that I. elegans has high dispersal ability, which is consistent with this species being an effective and early coloniser of new habitats.

Postglacial history of a widespread conifer produces inverse clines in selective neutrality tests

Deviations of the site frequency spectrum of mutations (SFS) from neutral expectations may be caused by natural selection or by demographic processes such as population subdivision or temporal changes in population size. As most widespread temperate and boreal tree species have expanded from glacial refugia in the past 13,000 years, colonization bottlenecks associated with this migration may have left variable demographic signatures among geographic populations corresponding to distance from the refugia. To determine whether the signature of postglacial re-colonization has skewed the SFS in the widely distributed conifer Sitka spruce (Picea sitchensis (Bong.) Carr.), we re-sequenced 153 nuclear genes in six populations from across the species range. We found that while the SFS for the pooled sample produced negative values for Tajima's D and Fay and Wu's H, these statistics exhibited strong clinal variation when populations were analysed separately (R(2) = 0.84, P = 0.007 for Tajima's D and R(2)  = 0.65, P = 0.033 for Fay and Wu's H). When historical bottlenecks of varying age were simulated using approximate Bayesian computation, distance of populations from the southern range limit explained most of the variation in bottleneck timing among populations (R(2) = 0.89, P = 0.003). These data suggest that sequential population bottlenecks during postglacial re-colonization have resulted in diverse among-population signatures within the contemporary SFS in Sitka spruce, with rare variants more common in the south, and medium-frequency variants more common in the north. Our results also emphasize the need to consider sampling strategy and to explore population-specific null demographic models in surveys of nucleotide variation in widely distributed species.

Pushing north one bottleneck at a time: site frequency spectra tell the history of Sitka spruce

Reconstructing the history of populations is a longstanding goal of molecular ecologists. In addition to a better understanding of the past, it is hoped that this knowledge would also facilitate predictions regarding species' responses to future events such as climate change. The traditional way of doing this is through the fossil record, but these historical records are often incomplete. Inferring historical demography from patterns of nucleotide variability can help to fill these gaps. In this issue of Molecular Ecology, Holliday et al. (2010) glimpse into the demographic past of Sitka spruce, Picea sitchensis, an economically and ecologically important species native to northwestern United States and Canada, by examining the site frequency spectrum (SFS) of 153 loci in six populations covering the species entire range.

Widespread, ecologically relevant genetic markers developed from association mapping of climate-related traits in Sitka spruce (Picea sitchensis)

• Genecological studies in widespread tree species have revealed steep genetic clines along environmental gradients for climate-related traits. In a changing climate, the ecological and economic importance of conifers necessitates an appraisal of how molecular genetic variation shapes quantitative trait variation, and one of the most promising approaches to answer this question is association mapping. • We phenotyped a wide collection of 410 individuals of the widely distributed conifer Sitka spruce rangewide (Picea sitchensis) for budset timing and autumn cold hardiness, and genotyped these individuals for a panel of 768 single nucleotide polymorphisms (SNPs) representing > 200 expressed nuclear genes. • After correcting for population structure, associations were detected in 28 of the candidate genes, which cumulatively explained 28 and 34% of the phenotypic variance in cold hardiness and budset, respectively. Most notable among the associations were five genes putatively involved in light signal transduction, the key pathway regulating autumn growth cessation in perennials. Many SNPs with phenotypic associations were also correlated with at least one climate variable. • This study represents a significant step toward the goal of characterizing the genomic basis of adaptation to local climate in conifers, and provides an important resource for breeding and conservation genetics in a changing climate.

The genetic structure of red raspberry (Rubus idaeus L.) populations in Lithuania

The red raspberry (Rubus idaeus L.) is widely distributed in Lithuania and occupies a range of habitats. The presence of coadapted gene pools in local populations of R. idaeus is a question of interest not only to plant scientists, but also to plant breeders. In this study, we investigated the genetic structure of R. idaeus and the influence of local habitats on the genetic diversity within and among populations. Nineteen populations of R. idaeus were sampled from different habitats in various agroclimatic subregions of Lithuania, and analyzed using RAPD markers. 113 RAPD bands were identified among 315 individuals; 84.31% of these were polymorphic. The mean values of Shannon’s information index for different populations ranged from 0.341 to 0.455. Nei’s gene diversity established within populations averaged 0.266. An AMOVA revealed 74% of genetic variation among individuals within populations of R. idaeus, and 23% among populations. The remaining genetic variation was distributed among populations from different agroclimatic subregions (3%). The results of this study suggest that the genetic structure of R. idaeus populations in Lithuania may be influenced partially by isolation by distance as well as by local environmental conditions.

Local adaptation at the range peripheries of Sitka spruce

High-dispersal rates in heterogeneous environments and historical rapid range expansion can hamper local adaptation; however, we often see clinal variation in high-dispersal tree species. To understand the mechanisms of the species' distribution, we investigated local adaptation and adaptive plasticity in a range-wide context in Sitka spruce, a wind-pollinated tree species that has recently expanded its range after glaciations. Phenotypic traits were observed using growth chamber experiments that mimicked temperature and photoperiodic regimes from the limits of the species realized niche. Bud phenology exhibited parallel reaction norms among populations; however, putatively adaptive plasticity and strong divergent selection were seen in bud burst and bud set timing respectively. Natural selection appears to have favoured genotypes that maximize growth rate during available frost-free periods in each environment. We conclude that Sitka spruce has developed local adaptation and adaptive plasticity throughout its range in response to current climatic conditions despite generally high pollen flow and recent range expansion.

Adaptation, migration or extirpation: climate change outcomes for tree populations

Species distribution models predict a wholesale redistribution of trees in the next century, yet migratory responses necessary to spatially track climates far exceed maximum post-glacial rates. The extent to which populations will adapt will depend upon phenotypic variation, strength of selection, fecundity, interspecific competition, and biotic interactions. Populations of temperate and boreal trees show moderate to strong clines in phenology and growth along temperature gradients, indicating substantial local adaptation. Traits involved in local adaptation appear to be the product of small effects of many genes, and the resulting genotypic redundancy combined with high fecundity may facilitate rapid local adaptation despite high gene flow. Gene flow with preadapted alleles from warmer climates may promote adaptation and migration at the leading edge, while populations at the rear will likely face extirpation. Widespread species with large populations and high fecundity are likely to persist and adapt, but will likely suffer adaptational lag for a few generations. As all tree species will be suffering lags, interspecific competition may weaken, facilitating persistence under suboptimal conditions. Species with small populations, fragmented ranges, low fecundity, or suffering declines due to introduced insects or diseases should be candidates for facilitated migration.

Global monitoring of autumn gene expression within and among phenotypically divergent populations of Sitka spruce (Picea sitchensis)

Cold acclimation in conifers is a complex process, the timing and extent of which reflects local adaptation and varies widely along latitudinal gradients for many temperate and boreal tree species. Despite their ecological and economic importance, little is known about the global changes in gene expression that accompany autumn cold acclimation in conifers. Using three populations of Sitka spruce (Picea sitchensis) spanning the species range, and a Picea cDNA microarray with 21,840 unique elements, within- and among-population gene expression was monitored during the autumn. Microarray data were validated for selected genes using real-time PCR. Similar numbers of genes were significantly twofold upregulated (1257) and downregulated (967) between late summer and early winter. Among those upregulated were dehydrins, pathogenesis-related/antifreeze genes, carbohydrate and lipid metabolism genes, and genes involved in signal transduction and transcriptional regulation. Among-population microarray hybridizations at early and late autumn time points revealed substantial variation in the autumn transcriptome, some of which may reflect local adaptation. These results demonstrate the complexity of cold acclimation in conifers, highlight similarities and differences to cold tolerance in annual plants, and provide a solid foundation for functional and genetic studies of this important adaptive process.

Ecological genetics and seed transfer guidelines for Pinus albicaulis (Pinaceae)

Whitebark pine (Pinus albicaulis Engelm.) has greatly declined throughout its range as a result of introduced disease, fire suppression, and other factors, and climate change is predicted to accelerate this decline. Restoration is needed; however, no information regarding the degree of local adaptation is available to guide these efforts. A seedling common-garden experiment was employed to assess genetic diversity and geographic differentiation (Q(ST)) of whitebark pine for traits involved in growth and adaptation to cold and to determine climatic variables revealing local adaptation. Seedlings from 48 populations were grown for two years and measured for height increment, biomass, root to shoot ratio, date of needle flush, fall and spring cold injury, and survival. Significant variation was observed among populations for most traits. The Q(ST) was low (0.07-0.14) for growth traits and moderate (0.36-0.47) for cold adaptation related traits, but varied by region. Cold adaptation traits were strongly correlated with mean temperature of the coldest month of population origins, while growth traits were generally correlated with growing season length. We recommend that seed transfer for restoration favor seed movement from milder to colder climates to a maximum of 1.9°C in mean annual temperature in the northern portion of the species range, and 1.0°C in the U. S. Rocky Mountains to avoid maladaptation to current conditions yet facilitate adaptation to future climates.