The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Ecological causes of morphological evolution in the three‐spined stickleback

The central assumption of evolutionary theory is that natural selection drives the adaptation of populations to local environmental conditions, resulting in the evolution of adaptive phenotypes. The three‐spined stickleback (Gasterosteus aculeatus) displays remarkable phenotypic variation, offering an unusually tractable model for understanding the ecological mechanisms underpinning adaptive evolutionary change. Using populations on North Uist, Scotland we investigated the role of predation pressure and calcium limitation on the adaptive evolution of stickleback morphology and behavior. Dissolved calcium was a significant predictor of plate and spine morph, while predator abundance was not. Stickleback latency to emerge from a refuge varied with morph, with populations with highly reduced plates and spines and high predation risk less bold. Our findings support strong directional selection in three‐spined stickleback evolution, driven by multiple selective agents.     

Population‐genomic approach reveals adaptive floral divergence in discrete populations of a hawk moth‐pollinated violet

Local adaptation to contrasting biotic or abiotic environments is an important evolutionary step that presumably precedes floral diversification at the species level, yet few studies have demonstrated the adaptive nature of intraspecific floral divergence in wild plant populations. We combine a population‐genomic approach with phenotypic information on floral traits to examine whether the differentiation in metric floral traits exhibited by 14 populations of the southern Spanish hawk moth‐pollinated violet Viola cazorlensis reflects adaptive divergence. Screening of many amplified fragment length polymorphism (AFLP) loci using a multiple‐marker‐based neutrality test identified nine outlier loci (2.6% of the total) that departed from neutral expectations and were potentially under selection. Generalized analysis of molecular variance revealed significant relationships between genetic distance and population divergence in three floral traits when genetic distance was based on outlier loci, but not when it was based on neutral ones. Population means of floral traits were closely correlated with population scores on the first principal coordinate axis of the genetic distance matrix using outlier loci, and with the allelic frequencies of four of the outlier loci. Results strongly support the adaptive nature of intraspecific floral divergence exhibited by V. cazorlensis and illustrate the potential of genome scans to identify instances of adaptive divergence when used in combination with phenotypic information.     

Evidence for a multi‐peak adaptive landscape in the evolution of trophic morphology in terapontid fishes

Using the Australian marine‐freshwater terapontid fishes as a model system, we examined the role of dietary phenotypic optima in an adaptive macro‐evolutionary landscape. Comparative modelling relying on both a priori and data‐driven identification of selective regimes suggested multi‐peak models as best describing much of the dietary phenotypic landscape of terapontids. Both approaches identified common phenotypic optima for different lineages of marine and freshwater herbivores, and minimal differentiation between carnivores and omnivores, irrespective of their phylogenetic relationships, as the model best describing morphological evolution. Significant correlations also existed between these phenotypic axes and proportions of non‐animal dietary items in species’ diets. While simulation results provided evidence for a multi‐peak adaptive landscape in the evolution of trophic morphology in terapontids, they could not rule out chance convergence in these adaptive peaks. However, they do provide scope for identifying areas for more detailed, functionally specific study of phenotypic convergence in herbivorous terapontid trophic habits. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 623–634.     

Symbiont‐mediated phenotypic variation without co‐evolution in an insect–fungus association

Recent studies have shown that symbionts can be a source of adaptive phenotypic variation for their hosts. It is assumed that co‐evolution between hosts and symbionts underlies these ecologically significant phenotypic traits. We tested this assumption in the ectosymbiotic fungal associate of the gall midge Asteromyia carbonifera. Phylogenetic analysis placed the fungal symbiont within a monophyletic clade formed by Botryosphaeria dothidea, a typically free‐living (i.e. not associated with an insect host) plant pathogen. Symbiont isolates from four divergent midge lineages demonstrated none of the patterns common to heritable microbial symbioses, including parallel diversification with their hosts, substitution rate acceleration, or A+T nucleotide bias. Amplified fragment length polymorphism genotyping of the symbiont revealed that within‐lineage genetic diversity was not clustered along host population lines. Culture‐based experiments demonstrated that the symbiont‐mediated variation in gall phenotype is not borne out in the absence of the midge. This study shows that symbionts can be important players in phenotypic variation for their hosts, even in the absence of a co‐evolutionary association.     

Genetic variance for behavioural ‘predictability’ of stress response

