Analysis of character divergence along environmental gradients and other covariates

Character displacement is typically identified by comparing phenotypic differences in sympatry and allopatry. Recently, however, Goldberg and Lande (2006) pointed out that when phenotypic characters vary along an environmental gradient, the standard approach may fail to identify sympatric character divergence. Here we present a general analytical procedure for identifying sympatric character divergence while accounting for phenotypic changes that covary with environmental variables. Our approach uses residual randomization from a generalized linear model, and allows the statistical comparison of sympatric phenotypic divergence to allopatric phenotypic divergence while accounting for phenotypic variation along a gradient. Through simulation we demonstrate that our approach correctly identifies patterns of sympatric character divergence when they are present, and does not identify such patterns when they are not. Our analytical approach complements and extends the suggestions of Goldberg and Lande (2006), by allowing a full statistical assessment of the varied patterns of character displacement along environmental gradients, or while accounting for other covariates and sources of variation.     

Divergent selection along elevational gradients promotes genetic and phenotypic disparities among small mammal populations

Species distributed along mountain slopes, facing contrasting habitats in short geographic scale, are of particular interest to test how ecologically based divergent selection promotes phenotypic and genetic disparities as well as to assess isolation‐by‐environment mechanisms. Here, we conduct the first broad comparative study of phenotypic variation along elevational gradients, integrating a large array of ecological predictors and disentangling population genetic driver processes. The skull form of nine ecologically distinct species distributed over a large altitudinal range (100–4200 m) was compared to assess whether phenotypic divergence is a common phenomenon in small mammals and whether it shows parallel patterns. We also investigated the relative contribution of biotic (competition and predation) and abiotic parameters on phenotypic divergence via mixed models. Finally, we assessed the population genetic structure of a rodent species (Niviventer confucianus) via analysis of molecular variance and F_ST along three mountain slopes and tested the isolation‐by‐environment hypothesis using Mantel test and redundancy analysis. We found a consistent phenotypic divergence and marked genetic structure along elevational gradients; however, the species showed mixed patterns of size and skull shape trends across mountain zones. Individuals living at lower altitudes differed greatly in both phenotype and genotype from those living at high elevations, while middle‐elevation individuals showed more intermediate forms. The ecological parameters associated with phenotypic divergence along elevation gradients are partly related to species' ecological and evolutionary constraints. Fossorial and solitary animals are mainly affected by climatic factors, while terrestrial and more gregarious species are influenced by biotic and abiotic parameters. A novel finding of our study is that predator richness emerged as an important factor associated with the intraspecific diversification of the mammalian skull along elevational gradients, a previously overlooked parameter. Population genetic structure was mainly driven by environmental heterogeneity along mountain slopes, with no or a week spatial effect, fitting the isolation‐by‐environment scenario. Our study highlights the strong and multifaceted effects of heterogeneous steep habitats and ecologically based divergent selective forces in small mammal populations.     

Interspecific Competition, Environmental Gradients, Gene Flow, and the Coevolution of Species' Borders

Darwin viewed species range limits as chiefly determined by an interplay between the abiotic environment and interspecific interactions. Haldane argued that species' ranges could be set intraspecifically when gene flow from a species' populous center overwhelms local adaptation at the periphery. Recently, Kirkpatrick and Barton have modeled Haldane's process with a quantitative genetic model that combines density-dependent local population growth with dispersal and gene flow across a linear environmental gradient in optimum phenotype. To address Darwin's ideas, we have extended the Kirkpatrick and Barton model to include interspecific competition and the frequency-dependent selection that it generates, as well as stabilizing selection on a quantitative character. Our model includes local population growth, movements over space, natural selection, and gene flow. It simultaneously addresses the evolution of character displacement and species borders. It reproduces the Kirkpatrick and Barton single-species result that limited ranges can be produced with sufficiently steep environmental gradients and strong dispersal. Further, in the absence of environmental gradients or barriers to dispersal, interspecific competition will not limit species ranges at evolutionary equilibrium. However, interspecific competition can interact with environmental gradients and gene flow to generate limited ranges with much less extreme gradient and dispersal parameters than in the single-species case. Species display character displacement in sympatry, yet the reduction in competition that results from this displacement does not necessarily allow the two species to become sympatric everywhere. When species meet, competition reduces population densities in the region of overlap, which, in turn, intensifies the asymmetry in gene flow from center to margin. This reduces the ability of each species to adapt to local physical conditions at their range limits. If environmental gradients are monotonic but not linear, the transition zone between species at coevolutionary equilibrium occurs where the environmental gradient is steepest. If productivity gradients are also introduced into the model, then patterns similar to Rapoport's rule emerge. Interacting species respond to climate change, as it affects the optimal phenotype over space, by a combination of range shifts and local evolution in mean phenotype, while solitary species respond solely by range shifts. Finally, we compare empirical estimates for intrinsic growth rates and diffusion coefficients for several species to those needed by the single-species model to produce a stable limited range. These empirical values are generally insufficient to produce limited ranges in the model suggesting a role for interspecific interactions.     

