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.     

Low neutral and immunogenetic diversity in northern fringe populations of the green toad Bufotes viridis: implications for conservation

Genetic variation is often lower at high latitudes, which may compromise the adaptability and hence survival of organisms. Here we show that genetic variability is negatively correlated with northern latitude in European green toads (Bufotes viridis). The result holds true for both putatively neutral microsatellite variation and supposedly adaptive MHC Class IIB variation. In particular, our findings have bearing on the conservation status of this species in Sweden, on the northern limit of its distribution where local populations are small and fragmented. These genetically impoverished populations are closely related to other populations found around the Baltic Sea basin. The low neutral and adaptive variation in these fringe populations compared to population at central ranges confirms a pattern shared across all other amphibians so far studied. In Sweden, the situation of green toads is of concern as the remaining populations may not have the evolutionary potential to cope with present and future environmental challenges.

     

ADAPTIVE GENETIC VARIATION IN GROWTH AND DEVELOPMENT OF TADPOLES OF THE HYBRIDOGENETIC RANA ESCULENTA COMPLEX

The distribution and proportion of the sexual species Rana lessonae to the hemiclonal hybrid R. esculenta among natural habitats suggests that these anurans may differ in adaptive abilities. I used a half-sib design to partition phenotypic and quantitative genetic variation in tadpole responses at two food levels into causal variance components. Rana lessonae displays strong phenotypic variation across food levels. Growth rate is strictly determined by environmental factors and includes weak maternal effects. Larval period and body size at metamorphosis both contain moderate levels of additive genetic variance. The sire x food interactions and the lack of environmental correlations indicate that adaptive phenotypic plasticity is present in both of these traits. In contrast, R. esculenta displays less phenotypic variation across food levels, especially for larval period. Variation in body size at metamorphosis is underlain by genetic variation as shown by high levels of additive genetic variance, yet growth rate and larval period are not. Significant environmental correlations between larval period at high food level and growth, larval period, and body size at low food, indicate phenotypic plasticity is absent. A positive phenotypic correlation between body size at metamorphosis and larval period for R. lessonae at both food levels suggests a trade-off between growing large and metamorphosing quickly to escape predation or pond drying. The lack of a similar correlation for R. esculenta at the high food level suggests it may be less constrained. Different levels of adaptive genetic variation among larval traits suggest that the sexual species and the hybridogenetic hemiclone differ in their abilities to cope with temporally and spatially heterogeneous environments.     

Testing local and global stressor impacts on a coastal foundation species using an ecologically realistic framework

Despite the abundance of literature on organismal responses to multiple environmental stressors, most studies have not matched the timing of experimental manipulations with the temporal pattern of stressors in nature. We test the interactive effects of diel‐cycling hypoxia with both warming and decreased salinities using ecologically realistic exposures. Surprisingly, we found no evidence of negative synergistic effects on Olympia oyster growth; rather, we found only additive and opposing effects of hypoxia (detrimental) and warming (beneficial). We suspect that diel‐cycling provided a temporal refuge that allowed physiological compensation. We also tested for latent effects of warming and hypoxia to low‐salinity tolerance using a seasonal delay between stressor events. However, we did not find a latent effect, rather a threshold survival response to low salinity that was independent of early life‐history exposure to warming or hypoxia. The absence of synergism is likely the result of stressor treatments that mirror the natural timing of environmental stressors. We provide environmental context for laboratory experimental data by examining field time series environmental data from four North American west coast estuaries and find heterogeneous environmental signals that characterize each estuary, suggesting that the potential stressor exposure to oysters will drastically differ over moderate spatial scales. This heterogeneity implies that efforts to conserve and restore oysters will require an adaptive approach that incorporates knowledge of local conditions. We conclude that studies of multiple environmental stressors can be greatly improved by integrating ecologically realistic exposure and timing of stressors found in nature with organismal life‐history traits.     

Fitness costs associated with low genetic variation are reduced in a harsher environment in amphibian island populations

