Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P_ST) which largely exceeds neutral genetic differentiation (F_ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.     

Contrasting patterns of morphological and neutral genetic divergence among geographically proximate populations of sockeye salmon Oncorhynchus nerka in Lake Aleknagik,Alaska

Levels of differentiation in morphological traits (age at maturity, body length at age, egg mass and body depth) and spawning time were examined in sockeye salmon Oncorhynchus nerka from three geographically proximate but physically distinct creeks in Lake Aleknagik, Alaska. Happy Creek fish had significantly greater values for most measured morphological traits, and Eagle Creek fish spawned significantly later than fish in the other creeks. Phenotypic differentiation between creeks, measured as P_ST, was then compared with microsatellite marker differentiation between creeks, measured as F_ST. No correlations were apparent between P_ST and F_ST values, and P_ST values were generally significantly larger than zero (P_ST= 0·0018–0·31) whereas F_ST values were not (F_ST=?0·0004 to 0·0016). The insignificant pair‐wise F_ST values between creek samples indicated that gene flow occurs between creeks, assuming the creek populations have reached migration–drift equilibrium. However, the strong homing behaviour of sockeye salmon precludes a scenario in which fish from the three creeks constitute a single population that segregates by body size. Rather, significant phenotypic differentiation suggests that strong divergent selection occurs on the phenotypic traits despite the homogenizing effects of gene flow.     

Risk‐taking and survival in the House Sparrow Passer domesticus: are plumage ornaments costly?

An important cost of sexual and social colour signals may be increased conspicuousness of the animals to visual predators. Although such predation costs have repeatedly been proposed for various ornaments of birds, especially for melanised and depigmented signals with low presumed production costs, tests of this hypothesis are rare and inconclusive. In this study we investigated whether individual variation in plumage ornamentation was associated with predator-related risk-taking behaviour and short-term survival in free-living House Sparrows Passer domesticus . In a large sample of birds we measured three aspects of coloration used in sexual and social signalling: the size of the melanised black throat patch in males, and the area and conspicuousness of the depigmented wing-bar in both sexes. We measured risk-taking by manipulating the distance of feeders from shelters, and recorded individually ringed birds feeding close to and far from shelter. Sparrows seemed to perceive the farther feeders as more risky, as indicated by the shorter time spent and smaller groups feeding far from rather than close to shelter. However, the use of the more risky (farther) feeders was not related to any of the colour traits we measured, suggesting that Sparrows do not adjust their risk-taking behaviour to their ornamentation. Males (the more ornamented sex) did not take less risk than females. Furthermore, we found no evidence that larger throat patches or more ornamented wing-bars reduced the probability of survival. Our findings were robust and consistent across multiple approaches, even when we controlled for several potential confounding effects. These results do not support the suggestion that melanised and depigmented plumage ornaments have significant predation costs in House Sparrows.     

Providing chicks with extra food lowers male but not female provisioning in the House Sparrow Passer domesticus

We assessed whether adult House Sparrows Passer domesticus adjusted their provisioning in response to an experimental increase in the nutritional condition of their nestlings. When we supplemented chicks directly with additional food, male parents, but not female parents, reduced their provisioning. The results for males, but not females, run contrary to a previous experiment in this species. In addition, female provisioning was positively associated with both brood size and the age of the brood. In contrast, whereas male provisioning was positively associated with brood size, males did not increase provisioning as their chicks grew older. Males, but not females, exhibited repeatability in their provisioning. Food supplementation had a larger positive effect upon nestling survival in smaller broods than in larger broods. Overall, there appear to be fundamental differences between males and females in how decisions regarding the level of parental investment in the current brood are made.     

Signals of local adaptation across an environmental gradient among highly connected populations of the Dead Sea Sparrow Passer moabiticus in Israel

Populations found at the edge of a species range often have decreased genetic diversity, which together with high gene flow may reduce the ability of a species to adapt to local environmental conditions. The Dead Sea Sparrow Passer moabiticus occupies a disjointed range, where the Israeli populations are considered peripheral and fragmented. The species is also thought to have undergone a recent range expansion. We aimed to describe the genetic and morphological variation of the Israeli populations and to determine the extent of gene flow among them. We expected that because of the small latitudinal gradient across Israel and the recent range expansion of the species that Dead Sea Sparrow populations would show no significant morphological adaptation to local environmental conditions, and that considerable gene flow would be taking place among populations. Our findings indicate the existence of gene flow, suggesting high connectivity among populations, but recovered no support for a recent range expansion, possibly due to insufficient time since expansion for mutations to have accumulated. However, despite recurrent gene flow among populations, latitudinal variation in wing length (male and female) and body mass (male) was indicative of local adaptation across Israel, in accordance with Bergmann's rule.     

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.     

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 Q_ST–F_ST 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. Q_ST was higher than F_ST in all but one comparison suggesting an important effect of divergent selection. Q_ST 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.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Natural selection and neutral evolutionary processes contribute to genetic divergence in leaf traits across a precipitation gradient in the tropical oak Quercus oleoides

The impacts of drought are expanding worldwide as a consequence of climate change. However, there is still little knowledge of how species respond to long‐term selection in seasonally dry ecosystems. In this study, we used Q_ST‐F_ST comparisons to investigate (i) the role of natural selection on population genetic differentiation for a set of functional traits related to drought resistance in the seasonally dry tropical oak Quercus oleoides and (ii) the influence of water availability at the site of population origin and in experimental treatments on patterns of trait divergence. We conducted a thorough phenotypic characterization of 1912 seedlings from ten populations growing in field and greenhouse common gardens under replicated watering treatments. We also genotyped 218 individuals from the same set of populations using eleven nuclear microsatellites. Q_ST distributions for leaf lamina area, specific leaf area, leaf thickness and stomatal pore index were higher than F_ST distribution. Results were consistent across growth environments. Genetic differentiation among populations for these functional traits was associated with the index of moisture at the origin of the populations. Together, our results suggest that drought is an important selective agent for Q. oleoides and that differences in length and severity of the dry season have driven the evolution of genetic differences in functional traits.     

Parallel invasions produce heterogenous patterns of life history adaptation: rapid divergence in an invasive insect

Biotic invasions provide a natural experiment in evolution: when invasive species colonize new ranges, they may evolve new clines in traits in response to environmental gradients. Yet it is not clear how rapidly such patterns can evolve and whether they are consistent between regions. We compare four populations of the invasive cabbage white butterfly (Pieris rapae) from North America and Japan, independently colonized by P. rapae 150 years ago and 300 years ago, respectively. On each continent, we employed a northern and southern population to compare the effects of latitude on body mass, development rate and immune function. For each population, we used a split‐sibling family design in which siblings were reared at either warm (26.7 °C) or cool (20 °C) temperatures to determine reaction norms for each trait. Latitudinal patterns in development time were similar between the two continents. In contrast, there were strong geographical differences in reaction norms for body size, but no consistent effects of latitude; there were no detectable effects of latitude or continent on immune function. These results imply that some life history traits respond consistently to selection along climatic gradients, whereas other traits may respond to local environmental factors, or not at all.     

A method for adjusting allele frequencies in the case of microsatellite allele drop‐out

A new method is proposed to adjust allele frequencies when allelic drop‐out is common. This method assumes Hardy–Weinberg equilibrium (HWE), and treats the problematic alleles as a one‐locus two‐allele system with dominance. By assuming that the homozygote frequency of the ‘recessive’ allele is measured correctly, we can back calculate the allele frequency of the ‘dominant’ allele, and adjust the heterozygote frequency accordingly. The drawback is that multilocus genotypes cannot be constructed and tests that use deviations from Hardy–Weinberg such as tests for bottlenecks become impossible. An example is given where a large homozygote excess (F_IS = 0.44) is adjusted to a reasonable level (F_IS = 0.046). The effect of scoring error was set in relation to sampling error and while F_IS values can be seriously biased, F_ST values are not necessarily so, if scoring error and sample size are both low. As sample size increases, the effect of scoring error increases.     

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.     

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of F_ST and to evaluate whether highly variable microsatellite loci resulted in deflated F_ST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (F_CT = 0.010) and within (F_ST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (P_CT > F_CT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (P_ST > F_ST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated F_ST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for P_ST–F_ST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.     

Phenotypic divergence during the invasion of Phyla canescens in Australia and France: evidence for selection-driven evolution

Rapid adaptive evolution has been advocated as a mechanism that promotes invasion. Demonstrating adaptive evolution in invasive species requires rigorous analysis of phenotypic shifts driven by selection. Here, we document selection-driven evolution of Phyla canescens , an Argentine weed, in two invaded regions (Australia and France). Invasive populations possessed similar or higher diversity than native populations, and displayed mixed lineages from different sources, suggesting that genetic bottlenecks in both countries might have been alleviated by multiple introductions. Compared to native populations, Australian populations displayed more investment in sexual reproduction, whereas French populations possessed enhanced vegetative reproduction and growth. We partitioned evolutionary forces (selection vs. stochastic events) using two independent methods. Results of both analyses suggest that the pattern of molecular and phenotypic variability among regions was consistent with selection-driven evolution, rather than stochastic events. Our findings indicate that selection has shaped the evolution of P . canescens in two different invaded regions.
Ecology Letters (2010) 13: 32–44     

Reaction norm variants for male calling song in populations of Achroia grisella (Lepidoptera: Pyralidae): toward a resolution of the lek paradox

Significant additive genetic variance often occurs for male advertisement traits in spite of the directional selection imposed by female choice, a problem generally known in evolutionary biology as the lek paradox. One hypothesis, which has limited support from recent studies, for the resolution of this paradox is the role of genotype x environment interaction in which no one genotype exhibits the superior performance in all environments--a crossover of reaction norms. However, these studies have not characterized the actual variation of reaction norms present in natural populations, and the extent to which crossover maintains genetic variance remains unknown. Here, we present a study of genotype x environment interaction for the male calling song in populations of Achroia grisella (Lepidoptera: Pyralidae; lesser waxmoth). We report significant variance among reaction norms for male calling song in two North American populations of A. grisella as measured along temperature, food availability, and density gradients, and there is a relatively high incidence of crossover of the temperature reaction norms. This range of reaction norm variants and their crossover may reflect the co-occurrence of plastic and canalized genotypes, and we argue that the different responses of these variants along environmental gradients may contribute toward the maintenance of genetic variance for male song.     

The search for Pleiades in trait constellations: functional integration and phenotypic selection in the complex flowers of Morrenia brachystephana (Apocynaceae)

Pollinator‐mediated natural selection on single traits, such as corolla tube or spur length, has been well documented. However, flower phenotypes are usually complex, and selection is expected to act on several traits that functionally interact rather than on a single isolated trait. Despite the fact that selection on complex phenotypes is expectedly widespread, multivariate selection modelling on such phenotypes still remains under‐explored in plants. Species of the subfamily Asclepiadoideae (Apocynaceae) provide an opportunity to study such complex flower contrivances integrated by fine‐scaled organs from disparate developmental origin. We studied the correlation structure among linear floral traits (i) by testing a priori morphological, functional or developmental hypotheses among traits and (ii) by exploring the organization of flower covariation, considering alternative expectations of modular organization or whole flower integration through conditional dependence analysis (CDA) and integration matrices. The phenotypic selection approach was applied to determine whether floral traits involved in the functioning of the pollination mechanism were affected by natural selection. Floral integration was low, suggesting that flowers are organized in more than just one correlation pleiad; our hypothetical functional correlation matrix was significantly correlated with the empirical matrix, and the CDA revealed three putative modules. Analyses of phenotypic selection showed significant linear and correlational gradients, lending support to expectations of functional interactions between floral traits. Significant correlational selection gradients found involved traits of different floral whorls, providing evidence for the existence of functional integration across developmental domains.     

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short‐range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19–29 individuals from each population. The sample includes two pairs of low‐ and high‐altitude populations from two different valleys. High‐altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low‐altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low‐ and high‐altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low‐ and high‐altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low‐ and high‐altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.     

Congruent population genetic structures and divergence histories in anther‐smut fungi and their host plants Silene italica and the Silene nutans species complex

Study of the congruence of population genetic structure between hosts and pathogens gives important insights into their shared phylogeographical and coevolutionary histories. We studied the population genetic structure of castrating anther‐smut fungi (genus Microbotryum) and of their host plants, the Silene nutans species complex, and the morphologically and genetically closely related Silene italica, which can be found in sympatry. Phylogeographical population genetic structure related to persistence in separate glacial refugia has been recently revealed in the S. nutans plant species complex across Western Europe, identifying several distinct lineages. We genotyped 171 associated plant–pathogen pairs of anther‐smut fungi and their host plant individuals using microsatellite markers and plant chloroplastic single nucleotide polymorphisms. We found clear differentiation between fungal populations parasitizing S. nutans and S. italica plants. The population genetic structure of fungal strains parasitizing the S. nutans plant species complex mirrored the host plant genetic structure, suggesting that the pathogen was isolated in glacial refugia together with its host and/or that it has specialized on the plant genetic lineages. Using random forest approximate Bayesian computation (ABC‐RF), we found that the divergence history of the fungal lineages on S. nutans was congruent with that previously inferred for the host plant and probably occurred with ancient but no recent gene flow. Genome sequences confirmed the genetic structure and the absence of recent gene flow between fungal genetic lineages. Our analyses of individual host–pathogen pairs contribute to a better understanding of co‐evolutionary histories between hosts and pathogens in natural ecosystems, in which such studies remain scarce.     

Disentangling the divergence and cladogenesis in the freshwater crab species (Potamonautidae: Potamonautes perlatus sensu lato) in the Cape Fold Mountains,South Africa,with the description of two novel cryptic lineages

River networks of major drainages can form barriers that shape the phylogeography of freshwater organisms, particularly those with low dispersal capabilities. Freshwater crab species' distributions can be used to examine hydrological patterns to expose historical drainage interconnectivity. We used molecular sequence data (mitochondrial and nuclear DNA) and divergence time estimations to determine the phylogeography of the freshwater crab, Potamonautes perlatus sensu lato, from six drainage systems along the Cape Fold Mountains, South Africa. Two major clades were detected: clade 1 comprised two geographically discrete haploclades occurring in southern flowing drainages, whereas clade 2 included specimens from western flowing drainages. Divergence time estimations suggested a Pleistocene (c. 2.61 Mya) divergence of P. perlatus s.l. The Pleistocene was associated with arid conditions and drainage contractions. However, it is likely that during the mesic conditions of the Pleistocene, P. perlatus s.l. migrated and diverged into contemporary patterns. We conclude that three lineages are nested within P. perlatus s.l., two representing novel species. Potamonautes perlatus sensu stricto is confined to western flowing drainages. The two novel species both occurring in southern flowing drainages are described here: P otamonautes barbarai sp. nov. occurs in the Gamtoos and Gourits Rivers and P otamonautes barnardi sp. nov. in the Breede River. © 2014 The Linnean Society of London