Spatial and environmental effects on Coho Salmon life history trait variation

Adult size, egg size, fecundity, and mass of gonads are affected by trade‐offs between reproductive investment and environmental conditions shaping the evolution of life history traits among populations for widely distributed species. Coho salmon Oncorhynchus kisutch have a large geographic distribution, and different environmental conditions are experienced by populations throughout their range. We examined the effect of environmental variables on female size, egg size, fecundity, and reproductive investment of populations of Coho Salmon from across British Columbia using an information theoretic approach. Female size increased with latitude and decreased with migration distance from the ocean to spawning locations. Egg size was lowest for intermediate intragravel temperature during incubation, decreased with migration distance, but increased in rivers below lakes. Fecundity increased with latitude, warmer temperature during the spawning period, and river size, but decreased in rivers below lakes compared with rivers with tributary sources. Relative gonad size increased with latitude and decreased with migration distance. Latitude of spawning grounds, migratory distance, and temperatures experienced by a population, but also hydrologic features—river size and headwater source—are influential in shaping patterns of reproductive investment, particularly egg size. Although, relative gonad size varied with latitude and migration distance, how gonadal mass was partitioned gives insight into the trade‐off between egg size and fecundity. The lack of an effect of latitude on egg size suggests that local optima for egg size related to intragravel temperature may drive the variation in fecundity observed among years.     

Disentangling complex genomic signals to understand population structure of an exploited,estuarine‐dependent flatfish

Interpreting contemporary patterns of population structure requires an understanding of the interactions among microevolutionary forces and past demographic events. Here, 4,122 SNP‐containing loci were used to assess structure in southern flounder (Paralichthys lethostigma) sampled across its range in the US Atlantic Ocean (Atlantic) and Gulf of Mexico (Gulf) and relationships among components of genomic variation and spatial and environmental variables were assessed across estuarine population samples in the Gulf. While hierarchical amova revealed significant heterogeneity within and between the Atlantic and Gulf, pairwise comparisons between samples within ocean basins demonstrated that all significant heterogeneity occurred within the Gulf. The distribution of Tajima''s D estimated at a genome‐wide scale differed significantly from equilibrium in all estuaries, with more negative values occurring in the Gulf. Components of genomic variation were significantly associated with environmental variables describing individual estuaries, and environment explained a larger component of variation than spatial proximity. Overall, results suggest that there is genetic spatial autocorrelation caused by shared larval sources for proximal nurseries (migration/drift), but that it is modified by environmentally driven differentiation (selection). This leads to conflicting signals in different parts of the genome and creates patterns of divergence that do not correspond to paradigms of strong local directional selection.     

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.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

Detection of environmental and morphological adaptation despite high landscape genetic connectivity in a pest grasshopper (Phaulacridium vittatum)

Widespread species that exhibit both high gene flow and the capacity to occupy heterogeneous environments make excellent models for examining local selection processes along environmental gradients. Here we evaluate the influence of temperature and landscape variables on genetic connectivity and signatures of local adaptation in Phaulacridium vittatum, a widespread agricultural pest grasshopper, endemic to Australia. With sampling across a 900‐km latitudinal gradient, we genotyped 185 P. vittatum from 19 sites at 11,408 single nucleotide polymorphisms (SNPs) using ddRAD sequencing. Despite high gene flow across sites (pairwise F_ST = 0.0003–0.08), landscape genetic resistance modelling identified a positive nonlinear effect of mean annual temperature on genetic connectivity. Urban areas and water bodies had a greater influence on genetic distance among sites than pasture, agricultural areas and forest. Together, F_ST outlier tests and environmental association analysis (EAA) detected 242 unique SNPs under putative selection, with the highest numbers associated with latitude, mean annual temperature and body size. A combination of landscape genetic connectivity analysis together with EAA identified mean annual temperature as a key driver of both neutral gene flow and environmental selection processes. Gene annotation of putatively adaptive SNPs matched with gene functions for olfaction, metabolic detoxification and ultraviolet light shielding. Our results imply that this widespread agricultural pest has the potential to spread and adapt under shifting temperature regimes and land cover change.     

