Local selection modifies phenotypic divergence amongRana temporariapopulations in the presence of gene flow

Food and light availability induce plastic responses in fire salamander larvae from contrasting environments

Phenotypic plasticity has been proposed as a mechanism facilitating the colonisation and adaptation to novel environments, such as caves. However, phenotypic plasticity in subterranean environments remains largely unexplored. Here, we test for plasticity in growth and development of fire salamander larvae (Salamandra salamandra) from subterranean and surface habitats, in response to contrasting food availability and light conditions. We hypothesized that: (i) low food availability and absence of light decrease larval growth and delay metamorphosis, (ii) light conditions mediate the effects of food availability on growth and time to metamorphosis, and (iii) larval response to contrasting light and food conditions is shaped by the habitat of origin. Our study showed that reduced food availability significantly delayed metamorphosis and slowed total length and body mass growth rates, while exposure to constant darkness slowed body mass growth rate. However, larvae slowed growth rates and increased time to metamorphosis without compromising size at metamorphosis. The effect of food availability on growth and time to metamorphosis did not change under different light conditions. Fire salamanders from subterranean and surface habitats responded differently only in relation to contrasting food availability conditions. Specifically, larvae from the surface habitat grew faster in high food conditions, while growth in larvae from the subterranean habitat was not influenced by food availability. Initial size also appeared to be an influential factor, since larger and heavier larvae grew slower, metamorphosed faster, and the size advantage was maintained in newly-metamorphosed juveniles. Overall, the results of our experiment suggest that plasticity and local adaptation favor the exploitation of aquatic subterranean habitats for breeding by fire salamanders, allowing successful development even under food shortage and day-length constraints, without compromising metamorphic size. Our findings have implications for conservation because they confirm that phenotypic plasticity plays a critical role in allowing fire salamanders to overcome altered environmental conditions.

Phenotypic variation in Xenopus laevis tadpoles from contrasting climatic regimes is the result of adaptation and plasticity

Phenotypic variations between populations often correlate with climatic variables. Determining the presence of phenotypic plasticity and local adaptation of a species to different environments over a large spatial scale can provide insight on the persistence of a species across its range. Amphibians, and in particular their larvae, are good models for studies of phenotypic variation as they are especially sensitive to their immediate environment. Few studies have attempted to determine the mechanisms that drive phenotypic variation between populations of a single amphibian species over a large spatial scale especially across contrasting climatic regimes. The African clawed frog, Xenopus laevis, occurs in two regions with contrasting rainfall regimes in southern Africa. We hypothesised that the phenotypic variation of life-history traits of X. laevis tadpoles emerges from a combination of plastic and genetic responses. We predicted that plasticity would allow the development of tadpoles from both regions in each environment. We also predicted that local adaptation of larval traits would drive the differentiation of reaction norms between populations and lower survival in tadpoles reared away from their home environment. We measured growth, time to metamorphosis, and survival in a reciprocal transplant experiment using outdoor mesocosms. Supporting our prediction, we found that the measured variation of all traits was explained by both adaptation and plasticity. However, the reaction norms differed between populations suggesting adaptive and asymmetric plasticity. All tadpoles experienced lower survival when translocated, but only translocated tadpoles from the winter rainfall region matched survival of local tadpoles. This has implications for the dynamics of translocated X. laevis into novel environments, especially from the winter rainfall region. Our discovery of their asymmetric capacity to overcome novel environmental conditions by phenotypic plasticity alone provides insight into their invasion success.

Thermal tolerance and preference are both consistent with the clinal distribution of house fly proto-Y chromosomes

Selection pressures can vary within localized areas and across massive geographical scales. Temperature is one of the best studied ecologically variable abiotic factors that can affect selection pressures across multiple spatial scales. Organisms rely on physiological (thermal tolerance) and behavioral (thermal preference) mechanisms to thermoregulate in response to environmental temperature. In addition, spatial heterogeneity in temperatures can select for local adaptation in thermal tolerance, thermal preference, or both. However, the concordance between thermal tolerance and preference across genotypes and sexes within species and across populations is greatly understudied. The house fly, Musca domestica, is a well-suited system to examine how genotype and environment interact to affect thermal tolerance and preference. Across multiple continents, house fly males from higher latitudes tend to carry the male-determining gene on the Y chromosome, whereas those from lower latitudes usually have the male determiner on the third chromosome. We tested whether these two male-determining chromosomes differentially affect thermal tolerance and preference as predicted by their geographical distributions. We identify effects of genotype and developmental temperature on male thermal tolerance and preference that are concordant with the natural distributions of the chromosomes, suggesting that temperature variation across the species range contributes to the maintenance of the polymorphism. In contrast, female thermal preference is bimodal and largely independent of congener male genotypes. These sexually dimorphic thermal preferences suggest that temperature-dependent mating dynamics within populations could further affect the distribution of the two chromosomes. Together, the differences in thermal tolerance and preference across sexes and male genotypes suggest that different selection pressures may affect the frequencies of the male-determining chromosomes across different spatial scales.

