Variation in the degree and costs of adaptive phenotypic plasticity among Rana temporaria populations

Adaptive phenotypic plasticity in the form of capacity to accelerate development as a response to pond drying risk is known from many amphibian species. However, very little is known about factors that might constrain the evolution of this type of plasticity, and few studies have explored to what degree plasticity might be constrained by trade-offs dictated by adaptation to different environmental conditions. We compared the ability of southern and northern Scandinavian common frog (Rana temporaria) larvae originating from 10 different populations to accelerate their development in response to simulated pond drying risk and the resulting costs in metamorphic size in a factorial laboratory experiment. We found that (i) northern larvae developed faster than the southern larvae in all treatments, (ii) a capacity to accelerate the response was present in all five southern and all five northern populations tested, but that the magnitude of the response was much larger (and less variable) in the southern than in the northern populations, and that (iii) significant plasticity costs in metamorphic size were present in the southern populations, the plastic genotypes having smaller metamorphic size in the absence of desiccation risk, but no evidence for plasticity costs was found in the northern populations. We suggest that the weaker response to pond drying risk in the northern populations is due to stronger selection on large metamorphic size as compared with southern populations. In other words, seasonal time constraints that have selected the northern larvae to be fast growing and developing, may also constrain their innate ability for adaptive phenotypic plasticity.     

The degree of adaptive phenotypic plasticity is correlated with the spatial environmental heterogeneity experienced by island populations of Rana temporaria

Although theoretical models have identified environmental heterogeneity as a prerequisite for the evolution of adaptive plasticity, this relationship has not yet been demonstrated experimentally. Because of pool desiccation risk, adaptation of development rate is important for many amphibians. In a simulated pool-drying experiment, we compared the development time and phenotypic plasticity in development time of populations of the common frog Rana temporaria, originating from 14 neighbouring islands off the coast of northern Sweden. Drying regime of pools used by frogs for breeding differed within and among the islands. We found that the degree of phenotypic plasticity in development time was positively correlated with the spatial variation in the pool-drying regimes present on each island. In addition, local adaptation in development time to the mean drying rate of the pools on each island was found. Hence, our study demonstrates the connection between environmental heterogeneity and developmental plasticity at the island population level, and also highlights the importance of the interplay between local specialization and phenotypic plasticity depending on the local selection pressures.     

Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow

Divergent natural selection, adaptive divergence and gene flow may interact in a number of ways. Recent studies have focused on the balance between selection and gene flow in natural populations, and empirical work has shown that gene flow can constrain adaptive divergence, and that divergent selection can constrain gene flow. A caveat is that phenotypic diversification may be under the direct influence of environmental factors (i.e. it may be due to phenotypic plasticity), in addition to partial genetic influence. In this case, phenotypic divergence may occur between populations despite high gene flow that imposes a constraint on genetic divergence. Plasticity may dampen the effects of natural selection by allowing individuals to rapidly adapt phenotypically to new conditions, thus slowing adaptive genetic divergence. On the other hand, plasticity may promote future adaptive divergence by allowing populations to persist in novel environments. Plasticity may promote gene flow between selective regimes by allowing dispersers to adapt to alternate conditions, or high gene flow may result in the selection for increased plasticity. Here I expand frameworks for understanding relationships among selection, adaptation and gene flow to include the effects of phenotypic plasticity in natural populations, and highlight its importance in evolutionary diversification.     

Proximate causes of adaptive growth rates: growth efficiency variation among latitudinal populations of Rana temporaria

In ectothermic organisms, declining season length and lower temperature towards higher latitudes often select for latitudinal variation in growth and development. However, the energetic mechanisms underlying this adaptive variation are largely unknown. We investigated growth, food intake and growth efficiency of Rana temporaria tadpoles from eight populations along a 1500 km latitudinal gradient across Sweden. To gain an insight into the mechanisms of adaptation at organ level, we also examined variation in tadpole gut length. The tadpoles were raised at two temperatures (16 and 20 degrees C) in a laboratory common garden experiment. We found increased growth rate towards higher latitudes, regardless of temperature treatment. This increase in growth was not because of a higher food intake rate, but populations from higher latitudes had higher growth efficiency, i.e. they were more efficient at converting ingested food into body mass. Low temperature reduced growth efficiency most strongly in southern populations. Relative gut length increased with latitude, and tadpoles at low temperature tended to have longer guts. However, variation in gut length was not the sole adaptive explanation for increased growth efficiency as latitude and body length still explained significant amounts of variation in growth efficiency. Hence, additional energetic adaptations are probably involved in growth efficiency variation along the latitudinal gradient.     