Landscape features along migratory routes influence adaptive genomic variation in anadromous steelhead (Oncorhynchus mykiss)

Organisms typically show evidence of adaptation to features within their local environment. However, many species undergo long‐distance dispersal or migration across larger geographic regions that consist of highly heterogeneous habitats. Therefore, selection may influence adaptive genetic variation associated with landscape features at residing sites and along migration routes in migratory species. We tested for genomic adaptation to landscape features at natal spawning sites and along migration paths to the ocean of anadromous steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin. Results from multivariate ordination, gene–environment association and outlier analyses using 24,526 single nucleotide polymorphisms (SNPs) provided evidence that adaptive allele frequencies were more commonly associated with landscape features along migration paths than features at natal sites (91.8% vs. 8.2% of adaptive loci, respectively). Among the 45 landscape variables tested, migration distance to the ocean and mean annual precipitation along migration paths were significantly associated with adaptive genetic variation in three distinct genetic groups. Additionally, variables such as minimum migration water temperature and mean migration slope were significant only in inland stocks of steelhead that migrate up to 1,200 km farther than those near the coast, indicating regional differences in migratory selective pressures. This study provides novel approaches for investigating migratory corridors and some of the first evidence that environment along migration paths can lead to substantial divergent selection. Consequently, our approach to understand genetic adaptation to migration conditions can be applied to other migratory species when migration or dispersal paths are generally known.     

Macroscale intraspecific variation and environmental heterogeneity: analysis of cold and warm zone abundance,mortality, and regeneration distributions of four eastern US tree species

I test for macroscale intraspecific variation of abundance, mortality, and regeneration of four eastern US tree species (Tsuga canadensis,Betula lenta,Liriodendron tulipifera, and Quercus prinus) by splitting them into three climatic zones based on plant hardiness zones (PHZs). The primary goals of the analysis are to assess the differences in environmental heterogeneity and demographic responses among climatic zones, map regional species groups based on decision tree rules, and evaluate univariate and multivariate patterns of species demography with respect to environmental variables. I use the Forest Inventory Analysis (FIA) data to derive abundance, mortality, and regeneration indices and split the range into three climatic zones based on USDA PHZs: (1) cold adapted, leading region; (2) middle, well‐adapted region; and (3) warm adapted, trailing region. I employ decision tree ensemble methods to assess the importance of environmental predictors on the abundance of the species between the cold and warm zones and map zonal variations in species groups. Multivariate regression trees are used to simultaneously explore abundance, mortality, and regeneration in tandem to assess species vulnerability. Analyses point to the relative importance of climate in the warm adapted, trailing zone (especially moisture) compared to the cold adapted, leading zone. Higher mortality and lower regeneration patterns in the warm trailing zone point to its vulnerability to growing season temperature and precipitation changes that could figure more prominently in the future. This study highlights the need to account for intraspecific variation of demography in order to understand environmental heterogeneity and differential adaptation. It provides a methodology for assessing the vulnerability of tree species by delineating climatic zones based on easily available PHZ data, and FIA derived abundance, mortality, and regeneration indices as a proxy for overall growth and fitness. Based on decision tree rules, ecologically meaningful variations in species abundance among the climatic zones can be related to environmental variability and mapped.     

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time‐consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome–environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described “islands of differentiation,” and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.     

Population genomics of marine fishes: identifying adaptive variation in space and time

Studies of adaptive evolution have experienced a recent revival in population genetics of natural populations and there is currently much focus on identifying genomic signatures of selection in space and time. Insights into local adaptation, adaptive response to global change and evolutionary consequences of selective harvesting can be generated through population genomics studies, allowing the separation of the effects invoked by neutral processes (drift-migration) from those due to selection. Such knowledge is important not only for improving our basic understanding of natural as well as human-induced evolutionary processes, but also for predicting future trajectories of biodiversity and for setting conservation priorities. Marine fishes possess a number of features rendering them well suited for providing general insights into adaptive genomic evolution in natural populations. These include well-described population structures, substantial and rapidly developing genomic resources and abundant archived samples enabling temporal studies. Furthermore, superior possibilities for conducting large-scale experiments under controlled conditions, due to the economic resources provided by the large and growing aquaculture industry, hold great promise for utilizing recent technological developments. Here, we review achievements in marine fish genomics to date and highlight potential avenues for future research, which will provide both general insights into evolution in high gene flow species, as well as specific knowledge which can lead to improved management of marine organisms.     