Ecological adaptation drives wood frog population divergence in life history traits

Phenotypic variation among populations is thought to be generated from spatial heterogeneity in environments that exert selection pressures that overcome the effects of gene flow and genetic drift. Here, we tested for evidence of isolation by distance or by ecology (i.e., ecological adaptation) to generate variation in early life history traits and phenotypic plasticity among 13 wood frog populations spanning 1200 km and 7° latitude. We conducted a common garden experiment and related trait variation to an ecological gradient derived from an ecological niche model (ENM) validated to account for population density variation. Shorter larval periods, smaller body weight, and relative leg lengths were exhibited by populations with colder mean annual temperatures, greater precipitation, and less seasonality in precipitation and higher population density (high-suitability ENM values). After accounting for neutral genetic variation, the Q_ST-F_ST analysis supported ecological selection as the key process generating population divergence. Further, the relationship between ecology and traits was dependent upon larval density. Specifically, high-suitability/high-density populations in the northern part of the range were better at coping with greater conspecific competition, evidenced by greater postmetamorphic survival and no difference in body weight when reared under stressful conditions of high larval density. Our results support that both climate and competition selection pressures drive clinal variation in larval and metamorphic traits in this species. Range-wide studies like this one are essential for accurate predictions of population's responses to ongoing ecological change.

Parental effects influence life history traits and covary with an environmental cline in common frog populations

Across latitudinal clines, the juvenile developmental rates of ectotherms often covary with the length of the growing season, due to life-history trade-offs imposed by the time-constrained environments. However, as the start of the growing season often varies substantially across years, adaptive parental effects on juvenile developmental rates may mediate the costs of a delayed season. By employing a meta-analysis, we tested whether larval developmental rates across a latitudinal cline of the common frog (Rana temporaria) are affected by fluctuating onsets of breeding, across years. We predicted that larval developmental rate will be inversely related to the onset of breeding, and that northern populations will be more prone to shorten their developmental rate in response to late breeding, as the costs of delayed metamorphosis should be highest in areas with a shorter growing season. We found that the larval period of both northern and southern populations responded to parental environmental conditions to a similar degree in absolute terms, but in different directions. In northern populations, a late season start correlated with decreased development time, suggesting that the evolution of parental effects aids population persistence in time-constrained environments. In southern populations, late season start correlated with increased development time, which could potentially be explained as a predator avoidance strategy. Our findings suggest that local ecological variables can induce adaptive parental effects, but responses are complex, and likely trade-off with other ecological factors.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Discovering divergence in the thermal physiology of intertidal crabs along latitudinal gradients using an integrated approach with machine learning

In intertidal marine crustaceans, phenotypic variation in physiological and life-history traits is pervasive along latitudinal clines. However, organisms have complex phenotypes, and their traits do not vary independently but rather interact differentially between them, effect that is caused by genetic and/or environmental forces. We evaluated the geographic variation in phenotypic integration of three marine crab species that inhabit different vertical thermal microhabitats of the intertidal zone. We studied seven populations of each species along a latitudinal gradient that spans more than 3000 km of the Chilean coast. Specifically we measured nine physiological traits that are highly related to thermal physiology. Of the nine traits, we selected four that contributed significantly to the observed geographical variation among populations; this variation was then evaluated using mixed linear models and an integrative approach employing machine learning. The results indicate that patterns of physiological variation depend on species vertical microhabitat, which may be subject to chronic or acute environmental variation. The species that inhabit the high- intertidal sites (i.e., exposed to chronic variation) better tolerated thermal stress compared with populations that inhabit the lower intertidal. While those in the low-intertidal only face conditions of acute thermal variation, using to a greater extent the plasticity to face these events. Our main results reflect that (1) species that inhabit the high-intertidal maintain a greater integration between their physiological traits and present lower plasticity than those that inhabit the low-intertidal. (2) Inverse relationship that exists between phenotypic plasticity and phenotypic integration of the physiological traits identified, which could help optimize energy resources. In general, the study of multiple physiological traits provides a more accurate picture of how the thermal traits of organisms vary along temperature gradients especially when exposed to conditions close to tolerance limits.

Deciphering the drivers of negative species-genetic diversity correlation in Alpine amphibians

The evolutionary and ecological importance of neutral and adaptive genetic diversity is widely recognized. Nevertheless, genetic diversity is rarely assessed for conservation planning, which often implicitly assumes a positive correlation between species and genetic diversity. Multiple drivers can cause the co-variation between the genetic diversity of one species and the richness of the whole communities, and explicit tests are needed to identify the processes that can determine species-genetic diversity correlations (SGDCs). Here, we tested whether intrapopulation genetic diversity (at neutral loci) and species richness co-vary in the amphibian communities of a southern Alpine region (Trentino, Italy), using the common frog (Rana temporaria) as focal species for the study of genetic diversity. We also analysed ecological similarity, niche overlap and interspecific interactions between the species, to unravel the processes determining SGDC. The neutral genetic diversity of common frogs was negatively related to species richness. The negative SGDC was probably due to an opposite influence of environmental gradients on the two levels of biodiversity, since the focal species and the other amphibians differ in ecological preferences, particularly in terms of thermal optimum. Conversely, we did not find evidence for a role of interspecific interactions in the negative SGDC. Our findings stress that species richness cannot be used as a universal proxy for genetic diversity, and only combining SGDC with analyses on the determinants of biodiversity can allow to identify the processes determining the relationships between genetic and species diversity.

