Spatial phenotypic and genetic structure of threespine stickleback (Gasterosteus aculeatus) in a heterogeneous natural system,Lake Mývatn,Iceland

Eco‐evolutionary responses of natural populations to spatial environmental variation strongly depend on the relative strength of environmental differences/natural selection and dispersal/gene flow. In absence of geographic barriers, as often is the case in lake ecosystems, gene flow is expected to constrain adaptive divergence between environments – favoring phenotypic plasticity or high trait variability. However, if divergent natural selection is sufficiently strong, adaptive divergence can occur in face of gene flow. The extent of divergence is most often studied between two contrasting environments, whereas potential for multimodal divergence is little explored. We investigated phenotypic (body size, defensive structures, and feeding morphology) and genetic (microsatellites) structure in threespine stickleback (Gasterosteus aculeatus) across five habitat types and two basins (North and South) within the geologically young and highly heterogeneous Lake Mývatn, North East Iceland. We found that (1) North basin stickleback were, on average, larger and had relatively longer spines than South basin stickleback, whereas (2) feeding morphology (gill raker number and gill raker gap width) differed among three of five habitat types, and (3) there was only subtle genetic differentiation across the lake. Overall, our results indicate predator and prey mediated phenotypic divergence across multiple habitats in the lake, in face of gene flow.     

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi‐interconnected ponds within a small geographic region (~7.5 km²) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.     

ADAPTIVE POPULATION DIVERGENCE IN CRYPTIC COLOR-PATTERN FOLLOWING A REDUCTION IN GENE FLOW

Adaptive population divergence is often driven by divergent natural selection, but can be constrained by the homogenizing effect of gene flow between populations. Indeed, a common pattern in nature is an inverse correlation between the degree of adaptive phenotypic divergence between populations and levels of gene flow between populations. However, there is essentially no experimental data on whether this correlation arises because gene flow constrains adaptation or, conversely, because adaptive divergence causes barriers to gene flow (ecological speciation). Here, I report increased adaptive divergence in cryptic color pattern between a pair of Timema insect populations following an experimental reduction in between-population gene flow. The reduction in gene flow arose due to a natural experiment, and thus was not replicated at a second site. However, temporal replication of the trends among six generations of data, coupled with a lack of increased adaptive divergence for two other population pairs where gene flow was not manipulated (i.e., control sites), argues that the results did not arise by chance. Estimates of dispersal ability and population size further support reduced gene flow, rather than increased genetic drift, as the cause of divergence. Thus, the findings provide experimental evidence that gene flow constrains adaptation in nature.     

Local adaptation within a hybrid species

Ecological divergence among populations may be strongly influenced by their genetic background. For instance, genetic admixture through introgressive hybridization or hybrid speciation is likely to affect the genetic variation and evolvability of phenotypic traits. We studied geographic variation in two beak dimensions and three other phenotypic traits of the Italian sparrow (Passer italiae), a young hybrid species formed through interbreeding between house sparrows (P. domesticus) and Spanish sparrows (P. hispaniolensis). We found that beak morphology was strongly influenced by precipitation regimes and that it appeared to be the target of divergent selection within Italian sparrows. Interestingly, however, the degree of parental genetic contribution in the hybrid species had no effect on phenotypic beak variation. Moreover, beak height divergence may mediate genetic differentiation between populations, consistent with isolation-by-adaptation within this hybrid species. The study illustrates how hybrid species may be relatively unconstrained by their admixed genetic background, allowing them to adapt rapidly to environmental variation.     

Toxic hydrogen sulfide and dark caves: phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicana

Divergent natural selection drives evolutionary diversification. It creates phenotypic diversity by favoring developmental plasticity within populations or genetic differentiation and local adaptation among populations. We investigated phenotypic and genetic divergence in the livebearing fish Poecilia mexicana along two abiotic environmental gradients. These fish typically inhabit nonsulfidic surface rivers, but also colonized sulfidic and cave habitats. We assessed phenotypic variation among a factorial combination of habitat types using geometric and traditional morphometrics, and genetic divergence using quantitative and molecular genetic analyses. Fish in caves (sulfidic or not) exhibited reduced eyes and slender bodies. Fish from sulfidic habitats (surface or cave) exhibited larger heads and longer gill filaments. Common-garden rearing suggested that these morphological differences are partly heritable. Population genetic analyses using microsatellites as well as cytochrome b gene sequences indicate high population differentiation over small spatial scale and very low rates of gene flow, especially among different habitat types. This suggests that divergent environmental conditions constitute barriers to gene flow. Strong molecular divergence over short distances as well as phenotypic and quantitative genetic divergence across habitats in directions classic to fish ecomorphology suggest that divergent selection is structuring phenotypic variation in this system.     

