Latitudinal divergence of common frog (Rana temporaria) life history traits by natural selection: evidence from a comparison of molecular and quantitative genetic data

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Genomic data and common garden experiments reveal climate-driven selection on ecophysiological traits in two Mediterranean oaks

Improving our knowledge of how past climate-driven selection has acted on present-day trait population divergence is essential to understand local adaptation processes and improve our predictions of evolutionary trajectories in the face of altered selection pressures resulting from climate change. In this study, we investigated signals of selection on traits related to drought tolerance and growth rates in two Mediterranean oak species (Quercus faginea and Q. lusitanica) with contrasting distribution ranges and climatic niches. We genotyped 182 individuals from 24 natural populations of the two species using restriction-site-associated DNA sequencing and conducted a thorough functional characterization in 1602 seedlings from 21 populations cultivated in common garden experiments under contrasting watering treatments. Our genomic data revealed that both Q. faginea and Q. lusitanica have very weak population genetic structure, probably as a result of high rates of pollen-mediated gene flow among populations and large effective population sizes. In contrast, common garden experiments showed evidence of climate-driven divergent selection among populations on traits related to leaf morphology, physiology and growth in both species. Overall, our study suggests that climate is an important selective factor for Mediterranean oaks and that ecophysiological traits have evolved in drought-prone environments even in a context of very high rates of gene flow among populations.

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.

Conservation genetics of yellow-bellied toads (Bombina variegata): a matter of geographical scale and isolation

Amphibian populations world-wide are threatened by declines and extinctions mainly due to habitat loss and fragmentation. Habitat fragmentation threatens the yellow-bellied toad Bombina variegata in the northern and western regions of its distribution where it is strictly protected. We studied the genetic structure and diversity of populations at three geographical scales using microsatellite loci to detect potential threats for population persistence. At the local scale, we sampled four neighbouring localities at 1–2.6 km distance to detect effects of short-term (decades) fragmentation on connectivity. At the regional scale, five additional localities in the mountains of the Westerwald (Rhineland-Palatinate, Germany) were studied at up to 50.1 km distance to analyse genetic diversity and population structure. At the continental scale, we included data from regions in the northern distribution with fragmented populations (Hesse and Lower Saxony, Germany) and more continuous populations in the South (Alsace, France; Geneva, Switzerland; Trentino, Italy) to evaluate variation of genetic diversity. At the local scale, short-term fragmentation caused significant genetic differentiation between breeding assemblages only 1.4 km apart from each other. At the regional scale, we found notable genetic distance among localities. At the continental scale, we identified Alsace, Trentino and Geneva in the South as regions with low genetic structuring and high allelic richness, and the northern remaining regions in Germany as deeply structured with reduced allelic richness. We suggest that reduced genetic diversity and habitat fragmentation in northern regions makes these populations particularly vulnerable to decline. In conclusion, informed conservation management of B. variegata should focus on measures maintaining or improving connectivity among neighbouring populations.

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.

Small-scale phenotypic differentiation along complex stream gradients in a non-native amphipod

Background

Selective landscapes in rivers are made up by an array of selective forces that vary from source to downstream regions or between seasons, and local/temporal variation in fitness maxima can result in gradual spatio-temporal variation of phenotypic traits. This study aimed at establishing freshwater amphipods as future model organisms to study adaptive phenotypic diversification (evolutionary divergence and/or adaptive plasticity) along stream gradients.

Methods

We collected Gammarus roeselii from 16 sampling sites in the Rhine catchment during two consecutive seasons (summer and winter). Altogether, we dissected n = 1648 individuals and quantified key parameters related to morphological and life-history diversification, including naturally selected (e.g., gill surface areas) as well as primarily sexually selected traits (e.g., male antennae). Acknowledging the complexity of selective regimes in streams and the interrelated nature of selection factors, we assessed several abiotic (e.g., temperature, flow velocity) and biotic ecological parameters (e.g., conspecific densities, sex ratios) and condensed them into four principal components (PCs).

Results

