Comparing the consequences of natural selection,adaptive phenotypic plasticity,and matching habitat choice for phenotype–environment matching,population genetic structure,and reproductive isolation in meta‐populations

Organisms commonly experience significant spatiotemporal variation in their environments. In response to such heterogeneity, different mechanisms may act that enhance ecological performance locally. However, depending on the nature of the mechanism involved, the consequences for populations may differ greatly. Building on a previous model that investigated the conditions under which different adaptive mechanisms (co)evolve, this study compares the ecological and evolutionary population consequences of three very different responses to environmental heterogeneity: matching habitat choice (directed gene flow), adaptive plasticity (associated with random gene flow), and divergent natural selection. Using individual‐based simulations, we show that matching habitat choice can have a greater adaptive potential than plasticity or natural selection: it allows for local adaptation while protecting genetic polymorphism despite global mating or strong environmental changes. Our simulations further reveal that increasing environmental fluctuations and unpredictability generally favor the emergence of specialist genotypes but that matching habitat choice is better at preventing local maladaptation by individuals. This confirms that matching habitat choice can speed up the genetic divergence among populations, cause indirect assortative mating via spatial clustering, and hence even facilitate sympatric speciation. This study highlights the potential importance of directed dispersal in local adaptation and speciation, stresses the difficulty of deriving its operation from nonexperimental observational data alone, and helps define a set of ecological conditions which should favor its emergence and subsequent detection in nature.     

Matching habitat choice promotes species persistence under climate change

Species may survive under contemporary climate change by either shifting their range or adapting locally to the warmer conditions. Theoretical and empirical studies recently underlined that dispersal, the central mechanism behind these responses, may depend on the match between an individuals’ phenotype and local environment. Such matching habitat choice is expected to induce an adaptive gene flow, but it now remains to be studied whether this local process could promote species’ responses to climate change. Here, we investigate this by developing an individual‐based model including either random dispersal or temperature‐dependent matching habitat choice. We monitored population composition and distribution through space and time under climate change. Relative to random dispersal, matching habitat choice induced an adaptive gene flow that lessened spatial range loss during climate warming by improving populations’ viability within the range (i.e. limiting range fragmentation) and by facilitating colonization of new habitats at the cold margin. The model even predicted range contraction under random dispersal but range expansion under optimal matching habitat choice. These benefits of matching habitat choice for population persistence mostly resulted from adaptive immigration decision and were greater for populations with larger dispersal distance and higher emigration probability. We also found that environmental stochasticity resulted in suboptimal matching habitat choice, decreasing the benefits of this dispersal mode under climate change. However population persistence was still better under suboptimal matching habitat choice than under random dispersal. Our results highlight the urgent need to implement more realistic mechanisms of dispersal such as matching habitat choice into models predicting the impacts of ongoing climate change on biodiversity.     

Conceptual issues in local adaptation

Studies of local adaptation provide important insights into the power of natural selection relative to gene flow and other evolutionary forces. They are a paradigm for testing evolutionary hypotheses about traits favoured by particular environmental factors. This paper is an attempt to summarize the conceptual framework for local adaptation studies. We first review theoretical work relevant for local adaptation. Then we discuss reciprocal transplant and common garden experiments designed to detect local adaptation in the pattern of deme × habitat interaction for fitness. Finally, we review research questions and approaches to studying the processes of local adaptation – divergent natural selection, dispersal and gene flow, and other processes affecting adaptive differentiation of local demes. We advocate multifaceted approaches to the study of local adaptation, and stress the need for experiments explicitly addressing hypotheses about the role of particular ecological and genetic factors that promote or hinder local adaptation. Experimental evolution of replicated populations in controlled spatially heterogeneous environments allow direct tests of such hypotheses, and thus would be a valuable way to complement research on natural populations.     

Local adaptation and ecological differentiation under selection,migration, and drift in Arabidopsis lyrata*

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole‐genome‐based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4N_em < 11) and have small effective population sizes (N_e < 6000), suggesting that both migration and drift can counteract local selection. However, the three‐year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.     