A basic premise of conservation geneticists is that low levels of genetic variation are associated with fitness costs in terms of reduced survival and fecundity. These fitness costs may frequently vary with environmental factors and should increase under more stressful conditions. However, there is no consensus on how fitness costs associated with low genetic variation change under natural conditions in relation to the stressfulness of the environment. On the Swedish west coast, natterjack toad Bufo calamita populations show a strong population genetic structure and large variation in the amount of within-population genetic variation. We experimentally examined the survival of natterjack larvae from six populations with different genetic variation in three thermal environments corresponding to (a) the mean temperature of natural ponds (stable, laboratory), (b) a high temperature environment occurring in desiccating ponds (stable, laboratory) and (c) an outdoor treatment mimicking the natural, variable thermal conditions (fluctuating, semi-natural). We found that larvae in the outdoor treatment had poorer survival than larvae in the stable environments suggesting that the outdoor treatment was more stressful. Overall, populations with higher genetic variation had higher larval survival. However, a significant interaction between treatments and genetic variation indicated that fitness costs associated with low genetic variation were less severe in the outdoor treatment. Thus, we found no support for the hypothesis that fitness costs associated with low genetic variation increase under more stressful conditions. Our results suggest that natural thermal stress may mask fitness losses associated with low genetic variation in these populations.     

The effects of selection,drift and genetic variation on life‐history trait divergence among insular populations of natterjack toad,Bufo calamita

Although loss of genetic variation is frequently assumed to be associated with loss of adaptive potential, only few studies have examined adaptation in populations with little genetic variation. On the Swedish west coast, the northern fringe populations of the natterjack toad Bufo calamita inhabit an atypical habitat consisting of offshore rock islands. There are strong among‐population differences in the amount of neutral genetic variation, making this system suitable for studies on mechanisms of trait divergence along a gradient of within‐population genetic variation. In this study, we examined the mechanisms of population divergence using Q_ST–F_ST comparisons and correlations between quantitative and neutral genetic variation. Our results suggest drift or weak stabilizing selection across the six populations included in this study, as indicated by low Q_ST–F_ST values, lack of significant population × temperature interactions and lack of significant differences among the islands in breeding pond size. The six populations included in this study differed in both neutral and quantitative genetic variation. Also, the correlations between neutral and quantitative genetic variation tended to be positive, however, the relatively small number of populations prevents any strong conclusions based on these correlations. Contrary to the majority of Q_ST–F_ST comparisons, our results suggest drift or weak stabilizing selection across the examined populations. Furthermore, the low heritability of fitness‐related traits may limit evolutionary responses in some of the populations.     

Contrasting patterns of quantitative and neutral genetic variation in locally adapted populations of the natterjack toad, Bufo calamita

The relative importance of natural selection and genetic drift in determining patterns of phenotypic diversity observed in nature is still unclear. The natterjack toad (Bufo calamita) is one of a few amphibian species capable of breeding in saline ponds, even though water salinity represents a considerable stress for them. Results from two common-garden experiments showed a pattern of geographic variation in embryonic salinity tolerance among populations from either fresh or brackish environments, consistent with the hypothesis of local adaptation. Full-sib analysis showed increased variation in survival among sibships within population for all populations as osmotic stress was increased (broad-sense heritability increased as salinity raised). Nevertheless, toads native to the brackish water environment had the highest overall survival under brackish conditions. Levels of population genetic differentiation for salinity tolerance were higher than those of neutral genetic differentiation, the latter obtained through the analysis of eight microsatellite loci. Microsatellite markers also revealed little population differentiation, lack of an isolation-by-distance pattern, and moderate gene flow connecting the populations. Therefore, environmental stress tolerance appears to have evolved in absence of geographic isolation, and consequently we reject the null hypothesis of neutral differentiation.     

Influence of density and salinity on larval development of salt‐adapted and salt‐naïve frog populations

Environmental change and habitat fragmentation will affect population densities for many species. For those species that have locally adapted to persist in changed or stressful habitats, it is uncertain how density dependence will affect adaptive responses. Anurans (frogs and toads) are typically freshwater organisms, but some coastal populations of green treefrogs (Hyla cinerea) have adapted to brackish, coastal wetlands. Tadpoles from coastal populations metamorphose sooner and demonstrate faster growth rates than inland populations when reared solitarily. Although saltwater exposure has adaptively reduced the duration of the larval period for coastal populations, increases in densities during larval development typically increase time to metamorphosis and reduce rates of growth and survival. We test how combined stressors of density and salinity affect larval development between salt‐adapted (“coastal”) and nonsalt‐adapted (“inland”) populations by measuring various developmental and metamorphic phenotypes. We found that increased tadpole density strongly affected coastal and inland tadpole populations similarly. In high‐density treatments, both coastal and inland populations had reduced growth rates, greater exponential decay of growth, a smaller size at metamorphosis, took longer to reach metamorphosis, and had lower survivorship at metamorphosis. Salinity only exaggerated the effects of density on the time to reach metamorphosis and exponential decay of growth. Location of origin affected length at metamorphosis, with coastal tadpoles metamorphosing slightly longer than inland tadpoles across densities and salinities. These findings confirm that density has a strong and central influence on larval development even across divergent populations and habitat types and may mitigate the expression (and therefore detection) of locally adapted phenotypes.     

