Widespread phenotypic and genetic divergence along altitudinal gradients in animals

Altitudinal gradients offer valuable study systems to investigate how adaptive genetic diversity is distributed within and between natural populations and which factors promote or prevent adaptive differentiation. The environmental clines along altitudinal gradients tend to be steep relative to the dispersal distance of many organisms, providing an opportunity to study the joint effects of divergent natural selection and gene flow. Temperature is one variable showing consistent altitudinal changes, and altitudinal gradients can therefore provide spatial surrogates for some of the changes anticipated under climate change. Here, we investigate the extent and patterns of adaptive divergence in animal populations along altitudinal gradients by surveying the literature for (i) studies on phenotypic variation assessed under common garden or reciprocal transplant designs and (ii) studies looking for signatures of divergent selection at the molecular level. Phenotypic data show that significant between‐population differences are common and taxonomically widespread, involving traits such as mass, wing size, tolerance to thermal extremes and melanization. Several lines of evidence suggest that some of the observed differences are adaptively relevant, but rigorous tests of local adaptation or the link between specific phenotypes and fitness are sorely lacking. Evidence for a role of altitudinal adaptation also exists for a number of candidate genes, most prominently haemoglobin, and for anonymous molecular markers. Novel genomic approaches may provide valuable tools for studying adaptive diversity, also in species that are not amenable to experimentation.     

Population genomic evidence for adaptive differentiation in the Baltic Sea herring

Detecting and estimating the degree of genetic differentiation among populations of highly mobile marine fish having pelagic larval stages is challenging because their effective population sizes can be large, and thus, little genetic drift and differentiation is expected in neutral genomic sites. However, genomic sites subject to directional selection stemming from variation in local environmental conditions can still show substantial genetic differentiation, yet these signatures can be hard to detect with low‐throughput approaches. Using a pooled RAD‐seq approach, we investigated genomewide patterns of genetic variability and differentiation within and among 20 populations of Atlantic herring in the Baltic Sea (and adjacent Atlantic sites), where previous low‐throughput studies and/or studies based on few populations have found limited evidence for genetic differentiation. Stringent quality control was applied in the filtering of 1 791 254 SNPs, resulting in a final data set of 68 182 polymorphic loci. Clear differentiation was identified between Atlantic and Baltic populations in many genomic sites, while differentiation within the Baltic Sea area was weaker and geographically less structured. However, outlier analyses – whether including all populations or only those within the Baltic Sea – uncovered hundreds of directionally selected loci in which variability was associated with either salinity, temperature or both. Hence, our results support the view that although the degree of genetic differentiation among Baltic Sea herring populations is low, there are many genomic regions showing elevated divergence, apparently as a response to temperature‐ and salinity‐related natural selection. As such, the results add to the increasing evidence of local adaptation in highly mobile marine organisms, and those in the young Baltic Sea in particular.     

Whither Pst? The approximation of Qst by Pst in evolutionary and conservation biology

Local adaptation through natural selection can be inferred in case the additive genetic divergence in a quantitative trait across populations (Q(st)) exceeds the neutral expectation based on differentiation of neutral alleles across these populations (e.g. F(st)). As such, measuring Q(st) in relation to neutral differentiation presents a first-line investigation applicable in evolutionary biology (selection on functional genes) and conservation biology (identification of locally adapted coding genes). However, many species, especially those in need of conservation actions, are not amenable for the kind of breeding design required to estimate either narrow- or broad-sense Q(st). In such cases, Q(st) has been approximated by the phenotypic divergence in a trait across populations (P(st)). I here argue that the critical aspect for how well P(st) approximates Q(st) depends on the extent that additive genetic effects determine variation between populations relative to within populations. I review how the sensitivity of conclusions regarding local adaptation based on P(st) have been evaluated in the literature and find that many studies make a anticonservative null assumption in estimating P(st) and/or use a nonconservative approach to explore sensitivity of their conclusions. Data from two studies that have provided a second, independent assessment of selection in their system suggest that P(st)-F(st) comparisons should be interpreted very conservatively. I conclude with recommendations for improving the robustness of the inferences drawn from comparing P(st) with neutral differentiation.     

Local adaptation and oceanographic connectivity patterns explain genetic differentiation of a marine diatom across the North Sea–Baltic Sea salinity gradient