Genetic factors underpinning phenotypic variation are required if natural selection is to result in adaptive evolution. However, evolutionary and behavioural ecologists typically focus on variation among individuals in their average trait values and seek to characterize genetic contributions to this. As a result, less attention has been paid to if and how genes could contribute towards within‐individual variance or trait ‘predictability’. In fact, phenotypic ‘predictability’ can vary among individuals, and emerging evidence from livestock genetics suggests this can be due to genetic factors. Here, we test this empirically using repeated measures of a behavioural stress response trait in a pedigreed population of wild‐type guppies. We ask (a) whether individuals differ in behavioural predictability and (b) whether this variation is heritable and so evolvable under selection. Using statistical methodology from the field of quantitative genetics, we find support for both hypotheses and also show evidence of a genetic correlation structure between the behavioural trait mean and individual predictability. We show that investigating sources of variability in trait predictability is statistically tractable and can yield useful biological interpretation. We conclude that, if widespread, genetic variance for ‘predictability’ will have major implications for the evolutionary causes and consequences of phenotypic variation.     

Within‐population Y‐linked genetic variation for lifespan in Drosophila melanogaster

The view that the Y chromosome is of little importance for phenotypic evolution stems from early studies of Drosophila melanogaster. This species’ Y chromosome contains only 13 protein‐coding genes, is almost entirely heterochromatic and is not necessary for male viability. Population genetic theory further suggests that non‐neutral variation can only be maintained at the Y chromosome under special circumstances. Yet, recent studies suggest that the D. melanogaster Y chromosome trans‐regulates hundreds to thousands of X and autosomal genes. This finding suggests that the Y chromosome may play a far more active role in adaptive evolution than has previously been assumed. To evaluate the potential for the Y chromosome to contribute to phenotypic evolution from standing genetic variation, we test for Y‐linked variation in lifespan within a population of D. melanogaster. Assessing variation for lifespan provides a powerful test because lifespan (i) shows sexual dimorphism, which the Y is primarily predicted to contribute to, (ii) is influenced by many genes, which provides the Y with many potential regulatory targets and (iii) is sensitive to heterochromatin remodelling, a mechanism through which the Y chromosome is believed to regulate gene expression. Our results show a small but significant effect of the Y chromosome and thus suggest that the Y chromosome has the potential to respond to selection from standing genetic variation. Despite its small effect size, Y‐linked variation may still be important, in particular when evolution of sexual dimorphism is genetically constrained elsewhere in the genome.     

Phenotypic variability can promote the evolution of adaptive plasticity by reducing the stringency of natural selection

Adaptive phenotypic plasticity is a potent but not ubiquitous solution to environmental heterogeneity, driving interest in what factors promote and limit its evolution. Here, a novel computational model representing stochastic information flow in development is used to explore evolution from a constitutive phenotype to an adaptively plastic response. Results show that populations tend to evolve robustness to developmental stochasticity, but that this evolved robustness limits evolvability; specifically, robust genotypes have less ability to evolve adaptive plasticity when presented with a mix of both the ancestral environment and a new environment. Analytic calculations and computational experiments confirm that this constraint occurs when the initial mutational steps towards plasticity are pleiotropic, such that mutant fitnesses decline in the environment to which their parents are well‐adapted. Greater phenotypic variability improves evolvability in the model by lessening this decline as well as by improving the fitness of partial adaptations to the new environment. By making initial plastic mutations more palatable to natural selection, phenotypic variability can increase the evolvability of an innovative, plastic response without improving evolvability to simpler challenges such as a shifted optimum in a single environment. Populations that evolved robustness by negative feedback between the trait and its rate of change show a particularly strong constraining effect on the evolvability of plasticity, revealing another mechanism by which evolutionary history can limit later innovation. These results document a novel mechanism by which weakening selection could actually stimulate the evolution of a major innovation.     

SPECIES‐SPECIFIC DIFFERENCES IN ADAPTIVE PHENOTYPIC PLASTICITY IN AN ECOLOGICALLY RELEVANT TROPHIC TRAIT: HYPERTROPHIC LIPS IN MIDAS CICHLID FISHES