Latitudinal variation in biophysical characteristics of avian eggshells to cope with differential effects of solar radiation

Solar radiation is an important driver of animal coloration, not only because of the effects of coloration on body temperature but also because coloration may protect from the deleterious effects of UV radiation. Indeed, dark coloration may protect from UV, but may increase the risk of overheating. In addition, the effect of coloration on thermoregulation should change with egg size, as smaller eggs have higher surface-volume ratios and greater convective coefficients than larger eggs, so that small eggs can dissipate heat quickly. We tested whether the reflectance of eggshells, egg spottiness, and egg size of the ground-nesting Kentish plover Charadrius alexandrinus is affected by maximum ambient temperature and solar radiation at breeding sites. We measured reflectance, both in the UV and human visible spectrum, spottiness, and egg size in photographs from a museum collection of plover eggshells. Eggshells of lower reflectance (darker) were found at higher latitudes. However, in southern localities where solar radiation is very high, eggshells are also of dark coloration. Eggshell coloration had no significant relationship with ambient temperature. Spotiness was site-specific. Small eggs tended to be light-colored. Thermal constraints may drive the observed spatial variation in eggshell coloration, which may be lighter in lower latitudes to diminish the risk of overheating as a result of higher levels of solar radiation. However, in southern localities with very high levels of UV radiation, eggshells are of dark coloration likely to protect embryos from more intense UV radiation. Egg size exhibited variation in relation to coloration, likely through the effect of surface area-to-volume ratios on overheating and cooling rates of eggs. Therefore, differential effects of solar radiation on functions of coloration and size of eggshells may shape latitudinal variations in egg appearance in the Kentish plover.

Temperature-driven selection on metabolic traits increases the strength of an algal-grazer interaction in naturally warmed streams

Trophic interactions are important determinants of the structure and functioning of ecosystems. Because the metabolism and consumption rates of ectotherms increase sharply with temperature, there are major concerns that global warming will increase the strength of trophic interactions, destabilizing food webs, and altering ecosystem structure and function. We used geothermally warmed streams that span an 11°C temperature gradient to investigate the interplay between temperature-driven selection on traits related to metabolism and resource acquisition, and the interaction strength between the keystone gastropod grazer, Radix balthica, and a common algal resource. Populations from a warm stream (~28°C) had higher maximal metabolic rates and optimal temperatures than their counterparts from a cold stream (~17°C). We found that metabolic rates of the population originating from the warmer stream were higher across all measurement temperatures. A reciprocal transplant experiment demonstrated that the interaction strengths between the grazer and its algal resource were highest for both populations when transplanted into the warm stream. In line with the thermal dependence of respiration, interaction strengths involving grazers from the warm stream were always higher than those with grazers from the cold stream. These results imply that increases in metabolism and resource consumption mediated by the direct, thermodynamic effects of higher temperatures on physiological rates are not mitigated by metabolic compensation in the long term, and suggest that warming could increase the strength of algal-grazer interactions with likely knock-on effects for the biodiversity and productivity of aquatic ecosystems.

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open- and closed-canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade-off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open- and closed-canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.

Microgeographic differentiation in thermal performance curves between rural and urban populations of an aquatic insect

The rapidly increasing rate of urbanization has a major impact on the ecology and evolution of species. While increased temperatures are a key aspect of urbanization ("urban heat islands"), we have very limited knowledge whether this generates differentiation in thermal responses between rural and urban populations. In a common garden experiment, we compared the thermal performance curves (TPCs) for growth rate and mortality in larvae of the damselfly Coenagrion puella from three urban and three rural populations. TPCs for growth rate shifted vertically, consistent with the faster-slower theoretical model whereby the cold-adapted rural larvae grew faster than the warm-adapted urban larvae across temperatures. In line with costs of rapid growth, rural larvae showed lower survival than urban larvae across temperatures. The relatively lower temperatures hence expected shorter growing seasons in rural populations compared to the populations in the urban heat islands likely impose stronger time constraints to reach a certain developmental stage before winter, thereby selecting for faster growth rates. In addition, higher predation rates at higher temperature may have contributed to the growth rate differences between urban and rural ponds. A faster-slower differentiation in TPCs may be a widespread pattern along the urbanization gradient. The observed microgeographic differentiation in TPCs supports the view that urbanization may drive life-history evolution. Moreover, because of the urban heat island effect, urban environments have the potential to aid in developing predictions on the impact of climate change on rural populations.