Drivers of population genetic structure are still poorly understood in marine micro‐organisms. We exploited the North Sea–Baltic Sea transition for investigating the seascape genetics of a marine diatom, Skeletonema marinoi. Eight polymorphic microsatellite loci were analysed in 354 individuals from ten locations to analyse population structure of the species along a 1500‐km‐long salinity gradient ranging from 3 to 30 psu. To test for salinity adaptation, salinity reaction norms were determined for sets of strains originating from three different salinity regimes of the gradient. Modelled oceanographic connectivity was compared to directional relative migration by correlation analyses to examine oceanographic drivers. Population genetic analyses showed distinct genetic divergence of a low‐salinity Baltic Sea population and a high‐salinity North Sea population, coinciding with the most evident physical dispersal barrier in the area, the Danish Straits. Baltic Sea populations displayed reduced genetic diversity compared to North Sea populations. Growth optima of low salinity isolates were significantly lower than those of strains from higher native salinities, indicating local salinity adaptation. Although the North Sea–Baltic Sea transition was identified as a barrier to gene flow, migration between Baltic Sea and North Sea populations occurred. However, the presence of differentiated neutral markers on each side of the transition zone suggests that migrants are maladapted. It is concluded that local salinity adaptation, supported by oceanographic connectivity patterns creating an asymmetric migration pattern between the Baltic Sea and the North Sea, determines genetic differentiation patterns in the transition zone.     

Habitat specialization and adaptive phenotypic divergence of anuran populations

We tested for adaptive population structure in the frog Rana temporaria by rearing tadpoles from 23 populations in a common garden experiment, with and without larval dragonfly predators. The goal was to compare tadpole phenotypes with the habitats of their source ponds. The choice of traits and habitat variables was guided by prior information about phenotypic function. There were large differences among populations in life history, behaviour, morphological shape, and the predator-induced plasticities in most of these. Body size and behaviour were correlated with predation risk in the source pond, in agreement with adaptive population divergence. Tadpoles from large sunny ponds were morphologically distinct from those inhabiting small woodland ponds, although here an adaptive explanation was unclear. There was no evidence that plasticity evolves in populations exposed to more variable environments. Much among-population variation in phenotype and plasticity was not associated with habitat, perhaps reflecting rapid changes in wetland habitats.     

Adaptation genomics: next generation sequencing reveals a shared haplotype for rapid early development in geographically and genetically distant populations of rainbow trout

Local adaptation occurs when a population evolves a phenotype that confers a selective advantage in its local environment, but which may not be advantageous in other habitats. Restricted gene flow and strong selection pressures are prerequisites for local adaptation. Fishes in the family Salmonidae are predicted to provide numerous examples of local adaptation because of the high fidelity of returning to spawn in their natal streams, which results in highly structured populations, and the wide diversity of environments that salmonids have colonized. These conditions are ideally suited for producing a set of specialist phenotypes, whose fitness is maximized for one specific habitat, rather than a generalist phenotype similarly viable in several environments. Understanding patterns and processes leading to local adaptations has long been a goal of evolutionary biology, but it is only recently that identifying the molecular basis for local adaptation has become feasible because of advances in genomic technologies. The study of shared adaptive phenotypes in populations that are both geographically distant and genetically distinct should reveal some of the fundamental molecular mechanisms associated with local adaptation. In this issue of Molecular Ecology, Miller et al. (2012) make a significant contribution to the development of adaptation genomics. This study suggests that salmonids use standing genetic variation to select beneficial alleles for local adaptations rather than de novo mutations in the same gene or alternative physiological pathways. Identifying the genetic basis for local adaptation has major implications for the management, conservation and potential restoration of salmonid populations.     

Identification of candidate loci for adaptive phenotypic plasticity in natural populations of spadefoot toads

Phenotypic plasticity allows organisms to alter their phenotype in direct response to changes in the environment. Despite growing recognition of plasticity's role in ecology and evolution, few studies have probed plasticity's molecular bases—especially using natural populations. We investigated the genetic basis of phenotypic plasticity in natural populations of spadefoot toads (Spea multiplicata). Spea tadpoles normally develop into an “omnivore” morph that is favored in long‐lasting, low‐density ponds. However, if tadpoles consume freshwater shrimp or other tadpoles, they can alternatively develop (via plasticity) into a “carnivore” morph that is favored in ephemeral, high‐density ponds. By combining natural variation in pond ecology and morph production with population genetic approaches, we identified candidate loci associated with each morph (carnivores vs. omnivores) and loci associated with adaptive phenotypic plasticity (adaptive vs. maladaptive morph choice). Our candidate morph loci mapped to two genes, whereas our candidate plasticity loci mapped to 14 genes. In both cases, the identified genes tended to have functions related to their putative role in spadefoot tadpole biology. Our results thereby form the basis for future studies into the molecular mechanisms that mediate plasticity in spadefoots. More generally, these results illustrate how diverse loci might mediate adaptive plasticity.     

