The genetic architecture of adaptation under migration-selection balance

Many ecologically important traits have a complex genetic basis, with the potential for mutations at many different genes to shape the phenotype. Even so, studies of local adaptation in heterogeneous environments sometimes find that just a few quantitative trait loci (QTL) of large effect can explain a large percentage of observed differences between phenotypically divergent populations. As high levels of gene flow can swamp divergence at weakly selected alleles, migration-selection-drift balance may play an important role in shaping the genetic architecture of local adaptation. Here, we use analytical approximations and individual-based simulations to explore how genetic architecture evolves when two populations connected by migration experience stabilizing selection toward different optima. In contrast to the exponential distribution of allele effect sizes expected under adaptation without migration (Orr 1998), we find that adaptation with migration tends to result in concentrated genetic architectures with fewer, larger, and more tightly linked divergent alleles. Even if many small alleles contribute to adaptation at the outset, they tend to be replaced by a few large alleles under prolonged bouts of stabilizing selection with migration. All else being equal, we also find that stronger selection can maintain linked clusters of locally adapted alleles over much greater map distances than weaker selection. The common empirical finding of QTL of large effect is shown to be expected with migration in a heterogeneous landscape, and these QTL may often be composed of several tightly linked alleles of smaller effect.     

The magnitude of local adaptation under genotype‐dependent dispersal

Dispersal moves individuals from patches where their immediate ancestors were successful to sites where their genotypes are untested. As a result, dispersal generally reduces fitness, a phenomenon known as “migration load.” The strength of migration load depends on the pattern of dispersal and can be dramatically lessened or reversed when individuals move preferentially toward patches conferring higher fitness. Evolutionary ecologists have long modeled nonrandom dispersal, focusing primarily on its effects on population density over space, the maintenance of genetic variation, and reproductive isolation. Here, we build upon previous work by calculating how the extent of local adaptation and the migration load are affected when individuals differ in their dispersal rate in a genotype‐dependent manner that alters their match to their environment. Examining a one‐locus, two‐patch model, we show that local adaptation occurs through a combination of natural selection and adaptive dispersal. For a substantial portion of parameter space, adaptive dispersal can be the predominant force generating local adaptation. Furthermore, genetic load may be largely averted with adaptive dispersal whenever individuals move before selective deaths occur. Thus, to understand the mechanisms driving local adaptation, biologists must account for the extent and nature of nonrandom, genotype‐dependent dispersal, and the potential for adaptation via spatial sorting of genotypes.     

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.     

Genetic variability,local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South‐East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South‐East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. F_ST and principal components analysis revealed striking divergence between South‐East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South‐East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South‐East Asia during mid‐Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the ‘genomic islands’ scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid‐Pleistocene.     

Migration-induced phenotypic divergence: the migration-selection balance of correlated traits

Genetically correlated traits are known to respond to indirect selection pressures caused by directional selection on other traits. It is however unclear how local adaptation in populations diverging along some phenotypic traits but not others is affected by the joint action of gene flow and genetic correlations among traits. This simulation study shows that although gene flow is a potent constraining mechanism of population adaptive divergence, it may induce phenotypic divergence in traits under homogeneous selection among habitats if they are genetically correlated with traits under divergent selection. This correlated phenotypic divergence is a nonmonotonous function of migration and increases with mutational correlation among traits. It also increases with the number of divergently selected traits provided their genetic autonomy relative to the uniformly selected trait is reduced by specific patterns of genetic covariances: populations with lower effective trait dimensionality are more likely to generate very large correlated divergence. The correlated divergence is likely to be picked up by Q(ST)-F(ST) analysis of population genetic differentiation and be erroneously ascribed to adaptive divergence under divergent selection. This study emphasizes the necessity to understand the interaction between selection and the genetic basis of adaptation in a multivariate rather than univariate context.     

Eco-evolutionary metapopulation dynamics and the spatial scale of adaptation

We construct a model that combines extinction-colonization dynamics with the dynamics of local adaptation in a network of habitat patches of dissimilar qualities. We derive a deterministic approximation for the stochastic model that allows the calculation of patch-specific incidences of occupancy and levels of adaptation at steady state. Depending on (i) the strength of local selection, (ii) the amount of genetic variance, (iii) the demographic cost of maladaptation, (iv) the spatial scale of gene flow, and (v) the amount of habitat heterogeneity, the model predicts adaptation at different spatial scales. Local adaptation is predicted when there is much genetic variance and strong selection, while network-level adaptation occurs when the demographic cost of maladaptation is low. For little genetic variance and high cost of maladaptation, the model predicts network-level habitat specialization in species with long-range migration but an intermediate scale of adaptation (mosaic specialization) in species with short-range migration. In fragmented landscapes, the evolutionary dynamics of adaptation may both decrease and enhance metapopulation viability in comparison with no evolution. The model can be applied to real patch networks with given sizes, qualities, and spatial positions of habitat patches.     

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait’s genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor “concentrated” genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.     

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.     

Cryptic lineage divergence in marine environments: genetic differentiation at multiple spatial and temporal scales in the widespread intertidal goby Gobiosoma bosc