A heuristic model on the role of plasticity in adaptive evolution: plasticity increases adaptation,population viability and genetic variation

An ongoing new synthesis in evolutionary theory is expanding our view of the sources of heritable variation beyond point mutations of fixed phenotypic effects to include environmentally sensitive changes in gene regulation. This expansion of the paradigm is necessary given ample evidence for a heritable ability to alter gene expression in response to environmental cues. In consequence, single genotypes are often capable of adaptively expressing different phenotypes in different environments, i.e. are adaptively plastic. We present an individual-based heuristic model to compare the adaptive dynamics of populations composed of plastic or non-plastic genotypes under a wide range of scenarios where we modify environmental variation, mutation rate and costs of plasticity. The model shows that adaptive plasticity contributes to the maintenance of genetic variation within populations, reduces bottlenecks when facing rapid environmental changes and confers an overall faster rate of adaptation. In fluctuating environments, plasticity is favoured by selection and maintained in the population. However, if the environment stabilizes and costs of plasticity are high, plasticity is reduced by selection, leading to genetic assimilation, which could result in species diversification. More broadly, our model shows that adaptive plasticity is a common consequence of selection under environmental heterogeneity, and hence a potentially common phenomenon in nature. Thus, taking adaptive plasticity into account substantially extends our view of adaptive evolution.     

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.     

ADAPTIVE RADIATION ALONG GENETIC LINES OF LEAST RESISTANCE

Are measurements of quantitative genetic variation useful for predicting long-term adaptive evolution? To answer this question, I focus on g_max, the multivariate direction of greatest additive genetic variance within populations. Original data on threespine sticklebacks, together with published genetic measurements from other vertebrates, show that morphological differentiation between species has been biased in the direction of g_max for at least four million years, despite evidence that natural selection is the cause of differentiation. This bias toward the direction of evolution tends to decay with time. Rate of morphological divergence between species is inversely proportional to θ, the angle between the direction of divergence and the direction of greatest genetic variation. The direction of greatest phenotypic variance is not identical with g_max, but for these data is nearly as successful at predicting the direction of species divergence. I interpret the findings to mean that genetic variances and covariances constrain adaptive change in quantitative traits for reasonably long spans of time. An alternative hypothesis, however, cannot be ruled out: that morphological differentiation is biased in the direction g_max because divergence and g_max are both shaped by the same natural selection pressures. Either way, the results reveal that adaptive differentiation occurs principally along “genetic lines of least resistance.”     

Bergmann's rule is maintained during a rapid range expansion in a damselfly

Climate‐induced range shifts result in the movement of a sample of genotypes from source populations to new regions. The phenotypic consequences of those shifts depend upon the sample characteristics of the dispersive genotypes, which may act to either constrain or promote phenotypic divergence, and the degree to which plasticity influences the genotype–environment interaction. We sampled populations of the damselfly Erythromma viridulum from northern Europe to quantify the phenotypic (latitude–body size relationship based on seven morphological traits) and genetic (variation at microsatellite loci) patterns that occur during a range expansion itself. We find a weak spatial genetic structure that is indicative of high gene flow during a rapid range expansion. Despite the potentially homogenizing effect of high gene flow, however, there is extensive phenotypic variation among samples along the invasion route that manifests as a strong, positive correlation between latitude and body size consistent with Bergmann's rule. This positive correlation cannot be explained by variation in the length of larval development (voltinism). While the adaptive significance of latitudinal variation in body size remains obscure, geographical patterns in body size in odonates are apparently underpinned by phenotypic plasticity and this permits a response to one or more environmental correlates of latitude during a range expansion.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Comparative landscape genetics and the adaptive radiation of Darwin's finches: the role of peripheral isolation

We use genetic divergence at 16 microsatellite loci to investigate how geographical features of the Galápagos landscape structure island populations of Darwin's finches. We compare the three most genetically divergent groups of Darwin's finches comprising morphologically and ecologically similar allopatric populations: the cactus finches (Geospiza scandens and Geospiza conirostris), the sharp-beaked ground finches (Geospiza difficilis) and the warbler finches (Certhidea olivacea and Certhidea fusca). Evidence of reduced genetic diversity due to drift was limited to warbler finches on small, peripheral islands. Evidence of low levels of recent interisland migration was widespread throughout all three groups. The hypothesis of distance-limited dispersal received the strongest support in cactus and sharp-beaked ground finches as evidenced by patterns of isolation by distance, while warbler finches showed a weaker relationship. Support for the hypothesis that gene flow constrains morphological divergence was only found in one of eight comparisons within these groups. Among warbler finches, genetic divergence was relatively high while phenotypic divergence was low, implicating stabilizing selection rather than constraint due to gene flow. We conclude that the adaptive radiation of Darwin's finches has occurred in the presence of ongoing but low levels of gene flow caused by distance-dependent interisland dispersal. Gene flow does not constrain phenotypic divergence, but may augment genetic variation and facilitate evolution due to natural selection. Both microsatellites and mtDNA agree in that subsets of peripheral populations of two older groups are genetically more similar to other species that underwent dramatic morphological change. The apparent decoupling of morphological and molecular evolution may be accounted for by a modification of Lack's two-stage model of speciation: relative ecological stasis in allopatry followed by secondary contact, ecological interactions and asymmetric phenotypic divergence.     