Investigating population differentiation in a major African agricultural pest: evidence from geometric morphometrics and connectivity suggests high invasion potential

The distribution, spatial pattern and population dynamics of a species can be influenced by differences in the environment across its range. Spatial variation in climatic conditions can cause local populations to undergo disruptive selection and ultimately result in local adaptation. However, local adaptation can be constrained by gene flow and may favour resident individuals over migrants-both are factors critical to the assessment of invasion potential. The Natal fruit fly (Ceratitis rosa) is a major agricultural pest in Africa with a history of island invasions, although its range is largely restricted to south east Africa. Across Africa, C. rosa is genetically structured into two clusters (R1 and R2), with these clusters occurring sympatrically in the north of South Africa. The spatial distribution of these genotypic clusters remains unexamined despite their importance for understanding the pest's invasion potential. Here, C. rosa, sampled from 22 South African locations, were genotyped at 11 polymorphic microsatellite loci and assessed morphologically using geometric morphometric wing shape analyses to investigate patterns of population structure and determine connectedness of pest-occupied sites. Our results show little to no intraspecific (population) differentiation, high population connectivity, high effective population sizes and only one morphological type (R2) within South Africa. The absence of the R1 morphotype at sites where it was previously found may be a consequence of differences in thermal niches of the two morphotypes. Overall, our results suggest high invasion potential of this species, that area-wide pest management should be undertaken on a country-wide scale, and that border control is critical to preventing further invasions.

Population-specific effects of developmental temperature on body condition and jumping performance of a widespread European frog

All physiological processes of ectotherms depend on environmental temperature. Thus, adaptation of physiological mechanisms to the thermal environments is important for achieving optimal performance and fitness. The European Common Frog, Rana temporaria, is widely distributed across different thermal habitats. This makes it an exceptional model for studying the adaptations to different thermal conditions. We raised tadpoles from Germany and Croatia at two constant temperature treatments (15°C, 20°C), and under natural temperature fluctuations (in outdoor treatments), and tested how different developmental temperatures affected developmental traits, that is, length of larval development, morphometrics, and body condition, as well as jumping performance of metamorphs. Our results revealed population‐specific differences in developmental time, body condition, and jumping performance. Croatian frogs developed faster in all treatments, were heavier, in better body condition, and had longer hind limbs and better jumping abilities than German metamorphs. The populations further differed in thermal sensitivity of jumping performance. While metamorphs from Croatia increased their jumping performance with higher temperatures, German metamorphs reached their performance maximum at lower temperatures. These population‐specific differences in common environments indicate local genetic adaptation, with southern populations being better adapted to higher temperatures than those from north of the Alps.

Genetic variation in plasticity of life-history traits between Atlantic cod (Gadus morhua) populations exposed to contrasting thermal regimes

We employed common-garden experiments to test for genetic variation in responses of larval life-history traits to temperature between two populations of Atlantic cod (Gadus morhua L., 1758) that naturally experience contrasting thermal environments during early life due to spatial and temporal differences in spawning. Southern Gulf of St. Lawrence cod larvae experienced faster growth in warmer water and low, uniform survival across all experimental temperatures (3, 7, 11 °C), consistent with previous studies on this spring-spawning population. In contrast, larvae from fall-spawning Southwestern Scotian Shelf cod collected near Sambro, Nova Scotia, lacked plasticity for growth but experienced much lower survival at higher temperatures. Phenotypes that are positively associated with fitness were observed at temperatures closest to those experienced in the wild, consistent with the hypothesis that these populations are adapted to local thermal regimes. The lack of growth plasticity observed in Sambro cod might be due to costly maintenance of plasticity in stable environments or energy savings at cold temperatures. However, additional experiments need to be conducted on Sambro cod and other fall-spawning marine fishes to determine to what extent responses to projected changes in climate will differ among populations.

Expression Differentiation Is Constrained to Low-Expression Proteins over Ecological Timescales

Protein expression level is one of the strongest predictors of protein sequence evolutionary rate, with high-expression protein sequences evolving at slower rates than low-expression protein sequences largely because of constraints on protein folding and function. Expression evolutionary rates also have been shown to be negatively correlated with expression level across human and mouse orthologs over relatively long divergence times (i.e., ∼100 million years). Long-term evolutionary patterns, however, often cannot be extrapolated to microevolutionary processes (and vice versa), and whether this relationship holds for traits evolving under directional selection within a single species over ecological timescales (i.e., <5000 years) is unknown and not necessarily expected. Expression is a metabolically costly process, and the expression level of a particular protein is predicted to be a tradeoff between the benefit of its function and the costs of its expression. Selection should drive the expression level of all proteins close to values that maximize fitness, particularly for high-expression proteins because of the increased energetic cost of production. Therefore, stabilizing selection may reduce the amount of standing expression variation for high-expression proteins, and in combination with physiological constraints that may place an upper bound on the range of beneficial expression variation, these constraints could severely limit the availability of beneficial expression variants. To determine whether rapid-expression evolution was restricted to low-expression proteins owing to these constraints on highly expressed proteins over ecological timescales, we compared venom protein expression levels across mainland and island populations for three species of pit vipers. We detected significant differentiation in protein expression levels in two of the three species and found that rapid-expression differentiation was restricted to low-expression proteins. Our results suggest that various constraints on high-expression proteins reduce the availability of beneficial expression variants relative to low-expression proteins, enabling low-expression proteins to evolve and potentially lead to more rapid adaptation.

