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.     

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.     

Strong signature of selection in seeder populations but not in resprouters of the fynbos heath Erica coccinea (Ericaceae)

A higher frequency of natural selection is expected in populations of organisms with shorter generation times. In fire‐prone ecosystems, populations of seeder plants behave as functionally semelparous populations, with short generation times compared to populations of resprouter plants, which are truly iteroparous. Therefore, a stronger signature of natural selection should be detected in seeder populations, favoured by their shorter generation times and higher rates of population turnover. Here we test this idea in Erica coccinea from the Cape Floristic Region, which is dimorphic for post‐fire regeneration mode. We measured three floral traits supposedly subject to natural selection in seeder and resprouter populations. We then compared phenotypic trait variation with neutral genetic variation in each group of populations using P_ST–F_ST comparisons to detect signatures of natural selection in seeders and resprouters. We found a strong signature of selection in seeder populations, but not in resprouters. Furthermore, anthers of seeders were more exserted (and larger) than those of resprouters. These differences were maintained at sites where seeders and resprouters co‐occurred, suggesting that phenotypic plasticity or adaptation to different growth environments are unlikely explanations for trait variation. These results provide empirical support for the hypothesis that the genetic signature of natural selection is certainly more intense in seeder than in resprouter populations, favoured by their comparatively faster generation turnovers. Increased frequency of natural selection would increase differentiation among populations, thus promoting speciation in pyrophytic seeder lineages of the Cape flora.     

Distribution of variation over populations

Understanding the significance of the distribution of genetic or phenotypic variation over populations is one of the central concerns of population genetic and ecological research. The import of the research decisively depends on the measures that are applied to assess the amount of variation residing within and between populations. Common approaches can be classified under two perspectives: differentiation and apportionment. While the former focuses on differences (distances) in trait distribution between populations, the latter considers the division of the overall trait variation among populations. Particularly when multiple populations are studied, the apportionment perspective is usually given preference (via F _ST/G _ST indices), even though the other perspective is also relevant. The differences between the two perspectives as well as their joint conceptual basis can be exposed by referring them to the association between trait states and population affiliations. It is demonstrated that the two directions, association of population affiliation with trait state and of trait state with population affiliation, reflect the differentiation and the apportionment perspective, respectively. When combining both perspectives and applying the suggested measure of association, new and efficient methods of analysis result, as is outlined for population genetic processes. In conclusion, the association approach to an analysis of the distribution of trait variation over populations resolves problems that are frequently encountered with the apportionment perspective and its commonly applied measures in both population genetics and ecology, suggesting new and more comprehensive methods of analysis that include patterns of differentiation and apportionment.     

Evidence for adaptive divergence of thermal responses among Bemisia tabaci populations from tropical Colombia following a recent invasion

