Epigenetic inheritance, genetic assimilation and speciation.

Epigenetic inheritance systems enable the environmentally induced phenotypes to be transmitted between generations. Jablonka and Lamb (1991, 1995) proposed that these systems have a substantial role during speciation. They argued that divergence of isolated populations may be first triggered by the accumulation of (heritable) phenotypic differences that are later followed and strengthened by genetic changes. The plausibility of this idea is examined in this paper. At first, we discuss the "exploratory" behaviour of an epigenetic inheritance system on a one peak adaptive landscape. If a quantitative trait is far from the optimum, then it is advantageous to induce heritable phenotypic variation. Conversely, if the genotypes get closer to the peak, it is more favorable to canalize the phenotypic expression of the character. This process would lead to genetic assimilation. Next we show that the divergence of heritable epigenetic marks acts to reduce or to eliminate the genetic barrier between two adaptive peaks. Therefore, an epigenetic inheritance system can increase the probability of transition from one adaptive state to another. Peak shift might be initiated by (i) slight changes in the inducing environment or by (ii) genetic drift of the genes controlling epigenetic variability. Remarkably, drift-induced transition is facilitated even if phenotypic variation is not heritable. A corollary of our thesis is that evolution can proceed through suboptimal phenotypic states, without passing through a deep adaptive valley of the genotype. We also consider the consequences of this finding on the dynamics and mode of reproductive isolation.     

FUNCTIONAL GENETIC DIVERGENCE IN HIGH CO2 ADAPTED EMILIANIA HUXLEYI POPULATIONS

Predicting the impacts of environmental change on marine organisms, food webs, and biogeochemical cycles presently relies almost exclusively on short‐term physiological studies, while the possibility of adaptive evolution is often ignored. Here, we assess adaptive evolution in the coccolithophore Emiliania huxleyi, a well‐established model species in biological oceanography, in response to ocean acidification. We previously demonstrated that this globally important marine phytoplankton species adapts within 500 generations to elevated CO₂. After 750 and 1000 generations, no further fitness increase occurred, and we observed phenotypic convergence between replicate populations. We then exposed adapted populations to two novel environments to investigate whether or not the underlying basis for high CO₂‐adaptation involves functional genetic divergence, assuming that different novel mutations become apparent via divergent pleiotropic effects. The novel environment “high light” did not reveal such genetic divergence whereas growth in a low‐salinity environment revealed strong pleiotropic effects in high CO₂ adapted populations, indicating divergent genetic bases for adaptation to high CO₂. This suggests that pleiotropy plays an important role in adaptation of natural E. huxleyi populations to ocean acidification. Our study highlights the potential mutual benefits for oceanography and evolutionary biology of using ecologically important marine phytoplankton for microbial evolution experiments.     

Locally adapted traits maintained in the face of high gene flow

Gene flow between phenotypically divergent populations can disrupt local adaptation or, alternatively, may stimulate adaptive evolution by increasing genetic variation. We capitalised on historical Trinidadian guppy transplant experiments to test the phenotypic effects of increased gene flow caused by replicated introductions of adaptively divergent guppies, which were translocated from high‐ to low‐predation environments. We sampled two native populations prior to the onset of gene flow, six historic introduction sites, introduction sources and multiple downstream points in each basin. Extensive gene flow from introductions occurred in all streams, yet adaptive phenotypic divergence across a gradient in predation level was maintained. Descendants of guppies from a high‐predation source site showed high phenotypic similarity with native low‐predation guppies in as few as ~12 generations after gene flow, likely through a combination of adaptive evolution and phenotypic plasticity. Our results demonstrate that locally adapted phenotypes can be maintained despite extensive gene flow from divergent populations.     

Molecular phenotyping of maternally mediated parallel adaptive divergence within Rana arvalis and Rana temporaria

When similar selection acts on the same traits in multiple species or populations, parallel evolution can result in similar phenotypic changes, yet the underlying molecular architecture of parallel phenotypic divergence can be variable. Maternal effects can influence evolution at ecological timescales and facilitate local adaptation, but their contribution to parallel adaptive divergence is unclear. In this study, we (i) tested for variation in embryonic acid tolerance in a common garden experiment and (ii) used molecular phenotyping of egg coats to investigate the molecular basis of maternally mediated parallel adaptive divergence in two amphibian species (Rana arvalis and Rana temporaria). Our results on three R. arvalis and two R. temporaria populations show that adaptive divergence in embryonic acid tolerance is mediated via maternally derived egg coats in both species. We find extensive polymorphism in egg jelly coat glycoproteins within both species and that acid‐tolerant clutches have more negatively charged egg jelly – indicating that the glycosylation status of the jelly coat proteins is under divergent selection in acidified environments, likely due to its impact on jelly water balance. Overall, these data provide evidence for parallel mechanisms of adaptive divergence in two species. Our study highlights the importance of studying intraspecific molecular variation in egg coats and, specifically, their glycoproteins, to increase understanding of underlying forces maintaining variation in jelly coats.     