Habitat specialization and adaptive phenotypic divergence of anuran populations

We tested for adaptive population structure in the frog Rana temporaria by rearing tadpoles from 23 populations in a common garden experiment, with and without larval dragonfly predators. The goal was to compare tadpole phenotypes with the habitats of their source ponds. The choice of traits and habitat variables was guided by prior information about phenotypic function. There were large differences among populations in life history, behaviour, morphological shape, and the predator-induced plasticities in most of these. Body size and behaviour were correlated with predation risk in the source pond, in agreement with adaptive population divergence. Tadpoles from large sunny ponds were morphologically distinct from those inhabiting small woodland ponds, although here an adaptive explanation was unclear. There was no evidence that plasticity evolves in populations exposed to more variable environments. Much among-population variation in phenotype and plasticity was not associated with habitat, perhaps reflecting rapid changes in wetland habitats.     

Adaptive phenotypic plasticity and genetics of larval life histories in two Rana temporaria populations

Phenotypic plasticity provides means for adapting to environmental unpredictability. In terms of accelerated development in the face of pond-drying risk, phenotypic plasticity has been demonstrated in many amphibian species, but two issues of evolutionary interest remain unexplored. First, the heritable basis of plastic responses is poorly established. Second, it is not known whether interpopulational differences in capacity to respond to pond-drying risk exist, although such differences, when matched with differences in desiccation risk would provide strong evidence for local adaptation. We investigated sources of within- and among-population variation in plastic responses to simulated pond-drying risk (three desiccation treatments) in two Rana temporaria populations originating from contrasting environments: (1) high desiccation risk with weak seasonal time constraint (southern population); and (2) low desiccation risk with severe seasonal time constraint (northern population). The larvae originating from the environment with high desiccation risk responded adaptively to the fast decreasing water treatment by accelerating their development and metamorphosing earlier, but this was not the case in the larvae originating from the environment with low desiccation risk. In both populations, metamorphic size was smaller in the high-desiccation-risk treatment, but the effect was larger in the southern population. Significant additive genetic variation in development rate was found in the northern and was nearly significant in the southern population, but there was no evidence for genetic variation in plasticity for development rates in either of the populations. No genetic variation for plasticity was found either in size at metamorphosis or growth rate. All metamorphic traits were heritable, and additive genetic variances were generally somewhat higher in the southern population, although significantly so in only one trait. Dominance variances were also significant in three of four traits, but the populations did not differ. Maternal effects in metamorphic traits were generally weak in both populations. Within-environment phenotypic correlations between larval period and metamorphic size were positive and genetic correlations negative in both populations. These results suggest that adaptive phenotypic plasticity is not a species-specific fixed trait, but evolution of interpopulational differences in plastic responses are possible, although heritability of plasticity appears to be low. The lack of adaptive response to desiccation risk in northern larvae is consistent with the interpretation that selection imposed by shorter growing season has favored rapid development in north (approximately 8% faster development in north as compared to south) or a minimum metamorphic size at the expense of phenotypic plasticity.     

Quantitative genetic analysis of larval life history traits in two alpine populations of Rana temporaria

We estimated genetic and maternal variance components of larval life history characters in alpine populations of Rana temporaria (the common frog) using a full-sib/half-sib breeding design. We studied trait plasticity by raising tadpoles at 14 or 20°C in the laboratory. Larval period and metamorphic mass were greater at 14°C. Larval period did not differ between populations, but high elevation metamorphs were larger than low elevation metamorphs. Significant additive variation for larval period was detected in the low altitude population. No significant additive variation was detected for mass at metamorphosis (MM), which instead displayed significant maternal effects. Plasticity in metamorphic mass of froglets was greater in the high altitude population. The plastic response of larval period to temperature did not differ between the populations. Evolution of metamorphic mass is likely constrained by lack of additive genetic variation. In contrast, significant heritability for larval period suggests this trait may evolve in response to environmental change. These results differ from other studies on R. temporaria, suggesting that populations of this broadly distributed species present substantial geographic variation in the genetic architecture and plasticity of tadpole life history traits.     

Population genetic structure and a possible role for selection in driving phenotypic divergence in a rainbowfish (Melanotaeniidae)

