Quantitative genetic variation in static allometry in the threespine stickleback

The common pattern of replicated evolution of a consistent shape-environment relationship might reflect selection acting in similar ways within each environment, but divergently among environments. However, phenotypic evolution depends on the availability of additive genetic variation as well as on the direction of selection, implicating a bias in the distribution of genetic variance as a potential contributor to replicated evolution. Allometry, the relationship between shape and size, is a potential source of genetic bias that is poorly understood. The threespine stickleback, Gasterosteus aculeatus, provides an ideal system for exploring the contribution of genetic variance in body shape allometry to evolutionary patterns. The stickleback system comprises marine populations that exhibit limited phenotypic variation, and young freshwater populations which, following independent colonization events, have often evolved similar phenotypes in similar environments. In particular, stickleback diversification has involved changes in both total body size and relative size of body regions (i.e., shape). In a laboratory-reared cohort derived from an oceanic Alaskan population that is phenotypically and genetically representative of the ancestor of the diverse freshwater populations in this region, we determined the phenotypic static allometry, and estimated the additive genetic variation about these population-level allometric functions. We detected significant allometry, with larger fish having relatively smaller heads, a longer base to their second dorsal fin, and longer, shallower caudal peduncles. There was additive genetic variance in body size and in size-independent body shape (i.e., allometric elevation), but typically not in allometric slopes. These results suggest that the parallel evolution of body shape in threespine stickleback is not likely to have been a correlated response to selection on body size, or vice versa. Although allometry is common in fishes, this study highlights the need for additional data on genetic variation in allometric functions to determine how allometry evolves and how it influences phenotypic evolution.     

Genetics of body shape and armour variation in threespine sticklebacks

Patterns of genetic variation and covariation can influence the rate and direction of phenotypic evolution. We explored the possibility that the parallel morphological evolution seen in threespine stickleback (Gasterosteus aculeatus) populations colonizing freshwater environments is facilitated by patterns of genetic variation and covariation in the ancestral (marine) population. We estimated the genetic (G) and phenotypic (P) covariance matrices and directions of maximum additive genetic (g(max) ) and phenotypic (p(max) ) covariances of body shape and armour traits. Our results suggest a role for the ancestral G in explaining parallel morphological evolution in freshwater populations. We also found evidence of genetic constraints owing to the lack of variance in the ancestral G. Furthermore, strong genetic covariances and correlations among traits revealed that selective factors responsible for threespine stickleback body shape and armour divergence may be difficult to disentangle. The directions of g(max) and p(max) were correlated, but the correlations were not high enough to imply that phenotypic patterns of trait variation and covariation within populations are very informative of underlying genetic patterns.     

Conserved phenotypic variation patterns, evolution along lines of least resistance, and departure due to selection in fossil rodents

Within a group of organisms, some morphologies are more readily generated than others due to internal developmental constraints. Such constraints can channel evolutionary changes into directions corresponding to the greatest intraspecific variation. Long-term evolutionary outputs, however, depend on the stability of these intraspecific patterns of variation over time and from the interplay between internal constraints and selective regimes. To address these questions, the relationship between the structure of phenotypic variance covariance matrices and direction of morphological evolution was investigated using teeth of fossil rodents. One lineage considered here leads to Stephanomys, a highly specialized genus characterized by a dental pattern supposedly favoring grass eating. Stephanomys evolved in the context of directional selection related to the climatic trend of global cooling causing an increasing proportion of grasslands in southwestern Europe. The initial divergence (up to approximately 6.5 mya) was channeled along the direction of greatest intraspecific variation, whereas after 6.5 mya, morphological evolution departed from the direction favored by internal constraints. This departure from the "lines of least resistance" was likely the consequence of an environmental degradation causing a selective gradient strong enough to overwhelm the constraints to phenotypic evolution. However, in a context of stabilizing selection, these constraints actually channel evolution, as exemplified by the lineage of Apodemus. This lineage retained a primitive diet and dental pattern over the last 10 myr. Limited morphological changes occurred nevertheless in accordance with the main patterns of intraspecific variation. The importance of these lines of least resistance directing long-term morphological evolution may explain parallel evolution of some dental patterns in murine evolution.     