QUANTITATIVE GENETICS OF BRYOZOAN PHENOTYPIC EVOLUTION. I. RATE TESTS FOR RANDOM CHANGE VERSUS SELECTION IN DIFFERENTIATION OF LIVING SPECIES

The possible roles of random genetic change and natural selection in bryozoan speciation were analyzed using quantitative genetic methods on breeding data for traits of skeletal morphology in two closely related species of the cheilostome Stylopoma. The hypothesis that morphologic differences between the species are caused entirely by mutation and genetic drift could not be rejected for reasonable rates of mutation maintained for as few as 10³ to 10⁴ generations. Divergence times this short or shorter are consistent with the abrupt appearances of many invertebrate species in the fossil record, commonly followed by millions of years of morphologic stasis. To produce these differences over 10³ generations or fewer, directional selection acting alone would require unrealistically high levels of minimum selective mortality throughout divergence. Thus, selection is unnecessary to explain the divergence of these species, except as a means of accelerating the effects of random genetic change on shorter time scales (directional selection), or decelerating them over longer ones (stabilizing selection). These results are consistent with a variety of models of phenotypic evolution involving random shifts between multiple adaptive peaks. Similar results were obtained by substituting trait heritabilities and genetic covariances reconstructed by partitioning within- and among-colony phenotypic variance in place of the values based on breeding data. Quantitative genetic analysis of speciation in fossil bryozoan lineages is thus justified.     

Exploiting genomic resources in studies of speciation and adaptive radiation of lizards in the genus Anolis

Lizards in the genus Anolis have radiated extensively within and among islands in the Caribbean. Here, I provide a prospectus for identifying genes underlying adaptive phenotypic traits in anoles. First I review patterns of diversification in Anolis and the important morphological axes along which divergence occurs. Then I discuss two features of anole diversification, the repeated, convergent evolution of ecomorphs, and phenotypic divergence among populations within species, that provide opportunities to identify genes underlying adaptive phenotypic variation. While small clutch size and difficulty with captive rearing currently limit the utility of quantitative trait locus analyses, comparative analyses of gene expression, and population genomic approaches are promising.     

A multiple peak adaptive landscape based on feeding strategies and roosting ecology shaped the evolution of cranial covariance structure and morphological differentiation in phyllostomid bats

We explored the evolution of morphological integration in the most noteworthy example of adaptive radiation in mammals, the New World leaf‐nosed bats, using a massive dataset and by combining phylogenetic comparative methods and quantitative genetic approaches. We demonstrated that the phenotypic covariance structure remained conserved on a broader phylogenetic scale but also showed a substantial divergence between interclade comparisons. Most of the phylogenetic structure in the integration space can be explained by splits at the beginning of the diversification of major clades. Our results provide evidence for a multiple peak adaptive landscape in the evolution of cranial covariance structure and morphological differentiation, based upon diet and roosting ecology. In this scenario, the successful radiation of phyllostomid bats was triggered by the diversification of dietary and roosting strategies, and the invasion of these new adaptive zones lead to changes in phenotypic covariance structure and average morphology. Our results suggest that intense natural selection preceded the invasion of these new adaptive zones and played a fundamental role in shaping cranial covariance structure and morphological differentiation in this hyperdiverse clade of mammals. Finally, our study demonstrates the power of combining comparative methods and quantitative genetic approaches when investigating the evolution of complex morphologies.     

Evaluating the role of natural selection in the evolution of gene regulation