Mating timing,dispersal and local adaptation in patchy environments

Dispersal is a life‐history trait that can evolve under various known selective pressures as identified by a multitude of theoretical and empirical studies. Yet only few of them are considering the succession of mating and dispersal. The sequence of these events influences gene flow and consequently affects the dynamics and evolution of populations. We use individual‐based simulations to investigate the evolution of the timing of dispersal and mating, i.e. mating before or after dispersal. We assume a discrete insect metapopulation in a heterogeneous environment, where populations may adapt to local conditions and only females are allowed to disperse. We run the model assuming different levels of species habitat tolerance, carrying capacity, and temporal environmental variability. Our results show that in species with narrow habitat tolerance, low to moderate dispersal evolves in combination with mating after dispersal (post‐dispersal mating). With such a strategy dispersing females benefit from mating with a resident male, as their offspring will be better adapted to the local habitat conditions. On the contrary, in species with wide habitat tolerance higher dispersal rates in combination with pre‐dispersal mating evolves. In this case individuals are adapted to the ‘average’ habitat where pre‐dispersal mating conveys the benefit of carrying relatives’ genes into a new population. With high dispersal rates and large population size, local adaptation and kin structure both vanish and the temporal sequence of dispersal and mating may become a (nearly) neutral trait.     

Local specialization and maladaptation in the Mediterranean blue tit (Parus caeruleus)

Although the evolutionary interplay between gene flow and local adaptation of organisms in heterogeneous environments has been widely discussed from a theoretical point of view, few empirical studies have been designed to test predictions on the consequences of habitat patchiness on the evolution of life history traits. Using blue tits in Mediterranean habitat mosaics as a model, we defined two nested levels of habitat heterogeneity: an inter-regional level which compares two isolated landscapes (mainland, southern France vs the island of Corsica), and an intra-regional level which compares two habitat types within each landscape (deciduous vs evergreen trees). Deciduous habitats are more common than evergreen habitats on the mainland whereas the opposite is true on the island. Results suggest that: (1) on a regional scale, each population is specialized to the more common habitat, i.e. life history traits have evolved in such a way that breeding success is high; (2) in the less common habitats within each landscape, birds are clearly mistimed because they mismatch the best period of food availability, and hence their breeding success is lower; and (3) the density of the populations and the morphometry of the birds support the model of ideal despotic distribution. These results, which are supported by preliminary data on the genetic variation and gene flow of populations in the mainland landscape, are consistent with a source-sink model of population structure within each landscape. They are discussed in the framework of metapopulation theory and habitat selection models.     

A multilocus perspective on colonization accompanied by selection and gene flow

The colonization of novel habitats involves complex interactions between founder events, selection, and ongoing migration, and can lead to diverse evolutionary outcomes from local extinction to adaptation to speciation. Although there have been several studies of the demography of colonization of remote habitats, less is known about the demographic consequences of colonization of novel habitats within a continuous species range. Populations of the Eastern Fence Lizard, Sceloporus undulatus, are continuously distributed across two dramatic transitions in substrate color in southern New Mexico and have undergone rapid adaptation following colonization of these novel environments. Blanched forms inhabit the gypsum sand dunes of White Sands and melanic forms are found on the black basalt rocks of the Carrizozo lava flow. Each of these habitats formed within the last 10,000 years, allowing comparison of genetic signatures of population history for two independent colonizations from the same source population. We present evidence on phenotypic variation in lizard color, environmental variation in substrate color, and sequence variation for mitochondrial DNA and 19 independent nuclear loci. To confirm the influence of natural selection and gene flow in this system, we show that phenotypic variation is best explained by environmental variation and that neutral genetic variation is related to distance between populations, not partitioned by habitat. The historical demography of colonization was inferred using an Approximate Bayesian Computation (ABC) framework that incorporates known geological information and allows for ongoing migration with the source population. The inferences differed somewhat between mtDNA and nuclear markers, but overall provided strong evidence of historical size reductions in both white sand and black lava populations at the time of colonization. Populations in both novel habitats appear to have undergone partial but incomplete recovery from the initial bottleneck. Both ABC analyses and measures of mtDNA sequence diversity also suggested that population reductions were more severe in the black lava compared to the white sands habitat. Differences observed between habitats may be explained by differences in colonization time, habitat geometry, and strength or response to natural selection for substrate matching. Finally, effective population size reductions in this system appear to be more dramatic when colonization is accompanied by a change in selection regime. Our analyses are consistent with a demographic cost of adaptation to novel environments and show that it is possible to infer aspects of the historical demography of local adaptation even in the presence of ongoing gene flow.     