The spectacular species richness of cichlids and their diversity in morphology, coloration, and behavior have made them an ideal model for the study of speciation and adaptive evolution. Hypertrophic lips evolved repeatedly and independently in African and Neotropical cichlid radiations. Cichlids with hypertrophic lips forage predominantly in rocky crevices and it has been hypothesized that mechanical stress caused by friction could result in larger lips through phenotypic plasticity. To test the influence of the environment on the size and development of lips, we conducted a series of breeding and feeding experiments on Midas cichlids. Full‐sibs of Amphilophus labiatus (thick‐lipped) and Amphilophus citrinellus (thin‐lipped) each were split into a control group which was fed food from the water column and a treatment group whose food was fixed to substrates. We found strong evidence for phenotypic plasticity on lip area in the thick‐lipped species, but not in the thin‐lipped species. Intermediate phenotypic values were observed in hybrids from thick‐ and thin‐lipped species reared under “control” conditions. Thus, both a genetic, but also a phenotypic plastic component is involved in the development of hypertrophic lips in Neotropical cichlids. Moreover, species‐specific adaptive phenotypic plasticity was found, suggesting that plasticity is selected for in recent thick‐lipped species.     

Genetic evolution,plasticity, and bet‐hedging as adaptive responses to temporally autocorrelated fluctuating selection: A quantitative genetic model

Adaptive responses to autocorrelated environmental fluctuations through evolution in mean reaction norm elevation and slope and an independent component of the phenotypic variance are analyzed using a quantitative genetic model. Analytic approximations expressing the mutual dependencies between all three response modes are derived and solved for the joint evolutionary outcome. Both genetic evolution in reaction norm elevation and plasticity are favored by slow temporal fluctuations, with plasticity, in the absence of microenvironmental variability, being the dominant evolutionary outcome for reasonable parameter values. For fast fluctuations, tracking of the optimal phenotype through genetic evolution and plasticity is limited. If residual fluctuations in the optimal phenotype are large and stabilizing selection is strong, selection then acts to increase the phenotypic variance (bet‐hedging adaptive). Otherwise, canalizing selection occurs. If the phenotypic variance increases with plasticity through the effect of microenvironmental variability, this shifts the joint evolutionary balance away from plasticity in favor of genetic evolution. If microenvironmental deviations experienced by each individual at the time of development and selection are correlated, however, more plasticity evolves. The adaptive significance of evolutionary fluctuations in plasticity and the phenotypic variance, transient evolution, and the validity of the analytic approximations are investigated using simulations.     

Adaptive processes drive ecomorphological convergent evolution in antwrens (Thamnophilidae)

Phylogenetic niche conservatism (PNC) and convergence are contrasting evolutionary patterns that describe phenotypic similarity across independent lineages. Assessing whether and how adaptive processes give origin to these patterns represent a fundamental step toward understanding phenotypic evolution. Phylogenetic model‐based approaches offer the opportunity not only to distinguish between PNC and convergence, but also to determine the extent that adaptive processes explain phenotypic similarity. The Myrmotherula complex in the Neotropical family Thamnophilidae is a polyphyletic group of sexually dimorphic small insectivorous forest birds that are relatively homogeneous in size and shape. Here, we integrate a comprehensive species‐level molecular phylogeny of the Myrmotherula complex with morphometric and ecological data within a comparative framework to test whether phenotypic similarity is described by a pattern of PNC or convergence, and to identify evolutionary mechanisms underlying body size and shape evolution. We show that antwrens in the Myrmotherula complex represent distantly related clades that exhibit adaptive convergent evolution in body size and divergent evolution in body shape. Phenotypic similarity in the group is primarily driven by their tendency to converge toward smaller body sizes. Differences in body size and shape across lineages are associated to ecological and behavioral factors.     

Parallel and non‐parallel morphological divergence among foraging specialists in European whitefish (Coregonus lavaretus)

Parallel phenotypic evolution occurs when independent populations evolve similar traits in response to similar selective regimes. However, populations inhabiting similar environments also frequently show some phenotypic differences that result from non‐parallel evolution. In this study, we quantified the relative importance of parallel evolution to similar foraging regimes and non‐parallel lake‐specific effects on morphological variation in European whitefish (Coregonus lavaretus). We found evidence for both lake‐specific morphological characteristics and parallel morphological divergence between whitefish specializing in feeding on profundal and littoral resources in three separate lakes. Foraging specialists expressed similar phenotypes in different lakes in both overall body shape and selected measured morphological traits. The morphology of the two whitefish specialists resembled that predicted from other fish species, supporting the conclusion of an adaptive significance of the observed morphological characteristics. Our results indicate that divergent natural selection resulting from foraging specialization is driving and/or maintaining the observed parallel morphological divergence. Whitefish in this study may represent an early stage of divergence towards the evolution of specialized morphs.     

Opportunities and challenges of next‐generation sequencing applications in ecological epigenetics