Surveys of gene expression reveal extensive variability both within and between a wide range of species. Compelling cases have been made for adaptive changes in gene regulation, but the proportion of expression divergence attributable to natural selection remains unclear. Distinguishing adaptive changes driven by positive selection from neutral divergence resulting from mutation and genetic drift is critical for understanding the evolution of gene expression. Here, we review the various methods that have been used to test for signs of selection in genomic expression data. We also discuss properties of regulatory systems relevant to neutral models of gene expression. Despite some potential caveats, published studies provide considerable evidence for adaptive changes in gene expression. Future challenges for studies of regulatory evolution will be to quantify the frequency of adaptive changes, identify the genetic basis of expression divergence and associate changes in gene expression with specific organismal phenotypes.     

ADAPTIVE RADIATION ALONG GENETIC LINES OF LEAST RESISTANCE

Are measurements of quantitative genetic variation useful for predicting long-term adaptive evolution? To answer this question, I focus on g_max, the multivariate direction of greatest additive genetic variance within populations. Original data on threespine sticklebacks, together with published genetic measurements from other vertebrates, show that morphological differentiation between species has been biased in the direction of g_max for at least four million years, despite evidence that natural selection is the cause of differentiation. This bias toward the direction of evolution tends to decay with time. Rate of morphological divergence between species is inversely proportional to θ, the angle between the direction of divergence and the direction of greatest genetic variation. The direction of greatest phenotypic variance is not identical with g_max, but for these data is nearly as successful at predicting the direction of species divergence. I interpret the findings to mean that genetic variances and covariances constrain adaptive change in quantitative traits for reasonably long spans of time. An alternative hypothesis, however, cannot be ruled out: that morphological differentiation is biased in the direction g_max because divergence and g_max are both shaped by the same natural selection pressures. Either way, the results reveal that adaptive differentiation occurs principally along “genetic lines of least resistance.”     

Concordance between stabilizing sexual selection,intraspecific variation,and interspecific divergence in Phymata

Empirical studies show that lineages typically exhibit long periods of evolutionary stasis and that relative levels of within‐species trait covariance often correlate with the extent of between‐species trait divergence. These observations have been interpreted by some as evidence of genetic constraints persisting for long periods of time. However, an alternative explanation is that both intra‐ and interspecific variation are shaped by the features of the adaptive landscape (e.g., stabilizing selection). Employing a genus of insects that are diverse with respect to a suite of secondary sex traits, we related data describing nonlinear phenotypic (sexual) selection to intraspecific trait covariances and macroevolutionary divergence. We found support for two key predictions (1) that intraspecific trait covariation would be aligned with stabilizing selection and (2) that there would be restricted macroevolutionary divergence in the direction of stabilizing selection. The observed alignment of all three matrices offers a point of caution in interpreting standing variability as metrics of evolutionary constraint. Our results also illustrate the power of sexual selection for determining variation observed at both short and long timescales and account for the apparently slow evolution of some secondary sex characters in this lineage.     

Multivariate Phenotypic Divergence Due to the Fixation of Beneficial Mutations in Experimentally Evolved Lineages of a Filamentous Fungus

The potential for evolutionary change is limited by the availability of genetic variation. Mutations are the ultimate source of new alleles, yet there have been few experimental investigations of the role of novel mutations in multivariate phenotypic evolution. Here, we evaluated the degree of multivariate phenotypic divergence observed in a long-term evolution experiment whereby replicate lineages of the filamentous fungus Aspergillus nidulans were derived from a single genotype and allowed to fix novel (beneficial) mutations while maintained at two different population sizes. We asked three fundamental questions regarding phenotypic divergence following approximately 800 generations of adaptation: (1) whether divergence was limited by mutational supply, (2) whether divergence proceeded in relatively many (few) multivariate directions, and (3) to what degree phenotypic divergence scaled with changes in fitness (i.e. adaptation). We found no evidence that mutational supply limited phenotypic divergence. Divergence also occurred in all possible phenotypic directions, implying that pleiotropy was either weak or sufficiently variable among new mutations so as not to constrain the direction of multivariate evolution. The degree of total phenotypic divergence from the common ancestor was positively correlated with the extent of adaptation. These results are discussed in the context of the evolution of complex phenotypes through the input of adaptive mutations.     

Divergence and evolution of assortative mating in a polygenic trait model of speciation with gene flow