There is an increasing evidence that populations of ectotherms can diverge genetically in response to different climatic conditions, both within their native range and (in the case of invasive species) in their new range. Here, we test for such divergence in invasive whitefly Bemisia tabaci populations in tropical Colombia, by considering heritable variation within and between populations in survival and fecundity under temperature stress, and by comparing population differences with patterns established from putatively neutral microsatellite markers. We detected significant differences among populations linked to mean temperature (for survival) and temperature variation (for fecundity) in local environments. A Q_ST ? F_ST analysis indicated that phenotypic divergence was often larger than neutral expectations (Q_ST > F_ST). Particularly, for survival after a sublethal heat shock, this divergence remained linked to the local mean temperature after controlling for neutral divergence. These findings point to rapid adaptation in invasive whitefly likely to contribute to its success as a pest species. Ongoing evolutionary divergence also provides challenges in predicting the likely impact of Bemisia in invaded regions.     

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 (Q_ST) and neutral genetic divergence (F_ST 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 Q_STs–F_STs 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.     

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

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

Morphometric differentiation across House Sparrow Passer domesticus populations in Finland in comparison with the neutral expectation for divergence

To understand the biology of organisms it is important to take into account the evolutionary forces that have acted on their constituent populations. Neutral genetic variation is often assumed to reflect variation in quantitative traits under selection, though with even low neutral divergence there can be substantial differentiation in quantitative genetic variation associated with locally adapted phenotypes. To study the relative roles of natural selection and genetic drift in shaping phenotypic variation, the levels of quantitative divergence based on phenotypes (P_ST) and neutral genetic divergence (F_ST) can be compared. Such a comparison was made between 10 populations of Finnish House Sparrows (n = 238 individuals) collected in 2009 across the whole country. Phenotypic variation in tarsus‐length, wing‐length, bill‐depth, bill‐length and body mass were considered and 13 polymorphic microsatellite loci were analysed to quantify neutral genetic variation. Calculations of P_ST were based on Markov‐Chain Monte Carlo Bayesian estimates of phenotypic variances across and within populations. The robustness of the conclusions of the P_ST–F_ST comparison was evaluated by varying the proportion of variation due to additive genetic effects within and across populations. Our results suggest that body mass is under directional selection, whereas the divergence in other traits does not differ from neutral expectations. These findings suggest candidate traits for considering gene‐based studies of local adaptation. The recognition of locally adapted populations may be of value in the conservation of this declining species.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

Adaptive responses and invasion: the role of plasticity and evolution in snail shell morphology

Invasive species often exhibit either evolved or plastic adaptations in response to spatially varying environmental conditions. We investigated whether evolved or plastic adaptation was driving variation in shell morphology among invasive populations of the New Zealand mud snail (Potamopyrgus antipodarum) in the western United States. We found that invasive populations exhibit considerable shell shape variation and inhabit a variety of flow velocity habitats. We investigated the importance of evolution and plasticity by examining variation in shell morphological traits 1) between the parental and F₁ generations for each population and 2) among populations of the first lab generation (F₁) in a common garden, full‐sib design using Canonical Variate Analyses (CVA). We compared the F₁ generation to the parental lineages and found significant differences in overall shell shape indicating a plastic response. However, when examining differences among the F₁ populations, we found that they maintained among‐population shell shape differences, indicating a genetic response. The F₁ generation exhibited a smaller shell morph more suited to the low‐flow common garden environment within a single generation. Our results suggest that phenotypic plasticity in conjunction with evolution may be driving variation in shell morphology of this widespread invasive snail.     

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

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

Recent local adaptation of sockeye salmon to glacial spawning habitats

Salmonids spawn in highly diverse habitats, exhibit strong genetic population structuring, and can quickly colonize newly created habitats with few founders. Spawning traits often differ among populations, but it is largely unknown if these differences are adaptive or due to genetic drift. To test if sockeye salmon (Oncorhynchus nerka) populations are adapted to glacial, beach, and tributary spawning habitats, we examined variation in heritable phenotypic traits associated with spawning in 13 populations of wild sockeye salmon in Lake Clark, Alaska. These populations were commonly founded between 100 and 400 hundred sockeye salmon generations ago and exhibit low genetic divergence at 11 microsatellite loci (F _ST < 0.024) that is uncorrelated with spawning habitat type. We found that mean P _ST (phenotypic divergence among populations) exceeded neutral F _ST for most phenotypic traits measured, indicating that phenotypic differences among populations could not be explained by genetic drift alone. Phenotypic divergence among populations was associated with spawning habitat differences, but not with neutral genetic divergence. For example, female body color was lighter and egg color was darker in glacial than non-glacial habitats. This may be due to reduced sexual selection for red spawning color in glacial habitats and an apparent trade-off in carotenoid allocation to body and egg color in females. Phenotypic plasticity is an unlikely source of phenotypic differences because Lake Clark sockeye salmon spend nearly all their lives in a common environment. Our data suggest that Lake Clark sockeye salmon populations are adapted to spawning in glacial, beach and tributary habitats and provide the first evidence of a glacial spawning ecotype in salmonids. Glacial spawning habitats are often young (i.e., <200 years old) and ephemeral. Thus, local adaptation of sockeye salmon to glacial habitats appears to have occurred recently.     