Adaptive traits are maintained on steep selective gradients despite gene flow and hybridization in the intertidal zone

Gene flow among hybridizing species with incomplete reproductive barriers blurs species boundaries, while selection under heterogeneous local ecological conditions or along strong gradients may counteract this tendency. Congeneric, externally-fertilizing fucoid brown algae occur as distinct morphotypes along intertidal exposure gradients despite gene flow. Combining analyses of genetic and phenotypic traits, we investigate the potential for physiological resilience to emersion stressors to act as an isolating mechanism in the face of gene flow. Along vertical exposure gradients in the intertidal zone of Northern Portugal and Northwest France, the mid-low shore species Fucus vesiculosus, the upper shore species Fucus spiralis, and an intermediate distinctive morphotype of F. spiralis var. platycarpus were morphologically characterized. Two diagnostic microsatellite loci recovered 3 genetic clusters consistent with prior morphological assignment. Phylogenetic analysis based on single nucleotide polymorphisms in 14 protein coding regions unambiguously resolved 3 clades; sympatric F. vesiculosus, F. spiralis, and the allopatric (in southern Iberia) population of F. spiralis var. platycarpus. In contrast, the sympatric F. spiralis var. platycarpus (from Northern Portugal) was distributed across the 3 clades, strongly suggesting hybridization/introgression with both other entities. Common garden experiments showed that physiological resilience following exposure to desiccation/heat stress differed significantly between the 3 sympatric genetic taxa; consistent with their respective vertical distribution on steep environmental clines in exposure time. Phylogenetic analyses indicate that F. spiralis var. platycarpus is a distinct entity in allopatry, but that extensive gene flow occurs with both higher and lower shore species in sympatry. Experimental results suggest that strong selection on physiological traits across steep intertidal exposure gradients acts to maintain the 3 distinct genetic and morphological taxa within their preferred vertical distribution ranges. On the strength of distributional, genetic, physiological and morphological differences, we propose elevation of F. spiralis var. platycarpus from variety to species level, as F. guiryi.     

Location-specific sympatric morphological divergence as a possible response to species interactions in West Virginia Plethodon salamander communities

