Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species

Global climate is rapidly changing, and the ability for tree species to adapt is dependent on standing genomic variation; however, the distribution and abundance of functional and adaptive variants are poorly understood in natural systems. We test key hypotheses regarding the genetics of adaptive variation in a foundation tree: genomic variation is associated with climate, and genomic variation is more likely to be associated with temperature than precipitation or aridity. To test these hypotheses, we used 9,593 independent, genomic single‐nucleotide polymorphisms (SNPs) from 270 individuals sampled from Corymbia calophylla's entire distribution in south‐western Western Australia, spanning orthogonal temperature and precipitation gradients. Environmental association analyses returned 537 unique SNPs putatively adaptive to climate. We identified SNPs associated with climatic variation (i.e., temperature [458], precipitation [75] and aridity [78]) across the landscape. Of these, 78 SNPs were nonsynonymous (NS), while 26 SNPs were found within gene regulatory regions. The NS and regulatory candidate SNPs associated with temperature explained more deviance (27.35%) than precipitation (5.93%) and aridity (4.77%), suggesting that temperature provides stronger adaptive signals than precipitation. Genes associated with adaptive variants include functions important in stress responses to temperature and precipitation. Patterns of allelic turnover of NS and regulatory SNPs show small patterns of change through climate space with the exception of an aldehyde dehydrogenase gene variant with 80% allelic turnover with temperature. Together, these findings provide evidence for the presence of adaptive variation to climate in a foundation species and provide critical information to guide adaptive management practices.     

Temporal stability and assignment power of adaptively divergent genomic regions between herring (Clupea harengus) seasonal spawning aggregations

Atlantic herring (Clupea harengus), a vital ecosystem component and target of the largest Northwest Atlantic pelagic fishery, undergo seasonal spawning migrations that result in elusive sympatric population structure. Herring spawn mostly in fall or spring, and genomic differentiation was recently detected between these groups. Here we used a subset of this differentiation, 66 single nucleotide polymorphisms (SNPs) to analyze the temporal dynamics of this local adaptation and the applicability of SNP subsets in stock assessment. We showed remarkable temporal stability of genomic differentiation corresponding to spawning season, between samples taken a decade apart (2005 N = 90 vs. 2014 N = 71) in the Gulf of St. Lawrence, and new evidence of limited interbreeding between spawning components. We also examined an understudied and overexploited herring population in Bras d'Or lake (N = 97); using highly reduced SNP panels (N_SNPs > 6), we verified little‐known sympatric spawning populations within this unique inland sea. These results describe consistent local adaptation, arising from asynchronous reproduction in a migratory and dynamic marine species. Our research demonstrates the efficiency and precision of SNP‐based assessments of sympatric subpopulations; and indeed, this temporally stable local adaptation underlines the importance of such fine‐scale management practices.     

The geographical and environmental determinants of genetic diversity for four alpine conifers of the European Alps

Climate is one of the most important drivers of local adaptation in forest tree species. Standing levels of genetic diversity and structure within and among natural populations of forest trees are determined by the interplay between climatic heterogeneity and the balance between selection and gene flow. To investigate this interplay, single nucleotide polymorphisms (SNPs) were genotyped in 24 to 37 populations from four subalpine conifers, Abies alba Mill., Larix decidua Mill., Pinus cembra L. and Pinus mugo Turra, across their natural ranges in the Italian Alps and Apennines. Patterns of population structure were apparent using a Bayesian clustering program, STRUCTURE, which identified three to five genetic groups per species. Geographical correlates with these patterns, however, were only apparent for P. cembra. Multivariate environmental variables [i.e. principal components (PCs)] were subsequently tested for association with SNPs using a Bayesian generalized linear mixed model. The majority of the SNPs, ranging from six in L. decidua to 18 in P. mugo, were associated with PC1, corresponding to winter precipitation and seasonal minimum temperature. In A. alba, four SNPs were associated with PC2, corresponding to the seasonal minimum temperature. Functional annotation of those genes with the orthologs in Arabidopsis revealed several genes involved in abiotic stress response. This study provides a detailed assessment of population structure and its association with environment and geography in four coniferous species in the Italian mountains.     

Novel trophic niches drive variable progress towards ecological speciation within an adaptive radiation of pupfishes