Generalized least squares models revealed pronounced phenotypic differentiation in most of the traits investigated herein, and components of the stream gradient (PCs) explained parts of the observed differences. Depending on the trait under investigation, phenotypic differentiation could be ascribed to variation in abiotic conditions, anthropogenic disturbance (influx of thermally polluted water), or population parameters. For example, female fecundity showed altitudinal variation and decreased with increasing conspecific densities, while sexual dimorphism in the length of male antennae—used for mate finding and assessment—increased with increasing population densities and towards female-biased sex ratios.

Conclusions

We provide a comprehensive protocol for comparative analyses of intraspecific variation in life history traits in amphipods. Whether the observed phenotypic differentiation over small geographical distances reflects evolutionary divergence or plasticity (or both) remains to be investigated in future studies. Independent of the mechanisms involved, variation in several traits is likely to have consequences for ecosystem functions. For example, leaf-shredding in G. roeselii strongly depends on body size, which varied in dependence of several ecological parameters.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0327-8) contains supplementary material, which is available to authorized users.

Mapping of quantitative trait loci for life history traits segregating within common frog populations

The evolution of complex traits is often shaped by adaptive divergence. However, very little is known about the number, effect size, and location of the genomic regions influencing the variation of these traits in natural populations. Based on a dense linkage map of the common frog, Rana temporaria, we have localized, for the first time in amphibians, three significant and nine suggestive quantitative trait loci (QTLs) for metabolic rate, growth rate, development time, and weight at metamorphosis, explaining 5.6-18.9% of the overall phenotypic variation in each trait. We also found a potential pleiotropic QTL between development time and size at metamorphosis that, if confirmed, might underlie the previously reported genetic correlation between these traits. Furthermore, we demonstrate that the genetic variation linked to fitness-related larval traits segregates within Rana temporaria populations. This study provides the first insight into the genomic regions that affect larval life history traits in anurans, providing a valuable resource to delve further into the genomic basis of evolutionary change in amphibians.

Regional variation in drivers of connectivity for two frog species (Rana pretiosa and R. luteiventris) from the U.S. Pacific Northwest

Comparative landscape genetics has uncovered high levels of variability in which landscape factors affect connectivity among species and regions. However, the relative importance of species traits versus environmental variation for predicting landscape patterns of connectivity is unresolved. We provide evidence from a landscape genetics study of two sister taxa of frogs, the Oregon spotted frog (Rana pretiosa) and the Columbia spotted frog (Rana luteiventris) in Oregon and Idaho, USA. Rana pretiosa is relatively more dependent on moisture for dispersal than R. luteiventris, so if species traits influence connectivity, we predicted that connectivity among R. pretiosa populations would be more positively associated with moisture than R. luteiventris. However, if environmental differences are important drivers of gene flow, we predicted that connectivity would be more positively related to moisture in arid regions. We tested these predictions using eight microsatellite loci and gravity models in two R. pretiosa regions and four R. luteiventris regions (n = 1,168 frogs). In R. pretiosa, but not R. luteiventris, connectivity was positively related to mean annual precipitation, supporting our first prediction. In contrast, connectivity was not more positively related to moisture in more arid regions. Various temperature metrics were important predictors for both species and in all regions, but the directionality of their effects varied. Therefore, the pattern of variation in drivers of connectivity was consistent with predictions based on species traits rather than on environmental variation.

The phenotypic variability in Rana temporaria decreases in response to drying habitats

In this study, we evaluated the diversity of skin coloration as a proxy for phenotypic diversity. The European common frog (Rana temporaria) populations from the Southern slope of central Pyrenees lie at the limit of the species distribution in latitude and altitude. We analysed the relationship of skin color typology with different environmental variables and found a large decrease in skin type variety in frogs developing in temporary water bodies when compared to those developing in permanent water bodies. Our results show that our method can be used as a non-invasive way to study phenotypic diversity and suggest that adaptation to an early metamorphosis in a rapidly-drying habitat can have negative effects on adult phenotypic diversity. In light of these results, we argue that access to permanent water bodies is important to prevent loss of diversity in anuran populations and reduce their vulnerability to environmental impacts as well as pathogens.

Genetic variability and structure of an isolated population of Ambystoma altamirani, a mole salamander that lives in the mountains of one of the largest urban areas in the world