Quantitative genetic divergence may be driven by drift or selection. The rainbowfish Melanotaenia australis exhibits phenotypic divergence among populations in Western Australia, although the mechanisms driving this divergence are unknown. We used microsatellites to assess neutral genetic divergence (F_ST), and found a hierarchical pattern of subdivision consistent with low divergence between upstream and downstream populations (within drainages), moderate divergence between drainages (within regions), and high divergence between regions. Using a common‐garden approach, we measured quantitative genetic divergence in phenotypic traits (Q_ST). By comparing this to expectations from neutral processes (F_ST), we concluded that the effect of selection varies depending on the spatial scale considered. Within drainages, selection may be causing divergence between upstream and downstream phenotypes but, between regions, selection appears to homogenize phenotypes. This highlights the importance of spatial scale in studies of this kind, and suggests that, because variance in selection pressures can drive speciation, polymorphism in M. australis may represent speciation in action. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 144–160.     

Latitudinal countergradient variation in the common frog (Rana temporaria) development rates--evidence for local adaptation

Adaptive genetic differentiation along a climatic gradient as a response to natural selection is not necessarily expressed at phenotypic level if environmental effects on population mean phenotypes oppose the genotypic effects. This form of cryptic evolution--called countergradient variation--has seldom been explicitly demonstrated for terrestrial vertebrates. We investigated the patterns of phenotypic and genotypic differentiation in developmental rates of common frogs (Rana temporaria) along a ca. 1600 km latitudinal gradient across Scandinavia. Developmental rates in the field were not latitudinally ordered, but displayed large variation even among different ponds within a given latitudinal area. In contrast, development rates assessed in the laboratory increased strongly and linearly with increasing latitude, suggesting a genetic capacity for faster development in the northern than the southern larvae. Experiments further revealed that environmental effects (temperature and food) could easily override the genetic effects on developmental rates, providing a possible mechanistic explanation as to why the genetic differentiation was not seen in the samples collected from the wild. Our results suggest that the higher developmental rates of the northern larvae are likely to be related to selection stemming from seasonal time constrains, rather than from selection dictated by low ambient temperatures per se. All in all, the results provide a demonstration of environmental effects concealing substantial latitudinally ordered genetic differentiation understandable in terms of adaptation to clinal variation in time constrains.     

Genetic and environmental components of phenotypic variation in body shape among populations of Atlantic cod (Gadus morhua L.)

A common-garden experiment was conducted on larvae to test for genetic differences in body shape among populations of Atlantic cod ( Gadus morhua ). Offspring from four north-west Atlantic regions were reared from hatching to postmetamorphosis at two temperatures (7 ± 1 °C and 11 ± 1 °C) and two food levels (1500 and 4500 prey L⁻¹). Body shape differed between populations and treatments. Population differences were greatest between south-west Scotian Shelf cod and those further north; the former were characterized by a deeper body, larger head, and longer caudal peduncle than cod from the other populations. Significant differences were also observed between two putative populations on the south-west Scotian Shelf, suggesting genetic divergence between spawning aggregations at small spatial scales (< 100 km). Temperature and food supply also influenced body shape, with the effect of the former being more pronounced. Individuals reared at the higher temperature or food level had a deeper body and a larger head than those reared at the lower temperature or food supply. Phenotypic responses to changes in the rearing environment also differed among populations, indicating genetic differences in phenotypic plasticity. Differences between populations in morphology and in phenotypic plasticity suggest genetic divergence at both large (> 1000 km) and small (< 100 km) spatial scales. The genetic differences at large spatial scales counteracted the expected effects of temperature differences in the wild, suggesting countergradient variation in morphology among these populations.  © 2006 Her Majesty the Queen in Right of Canada. Journal compilation © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 351–365.     

Maternal and genetic contributions to geographical variation in Rana temporaria larval life-history traits

The relative importance of genetic, environmental, and maternal effects as determinants of geographical variation in vertebrate life-histories has not often been explored. We examined the role of genetic and maternal effects as determinants of population divergence in survival and three important larval life-history traits (growth rate, age, and size at metamorphosis) using reciprocal crosses between two latitudinally separated populations of the common frog ( Rana temporaria Linnaeus). Genetic effects were important in all three traits as indicated by the significant effect of male origin, but there was also evidence for nonadditive genetic contributions on metamorphic size and growth rate. Likewise, maternal effect contributions to population divergence were large, partially environment dependent, and apparently acting primarily through egg size in two of three traits. These results suggest that both genetic and maternal effects are important determinants of geographical variation in amphibian life-histories, and that much of the differentiation resulting from maternal effects is mediated through variation in egg size. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 76 , 61–70.     

Phenotypic plasticity in gene expression contributes to divergence of locally adapted populations of Fundulus heteroclitus