Shape variation in the least killifish: ecological associations of phenotypic variation and the effects of a common garden

Studies of the adaptive significance of variation among conspecific populations often focus on a single ecological factor. However, habitats rarely differ in only a single ecological factor, creating a challenge for identifying the relative importance of the various ecological factors that might be maintaining local adaptation. Here we investigate the ecological factors associated with male body shape variation among nine populations of the poeciliid fish, Heterandria formosa, from three distinct habitats and combine those results with a laboratory study of three of those populations to assess the contributions of genetic and environmental influences to shape variation. Field‐collected animals varied principally in three ways: the orientation of the gonopodium, the intromittent organ; the degree of body depth and streamlining; and the shape of the tail musculature. Fish collected in the spring season were larger and had a more anteriorly positioned gonopodium than fish collected in autumn. Fish collected from lotic springs were larger and more streamlined than those collected from lentic ponds or tidal marshes. Some of the variation in male shape among populations within habitats was associated with population‐level variation in species richness, adult density, vegetative cover, predation risk, and female standard length. Population‐level differences among males in body size, position of the gonopodium, and shape of the tail musculature were maintained among males reared in a common environment. In contrast, population variation in the degree of streamlining was eliminated when males were reared in a common environment. These results illustrate the complicated construction of multivariate phenotypic variation and suggest that different agents of selection have acted on different components of shape.     

The genetic basis of adaptive population differentiation: a quantitative trait locus analysis of fitness traits in two wild barley populations from contrasting habitats

We used a quantitative trait locus (QTL) approach to study the genetic basis of population differentiation in wild barley, Hordeum spontaneum. Several ecotypes are recognized in this model species, and population genetic studies and reciprocal transplant experiments have indicated the role of local adaptation in shaping population differences. We derived a mapping population from a cross between a coastal Mediterranean population and a steppe inland population from Israel and assessed F3 progeny fitness in the natural growing environments of the two parental populations. Dilution of the local gene pool, estimated as the proportion of native alleles at 96 marker loci in the recombinant lines, negatively affected fitness traits at both sites. QTLs for fitness traits tended to differ in the magnitude but not in the direction of their effects across sites, with beneficial alleles generally conferring a greater fitness advantage at their native site. Several QTLs showed fitness effects at one site only, but no opposite selection on individual QTLs was observed across the sites. In a common-garden experiment, we explored the hypothesis that the two populations have adapted to divergent nutrient availabilities. In the different nutrient environments of this experiment, but not under field conditions, fitness of the F3 progeny lines increased with the number of heterozygous marker loci. Comparison of QTL-effects that underlie genotype x nutrient interaction in the common-garden experiment and genotype x site interaction in the field suggested that population differentiation at the field sites may have been driven by divergent nutrient availabilities to a limited extent. Also in this experiment no QTLs were observed with opposite fitness effects in contrasting environments. Our data are consistent with the view that adaptive differentiation can be based on selection on multiple traits changing gradually along ecological gradients. This can occur without QTLs showing opposite fitness effects in the different environments, that is, in the absence of genetic trade-offs in performance between environments.     

Low genome‐wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome‐wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with F_ST statistics significantly higher than those expected on the basis of overall, genome‐wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome‐wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with F_ST‐based outlier detection tests. Instead, our analyses confirmed the lack of genome‐wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.     

The genetic differentiation at quantitative trait loci under local adaptation