Local divergence of thermal reaction norms among amphibian populations is affected by pond temperature variation

Although temperature variation is known to cause large-scale adaptive divergence, its potential role as a selective factor over microgeographic scales is less well-understood. Here, we investigated how variation in breeding pond temperature affects divergence in multiple physiological (thermal performance curve and critical thermal maximum [CTmax]) and life-history (thermal developmental reaction norms) traits in a network of Rana arvalis populations. The results supported adaptive responses to face two main constraints limiting the evolution of thermal adaptation. First, we found support for the faster-slower model, indicating an adaptive response to compensate for the thermodynamic constraint of low temperatures in colder environments. Second, we found evidence for the generalist-specialist trade-off with populations from colder and less thermally variable environments exhibiting a specialist phenotype performing at higher rates but over a narrower range of temperatures. By contrast, the local optimal temperature for locomotor performance and CTmax did not match either mean or maximum pond temperatures. These results highlight the complexity of the adaptive multiple-trait thermal responses in natural populations, and the role of local thermal variation as a selective force driving diversity in life-history and physiological traits in the presence of gene flow.

Variation in spawning time promotes genetic variability in population responses to environmental change in a marine fish

The level of phenotypic plasticity displayed within a population (i.e. the slope of the reaction norm) reflects the short-term response of a population to environmental change, while variation in reaction norm slopes among populations reflects spatial variation in these responses. Thus far, studies of thermal reaction norm variation have focused on geographically driven adaptation among different latitudes, altitudes or habitats. Yet, thermal variability is a function of both space and time. For organisms that reproduce at different times of year, such variation has the potential to promote adaptive variability in thermal responses for critical early life stages. Using common-garden experiments, we examined the spatial scale of genetic variation in thermal plasticity for early life-history traits among five populations of endangered Atlantic cod (Gadus morhua) that spawn at different times of year. Patterns of plasticity for larval growth and survival suggest that population responses to climate change will differ substantially, with increasing water temperatures posing a considerably greater threat to autumn-spawning cod than to those that spawn in winter or spring. Adaptation to seasonal cooling or warming experienced during the larval stage is suggested as a possible cause. Furthermore, populations that experience relatively cold temperatures during early life might be more sensitive to changes in temperature. Substantial divergence in adaptive traits was evident at a smaller spatial scale than has previously been shown for a marine fish with no apparent physical barriers to gene flow (∼200 km). Our findings highlight the need to consider the impact of intraspecific variation in reproductive timing on thermal adaptation when forecasting the effects of climate change on animal populations.

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age-specific selection is constrained to act on the same genome by altering age-specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype-phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.

Rapid evolution of constitutive and inducible defenses against an invasive predator

Invasive alien predators can impose strong selection on native prey populations and induce rapid evolutionary change in the invaded communities. However, studies on evolutionary responses to invasive predators are often complicated by the lack of replicate populations differing in coexistence time with the predator, which would allow the determination of how prey traits change during the invasion. The red swamp crayfish Procambarus clarkii has invaded many freshwater areas worldwide, with negative impacts for native fauna. Here, we examined how coexistence time shapes antipredator responses of the Iberian waterfrog (Pelophylax perezi) to the invasive crayfish by raising tadpoles from five populations differing in historical exposure to P. clarkii (30 years, 20 years, or no coexistence). Tadpoles from non-invaded populations responded to the presence of P. clarkii with behavioral plasticity (reduced activity), whereas long-term invaded populations showed canalized antipredator behavior (constant low activity level). Tadpoles from one of the long-term invaded populations responded to the crayfish with inducible morphological defenses (deeper tails), reflecting the use of both constitutive and inducible antipredator defenses against the exotic predator by this population. Our results suggest that, while naive P. perezi populations responded behaviorally to P. clarkii, the strong predation pressure imposed by the crayfish has induced the evolution of qualitatively different antipredator defenses in populations with longer coexistence time. These responses suggest that strong selection by invasive predators may drive rapid evolutionary change in invaded communities. Examining responses of prey species to biological invasions using multiple populations will help us better forecast the impact of invasive predators in natural communities.

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

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

How populations differentiate despite gene flow: sexual and natural selection drive phenotypic divergence within a land fish, the Pacific leaping blenny

Background