Evolutionary theory posits that adaptation can result when populations harbour heritable phenotypic variation for traits that increase tolerance to local conditions. However, the actual mechanisms that underlie heritable phenotypic variation are not completely understood (Keller 2014 ). Recently, the potential role of epigenetic mechanisms in the process of adaptive evolution has been the subject of much debate (Pigliucci & Finkelman 2014 ). Studies of variation in DNA methylation in particular have shown that natural populations harbour high amounts of epigenetic variation, which can be inherited across generations and can cause heritable trait variation independently of genetic variation (Kilvitis et al. 2014 ). While we have made some progress addressing the importance of epigenetics in ecology and evolution using methylation‐sensitive AFLP (MS‐AFLP), this approach provides relatively few anonymous and dominant markers per individual. MS‐AFLP are difficult to link to functional genomic elements or phenotype and are difficult to compare directly to genetic variation, which has limited the insights drawn from studies of epigenetic variation in natural nonmodel populations (Schrey et al. 2013 ). In this issue, Platt et al. provide an example of a promising approach to address this problem by applying a reduced representation bisulphite sequencing (RRBS) approach based on next‐generation sequencing methods in an ecological context.     

Genetic component of flammability variation in a Mediterranean shrub

Recurrent fires impose a strong selection pressure in many ecosystems worldwide. In such ecosystems, plant flammability is of paramount importance because it enhances population persistence, particularly in non‐resprouting species. Indeed, there is evidence of phenotypic divergence of flammability under different fire regimes. Our general hypothesis is that flammability‐enhancing traits are adaptive; here, we test whether they have a genetic component. To test this hypothesis, we used the postfire obligate seeder Ulex parviflorus from sites historically exposed to different fire recurrence. We associated molecular variation in potentially adaptive loci detected with a genomic scan (using AFLP markers) with individual phenotypic variability in flammability across fire regimes. We found that at least 42% of the phenotypic variation in flammability was explained by the genetic divergence in a subset of AFLP loci. In spite of generalized gene flow, the genetic variability was structured by differences in fire recurrence. Our results provide the first field evidence supporting that traits enhancing plant flammability have a genetic component and thus can be responding to natural selection driven by fire. These results highlight the importance of flammability as an adaptive trait in fire‐prone ecosystems.     

Maternal source of variability in the embryo development of an annual killifish

Organisms inhabiting unpredictable environments often evolve diversified reproductive bet‐hedging strategies, expressed as production of multiple offspring phenotypes, thereby avoiding complete reproductive failure. To cope with unpredictable rainfall, African annual killifish from temporary savannah pools lay drought‐resistant eggs that vary widely in the duration of embryo development. We examined the sources of variability in the duration of individual embryo development, egg production and fertilization rate in Nothobranchius furzeri. Using a quantitative genetics approach (North Carolina type II design), we found support for maternal effects rather than polyandrous mating as the primary source of the variability in the duration of embryo development. The number of previously laid eggs appeared to serve as an internal physiological cue initiating a shift from rapid‐to‐slow embryo developmental mode. In annual killifish, extensive phenotypic variability in progeny traits is adaptive, as the conditions experienced by parents have limited relevance to the offspring generation. In contrast to genetic control, with high phenotypic expression and heritability, maternal control of traits under natural selection prevents standing genetic diversity from potentially detrimental effects of selection in fluctuating environments.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole‐genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.     

Short‐term insurance versus long‐term bet‐hedging strategies as adaptations to variable environments

Understanding how organisms adapt to environmental variation is a key challenge of biology. Central to this are bet‐hedging strategies that maximize geometric mean fitness across generations, either by being conservative or diversifying phenotypes. Theoretical models have identified environmental variation across generations with multiplicative fitness effects as driving the evolution of bet‐hedging. However, behavioral ecology has revealed adaptive responses to additive fitness effects of environmental variation within lifetimes, either through insurance or risk‐sensitive strategies. Here, we explore whether the effects of adaptive insurance interact with the evolution of bet‐hedging by varying the position and skew of both arithmetic and geometric mean fitness functions. We find that insurance causes the optimal phenotype to shift from the peak to down the less steeply decreasing side of the fitness function, and that conservative bet‐hedging produces an additional shift on top of this, which decreases as adaptive phenotypic variation from diversifying bet‐hedging increases. When diversifying bet‐hedging is not an option, environmental canalization to reduce phenotypic variation is almost always favored, except where the tails of the fitness function are steeply convex and produce a novel risk‐sensitive increase in phenotypic variance akin to diversifying bet‐hedging. Importantly, using skewed fitness functions, we provide the first model that explicitly addresses how conservative and diversifying bet‐hedging strategies might coexist.