We examine the interaction between phenotypic plasticity and evolutionary adaptation using muscle gene expression levels among populations of the fish Fundulus heteroclitus acclimated to three temperatures. Our analysis reveals shared patterns of phenotypic plasticity due to thermal acclimation as well as non‐neutral patterns of variation among populations adapted to different thermal environments. For the majority of significant differences in gene expression levels, phenotypic plasticity and adaptation operate on different suites of genes. The subset of genes that demonstrate both adaptive differences and phenotypic plasticity, however, exhibit countergradient variation of expression. Thus, expression differences among populations counteract environmental effects, reducing the phenotypic differentiation between populations. Finally, gene‐by‐environment interactions among genes with non‐neutral patterns of expression suggest that the penetrance of adaptive variation depends on the environmental conditions experienced by the individual.     

Local adaptation and cogradient selection in the alpine plant, Poa hiemata, along a narrow altitudinal gradient

Alpine environments are particularly susceptible to environmental changes associated with global warming but there is potential for alpine plants to adapt to warming if local adaptation occurs and gene flow allows genotypes adapted to low altitudes to colonize higher altitude sites. Here we examine the adaptive potential of a common alpine grass, Poa hiemata, within the restricted alpine habitat of Australian mountains, across a narrow altitudinal gradient replicated in three areas. Grasses at high altitude sites had shorter leaf lengths and larger circumferences than those at lower sites. Transplant experiments with clonal material and plants grown from seed indicated that these differences were partly genetic, with environmental and genetic factors both contributing to the differences between altitudes. Differences in altitudinal forms were also evident in a common garden experiment. Plants showed a home-site advantage in terms of survival. A fitness analysis indicated that at high altitude sites, selection favored plants with short leaves and larger circumferences, whereas these traits were selected in the opposite direction at the low altitude sites. These findings indicate cogradient selection and potential for both plastic and genotypic shifts in response to climate change in P. hiemata.     

Interpopulation variation in a fish predator drives evolutionary divergence in prey in lakes

Ecological factors are known to cause evolutionary diversification. Recent work has shown that evolution in strongly interacting predator species has reciprocal impacts on ecosystems. These divergent impacts of predators may alter the selective landscape and cause the evolution of prey. Yet, this link between intraspecific variation and evolution is unexplored. We compared the life history of a species of zooplankton (Daphnia ambigua) from lakes in New England in which the dominant planktivorous predator, the alewife (Alosa pseudoharengus), differs in feeding traits and migratory behaviour. Anadromous alewife (seasonal migrants) exhibit larger gapes, gill-raker spacing and target larger prey than landlocked alewife (year-round freshwater resident). In 'anadromous' lakes, Daphnia are abundant in the spring but extirpated by alewife predation in summer. Daphnia are rare year-round in 'landlocked' lakes. We show that Daphnia from lakes with anadromous alewife grew faster, matured earlier but at the same size and produced more offspring than Daphnia from lakes with landlocked or no alewife across multiple temperature and resource treatments. Our results are inconsistent with a response to size-selective predation but are better explained as an adaptation to colder temperatures and shorter periods of development (countergradient variation) mediated by seasonal alewife predation.     

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.     

Population differentiation in G matrix structure due to natural selection in Rana temporaria

The additive genetic variance-covariance matrix (G) is a concept central to discussions about evolutionary change over time in a suite of traits. However, at the moment we do not know how fast G itself changes as a consequence of selection or how sensitive it is to environmental influences. We investigated possible evolutionary divergence and environmental influences on G using data from a factorial common-garden experiment where common frog (Rana temporaria) tadpoles from two divergent populations were exposed to three different environmental treatments. G-matrices were estimated using an animal model approach applied to data from a NCII breeding design. Matrix comparisons using both Flury and multivariate analysis of variance methods revealed significant differences in G matrices both between populations and between treatments within populations, the former being generally larger than the latter. Comparison of levels of population differentiation in trait means using Q(ST) indices with that observed in microsatellite markers (F(ST)) revealed that the former values generally exceeded the neutral expectation set by F(ST). Hence, the results suggest that intraspecific divergence in G matrix structure has occurred mainly due to natural selection.     

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

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

Establishment and maintenance of adaptive genetic divergence under migration, selection, and drift

There is a long tradition in population genetics of exploring the maintenance of variation under migration-selection balance using deterministic models that assume infinite population size. With finite population size, stochastic dynamics can greatly reduce the potential for the maintenance of polymorphism, but this has yet to be explored in detail. Here, classical two-patch models are extended to predict: (1) the probability of a locally beneficial mutation rising in frequency in the patch where it is favored and (2) the critical threshold migration rate above which the maintenance of polymorphism is much less likely. Individual-based simulations show that these approximations provide accurate predictions across a wide range of parameter space.