The adaptive radiation of the seven‐spined gobies (Gobiidae: Gobiosomatini) represents a classic example of how ecological specialization and larval retention can drive speciation through local adaptation. However, geographically widespread and phenotypically uniform species also do occur within Gobiosomatini. This lack of phenotypic variation across large geographic areas could be due to recent colonization, widespread gene flow, or stabilizing selection acting across environmental gradients. We use a phylogeographic approach to test these alternative hypotheses in the naked goby Gobiosoma bosc, a widespread and phenotypically invariable intertidal fish found along the Atlantic Coast of North America. Using DNA sequence from 218 individuals sampled at 15 localities, we document marked intraspecific genetic structure in mitochondrial and nuclear genes at three main geographic scales: (i) between Gulf of Mexico and Atlantic Coast, (ii) between the west coast of the Florida peninsula and adjacent Gulf of Mexico across the Apalachicola Bay, and (iii) at local scales of a few hundred kilometers. Clades on either side of Florida diverged about 8 million years ago, whereas some populations along the East Cost show divergent phylogroups that have differentiated within the last 200,000 years. The absence of noticeable phenotypic or ecological differentiation among lineages suggests the role of stabilizing selection on ancestral phenotypes, together with isolation in allopatry due to reduced dispersal and restricted gene flow, as the most likely explanation for their divergence. Haplotype phylogenies and spatial patterns of genetic diversity reveal frequent population bottlenecks followed by rapid population growth, particularly along the Gulf of Mexico. The magnitude of the genetic divergence among intraspecific lineages suggests the existence of cryptic species within Gobiosoma and indicates that modes of speciation can vary among lineages within Gobiidae.     

The genomic footprint of climate adaptation in Chironomus riparius

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool‐Seq data and population genetic modelling. Common‐garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an F_ST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.     

When can ecological speciation be detected with neutral loci?

It is not yet clear under what conditions empirical studies can reliably detect progress toward ecological speciation through the analysis of allelic variation at neutral loci. We use a simulation approach to investigate the range of parameter space under which such detection is, and is not, likely. We specifically test for the conditions under which divergent natural selection can cause a ‘generalized barrier to gene flow’ that is present across the genome. Our individual‐based numerical simulations focus on how population divergence at neutral loci varies in relation to recombination rate with a selected locus, divergent selection on that locus, migration rate and population size. We specifically test whether genetic differences at neutral markers are greater between populations in different environments than between populations in similar environments. We find that this expected signature of ecological speciation can be detected under part of the parameter space, most consistently when divergent selection is strong and migration is intermediate. By contrast, the expected signature of ecological speciation is not reliably detected when divergent selection is weak or migration is low or high. These findings provide insights into the strengths and weaknesses of using neutral markers to infer ecological speciation in natural systems.     

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80‐year period. Screening of >1000 gene‐associated single‐nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.     

The conditions for speciation through intraspecific competition

Abstract It has been shown theoretically that sympatric speciation can occur if intraspecific competition is strong enough to induce disruptive selection. However, the plausibility of the involved processes is under debate, and many questions on the conditions for speciation remain unresolved. For instance, is strong disruptive selection sufficient for speciation? Which roles do genetic architecture and initial composition of the population play? How strong must assortative mating be before a population can split in two? These are some of the issues we address here. We investigate a diploid multilocus model of a quantitative trait that is under frequency‐dependent selection caused by a balance of intraspecific competition and frequency‐independent stabilizing selection. This trait also acts as mating character for assortment. It has been established previously that speciation can occur only if competition is strong enough to induce disruptive selection. We find that speciation becomes more difficult for very strong competition, because then extremely strong assortment is required. Thus, speciation is most likely for intermediate strengths of competition, where it requires strong, but not extremely strong, assortment. For this range of parameters, however, it is not obvious how assortment can evolve from low to high levels, because with moderately strong assortment less genetic variation is maintained than under weak or strong assortment sometimes none at all. In addition to the strength of frequency‐dependent competition and assortative mating, the roles of the number of loci, the distribution of allelic effects, the initial conditions, costs to being choosy, the strength of stabilizing selection, and the particular choice of the fitness function are explored. A multitude of possible evolutionary outcomes is observed, including loss of all genetic variation, splitting in two to five species, as well as very short and extremely long stable limit cycles. On the methodological side, we propose quantitative measures for deciding whether a given distribution reflects two (or more) reproductively isolated clusters.     

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra‐ and interspecific genetic variation at 25 climate‐related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two‐thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. F_ST outlier analysis identified six SNPs in five climate‐related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation‐with‐migration model indicated that gene flow between the species has been more restricted for climate‐related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species‐specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.     

Fin whale MDH‐1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH‐1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH‐1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.     

Local adaptation, patterns of selection, and gene flow in the Californian serpentine sunflower (Helianthus exilis)