Divergence between populations in reproductively important features is often vital for speciation. Many studies attempt to identify the cause of population differentiation in phenotype through the study of a specific selection pressure. Holistic studies that consider the interaction of several contrasting forms of selection are more rare. Most studies also fail to consider the history of connectivity among populations and the potential for genetic drift or gene flow to facilitate or limit phenotypic divergence. We examined the interacting effects of natural selection, sexual selection and the history of connectivity on phenotypic differentiation among five populations of the Pacific leaping blenny (Alticus arnoldorum), a land fish endemic to the island of Guam.

Results

We found key differences among populations in two male ornaments—the size of a prominent head crest and conspicuousness of a coloured dorsal fin—that reflected a trade-off between the intensity of sexual selection (male biased sex ratios) and natural selection (exposure to predators). This differentiation in ornamentation has occurred despite evidence suggesting extensive gene flow among populations, which implies that the change in ornament expression has been recent (and potentially plastic).

Conclusions

Our study provides an early snapshot of divergence in reproductively important features that, regardless of whether it reflects genetic or plastic changes in phenotype, could ultimately form a reproductive barrier among populations.

Sensory trait variation in an echolocating bat suggests roles for both selection and plasticity

Background

Across heterogeneous environments selection and gene flow interact to influence the rate and extent of adaptive trait evolution. This complex relationship is further influenced by the rarely considered role of phenotypic plasticity in the evolution of adaptive population variation. Plasticity can be adaptive if it promotes colonization and survival in novel environments and in doing so may increase the potential for future population differentiation via selection. Gene flow between selectively divergent environments may favour the evolution of phenotypic plasticity or conversely, plasticity itself may promote gene flow, leading to a pattern of trait differentiation in the presence of gene flow. Variation in sensory traits is particularly informative in testing the role of environment in trait and population differentiation. Here we test the hypothesis of ‘adaptive differentiation with minimal gene flow’ in resting echolocation frequencies (RF) of Cape horseshoe bats (Rhinolophus capensis) across a gradient of increasingly cluttered habitats.

Results

Our analysis reveals a geographically structured pattern of increasing RF from open to highly cluttered habitats in R. capensis; however genetic drift appears to be a minor player in the processes influencing this pattern. Although Bayesian analysis of population structure uncovered a number of spatially defined mitochondrial groups and coalescent methods revealed regional-scale gene flow, phylogenetic analysis of mitochondrial sequences did not correlate with RF differentiation. Instead, habitat discontinuities between biomes, and not genetic and geographic distances, best explained echolocation variation in this species. We argue that both selection for increased detection distance in relatively less cluttered habitats and adaptive phenotypic plasticity may have influenced the evolution of matched echolocation frequencies and habitats across different populations.

Conclusions

Our study reveals significant sensory trait differentiation in the presence of historical gene flow and suggests roles for both selection and plasticity in the evolution of echolocation variation in R. capensis. These results highlight the importance of population level analyses to i) illuminate the subtle interplay between selection, plasticity and gene flow in the evolution of adaptive traits and ii) demonstrate that evolutionary processes may act simultaneously and that their relative influence may vary across different environments.

Giants, dwarfs and the environment - metamorphic trait plasticity in the common frog

In order to understand adaptation processes and population dynamics, it is central to know how environmental parameters influence performance of organisms within populations, including their phenotypes. The impact of single or few particular parameters in concert was often assessed in laboratory and mesocosm experiments. However, under natural conditions, with many biotic and abiotic factors potentially interacting, outcomes on phenotypic changes may be different. To study the potential environmental impact on realized phenotypic plasticity within a natural population, we assessed metamorphic traits (developmental time, size and body mass) in an amphibian species, the European common frog Rana temporaria, since a) larval amphibians are known to exhibit high levels of phenotypic plasticity of these traits in response to habitat parameters and, b) the traits' features may strongly influence individuals' future performance and fitness. In 2007 we studied these metamorphic traits in 18 ponds spread over an area of 28 km2. A subset of six ponds was reinvestigated in 2009 and 2010. This study revealed locally high variances in metamorphic traits in this presumed generalist species. We detected profound differences between metamorphing froglets (up to factor ten); both between and within ponds, on a very small geographic scale. Parameters such as predation and competition as well as many other pond characteristics, generally expected to have high impact on development, could not be related to the trait differences. We observed high divergence of patterns of mass at metamorphosis between ponds, but no detectable pattern when metamorphic traits were compared between ponds and years. Our results indicate that environment alone, i.e. as experienced by tadpoles sharing the same breeding pond, can only partly explain the variability of metamorphic traits observed. This emphasizes the importance to assess variability of reaction norms on the individual level to explain within-population variability.

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high- and low-altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high- and low-altitude sites using three common temperature treatments. Local adaptation was assessed using Q(ST)-F(ST) analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high-altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.

Compensatory growth strategies are affected by the strength of environmental time constraints in anuran larvae