GENE FLOW AND SELECTION ON PHENOTYPIC PLASTICITY IN AN ISLAND SYSTEM OF RANA TEMPORARIA

Gene flow is often considered to be one of the main factors that constrains local adaptation in a heterogeneous environment. However, gene flow may also lead to the evolution of phenotypic plasticity. We investigated the effect of gene flow on local adaptation and phenotypic plasticity in development time in island populations of the common frog Rana temporaria which breed in pools that differ in drying regimes. This was done by investigating associations between traits (measured in a common garden experiment) and selective factors (pool drying regimes and gene flow from other populations inhabiting different environments) by regression analyses and by comparing pairwise F_ST values (obtained from microsatellite analyses) with pairwise Q_ST values. We found that the degree of phenotypic plasticity was positively correlated with gene flow from other populations inhabiting different environments (among‐island environmental heterogeneity), as well as with local environmental heterogeneity within each population. Furthermore, local adaptation, manifested in the correlation between development time and the degree of pool drying on the islands, appears to have been caused by divergent selection pressures. The local adaptation in development time and phenotypic plasticity is quite remarkable, because the populations are young (less than 300 generations) and substantial gene flow is present among islands.     

Are QST–FST comparisons for natural populations meaningful?

Comparisons between putatively neutral genetic differentiation amongst populations, F_ST, and quantitative genetic variation, Q_ST, are increasingly being used to test for natural selection. However, we find that approximately half of the comparisons that use only data from wild populations confound phenotypic and genetic variation. We urge the use of a clear distinction between narrow‐sense Q_ST, which can be meaningfully compared with F_ST, and phenotypic divergence measured between populations, P_ST, which is inadequate for comparisons in the wild. We also point out that an unbiased estimate of Q_ST can be found using the so‐called ‘animal model’ of quantitative genetics.     

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 Q_ST–F_ST 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 Q_ST values with theoretical F_ST 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.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.     

Population size is weakly related to quantitative genetic variation and trait differentiation in a stream fish

How population size influences quantitative genetic variation and differentiation among natural, fragmented populations remains unresolved. Small, isolated populations might occupy poor quality habitats and lose genetic variation more rapidly due to genetic drift than large populations. Genetic drift might furthermore overcome selection as population size decreases. Collectively, this might result in directional changes in additive genetic variation (V_A) and trait differentiation (Q_ST) from small to large population size. Alternatively, small populations might exhibit larger variation in V_A and Q_ST if habitat fragmentation increases variability in habitat types. We explored these alternatives by investigating V_A and Q_ST using nine fragmented populations of brook trout varying 50‐fold in census size N (179–8416) and 10‐fold in effective number of breeders, N_b (18–135). Across 15 traits, no evidence was found for consistent differences in V_A and Q_ST with population size and almost no evidence for increased variability of V_A or Q_ST estimates at small population size. This suggests that (i) small populations of some species may retain adaptive potential according to commonly adopted quantitative genetic measures and (ii) populations of varying sizes experience a variety of environmental conditions in nature, however extremely large studies are likely required before any firm conclusions can be made.     

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.Subject terms: Biogeography, Ecological modelling, Evolutionary ecology, Evolution     

Selection on life‐history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation – the Monopolization hypothesis – has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F_ST) and in quantitative traits (Q_ST) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life‐history traits. Five life‐history traits were measured in four temperature × salinity combinations in common‐garden experiments. Population differentiation for neutral genetic markers – 11 microsatellite loci – was very high (F_ST = 0.482). Differentiation in life‐history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q_ST values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F_ST in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.     

GENETIC AND PHENOTYPIC DIVERGENCE BETWEEN LOW‐ AND HIGH‐ALTITUDE POPULATIONS OF TWO RECENTLY DIVERGED CINNAMON TEAL SUBSPECIES

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species’ range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low‐ and high‐altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P_ST= 0.56) and exhibited significant frequency differences in a single nonsynonymous α‐hemoglobin amino acid polymorphism (Asn/Ser‐α9; F_ST= 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F_ST= 0.004–0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high‐altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.