Amphibians are globally threatened by habitat loss and fragmentation; species within the order Ambystoma are not the exception, as there are 18 species of mole salamanders in México, of which 16 are endemic and all species are under some national or international status of protection. The mole salamander, Ambystoma altamirani is a microendemic species, which is distributed in central México, within the trans-Mexican volcanic belt, and is one of the most threatened species due to habitat destruction and the introduction of exotic species. Nine microsatellite markers were used to determine the genetic structure, genetic variability, effective population size, presence of bottlenecks and inbreeding coefficient of one population of A. altamirani to generate information which might help to protect and conserve this threatened species. We found two genetic subpopulations with significant level of genetic structure (F_ST = 0.005) and high levels of genetic variability (H_o = 0.883; H_e = 0.621); we also found a small population size (N_e = 8.8), the presence of historical (M = 0.486) and recent bottlenecks under IAM and TPM models, with a low, but significant coefficient of inbreeding (F_IS = -0.451). This information will help us to raise conservation strategies of this microendemic mole salamander species.

Genetic differentiation of regional populations of the widespread Asiatic toad (Bufo gargarizans), as revealed by development of novel microsatellite markers

Dispersal is a key component of a species’ life history, by influencing population persistence, genetic structure, adaptation and maintenance of genetic diversity. The Asiatic toad (Bufo gargarizans) is a widespread species in east Asia. However, we still have no knowledge of what kind of geographical scale equates to genetic differentiation within B. gargarizans. In this study, the population genetics of B. gargarizans was studied at five localities, with the Yangtze River running through the sampling area, in order to detect the level of genetic differentiation and the natural barriers to the species’ dispersal on a small geographic scale, by means of the development and use of novel microsatellite loci. These markers revealed a relatively high level of genetic diversity. Distinct genetic structure among populations in B. gargarizans was observed, as described by genetic distance, AMOVA, PCA and Geneland results. A weak but significant positive correlation between genetic distance and geographical distance. The combination of these findings suggests that the Yangtze River and geographic distance may act as effective barriers for B. gargarizans. These results serve as benchmark data for understanding the impacts of dispersal barriers and continued landscape research on B. gargarizans.

Shared and unique patterns of phenotypic diversification along a stream gradient in two congeneric species

Stream ecosystems show gradual variation of various selection factors, which can result in a zonation of species distributions and gradient evolution of morphological and life-history traits within species. Identifying the selective agents underlying such phenotypic evolution is challenging as different species could show shared and/or unique (species-specific) responses to components of the river gradient. We studied a stream gradient inhabited by two mosquitofishes (genus Gambusia) in the Río Grijalva basin in southern Mexico and found a patchy distribution pattern of both congeners along a stretch of 100 km, whereby one species was usually dominant at a given site. We uncovered both shared and unique patterns of diversification: some components of the stream gradient, including differences in piscine predation pressure, drove shared patterns of phenotypic divergence, especially in females. Other components of the gradient, particularly abiotic factors (max. annual temperature and temperature range) resulted in unique patterns of divergence, especially in males. Our study highlights the complexity of selective regimes in stream ecosystems. It exemplifies that even closely related, congeneric species can respond in unique ways to the same components of the river gradient and shows how both sexes can exhibit quite different patterns of divergence in multivariate phenotypic character suites.

Adaptive divergence in embryonic thermal plasticity among Atlantic salmon populations

In the context of global changes, the long-term viability of populations of endangered ectotherms may depend on their adaptive potential and ability to cope with temperature variations. We measured responses of Atlantic salmon embryos from four populations to temperature variations and used a QST -FST approach to study the adaptive divergence among these populations. Embryos were reared under two experimental conditions: a low temperature regime at 4 °C until eyed-stage and 10 °C until the end of embryonic development and a high temperature regime with a constant temperature of 10 °C throughout embryonic development. Significant variations among populations and population × temperature interactions were observed for embryo survival, incubation time and length. QST was higher than FST in all but one comparison suggesting an important effect of divergent selection. QST was also higher under the high-temperature treatment than at low temperature for length and survival due to a higher variance among populations under the stressful warmer treatment. Interestingly, heritability was lower for survival under high temperature in relation to a lower additive genetic variance under that treatment. Overall, these results reveal an adaptive divergence in thermal plasticity in embryonic life stages of Atlantic salmon suggesting that salmon populations may differentially respond to temperature variations induced by climate change. These results also suggest that changes in temperature may alter not only the adaptive potential of natural populations but also the selection regimes among them.