Organisms normally grow at a sub-maximal rate. After experiencing a period of arrested growth, individuals often show compensatory growth responses by modifying their life-history, behaviour and physiology. However, the strength of compensatory responses may vary across broad geographic scales as populations differ in their exposition to varying time constraints. We examined differences in compensatory growth strategies in common frog (Rana temporaria) populations from southern and northern Sweden. Tadpoles from four populations were reared in the laboratory and exposed to low temperature to evaluate the patterns and mechanisms of compensatory growth responses. We determined tadpoles' growth rate, food intake and growth efficiency during the compensation period. In the absence of arrested growth conditions, tadpoles from all the populations showed similar (size-corrected) growth rates, food intake and growth efficiency. After being exposed to low temperature for 1 week, only larvae from the northern populations increased growth rates by increasing both food intake and growth efficiency. These geographic differences in compensatory growth mechanisms suggest that the strategies for recovering after a period of growth deprivation may depend on the strength of time constraints faced by the populations. Due to the costs of fast growth, only populations exposed to the strong time constraints are prone to develop fast recovering strategies in order to metamorphose before conditions deteriorate. Understanding how organisms balance the cost and benefits of growth strategies may help in forecasting the impact of fluctuating environmental conditions on life-history strategies of populations likely to be exposed to increasing environmental variation in the future.

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.

Using genetic variation to infer associations with climate in the common frog, Rana temporaria

Recent and historical species' associations with climate can be inferred using molecular markers. This knowledge of population and species-level responses to climatic variables can then be used to predict the potential consequences of ongoing climate change. The aim of this study was to predict responses of Rana temporaria to environmental change in Scotland by inferring historical and contemporary patterns of gene flow in relation to current variation in local thermal conditions. We first inferred colonization patterns within Europe following the last glacial maximum by combining new and previously published mitochondrial DNA sequences. We found that sequences from our Scottish samples were identical to (92%), or clustered with, the common haplotype previously identified from Western Europe. This clade showed very low mitochondrial variation, which did not allow inference of historical colonization routes but did allow interpretation of patterns of current fine-scale population structure without consideration of confounding historical variation. Second, we assessed fine-scale microsatellite-based patterns of genetic variation in relation to current altitudinal temperature gradients. No population structure was found within altitudinal gradients (average FST=0.02), despite a mean annual temperature difference of 4.5 °C between low- and high-altitude sites. Levels of genetic diversity were considerable and did not vary between sites. The panmictic population structure observed, even along temperature gradients, is a potentially positive sign for R. temporaria persistence in Scotland in the face of a changing climate. This study demonstrates that within taxonomic groups, thought to be at high risk from environmental change, levels of vulnerability can vary, even within species.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Connectivity in a pond system influences migration and genetic structure in threespine stickleback

Neutral genetic structure of natural populations is primarily influenced by migration (the movement of individuals and, subsequently, their genes) and drift (the statistical chance of losing genetic diversity over time). Migration between populations is influenced by several factors, including individual behavior, physical barriers, and environmental heterogeneity among populations. However, drift is expected to be stronger in populations with low immigration rate and small effective population size. With the technological advancement in geological information systems and spatial analysis tools, landscape genetics now allows the development of realistic migration models and increased insight to important processes influencing diversity of natural populations. In this study, we investigated the relationship between landscape connectivity and genetic distance of threespine stickleback (Gasterosteus aculeatus) inhabiting a pond complex in Belgjarskógur, Northeast Iceland. We used two landscape genetic approaches (i.e., least-cost-path and isolation-by-resistance) and asked whether gene flow, as measured by genetic distance, was more strongly associated with Euclidean distance (isolation-by-distance) or with landscape connectivity provided by areas prone to flooding (as indicated by Carex sp. cover)? We found substantial genetic structure across the study area, with pairwise genetic distances among populations (DPS) ranging from 0.118 to 0.488. Genetic distances among populations were more strongly correlated with least-cost-path and isolation-by-resistance than with Euclidean distance, whereas the relative contribution of isolation-by-resistance and Euclidian distance could not be disentangled. These results indicate that migration among stickleback populations occurs via periodically flooded areas. Overall, this study highlights the importance of transient landscape elements influencing migration and genetic structure of populations at small spatial scales.

Discordant patterns of nuclear and mitochondrial introgression in Iberian populations of the European common frog (Rana temporaria)