PLASTICITY, CANALIZATION, AND DEVELOPMENTAL STABILITY OF THE DROSOPHILA WING: JOINT EFFECTS OF MUTATIONS AND DEVELOPMENTAL TEMPERATURE

The phenotypic effects of genetic and environmental manipulations have been rarely investigated simultaneously. In addition to phenotypic plasticity, their effect on the amount and directions of genetic and phenotypic variation is of particular evolutionary importance because these constitute the material for natural selection. Here, we used heterozygous insertional mutations of 16 genes involved in the formation of the Drosophila wing. The flies were raised at two developmental temperatures (18°C and 28°C). Landmark-based geometric morphometrics was used to analyze the variation of the wing size and shape at different hierarchical levels: among genotypes and temperatures; among individuals within group; and fluctuating asymmetry (FA). Our results show that (1) the phenotypic effects of the mutations depend on temperature; (2) reciprocally, most mutations affect wing plasticity; (3) both temperature and mutations modify the levels of FA and of among individuals variation within lines. Remarkably, the patterns of shape FA seem unaffected by temperature whereas those associated with individual variation are systematically altered. By modifying the direction of available phenotypic variation, temperature might thus directly affect the potential for further evolution. It suggests as well that the developmental processes responsible for developmental stability and environmental canalization might be partially distinct.     

Developmental constraints revealed by co-variation within and among molar rows in two murine rodents

SUMMARY Morphological integration corresponds to interdependency between characters that can arise from several causes. Proximal causes of integration include that different phenotypic features may share common genetic sets and/or interact during their development. Ultimate causes may be the prolonged effect of selection favoring integration of functionally interacting characters, achieved by the molding of these proximal causes. Strong and direct interactions among successive teeth of a molar row are predicted by genetic and developmental evidences. Functional constraints related to occlusion, however, should have selected more strongly for a morphological integration of occluding teeth and a corresponding evolution of the underlying developmental and genetic pathways. To investigate how these predictions match the patterns of phenotypic integration, we studied the co‐variation among the six molars of the murine molar row, focusing on two populations of house mice (Mus musculus domesticus) and wood mice (Apodemus sylvaticus). The size and shape of the three upper and lower molars were quantified and compared. Our results evidenced similar patterns in both species, size being more integrated than shape among all the teeth, and both size and shape co‐varying strongly between adjacent teeth, but also between occluding teeth. Strong co‐variation within each molar row is in agreement with developmental models showing a cascade influence of the first molar on the subsequent molars. In contrast, the strong co‐variation between molars of the occluding tooth rows confirms that functional constraints molded patterns of integration and probably the underlying developmental pathways despite the low level of direct developmental interactions occurring among molar rows. These patterns of co‐variation are furthermore conserved between the house mouse and the wood mouse that diverged >10 Ma, suggesting that they may constitute long‐running constraints to the diversification of the murine rodent dentition.     

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.     

Developmental plasticity and the origin of novel forms: unveiling cryptic genetic variation via "use and disuse"

Natural selection eliminates phenotypic variation from populations, generation after generation-an observation that haunted Darwin. So, how does new phenotypic variation arise, and is it always random with respect to fitness? Repeated behavioral responses to a novel environment-particularly those that are learned-are typically advantageous. If those behaviors yield more extreme or novel morphological variants via developmental plasticity, then previously cryptic genetic variation may be exposed to natural selection. Significantly, because the mean phenotypic effect of "use and disuse" is also typically favorable, previously cryptic genetic variation can be transformed into phenotypic variation that is both visible to selection and biased in an adaptive direction. Therefore, use-induced developmental plasticity in a very real sense "creates" new phenotypic variation that is nonrandom with respect to fitness, in contrast to the random phenotypic effects of mutation, recombination, and "direct effects" of environment (stress, nutrition). I offer here (a) a brief review of the immense literature on the effects of "use and disuse" on morphology, (b) a simple yet general model illustrating how cryptic genetic variation may be exposed to selection by developmentally plastic responses that alter trait performance in response to "use and disuse," and (c) a more detailed model of a positive feedback loop between learning (handed behavior) and morphological plasticity (use-induced morphological asymmetry) that may rapidly generate novel phenotypic variation and facilitate the evolution of conspicuous morphological asymmetries. Evidence from several sources suggests that handed behaviors played an important role both in the origin of novel forms (asymmetries) and in their subsequent evolution.     