Phenotype‐dependent selection underlies patterns of sorting across habitats: the case of stream‐fishes

Spatial and temporal heterogeneity within landscapes influences the distribution and phenotypic diversity of individuals both within and across populations. Phenotype–habitat correlations arise either through phenotypes within an environment altering through the process of natural selection or plasticity, or phenotypes remaining constant but individuals altering their distribution across environments. The mechanisms of non‐random movement and phenotype‐dependent habitat choice may account for associations within highly heterogeneous systems, such as streams, where local adaptation may be negated, plasticity too costly and movement is particularly important. Despite growing attention, however, few empirical tests have yet to be conducted. Here we provide a test of phenotype‐dependent habitat choice and ask: 1) if individuals collected from a single habitat type continue to select original habitat; 2) if decisions are phenotype‐dependent and functionally related to habitat requirements; and 3) if phenotypic‐sorting continues despite increasing population density. To do so we both conducted experimental trials manipulating the density of four stream‐fish species collected from either a single riffle or pool and developed a game‐theoretical model exploring the influence of individuals’ growth rate, sampling and competitive abilities as well as interference on distribution across two habitats as a function of density. Our experimental trials show individuals selecting original versus alternative habitats differed in their morphologies, that morphologies were functionally related to habitat‐type swimming demands, and that phenotypic‐sorting remained significant (although decreased) as density increased. According to our model this only occurs when phenotypes have contrasting habitat preferences and only one phenotype disperses (i.e. selects alternatives) in response to density pressures. This supports our explanation that empirical habitat selection was due to a combination of collecting a fraction of mobile individuals with different habitat preferences and the exclusion of individuals via scramble competition at increased densities. Phenotype‐dependent habitat choice can thereby account for observed patterns of natural stream‐fish distribution.     

GENE FLOW VERSUS LOCAL ADAPTATION IN THE NORTHERN ACORN BARNACLE, SEMIBALANUS BALANOIDES: INSIGHTS FROM MITOCHONDRIAL DNA VARIATION

Abstract In reciprocal transplant experiments, Bertness and Gaines (1993) found that Semibalanus balanoides juveniles that had settled in an upper Narragansett Bay estuary survived better in that estuary that did juveniles from coastal localities. The observed pattern of survivorship led to the claim that local adaptation may result from a combination of limited gene flow between and strong selection within these habitats. Here we test the hypothesis that limited gene flow has led to habitat‐specific population differentiation using sequence and restriction fragment length polymorphism analyses of the mitochondrial DNA D‐loop region of S. balanoides. Samples were analyzed from replicated coastal and estuary localities in both Narragansett Bay, Rhode Island, and Damariscotta River, Maine. The patterns of F_ST indicate that gene flow between coast and estuary is extensive (Nm > 100) and is not lower in the estuary with lower flushing rates (Narragansett Bay). Given the high estimate of genetic exchange, adaptations for unpredictable environments seem more likely than local adaptation in this species because loci that respond to selection in one generation are essentially homogenized by the next seasons' settlement. Nevertheless, these estimates of neutral gene flow can help identify the strength of selection necessary for local adaptation to accumulate in Semibalanus.     

LOCAL ADAPTATION WHEN COMPETITION DEPENDS ON PHENOTYPIC SIMILARITY