Assortative mating is an important driver of speciation in populations with gene flow and is predicted to evolve under certain conditions in few‐locus models. However, the evolution of assortment is less understood for mating based on quantitative traits, which are often characterized by high genetic variability and extensive linkage disequilibrium between trait loci. We explore this scenario for a two‐deme model with migration, by considering a single polygenic trait subject to divergent viability selection across demes, as well as assortative mating and sexual selection within demes, and investigate how trait divergence is shaped by various evolutionary forces. Our analysis reveals the existence of sharp thresholds of assortment strength, at which divergence increases dramatically. We also study the evolution of assortment via invasion of modifiers of mate discrimination and show that the ES assortment strength has an intermediate value under a range of migration‐selection parameters, even in diverged populations, due to subtle effects which depend sensitively on the extent of phenotypic variation within these populations. The evolutionary dynamics of the polygenic trait is studied using the hypergeometric and infinitesimal models. We further investigate the sensitivity of our results to the assumptions of the hypergeometric model, using individual‐based simulations.     

Diversification on an ecologically constrained adaptive landscape

We used phylogenetic analysis of body-size ecomorphs in a crustacean species complex to gain insight into how spatial complexity of ecological processes generates and maintains biological diversity. Studies of geographically widespread species of Hyalella amphipods show that phenotypic evolution is tightly constrained in a manner consistent with adaptive responses to alternative predation regimes. A molecular phylogeny indicates that evolution of Hyalella ecomorphs is characterized by parallel evolution and by phenotypic stasis despite substantial levels of underlying molecular change. The phylogeny suggests that species diversification sometimes occurs by niche shifts, and sometimes occurs without a change in niche. Moreover, diversification in the Hyalella ecomorphs has involved the repeated evolution of similar phenotypic forms that exist in similar ecological settings, a hallmark of adaptive evolution. The evolutionary stasis observed in clades separated by substantial genetic divergence, but existing in similar habitats, is also suggestive of stabilizing natural selection acting to constrain phenotypic evolution within narrow bounds. We interpret the observed decoupling of genetic and phenotypic diversification in terms of adaptive radiation on an ecologically constrained adaptive landscape, and suggest that ecological constraints, perhaps acting together with genetic and functional constraints, may explain the parallel evolution and evolutionary stasis inferred by the phylogeny.     

HOW STABLE ‘SHOULD’ EPIGENETIC MODIFICATIONS BE? INSIGHTS FROM ADAPTIVE PLASTICITY AND BET HEDGING

Although there is keen interest in the potential adaptive value of epigenetic variation, it is unclear what conditions favor the stability of these variants either within or across generations. Because epigenetic modifications can be environmentally sensitive, existing theory on adaptive phenotypic plasticity provides relevant insights. Our consideration of this theory suggests that stable maintenance of environmentally induced epigenetic states over an organism's lifetime is most likely to be favored when the organism accurately responds to a single environmental change that subsequently remains constant, or when the environmental change cues an irreversible developmental transition. Stable transmission of adaptive epigenetic states from parents to offspring may be selectively favored when environments vary across generations and the parental environment predicts the offspring environment. The adaptive value of stability beyond a single generation of parent–offspring transmission likely depends on the costs of epigenetic resetting. Epigenetic stability both within and across generations will also depend on the degree and predictability of environmental variation, dispersal patterns, and the (epi)genetic architecture underlying phenotypic responses to environment. We also discuss conditions that favor stability of random epigenetic variants within the context of bet hedging. We conclude by proposing research directions to clarify the adaptive significance of epigenetic stability.     

Population mixing and the adaptive divergence of quantitative traits in discrete populations: a theoretical framework for empirical tests

Empirical tests for the importance of population mixing in constraining adaptive divergence have not been well grounded in theory for quantitative traits in spatially discrete populations. We develop quantitative-genetic models to examine the equilibrium difference between two populations that are experiencing different selective regimes and exchanging individuals. These models demonstrate that adaptive divergence is negatively correlated with the rate of population mixing (m, most strongly so when m is low), positively correlated with the difference in phenotypic optima between populations, and positively correlated with the amount of additive genetic variance (G, most strongly so when G is low). The approach to equilibrium is quite rapid (fewer than 50 generations for two populations to evolve 90% of the distance to equilibrium) when either heritability or mixing are not too low (h2 > 0.2 or m > 0.05). The theory can be used to aid empirical tests that: (1) compare observed divergence to that predicted using estimates of population mixing, additive genetic variance/covariance, and selection; (2) test for a negative correlation between population mixing and adaptive divergence across multiple independent population pairs; and (3) experimentally manipulate the rate of mixing. Application of the first two of these approaches to data from two well-studied natural systems suggests that population mixing has constrained adaptive divergence for color patterns in Lake Erie water snakes (Nerodia sipedon), but not for trophic traits in sympatric pairs of benthic and limnetic stickleback (Gasterosteus aculeatus). The theoretical framework we outline should provide an improved basis for future empirical tests of the role of population mixing in adaptive divergence.     