Alternative life histories in the Atlantic salmon: genetic covariances within the sneaker sexual tactic in males

Alternative reproductive tactics are ubiquitous in many species. Tactic expression often depends on whether an individual's condition surpasses thresholds that are responsible for activating particular developmental pathways. Two central goals in understanding the evolution of reproductive tactics are quantifying the extent to which thresholds are explained by additive genetic effects, and describing their covariation with condition-related traits. We monitored the development of early sexual maturation that leads to the sneaker reproductive tactic in Atlantic salmon (Salmo salar L.). We found evidence for additive genetic variance in the timing of sexual maturity (which is a measure of the surpassing of threshold values) and body-size traits. This suggests that selection can affect the patterns of sexual development by changing the timing of this event and/or body size. Significant levels of covariation between these traits also occurred, implying a potential for correlated responses to selection. Closer examination of genetic covariances suggests that the detected genetic variation is distributed along at least five directions of phenotypic variation. Our results show that the potential for evolution of the life-history traits constituting this reproductive phenotype is greatly influenced by their patterns of genetic covariance.     

METAMODELS AND PHYLOGENETIC REPLICATION: A SYSTEMATIC APPROACH TO THE EVOLUTION OF DEVELOPMENTAL PATHWAYS

Molecular genetic analysis of phenotypic variation has revealed many examples of evolutionary change in the developmental pathways that control plant and animal morphology. A major challenge is to integrate the information from diverse organisms and traits to understand the general patterns of developmental evolution. This integration can be facilitated by evolutionary metamodels—traits that have undergone multiple independent changes in different species and whose development is controlled by well-studied regulatory pathways. The metamodel approach provides the comparative equivalent of experimental replication, allowing us to test whether the evolution of each developmental pathway follows a consistent pattern, and whether different pathways are predisposed to different modes of evolution by their intrinsic organization. A review of several metamodels suggests that the structure of developmental pathways may bias the genetic basis of phenotypic evolution, and highlights phylogenetic replication as a value-added approach that produces deeper insights into the mechanisms of evolution than single-species analyses.     

SIZE AS A LINE OF LEAST RESISTANCE II: DIRECT SELECTION ON SIZE OR CORRELATED RESPONSE DUE TO CONSTRAINTS?

Evolutionary change in New World Monkey (NWM) skulls occurred primarily along the line of least resistance defined by size (including allometric) variation ( g_max ). Although the direction of evolution was aligned with this axis, it was not clear whether this macroevolutionary pattern results from the conservation of within population genetic covariance patterns (long‐term constraint) or long‐term selection along a size dimension, or whether both, constraints and selection, were inextricably involved. Furthermore, G ‐matrix stability can also be a consequence of selection, which implies that both, constraints embodied in g_max and evolutionary changes observed on the trait averages, would be influenced by selection. Here, we describe a combination of approaches that allows one to test whether any particular instance of size evolution is a correlated by‐product due to constraints ( g_max ) or is due to direct selection on size and apply it to NWM lineages as a case study. The approach is based on comparing the direction and amount of evolutionary change produced by two different simulated sets of net‐selection gradients ( β ), a size (isometric and allometric size) and a nonsize set. Using this approach it is possible to distinguish between the two hypotheses (indirect size evolution due to constraints or direct selection on size), because although both may produce an evolutionary response aligned with g_max , the amount of change produced by random selection operating through the variance/covariance patterns (constraints hypothesis) will be much smaller than that produced by selection on size (selection hypothesis). Furthermore, the alignment of simulated evolutionary changes with g_max when selection is not on size is not as tight as when selection is actually on size, allowing a statistical test of whether a particular observed case of evolution along the line of least resistance is the result of selection along it or not. Also, with matrix diagonalization (principal components [PC]) it is possible to calculate directly the net‐selection gradient on size alone (first PC [PC1]) by dividing the amount of phenotypic difference between any two populations by the amount of variation in PC1, which allows one to benchmark whether selection was on size or not.     

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.     

Selection on Wing Allometry in Drosophila Melanogaster