Role of selection and gene flow in population differentiation at the edge vs. interior of the species range differing in climatic conditions

Evaluating the relative importance of neutral and adaptive processes as determinants of population differentiation across environments is a central theme of evolutionary biology. We applied the QST-FST comparison flanked by a direct test for local adaptation to infer the role of climate-driven selection and gene flow in population differentiation of an annual grass Avena sterilis in two distinct parts of the species range, edge and interior, which represent two globally different climates, desert and Mediterranean. In a multiyear reciprocal transplant experiment, the plants of desert and Mediterranean origin demonstrated home advantage, and population differentiation in several phenotypic traits related to reproduction exceeded neutral predictions, as determined by comparisons of QST values with theoretical FST distributions. Thus, variation in these traits likely resulted from local adaptation to desert and Mediterranean environments. The two separate common garden experiments conducted with different experimental design revealed that two population comparisons, in contrast to multi-population comparisons, are likely to detect population differences in virtually every trait, but many of these differences reflect effects of local rather than regional environment. We detected a general reduction in neutral (SSR) genetic variation but not in adaptive quantitative trait variation in peripheral desert as compared with Mediterranean core populations. On the other hand, the molecular data indicated intensive gene flow from the Mediterranean core towards desert periphery. Although species range position in our study (edge vs. interior) was confounded with climate (desert vs. Mediterranean), the results suggest that the gene flow from the species core does not have negative consequences for either performance of the peripheral plants or their adaptive potential.

Postglacial recolonization in a cold climate specialist in western Europe: patterns of genetic diversity in the adder (Vipera berus) support the central-marginal hypothesis

Understanding the impact of postglacial recolonization on genetic diversity is essential in explaining current patterns of genetic variation. The central-marginal hypothesis (CMH) predicts a reduction in genetic diversity from the core of the distribution to peripheral populations, as well as reduced connectivity between peripheral populations. While the CMH has received considerable empirical support, its broad applicability is still debated and alternative hypotheses predict different spatial patterns of genetic diversity. Using microsatellite markers, we analysed the genetic diversity of the adder (Vipera berus) in western Europe to reconstruct postglacial recolonization. Approximate Bayesian Computation (ABC) analyses suggested a postglacial recolonization from two routes: a western route from the Atlantic Coast up to Belgium and a central route from the Massif Central to the Alps. This cold-adapted species likely used two isolated glacial refugia in southern France, in permafrost-free areas during the last glacial maximum. Adder populations further from putative glacial refugia had lower genetic diversity and reduced connectivity; therefore, our results support the predictions of the CMH. Our study also illustrates the utility of highly variable nuclear markers, such as microsatellites, and ABC to test competing recolonization hypotheses.

Commentary: When does understanding phenotypic evolution require identification of the underlying genes?

Adaptive evolution is fundamentally a genetic process. Over the past three decades, characterizing the genes underlying adaptive phenotypic change has revealed many important aspects of evolutionary change. At the same time, natural selection is often fundamentally an ecological process that can often be studied without identifying the genes underlying the variation on which it acts. This duality has given rise to disagreement about whether, and under what circumstances, it is necessary to identify specific genes associated with phenotypic change. This issue is of practical concern, especially for researchers who study nonmodel organisms, because of the often enormous cost and labor required to "go for the genes." We here consider a number of situations and questions commonly addressed by researchers. Our conclusion is that although gene identification can be crucial for answering some questions, there are others for which definitive answers can be obtained without finding underlying genes. It should thus not be assumed that considerations of "empirical completeness" dictate that gene identification is always desirable.

Molecular and quantitative genetic differentiation in Sitobion avenae populations from both sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between FST and QST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.

Behavioural adaptations of Rana temporaria to cold climates