Recent work incorporating demographic–genetic interactions indicates the importance of population size, gene flow, and selection in influencing local adaptation. This work typically assumes that density‐dependent survival affects individuals equally, but individuals in natural population rarely compete equally. Among‐individual differences in resource use generate stronger competition between more similar phenotypes (frequency‐dependent competition) but it remains unclear how this additional form of selection changes the interactions between population size, gene flow, and local stabilizing selection. Here, we integrate migration–selection dynamics with frequency‐dependent competition. We developed a coupled demographic‐quantitative genetic model consisting of two patches connected by dispersal and subject to local stabilizing selection and competition. Our model shows that frequency‐dependent competition slightly increases local adaptation, greatly increases genetic variance within patches, and reduces the amount that migration depresses population size, despite the increased genetic variance load. The effects of frequency‐dependence depend on the strength of divergent selection, trait heritability, and when mortality occurs in the life cycle in relation to migration and reproduction. Essentially, frequency‐dependent competition reduces the density‐dependent interactions between migrants and residents, the extent to which depends on how different and common immigrants are compared to residents. Our results add new dynamics that illustrate how competition can alter the effects of gene flow and divergent selection on local adaptation and population carrying capacities.     

DEMOGRAPHIC SOURCE-SINK DYNAMICS RESTRICT LOCAL ADAPTATION IN ELLIOTT'S BLUEBERRY (VACCINIUM ELLIOTTII)

In heterogeneous landscapes, divergent selection can favor the evolution of locally adapted ecotypes, especially when interhabitat gene flow is minimal. However, if habitats differ in size or quality, source-sink dynamics can shape evolutionary trajectories. Upland and bottomland forests of the southeastern USA differ in water table depth, light availability, edaphic conditions, and plant community. We conducted a multiyear reciprocal transplant experiment to test whether Elliott's blueberry ( Vaccinium elliottii ) is locally adapted to these contrasting environments. Additionally, we exposed seedlings and cuttings to prolonged drought and flooding in the greenhouse to assess fitness responses to abiotic stress. Contrary to predictions of local adaptation, V. elliottii families exhibited significantly higher survivorship and growth in upland than in bottomland forests and under drought than flooded conditions, regardless of habitat of origin. Neutral population differentiation was minimal, suggesting widespread interhabitat migration. Population density, reproductive output, and genetic diversity were all significantly greater in uplands than in bottomlands. These disparities likely result in asymmetric gene flow from uplands to bottomlands. Thus, adaptation to a marginal habitat can be constrained by small populations, limited fitness, and immigration from a benign habitat. Our study highlights the importance of demography and genetic diversity in the evolution of local (mal)adaptation.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Local adaptation and matching habitat choice in female barn owls with respect to melanic coloration

Local adaptation is a major mechanism underlying the maintenance of phenotypic variation in spatially heterogeneous environments. In the barn owl (Tyto alba), dark and pale reddish-pheomelanic individuals are adapted to conditions prevailing in northern and southern Europe, respectively. Using a long-term dataset from Central Europe, we report results consistent with the hypothesis that the different pheomelanic phenotypes are adapted to specific local conditions in females, but not in males. Compared to whitish females, reddish females bred in sites surrounded by more arable fields and less forests. Colour-dependent habitat choice was apparently beneficial. First, whitish females produced more fledglings when breeding in wooded areas, whereas reddish females when breeding in sites with more arable fields. Second, cross-fostering experiments showed that female nestlings grew wings more rapidly when both their foster and biological mothers were of similar colour. The latter result suggests that mothers should particularly produce daughters in environments that best match their own coloration. Accordingly, whiter females produced fewer daughters in territories with more arable fields. In conclusion, females displaying alternative melanic phenotypes bred in habitats providing them with the highest fitness benefits. Although small in magnitude, matching habitat selection and local adaptation may help maintain variation in pheomelanin coloration in the barn owl.     

Quantifying the effects of migration and mutation on adaptation and demography in spatially heterogeneous environments