Microgeographic variation in temperature-induced plasticity in an isolated amphibian metapopulation

Variation in local environments may lead to variation in the selection pressures and differentiation among local populations even at microgeographic scale. We investigated variation in temperature-induced plasticity in larval life-history traits among populations of an isolated pool frog (Rana lessonae) metapopulation in Central Sweden. Successful breeding of this northern fringe metapopulation is highly dependent on early summer temperature, however, the metapopulation shows very little variation in molecular genetic markers suggesting limited potential for local differentiation. We exposed larvae from three closely-located populations to two temperatures (20 and 25°C) in laboratory to investigate their growth and development responses to temperature variation. In general, larvae exposed to warmer temperature experienced higher survival and metamorphosed faster, but at a smaller size than those at low temperature. We found differences among the populations in both trait mean values and in the plastic responses. Among-family variation within populations was found in growth rate and time to metamorphosis, as well as in plasticity suggesting that these traits have a capacity to evolve. Our results indicate ample phenotypic variation within and among these closely-located populations despite the low molecular genetic variation. The differences in pond temperature characteristics detected in the study in the three localities may suggest that differential selection is acting in the populations. The strong differentiation found in the larval traits implies that understanding the factors that influence the potential of the populations to adapt to environmental changes may be essential for successful conservation strategies.     

Genetic structure in peripheral populations of the natterjack toad, <Emphasis Type="Italic">Bufo calamita</Emphasis>, as revealed by AFLP

Decreased fitness due to loss of genetic variation is a well recognised issue in conservation biology. Along the Swedish west coast, the endangered natterjack toad (Bufo calamita) occurs on, for the species, highly unusual habitat of rocky islands. Although the toads inhabit a restricted geographical area (maximum distance between the populations is 71 km), the fragmented nature of the landscape makes the genetic properties of the populations of conservation interest. However, lack of genetic variation found using conventional methods (microsatellites) has impeded genetic studies within these peripheral populations so far. In this study we assess population structure and genetic variation among seven of these fringe populations using 105 polymorphic Amplified Fragment Length Polymorphism (AFLP) loci. We found a well-defined population structure without evidence for isolation by distance, implying restricted gene flow between populations. Additionally, the populations differed in their amount of genetic variation, emphasizing the need to monitor genetically impoverished populations for possible declines mediated by inbreeding depression and reduced adaptive potential. Conservation implications for these unique populations are discussed in the light of our results.     

Biocomplexity in Populations of European Anchovy in the Adriatic Sea

The sustained exploitation of marine populations requires an understanding of a species'' adaptive seascape so that populations can track environmental changes from short- and long-term climate cycles and from human development. The analysis of the distributions of genetic markers among populations, together with correlates of life-history and environmental variability, can provide insights into the extent of adaptive variation. Here, we examined genetic variability among populations of mature European anchovies (n = 531) in the Adriatic (13 samples) and Tyrrhenian seas (2 samples) with neutral and putative non-neutral microsatellite loci. These genetic markers failed to confirm the occurrence of two anchovy species in the Adriatic Sea, as previously postulated. However, we found fine-scale population structure in the Adriatic, especially in northern areas, that was associated with four of the 13 environmental variables tested. Geographic gradients in sea temperature, salinity and dissolved oxygen appear to drive adaptive differences in spawning time and early larval development among populations. Resolving adaptive seascapes in Adriatic anchovies provides a means to understand mechanisms underpinning local adaptation and a basis for optimizing exploitation strategies for sustainable harvests.     

Antiporter NHX2 differentially induced in Mesembryanthemum crystallinum natural genetic variant under salt stress