Adaptive radiation is recognized by a rapid burst of phenotypic, ecological and species diversification. However, it is unknown whether different species within an adaptive radiation evolve reproductive isolation at different rates. We compared patterns of genetic differentiation between nascent species within an adaptive radiation of Cyprinodon pupfishes using genotyping by sequencing. Similar to classic adaptive radiations, this clade exhibits rapid morphological diversification rates and two species are novel trophic specialists, a scale‐eater and hard‐shelled prey specialist (durophage), yet the radiation is <10 000 years old. Both specialists and an abundant generalist species all coexist in the benthic zone of lakes on San Salvador Island, Bahamas. Based on 13 912 single‐nucleotide polymorphisms (SNPs), we found consistent differences in genetic differentiation between each specialist species and the generalist across seven lakes. The scale‐eater showed the greatest genetic differentiation and clustered by species across lakes, whereas durophage populations often clustered with sympatric generalist populations, consistent with parallel speciation across lakes. However, we found strong evidence of admixture between durophage populations in different lakes, supporting a single origin of this species and genome‐wide introgression with sympatric generalist populations. We conclude that the scale‐eater is further along the speciation‐with‐gene‐flow continuum than the durophage and suggest that different adaptive landscapes underlying these two niche environments drive variable progress towards speciation within the same habitat. Our previous measurements of fitness surfaces in these lakes support this conclusion: the scale‐eating fitness peak may be more distant than the durophage peak on the complex adaptive landscape driving adaptive radiation.     

Signatures of local adaptation in candidate genes of oaks (Quercus spp.) with respect to present and future climatic conditions

Testing how populations are locally adapted and predicting their response to their future environment is of key importance in view of climate change. Landscape genomics is a powerful approach to investigate genes and environmental factors involved in local adaptation. In a pooled amplicon sequencing approach of 94 genes in 71 populations, we tested whether >3500 single nucleotide polymorphisms (SNPs) in the three most common oak species in Switzerland (Quercus petraea, Q. pubescens, Q. robur) show an association with abiotic factors related to local topography, historical climate and soil characteristics. In the analysis including all species, the most frequently associated environmental factors were those best describing the habitats of the species. In the species‐specific analyses, the most important environmental factors and associated SNPs greatly differed among species. However, we identified one SNP and seven genes that were associated with the same environmental factor across all species. We finally used regressions of allele frequencies of the most strongly associated SNPs along environmental gradients to predict the risk of nonadaptedness (RONA), which represents the average change in allele frequency at climate‐associated loci theoretically required to match future climatic conditions. RONA is considerable for some populations and species (up to 48% in single populations) and strongly differs among species. Given the long generation time of oaks, some of the required allele frequency changes might not be realistic to achieve based on standing genetic variation. Hence, future adaptedness requires gene flow or planting of individuals carrying beneficial alleles from habitats currently matching future climatic conditions.     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Climate‐associated genetic variation in Fagus sylvatica and potential responses to climate change in the French Alps

Local adaptation patterns have been found in many plants and animals, highlighting the genetic heterogeneity of species along their range of distribution. In the next decades, global warming is predicted to induce a change in the selective pressures that drive this adaptive variation, forcing a reshuffling of the underlying adaptive allele distributions. For species with low dispersion capacity and long generation time such as trees, the rapidity of the change could impede the migration of beneficial alleles and lower their capacity to track the changing environment. Identifying the main selective pressures driving the adaptive genetic variation is thus necessary when investigating species capacity to respond to global warming. In this study, we investigate the adaptive landscape of Fagus sylvatica along a gradient of populations in the French Alps. Using a double‐digest restriction‐site‐associated DNA (ddRAD) sequencing approach, we identified 7,000 SNPs from 570 individuals across 36 different sites. A redundancy analysis (RDA)‐derived method allowed us to identify several SNPs that were strongly associated with climatic gradients; moreover, we defined the primary selective gradients along the natural populations of F. sylvatica in the Alps. Strong effects of elevation and humidity, which contrast north‐western and south‐eastern site, were found and were believed to be important drivers of genetic adaptation. Finally, simulations of future genetic landscapes that used these findings allowed identifying populations at risk for F. sylvatica in the Alps, which could be helpful for future management plans.     

Evidence of genomic adaptation to climate in Eucalyptus microcarpa: Implications for adaptive potential to projected climate change