Environmental conditions at the edge of a species' ecological optimum can exert great ecological or evolutionary pressure at local populations. For ectotherms like amphibians temperature is one of the most important abiotic factors of their environment as it influences directly their metabolism and sets limits to their distribution. Amphibians have evolved three ways to cope with sub-zero temperatures: freeze tolerance, freeze protection, freeze avoidance. The aim of this study was to assess which strategy common frogs at mid and high elevation use to survive and thrive in cold climates. In particular we (1) tested for the presence of physiological freeze protection, (2) evaluated autumnal activity and overwintering behaviour with respect to freeze avoidance and (3) assessed the importance of different high-elevation microhabitats for behavioural thermoregulation. Common frogs did not exhibit any signs of freeze protection when experiencing temperatures around 0 °C. Instead they retreated to open water for protection and overwintering. High elevation common frogs remained active for around the same period of time than their conspecifics at lower elevation. Our results suggest that at mid and high elevation common frogs use freeze avoidance alone to survive temperatures below 0 °C. The availability of warm microhabitats, such as rock or pasture, provides high elevation frogs with the opportunity of behavioural thermoregulation and thus allows them to remain active at temperatures at which common frogs at lower elevation cease activity.

Plasticity in metamorphic traits of Chinese brown frog (Rana chensinensis) tadpoles: the interactive effects of food level and rearing temperature

Larval amphibians are especially likely to encounter variation in temperature and resource availability because they live in a variety of aquatic habitats. In this study, plasticity in growth rates, survivorship, age and size at metamorphosis were examined in the Chinese brown frog (Rana chensinensis) under different combinations of rearing temperature and food level. Tadpoles reared at 23.2°C had larger mass at metamorphosis and a longer larval period than tadpoles reared at 25.3 and 28.1°C. High food level shortened the larval period and produced a larger size at metamorphosis. The interaction of rearing temperature and food level significantly influenced the growth and survival of Chinese brown frog tadpoles. At a low food level, tadpoles reared at 23.2°C displayed faster growth than tadpoles reared at the higher temperatures, and tadpoles had lower survival at 28.1°C than at the lower temperatures. Therefore, global warming or local manipulations of the environment could limit growth and development of Chinese brown frog tadpoles, resulting in low survival, but food availability may mediate effects of temperature. We suggest that Chinese brown frog tadpoles prefer cool a temperature and high food level for longer larval periods to capitalize on the opportunity of entering the terrestrial habitat with a larger body size.

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation--or clines--have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by-product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.

Environmental versus anthropogenic effects on population adaptive divergence in the freshwater snail Lymnaea stagnalis

Repeated pesticide contaminations of lentic freshwater systems located within agricultural landscapes may affect population evolution in non-target organisms, especially in species with a fully aquatic life cycle and low dispersal ability. The issue of evolutionary impact of pollutants is therefore conceptually important for ecotoxicologists. The impact of historical exposure to pesticides on genetic divergence was investigated in the freshwater gastropod Lymnaea stagnalis, using a set of 14 populations from contrasted environments in terms of pesticide and other anthropogenic pressures. The hypothesis of population adaptive divergence was tested on 11 life-history traits, using Q(ST)-F(ST) comparisons. Despite strong neutral differentiation (mean F(ST) = 0.291), five adult traits or parameters were found to be under divergent selection. Conversely, two early expressed traits showed a pattern consistent with uniform selection or trait canalization, and four adult traits appeared to evolve neutrally. Divergent selection patterns were mostly consistent with a habitat effect, opposing pond to ditch and channel populations. Comparatively, pesticide and other human pressures had little correspondence with evolutionary patterns, despite hatching rate impairment associated with global anthropogenic pressure. Globally, analyses revealed high genetic variation both at neutral markers and fitness-related traits in a species used as model in ecotoxicology, providing empirical support for the need to account for genetic and evolutionary components of population response in ecological risk assessment.

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.

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.

Molecular phylogeography and population genetic structure of an endangered species Pachyhynobius shangchengensis (hynobiid Salamander) in a fragmented habitat of southeastern China