Temporal and environmental influences on the variation in Atlantic salmon smolt migration in the Burrishoole system 1970–2000

The relationships between a number of environmental variables and the number of Atlantic salmon Salmo salar smolts migrating in the Burrishoole system, western Ireland, were examined over a 30 year period from 1970 to 2000. The number of Atlantic salmon smolts recorded migrating downstream decreased significantly from an annual mean of 11 579 in the 1970s to a mean of 6272 in the 1990s. The primary factor in the decline in Atlantic salmon numbers was a consistent decline in the number of returning adults from the 1970s until the mid 1990s. Timing of the smolt migrations (runs) was consistent throughout the three decades. The mean durations of the smolt runs were 102·1, 92·6 and 103·2 days for the 1970s, 1980s and 1990s respectively. Mean surface water temperatures at key points in the smolt run for the three decades were also similar. Mean water temperature values were 5·4, 5·1 and 5·3° C at the start of the smolt runs and 15·2, 14·9 and 15·3° C at the end of the smolt runs in the 1970s, 1980s and 1990s respectively. Multivariate analysis identified two groups of environmental variables which had a significant influence on the daily smolt catch. One group of variables dominated by photoperiod and temperature operated prior to the smolt run and was considered to regulate the development of smoltification. The second group of variables dominated by total light and water level operated within the smolt run and was considered to control daily smolt migration.     

Exome sequencing of deer mice on two California Channel Islands identifies potential adaptation to strongly contrasting ecological conditions

Understanding the forces that drive genotypic and phenotypic change in wild populations is a central goal of evolutionary biology. We examined exome variation in populations of deer mice from two of the California Channel Islands: Peromyscus maniculatus elusus from Santa Barbara Island and P. m. santacruzae from Santa Cruz Island exhibit significant differences in olfactory predator recognition, activity timing, aggressive behavior, morphology, prevalence of Sin Nombre virus, and population densities. We characterized variation in protein‐coding regions using exome capture and sequencing of 25 mice from Santa Barbara Island and 22 mice from Santa Cruz Island. We identified and examined 386,256 SNPs using three complementary methods (BayeScan, pcadapt, and LFMM). We found strong differences in molecular variation between the two populations and 710 outlier SNPs in protein‐coding genes that were detected by all three methods. We identified 35 candidate genes from this outlier set that were related to differences in phenotypes between island populations. Enrichment analyses demonstrated that patterns of molecular variation were associated with biological processes related to response to chemical stimuli and regulation of immune processes. Candidate genes associated with olfaction (Gfy, Tlr2, Vmn13r2, numerous olfactory receptor genes), circadian activity (Cry1), anxiety (Brca1), immunity (Cd28, Eif2ak4, Il12a, Syne1), aggression (Cyp19a, Lama2), and body size (Bc16, Syne1) exhibited non‐synonymous mutations predicted to have moderate to large effects. Variation in olfaction‐related genes, including a stop codon in the Santa Barbara Island population, suggests loss of predator‐recognition traits at the molecular level, consistent with a lack of behavioral aversion to fox feces. These findings also suggest that divergent pathogen prevalence and population density may have influenced adaptive immunity and behavioral phenotypes, such as reduced aggression. Overall, our study indicates that ecological differences between islands are associated with signatures of selection in protein‐coding genes underlying phenotypes that promote success in those environments.     

Adaptation to abiotic conditions drives local adaptation in bacteria and viruses coevolving in heterogeneous environments