Understanding whether populations can adapt in situ or whether interventions are required is of key importance for biodiversity management under climate change. Landscape genomics is becoming an increasingly important and powerful tool for rapid assessments of climate adaptation, especially in long‐lived species such as trees. We investigated climate adaptation in Eucalyptus microcarpa using the DArTseq genomic approach. A combination of F_ST outlier and environmental association analyses were performed using >4200 genomewide single nucleotide polymorphisms (SNPs) from 26 populations spanning climate gradients in southeastern Australia. Eighty‐one SNPs were identified as putatively adaptive, based on significance in F_ST outlier tests and significant associations with one or more climate variables related to temperature (70/81), aridity (37/81) or precipitation (35/81). Adaptive SNPs were located on all 11 chromosomes, with no particular region associated with individual climate variables. Climate adaptation appeared to be characterized by subtle shifts in allele frequencies, with no consistent fixed differences identified. Based on these associations, we predict adaptation under projected changes in climate will include a suite of shifts in allele frequencies. Whether this can occur sufficiently rapidly through natural selection within populations, or would benefit from assisted gene migration, requires further evaluation. In some populations, the absence or predicted increases to near fixation of particular adaptive alleles hint at potential limits to adaptive capacity. Together, these results reinforce the importance of standing genetic variation at the geographic level for maintaining species’ evolutionary potential.     

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.     

Landscape genomics and a common garden trial reveal adaptive differentiation to temperature across Europe in the tree species Alnus glutinosa

The adaptive potential of tree species to cope with climate change has important ecological and economic implications. Many temperate tree species experience a wide range of environmental conditions, suggesting high adaptability to new environmental conditions. We investigated adaptation to regional climate in the drought‐sensitive tree species Alnus glutinosa (Black alder), using a complementary approach that integrates genomic, phenotypic and landscape data. A total of 24 European populations were studied in a common garden and through landscape genomic approaches. Genotyping‐by‐sequencing was used to identify SNPs across the genome, resulting in 1990 SNPs. Although a relatively low percentage of putative adaptive SNPs was detected (2.86% outlier SNPs), we observed clear associations among outlier allele frequencies, temperature and plant traits. In line with the typical drought avoiding nature of A. glutinosa, leaf size varied according to a temperature gradient and significant associations with multiple outlier loci were observed, corroborating the ecological relevance of the observed outlier SNPs. Moreover, the lack of isolation by distance, the very low genetic differentiation among populations and the high intrapopulation genetic variation all support the notion that high gene exchange combined with strong environmental selection promotes adaptation to environmental cues.     

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.     

Time‐specific and pleiotropic quantitative trait loci coordinately modulate stem growth in Populus

In perennial woody plants, the coordinated increase of stem height and diameter during juvenile growth improves competitiveness (i.e. access to light); however, the factors underlying variation in stem growth remain unknown in trees. Here, we used linkage‐linkage disequilibrium (linkage‐LD) mapping to decipher the genetic architecture underlying three growth traits during juvenile stem growth. We used two Populus populations: a linkage mapping population comprising a full‐sib family of 1,200 progeny and an association mapping panel comprising 435 unrelated individuals from nearly the entire natural range of Populus tomentosa. We mapped 311 quantitative trait loci (QTL) for three growth traits at 12 timepoints to 42 regions in 17 linkage groups. Of these, 28 regions encompassing 233 QTL were annotated as 27 segmental homology regions (SHRs). Using SNPs identified by whole‐genome re‐sequencing of the 435‐member association mapping panel, we identified significant SNPs (P ≤ 9.4 × 10^?7) within 27 SHRs that affect stem growth at nine timepoints with diverse additive and dominance patterns, and these SNPs exhibited complex allelic epistasis over the juvenile growth period. Nineteen genes linked to potential causative alleles that have time‐specific or pleiotropic effects, and mostly overlapped with significant signatures of selection within SHRs between climatic regions represented by the association mapping panel. Five genes with potential time‐specific effects showed species‐specific temporal expression profiles during the juvenile stages of stem growth in five representative Populus species. Our observations revealed the importance of considering temporal genetic basis of complex traits, which will facilitate the molecular design of tree ideotypes.     

Environmental adaptation in Chinook salmon (Oncorhynchus tshawytscha) throughout their North American range