How do mutation and gene flow influence population persistence, niche expansion and local adaptation in spatially heterogeneous environments? In this article, we analyse a demographic and evolutionary model of adaptation to an environment containing two habitats in equal frequencies, and we bridge the gap between different theoretical frameworks. Qualitatively, our model yields four qualitative types of outcomes: (i) global extinction of the population, (ii) adaptation to one habitat only, but also adaptation to both habitats with, (iii) specialized phenotypes or (iv) with generalized phenotypes, and we determine the conditions under which each equilibrium is reached. We derive new analytical approximations for the local densities and the distributions of traits in each habitat under a migration–selection–mutation balance, compute the equilibrium values of the means, variances and asymmetries of the local distributions of phenotypes, and contrast the effects of migration and mutation on the evolutionary outcome. We then check our analytical results by solving our model numerically, and also assess their robustness in the presence of demographic stochasticity. Although increased migration results in a decrease in local adaptation, mutation in our model does not influence the values of the local mean traits. Yet, both migration and mutation can have dramatic effects on population size and even lead to metapopulation extinction when selection is strong. Niche expansion, the ability for the population to adapt to both habitats, can also be prevented by small migration rates and a reduced evolutionary potential characterized by rare mutation events of small effects; however, niche expansion is otherwise the most likely outcome. Although our results are derived under the assumption of clonal reproduction, we finally show and discuss the links between our model and previous quantitative genetics models.     

A Meta-Analysis of Local Adaptation in Plants

Local adaptation is of fundamental importance in evolutionary, population, conservation, and global-change biology. The generality of local adaptation in plants and whether and how it is influenced by specific species, population and habitat characteristics have, however, not been quantitatively reviewed. Therefore, we examined published data on the outcomes of reciprocal transplant experiments using two approaches. We conducted a meta-analysis to compare the performance of local and foreign plants at all transplant sites. In addition, we analysed frequencies of pairs of plant origin to examine whether local plants perform better than foreign plants at both compared transplant sites. In both approaches, we also examined the effects of population size, and of the habitat and species characteristics that are predicted to affect local adaptation. We show that, overall, local plants performed significantly better than foreign plants at their site of origin: this was found to be the case in 71.0% of the studied sites. However, local plants performed better than foreign plants at both sites of a pair-wise comparison (strict definition of local adaption) only in 45.3% of the 1032 compared population pairs. Furthermore, we found local adaptation much more common for large plant populations (>1000 flowering individuals) than for small populations (<1000 flowering individuals) for which local adaptation was very rare. The degree of local adaptation was independent of plant life history, spatial or temporal habitat heterogeneity, and geographic scale. Our results suggest that local adaptation is less common in plant populations than generally assumed. Moreover, our findings reinforce the fundamental importance of population size for evolutionary theory. The clear role of population size for the ability to evolve local adaptation raises considerable doubt on the ability of small plant populations to cope with changing environments.     

A genomic selection component analysis characterizes migration‐selection balance

The genetic differentiation of populations in response to local selection pressures has long been studied by evolutionary biologists, but key details about the process remain obscure. How rapidly can local adaptation evolve, how extensive is the process across the genome, and how strong are the opposing forces of natural selection and gene flow? Here, we combine direct measurement of survival and reproduction with whole‐genome genotyping of a plant species (Mimulus guttatus) that has recently invaded a novel habitat (the Quarry population). We renovate the classic selection component method to accommodate genomic data and observe selection at SNPs throughout the genome. SNPs showing viability selection in Quarry exhibit elevated divergence from neighboring populations relative to neutral SNPs. We also find that nonsignificant SNPs exhibit a subtle, but still significant, change in allele frequency toward neighboring populations, a predicted effect of gene flow. Given that the Quarry population is most probably only 30–40 generations old, the alleles conferring local advantage are almost certainly older than the population itself. Thus, local adaptation owes to the recruitment of standing genetic variation.     

Developments of the Price equation and natural selection under uncertainty

Many approaches to the study of adaptation, following Darwin, centre on the number of offspring of individuals. Population genetics theory makes clear that predicting gene frequency changes requires more detailed knowledge, for example of linkage and linkage disequilibrium and mating systems. Because gene frequency changes underlie adaptation, this can lead to a suspicion that approaches ignoring these sophistications are approximate or tentative or wrong. Stochastic environments and sexual selection are two topics in which there are widespread views that focusing on number of offspring of individuals is not enough, and that proper treatments require the introduction of further details, namely variability in offspring number and linkage disequilibrium, respectively. However, the bulk of empirical research on adaptation and a great deal of theoretical work continue to employ these approaches. Here, a new theoretical development arising from the Price equation provides a formal justification in very general circumstances for focusing on the arithmetic average of the relative number of offspring of individuals.     