Parasite local adaptation, the greater performance of parasites on their local compared with foreign hosts, has important consequences for the maintenance of diversity and epidemiology. While the abiotic environment may significantly affect local adaptation, most studies to date have failed either to incorporate the effects of the abiotic environment, or to separate them from those of the biotic environment. Here, we tease apart biotic and abiotic components of local adaptation using the bacterium Pseudomonas fluorescens and its viral parasite bacteriophage Φ2. We coevolved replicate populations of bacteria and phages at three different temperatures, and determined their performance against coevolutionary partners from the same and different temperatures. Crucially, we measured performance at different assay temperatures, which allowed us to disentangle adaptation to biotic and abiotic habitat components. Our results show that bacteria and phages are more resistant and infectious, respectively, at the temperature at which they previously coevolved, confirming that local adaptation to abiotic conditions can play a crucial role in determining parasite infectivity and host resistance. Our work underlines the need to assess host–parasite interactions across multiple relevant abiotic environments, and suggests that microbial adaption to local temperatures can create ecological barriers to dispersal across temperature gradients.     

Coevolution of adaptive technology,maladaptive culture and population size in a producer–scrounger game

Technology (i.e. tools, methods of cultivation and domestication, systems of construction and appropriation, machines) has increased the vital rates of humans, and is one of the defining features of the transition from Malthusian ecological stagnation to a potentially perpetual rising population growth. Maladaptations, on the other hand, encompass behaviours, customs and practices that decrease the vital rates of individuals. Technology and maladaptations are part of the total stock of culture carried by the individuals in a population. Here, we develop a quantitative model for the coevolution of cumulative adaptive technology and maladaptive culture in a ‘producer–scrounger’ game, which can also usefully be interpreted as an ‘individual–social’ learner interaction. Producers (individual learners) are assumed to invent new adaptations and maladaptations by trial-and-error learning, insight or deduction, and they pay the cost of innovation. Scroungers (social learners) are assumed to copy or imitate (cultural transmission) both the adaptations and maladaptations generated by producers. We show that the coevolutionary dynamics of producers and scroungers in the presence of cultural transmission can have a variety of effects on population carrying capacity. From stable polymorphism, where scroungers bring an advantage to the population (increase in carrying capacity), to periodic cycling, where scroungers decrease carrying capacity, we find that selection-driven cultural innovation and transmission may send a population on the path of indefinite growth or to extinction.     

A practical guide to environmental association analysis in landscape genomics

Landscape genomics is an emerging research field that aims to identify the environmental factors that shape adaptive genetic variation and the gene variants that drive local adaptation. Its development has been facilitated by next‐generation sequencing, which allows for screening thousands to millions of single nucleotide polymorphisms in many individuals and populations at reasonable costs. In parallel, data sets describing environmental factors have greatly improved and increasingly become publicly accessible. Accordingly, numerous analytical methods for environmental association studies have been developed. Environmental association analysis identifies genetic variants associated with particular environmental factors and has the potential to uncover adaptive patterns that are not discovered by traditional tests for the detection of outlier loci based on population genetic differentiation. We review methods for conducting environmental association analysis including categorical tests, logistic regressions, matrix correlations, general linear models and mixed effects models. We discuss the advantages and disadvantages of different approaches, provide a list of dedicated software packages and their specific properties, and stress the importance of incorporating neutral genetic structure in the analysis. We also touch on additional important aspects such as sampling design, environmental data preparation, pooled and reduced‐representation sequencing, candidate‐gene approaches, linearity of allele–environment associations and the combination of environmental association analyses with traditional outlier detection tests. We conclude by summarizing expected future directions in the field, such as the extension of statistical approaches, environmental association analysis for ecological gene annotation, and the need for replication and post hoc validation studies.     

Drivers and consequences of partial migration in an alpine bird species

1. Partial migration, where a portion of the population migrates between winter and summer (breeding) areas and the rest remain year‐round resident, is a common phenomenon across several taxonomic groups. Several hypotheses have been put forward to explain why some individuals migrate while others stay resident, as well as the fitness consequences of the different strategies. Yet, the drivers and consequences of the decision to migrate or not are poorly understood.
2. We used data from radio‐tagged female (n = 73) willow ptarmigan Lagopus lagopus in an alpine study area in Central Norway to test if (i) the decision to migrate was dependent on individual state variables (age and body weight), (ii) individuals repeated migratory decisions between seasons, and (iii) the choice of migratory strategy was related to reproductive success.
3. Partially supporting our prediction that migratory strategy depends on individual state, we found that juvenile birds with small body sizes were more likely to migrate, whereas large juveniles remained resident. For adult females, we found no relationship between the decision to migrate or stay resident and body weight. We found evidence for high individual repeatability of migratory decision between seasons. Migratory strategy did not explain variation in clutch size or nest fate among individuals, suggesting no direct influence of the chosen strategy on reproductive success.
4. Our results indicate that partial migration in willow ptarmigan is related to juvenile body weight, and that migratory behavior becomes a part of the individual life history as a fixed strategy. Nesting success was not affected by migratory strategy in our study population, but future studies should assess other traits to further test potential fitness consequences.

     