1. The competitive interactions of closely related species have long been considered important determinants of community composition and a major cause of phenotypic diversification. However, while patterns such as character displacement are well documented, less is known about how local adaptation influences diversifying selection from interspecific competition. 2. We examined body size and head shape variation among allopatric and sympatric populations of two salamander species, the widespread Plethodon cinereus and the geographically restricted P. nettingi. We quantified morphology from 724 individuals from 20 geographical localities throughout the range of P. nettingi. 3. Plethodon nettingi was more robust in cranial morphology relative to P. cinereus, and sympatric localities were more robust relative to allopatric localities. Additionally, there was significantly greater sympatric head shape divergence between species relative to allopatric communities, and sympatric localities of P. cinereus exhibited greater morphological variation than sympatric P. nettingi. 4. The sympatric morphological divergence and increase in cranial robustness of one species (P. nettingi) were similar to observations in other Plethodon communities, and were consistent with the hypothesis of interspecific competition. These findings suggest that interspecific competition in Plethodon may play an important role in phenotypic diversification in this group. 5. The increase in among-population variance in sympatric P. cinereus suggests a species-specific response to divergent natural selection that is influenced in part by other factors. We hypothesize that enhanced morphological flexibility and ecological tolerance allow P. cinereus to more rapidly adapt to local environmental conditions, and initial differences among populations have allowed the evolutionary response of P. cinereus to vary across replicate sympatric locations, resulting in distinct evolutionary trajectories of morphological change.     

Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P_ST) which largely exceeds neutral genetic differentiation (F_ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.     

Molecular divergence in tropical tree populations occupying environmental mosaics

Unveiling the genetic basis of local adaptation to environmental variation is a major goal in molecular ecology. In rugged landscapes characterized by environmental mosaics, living populations and communities can experience steep ecological gradients over very short geographical distances. In lowland tropical forests, interspecific divergence in edaphic specialization (for seasonally flooded bottomlands and seasonally dry terra firme soils) has been proven by ecological studies on adaptive traits. Some species are nevertheless capable of covering the entire span of the gradient; intraspecific variation for adaptation to contrasting conditions may explain the distribution of such ecological generalists. We investigated whether local divergence happens at small spatial scales in two stands of Eperua falcata (Fabaceae), a widespread tree species of the Guiana Shield. We investigated Single Nucleotide Polymorphisms (SNP) and sequence divergence as well as spatial genetic structure (SGS) at four genes putatively involved in stress response and three genes with unknown function. Significant genetic differentiation was observed among sub‐populations within stands, and eight SNP loci showed patterns compatible with disruptive selection. SGS analysis showed genetic turnover along the gradients at three loci, and at least one haplotype was found to be in repulsion with one habitat. Taken together, these results suggest genetic differentiation at small spatial scale in spite of gene flow. We hypothesize that heterogeneous environments may cause molecular divergence, possibly associated to local adaptation in E. falcata.     

Evolution of phenotypic clusters through competition and local adaptation along an environmental gradient

We have analyzed the evolution of a quantitative trait in populations that are spatially extended along an environmental gradient, with gene flow between nearby locations. In the absence of competition, there is stabilizing selection toward a locally best-adapted trait that changes gradually along the gradient. According to traditional ideas, gradual spatial variation in environmental conditions is expected to lead to gradual variation in the evolved trait. A contrasting possibility is that the trait distribution instead breaks up into discrete clusters. Doebeli and Dieckmann (2003) argued that competition acting locally in trait space and geographical space can promote such clustering. We have investigated this possibility using deterministic population dynamics for asexual populations, analyzing our model numerically and through an analytical approximation. We examined how the evolution of clusters is affected by the shape of competition kernels, by the presence of Allee effects, and by the strength of gene flow along the gradient. For certain parameter ranges clustering was a robust outcome, and for other ranges there was no clustering. Our analysis shows that the shape of competition kernels is important for clustering: the sign structure of the Fourier transform of a competition kernel determines whether the kernel promotes clustering. Also, we found that Allee effects promote clustering, whereas gene flow can have a counteracting influence. In line with earlier findings, we could demonstrate that phenotypic clustering was favored by gradients of intermediate slope.     

An intraspecific comparative analysis of character divergence between sympatric species

Although sympatric character divergence between closely related species has been described in a wide variety of taxa, the evolutionary processes responsible for generating these patterns are difficult to identify. One hypothesis that can explain sympatric differences is ecological character displacement: the sympatric origin of morphologically divergent phenotypes in response to selection caused by interspecific competition. Alternatively, populations may adapt to different conditions in allopatry, with sympatric distributions evolving through selective colonization and proliferation of ecologically compatible phenotypes. In this study, I characterize geographic variation within two sibling species of rocky-shore gastropods that have partially overlapping distributions in central California. In sympatry, both Nucella emarginata and N. ostrina show significant differences in shell shape and shell ornamentation that together suggest that where the two species co-exist, divergent phenotypes arose as an evolutionary consequence of competition. To examine the evolutionary origins of divergent characters in sympatry, I used a comparative method based on spatial autocorrelation to remove the portion of the phenotypic variance among populations that is explained by genetic distance (using mitochondrial DNA sequences and allozyme frequency data). Because the remaining portion of the phenotypic variance represents the independent divergence of individual populations, a significant sympatric difference in the corrected dataset provides evidence of true character displacement: significant sympatric character evolution that is independent of population history. After removal of genetic distance effects in Nucella, shell shape differences remain statistically significant in N. emarginata, providing evidence of significant sympatric character divergence. However, for external shell ornamentation in both species and shell shape in N. ostrina, the significance of sympatric differences is lost in the corrected dataset, indicating that colonization events and gene flow have played important roles in the evolutionary history of character divergence in sympatry. Although the absence of a widely dispersing planktonic larva in the life cycle of Nucella will promote local adaptation, the results here indicate that once advantageous traits arise, demographic processes, such as recurrent gene flow between established populations and extinction and recolonization, are important factors contributing to the geographic pattern of sympatric character divergence.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

The evolution of conditional dispersal and reproductive isolation along environmental gradients

Dispersal modulates gene flow throughout a population's spatial range. Gene flow affects adaptation at local spatial scales, and consequently impacts the evolution of reproductive isolation. A recent theoretical investigation has demonstrated that local adaptation along an environmental gradient, facilitated by the evolution of limited dispersal, can lead to parapatric speciation even in the absence of assortative mating. This and other studies assumed unconditional dispersal, so individuals start dispersing without regard to local environmental conditions. However, many species disperse conditionally; their propensity to disperse is contingent upon environmental cues, such as the degree of local crowding or the availability of suitable mates. Here, we use an individual-based model in continuous space to investigate by numerical simulation the relationship between the evolution of threshold-based conditional dispersal and parapatric speciation driven by frequency-dependent competition along environmental gradients. We find that, as with unconditional dispersal, parapatric speciation occurs under a broad range of conditions when reproduction is asexual, and under a more restricted range of conditions when reproduction is sexual. In both the asexual and sexual cases, the evolution of conditional dispersal is strongly influenced by the slope of the environmental gradient: shallow environmental gradients result in low dispersal thresholds and high dispersal distances, while steep environmental gradients result in high dispersal thresholds and low dispersal distances. The latter, however, remain higher than under unconditional dispersal, thus undermining isolation by distance, and hindering speciation in sexual populations. Consequently, the speciation of sexual populations under conditional dispersal is triggered by a steeper gradient than under unconditional dispersal. Enhancing the disruptiveness of frequency-dependent selection, more box-shaped competition kernels dramatically lower the speciation-enabling slope of the environmental gradient.     

Weak disruptive selection and incomplete phenotypic divergence in two classic examples of sympatric speciation: cameroon crater lake cichlids

Abstract Recent documentation of a few compelling examples of sympatric speciation led to a proliferation of theoretical models. Unfortunately, plausible examples from nature have rarely been used to test model predictions, such as the initial presence of strong disruptive selection. Here I estimated the form and strength of selection in two classic examples of sympatric speciation: radiations of Cameroon cichlids restricted to Lakes Barombi Mbo and Ejagham. I measured five functional traits and relative growth rates in over 500 individuals within incipient species complexes from each lake. Disruptive selection was prevalent in both groups on single and multivariate trait axes but weak relative to stabilizing selection on other traits and most published estimates of disruptive selection. Furthermore, despite genetic structure, assortative mating, and bimodal species-diagnostic coloration, trait distributions were unimodal in both species complexes, indicating the earliest stages of speciation. Long waiting times or incomplete sympatric speciation may result when disruptive selection is initially weak. Alternatively, I present evidence of additional constraints in both species complexes, including weak linkage between coloration and morphology, reduced morphological variance aligned with nonlinear selection surfaces, and minimal ecological divergence. While other species within these radiations show complete phenotypic separation, morphological and ecological divergence in these species complexes may be slow or incomplete outside optimal parameter ranges, in contrast to rapid divergence of their sexual coloration.     

Sexual selection at the protein level drives the extraordinary divergence of sex-related genes during sympatric speciation

An increasing number of molecular studies are indicating that, in a wide variety of species, genes directly related to fertilization evolve at extraordinarily high rates. We try to gain insight into the dynamics of this rapid evolution and its underlying mechanisms by means of a simple theoretical model. In the model, sexual selection and sympatric speciation act together in order to drive rapid divergence of gamete recognition proteins. In this process, intraspecific competition for fertilizations enlarges male gamete protein variation by means of evolutionary branching, which initiates sympatric speciation. In addition, avoidance of competition for fertilizations between the incipient species drives the rapid evolution of gamete recognition proteins. This mechanism can account for both strong stabilizing selection on gamete recognition proteins within species and rapid divergence between species. Moreover, it can explain the empirical finding that the rate of divergence of fertilization genes is not constant, but highest between closely related species.     

Niche diversity in crustacean cryptic species: complementarity in spatial distribution and predation risk

Recent genetic studies indicate that species with very close phenotypic similarity (“cryptic species”) are a common feature of nature, and that such cryptic species often coexist in communities. Because traditional views of species coexistence demand that species differ in phenotype to coexist stably, the existence of sympatric cryptic species appears to challenge traditional perspectives of coexistence. We evaluated niche diversity in three recently discovered species of Hyalella amphipods that occur sympatrically in lakes and share close phenotypic similarity. We found that, in some cases, these species exhibited strong complementary spatial distributions within the littoral zone of lakes, both across a distance-from-shore gradient, and a vertical depth gradient. Additionally, we compared fish stomach contents with habitat samples and found that species differed in their vulnerability to predation from sunfish (Lepomis spp.). Complementarity among species across axes of spatial distribution and predation risk, two important niche components, suggests that species with close phenotypic similarity may differ appreciably along ecologically relevant axes. Our results, considered in the light of previous studies, suggest a community structured by predator-mediated coexistence or sequential dominance across environmental gradients in the littoral zone.     

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (Q_ST) and neutral genetic divergence (F_ST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise Q_STs–F_STs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

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.     

Natural selection drives altitudinal divergence at the albumin locus in deer mice, Peromyscus maniculatus

In populations that are distributed across steep environmental gradients, the potential for local adaptation is largely determined by the spatial scale of fitness variation relative to dispersal distance. Since altitudinal gradients are generally characterized by dramatic ecological transitions over relatively short linear distances, adaptive divergence across such gradients will typically require especially strong selection to counterbalance the homogenizing effect of gene flow. Here we report the results of a study that was designed to test for evidence of adaptive divergence across an altitudinal gradient in a natural population of deer mice, Peromyscus maniculatus. We conducted a multilocus survey of allozyme variation across a steep altitudinal gradient in the southern Rocky Mountains that spanned several distinct biomes, from prairie grassland to alpine tundra. As a control for the effects of altitude, we also surveyed the same loci in mice sampled along a latitudinal transect through the prairie grassland that ran perpendicular to the east-west altitudinal transect. We used a coalescent-based simulation model to identify loci that deviated from neutral expectations, and we then assessed whether locus-specific patterns of variation were nonrandom with respect to altitude. Results indicated that the albumin locus (Alb) reflects a history of diversifying selection across the altitudinal gradient. This conclusion is supported by two main lines of evidence: (1) Alb was characterized by levels of divergence across the altitudinal transect that exceeded neutral expectations in two consecutive years of sampling (in contrast to the spatial pattern of variation across the latitudinal transect), and (2) levels of divergence at the Alb locus exhibited a positive association with altitudinal distance in both years (in contrast to the pattern observed at unlinked loci). We conclude that clinal variation at the Alb locus reflects a balance between gene flow and diversifying selection that results from elevational changes in fitness rankings among alternative genotypes.     

Species range expansion by beneficial mutations

A species' range can be limited when there is no genetic variation for a trait that allows for adaptation to more extreme environments. We study how range expansion occurs by the establishment of a new mutation that affects a quantitative trait in a spatially continuous population. The optimal phenotype for the trait varies linearly in space. The survival probabilities of new mutations affecting the trait are found by simulation. Shallow environmental gradients favour mutations that arise nearer to the range margin and that have smaller phenotypic effects than do steep gradients. Mutations that become established in shallow environmental gradients typically result in proportionally larger range expansions than those that establish in steep gradients. Mutations that become established in populations with high maximum growth rates tend to originate nearer to the range edge and to cause relatively smaller range expansion than mutations that establish in populations with low maximum growth rates. Under plausible parameter values, mutations that allow for range expansion tend to have large phenotypic effects (more than one phenotypic standard deviation) and cause substantial range expansions (15% or more). Sexual reproduction allows for larger range expansions and adaptation to more extreme environments than asexual reproduction.