Natural or artificial selection that favors higher values of a particular trait within a given population should engender an evolutionary response that increases the mean value of the trait. For this prediction to hold, the phenotypic variance of the trait must be caused in part by additive effects of alleles segregating in the population, and also the trait must not be too strongly genetically correlated with other traits that are under selection. Another prediction, rarely discussed in the literature, is that directional selection should favor alleles that increase phenotypic plasticity in the direction of selection, where phenotypic plasticity is defined as the ability of one genotype to produce more than one phenotype when exposed to different environments. This prediction has received relatively little empirical attention. Nonetheless, many laboratory experiments impose selection regimes that could allow for the evolution of enhanced plasticity (e.g. desiccation trials with Drosophila that last for several hours or days). We review one example that involved culturing of Drosophila on lemon for multiple generations and then tested for enhanced plasticity of detoxifying enzymes. We also review an example with vertebrates that involves selective breeding for high voluntary activity levels in house mice, targeting wheel-running behavior on days 5+6 of a 6-day wheel exposure. This selection regime allows for the possibility of wheel running itself or subordinate traits that support such running to increase in plasticity over days 1–4 of wheel access. Indeed, some traits, such as the concentration of the glucose transporter GLUT4 in gastrocnemius muscle, do show enhanced plasticity in the selected lines over a 5–6 day period. In several experiments we have housed mice from both the Selected (S) and Control (C) lines with or without wheel access for several weeks to test for differences in plasticity (training effects). A variety of patterns were observed, including no training effects in either S or C mice, similar changes in both the S and C lines, greater changes in the S lines but in the same direction in the C lines, and even opposite directions of change in the S and C lines. For some of the traits that show a greater training effect in the S lines, but in the same direction as in C lines, the greater effect can be explained statistically by the greater wheel running exhibited by S lines (more pain, more gain'). For others, however, the differences seem to reflect inherently greater plasticity in the S lines (i.e. for a given amount of stimulus, such as wheel running/day, individuals in the S lines show a greater response as compared with individuals in the C lines). We suggest that any selection experiment in which the selective event is more than instantaneous should explore whether plasticity in the appropriate (adaptive) direction has increased as a component of the response to selection.

Natural selection tends to act most strongly on aspects of the phenotype(traits) at relatively high levels of biological organization because they are the most strongly correlated with Darwinian fitness (e.g. lifetime reproductive success). Components of life history, behaviors and aspects of organismal performance (for reviews, see Ketterson and Nolan, Jr, 1999; Irschick and Garland, Jr,2001; Kingsolver and Huey,2003; Costa and Sinervo,2004) are complex traits' in that they are composed of many subordinate traits at lower levels of biological organization(Swallow and Garland, Jr,2005) (Fig. 1). Thus, the evolutionary response to selection on such complex traits necessarily entails associated changes in aspects of morphology, physiology and biochemical pathways (Ghalambor et al., 2003; Sinervo and Calsbeek, 2003). In addition, complex patterns of trade-offs and constraints are expected to occur, and the genetic architecture underlying these may itself evolve in response to selection (e.g. Chippindale et al., 2003; Rose et al., 2005).

Regardless of position in the biological hierarchy, most if not all traits are affected by alleles segregating in the population at many loci, whose expression is affected by numerous environmental factors, including both abiotic (e.g. temperature) and biotic (e.g. social interactions). Environmental factors can influence development by acting at any time after formation of the zygote, or in some cases even before (e.g. maternal effects acting on the unfertilized egg). Moreover, organisms often exert some level of choice with respect to environmental conditions that may affect their own development, such as through selection of habitat or diet (e.g. Geiser et al., 1997; Kupferberg, 1997). Whenever they act, the consequences of environmental effects are often termed developmental or phenotypic plasticity.

Fig. 1.

Complex traits, such as behavior, are composed of numerous lower-level(subordinate) traits, themselves interrelated in a strongly hierarchical fashion. In general, natural and sexual selection will tend to act more strongly at higher levels of biological organization, as indicated by the relative thickness of the black arrows. As typically viewed by organismal and evolutionary biologists, selection acts on phenotypic variation (which reflects variation in gene expression), but does not generally act directly on genetic variation (e.g. at the level of DNA sequences). Exceptions to this point can occur via such phenomena as genomic conflict (e.g. Stearns and Hoekstra,2005).

Fig. 1.

Complex traits, such as behavior, are composed of numerous lower-level(subordinate) traits, themselves interrelated in a strongly hierarchical fashion. In general, natural and sexual selection will tend to act more strongly at higher levels of biological organization, as indicated by the relative thickness of the black arrows. As typically viewed by organismal and evolutionary biologists, selection acts on phenotypic variation (which reflects variation in gene expression), but does not generally act directly on genetic variation (e.g. at the level of DNA sequences). Exceptions to this point can occur via such phenomena as genomic conflict (e.g. Stearns and Hoekstra,2005).

Directional natural selection is predicted to have various effects, some fairly obvious but others less so. If natural selection in a given population favors individuals with higher values of a particular trait, then the population mean value of that trait is predicted to increase from generation to generation (Fig. 2A),assuming that some additive genetic variance exists and that the trait is not too strongly genetically correlated with other traits under selection. Evolutionary biology is replete with empirical examples illustrating the validity of this prediction (Endler,1986). Moving from phenotype to genotype, a second prediction is that alleles with appropriate' pleiotropic effects will be favored, which would facilitate the coordinated evolution of components of complex phenotypes. For example, if selection were to favor individuals that foraged widely to find food, then alleles that increased motivation for high locomotor activity might be favored most directly, and the subset of those alleles that also tended to increase ability for high activity would be particularly favored. This sort of process, in which the genetic architecture of the traits involved (especially the additive genetic variance–covariance matrix)evolves to become more consistent with the prevailing pattern of multivariate selection, could facilitate further evolution and adaptive radiation (e.g. Garland, 1994; Schluter, 1996). A somewhat more subtle genetic prediction is that directional selection should tend to favor alleles that exhibit phenotypic dominance in the direction of selection,and this has also received empirical support (e.g. Broadhurst and Jinks, 1974; Henderson, 1981; Hewitt et al., 1981; Mather and Jinks, 1982; Falconer, 1989; Garland et al., 1990; Lynch, 1994; Lynch and Walsh, 1998). For example, under a selective regime that favored high activity levels, alleles that promoted high activity and were dominant to alleles with neutral or negative effects on activity would be the most favored among the spectrum ofhigh-activity alleles'.

A fourth hypothesis, not mutually exclusive with the previous three, is that the average plasticity of the population should also evolve if the selective agent imposes more than instantaneous stress' (sensuHarshman et al., 1999; Wilson and Franklin, 2002; Gabriel, 2005) on the population. More specifically, plasticity of the trait under selection or of subordinate traits that contribute to that trait should increase in the direction that would be adaptive (tend to confer higher organismal performance and/or higher Darwinian fitness) under the prevailing selective regime(Fig. 2B). To clarify this hypothesis, we need to distinguish between selective agents that (1) impose selective events that are virtually instantaneous relative to the time course of possible plastic responses and (2) impose relatively prolonged selective events. As an example of an instantaneous' selective event, a cheetah(selective agent) chases an antelope (selective event) and the outcome (life or death for the antelope) occurs so quickly (a matter of seconds) that the exercise physiology of the antelope has no chance to train'. In addition to its own behavior, the motivation and abilities of the cheetah, and some element of chance (e.g. tripping over a rock), what determines the outcome of the selective event are the innate' (constitutive or intrinsic) exercise abilities of the antelope at the instant the cheetah began its pursuit. As an example of a prolonged' selective event, global warming trends that occur over years would allow long-lived organisms to acclimatize to the higher temperatures in a way that might be beneficial (e.g. upregulation of heat-shock proteins). Thus, the Darwinian fitness of individuals in a crocodilian population experiencing environmental warming might depend on both innate and induced components of heat tolerance. Of course, the distinction between instantaneous and prolonged selective events is not always clean. For example, an antelope might learn from a failed close encounter with a cheetah and became better at avoiding capture in subsequent encounters by altering its own behavior. The classic selective event of an ice storm that killed some(but not all) sparrows in a fortuitously collected sample(Bumpus, 1899; Lande and Arnold, 1983; Endler, 1986) lasted overnight, such that individual sparrows may have varied in thermal tolerance because of both innate and induced individual variation in components of heat tolerance.

Fig. 2.

Hypothetical example of the effects of positive directional selection favoring individuals with higher values for a particular trait on the mean value of that trait (A) and on the plasticity of that trait or of a subordinate trait (B). (A) The standard expectation for the effects of positive directional selection on the distribution of a trait (for example,heat tolerance) across several generations. During generation one, a selective event – high temperature lasting for several days – kills a majority of the individuals in the population (G1) before they can breed. The survivors (S1) of this selective event then breed and the mean heat tolerance in their offspring (G2) is somewhat higher than for their parents (G1). The difference in population mean phenotype between generations one (G1) and two (G2)indicates that evolution has occurred (assuming that the environment in which the organisms are living has not changed in a way that causes the altered phenotypes via direct environmental effects). This process continues for several generations such that the mean value of the trait in generation five (G5) is substantially higher than in generation one. (B) A hypothesis regarding the correlated evolution of the plasticity of heat tolerance or of a subordinate trait that supports heat tolerance (e.g. expression of heat shock proteins). In the original population, exposure to high temperatures for a few hours or days causes some individuals to increase in heat tolerance (which would probably be adaptive if the high temperatures continued) while an equal number of other individuals actually exhibit a decrease in heat tolerance, which would be maladaptive (inappropriate) if high temperatures persisted. For the population as a whole, the average plastic response is zero. Following a selective event and subsequent breeding of the survivors (S1), which produces the next generation (G2),the average plastic response in this new generation tends to be an increase in heat tolerance. Thus, natural selection has caused an evolutionary increase in both the average innate' (or constitutive' or intrinsic') heat tolerance(A) and a shift in the average plasticity of individuals (B) such that, on average, they become more heat tolerant when exposed (acutely) to high temperatures. This constitutes the evolution of adaptive plasticity. See text for discussion of possible genetic mechanisms of such a correlated response to selection.

Fig. 2.

Hypothetical example of the effects of positive directional selection favoring individuals with higher values for a particular trait on the mean value of that trait (A) and on the plasticity of that trait or of a subordinate trait (B). (A) The standard expectation for the effects of positive directional selection on the distribution of a trait (for example,heat tolerance) across several generations. During generation one, a selective event – high temperature lasting for several days – kills a majority of the individuals in the population (G1) before they can breed. The survivors (S1) of this selective event then breed and the mean heat tolerance in their offspring (G2) is somewhat higher than for their parents (G1). The difference in population mean phenotype between generations one (G1) and two (G2)indicates that evolution has occurred (assuming that the environment in which the organisms are living has not changed in a way that causes the altered phenotypes via direct environmental effects). This process continues for several generations such that the mean value of the trait in generation five (G5) is substantially higher than in generation one. (B) A hypothesis regarding the correlated evolution of the plasticity of heat tolerance or of a subordinate trait that supports heat tolerance (e.g. expression of heat shock proteins). In the original population, exposure to high temperatures for a few hours or days causes some individuals to increase in heat tolerance (which would probably be adaptive if the high temperatures continued) while an equal number of other individuals actually exhibit a decrease in heat tolerance, which would be maladaptive (inappropriate) if high temperatures persisted. For the population as a whole, the average plastic response is zero. Following a selective event and subsequent breeding of the survivors (S1), which produces the next generation (G2),the average plastic response in this new generation tends to be an increase in heat tolerance. Thus, natural selection has caused an evolutionary increase in both the average innate' (or constitutive' or intrinsic') heat tolerance(A) and a shift in the average plasticity of individuals (B) such that, on average, they become more heat tolerant when exposed (acutely) to high temperatures. This constitutes the evolution of adaptive plasticity. See text for discussion of possible genetic mechanisms of such a correlated response to selection.

The actual mechanisms for the evolution of increased plasticity could be several, of which we will mention two. First, referring to Fig. 2A, in addition to individuals whose phenotypes are intrinsically high, individuals that exhibit plasticity in the direction favored by selection will tend to be among thesurvivors' each generation; thus, appropriately plastic individuals will be favored by phenotypic selection [see p. 67(Falconer, 1990)]. If plasticity is heritable, then it will evolve in response to such selection(assuming genetic correlations with other traits under selection are not too strong). Second, the genes that affect the constitutive value of the phenotype might also have pleiotropic effects (in the appropriate direction) on the plastic response of that phenotype when the organism experiences chronic (more than instantaneous) exposure to the selective agent (e.g. gradual warming,repeated encounters with predators). The genetics of plasticity are discussed further elsewhere (Scheiner,1993; Pigliucci,2005). Here it is also worth noting that the evolution of plasticity is related to the concept of genetic assimilation, a process in which environmentally induced phenotypic variation that is favored by selection (natural or artificial) gradually (across many generations) comes to be constitutively produced [reviewed elsewhere in this issue(Pigliucci et al., 2006)].

The foregoing ideas about evolutionary processes would seem to be implied by the beneficial acclimation hypothesis' (see below), but to our knowledge they have not been discussed so explicitly in the literature. In any case, we hypothesize that the mean plasticity of a population under directional selection (Fig. 2A) should evolve from being neutral (or possibly deleterious) to being beneficial or adaptive (Fig. 2B). This evolutionary hypothesis has received little direct empirical attention (but see Falconer, 1990; Scheiner, 2002). Nonetheless,as discussed below, many laboratory experiments impose selective regimes that could allow for the evolution of enhanced plasticity (e.g. desiccation trials with Drosophila that last for several days), and the few that have tested for evolutionary changes in plasticity have found some evidence for it.

When behavioral performance traits are the subject of directional selection, the role of phenotypic plasticity in evolutionary response may be particularly interesting [other perspectives on the importance of behavior have been discussed elsewhere (Huey et al., 2003; Price et al.,2003; Price,2006)]. The term self-induced adaptive plasticity' was proposed(Swallow et al., 2005) for situations in which a behavior induces plastic changes in morphological or physiological traits that in turn enhance the ability to perform the behavior. For example, animals that migrate altitudinally might make trial runs' that would cause cardiovascular, pulmonary or metabolic changes that would improve their ability to function at high altitude. Similarly, animals that begin feeding on a new type of food may experience changes in digestive enzymes that increase efficiency of nutrient extraction and/or detoxification [examples of related effects of diet are reported elsewhere, including references therein(Geiser et al., 1997; Kupferberg, 1997)].

The first purpose of this paper is to provide a brief introduction to phenotypic plasticity from an ecological and evolutionary perspective [see also elsewhere in this issue (Fordyce,2006; Pigliucci et al.,2006; Price,2006)]. Second, we discuss how the evolution of plasticity can be studied, with an emphasis on the experimental evolution approach. Finally, we review some results from a study on the experimental evolution of high voluntary activity levels in house mice, including examples of self-induced adaptive plasticity.

From the perspective of evolutionary biology, classic and dramatic examples of phenotypic plasticity in animals include wing polymorphisms in some insects, the timing of metamorphosis in amphibians, and alternative reproductive tactics in male vertebrates – all of which exhibit complex neuro-endocrine control mechanisms that are sensitive to various environmental factors (Ketterson and Nolan, Jr,1999; Sinervo and Calsbeek,2003; Boorse and Denver,2004; Knapp, 2004; Zera, 2004). From the biomedical perspective, well-known examples of plasticity include effects of intentional physical conditioning (exercise training)(Flück, 2006) such as weight lifting, on human morphology and physiology. Various biomedical subfields use additional terminology, such as metabolic plasticity' orcardiac remodeling', and the molecular mechanisms underlying such processes as muscular and neuronal plasticity are the subject of intensive study [for reviews, see other articles in this issue(Flück, 2006; Hood et al., 2006; Johnston, 2006; Magistretti, 2006; Swynghedauw, 2006)]. (Many environmental insults, e.g. excessive alcohol consumption, smoking, inhalation of coal dust, can lead to plastic' changes in organs and organ systems, but when such changes are clearly pathological they are not typically included under the rubric of phenotypic plasticity.) In plants, basic growth form is notoriously plastic, and many readers will be familiar with the differences between dandelions growing in shade versus sun [although genetic differences among clones may also be involved(Collier and Rogstad,2004)].

As with the term adaptation' (see below), phenotypic plasticity can refer both to a process and to the outcome of that process. Phenotypic plasticity can be defined formally as the ability of one genotype to produce more than one phenotype when exposed to different environments, as the modification of developmental events by the environment, or as the ability of an individual organism to alter its phenotype in response to changes in environmental conditions (Gordon, 1992; Scheiner, 1993; Via et al., 1995; Futuyma, 1998; Freeman and Herron, 2004; Pigliucci, 2005; Rezende et al., 2005; Stearns and Hoekstra, 2005; Pigliucci et al., 2006). The range of phenotypes that a given genotype (possessed by an individual organism or by an entire clone or inbred line) may produce when exposed to a range of environmental conditions is termed its norm of reaction, and non-parallel reaction norms of different genotypes indicate the presence of genotype-by-environment interaction.

The sequence of events involved in phenotypic plasticity often includes the following components: (1) something in the environment changes; (2) the organism senses that change; (3) the organism alters gene expression; and (4),usually, the altered gene expression yields additional observable phenotypes[e.g. see fig. 8 in Flück's paper in this issue(Flück, 2006)]. Several aspects of this scenario require amplification. With respect to (1), we may attempt to draw a distinction between environmental factors that are external or internal to an organism. Changes in ambient temperature, humidity or oxygen concentration would constitute external environmental factors, and many organisms respond to these with phenotypic plasticity that involves multiple organ systems and multiple levels of biological organization. Mechanical overload of the heart is an example of an environmental change that occurs within an organism, and it leads mainly to organ-specific changes that necessarily involve fewer levels of biological organization(Swynghedauw, 2006). Of course, external environmental stresses' can also lead to tissue-specific responses (e.g. Cossins et al.,2006). Nonetheless, we may predict that, in general, external environmental changes will lead to more and more pervasive plastic responses as compared with internal changes. With respect to (2), some changes may occur without any formal sensing by the organism, e.g. as a result of direct (and possibly differential) effects of temperature on the rates of ongoing biochemical and physiological processes. With respect to (3), it is important to note that some plastic responses need not involve changes in gene expression (transcription) but instead could occur viaphosphorylation of existing proteins, changes in protein levels caused by variation in protein ubiquitination, or stimulation of existing microRNAs(Nelson et al., 2003; Schratt et al., 2006). For point (4) we emphasize the word usually' because it is possible that lower level traits might change in offsetting ways such that higher level traits could show little or no apparent change. For example, it would be theoretically possible (though perhaps unlikely) for exercise training to cause an increase in maximal heart rate but a reduction in stroke volume such that cardiac output was unchanged.

Acclimation and acclimatization(Wilson and Franklin, 2002),as well as learning and memory (e.g. Magistretti, 2006), are encompassed by the most inclusive definitions of phenotypic plasticity. Therefore, environmentally induced changes may or may not be reversible,depending on the organism, trait, and when in the lifecycle and for how long the environmental exposure occurs (Hatle,2004; Johnston,2006). If the capacity for change is more-or-less fully reversible, then it may be termed phenotypic flexibility(Piersma and Lindstrom,1997).

Whether reversible or not, it is generally assumed that environmentally induced modifications are adaptive in the sense that they improve organismal function and/or enhance Darwinian fitness of the individual organisms that exhibit such effects (Nunney and Cheung,1997). In fact, this may or may not be true, and the claim that such changes will aid the organism has been termed the beneficial acclimation hypothesis (Leroi et al.,1994; Huey and Berrigan,1996; Huey et al.,1999; Wilson and Franklin,2002). In some cases, behavioral plasticity can shield lower level traits from selection (Huey et al.,2003; Price et al.,2003). At the population level, phenotypic plasticity in behavior and other traits can facilitate invasions of new habitats(Price et al., 2003; Price, 2006; Pigliucci et al., 2006). As reviewed elsewhere in this issue (Fordyce,2006), many ecological (cross-species) interactions are mediated by the phenotypically plastic responses of one or more species involved in the interaction. Some of these ecological interactions can be quite complex and difficult to predict, as when an herbivore induces a plant phenotype that in turn affects the performance of other herbivores(Fordyce, 2006)

At this point it is worth remembering that the word adaptation has numerous meanings in biology (Garland and Carter,1994; Bennett,1997). Most generally, we should keep in mind the distinction between what is often called physiological adaptation' (environmentally induced changes that occur within individual organisms during their lifetimes,including acclimation and acclimatization) and evolutionary adaptation'(cross-generational changes in the genetic composition of a population in response to natural selection). Physiological adaptation is one type of phenotypic plasticity, but the ability to be plastic for any particular trait may also be an evolutionary adaptation whose details vary among organisms.

As noted above, although biologists have usually assumed that physiological adaptation is adaptive in the evolutionary sense, this is not always a safe assumption because some changes will be simply the result of activation of control systems designed to do something else, and they can even be maladaptive, including various human pathologies(Nesse, 2005; Swynghedauw, 2006). In general, non-adaptive plasticity might be expected to occur any time that an organism is exposed to environmental conditions with which it is unfamiliar'in terms of its evolutionary history. This follows from the general evolutionary principle that organisms gradually lose abilities and traits that are no longer under positive selection, well-illustrated by things like blind cave fish or flightless birds on islands that lack predators(Diamond, 1986). Thus, imagine a species that has inhabited low-elevation environments for millions of years,adapting evolutionarily to function (reasonably) well in normal' levels of atmospheric oxygen (∼21%). If one were to expose individuals of this species to high altitude, then they might be expected to exhibit inappropriate physiological responses to reduced atmospheric oxygen. The literature on human physiological responses to high altitude, both acute and chronic, is interesting in this context because it offers conflicting views on whether and to what extent various changes are adaptive versus maladaptive, and whether long-term, high-altitude native populations exhibit evolutionary adaptations to hypoxia (e.g. Winslow et al., 1989; Beall,2001; Brutsaert et al.,2005; Norcliffe et al.,2005; Wu et al.,2005). More generally, it is worth noting that the environment that many human beings experience (including aspects of nutrition, sanitation,medicine and the so-called built environment) has changed very rapidly relative to our generation time. Concomitantly, average lifespan has increased in many countries and diseases associated with aging have become much more common (e.g. Swynghedauw,2006). Therefore, it may be expected that at least some aspects of our phenotypic plasticity may not be adaptive.

Intuitively, plasticity might be good or bad, depending on the amount of spatial heterogeneity in the environment, the speed of temporal environmental changes, the predictability of spatial and temporal heterogeneity, and the size or duration of heterogeneity relative to an organism's mobility and lifespan. From a formal theoretical perspective, the evolution of plasticity has been studied with optimality models, quantitative genetic models, and gametic models (Scheiner,1993). Generally, all of these models suggest that adaptive plasticity will evolve when environmental heterogeneity exists, environmental cues about that heterogeneity are somewhat reliable, plastic responses confer a net fitness benefit, and the population contains some additive genetic variance for the plastic response (Berrigan and Scheiner, 2004). With regard to spatial variability,optimality, quantitative genetic and gametic models all predict further that plasticity is most favored when (1) inter-habitat variability is high, (2) all habitats are equally regular, (3) selection acts equally strongly across habitats, (4) the environmental cue-dependent phenotype is correlated with the environment of selection, (5) habitat selection is correlated with trait plasticity [for specific references, see elsewhere(Scheiner, 1993)].

Phenotypic plasticity is typically induced by environmental heterogeneity or environmental stress (Harshman et al.,1999; Wilson and Franklin,2002; Berrigan and Scheiner,2004; Gabriel,2005). In this context, stress' is generally taken to mean anything that threatens physiological homeostasis (e.g. Sapolsky, 2005) and/or reduces Darwinian fitness. Environmental stress can be categorized into biotic(e.g. predator presence) versus abiotic (e.g. ambient temperature),and either type may cause changes in behavior, morphology, and/or physiology(Gabriel, 2005). If the mean fitness of individuals with plastic strategies exceeds the mean fitness of those with fixed strategies, then phenotypic plasticity or flexibility will tend to evolve (Scheiner,1993; Berrigan and Scheiner,2004; Gabriel,2005). Environment tolerance curves have been defined as the response of a genotype's total fitness over an environmental gradient'(Lynch and Gabriel, 1987),distinguishing this as a special case of the norm of reaction, and using them to predict when irreversible plasticity will tend to evolve.

However, as noted elsewhere, If stress periods are short compared to the life-time of an organism, then irreversible phenotypic plasticity is unlikely to be a favorable response'(Gabriel, 2005). Therefore,Gabriel proposed models predicting the selective advantage of reversible plasticity (phenotypic flexibility)(Gabriel, 1999; Gabriel, 2005). He concluded(Gabriel, 2005) that ... reversible phenotypic plasticity would be expected for all organisms under the following conditions: they are exposed to stress periods that last shorter than life span; stress appears in the long run with some regularity so that natural selection can shape non-induced and induced values of adaptive plastic traits.' In these models, he assumed that plasticity was not costly, with the rationale that Plasticity costs would usually enter as constant factors that do not alter the optimal values of mode and breadth' [see p. 875 (Gabriel, 2005)]. He added the caveat that if plasticity costs depend significantly on the amount of performed phenotypic change, then costs might become a function of the environmental state during stress in a way that the optimal values of mode and breadth are affected' (p. 875), but concluded by arguing that given the predicted huge fitness advantages, the cost of plasticity would have to be unexpectedly high in order to counteract selection for reversible phenotypic plasticity' (pp. 880–881). Thus, it is important to remember Pigliucci's point on p. 483 (Pigliucci,2005) that Research of costs of plasticity is still in its infancy, but is both theoretically important and empirically challenging, and should become a major area of future inquiry.'

As discussed above, natural selection ought to affect plasticity, and organisms ought to vary in plasticity. How can we test such theoretical predictions? In general, the same way that we may seek to study adaptation in any sort of trait. Four general approaches to studying adaptation are commonly used by evolutionary biologists (e.g. see Huey and Kingsolver, 1993; Garland and Carter, 1994; Bennett, 1997; Futuyma, 1998; Schlichting and Pigliucci,1998; Feder et al.,2000; Orzack and Sober,2001; Pigliucci,2001; Freeman and Herron,2004; Stearns and Hoekstra,2005). First, as outlined in the previous section, real organisms can be compared with predictions of theoretical models, such as those based on optimality (e.g. Garland,1998; Orzack and Sober,2001). Second, examinations of the biology of natural populations can determine what sorts of traits vary, are heritable, and are currently under sexual or natural selection (e.g. Young et al., 2004). Experimental manipulations of putatively adaptive traits are often employed in such studies (e.g. Sinervo and Basolo,1996; Ketterson and Nolan, Jr,1999) and pp. 224-229 (Costa and Sinervo, 2004). Although several studies have attempted to quantify how natural selection acts on plasticity in the field (e.g. Trussell, 1997; Donohue et al., 2000; Nussey et al., 2005), this area of investigation will not be covered here. Third, one can compare species(or populations) that vary with respect to ecological factors that might cause variation in how selection views' plasticity [overviews of studying adaptation via the comparative method' and with a phylogenetic perspective have been published previously(Garland and Adolph, 1994; Garland et al., 2005)]. In the following subsection, we provide a brief summary of some comparative studies of plasticity. Finally, selection experiments(Bennett, 2003; Garland, 2003; Swallow and Garland, 2005)can be used to study adaptation, and this is our main focus, with emphasis on those that would qualify as experimental evolution' (e.g. Rose et al., 1996; Rose et al., 2004; Ebert, 1998; Bennett, 2002; Bennett, 2003; Swallow and Garland, 2005)(http://en.wikipedia.org/wiki/Experimental_evolution).

### Comparative studies

Vertebrate morphology and physiology provide dramatic examples of both inter-specific and inter-trait variation in plasticity [plasticity of the water barrier in vertebrate integument is reviewed elsewhere(Lillywhite, 2004)]. With respect to variation among traits, vertebrate skeletal muscle(Flück, 2006) and the gastrointestinal tract (Secor,2005) are very responsive to use and disuse (training' anddetraining' effects). Bone size, shape and architecture also change in response to variation in loading conditions, but to a much smaller extent than for muscle [for example, compare(Houle-Leroy et al., 2000)with (Kelly et al., 2006)]. Adult vertebrate lung also seems to have relatively low plasticity (e.g. Hoppeler et al., 1995; Weibel, 2000; Hsia, 2001; Henderson et al., 2002). In plants, one study shows that aspects of gas exchange may be more plastic than structural traits (Valladares et al.,2000).

Interspecific variation in plasticity has also been documented. In vertebrates, for example, attempts at aerobic exercise training (to improve cardiopulmonary and/or muscular function) of lizards have generally not been successful, even when patterned after those that cause large changes in mammals (Garland and Else,1987; Conley et al.,1995) (A. Szucsik, personal communication). In amphibians and squamates, species differences in gut plasticity seem to be related to their feeding ecology, in particular the frequency and/or regularity of feeding(Secor, 2005). In Burmese pythons, ventricular mass can increase 40% within 48 h after feeding, a change that is fully reversible (Andersen et al.,2005). Among species of fishes, carp and goldfish seem to be especially plastic (Cossins et al.,2006; Johnston,2006). In plants, a common-garden study of 16 shrubs in the genus Psychotria showed that species found in the understory, where light is less variable, showed less plasticity for a variety of traits as compared with species that generally occur in forest gaps, where light is more variable(Valladares et al., 2000). Population differences in plasticity have also received considerable attention in plants, with several studies suggesting that the they are indeed adaptive(e.g. Cook and Johnson, 1968; Donohue et al., 2000).

### Selection experiments and experimental evolution

Selection experiments have provided valuable insights into central questions surrounding the evolution of phenotypic plasticity (see Scheiner, 2002). At their most basic, they have demonstrated that the plasticity of a trait is often heritable, capable of responding rapidly to selection, and determined by multiple genetic loci. [A genetic basis for the response to aerobic exercise training has also been demonstrated in human twin studies (references in Koch et al., 2005).] In addition, selection experiments have shown that plasticity (environmental sensitivity) of a given trait can evolve independently of the population mean value for that trait. More specifically, experiments with plants and invertebrates have shown plasticity to evolve in response to selection (1)directly on the reaction norm, (2) on a single trait in one environment, and(3) on a single trait across multiple environments. With respect to vertebrates, although many selection experiments have been performed, very few have focused on phenotypic plasticity as a component of the response to selection (Falconer, 1990; Scheiner, 2002).

The reaction norm has been directly selected upon in Drosophila melanogaster, the butterfly Bicyclus anynana, and the tobacco plant Nicotina rustica(Scheiner, 2002). We will only highlight the experiment performed by Scheiner and Lyman(Scheiner and Lyman, 1991),as it appears to be the most comprehensive and has also been reviewed in detail (Scheiner, 2002).

The stated purpose in Scheiner and Lyman's experiment(Scheiner and Lyman, 1991)was to determine if plasticity could respond to selection that was imposed under controlled and reproducible conditions. They began by capturing 301 individual D. melanogaster from the wild. These flies founded a stock that was maintained in the laboratory by mass culture at 21°C for 2–3 months (several generations), thus establishing a genetically heterogeneous base population. They then used 50 randomly chosen pairs to establish each of 14 separate experimental lines, which comprised two replicates of each of six selection regimes (increased thorax size at 19°C, decreased thorax size at 19°C, increased thorax size at 25°C, decreased thorax size at 25°C, increased plasticity, decreased plasticity) plus a control line that was not intentionally selected. To impose selection, plasticity was defined as the difference in average thorax length for sets of full-sibs raised at 19°C and 25°C. Plasticity did indeed respond to selection, but with a rather low realized heritability of 0.088±0.027 (mean ± s.e.m.). The authors concluded that the plasticity was not the result of overdominance, but rather a genetic interaction among multiple loci.

As noted above, the evolution of phenotypic plasticity has also been examined as a correlated response to selection on a specific trait in a single environment. For example, Harshman et al.(Harshman et al., 1991)studied detoxification enzymes in D. melanogaster. After establishing a base population from wild-caught flies, three Control (C) lines were reared on standard medium and three Selected (S) lines on lemon for 20 generations. For the lemon-cultured lines, the selection process was as follows: (1) flies were placed in bottles with freshly cut lemon (10 g, pesticide free) at room temperature for 7–10 days; (2) approximately 50% mortality occurred; (3)survivors were placed into a new bottle of freshly cut lemon (30 g) and vermiculite to produce the next generation. According to Harshman et al., the 50% mortality (during the lemon selection episodes) may have been caused by natural insecticide activity in lemons, or by toxins produced by bacteria growing on the fruit (Harshman et al.,1991). Flies (35–70) were randomly mated to produce the subsequent generation in all six lines. In the Control lines, flies were transferred to fresh medium for mating.

After 20 generations, all flies to be tested were reared on ordinary medium for one generation to standardize environmental conditions. They were then transferred to either lemon (which may induce the expression of detoxification enzymes) or fresh medium (to allow determination of enzyme activities under baseline conditions) for 24 h prior to sacrifice. Activities of epoxide hydrolases and glutathione S-tranferase (GST) were then measured. For GST measured using trans-stilbene oxide (TSO) as a substrate, S and C lines showed no significant difference for the sample exposed to fresh medium for 24 h, but the S lines showed substantially higher enzyme activities than C lines when lemon-exposed for 24 h (Fig. 3). Harshman et al. concluded(Harshman et al., 1991):After 20 generations on lemon there was a pronounced change in environment-dependent expression... The response appeared independently in all three lines on lemon.' They also noted that: In the present study the culturing regime used was ostensibly continuous, unless the process of lemon rotting every generation constitutes temporal variation. Normally, one would anticipate selection for change in environment-dependent enzyme expression to occur in variable environments but the results of the present study suggest it can evolve in a relatively constant regime.' Harshman et al. give additional examples in which the plasticity of an enzyme activity seems to have evolved as a correlated response(Harshman et al., 1999).

Fig. 3.

Example of a selection experiment (20 generations) with Drosophila melanogaster (Harshman et al.,1991) in which plasticity evolved to be higher in the Selected lines (S; N=3) as compared with the Control lines (C; N=3). C flies (left) were exposed to either standard medium or lemon for 24 h prior to sacrifice for measurement of a detoxification enzyme, and S flies (right)were similarly treated. Values are means ± s.d. For C lines, the magnitude of the induction caused by lemon exposure, as indicated by glutathionase S-transferase activity on lemon divided by the value on standard medium, was 1.16, whereas for S lines the value was 2.58. The greater induction caused by lemon exposure in the S lines relative to the C lines is an example of a genotype-by-environment interaction. TSO, trans-stilbene oxide. See text for further details.

Fig. 3.

Example of a selection experiment (20 generations) with Drosophila melanogaster (Harshman et al.,1991) in which plasticity evolved to be higher in the Selected lines (S; N=3) as compared with the Control lines (C; N=3). C flies (left) were exposed to either standard medium or lemon for 24 h prior to sacrifice for measurement of a detoxification enzyme, and S flies (right)were similarly treated. Values are means ± s.d. For C lines, the magnitude of the induction caused by lemon exposure, as indicated by glutathionase S-transferase activity on lemon divided by the value on standard medium, was 1.16, whereas for S lines the value was 2.58. The greater induction caused by lemon exposure in the S lines relative to the C lines is an example of a genotype-by-environment interaction. TSO, trans-stilbene oxide. See text for further details.

In numerous other selection experiments where plasticity could potentially evolve (because the selection regime is more than instantaneous), plasticity of the selected trait or of potentially related or subordinate traits does not appear to have been examined. In one such example(Bubliy and Loeschcke, 2005),correlated responses to selection for stress resistance and longevity in a laboratory population of D. melanogaster were examined. Several selection regimes were imposed in this experiment: cold-shock resistance selection, heat-shock resistance selection, heat knockdown resistance selection, desiccation resistance selection, starvation resistance selection and longevity selection. Here we discuss the cold-shock resistance selection line, for which selection was imposed for 21 generations. The selection regime was as follows. Flies were maintained on standard medium for 5 days at 11°C for acclimation, then placed in empty vials and exposed to 0.5°C for 27–50 h with relative humidity near 100%. Surviving flies were allowed to recover for 24 h at 25°C in vials with standard medium, then allowed to reproduce. We would argue that plasticity of traits that may support cold resistance in D. melanogaster could potentially be altered during the acclimation phase of this selection protocol. Furthermore,plasticity of traits that may support cold resistance may be evolving across generations. Although phenotypic plasticity apparently has not been investigated in these flies, we contend that selection experiments of this type should explore whether plasticity has increased as a component of the response to selection. In the following section, we discuss an ongoing experiment with house mice that has begun to examine plasticity in various traits as a potential correlated response to artificial selection for high voluntary wheel running, as expressed during days 5 and 6 of a 6-day exposure to wheels.

Since 1993 our laboratory has been conducting a replicated selection experiment for high voluntary wheel-running behavior on days 5+6 of a 6-day wheel exposure. By housing mice from each of the four replicate S lines and from each of the four replicate C lines with or without wheel access for several days or weeks, we can test for differences in plasticity (training effects) in various traits. As outlined in the remainder of this section, we have found several traits that show greater differences between S and C lines when they are housed with wheel access than when they are housed without wheel access (or, in some cases, housed with access to wheels that are locked to prevent rotation). For some of these traits, the greater differences can be explained statistically by the greater wheel running exhibited by mice from S lines. For others, however, the differences seem to reflect greater plasticity in the S lines [i.e. for a given amount of stimulus (wheel running/day),individuals in the S lines show a greater response than in the C lines].

### Animals and experimental protocol

The original progenitors (founding population) were outbred, genetically variable house mice (Mus domesticus) of the Hsd:ICR strain(Harlan-Sprague-Dawley, Indianapolis, IN, USA). After purchase from HSD, mice were randomly mated for two generations, paired, and then assigned randomly to eight closed lines (10 pairs in each). Four of these lines have been designated to experience selective breeding for high voluntary activity (lab designation, lines 3, 6, 7, 8) and four serve as controls (lines 1, 2, 4,5).

The selection protocol has been described in detail elsewhere(Swallow et al., 1998a). In brief, when each generation of mice are 6–8 weeks old, they are housed individually with access to running wheels (circumference=1.12 m) for 6 days. Daily wheel-running activity is monitored with photocell counters linked to a computer-automated system. Wheel activity is recorded in 1-min bins for 23–24 h of each of the 6 days of wheel access. For purposes of selection, wheel running is quantified as the total number of revolutions on days 5 and 6 of the 6-day test. After accounting statistically for any variation related to measurement block, age, wheel resistance, and sex,breeders are chosen. In the four S lines, the highest running male and female are chosen from each family as breeders to propagate the lines of the next generation. Within-family selection is performed to increase the effective population size (Ne), while reducing maternal and environmental variances, including effects of genotype-environment interactions (Henderson,1989). In the four C lines, breeders are randomly chosen from each family. Within all lines, sibling matings are disallowed.

By generation 16, the high-activity lines exhibited a 170% increase in total revolutions/day as compared with the C lines. This was caused primarily by S mice running faster rather than for more minutes each day, but the relative importance of the two components differs between the sexes, with females from the S lines typically showing little or no increase in amount of time running whereas males do show an increase in time running(Swallow et al., 1998a; Koteja et al., 1999a; Koteja et al., 1999b; Rhodes et al., 2000; Girard et al., 2001). This increase in wheel running greatly exceeds that of wild house mice born and raised under the same conditions (Dohm et al., 1994), and comes close to spanning the range of variation that has been reported among 13 species of wild murid rodents(Garland, 2003). Therefore, it seems that we have an evolutionarily important' amount of divergence in wheel running between the S and C mice. Additionally, based on high-speed video analyses, estimates of instantaneous running speeds have shown that S line females run twice as fast as C line females, as well as more intermittently(Girard et al., 2001). However, since approximately generation 16, the differential in wheel-running distances has remained relatively constant, indicating that a selection limit or plateau may have been attained.

### Plasticity of wheel running

Because our wheel-testing protocol is prolonged (6 days) rather than instantaneous (e.g. a few minutes), it is possible that the S lines may have evolved greater plasticity in this behavior. In other words, as compared with the C lines, mice from S lines might now exhibit a greater increase in wheel running across the 6 days of wheel access, given that only their performance on days 5+6 affects their probability of reproducing. Fig. 4 shows some hypothetical examples of how plasticity might be greater in the S lines. In Fig. 4A, Control lines show constancy of wheel running, whereas S lines increase monotonically across days 1–6; clearly, plasticity is greater in S lines. In Fig. 4B, wheel running is identical and increases gradually in both S and C lines over the first four days of testing. Selected lines then show a much greater increase between days 4 and 5, thus indicating greater plasticity during this time period. In these two cases, the greater plasticity of S as compared with C lines is reflected in the ratio of S/C (see Fig. 4A,B).

Fig. 4.

Hypothetical trajectories for the amount of voluntary wheel running across a 6-day test period, as implemented in the selective breeding experiment with house mice. (A,B) Examples of greater plasticity in the Selected lines (S;black circles) than in the Control lines (C; gray squares). (C) Similar relative increases in wheel running on a day-to-day basis, but greater absolute increases in the S lines. See text for discussion.

Fig. 4.

Hypothetical trajectories for the amount of voluntary wheel running across a 6-day test period, as implemented in the selective breeding experiment with house mice. (A,B) Examples of greater plasticity in the Selected lines (S;black circles) than in the Control lines (C; gray squares). (C) Similar relative increases in wheel running on a day-to-day basis, but greater absolute increases in the S lines. See text for discussion.

In the case of Fig. 4C, the interpretation is more complicated. Both S and C lines increase monotonically across days 1–6. On an absolute basis, S lines increase more (2000) than C lines (667) on each day. Relative to their own starting values, however, S and C lines increase by the same percentage each day, although this increase becomes smaller each day (40, 29, 22, 18 and 15%, respectively). As a result,the S/C ratio is a constant. Thus, whether one considers the S and C lines to differ in plasticity depends on whether absolute or relative values are considered.

Fig. 5 shows example data from our selection experiment, and the pattern resembles the one shown in Fig. 4C. As reported elsewhere(Belter et al., 2004), 48 female mice from generation 23 were studied. As shown in Table 1, S lines ran significantly more than C on every day. A repeated-measures ANOVA (SAS Proc Mixed with autoregressive error structure) indicated highly significant effects of day (P<0.0001) and line type (P=0.0001), but no significant day-by-line type interaction (P=0.7184). The foregoing results suggest that S lines do not exhibit a greater plasticity in wheel running.

Table 1.

Total revolutions during 6 days of wheel access for Control and Selected female mice from generation 23(EF3Belter et al.,2004)

Total no. of revolutions
TraitControlSelectedS/C RatioP
Day 1 3071±822 8509±746 2.77 0.0017
Day 2 4965±1234 12128±1153 2.44 0.0041
Day 3 3822±800 11379±678 2.98 0.0001
Day 4 4752±1219 13063±1074 2.75 0.0012
Day 5 5105±994 13277±819 2.60 0.0003
Day 6 5153±1031 14141±878 2.74 0.0002
Day 6—Day 1 2112±1506 5658±1365 2.68 0.1047*
log10 (Day 6/Day 1) 0.2525±0.0723 0.2236±0.0649  0.7541
Total no. of revolutions
TraitControlSelectedS/C RatioP
Day 1 3071±822 8509±746 2.77 0.0017
Day 2 4965±1234 12128±1153 2.44 0.0041
Day 3 3822±800 11379±678 2.98 0.0001
Day 4 4752±1219 13063±1074 2.75 0.0012
Day 5 5105±994 13277±819 2.60 0.0003
Day 6 5153±1031 14141±878 2.74 0.0002
Day 6—Day 1 2112±1506 5658±1365 2.68 0.1047*
log10 (Day 6/Day 1) 0.2525±0.0723 0.2236±0.0649  0.7541

Values are least-squares (adjusted) means ± s.e.m.(N=48).

P-values are for two-tailed tests comparing Control and Selected lines.

Data were analyzed by mixed-model nested ANCOVA in SAS Proc Mixed version 8. The effect of line type (selected vs control) was tested over replicate lines nested within line types (d.f.=1,6), and replicate lines were considered a random factor nested within line type. Age was entered as a covariate and mini-muscle' was entered as an additional factor in all analyses; neither was ever statistically significant (all P>0.05).

*

If one statistical outlier is removed, the line type effect becomes significant (P=0.0325).

On the other hand, the difference between total revolutions on day 6 and day 1 was considerably higher, on average, for S lines (5658) than for C lines(2112). This greater absolute increase in wheel running across 6 days is not statistically significant (P=0.1047), but becomes significant(P=0.0325) when one outlier is removed. This was an S individual whose wheel running declined anomalously from 14 375 on day 1 to 7603 on day 6, the greatest decline for any mouse in the sample of 48. This may represent a real phenomenon, or it might indicate a problem with the wheel on day 6. We intend to explore the plasticity of wheel running more in future studies, with larger sample sizes. In any case, we believe that the greater increase in wheel running across the 6-day trial may well have biological significance,and may well have required coadaptational changes in one or more subordinate traits that support wheel running.

Fig. 5.

Wheel running of 48 female mice from generation 23(Belter et al., 2004), with Selected lines (S) depicted by black circles and Control lines (C) as gray squares. Values are least-squares (adjusted) means from ANOVA, as shown in Table 1. Compare with Fig. 4C.

Fig. 5.

Wheel running of 48 female mice from generation 23(Belter et al., 2004), with Selected lines (S) depicted by black circles and Control lines (C) as gray squares. Values are least-squares (adjusted) means from ANOVA, as shown in Table 1. Compare with Fig. 4C.

### Apparent exercise adaptations in the high-activity lines

A main goal of the selection experiment was to identify traits that have evolved in concert with increased activity levels and that may be necessary for them, i.e. evolutionary adaptations for the high wheel running. Several considerations make this goal non-trivial. First, exercise physiology is complicated, and we have not examined all possible subordinate traits that could be key in terms of allowing high wheel running. Second, of those traits that have been examined, not all have been examined in the same generation. Some adaptations may have occurred in earlier generations and others in later ones, and indeed those occurring in later generations may even have supplanted some that occurred earlier. Third, adaptations may only exist, or at least be more developed, around the age at which wheel testing normally occurs, which is 6–8 weeks of age. Fourth, adaptations may only exist on days 5 and 6 of wheel testing, i.e. they require some days of wheel access to develop. Fifth, adaptations may to some extent be sex-specific, especially given that females in the S lines have increased total activity almost entirely by running faster, whereas males also show an increase in amount of time spent running. Given that we have not studied both sexes, at all ages, under all possible housing conditions (e.g. with or without wheel access), let alone in every generation, we may well have missed some key adaptations. With those cautions in mind, we have discovered a number of traits that seem to represent adaptations for high wheel running in the S lines. We review the motivational basis for high wheel running elsewhere(Rhodes et al., 2005).

Mice from the selected lines have higher maximal oxygen consumption during forced treadmill exercise (VO2max), especially in males (Swallow et al.,1998b; Rezende et al.,2006a; Rezende et al.,2006b) and higher insulin-stimulated glucose uptake in the extensor digitorum longus muscle [located in the hindlimb(Dumke et al., 2001)]. Mice from S lines have larger femoral heads and more symmetrical hindlimb bone lengths (Garland and Freeman,2005; Kelly et al.,2006). Interestingly, S lines exhibit reduced hindlimb muscle mass, especially in two lines that contain a Mendelian recessive allele that halves hindlimb muscle mass while increasing mass-specific aerobic capacity and having a variety of other pleiotropic effects(Garland et al., 2002; Houle-Leroy et al., 2003; Swallow et al., 2005; Syme et al., 2005; Kelly et al., 2006). The S and C lines differ with respect to many other traits as well, such as higher plasma corticosterone levels (Girard and Garland, 2002) and reduced body fat in S lines(Swallow et al., 2001; Dumke et al., 2001). We are currently attempting to determine which of these are adaptations that enhance wheel-running ability, as opposed to non-adaptive (and possibly even maladaptive) correlated responses.

### Plasticity of exercise-related traits

Many traits (e.g. heart mass, VO2max) that one might expect to evolve as a correlated response to selection for high activity levels are also known to respond to the amount of exercise that an individual organism exhibits. Indeed, the literature on mammalian training effects is immense, in part because of our interest in competitive athletics but also because many exercise-related traits are known or thought to be important in promoting physical and/or psychological health(Booth et al., 2002; Castaneda et al., 2005)(Health Activity Center: www.cvm.missouri.edu/hac/index.html). Given that mice from the S lines run more than C when given wheel access, they might also be expected to exhibit greater training responses (physical conditioning) over a given period of time, such as several weeks.

Imagine that groups of both S and C mice were housed either without(sedentary group) or with (active group) access to a running wheel for 8 weeks(e.g. Swallow et al., 2005). Imagine further that for mice housed without wheels, we observed no difference in some phenotype, such as hematocrit. For the mice housed with running wheels, consider a hypothetical phenotype for which values are higher in the S lines (Fig. 6A); this can be explained, statistically at least, by their higher wheel running: a single regression line adequately describes the relation. In this case, we would interpret the data as indicating that mice from S and C lines are equally plastic: it seems to be a simple case of more pain, more gain.' One real example of this pattern involves the level of brain-derived neurotrophic factor (BDNF) in the hippocampus of S and C mice after one week of access to running wheels [see fig. 2 (Johnson et al., 2003)].

Fig. 6.

Hypothetical relations between a phenotypic trait and the amount of running exhibited during the final week of a multi-week exposure to wheels. In both A and B, it is assumed that mice housed without wheels (not shown) would have values of the phenotype lower than or equal to (about 24) those exhibited by Control mice housed with wheels. (A) The greater phenotypic values for the Selected lines (S; black circles) as compared with Control lines (C; gray squares) are explainable statistically by their greater amount of running(more pain, more gain'). A real example of this pattern involves the level of brain-derived neurotrophic factor (BDNF) in the hippocampus of S and C mice after 1 week of access to running wheels [see fig. 2(Johnson et al., 2003)]. (B)There is no relation between the amount of running and phenotype within either group and for a given amount of running the increase in phenotype (relative to the values when mice do not have wheel access) is greater for selected lines than for C lines. Hence, S lines exhibit a greater plastic response. See text for further explanation.

Fig. 6.

Hypothetical relations between a phenotypic trait and the amount of running exhibited during the final week of a multi-week exposure to wheels. In both A and B, it is assumed that mice housed without wheels (not shown) would have values of the phenotype lower than or equal to (about 24) those exhibited by Control mice housed with wheels. (A) The greater phenotypic values for the Selected lines (S; black circles) as compared with Control lines (C; gray squares) are explainable statistically by their greater amount of running(more pain, more gain'). A real example of this pattern involves the level of brain-derived neurotrophic factor (BDNF) in the hippocampus of S and C mice after 1 week of access to running wheels [see fig. 2(Johnson et al., 2003)]. (B)There is no relation between the amount of running and phenotype within either group and for a given amount of running the increase in phenotype (relative to the values when mice do not have wheel access) is greater for selected lines than for C lines. Hence, S lines exhibit a greater plastic response. See text for further explanation.

Fig. 6B shows a different situation. When housed with wheel access, mice from S lines again have higher values for the phenotype, but we see no relation with the amount of running within either group. If we imagine further that S and C mice housed without wheels showed no difference (or at least values similar to those of C mice housed with wheels), then the S mice seem to be more responsive to wheel exposure, i.e. they are more plastic. For a given amount of exercise (wheel running), S mice experience a greater training response. Remember also that phenotypic differences between genotypes (e.g. S versus C mice) that appear only in some environments are termed genotype-by-environment interactions.

As in the hypothetical scenarios just discussed, we have published several papers that involved groups of both S and C mice housed with or without access to functional wheels for several weeks. We have studied various traits,including body mass, VO2max, organ masses, bone properties and enzyme activities (Swallow et al., 1999; Houle-Leroy et al., 2000; Thomson et al.,2002; Belter et al.,2004; Swallow et al.,2005; Kelly et al.,2006). We found a variety of responses in these phenotypes,including some that differ between the sexes. Some traits do not differ between S and C mice regardless of housing conditions [e.g. tail length,adjusted for variation in body mass, in both sexes(Swallow et al., 2005)]. Some traits were found to differ between S and C mice regardless of housing conditions [e.g. S mice are smaller in body mass but have relatively larger kidneys (Swallow et al.,2001; Swallow et al.,2005; Kelly et al.,2006)]. Others showed a difference between S and C lines when housed with wheels but not when housed without [e.g. hematocrit and blood hemoglobin content (Swallow et al.,2005)].

For traits that differ more between S and C lines when they are housed with wheel access, we can examine statistically which of the competing patterns shown in Fig. 6A,B better describes the data (see also Swallow et al., 2005). The general strategy is as follows. First, we identify a trait that shows a statistical interaction between the effects of line type(S versus C lines) and wheel access (the environmental factor). For these analyses, we use SAS Proc Mixed to implement a mixed-model, nested ANOVA(or ANCOVA if such covariates as age or body mass are included in the model),in which replicate line is a random effect nested within line type (S or C). Degrees of freedom for testing the effect of line type, the effect of wheel access, and the line type × wheel access interaction are all 1 and 6. With this type of analysis, one trait that showed a statistically significant interaction is hematocrit in a sample of 81 female mice housed with or without wheels for 8 weeks (Swallow et al.,2005), as repeated here in Table 2.

Table 2.

Comparison of all (female) mice for interaction between line type and wheel access to determine whether Selected lines show greater phenotypic plasticity

P
Trait (± wheels)In likelihood (larger is better)Likelihood ratio testaNo. parametersbAIC (smaller is better)cNS vs CdWheel access or revs in final weekdInteractiondMini-muscled
Hematocrit
±     81 0.1756+ 0.0211+ 0.0465 0.0044-
+ only -77.9  10 175.8 40 0.1559+ 0.6449+  0.0044-
-78.1 0.3 174.1 40 0.0472+   0.0039-
Hemoglobin
±     81 0.0520+ 0.0078+ 0.0146 0.0892-
+ only -34.5  10 88.9 40 0.1375+ 0.5450+  0.0154-
-34.8 0.7 87.6 40 0.0424+   0.0145-
Citrate synthase
±     79 0.0877+ <0.0001+ 0.0039 <0.0001+
+ only -132.3  10 284.6 39 0.4062+ 0.0224+  <0.0001+
-135.8 6.9* 289.6 39 0.0546+   <0.0001+
Cytochrome-c oxidase
±     81 0.0573+ 0.0081+ 0.0224 0.0001+
+ only -162.3  10 344.7 40 0.0490+ 0.8552-  0.0004+
-162.4 0.0 342.7 40 0.0316+   0.0003+
Pyruvate dehydrogenase
±     81 0.0296+ 0.0060+ 0.0704 0.1714+
+ only -51.7  10 123.5 40 0.2616+ 0.1166+  0.6657+
-53.3 3.2 124.6 40 0.0173+   0.3476+
P
Trait (± wheels)In likelihood (larger is better)Likelihood ratio testaNo. parametersbAIC (smaller is better)cNS vs CdWheel access or revs in final weekdInteractiondMini-muscled
Hematocrit
±     81 0.1756+ 0.0211+ 0.0465 0.0044-
+ only -77.9  10 175.8 40 0.1559+ 0.6449+  0.0044-
-78.1 0.3 174.1 40 0.0472+   0.0039-
Hemoglobin
±     81 0.0520+ 0.0078+ 0.0146 0.0892-
+ only -34.5  10 88.9 40 0.1375+ 0.5450+  0.0154-
-34.8 0.7 87.6 40 0.0424+   0.0145-
Citrate synthase
±     79 0.0877+ <0.0001+ 0.0039 <0.0001+
+ only -132.3  10 284.6 39 0.4062+ 0.0224+  <0.0001+
-135.8 6.9* 289.6 39 0.0546+   <0.0001+
Cytochrome-c oxidase
±     81 0.0573+ 0.0081+ 0.0224 0.0001+
+ only -162.3  10 344.7 40 0.0490+ 0.8552-  0.0004+
-162.4 0.0 342.7 40 0.0316+   0.0003+
Pyruvate dehydrogenase
±     81 0.0296+ 0.0060+ 0.0704 0.1714+
+ only -51.7  10 123.5 40 0.2616+ 0.1166+  0.6657+
-53.3 3.2 124.6 40 0.0173+   0.3476+

Values represent nested ANCOVAs comparing all (female) mice (S vsC lines housed with vs without wheel access for 8 weeks). Selected(S) and Control (C) lines for mice housed only with wheel access are then compared to determine whether S lines show greater phenotypic plasticity; see text for further explanation.

AIC, Akaike Information Criterion.

*

P<0.05. For cytrate synthase, this statistical difference indicates that the higher values exhibited by mice from the S lines can be explained as being a simple linear function of their greater amount of wheel running during the final week of wheel access. In contrast, for the other four traits, mice from S lines seem to exhibit greater phenotypic plasticity(training effects).

a

Twice the difference in likelihood is distributed as a χ2with 1 d.f., i.e. 3.841 for P=0.05. Values larger than this indicate that the model including amount of running during final week as a covariate(full model) fits the data significantly better than a model that does not include this covariate (reduced model).

b

Including intercept, fixed effects, random effects (even if estimated as zero), residual variance, and all covariates. Other covariates (not shown) in all analyses included log body mass (except for analysis of log body mass itself), age, time of day and (z-transformed time of day)2.

c

AIC calculated using all parameters in model, as indicated in previous column.

d

Degrees of freedom for S vs C, for wheel access vssedentary, and for the interaction of those two factors were 1 and 6 in all models. For effect of mini-muscle phenotype, d.f. were approximately 1 and 26-28, depending on sample size. In the analyses of wheel-access mice only,d.f. for amount of running during final week were approximately 1 and 26,depending on sample size; see Swallow et al.(2005) and text for further details. Signs indicate direction of effect: + indicates S>C, wheel access>sedentary, or positive effect of body mass.

P values for S vs C, Wheel access vs sedentary,Interaction of those two effects, and Mini-muscle are for two-tailed tests;values <0.05 effects are in bold type. P values for amount of running in final week are also for two-tailed tests.

Second, we examine the mean values for the four subgroups. In the case of hematocrit in females, adjusted means (SAS Proc Mixed) were 48.51, 48.76,48.66, and 50.90 for Control Sedentary, Control Active, Selected Sedentary and Selected Active, respectively [see table 3(Swallow et al., 2005)]. Thus, the line type effect is greater when mice are housed with wheel access. Indeed, separate ANCOVAs reveal no significant effect of line type(P=0.8502) for sedentary mice but a significant effect(P=0.0472) for the active group. Third, within the active group(Fig. 7), we can ask whether the data are better fit by a model that does or does not include the amount of wheel running as an additional covariate. For hematocrit in females, Table 2 shows that the ln likelihood of the nested ANCOVA model without wheel running (–75.7) is larger (less negative, in this case) than for the model with wheel running(–83.7). As the latter model contains one additional parameter(estimating the effect of wheel running), twice the difference in ln likelihoods (16.0, in this case) can be compared with a χ2distribution with one degree of freedom, for which the critical value for P=0.05 is 3.841. Therefore, the model with wheel running as an additional covariate yields a significantly worse fit to the data, and we conclude that the difference in hematocrit between S and C mice when housed with wheel access is not best explained as a simple function of the greater running by S mice. Instead, the greater training effect experienced by S mice seems to indicate that they have greater plasticity for this trait when given wheel access. Results are similar for blood hemoglobin content.

Fig. 7.

Hematocrit of female mice from a study reported elsewhere(Swallow et al., 2005). Mice were given wheel access for 8 weeks beginning at weaning. Additional mice (not shown) were housed without access to running wheels, and they showed no Selected (S) versus Control (C) difference in hematocrit. Thus, S mice are more responsive (more plastic) when granted wheel access, but not by virtue of a simple linear relation with amount of running that applies to all animals. Compare with Fig. 6B.

Fig. 7.

Hematocrit of female mice from a study reported elsewhere(Swallow et al., 2005). Mice were given wheel access for 8 weeks beginning at weaning. Additional mice (not shown) were housed without access to running wheels, and they showed no Selected (S) versus Control (C) difference in hematocrit. Thus, S mice are more responsive (more plastic) when granted wheel access, but not by virtue of a simple linear relation with amount of running that applies to all animals. Compare with Fig. 6B.

We are currently in the process of examining or reexamining a variety of traits, from several different published and unpublished studies, by the approach outlined in the previous paragraph. Table 2 shows some additional examples for the female mice studied elsewhere(Houle-Leroy et al., 2000; Swallow et al., 2005). We are finding a number of traits that fit the pattern described for hematocrit in females, including hemoglobin content of the blood as well as cytochrome c oxidase and pyruvate dehydrogenase activity in mixed hindlimb muscle in females. Thus, several traits seem to show greater plasticity in response to wheel access in S lines as compared with C lines. In addition, the studies by Houle-Leroy et al. and Swallow et al.(Houle-Leroy et al., 2000; Swallow et al., 2005) housed the sedentary' mice with access to locked wheels (unable to rotate), and a subsequent study revealed that mice from S lines climb more than those from C lines in locked wheels. Therefore, the relative magnitude of training effects in S and C lines might differ if sedentary' mice were housed in ordinary cages with no access to even locked wheels.

The results shown in Table 2are for (female) mice housed with or without access to a functional wheel for 8 weeks. However, the selection regime involves only 6 days of wheel access,so it will be crucial in future studies to see if a similar pattern emerges for shorter periods of wheel access. In fact, we have already found that, for some traits, the effects of wheel access can be dramatic even over a matter of days. For example, the amount of GLUT-4 glucose transporter in gastrocnemius muscle did not differ between S and C females when they were housed without wheels (Gomes et al., 2004). After 5 days, both groups exhibited an increase in GLUT-4, but the increase was much greater in S mice, such that S and C showed no overlap in values. When the amount of GLUT-4 was plotted against the amount of wheel running on day 5, the relation was similar to that shown in Fig. 6B. Thus, the difference between S and C lines in gastrocnemius GLUT-4 expression is not a simple linear function of the amount of wheel running; rather, mice from the S lines seem to have greater plasticity for this trait, and this greater plasticity can have large effects even in as few as 5 days.

Some traits show altered plasticity in the S lines, but in a complicated way. For example, the amount of neurogenesis in the hippocampus [see fig. 2D(Rhodes et al., 2003)] shows a relation similar to that depicted in Fig. 8. Mice from C lines (gray squares) exhibit a positive and quantitative relation with the amount of wheel running exhibited over several weeks, but this relation is lost in the S lines (black circles). Finally, we have also observed some traits that show an actual reversal of the direction of plasticity in S lines. For instance, relative ovary mass was found to be larger in S mice than in C when both were housed without wheels, but the opposite was true for mice housed with wheels for 8 weeks(Swallow et al., 2005).

Fig. 8.

Complicated hypothetical relations between a phenotypic trait and the amount of running exhibited during the final week of a multi-week exposure to wheels. Mice from Control lines (gray squares) exhibit a positive and linear quantitative relation with the amount of wheel running, but this relation is lost in the Selected lines (black circles). Some traits, such as the amount of neurogenesis in the hippocampus [see fig. 2D(Rhodes et al., 2003)],actually show this sort of complicated pattern.

Fig. 8.

Complicated hypothetical relations between a phenotypic trait and the amount of running exhibited during the final week of a multi-week exposure to wheels. Mice from Control lines (gray squares) exhibit a positive and linear quantitative relation with the amount of wheel running, but this relation is lost in the Selected lines (black circles). Some traits, such as the amount of neurogenesis in the hippocampus [see fig. 2D(Rhodes et al., 2003)],actually show this sort of complicated pattern.

### Evolutionary change versus phenotypic plasticity

As noted above, some traits do not show a significant interaction between line type and wheel access. For these traits, the magnitude of any S versus C difference is relatively constant, regardless of housing conditions, and the magnitude of any training effect is similar in both S and C. Therefore, we can compare these two effects in a straightforward way plotting one versus the other.

Kelly et al. report hindlimb bone properties for male mice from generation 21 that were given wheel access for 8 weeks and compared with counterparts housed in ordinary cages with no wheels(Kelly et al., 2006). As shown in Fig. 9, bone lengths were not affected by either selective breeding or chronic wheel access. Diameters,in contrast, tended to be increased by both factors, with the magnitude of the evolutionary effect being somewhat greater than the training effect. In spite of the fact that mice from S lines ran considerably more than C, the magnitude of training effects was similar in S and C lines, indicating no genotype-by-environment interactions (see Kelly et al., 2006). Similarly, plastic and evolved (5–6 generations) responses of larval tracheae were in the same direction in Drosophila melanogasterexposed to different atmospheric oxygen content(Henry and Harrison, 2004)[for other examples in Drosophila, see(Price et al., 2003) and references therein]. An important area for future research will be to determine how often the underlying mechanisms of change are similar versus different for plastic and evolved responses [e.g.(Hodin and Riddiford, 2000)and references therein].

Fig. 9.

Comparison of the effects of 21 generations of selective breeding for high voluntary wheel running (Evolution') with the effects of 8 weeks of wheel access (Phenotypic plasticity'), based on least-squares (adjusted) means for hindlimb bones [as presented in table 3(Kelly et al., 2006)]. Bone lengths were unaffected by either selective breeding or chronic wheel access,but both factors increased bone diameters (one exception), with the former generally having somewhat larger effects. The magnitude of training effects is similar to what has been reported previously for mammalian bone in response to either forced or voluntary exercise (see Kelly et al., 2006).

Fig. 9.

Comparison of the effects of 21 generations of selective breeding for high voluntary wheel running (Evolution') with the effects of 8 weeks of wheel access (Phenotypic plasticity'), based on least-squares (adjusted) means for hindlimb bones [as presented in table 3(Kelly et al., 2006)]. Bone lengths were unaffected by either selective breeding or chronic wheel access,but both factors increased bone diameters (one exception), with the former generally having somewhat larger effects. The magnitude of training effects is similar to what has been reported previously for mammalian bone in response to either forced or voluntary exercise (see Kelly et al., 2006).

Whenever a selective event is more than instantaneous, individuals exposed to it have the opportunity to exhibit phenotypic change (plasticity) even as that event is in progress. If the direction of plastic change is such that it increases the probability of survival or other components of Darwinian fitness, then those individuals will be at a selective advantage. Hence, we may expect that directional selection will often lead to the evolution of increased plasticity in the appropriate direction(Fig. 2). This prediction would apply to organisms in nature or in laboratory or other experimental systems. Hence, the tools of selective breeding (sensu lato) can be used to study the evolution of phenotypic plasticity, just as they can be used to study any other complex phenotype(Swallow and Garland, 2005). Even if plasticity per se is not the intentional target of a selection experiment, we suggest that any experiment in which the selective event is more than instantaneous should explore whether plasticity has increased as a component of the response to selection. Our own experiment that involves selective breeding for high voluntary activity levels in mice has uncovered some examples of increased and apparently adaptive plasticity in the selected lines, but not for all traits (e.g. not for hindlimb skeletal dimensions). We believe that experimental evolution approaches have much to offer with respect to clarifying the evolution of plasticity, including self-induced adaptive plasticity(Houle-Leroy et al., 2000; Swallow et al., 2005) and the related concept of genetic assimilation(Price et al., 2003; Pigliucci et al., 2006).

Although we have argued here that plasticity might be predicted to increase in the appropriate direction as a correlated response to directional selection(e.g. see Fig. 2), the opposite prediction has also been offered. For example, in our original paper on enzyme activities, we noted that genetic selection for increased voluntary wheel running did not reduce the capability of muscle aerobic capacity to respond to training' [see pp. 1608 and 1613(Houle-Leroy et al., 2000)]because, a priori, we had thought that organisms innately higher for a physiological function might exhibit relativity less ability to increase further in that direction via physical conditioning, perhaps because of some sort of ceiling effect'. Similarly, Koch et al. in a training study of rat exercise capacity (Koch et al.,2005) noted that According to the principle of initial values, if capacity is low then the percentage gain in capacity in response to training will be high, and vice versa.' However, their results for two inbred strains that differed in intrinsic aerobic running capacity were not in accord with the principle of initial value.

We thank A. F. Bennett, T. D. Brutsaert, A. G. Gibbs, R. B. Huey and E. R. Weibel for helpful discussions. We also thank L. G. Harshman, M. Pigliucci, T. Price, M. Van Sant and an anonymous reviewer for comments on the manuscript. Supported by US National Science Foundation grant IBN-0212567 to T.G.

Andersen, J. B., Rourke, B. C., Caiozzo, V. J., Bennett, A. F. and Hicks, J. W. (
2005
). Postprandial cardiac hypertrophy in pythons.
Nature
434
,
37
-38.
Beall, C. M. (
2001
). Adaptations to altitude: a current assessment.
Annu. Rev. Anthropol.
30
,
423
-456.
Belter, J. G., Carey, H. V. and Garland, T., Jr(
2004
). Effects of voluntary exercise and genetic selection for high activity levels on HSP72 expression in house mice.
J. Appl. Physiol.
96
,
1270
-1276.
Bennett, A. F. (
1997
). Adaptation and the evolution of physiological characters. In
Handbook of Physiology:Section 13, Comparative Physiology.
Vol.
I
(ed. W. H. Dantzler), pp.
3
-16. New York: Oxford University Press.
Bennett, A. F. (
2002
). Experimental evolution. In
Encyclopedia of Evolution
(ed. M. Pagel), pp.
339
-342. New York: Oxford University Press.
Bennett, A. F. (
2003
). Experimental evolution and the Krogh Principle: generating biological novelty for functional and genetic analyses.
Physiol. Biochem. Zool.
76
,
1
-11.
Berrigan, D. and Scheiner, S. M. (
2004
). Modeling the evolution of phenotypic plasticity. In
Phenotypic Plasticity: Functional and Conceptual Approaches
(ed. T. J. Dewitt and S. M. Scheiner), pp.
82
-97. New York:Oxford University Press.
Boorse, G. C. and Denver, R. J. (
2004
). Endocrine mechanisms underlying plasticity in metamorphic timing in spadefoot toads.
Integr. Comp. Biol.
43
,
646
-657.
Booth, F. W., Chakravathy, M. U., Gordon, S. E. and Spangenburg,E. E. (
2002
). Waging war on physical inactivity: using modern molecular ammunition against an ancient enemy.
J. Appl. Physiol.
93
,
3
-30.
Broadhurst, P. L. and Jinks, J. L. (
1974
). What genetical architecture can tell us about the natural selection of behavioral traits. In
The Genetics of Behaviour
(ed. J. H. F. van Abeelen), pp.
43
-63. Amsterdam:North-Holland.
Brutsaert, T. D., Parra, E. J., Shriver, M. D., Gamboa, A.,Rivera-Ch, M. and Leon-Velarde, F. (
2005
). Ancestry explains the blunted ventilatory response to sustained hypoxia and lower exercise ventilation of Quechua altitude natives.
Am. J. Physiol.
289
,
R225
-R234.
Bubliy, O. A. and Loeschcke, V. (
2005
). Correlated responses to selection for stress resistance and longevity in a laboratory population of Drosophila melanogaster.
J. Evol. Biol.
18
,
789
-803.
Bumpus, H. C. (
1899
). The elimination of the unfit as illustrated by the introduced sparrow, Passer domesticus. In
Biological Lectures Delivered at the Marine Biological Laboratory of Woods Hole, 1896-7
. pp.
209
-226. Boston: Ginn & Co.
Castaneda, T. R., Jurgens, H., Wiedmer, P., Pfluger, P., Diano,S., Horvath, T. L., Tang-Christensen, M. and Tschop, M. H.(
2005
). Obesity and the neuroendocrine control of energy homeostasis: the role of spontaneous locomotor activity.
J. Nutr.
135
,
1314
-1319.
Chippindale, A. K., Ngo, A. L. and Rose, M. R.(
2003
). The devil in the details of life-history evolution:instability and reversal of genetic correlations during selection on Drosophila development.
J. Genet.
82
,
133
-145.
Collier, M. H. and Rogstad, S. H. (
2004
). Clonal variation in floral stage timing in the common dandelion Taraxacum officinale (Asteraceae).
Am. J. Bot.
91
,
1828
-1833.
Conley, K. E., Christian, K. A., Hoppeler, H. and Weibel, E. R. (
1995
). Heart mitochondrial properties and aerobic capacity are similarly related in a mammal and a reptile.
J. Exp. Biol.
198
,
739
-746.
Cook, S. A. and Johnson, M. P. (
1968
). Adaptation to heterogeneous environments. I. Variation in heterophylly in Ranunculus flammula L.
Evolution
22
,
496
-516.
Cossins, A., Fraser, J. and Gracey, A. (
2006
). Post-genomic approaches to understanding the mechanisms of environmentally induced phenotypic plasticity.
J. Exp. Biol.
209
,
2328
-2336.
Costa, D. P. and Sinervo, B. (
2004
). Field physiology: physiological insights from animals in nature.
Annu. Rev. Physiol.
66
,
209
-238.
Diamond, J. M. (
1986
). Why do disused proteins become genetically lost or repressed?
Nature
321
,
565
-566.
Dohm, M. R., Richardson, C. S. and Garland, T., Jr(
1994
). Exercise physiology of wild and random-bred laboratory house mice and their reciprocal hybrids.
Am. J. Physiol.
267
,
R1098
-R1108.
Donohue, K., Messiqua, D., Pyle, E. H., Heschel, M. S. and Schmitt, J. (
2000
). Evidence of adaptive divergence in plasticity: density-and-site-dependent selection on shade-avoidance responses in Impatiens capensis.
Evolution
54
,
1956
-1968.
Dumke, C. L., Rhodes, J. S., Garland, T., Jr, Maslowski, E.,Swallow, J. G., Wetter, A. C. and Cartee, G. D.(
2001
).Genetic selection of mice for high voluntary wheel running: effect on skeletal muscle glucose uptake.
J. Appl. Physiol.
91
,
1289
-1297.
Ebert, D. (
1998
). Experimental evolution of parasites.
Science
282
,
1432
-1435.
Endler, J. A. (
1986
).
Natural Selection in the Wild
. Princeton, NJ: Princeton University Press.
Falconer, D. S. (
1989
).
Introduction to Quantitative Genetics
(3rd edn). London: Longman.
Falconer, D. S. (
1990
). Selection in different environments: effects on environmental sensitivity (reaction norm) and on mean performance.
Genet. Res., Camb.
56
,
57
-70.
Feder, M. E., Bennett, A. F. and Huey, R. B.(
2000
). Evolutionary physiology.
Annu. Rev. Ecol. Syst.
31
,
315
-341.
Flück, M. (
2006
). Functional, structural and molecular plasticity of mammalian skeletal muscle in response to exercise stimuli.
J. Exp. Biol.
209
,
2239
-2248.
Fordyce, J. A. (
2006
). The evolutionary consequences of ecological interactions mediated through phenotypic plasticity.
J. Exp. Biol.
209
,
2377
-2383.
Freeman, S. and Herron, J. C. (
2004
).
Evolutionary Analysis
(3rd edn). Upper Saddle River,NJ: Prentice Hall.
Futuyma, D. J. (
1998
).
Evolutionary Biology
(3rd edn). Sunderland, MA: Sinauer Associates.
Gabriel, W. (
1999
). Evolution of reversible plastic responses: inducible defenses and environmental tolerance. In
The Ecology and Evolution of Inducible Defenses
(ed. R. Tollrian and C. D. Harvell), pp.
286
-305. Princeton, NJ: Princeton University Press.
Gabriel, W. (
2005
). How stress selects for reversible phenotypic plasticity.
J. Evol. Biol.
18
,
873
-883.
Garland, T., Jr (
1994
). Quantitative genetics of locomotor behavior and physiology in a garter snake. In
Quantitative Genetic Studies of Behavioral Evolution
(ed. C. R. B. Boake), pp.
251
-277. Chicago:University of Chicago Press.
Garland, T., Jr (
1998
). Testing the predictions of symmorphosis: conceptual and methodological issues. In
Principles of Animal Design: the Optimization and Symmorphosis Debate
(ed. E. R. Weibel, L. Bolis and C. R. Taylor), pp.
40
-47. New York: Cambridge University Press.
Garland, T., Jr (
1999
). Artificial selection for increased wheel-running activity in house mice results in decreased body mass at maturity.
J. Exp. Biol.
202
,
2513
-2520.
Garland, T., Jr (
2003
). Selection experiments:an under-utilized tool in biomechanics and organismal biology. In
Vertebrate Biomechanics and Evolution
(ed. V. L. Bels,J. P. Gasc and A. Casinos), pp.
23
-56. Oxford: BIOS Scientific Publishers.
Garland, T., Jr and Adolph, S. C. (
1994
). Why not to do two-species comparative studies: limitations on inferring adaptation.
Physiol. Zool.
67
,
797
-828.
Garland, T., Jr and Carter, P. A. (
1994
). Evolutionary physiology.
Annu. Rev. Physiol.
56
,
579
-621.
Garland, T., Jr and Else, P. L. (
1987
). Seasonal, sexual, and individual variation in endurance and activity metabolism in lizards.
Am. J. Physiol.
252
,
R439
-R449.
Garland, T., Jr and Freeman, P. A. (
2005
). Selective breeding for high endurance running increases hindlimb symmetry.
Evolution
59
,
1851
-1854.
Garland, T., Jr, Bennett, A. F. and Daniels, C. B.(
1990
). Heritability of locomotor performance and its correlates in a natural population of vertebrates.
Experientia
46
,
530
-533.
Garland, T., Jr, Morgan, M. T., Swallow, J. G., Rhodes, J. S.,Girard, I., Belter, J. G. and Carter, P. A. (
2002
). Evolution of a small-muscle polymorphism in lines of house mice selected for high activity levels.
Evolution
56
,
1267
-1275.
Garland, T., Jr, Bennett, A. F. and Rezende, E. L.(
2005
). Phylogenetic approaches in comparative physiology.
J. Exp. Biol.
218
,
3015
-3035.
Geiser, F., Kenagy, G. J. and Wingfield, J. C.(
1997
). Dietary cholesterol enhances torpor in a rodent hibernator.
J. Comp. Physiol. B
167
,
416
-422.
Ghalambor, C. K., Walker, J. A. and Reznick, D. N.(
2003
). Multi-trait selection, adaptation, and constraints on the evolution of burst swimming performance.
Integr. Comp. Biol.
43
,
431
-438.
Girard, I. and Garland, T., Jr (
2002
). Plasma corticosterone response to acute and chronic voluntary exercise in female house mice.
J. Appl. Physiol.
92
,
1553
-1561.
Girard, I., McAleer, M. C., Rhodes, J. S. and Garland, T.,Jr (
2001
). Selection for high voluntary wheel-running increases speed and intermittency in house muice (Mus domesticus).
J. Exp. Biol.
204
,
4311
-4320.
Gomes, F. R., Rezende, E. L., Bunkers, J. L., Rivas, D. A.,Yaspelkis, B. B., III and Garland, T., Jr (
2004
). Muscle glucose transporters (GLUT-4) and glycogen storage of mice selectively bred for high activity levels.
Int. Comp. Biol.
44
,
984
.
Gordon, D. M. (
1992
). Phenotypic plasticity. In
Keywords in Evolutionary Biology
(ed. E. F. Keller and E. A. Loyd), pp.
255
-262. Cambridge, MA: Harvard University Press.
Harshman, L. G., Ottea, J. A. and Hammock, B. D.(
1991
). Evolved environment-dependent expression of detoxication enzyme activity in Drosophila melanogaster.
Evolution
45
,
791
-795.
Harshman, L. G., Hoffmann, A. A. and Clark, A. G.(
1999
). Selection for starvation resistance in Drosophila melanogaster: physiological correlates, enzyme activities and multiple stress responses.
J. Evol. Biol.
12
,
370
-379.
Hatle, J. D. (
2004
). Physiology underlying phenotypic plasticity and polyphenisms: introduction to the symposium.
Integr. Comp. Biol.
43
,
605
-606.
Henderson, K. K., Wagner, H., Favret, F., Britton, S. L., Koch,L. G., Wagner, P. D. and Gonzalez, N. C. (
2002
). Determinants of maximal O2 uptake in rats selectively bred for endurance running capacity.
J. Appl. Physiol.
93
,
1265
-1274.
Henderson, N. D. (
1981
). Genetic influences on locomotor activity in 11-day-old housemice.
Behav. Genet.
11
,
209
-225.
Henderson, N. D. (
1989
). Interpreting studies that compare high- and low-selected lines on new characters.
Behav. Genet.
19
,
473
-502.
Henry, J. R. and Harrison, J. F. (
2004
). Plastic and evolved responses of larval tracheae and mass to varying atmospheric oxygen content in Drosophila melanogaster.
J. Exp. Biol.
207
,
3559
-3567.
Hewitt, J. K., Fulker, D. W. and Broadhurst, P. L.(
1981
). Genetics of escape-avoidance conditioning in laboratory and wild populations of rats: a biometrical approach.
Behav. Genet.
11
,
533
-544.
Hodin, J. and Riddiford, L. M. (
2000
). Different mechanisms underlie phenotypic plasticity and interspecific variaiton for a reproductive character in drosophilids (Insecta: Diptera).
Evolution
54
,
1638
-1653.
Hood, D. A., Irrcher, I., Ljubicic, V. and Joseph, A.-M.(
2006
). Coordination of metabolic plasticity in skeletal muscle.
J. Exp. Biol.
209
,
2265
-2275.
Hoppeler, H., Altpeter, E., Wagner, M., Turner, D. L., Hokanson,J., König, M., Stalder-Navarro, V. P. and Weibel, E. R.(
1995
). Cold acclimation and endurance training in guinea pigs:changes in lung, muscle and brown fat tissue.
Respir. Physiol.
101
,
189
-198.
Houle-Leroy, P., Garland, T., Jr, Swallow, J. G. and Guderley,H. (
2000
). Effects of voluntary activity and genetic selection on muscle metabolic capacitites in house mice Mus domesticus.
J. Appl. Physiol.
89
,
1608
-1616.
Houle-Leroy, P., Guderley, H., Swallow, J. G. and Garland, T.,Jr (
2003
). Artificial selection for high activity favors mighty mini-muscles in house mice.
Am. J. Physiol.
284
,
R433
-R443.
Hsia, C. C. (
2001
). Coordinated adaptation of oxygen transport in cardiopulmonary disease.
Circulation
104
,
963
-969.
Huey, R. B. and Berrigan, D. (
1996
). Testing evolutionary hypotheses of acclimation. In
Phenotypic and Evolutionary Adaptation to Temperature
(ed. I. A. Johnston and A. F. Bennett), pp.
205
-237. Cambridge: Cambridge University Press.
Huey, R. B. and Kingsolver, J. G. (
1993
). Evolution of resistance to high temperature in ectotherms.
Am. Nat.
142
,
S21
-S46.
Huey, R. B., Berrigan, D., Gilchrist, G. W. and Herron, J. C. (
1999
). Testing the adaptive significance of acclimation:a strong inference approach.
Am. Zool.
39
,
323
-336.
Huey, R. B., Hertz, P. E. and Sinervo, B.(
2003
). Behavioral drive versus behavioral inertia in evolution:a null model approach.
Am. Nat.
161
,
357
-366.
Irschick, D. J. and Garland, T., Jr (
2001
). Integrating function and ecology in studies of adaptation: investigations of locomotor capacity as a model system.
Annu. Rev. Ecol. Syst.
32
,
367
-396.
Johnson, R. A., Rhodes, J. S., Jeffrey, S. L., Garland, T., Jr and Mitchell, G. S. (
2003
). Hippocampal brain-derived neurotrophic factor but not neurotrophin-3 increases more in mice selected for increased voluntary wheel running.
Neuroscience
121
,
1
-7.
Johnston, I. (
2006
). Environment and plasticity of myogenesis in teleost fish.
J. Exp. Biol.
209
,
2249
-2264.
Kelly, S. A., Czech, P. P., Wight, J. T., Blank, K. M. and Garland, T., Jr (
2006
). Experimental evolution and phenotypic plasticity of hindlimb bones in high-activity house mice.
J. Morphol.
267
,
360
-374.
Ketterson, E. D. and Nolan, V., Jr (
1999
). Adaptation, exaptation, and constraint: a hormonal perspective.
Am. Nat.
154
,
S4
-S25.
Kingsolver, J. G. and Huey, R. B. (
2003
). Introduction: the evolution of morphology, performance, and fitness.
Integr. Comp. Biol.
43
,
361
-366.
Knapp, R. (
2004
). Endocrine mediation of vertebrate male alternative reproductive tactics: the next generation of studies.
Integr. Comp. Biol.
43
,
658
-668.
Koch, L. G., Green, C. L., Lee, A. D., Hornyak, J. E., Cicila,G. T. and Britton, S. L. (
2005
). Test of the principle of initial value in rat genetic models of exercise capacity.
Am. J. Physiol.
288
,
R466
-R472.
Koteja, P., Swallow, J. G., Carter, P. A. and Garland, T.,Jr (
1999a
). Energy cost of wheel running in house mice:implications for coadaptation of locomotion energy budgets.
Physiol. Biochem. Zool.
72
,
238
-249.
Koteja, P., Garland, T., Jr, Sax, J. N., Swallow, J. G. and Carter, P. A. (
1999b
). Behaviour of house mice artificially selected for high levels of voluntary wheel running.
Anim. Behav.
58
,
1307
-1318.
Kupferberg. S. J. (
1997
). The role of larval diet in anuran metamorphosis.
Am. Zool.
37
,
146
-159.
Lande, R. and Arnold, S. J. (
1983
). The measurement of selection on correlated characters.
Evolution
37
,
1210
-1226.
Leroi, A., Bennett, A. F. and Lenski, R. E.(
1994
). Temperature acclimation and competitive fitness: an experimental test of the beneficial acclimation assumption.
Proc. Natl. Acad. Sci. USA
91
,
1917
-1921.
Lillywhite, H. B. (
2004
). Plasticity of the water barrier in vertebrate integument. 3rd International Conference of Comparative Physiology and Biochemistry, August 7-13, 2004, KwaZulu Natal,South Africa.
Int. Congr. Ser.
1275
,
283
-290.
Lynch, C. B. (
1994
). Evolutionary inferences from genetic analyses of cold adaptation in laboratory and wild populations of the house mouse. In
Quantitative Genetic Studies of Behavioral Evolution
(ed. C. R. B. Boake), pp.
278
-301. Chicago: University of Chicago Press.
Lynch, M. and Gabriel, W. (
1987
). Environmental tolerance.
Am. Nat.
129
,
283
-303.
Lynch, M. and Walsh, J. B. (
1998
).
Genetics and Analysis of Quantitative Traits
. Sunderland, MA: Sinauer Associates.
Magistretti, P. J. (
2006
). Neuron-glia metabolic coupling and plasticity.
J. Exp. Biol.
209
,
2304
-2311.
Mather, K. and Jinks, J. L. (
1982
).
Biometrical Genetics: The Study of Continuous Variation
(3rd edn). London: Chapman & Hall.
Nelson, P., Kiriakidou, M., Sharma, A., Maniataki, E. and Mourelatos, Z. (
2003
). The microRNA world: small is mighty.
Trends Biochem. Sci.
28
,
534
-540.
Nesse, R. M. (
2005
). Maladaptation and natural selection.
Q. Rev. Biol.
80
,
62
-70.
Norcliffe, L. J., Rivera-Ch, M., Claydon, V. E., Moore, J. P.,Leon-Velarde, F., Appenzeller, O. and Hainsworth, R. (
2005
). Cerebrovascular responses to hypoxia and hypocapnia in high-altitude dwellers.
J. Physiol.
566
,
287
-294.
Nunney, L. and Cheung, W. (
1997
). The effect of temperature on body size and fecundity in female Drosophila melanogaster: evidence for adaptive plasticity.
Evolution
51
,
1529
-1535.
Nussey, D. H., Postma, E., Gienapp, P. and Visser, M. E.(
2005
). Selection on heritable phenotypic plasticity in a wild bird population.
Science
310
,
304
-306.
Orzack, S. and Sober, E. (ed.) (
2001
).
. New York: Cambridge University Press.
Piersma, T. and Lindstrom, A. (
1997
). Rapid reversible changes in organ size as a component of adaptive behavior.
Trends Ecol. Evol.
12
,
134
-138.
Pigliucci, M. (
2001
).
Phenotypic Plasticity: Beyond Nature and Nurture
. Baltimore, MD: The Johns Hopkins University Press.
Pigliucci, M. (
2005
). Evolution of phenotypic plasticity: where are we going now?
Trends Ecol. Evol.
20
,
481
-486.
Pigliucci, M., Murran, C. J. and Schlichting, C. D.(
2006
). Phenotypic plasticity and evolution by genetic assimilation.
J. Exp. Biol.
209
,
2362
-2367.
Price, T. D. (
2006
). Phenotypic plasticity,sexual selection, and the evolution of colour patterns.
J. Exp. Biol.
209
,
2368
-2376.
Price, T. D., Qvarnstrom, A. and Irwin, D. E.(
2003
). The role of phenotypic plasticity in driving genetic evolution.
Proc. R. Soc. Lond. B Biol. Sci.
270
,
1433
-1440.
Rezende, E. L., Gomes, F. R., Ghalambor, C. K., Russell, G. A. and Chappell, M. A. (
2005
). An evolutionary frame of work to study physiological adaptation to high altitudes.
Rev. Chil. Hist. Nat.
78
,
323
-336.
Rezende, E. L., Garland, T., Jr, Chappell, M. A., Malisch, J. L. and Gomes, F. R. (
2006a
). Maximum aerobic performance in lines of Mus selected for high wheel-running activity: effects of selection, oxygen availability, and the mini-muscle phenotype.
J. Exp. Biol.
209
,
115
-127.
Rezende, E. L., Kelly, S. A., Gomes, F. R., Chappell, M. A. and Garland, T., Jr (
2006b
). Effects of size, sex, and voluntary running speeds on costs of locomotion in lines of laboratory mice selectively bred for high wheel-running activity.
Physiol. Biochem. Zool.
79
,
83
-99.
Rhodes, J. S., Koteja, P., Swallow, J. G., Carter, P. A. and Garland, T., Jr (
2000
). Body temperatures of house mice artificially selected for high voluntary wheel-running behavior:repeatability and effect of genetic selection.
J. Therm. Biol.
25
,
391
-400.
Rhodes, J. S., van Praag, H., Jeffrey, S., Girard, I., Mitchell,G. S., Garland, T., Jr and Gage, F. H. (
2003
). Exercise increases hippocampal neurogenesis to high levels but does not improve spatial learning in mice bred for increased voluntary wheel running.
Behav. Neurosci.
117
,
1006
-1016.
Rhodes, J. S., Gammie, S. C. and Garland, T., Jr(
2005
). Neurobiology of mice selected for high voluntary wheel-running activity.
Integr. Comp. Biol.
45
,
438
-455.
Rose, M. R., Chippindale, A. K. and Nusbaum, T. J.(
1996
). Laboratory evolution: the experimental wonderland and the Cheshire Cat Syndrome. In
(ed. M. R. Rose and G. V. Lauder), pp.
221
-241. New York: Academic Press.
Rose, M. R., Passananti, H. B. and Matos, M. (ed.)(
2004
).
Methuselah Flies: A Case Study in the Evolution of Aging
. Singapore: World Scientific Publishing.
Rose, M. R., Passananti, H. B., Chippindale, A. K., Phelan J. P., Matos, M., Teotonio, H. and Mueller, L. D. (
2005
). The effects of evolution are local: Evidence from experimental evolution in Drosophila.
Int. Comp. Biol.
45
,
486
-491.
Sapolsky, R. M. (
2005
). The influence of social hierarchy on primate health.
Science
308
,
648
-652.
Scheiner, S. M. (
1993
). Genetics and evolution of phenotypic plasticity.
Annu. Rev. Ecol. Syst.
24
,
35
-68.
Scheiner, S. M. (
2002
). Selection experiments and the study of phenotypic plasticity.
J. Evol. Biol.
15
,
889
-898.
Scheiner, S. M. and Lyman, R. F. (
1991
). The genetics of phenotypic plasticity. II. Responses to selection.
J. Evol. Biol.
4
,
23
-50.
Schlichting, C. D. and Pigliucci, M. (
1998
).
Phenotypic Evolution: A Reaction Norm Perspective
. Sunderland, MA: Sinauer Associates.
Schluter, D. (
1996
). Adaptive radiation along genetic lines of least resistance.
Evolution
50
,
1766
-1774.
Schratt, G. M., Tuebing, F., Nigh, E. A., Kane, C. G., Sabatini,M. E., Kiebler, M. and Greenberg, M. E. (
2006
). A brain-specific microRNA regulates dendritic spine development.
Nature
439
,
283
-289.
Secor, S. M. (
2005
). Evolutionary and cellular mechanisms regulating intestinal performance of amphibians and reptiles.
Integr. Comp. Biol.
45
,
282
-294.
Sinervo, B. and Basolo, A. L. (
1996
). Testing adaptation using phenotypic manipulations. In
(ed. M. R. Rose and G. V. Lauder), pp.
149
-185. San Diego: Academic Press.
Sinervo, B. and Calsbeek, R. (
2003
). Physiological epistasis, ontogenetic conflict and natural selection on physiology and life history.
Integr. Comp. Biol.
43
,
419
-430.
Stearns, S. C. and Hoekstra, R. F. (
2005
).
Evolution: An Introduction
(2nd edn). Oxford: Oxford University Press.
Swallow, J. G. and Garland, T., Jr (
2005
). Selection experiments as a tool in evolutionary and comparative physiology:insights into complex traits – an introduction to the symposium.
Integr. Comp. Biol.
45
,
387
-390.
Swallow, J. G., Carter, P. A. and Garland, T., Jr(
1998a
). Artificial selection for increased wheel-running behavior in house mice.
Behav. Genet.
28
,
227
-237.
Swallow, J. G., Garland, T., Jr, Carter, P. A., Zhan, W.-Z. and Sieck, G. (
1998b
). Effects of voluntary activity and genetic selection on aerobic capacity in house mice.
J. Appl. Physiol.
84
,
69
-76.
Swallow, J. G., Koteja, P., Carter, P. A. and Garland, T.,Jr (
1999
). Artificial selection for increased wheel-running activity in house mice results in decreased body mass at maturity.
J. Exp. Biol.
202
,
2513
-2520.
Swallow, J. G., Garland, T., Jr, Koteja, P. and Carter, P. A. (
2001
). Food consumption and body composition in mice selected for high wheel-running activity.
J. Comp. Physiol. B
171
,
651
-659.
Swallow, J. G., Rhodes, J. S. and Garland, T., Jr(
2005
). Phenotypic and evolutionary plasticity of organ masses in response to voluntary exercise in house mice.
Integr. Comp. Biol.
45
,
426
-437.
Swynghedauw, B. (
2006
). Phenotypic plasticity of adult myocardium. Molecular mechanisms.
J. Exp. Biol.
209
,
2320
-2327.
Syme, D. A., Evashuk, K., Grintuch, B., Rezende, E. L. and Garland, T., Jr (
2005
). Contractile abilities of normal and `mini' triceps surae muscles from mice (Mus domesticus)selectively bred for high voluntary wheel running.
J. Appl. Physiol.
99
,
1308
-1316.
Thomson, S. L., Garland, T., Jr, Swallow, J. G. and Carter, P. A. (
2002
). Response of Sod-2 enzyme activity to selection for high voluntary wheel running.
Heredity
88
,
52
-61.
Trussell, G. C. (
1997
). Phenotypic plasticity in the foot size of an intertidal snail.
Ecology
78
,
1033
-1048.
Valladares, F., Wright, S. J., Lasso, E., Kitajima, K. and Pearcy, R. W. (
2000
). Plastic phenotypic response to light of 16 congeneric shrubs from a Panamanian rainforest.
Ecology
81
,
1925
-1936.
Via, S., Gomulkiewicz, R., de Jong, G., Scheiner, S.,Schlichting, C. D. and van Tienderen, P. H. (
1995
). Adaptive phenotypic plasticity: concensus and controversy.
Trends Ecol. Evol.
10
,
212
-217.
Weibel, E. R. (
2000
).
Symmorphosis:On Form and Function in Shaping Life.
Cambridge, MA: Harvard University Press.
Wilson, R. S. and Franklin, C. E. (
2002
). Testing the beneficial acclimation hypothesis.
Trends Ecol. Evol.
17
,
66
-70.
Winslow, R. M., Chapman, K. W., Gibson, C. C., Samaja, M.,Monge, C. C., Goldwasser, E., Sherpa, M., Blume, F. D. and Santolaya,R. (
1989
). Different hematologic responses to hypoxia in Sherpas and Quechua Indians.
J. Appl. Physiol.
66
,
1561
-1569.
Wu, T., Wang, X., Wei, C., Cheng, H., Wang, X., Li, Y., Ge-Dong Zhao, H., Young, P., Li, G. and Wang, Z. (
2005
). Hemoglobin levels in Qinghai-Tibet: different effects of gender for Tibetans vs. Han.
J. Appl. Physiol.
98
,
598
-604.
Young, K. V., Brodie, E. D., Jr and Brodie, E. D., III(
2004
). How the horned lizard got its horns.
Science
304
,
65
.
Zera, A. J. (
2004
). The endocrine regulation of wing polymorphism in insects: state of the art, recent surprises, and future directions.
Integr. Comp. Biol.
43
,
607
-616.