Landscape genomics is a rapidly growing field with recent advances in both genotyping efficiency and statistical analyses that provide insight towards local adaptation of populations under varying environmental and selective pressure. Chinook salmon (Oncorhynchus tshawytscha) are a broadly distributed Pacific salmon species, occupying a diversity of habitats throughout the northeastern Pacific with pronounced variation in environmental and climate features but little is understood regarding local adaptation in this species. We used a multivariate method, redundancy analysis (RDA), to identify polygenic correlations between 19 703 SNP loci and a suite of environmental variables in 46 collections of Chinook salmon (1956 total individuals) distributed throughout much of its North American range. Models in RDA were conducted on both rangewide and regional scales by hierarchical partitioning of the populations into three distinct genetic lineages. Our results indicate that between 5.8 and 21.8% of genomic variation can be accounted for by environmental features, and 566 putatively adaptive loci were identified as targets of environmental adaptation. The most influential drivers of adaptive divergence included precipitation in the driest quarter of the year (Rangewide and North Coastal Lineage, anova P = 0.002 and 0.01, respectively), precipitation in the wettest quarter of the year (Interior Columbia River Stream‐Type Lineage, anova P = 0.03), variation in mean diurnal range in temperature (South Coastal Lineage, anova P = 0.005), and migration distance (Rangewide, anova P = 0.001). Our results indicate that environmental features are strong drivers of adaptive genomic divergence in this species, and provide a foundation to investigate how Chinook salmon might respond to global environmental change.     

Lineage‐specific adaptation to climate involves flowering time in North American Arabidopsis lyrata

Adaptation to local climatic conditions is commonly found within species, but whether it involves the same intraspecific genomic variants is unknown. We studied this question in North American Arabidopsis lyrata, whose current distribution is shaped by post‐glacial range expansion from two refugia, resulting in two distinct genetic clusters covering comparable climatic gradients. Using pooled whole‐genome sequence data of 41 outcrossing populations, we identified loci associated with three niche‐determining climatic variables in the two clusters and compared these outliers. Little evidence was found for parallelism in climate adaptation for single nucleotide polymorphisms (SNPs) and for genes with an accumulation of outlier SNPs. Significantly increased selection coefficients supported them as candidates of climate adaptation. However, the fraction of gene ontology (GO) terms shared between clusters was higher compared to outlier SNPs and outlier genes, suggesting that selection acts on similar pathways but not necessarily the same genes. Enriched GO terms involved responses to abiotic and biotic stress, circadian rhythm and development, with flower development and reproduction being among the most frequently detected. In line with GO enrichment, regulators of flowering time were detected as outlier genes. Our results suggest that while adaptation to environmental gradients on the genomic level are lineage‐specific in A. lyrata, similar biological processes seem to be involved. Differential loss of standing genetic variation, probably driven by genetic drift, can in part account for the lack of parallel evolution on the genomic level.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

Forest tree growth response to hydroclimate variability in the southern Appalachians

Climate change will affect tree species growth and distribution; however, under the same climatic conditions species may differ in their response according to site conditions. We evaluated the climate‐driven patterns of growth for six dominant deciduous tree species in the southern Appalachians. We categorized species into two functional groups based on their stomatal regulation and xylem architecture: isohydric, diffuse porous and anisohydric, ring porous. We hypothesized that within the same climatic regime: (i) species‐specific differences in growth will be conditional on topographically mediated soil moisture availability; (ii) in extreme drought years, functional groups will have markedly different growth responses; and (iii) multiple hydroclimate variables will have direct and indirect effects on growth for each functional group. We used standardized tree‐ring chronologies to examine growth of diffuse‐porous (Acer, Liriodendron, and Betula) and ring‐porous (Quercus) species vs. on‐site climatic data from 1935 to 2003. Quercus species growing on upslope sites had higher basal area increment (BAI) than Quercus species growing on mesic, cove sites; whereas, Acer and Liriodendron had lower BAI on upslope compared to cove sites. Diffuse‐porous species were more sensitive to climate than ring porous, especially during extreme drought years. Across functional groups, radial growth was more sensitive to precipitation distribution, such as small storms and dry spell length (DSL), rather than the total amount of precipitation. Based on structural equation modeling, diffuse‐porous species on upslope sites were the most sensitive to multiple hydroclimate variables (r² = 0.46), while ring‐porous species on upslope sites were the least sensitive (r² = 0.32). Spring precipitation, vapor pressure deficit, and summer storms had direct effects on summer AET/P, and summer AET/P, growing season small storms and DSL partially explained growth. Decreasing numbers of small storms and extending the days between rainfall events will result in significant growth reduction, even in regions with relatively high total annual rainfall.