Five bivariate distributions of wing dimensions of Drosophila melanogaster were measured, in flies 1) subjected to four defined environmental regimes during development, 2) taken directly from nature in seven U.S. states, 3) selected in ten populations for change in wing form, and 4) sampled from 21 long inbred wild-type lines. Environmental stresses during development altered both wing size and the ratios of wing dimensions, but regardless of treatment all wing dimensions fell near a common allometric baseline in each bivariate distribution. The wings of wild-caught flies from seven widely separated localities, and of their laboratory-reared offspring, also fell along the same baselines. However, when flies were selected divergently for lateral offset from these developmental baselines, response to selection was rapid in every case. The mean divergence in offset between oppositely selected lines was 14.68 SD of the base population offset, after only 15 generations of selection at 20%. Measurements of 21 isofemale lines, founded from wild-caught flies and maintained in small populations for at least 22 years, showed large reductions in phenotypic variance of offsets within lines, but a large increase in the variance among lines. The variance of means of isofemale lines within collection localities was ten times the variance of means among localities of newly established wild lines. These observations show that much additive genetic variance exists for individual dimensions within the wing, such that bivariate developmental patterns can be changed in any direction by selection or by drift. The relative invariance of the allometric baselines of wing morphology in nature is most easily explained as the result of continuous natural selection around a local optimum of functional design.     

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.     

Measuring and comparing evolvability and constraint in multivariate characters

The Lande equation forms the basis for our understanding of the short-term evolution of quantitative traits in a multivariate context. It predicts the response to selection as the product of an additive genetic variance matrix and a selection gradient. The selection gradient approximates the force and direction of selection, and the genetic variance matrix quantifies the role of the genetic system in evolution. Attempts to understand the evolutionary significance of the genetic variance matrix are hampered by the fact that the majority of the methods used to characterize and compare variance matrices have not been derived in an explicit theoretical context. We use the Lande equation to derive new measures of the ability of a variance matrix to allow or constrain evolution in any direction in phenotype space. Evolvability captures the ability of a population to evolve in the direction of selection when stabilizing selection is absent. Conditional evolvability captures the ability of a population to respond to directional selection in the presence of stabilizing selection on other trait combinations. We then derive measures of character autonomy and integration from these evolvabilities. We study the properties of these measures and show how they can be used to interpret and compare variance matrices. As an illustration, we show that divergence of wing shape in the dipteran family Drosophilidae has proceeded in directions that have relatively high evolvabilities.     

Change of Genetic Architecture in Response to Sex

A traditional view is that sexual reproduction increases the potential for phenotypic evolution by expanding the range of genetic variation upon which natural selection can act. However, when nonadditive genetic effects and genetic disequilibria underlie a genetic system, genetic slippage (a change in the mean genotypic value contrary to that promoted by selection) in response to sex may occur. Additionally, depending on whether natural selection is predominantly stabilizing or disruptive, recombination may either enhance or reduce the level of expressed genetic variance. Thus, the role of sexual reproduction in the dynamics of phenotypic evolution depends heavily upon the nature of natural selection and the genetic system of the study population. In the present study, on a permanent lake Daphnia pulicaria population, sexual reproduction resulted in significant genetic slippage and a significant increase in expressed genetic variance for several traits. These observations provide evidence for substantial genetic disequilibria and nonadditive genetic effects underlying the genetic system of the study population. From these results, the fitness function of the previous clonal selection phase is inferred to be directional and/or stabilizing. The data are also used to infer the effects of natural selection on the mean and the genetic variance of the population.     

High evolutionary constraints limited adaptive responses to past climate changes in toad skulls