Most adaptive traits are controlled by large number of genes that may all together be the targets of selection. Adaptation may thus involve multiple but not necessarily substantial allele frequency changes. This has important consequences for the detection of selected loci and implies that a quantitative genetics framework may be more appropriate than the classical 'selective sweep' paradigm. Preferred methods to detect loci involved in local adaptation are based on the detection of 'outlier' values of the allelic differentiation F(ST) . A quantitative genetics framework is adopted here to review theoretical expectations for how allelic differentiation at quantitative trait loci (F(STQ) ) relates to (i), neutral genetic differentiation (F(ST) ) and (ii), phenotypic differentiation (Q(ST) ). We identify cases where results of outlier-based methods are likely to be poor and where differentiation at selected loci conveys little information regarding local adaptation. A first case is when neutral differentiation is high, so that local adaptation does not necessitate increased differentiation. A second case is when local adaptation is reached via an increased covariance of allelic effects rather than via allele frequency changes, which is more likely under high gene flow when the number of loci is high and selection is recent. The comparison of theoretical predictions with observed data from the literature suggests that polygenic local adaptation involving only faint allele frequency changes are very likely in some species such as forest trees and for climate-related traits. Recent methodological improvements that may alleviate the weakness of F(ST) -based detection methods are presented.     

A Pleiotropic Flowering Time QTL Exhibits Gene-by-Environment Interaction for Fitness in a Perennial Grass

Appropriate flowering time is a crucial adaptation impacting fitness in natural plant populations. Although the genetic basis of flowering variation has been extensively studied, its mechanisms in nonmodel organisms and its adaptive value in the field are still poorly understood. Here, we report new insights into the genetic basis of flowering time and its effect on fitness in Panicum hallii, a native perennial grass. Genetic mapping in populations derived from inland and coastal ecotypes identified flowering time quantitative trait loci (QTL) and many exhibited extensive QTL-by-environment interactions. Patterns of segregation within recombinant hybrids provide strong support for directional selection driving ecotypic divergence in flowering time. A major QTL on chromosome 5 (q-FT5) was detected in all experiments. Fine-mapping and expression studies identified a gene with orthology to a rice FLOWERING LOCUS T-like 9 (PhFTL9) as the candidate underlying q-FT5. We used a reciprocal transplant experiment to test for local adaptation and the specific impact of q-FT5 on performance. We did not observe local adaptation in terms of fitness tradeoffs when contrasting ecotypes in home versus away habitats. However, we observed that the coastal allele of q-FT5 conferred a fitness advantage only in its local habitat but not at the inland site. Sequence analyses identified an excess of low-frequency polymorphisms at the PhFTL9 promoter in the inland lineage, suggesting a role for either selection or population expansion on promoter evolution. Together, our findings demonstrate the genetic basis of flowering variation in a perennial grass and provide evidence for conditional neutrality underlying flowering time divergence.     

Genetic basis for rapidly evolved tolerance in the wild: adaptation to toxic pollutants by an estuarine fish species

Atlantic killifish (Fundulus heteroclitus) residing in some urban and industrialized estuaries of the US eastern seaboard demonstrate recently evolved and extreme tolerance to toxic aryl hydrocarbon pollutants, characterized as dioxin‐like compounds (DLCs). Here, we provide an unusually comprehensive accounting (69%) through quantitative trait locus (QTL) analysis of the genetic basis for DLC tolerance in killifish inhabiting an urban estuary contaminated with PCB congeners, the most toxic of which are DLCs. Consistent with mechanistic knowledge of DLC toxicity in fish and other vertebrates, the aryl hydrocarbon receptor (ahr2) region accounts for 17% of trait variation; however, QTL on independent linkage groups and their interactions have even greater explanatory power (44%). QTL interpreted within the context of recently available Fundulus genomic resources and shared synteny among fish species suggest adaptation via interacting components of a complex stress response network. Some QTL were also enriched in other killifish populations characterized as DLC‐tolerant and residing in distant urban estuaries contaminated with unique mixtures of pollutants. Together, our results suggest that DLC tolerance in killifish represents an emerging example of parallel contemporary evolution that has been driven by intense human‐mediated selection on natural populations.     

Evidence of adaptive divergence in plasticity: density- and site-dependent selection on shade-avoidance responses in Impatiens capensis

We investigated the conditions under which plastic responses to density are adaptive in natural populations of Impatiens capensis and determined whether plasticity has evolved differently in different selective environments. Previous studies showed that a population that evolved in a sunny site exhibited greater plasticity in response to density than did a population that evolved in a woodland site. Using replicate inbred lines in a reciprocal transplant that included a density manipulation, we asked whether such population differentiation was consistent with the hypothesis of adaptive divergence. We hypothesized that plasticity would be more strongly favored in the sunny site than in the woodland site; consequently, we predicted that selection would be more strongly density dependent in the sunny site, favoring the phenotype that was expressed at each density. Selection on internode length and flowering date was consistent with the hypothesis of adaptive divergence in plasticity. Few costs or benefits of plasticity were detected independently from the expressed phenotype, so plasticity was selected primarily through selection on the phenotype. Correlations between phenotypes and their plasticity varied with the environment and would cause indirect selection on plasticity to be environment dependent. We showed that an appropriate plastic response even to a rare environment can greatly increase genotypic fitness when that environment is favorable. Selection on the measured characters contributed to local adaptation and fully accounted for fitness differences between populations in all treatments except the woodland site at natural density.     