Candidate gene SNP variation in floodplain populations of pedunculate oak (Quercus robur L.) near the species' southern range margin: Weak differentiation yet distinct associations with water availability

Populations residing near species' low‐latitude range margins (LLMs) often occur in warmer and drier environments than those in the core range. Thus, their genetic composition could be shaped by climatic drivers that differ from those occurring at higher latitudes, resulting in potentially adaptive variants of conservation value. Such variants could facilitate the adaptation of populations from other portions of the geographical range to similar future conditions anticipated under ongoing climate change. However, very few studies have assessed standing genetic variation at potentially adaptive loci in natural LLM populations. We investigated standing genetic variation at single nucleotide polymorphisms (SNPs) located within 117 candidate genes and its links to putative climatic selection pressures across 19 pedunculate oak (Quercus robur L.) populations distributed along a regional climatic gradient near the species' southern range margin in southeastern Europe. These populations are restricted to floodplain forests along large lowland rivers, whose hydric regime is undergoing significant shifts under modern rapid climate change. The populations showed very weak geographical structure, suggesting extensive genetic connectivity and gene flow or shared ancestry. We identified eight (6.2%) positive F_ST‐outlier loci, and genotype–environment association analyses revealed consistent associations between SNP allele frequencies and several climatic variables linked to water availability. A total of 61 associations involving 37 SNPs (28.5%) from 35 annotated genes provided important insights into putative functional mechanisms in our system. Our findings provide empirical support for the role of LLM populations as sources of potentially adaptive variation that could enhance species’ resilience to climate change‐related pressures.     

Signatures of local adaptation along environmental gradients in a range‐expanding damselfly (Ischnura elegans)

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range‐expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site‐associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V‐ATPase) and visual processes (long‐wavelength‐sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate‐mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.     

Effects of spawning Pacific salmon on terrestrial invertebrates: Insects near spawning habitat are isotopically enriched with nitrogen‐15 but display no differences in body size

When Pacific salmon (Oncorhynchus spp.) spawn and die, they deliver marine‐derived nutrient subsidies to freshwater and riparian ecosystems. These subsidies can alter the behavior, productivity, and abundance of recipient species and their habitats. Isotopes, such as nitrogen‐15 (¹⁵N), are often used to trace the destination of marine‐derived nutrients in riparian habitats. However, few studies have tested for correlations between stable isotopes and physiological responses of riparian organisms. We examined whether increases in δ ¹⁵N in terrestrial insect bodies adjacent to salmon spawning habitat translate to changes in percent nitrogen content and body size. This involved comparisons between distance from a salmon‐bearing river, marine‐derived nutrients in soils and insects, soil moisture content, and body size and nitrogen content in two common beetle families (Coleoptera: Curculionidae, Carabidae). As predicted, δ¹⁵N in riparian soils attenuated with distance from the river but was unaffected by soil moisture. This gradient was mirrored by δ¹⁵N in the herbivorous curculionid beetles, whereas carabid beetles, which feed at a higher trophic level and are more mobile, did not show discernable patterns in their δ¹⁵N content. Additionally, neither distance from the river nor body δ¹⁵N content was related to beetle body size. We also found that nitrogen‐15 was not correlated with total percent nitrogen in insect bodies, meaning that the presence of spawning salmon did not increase the percent nitrogen content of these insects. We conclude that while salmon‐derived nutrients had entered terrestrial food webs, the presence of δ¹⁵N alone did not indicate meaningful physiological changes in these insects in terms of percent nitrogen nor body size. While stable isotopes may be useful tracers of marine‐derived nutrients, they cannot necessarily be used as a proxy for physiologically important response variables.     

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.