Plants exhibit several mechanisms to survive under high salinity conditions. The uptake and compartmentalization of Na⁺ ion by the NHX antiporter is a crucial mechanism in homeostasis maintenance. Therefore, we evaluated McNHX2 gene expression and several physiological responses induced in three natural genetic variants of ice plants under salt stress. Based on morphology and growth behavior of wild type populations from an arid region of northwestern Mexico, we identified three ice plant natural genetic variants and called P0, P9, and P11. Several physiological parameters, such as water potential, relative water content, chlorophyll, and Na⁺ and K⁺ ion contents from all natural genetic variants exhibited a differential response under high salinity conditions. Specifically, the P0 variant showed lower water potential changes and least perturbation of Na⁺/K⁺ ratio than those of the P9 and P11 variants under saline conditions, suggesting that the P0 variant is the most salt tolerant. Unexpectedly, McNHX2 expression was repressed in the P11 variant while it was upregulated in the P0 and P9 variants under saline treatments. The McNHX2 gene was sequenced showing 15 introns and 16 exons; polymorphisms were found among the cDNAs sequences from the three natural genetic variants. All these data suggest that differential responses to high salinity involve different mechanisms operating in each variant for counteracting saline stress effects.

     

Epigenetic variation associated with responses to different habitats in the context of genetic divergence in Phragmites australis

The mechanisms underlying heritable phenotypic divergence associated with adaptation in response to environmental stresses may involve both genetic and epigenetic variations. Several prior studies have revealed even higher levels of epigenetic variation than genetic variation. However, few population‐level studies have explored the effects of epigenetic variation on species with high levels of genetic diversity distributed across different habitats. Using AFLP and methylation‐sensitive AFLP markers, we tested the hypothesis that epigenetic variation may contribute to differences in plants occupying different habitats when genetic variation alone cannot fully explain adaptation. As a cosmopolitan invasive species, Phragmites australis (common reed) together with high genetic diversity and remarkable adaptability has been suggested as a model for responses to global change and indicators of environmental fluctuations. We found high levels of genetic and epigenetic diversity and significant genetic/epigenetic structure within each of 12 studied populations sampled from four natural habitats of P. australis. Possible adaptive epigenetic variation was suggested by significant correlations between DNA methylation‐based epigenetic differentiation and adaptive genetic divergence in populations across the habitats. Meanwhile, various AMOVAs indicated that some epigenetic differences may respond to various local habitats. A partial Mantel test was used to tease out the correlations between genetic/epigenetic variation and habitat after controlling for the correlation between genetic and epigenetic variations. We found that epigenetic diversity was affected mostly by soil nutrient availability, suggesting that at least some epigenetic differentiation occurred independently of genetic variation. We also found stronger correlations between epigenetic variation and phenotypic traits than between genetic variation and such traits. Overall, our findings indicate that genetically based differentiation correlates with heterogeneous habitats, while epigenetic variation plays an important role in ecological differentiation in natural populations of P. australis. In addition, our results suggest that when assessing global change responses of plant species, intraspecific variation needs to be considered.     

Genetic variation in life-history reaction norms in a marine fish

Neither the scale of adaptive variation nor the genetic basis for differential population responses to the environment is known for broadcast-spawning marine fishes. Using a common-garden experimental protocol, we document how larval growth, survival and their norms of reaction differ genetically among four populations of Atlantic cod (Gadus morhua). These traits, and their plastic responses to food and temperature, differed across spatial scales at which microsatellite DNA failed to detect population structure. Divergent survival reaction norms indicate that warm-water populations are more sensitive to changes in food, whereas cold-water populations are more sensitive to changes in temperature. Our results suggest that neither the direction nor the magnitude of demographic responses to environmental change need be the same among populations. Adaptive phenotypic plasticity, previously undocumented in marine fishes, can significantly influence the probability of recovery and persistence of collapsed populations by affecting their ability to respond to natural and anthropogenic environmental change.     

Population structure of flounder (Platichthys flesus) in the Baltic Sea: differences among demersal and pelagic spawners