The traditional view of the species as the fundamental unit of evolution has been challenged by observations that in heterogeneous environments, gene flow may be too restricted to overcome the effects of local selection. Whether a species evolves as a cohesive unit depends critically on the dynamic balance between homogenizing gene flow among populations and potentially disruptive local adaptation. To examine this evolutionary balance between "global" gene flow and local selection, we studied northern Californian populations of Helianthus exilis, the serpentine sunflower, within a mosaic of contrasting serpentine and nonserpentine areas that differ considerably in soil chemistry and water availability. Local adaptation to riparian and serpentine habitats was studied in Helianthus exilis along with an analysis of gene flow patterns among populations within these habitats. Local adaptation was assessed in H. exilis during 2002 and 2003 using reciprocal transplant experiments at multiple locations within serpentine and riparian habitats. Effects of competition and germination date on the expression of local adaptation were also examined within the reciprocal transplant experiments. Local adaptation was detected in both years at the local site level and at the level of habitat. The analysis of the transplanted populations indicated that the patterns of selection differed considerably between riparian and serpentine sites. Differential survivorship occurred in serpentine habitats, whereas selection on reproductive output predominated in riparian habitats. Local adaptation was expressed only in the absence of competition. Local adaptation in terms of survivorship was most strongly expressed in treatments with delayed seed germination. Microsatellite markers were used to quantify population genetic parameters and examine the patterns of gene flow among sampled populations. Analysis of molecular markers revealed a system of population patches that freely exchange genes with each other. Strong selection seems to maintain ecotypic variation within this endemic sunflower species, while extensive gene flow among populations prevents local speciation between serpentine and riparian ecotypes.     

Competition and evolution along environmental gradients: patterns, boundaries and sympatric divergence

The present study examined how competitive interactions and environmental conditions generate species boundaries and determine species distributions. A spatially explicit, quantitative genetic, two-species competition model was used to manipulate the strengths of competition, gene flow and local adaptation along environmental gradients. This allowed us to assess the long-term persistence of each species and whether the ranges they inhabited had boundaries in space or were unlimited. We found that a species boundary arises along less steep environmental gradients when the strength of stabilizing selection and diversifying selection are similar. We also found that a species boundary may arise along shallow environmental gradients if interspecific competition is more intense than intraspecific, which relaxes previous requirements for steep gradients for generating range limits. We determined an analytical form for the critical environmental gradient as a function of ecological and genetic parameters at which a species boundary is expected to arise by competition. Results suggest an alternative to resource competition as an explanation for phenotypic divergence between sympatric competitors. Competitors sharing a trait that is under stabilizing selection along an environmental gradient may segregate spatially and evolve in different regions, with phenotypic sympatric divergence reflecting the resulting clines. Along various types of environmental gradients, variation in stabilizing selection intensities could lead to contrasting patterns in the distribution of species. For stabilizing selection strengths in accord with field data estimates, this study predicts that the level of sympatric character divergence would be limited along environmental gradients.     

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site‐specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole‐genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree‐based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.     

Convergent evolution and divergent selection: lizards at the White Sands ecotone

Ecological transition zones, where organismal phenotypes result from a delicate balance between selection and migration, highlight the interplay of local adaptation and gene flow. Here, I study the response of an entire species assemblage to natural selection across a common ecotone. Three lizard species, distributed along a dramatic environmental gradient in substrate color, display convergent adaptation of blanched coloration on the gypsum dunes of White Sands National Monument. I investigate the role of gene flow in modulating phenotypic response to selection by quantifying color variation and genetic variation across the ecotone. I find species differences in degree of background matching and in genetic connectivity of populations across the ecotone. Differences among species in phenotypic response to selection scale precisely to levels of genetic isolation. Species with higher levels of gene flow across the ecotone exhibit less dramatic responses to selection. Results also reveal a strong signal of ecologically mediated divergence for White Sands lizards. For all species, phenotypic variation is better explained by habitat similarity than genetic similarity. Convergent evolution of blanched coloration at White Sands clearly reflects the action of strong divergent selection; however, adaptive response appears to be modulated by gene flow and demographic history and can be predicted by divergence-with-gene-flow models.     

Gene flow in the anemone Anthopleura elegantissima limits signatures of local adaptation across an extensive geographic range

Species inhabiting marine environments face a wide range of environmental conditions that vary spatially across several orders of magnitude. The selective pressures that these conditions impose on marine organisms, in combination with potentially high rates of gene flow between distant populations, make it difficult to predict the extent to which these populations can locally adapt. Here, I identify how selection and gene flow influence the population genetic structure of the anemone Anthopleura elegantissima along the Pacific coast of North America. Isolation by distance is the dominant pattern across the range of this species, with a genetic break near Pt. Conception, CA. Furthermore, demographic modelling suggests that this species was historically confined to southerly latitudes before expanding northward. Outlier analyses identify 24 loci under selection (out of ~1,100), but the same analysis on simulated genetic data generated using the most likely demographic model erroneously identified the same number of loci under selection, if not more. Taken together, these results suggest that demographic processes are the dominant force shaping population genetic patterns in A. elegantissima along the Pacific coast of North America. I discuss these patterns in terms of the evolutionary history of A. elegantissima, the potential for local adaptation, and their consequences with respect to interactions with the endosymbiont Breviolum muscatinei across their geographic range.