The consequences of phenotypic plasticity for ecological speciation

We use an individual-based numerical simulation to study the effects of phenotypic plasticity on ecological speciation. We find that adaptive plasticity evolves readily in the presence of dispersal between populations from different ecological environments. This plasticity promotes the colonization of new environments but reduces genetic divergence between them. We also find that the evolution of plasticity can either enhance or degrade the potential for divergent selection to form reproductive barriers. Of particular importance here is the timing of plasticity in relation to the timing of dispersal. If plasticity is expressed after dispersal, reproductive barriers are generally weaker because plasticity allows migrants to be better suited for their new environment. If plasticity is expressed before dispersal, reproductive barriers are either unaffected or enhanced. Among the potential reproductive barriers we considered, natural selection against migrants was the most important, primarily because it was the earliest-acting barrier. Accordingly, plasticity had a much greater effect on natural selection against migrants than on sexual selection against migrants or on natural and sexual selection against hybrids. In general, phenotypic plasticity can strongly alter the process of ecological speciation and should be considered when studying the evolution of reproductive barriers.     

Potential local adaptation in populations of invasive reed canary grass (Phalaris arundinacea) across an urbanization gradient

Urban stressors represent strong selective gradients that can elicit evolutionary change, especially in non‐native species that may harbor substantial within‐population variability. To test whether urban stressors drive phenotypic differentiation and influence local adaptation, we compared stress responses of populations of a ubiquitous invader, reed canary grass (Phalaris arundinacea). Specifically, we quantified responses to salt, copper, and zinc additions by reed canary grass collected from four populations spanning an urbanization gradient (natural, rural, moderate urban, and intense urban). We measured ten phenotypic traits and trait plasticities, because reed canary grass is known to be highly plastic and because plasticity may enhance invasion success. We tested the following hypotheses: (a) Source populations vary systematically in their stress response, with the intense urban population least sensitive and the natural population most sensitive, and (b) plastic responses are adaptive under stressful conditions. We found clear trait variation among populations, with the greatest divergence in traits and trait plasticities between the natural and intense urban populations. The intense urban population showed stress tolerator characteristics for resource acquisition traits including leaf dry matter content and specific root length. Trait plasticity varied among populations for over half the traits measured, highlighting that plasticity differences were as common as trait differences. Plasticity in root mass ratio and specific root length were adaptive in some contexts, suggesting that natural selection by anthropogenic stressors may have contributed to root trait differences. Reed canary grass populations in highly urbanized wetlands may therefore be evolving enhanced tolerance to urban stressors, suggesting a mechanism by which invasive species may proliferate across urban wetland systems generally.     

Integrating patterns of thermal tolerance and phenotypic plasticity with population genetics to improve understanding of vulnerability to warming in a widespread copepod

Differences in population vulnerability to warming are defined by spatial patterns in thermal adaptation. These patterns may be driven by natural selection over spatial environmental gradients, but can also be shaped by gene flow, especially in marine taxa with high dispersal potential. Understanding and predicting organismal responses to warming requires disentangling the opposing effects of selection and gene flow. We begin by documenting genetic divergence of thermal tolerance and developmental phenotypic plasticity. Ten populations of the widespread copepod Acartia tonsa were collected from sites across a large thermal gradient, ranging from the Florida Keys to Northern New Brunswick, Canada (spanning over 20° latitude). Thermal performance curves (TPCs) from common garden experiments revealed local adaptation at the sampling range extremes, with thermal tolerance increasing at low latitudes and decreasing at high latitudes. The opposite pattern was observed in phenotypic plasticity, which was strongest at high latitudes. No relationship was observed between phenotypic plasticity and environmental variables. Instead, the results are consistent with the hypothesis of a trade‐off between thermal tolerance and the strength of phenotypic plasticity. Over a large portion of the sampled range, however, we observed a remarkable lack of differentiation of TPCs. To examine whether this lack of divergence is the result of selection for a generalist performance curve or constraint by gene flow, we analyzed cytochrome oxidase I mtDNA sequences, which revealed four distinct genetic clades, abundant genetic diversity, and widely distributed haplotypes. Strong divergence in thermal performance within genetic clades, however, suggests that the pace of thermal adaptation can be relatively rapid. The combined insight from the laboratory physiological experiments and genetic data indicate that gene flow constrains differentiation of TPCs. This balance between gene flow and selection has implications for patterns of vulnerability to warming. Taking both genetic differentiation and phenotypic plasticity into account, our results suggest that local adaptation does not increase vulnerability to warming, and that low‐latitude populations in general may be more vulnerable to predicted temperature change over the next century.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