Amphibians often show complex histories of intraspecific and interspecific genetic introgression, which might differ in mitochondrial and nuclear genes. In our study of the genetic differentiation of the European common frog, Rana temporaria (159 specimens from 23 populations were analyzed for 24 presumptive allozyme loci; 82 specimens were sequenced for a 540-bp fragment of the mitochondrial 16S rRNA gene), multilocus correspondence analysis (CA) and Bayesian assignment tests of the nuclear data were concordant in identifying 2 population groups corresponding to 1) the Pyrenees in the east and 2) the Galicia and Asturias regions in the west, the latter corresponding to the subspecies R. temporaria parvipalmata. Geographically intermediate populations were genetically intermediate in the allozyme CA and, less clearly in the Bayesian assignment, with mitochondrial haplotypes exclusively belonging to the parvipalmata group. This indicates different degrees of introgression in the mitochondrial and nuclear genomes. Although Pyrenean high-altitude populations are morphologically distinct from low-altitude populations, these 2 groups were not separate clusters in any analysis. This suggests that the morphological differences may be due to fast adaptation to elevational gradients, likely under maintenance of gene flow, and that the underlying genetic changes are not detectable by the analyzed markers. We argue that a parsimonious explanation for the observed pattern along the east-west axis in northern Spain may be competition between invading and resident populations, with no need to invoke selection. However, in order to conclusively rule out selective processes, additional and finer scale data are required to test for asymmetric mating preference/behaviour, sex-biased gene flow, or sex-biased survival of potential hybrids.

Adaptive divergence in moor frog (Rana arvalis) populations along an acidification gradient: inferences from Q(st) -F(st) correlations

Microevolutionary responses to spatial variation in the environment seem ubiquitous, but the relative role of selection and neutral processes in driving phenotypic diversification remain often unknown. The moor frog (Rana arvalis) shows strong phenotypic divergence along an acidification gradient in Sweden. We here used correlations among population pairwise estimates of quantitative trait (P(ST) or Q(ST) from common garden estimates of embryonic acid tolerance and larval life-history traits) and neutral genetic divergence (F(ST) from neutral microsatellite markers), as well as environmental differences (pond pH, predator density, and latitude), to test whether this phenotypic divergence is more likely due to divergent selection or neutral processes. We found that trait divergence was more strongly correlated with environmental differences than the neutral marker divergence, suggesting that divergent natural selection has driven phenotypic divergence along the acidification gradient. Moreover, pairwise P(ST) s of embryonic acid tolerance and Q(ST) s of metamorphic size were strongly correlated with breeding pond pH, whereas pairwise Q(ST) s of larval period and growth rate were more strongly correlated with geographic distance/latitude and predator density, respectively. We suggest that incorporating measurements of environmental variation into Q(ST) -F(ST) studies can improve our inferential power about the agents of natural selection in natural populations.

Comparisons between Q(ST) and F(ST) --how wrong have we been?

The comparison between quantitative genetic divergence (Q(ST) ) and neutral genetic divergence (F(ST) ) among populations has become the standard test for historical signatures of selection on quantitative traits. However, when the mutation rate of neutral markers is relatively high in comparison with gene flow, estimates of F(ST) will decrease, resulting in upwardly biased comparisons of Q(ST) vs. F(ST) . Reviewing empirical studies, the difference between Q(ST) and F(ST) is positively related to marker heterozygosity. After refuting alternative explanations for this pattern, we conclude that marker mutation rate indeed has had a biasing effect on published Q(ST) -F(ST) comparisons. Hence, it is no longer clear that populations have commonly diverged in response to divergent selection. We present and discuss potential solutions to this bias. Comparing Q(ST) with recent indices of neutral divergence that statistically correct for marker heterozygosity (Hedrick's G'st and Jost's D) is not advised, because these indices are not theoretically equivalent to Q(ST) . One valid solution is to estimate F(ST) from neutral markers with mutation rates comparable to those of the loci underlying quantitative traits (e.g. SNPs). Q(ST) can also be compared to Φ(ST) (Phi(ST) ) of amova, as long as the genetic distance among allelic variants used to estimate Φ(ST) reflects evolutionary history: in that case, neutral divergence is independent of mutation rate. In contrast to their common usage in comparisons of Q(ST) and F(ST) , microsatellites typically have high mutation rates and do not evolve according to a simple evolutionary model, so are best avoided in Q(ST) -F(ST) comparisons.

Conserving adaptive genetic diversity in dynamic landscapes

Genetic variation supplies the raw material for adaptation, evolution and survival of populations and has therefore been a key focus of conservation biology ever since its foundation (Soulé 1985). In previous decades, the neutral component of genetic diversity (generated by mutation and shaped by drift) has been the subject of intense scientific research, fuelled by the increasing availability of molecular markers. On the other hand, the adaptive component of genetic diversity, which is shaped by the action of natural selection, has long remained elusive and difficult to assess, especially at small spatial or temporal scales (Ouborg et al. 2010). Fortunately, new technological and methodological developments now make it possible to identify loci in the genome that are influenced by selection, and thus to get a more complete view of genetic diversity. One article featured in this issue of Molecular Ecology is a good example of this recent breakthrough. Richter-Boix et al. (2011) examined a network of moor frog populations breeding in contrasting habitats in order to understand how landscape features influence patterns of genetic variation. They combined information from both neutral markers and loci putatively under selection to quantify the relative roles of selection and isolation in the evolution of fine-scale local adaptations in these populations. This study nicely illustrates how data on polymorphisms of neutral and adaptive loci can now be judiciously synthesized to help identify the best strategies for preserving adaptive variation, and more generally to enlighten conservation and population-management plans.