Is local adaptation in Mimulus guttatus caused by trade‐offs at individual loci?

Local adaptation is considered to be the result of fitness trade‐offs for particular phenotypes across different habitats. However, it is unclear whether such phenotypic trade‐offs exist at the level of individual genetic loci. Local adaptation could arise from trade‐offs of alternative alleles at individual loci or by complementary sets of loci with different fitness effects of alleles in one habitat but selective neutrality in the alternative habitat. To evaluate the genome‐wide basis of local adaptation, we performed a field‐based quantitative trait locus (QTL) mapping experiment on recombinant inbred lines (RILs) created from coastal perennial and inland annual races of the yellow monkeyflower (Mimulus guttatus) grown reciprocally in native parental habitats. Overall, we detected 19 QTLs affecting one or more of 16 traits measured in two environments, most of small effect. We identified 15 additional QTL effects at two previously identified candidate QTLs [DIV ERGENCE (DIV)]. Significant QTL by environment interactions were detected at the DIV loci, which was largely attributable to genotypic differences at a single field site. We found no detectable evidence for trade‐offs for any one component of fitness, although DIV2 showed a trade‐off involving different fitness traits between sites, suggesting that local adaptation is largely controlled by non‐overlapping loci. This is surprising for an outcrosser, implying that reduced gene flow prevents the evolution of individuals adapted to multiple environments. We also determined that native genotypes were not uniformly adaptive, possibly reflecting fixed mutational load in one of the populations.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Relaxed predation results in reduced phenotypic integration in a suite of dragonflies

Although changes in magnitude of single traits responding to selective agents have been studied intensively, little is known about selection shaping networks of traits and their patterns of covariation. However, this is central for our understanding of phenotypic evolution as traits are embedded in a multivariate environment with selection affecting a multitude of traits simultaneously rather than individually. Here, we investigate inter‐ and intraspecific patterns of trait integration (trait correlations) in the larval abdomen of dragonflies as a response to a change in predator selection. Species of the dragonfly genus Leucorrhinia underwent a larval habitat shift from predatory fish to predatory dragonfly‐dominated lakes with an associated relaxation in selection pressure from fish predation. Our results indicate that the habitat‐shift‐induced relaxed selection pressure caused phenotypic integration of abdominal traits to be reduced. Intraspecific findings matched patterns comparing species from both habitats with higher abdominal integration in response to predatory fish. This higher integration is probably a result of faster burst swimming speed. The abdomen holds the necessary morphological machinery to successfully evade predatory fish via burst swimming. Hence, abdominal traits have to function in a tight coordinated manner, as maladaptive variation and consequently nonoptimal burst swimming would cause increased mortality. In predatory dragonfly‐dominated lakes, no such strong link between burst swimming and mortality is present. Our findings highlight the importance of studying multivariate trait relationships as a response to selection for understanding patterns of phenotypic diversification.     

The degree of adaptive phenotypic plasticity is correlated with the spatial environmental heterogeneity experienced by island populations of Rana temporaria

Although theoretical models have identified environmental heterogeneity as a prerequisite for the evolution of adaptive plasticity, this relationship has not yet been demonstrated experimentally. Because of pool desiccation risk, adaptation of development rate is important for many amphibians. In a simulated pool-drying experiment, we compared the development time and phenotypic plasticity in development time of populations of the common frog Rana temporaria, originating from 14 neighbouring islands off the coast of northern Sweden. Drying regime of pools used by frogs for breeding differed within and among the islands. We found that the degree of phenotypic plasticity in development time was positively correlated with the spatial variation in the pool-drying regimes present on each island. In addition, local adaptation in development time to the mean drying rate of the pools on each island was found. Hence, our study demonstrates the connection between environmental heterogeneity and developmental plasticity at the island population level, and also highlights the importance of the interplay between local specialization and phenotypic plasticity depending on the local selection pressures.