Natural selection contributes to geographic patterns of thermal plasticity in Plantago lanceolata

A long‐standing debate in evolutionary biology concerns the relative importance of different evolutionary forces in explaining phenotypic diversification at large geographic scales. For example, natural selection is typically assumed to underlie divergence along environmental gradients. However, neutral evolutionary processes can produce similar patterns. We collected molecular genetic data from 14 European populations of Plantago lanceolata to test the contributions of natural selection versus neutral evolution to population divergence in temperature‐sensitive phenotypic plasticity of floral reflectance. In P. lanceolata, reflectance plasticity is positively correlated with latitude/altitude. We used population pairwise comparisons between neutral genetic differentiation (F_ST and Jost's D) and phenotypic differentiation (P_ST) to assess the contributions of geographic distance and environmental parameters of the reproductive season in driving population divergence. Data are consistent with selection having shaped large‐scale geographic patterns in thermal plasticity. The aggregate pattern of P_ST versus F_ST was consistent with divergent selection. F_ST explained thermal plasticity differences only when geographic distance was not included in the model. Differences in the extent of cool reproductive season temperatures, and not overall temperature variation, explained plasticity differences independent of distance. Results are consistent with the hypothesis that thermal plasticity is adaptive where growing seasons are shorter and cooler, that is, at high latitude/altitude.     

Patterns of phenotypic and genetic variation for the plasticity of diapause incidence

Phenotypic plasticity describes an organism's ability to produce multiple phenotypes in direct response to its environmental conditions. Over the past 15 years empiricists have found that this plasticity frequently exhibits geographic variation and often possesses a significant heritable genetic basis. However, few studies have examined both of these aspects of plasticity simultaneously. Here, we examined both the geographic and genetic variations of the plasticity for diapause incidence (the proportion of eggs that enter an arrested state of development capable of surviving over the winter) relative to temperatures and photoperiods associated with long and short season environments across six populations of the striped ground cricket, Allonemobius socius, using a half-sibling split brood quantitative genetic design. We found that plasticity, as measured by the slope of the reaction norm, was greater in the southern-low altitude region (where populations are bivoltine) relative to the southern-high and northern-low altitude regions (where populations are univoltine). However, the heritability of plasticity was only significantly different from zero in univoltine populations that experienced "intermediate" natal season lengths. These patterns suggest that selection may favor the plasticity of diapause incidence in bivoltine regions, but act against plasticity in regions in which populations are univoltine. Furthermore, our data suggest that under "intermediate" natal season length conditions, the interplay between local adaptation and gene flow may keep the plasticity of diapause incidence low (but still significant) while maintaining its genetic variation. As such, this study not only provides a novel observation into the geographic variation of phenotypic plasticity, but also provides much needed groundwork for tests of its adaptive significance.     

Geographic variation in phenotypic plasticity in response to dissolved oxygen in an African cichlid fish

Genetic adaptation and phenotypic plasticity are two ways in which organisms can adapt to local environmental conditions. We examined genetic and plastic variation in gill and brain size among swamp (low oxygen; hypoxic) and river (normal oxygen; normoxic) populations of an African cichlid fish, Pseudocrenilabrus multicolor victoriae. Larger gills and smaller brains should be advantageous when oxygen is low, and we hypothesized that the relative contribution of local genetic adaptation vs. phenotypic plasticity should be related to potential for dispersal between environments (because of gene flow’s constraint on local genetic adaptation). We conducted a laboratory‐rearing experiment, with broods from multiple populations raised under high‐oxygen and low‐oxygen conditions. We found that most of the variation in gill size was because of plasticity. However, both plastic and genetic effects on brain mass were detected, as were genetic effects on brain mass plasticity. F₁ offspring from populations with the highest potential for dispersal between environments had characteristically smaller and more plastic brains. This phenotypic pattern might be adaptive in the face of gene flow, if smaller brains and increased plasticity confer higher average fitness across environment types.