Interactions among traits that build a complex structure may be represented as genetic covariation and correlation. Genetic correlations may act as constraints, deflecting the evolutionary response from the direction of natural selection. We investigated the relative importance of drift, selection, and constraints in driving skull divergence in a group of related toad species. The distributional range of these species encompasses very distinct habitats with important climatic differences and the species are primarily distinguished by differences in their skulls. Some parts of the toad skull, such as the snout, may have functional relevance in reproductive ecology, detecting water cues. Thus, we hypothesized that the species skull divergence was driven by natural selection associated with climatic variation. However, given that all species present high correlations among skull traits, our second prediction was of high constraints deflecting the response to selection. We first extracted the main morphological direction that is expected to be subjected to selection by using within- and between-species covariance matrices. We then used evolutionary regressions to investigate whether divergence along this direction is explained by climatic variation between species. We also used quantitative genetics models to test for a role of random drift versus natural selection in skull divergence and to reconstruct selection gradients along species phylogeny. Climatic variables explained high proportions of between-species variation in the most selected axis. However, most evolutionary responses were not in the direction of selection, but aligned with the direction of allometric size, the dimension of highest phenotypic variance in the ancestral population. We conclude that toad species have responded to selection related to climate in their skulls, yet high evolutionary constraints dominated species divergence and may limit species responses to future climate change.     

How many processes are responsible for phenotypic evolution?

SUMMARY In addressing phenotypic evolution, this article reconsiders natural selection, random drift, developmental constraints, and internal selection in the new extended context of evolutionary developmental biology. The change of perspective from the "evolution of phenotypes" toward an "evolution of ontogenies" (evo-devo perspective) affects the reciprocal relationships among these different processes. Random drift and natural selection are sibling processes: two forms of post-productional sorting among alternative developmental trajectories, the former random, the latter nonrandom. Developmental constraint is a compound concept; it contains even some forms of natural ("external" and "internal") selection. A narrower definition ("reproductive constraints") is proposed. Internal selection is not a selection caused by an internal agent. It is a form of environment-independent selection depending on the level of the organism's internal developmental or functional coordination. Selection and constraints are the main deterministic processes in phenotypic evolution but they are not opposing forces. Indeed, they are continuously interacting processes of evolutionary change, but with different roles that should not be confused.     

The Baldwin effect and genetic assimilation: revisiting two mechanisms of evolutionary change mediated by phenotypic plasticity

Two different, but related, evolutionary theories pertaining to phenotypic plasticity were proposed by James Mark Baldwin and Conrad Hal Waddington. Unfortunately, these theories are often confused with one another. Baldwin's notion of organic selection posits that plasticity influences whether an individual will survive in a new environment, thus dictating the course of future evolution. Heritable variations can then be selected upon to direct phenotypic evolution (i.e., "orthoplasy"). The combination of these two processes (organic selection and orthoplasy) is now commonly referred to as the "Baldwin effect." Alternately, Waddington's genetic assimilation is a process whereby an environmentally induced phenotype, or "acquired character," becomes canalized through selection acting upon the developmental system. Genetic accommodation is a modern term used to describe the process of heritable changes that occur in response to a novel induction. Genetic accommodation is a key component of the Baldwin effect, and genetic assimilation is a type of genetic accommodation. I here define both the Baldwin effect and genetic assimilation in terms of genetic accommodation, describe cases in which either should occur in nature, and propose that each could play a role in evolutionary diversification.     

Geometric morphometrics and leaf phenotypic plasticity: assessing fluctuating asymmetry and allometry in European white oaks (Quercus)

Because plants are unable to move away from unfavourable habitats and environmental perturbations, leaf phenotypic plasticity facilitates light absorption and gas exchange. Oaks (Quercus spp.) are particularly known for their adaptability and plastic phenotypes, and leaf allometry and developmental instability may represent important mechanisms for their adaptation to environments and evolution. Because of its important role in the adaptation of plant populations to different environments, allometry can be involved in diversifying selection. Developmental instability is related to environmental perturbations and stresses by producing random deviations in structures characterized by bilateral symmetry, such as oak leaves. In addition, developmental instability can also arise from genetic bottlenecks or as a result of hybridization. The splitting of symmetric and asymmetric components of variation and their separate analysis allows the variability in leaf shape traits to be summarized, reducing the variation produced by developmental instability. The geometric morphometric approach is a useful method for the study of leaf asymmetry and allometric patterns. This method provides an important tool for the visualization of shape attributes that characterize species with highly variable leaf phenotypic patterns. In this study, leaf shape and size variability of three white oak species was investigated by means of a two‐dimensional landmark‐based method providing improved knowledge of variance partitioning, species discrimination, fluctuating asymmetry and allometric patterns of variation resulting from the different analyses. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 335–348.