The salamander Pachyhynobius shangchengensis (Hynobiidae) is a vulnerable species restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes in the Mount Dabieshan region in southeastern China. However, molecular phylogeography and population genetic structure of P. shangchengensis remain poorly investigated. In this study, we explored the genetic structure and phylogeography of P. shangchengensis based on partial sequences of the mitochondrial DNA (mtDNA) cytochrome b and cytochrome c oxidase subunit I genes. Fifty-one haplotypes and four clades were found among 93 samples. Phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographic regions separated by complicated topography and long distances. The distinct geographic distributions of all lineages and the estimated divergence time suggest spatial and temporal separation coinciding with climatic changes during the Pleistocene. Analysis of molecular variance indicated that most of the observed genetic variation occurred among the four groups, implying long-term interruption of gene flow, and the possible separation of P. shangchengensis into four management units for conservation.

Potential effects of climate change on the distribution of the common frog Rana temporaria at its northern range margin

Climate change is projected to be particularly strong in northern latitudes, and subarctic species are thus likely to be especially susceptible to the effects of climate warming. We forecast potential effects of climate change on the extent of the suitable habitat of the common frog,Rana temporaria,at the margin of its northern range. We investigated 179 potential breeding sites in subarctic Finland and subjected the data to detailed bioclimate envelope modelling using three state-of-the-art techniques: generalized additive models, maximum entropy and generalized boosting methods. Moreover, we included local environmental factors in the models to investigate whether they improve model performance. Under all tested climate change projections and irrespective of the modelling method, the suitable habitat forR. temporariaincreased in warming climate. The inclusion of local abiotic variables significantly improved the performance of the models. However, June temperature appeared to be the most informative variable in all modelling approaches: a major increase in the extent of suitable habitat occurred when it increased by 1°C. Overall, the modelling results indicate that the distribution of northernR. temporariais likely to be very sensitive to climate warming. The results also highlight the fact that overlooking local abiotic variation can significantly bias bioclimatic modelling results.

Geographically multifarious phenotypic divergence during speciation

Speciation is an important evolutionary process that occurs when barriers to gene flow evolve between previously panmictic populations. Although individual barriers to gene flow have been studied extensively, we know relatively little regarding the number of barriers that isolate species or whether these barriers are polymorphic within species. Herein, we use a series of field and lab experiments to quantify phenotypic divergence and identify possible barriers to gene flow between the butterfly species Lycaeides idas and Lycaeides melissa. We found evidence that L. idas and L. melissa have diverged along multiple phenotypic axes. Specifically, we identified major phenotypic differences in female oviposition preference and diapause initiation, and more moderate divergence in mate preference. Multiple phenotypic differences might operate as barriers to gene flow, as shown by correlations between genetic distance and phenotypic divergence and patterns of phenotypic variation in admixed Lycaeides populations. Although some of these traits differed primarily between species (e.g., diapause initiation), several traits also varied among conspecific populations (e.g., male mate preference and oviposition preference).

Phenotypic landscapes: phenological patterns in wild and cultivated barley

Phenotypic variation in natural populations is the outcome of the joint effects of environmentally induced adaptations and neutral processes on the genetic architecture of quantitative traits. In this study, we examined the role of adaptation in shaping wild barley phenotypic variation along different environmental gradients. Detailed phenotyping of 164 wild barley (Hordeum spontaneum) accessions from Israel (of the Barley1K collection) and 18 cultivated barley (H. vulgare) varieties was conducted in common garden field trials. Cluster analysis based on phenotypic data indicated that wild barley in this region can be differentiated into three ecotypes in accordance with their ecogeographical distribution: north, coast and desert. Population differentiation (Q(ST) ) for each trait was estimated using a hierarchical Bayesian model and compared to neutral differentiation (F(ST) ) based on 42 microsatellite markers. This analysis indicated that the three clusters diverged in morphological but not in reproductive characteristics. To address the issue of phenotypic variation along environmental gradients, climatic and soil gradients were compared with each of the measured traits given the geographical distance between sampling sites using a partial Mantel test. Flowering time and plant growth were found to be differentially correlated with climatic and soil characteristic gradients, respectively. The H. vulgare varieties were superior to the H. spontaneum accessions in yield components, yet resembled the Mediterranean types in vegetative characteristics and flowering time, which may indicate the geographical origin of domesticated barley.