Evolution of adaptive diversity and genetic connectivity in Arctic charr (Salvelinus alpinus) in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.     

Testing limits to adaptation along altitudinal gradients in rainforest Drosophila

Given that evolution can generate rapid and dramatic shifts in the ecological tolerance of a species, what prevents populations adapting to expand into new habitat at the edge of their distributions? Recent population genetic models have focused on the relative costs and benefits of migration between populations. On the one hand, migration may limit adaptive divergence by preventing local populations from matching their local selective optima. On the other hand, migration may also contribute to the genetic variance necessary to allow populations to track these changing optima. Empirical evidence for these contrasting effects of gene flow in natural situations are lacking, largely because it remains difficult to acquire. Here, we develop a way to explore theoretical models by estimating genetic divergence in traits that confer stress resistance along similar ecological gradients in rainforest Drosophila. This approach allows testing for the coupling of clinal divergence with local density, and the effects of genetic variance and the rate of change of the optimum on the response to selection. In support of a swamping effect of migration on phenotypic divergence, our data show no evidence for a cline in stress-related traits where the altitudinal gradient is steep, but significant clinal divergence where it is shallow. However, where clinal divergence is detected, sites showing trait means closer to the presumed local optimum have more genetic variation than sites with trait means distant from their local optimum. This pattern suggests that gene flow also aids a sustained response to selection.     

Measuring the contribution of evolution to community trait structure in freshwater zooplankton

There are currently few predictions about when evolutionary processes are likely to play an important role in structuring community features. Determining predictors that indicate when evolution is expected to impact ecological processes in natural landscapes can help researchers identify eco-evolutionary ‘hotspots', where eco-evolutionary interactions are more likely to occur. Using data collected from a survey in freshwater cladoceran communities, landscape population genetic data and phenotypic trait data measured in a common garden, we applied a Bayesian linear model to assess whether the impact of local trait evolution in the keystone species Daphnia magna on cladoceran community trait values could be predicted by population genetic properties (within-population genetic diversity, genetic distance among populations), ecological properties (Simpson's diversity, phenotypic divergence) or environmental divergence. We found that the impact of local trait evolution varied among communities. Moreover, community diversity and phenotypic divergence were found to be better predictors of the contribution of evolution to community trait values than environmental features or genetic properties of the evolving species. Our results thus indicate the importance of ecological context for the impact of evolution on community features. Our study also demonstrates one way to detect signatures of eco-evolutionary interactions in communities inhabiting heterogeneous landscapes using survey data of contemporary ecological and evolutionary structure.     

The long-term evolution of multilocus traits under frequency-dependent disruptive selection

Frequency-dependent disruptive selection is widely recognized as an important source of genetic variation. Its evolutionary consequences have been extensively studied using phenotypic evolutionary models, based on quantitative genetics, game theory, or adaptive dynamics. However, the genetic assumptions underlying these approaches are highly idealized and, even worse, predict different consequences of frequency-dependent disruptive selection. Population genetic models, by contrast, enable genotypic evolutionary models, but traditionally assume constant fitness values. Only a minority of these models thus addresses frequency-dependent selection, and only a few of these do so in a multilocus context. An inherent limitation of these remaining studies is that they only investigate the short-term maintenance of genetic variation. Consequently, the long-term evolution of multilocus characters under frequency-dependent disruptive selection remains poorly understood. We aim to bridge this gap between phenotypic and genotypic models by studying a multilocus version of Levene's soft-selection model. Individual-based simulations and deterministic approximations based on adaptive dynamics theory provide insights into the underlying evolutionary dynamics. Our analysis uncovers a general pattern of polymorphism formation and collapse, likely to apply to a wide variety of genetic systems: after convergence to a fitness minimum and the subsequent establishment of genetic polymorphism at multiple loci, genetic variation becomes increasingly concentrated on a few loci, until eventually only a single polymorphic locus remains. This evolutionary process combines features observed in quantitative genetics and adaptive dynamics models, and it can be explained as a consequence of changes in the selection regime that are inherent to frequency-dependent disruptive selection. Our findings demonstrate that the potential of frequency-dependent disruptive selection to maintain polygenic variation is considerably smaller than previously expected.     

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.