We found significant population structure and isolation by distance among samples of flounder (Platichthys flesus) in the Baltic, Kattegat and Skagerrak seas using microsatellite genetic markers. This pattern was almost entirely due to a difference between flounder that have demersal spawning in the northern Baltic, as compared to pelagic spawners in the southern Baltic and on the west coast of Sweden. Among demersal spawners we found neither genetic differentiation nor any isolation by distance among sampling sites. We speculate that demersal flounder are descendants of a population that colonized the Baltic previous to pelagic spawners. The demersal flounder may thus have had longer time to adapt to the low salinity in the Baltic, and accordingly display egg characteristics that make it possible to reproduce at the low salinity levels in the northern Baltic. Among pelagic spawners significant isolation by distance was detected. Pelagic spawners have previously been shown to display clinal variation in egg size, which allows them to float also at the moderate salinity levels up to the region north of the island Bornholm. Management units for harvesting should ideally be based on true biological populations, and for the commercially important flounder up to 15 different management stocks in the Baltic have been suggested. We could not find a population genetic foundation for such a high number of management units, and our data suggest three management units: the northern Baltic (demersal populations), southern Baltic with the Oresund straits and the most northwestern sampling sites (Skagerrak, Kattegat and North Sea).     

Adrenocortical stress responses influence an invasive vertebrate's fitness in an extreme environment

Continued range expansion into physiologically challenging environments requires invasive species to maintain adaptive phenotypic performance. The adrenocortical stress response, governed in part by glucocorticoid hormones, influences physiological and behavioural responses of vertebrates to environmental stressors. However, any adaptive role of this response in invasive populations that are expanding into extreme environments is currently unclear. We experimentally manipulated the adrenocortical stress response of invasive cane toads (Rhinella marina) to investigate its effect on phenotypic performance and fitness at the species'' range front in the Tanami Desert, Australia. Here, toads are vulnerable to overheating and dehydration during the annual hot–dry season and display elevated plasma corticosterone levels indicative of severe environmental stress. By comparing unmanipulated control toads with toads whose adrenocortical stress response was manipulated to increase acute physiological stress responsiveness, we found that control toads had significantly reduced daily evaporative water loss and higher survival relative to the experimental animals. The adrenocortical stress response hence appears essential in facilitating complex phenotypic performance and setting fitness trajectories of individuals from invasive species during range expansion.     

Climatic adaptation in an isolated and genetically impoverished amphibian population

The capacity of populations to respond adaptively to environmental change is essential for their persistence. Isolated populations often harbour reduced genetic variation, which is predicted to decrease adaptive potential, and can be detrimental under the current scenarios of global change. In this study, we examined climatic adaptation in larval life history traits in the pool frog Rana lessonae along a latitudinal gradient across the northern distribution area of the species, paying special attention to the isolated and genetically impoverished fringe populations in central Sweden. Larvae from eight populations within three geographic areas (Poland, Latvia and Sweden) were reared under three temperatures (19, 22 and 26°C) in a common garden laboratory experiment. We found clear evidence for latitudinal adaptation in R. lessonae populations, larvae from higher latitudes growing and developing faster than low‐latitude ones. Larvae from the Swedish populations were able to compensate for the effects of cooler temperatures and a shorter growth season with genetically higher growth and development rates (i.e. countergradient variation) in the two higher temperature treatments, but there was no difference among the populations at the lowest temperature treatment, which is likely to be close to the temperature limiting growth in R. lessonae. Our results demonstrate that isolated and genetically impoverished populations can be locally adapted, and identify the Swedish fringe populations as a significant conservation unit adapted to the northern environmental conditions.     

Larval Tolerance to Salinity in Three Species of Australian Anuran: An Indication of Saline Specialisation in Litoria aurea

Recent anthropogenic influences on freshwater habitats are forcing anuran populations to rapidly adapt to high magnitude changes in environmental conditions or face local extinction. We examined the effects of ecologically relevant elevated salinity levels on larval growth, metamorphosis and survival of three species of Australian anuran; the spotted marsh frog (Limnodynastes tasmaniensis), the painted burrowing frog (Neobatrachus sudelli) and the green and golden bell frog (Litoria aurea), in order to better understand the responses of these animals to environmental change. Elevated salinity (16% seawater) negatively impacted on the survival of L. tasmaniensis (35% survival) and N sudelli (0% survival), while reduced salinity had a negative impact on L. aurea. (16% seawater: 85% survival; 0.4% seawater: 35% survival). L. aurea tadpoles survived in salinities much higher than previously reported for this species, indicating the potential for inter-populations differences in salinity tolerance. In L. tasmaniensis and L. aurea, development to metamorphosis was fastest in low and high salinity treatments suggesting it is advantageous for tadpoles to invest energy in development in both highly favourable and developmentally challenging environments. We propose that this response might either maximise potential lifetime fecundity when tadpoles experience favourable environments, or, facilitate a more rapid escape from pond environments where there is a reduced probability of survival.     

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.