Role of selection and gene flow in population differentiation at the edge vs. interior of the species range differing in climatic conditions

Evaluating the relative importance of neutral and adaptive processes as determinants of population differentiation across environments is a central theme of evolutionary biology. We applied the Q_ST–F_ST comparison flanked by a direct test for local adaptation to infer the role of climate‐driven selection and gene flow in population differentiation of an annual grass Avena sterilis in two distinct parts of the species range, edge and interior, which represent two globally different climates, desert and Mediterranean. In a multiyear reciprocal transplant experiment, the plants of desert and Mediterranean origin demonstrated home advantage, and population differentiation in several phenotypic traits related to reproduction exceeded neutral predictions, as determined by comparisons of Q_ST values with theoretical F_ST distributions. Thus, variation in these traits likely resulted from local adaptation to desert and Mediterranean environments. The two separate common garden experiments conducted with different experimental design revealed that two population comparisons, in contrast to multi‐population comparisons, are likely to detect population differences in virtually every trait, but many of these differences reflect effects of local rather than regional environment. We detected a general reduction in neutral (SSR) genetic variation but not in adaptive quantitative trait variation in peripheral desert as compared with Mediterranean core populations. On the other hand, the molecular data indicated intensive gene flow from the Mediterranean core towards desert periphery. Although species range position in our study (edge vs. interior) was confounded with climate (desert vs. Mediterranean), the results suggest that the gene flow from the species core does not have negative consequences for either performance of the peripheral plants or their adaptive potential.     

Carryover effects from natal habitat type upon competitive ability lead to trait divergence or source–sink dynamics

Local adaptation to rare habitats is difficult due to gene flow, but can occur if the habitat has higher productivity. Differences in offspring phenotypes have attracted little attention in this context. We model a scenario where the rarer habitat improves offspring's later competitive ability – a carryover effect that operates on top of local adaptation to one or the other habitat type. Assuming localised dispersal, so the offspring tend to settle in similar habitat to the natal type, the superior competitive ability of offspring remaining in the rarer habitat hampers immigration from the majority habitat. This initiates a positive feedback between local adaptation and trait divergence, which can thereafter be reinforced by coevolution with dispersal traits that match ecotype to habitat type. Rarity strengthens selection on dispersal traits and promotes linkage disequilibrium between locally adapted traits and ecotype‐habitat matching dispersal. We propose that carryover effects may initiate isolation by ecology.     

MIGRATION-SELECTION BALANCE AND LOCAL ADAPTATION OF MITOCHONDRIAL HAPLOTYPES IN RUFOUS-COLLARED SPARROWS (ZONOTRICHIA CAPENSIS) ALONG AN ELEVATIONAL GRADIENT

Variable selection pressures across heterogeneous landscapes can lead to local adaptation of populations. The extent of local adaptation depends on the interplay between natural selection and gene flow, but the nature of this relationship is complex. Gene flow can constrain local adaptation by eroding differentiation driven by natural selection, or local adaptation can itself constrain gene flow through selection against maladapted immigrants. Here we test for evidence that natural selection constrains gene flow among populations of a widespread passerine bird ( Zonotrichia capensis ) that are distributed along an elevational gradient in the Peruvian Andes. Using multilocus sequences and microsatellites screened in 142 individuals collected along a series of replicate transects, we found that mitochondrial gene flow was significantly reduced along elevational transects relative to latitudinal control transects. Nuclear gene flow, however, was not similarly reduced. Clines in mitochondrial haplotype frequency were strongly associated with transitions in environmental variables along the elevational transects, but this association was not observed for the nuclear markers. These results suggest that natural selection constrains mitochondrial gene flow along elevational gradients and that the mitonuclear discrepancy may be due to local adaptation of mitochondrial haplotypes.