Skip Navigation

Integrative and Comparative Biology 2003 43(3):367-375; doi:10.1093/icb/43.3.367
© 2003 by The Society for Integrative and Comparative Biology
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (22)
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Arnold, S. J.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?


Performance Surfaces and Adaptive Landscapes1

Stevan J. Arnold2,1
1 Department of Zoology, Oregon State University, Corvallis, Oregon 97331


    SYNOPSIS
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
In an earlier characterization of the relationship between morphology, performance and fitness, I focused only on directional selection (Arnold, 1983Go). The aim of this article is to extend that characterization to include stabilizing and other forms of nonlinear selection. As in the earlier characterization, this more general description of the morphology-performance-fitness relationship splits empirical analysis into two parts: the study of the relationship between morpholgy and performance, and the study of the relationship between performance and fitness. From a conceptual standpoint, my goal is to specify the relationship of performance studies to the adaptive landscape. I begin by reviewing the adaptive landscape concept and its importance in evolutionary biology. A central point emerging from that review is that that key descriptors of the adaptive landscape can be estimated by measuring the impact of selection on the means, variances and covariances of phenotypic traits. Those descriptors can be estimated by making a quadratic (regression) approximation to the selection surface that describes the relationship between the phenotypic traits of individuals and their fitness. Analysis of the effects of morphology on performance follows an analogous procedure: making a quadratic approximation to the individual performance surface and then using that approximation to solve for the descriptors of the performance landscape. I conclude by discussing the evolution of performance and adaptive landscapes. One possibility with biomechanical justification is that the performance landscape evolves along performance lines of least resistance.


    INTRODUCTION
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
An adaptive landscape for phenotypic characters was first proposed by Simpson (1944Go, 1952Go). This landscape characterization of evolution was derived from Wright's concept of an adaptive landscape in which population fitness is a function of gene frequencies (Wright, 1932Go). In Wright's landscape, mean fitness in the population, W, is a function of the frequency, p, of an allele. Simpson (1944)Go used phenotypic adaptive landscapes to successfully argue that the results of paleontology could be reconciled with geneticist's view of population evolution (Wright, 1945Go). Lande (1976Go, 1979Go) showed that properties of an adaptive landscape for phenotypic traits specified the force of selection in equations for change in the means, variances and covariances of the traits. In Lande's landscape, the natural log of population mean fitness, lnW, is a function of the average values of one or more continuously distributed phenotypic traits, z. As in Wright's landscape, the population mean tends to evolve in an uphill direction on this phenotypic landscape (Lande, 1979Go). These theoretical developments merged Simpson's idea of an adaptive landscape with concepts of inheritance that had been developed in the field of quantitative genetics.

The importance of the adaptive landscape as an integrating concept has long been appreciated. This importance is reflected in the long history of the landscape and related concepts in the evolutionary literature (Pearson, 1903Go; Fisher, 1930Go; Wright, 1931Go; Dobzhansky, 1937Go; Simpson, 1944Go; Schmalhausen, 1949Go). The power of the landscape idea is that it ties together many of the fundamental ideas in evolutionary biology (modes of selection, responses to selection, adaptation, speciation, adaptive radiation). Because of their power, landscapes have been sought for various kinds of traits (e.g., proteins, RNA) with varying degrees of success (Govindarajan and Goldstein, 1997Go; Fontana and Schuster, 1998Go; Bornberg-Bauer and Chan, 1999Go). The most common kind of failure in landscape pursuit is inability to find empirically accessible aspects of traits that can in turn be related to a property, such a population mean fitness, that governs evolutionary trajectories. Against the background of these unsuccessful quests, the tangible existence of an adaptive landscape for phenotypic traits shines like a beacon across the field of evolutionary biology. We have succeeded here, so perhaps we can succeed elsewhere.


    THE ADAPTIVE LANDSCAPE FOR PHENOTYPIC TRAITS
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
Success in specifying a useful adaptive landscape for phenotypic traits rests on a foundation of three assumptions (Lande, 1976Go, 1979Go; Lande and Arnold, 1983Go). (a) Phenotypic traits are normally distributed within populations, or are approximately normal after suitable transformation. (b) The function that relates the trait values of individuals to their expected individual fitness can plausibly assumed to be smooth and continuous. In other words, we can specify its first and second derivatives. (c) Fitness is frequency-independent in the sense that the fitness of individuals is affected by their own trait values, but not by the trait values of other individuals in the population. Although these assumptions somewhat restrict our domain of application, they grant enormous conceptual power.

How can we tap into the conceptual power of the adaptive landscape? Many authors have assumed that the adaptive landscape is merely a metaphor (Eldredge and Cracraft, 1980Go; Provine, 1986Go; Dawkins, 1996Go). These authors have ignored the possibility that the adaptive landscape—like mutation, drift, and inheritance—is an object for empirical study. The adaptive landscape for phenotypic traits is accessible to us because of the immediate changes that it induces in the means, variances and covariances of traits within a generation. The features of the landscape that cause changes in means, variances and covariances are its slope and curvature, its first and second derivatives. The first and second derivatives of the landscape also play a fundamental role in predicting how the population mean, z, will change from one generation to the next in respond to selection, and how the inheritance matrix, G, will change within a generation in response to selection;


{i1540-7063-043-03-0367-e1}

(Lande, 1979Go, 1980Go; Lande and Arnold, 1983Go). In these equations z and P are, respectively, a column vector of phenotypic means and their variances and covariances (P is a matrix; P–1 is its inverse) before selection; asterisks denote these same parameters after selection within a generation. G is the inheritance matrix (additive genetic variance-covariance matrix) that translates the effects of phenotypic selection into genetic values, either across generations (1) or within a generation (2). These equations also show that the first and second derivatives of the adaptive landscape can be estimated from the changes that selection causes in the first and second moments of the phenotypic distribution. In other words, the slope of the adaptive landscape can be inferred from the shift that selection causes in the phenotypic mean. The curvature and orientation of the adaptive landscape can be inferred from the expansions and contractions that selection causes in phenotypic variances and covariances. Quadratic regression is a useful tool for accomplishing these inferences (Lande and Arnold, 1983Go; Phillips and Arnold, 1989Go; Brodie et al., 1995Go; Stinchcombe et al., 2002Go).

In our empirical quest for the adaptive landscape it will be useful to recognize an important way station known as the individual selection surface. Our approach will be to first characterize this individual selection surface and then use theoretically-derived connections between this surface and the adaptive landscape to make inferences about the adaptive landscape. What is the individual selection surface?

The individual selection surface (Pearson, 1903Go) is the relationship between the relative fitness of individuals in a population and the values of a set of phenotypic traits. The adaptive landscape is a function of mean trait values; the selection surface is a function of individual trait values (Table 1). Using data on the fitness of individuals and their phenotypic traits, we could use any of a variety of curve-fitting approaches to estimate the individual selection surface. A popular approach is to use the method of cubic splines, which is capable of matching the shape of a wide variety of surfaces in one or more character dimensions (Schluter, 1988Go; Schluter and Nychka, 1994Go). The goal in such an exercise is to accurately describe the shape of the selection surface. Our goal, however, is different and so demands another approach. Our goal is to characterize the selection surface in a way that will give us straightforward insights about the adaptive landscape. The cubic spline approach will not accomplish this goal because the connections between cubic splines and the adaptive landscape are unknown and may not exist. In other words, it has not been established that the first and second derivatives in equations (1) and (2) can be estimated with cubic splines. Although cubic splines will not get us to our goal, another approach will. That approach is called quadratic approximation. In this approach, second-order (quadratic) rather than third-order (cubic) polynomials are used to approximate the individual selection surface.


View this table:
[in this window]
[in a new window]
 
TABLE 1. Nomenclature for surfaces and landscapes

 
A quadratic approximation to the individual selection surface can be made by using a least-squares regression approach (Lande and Arnold, 1983Go). For convenience we will assume from here on that the values of individual values for traits, z, have been standardized so that their means are zero, and relative fitness, w, has been standardized so that its mean is one. We will use the following multiple regression model to relate the individual values of multiple traits, z, to relative fitness, w


{i1540-7063-043-03-0367-e3}

using the notation of Phillips and Arnold (1989)Go. In this multivariate version of the model, z is a column vector of individual traits values, zT is the row vector version of that same vector (the transpose of z). The intercepts in the model are {alpha} (a column vector) and errors or residuals are {varepsilon} (a column vector). The point of this model is to estimate the linear directional selection gradients ß (a column vector) and the nonlinear selection gradients {gamma} (a matrix). The factors of 1/2 in (3) serve the role of making the {gamma} coefficients equivalent to second derivatives.

Why are we interested in ß and {gamma}? One point of interest is that these coefficients correct for the effect of trait correlation and so can help us identify the targets of selection. We shall return to that point shortly. The other point of interest—the primary one in our present discussion—is that we can infer the first and second derivatives of the adaptive landscape from ß and {gamma}. The relationship of ß to the adaptive landscape is straightforward. ß is the average slope of the individual selection surface, w(z), but it is also the direction of steepest uphill slope on the adaptive landscape, ln W, in the vicinity of the phenotypic mean, z,


{i1540-7063-043-03-0367-e4}

where p(z) are the frequencies of the phenotypic values z in the population before selection (Lande and Arnold, 1983Go). The relationship of {gamma} to the adaptive landscape is slightly more complicated. In parallel with the first expression in (4), {gamma} is the average curvature of the individual selection surface,


{i1540-7063-043-03-0367-e5}

(Lande and Arnold, 1983Go). The curvature of the adaptive landscape, however, is a function of both {gamma} and ß,


{i1540-7063-043-03-0367-e6}

(Lande, 1979Go; Phillips and Arnold, 1989Go). One way to understand this more complicated relationship, is to notice that, whereas the elements in {gamma} represent just the effects of nonlinear selection on the variance or covariance of characters (Lande and Arnold, 1983Go, eq. 14a), the elements in the matrix {partial}2lnW/{partial}z2 represent both the effects of linear and nonlinear selection on variances and covariances. The matrix ßßT represents just the effects of linear selection on variances and covariances.


    SELECTION ON MORPHOLOGY VIA EFFECTS ON PERFORMANCE AND FITNESS
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
Performance is a measure of whole organism capacity (Bartholomew, 1958Go; Bennett, 1980Go). Depending on the traits and ecological circumstance of interest, performance might take the form of seed dispersal ability, shell-crushing capacity, swimming speed, or tree-climbing ability. Performance in its various guises is a central issue in ecological morphology, physiological ecology, functional morphology and other disciplines poised at the interface of phenotype and ecology (Wainwright and Reilly, 1994Go).

The relationship of performance to selection theory can be visualized with path diagrams that portray selection surfaces (Arnold, 1983Go). In such a diagram (Fig. 1) we see performance interposed between morphology and fitness. In the discussion that follow ‘morphology’ will be used as a shorthand for biochemical, physiological, behavioral and morphological attributes that affect performance. Turning to the regression equations, the interposition of performance means that the equation that relates morphology to fitness can be decomposed into two equations: one that relates morphology to performance and one that relates performance to fitness:


{i1540-7063-043-03-0367-e7}

where f denotes various measures of performance (a column vector) and fi denotes an element in that vector. The total selection gradient ß in the linear version of equation (3),


{i1540-7063-043-03-0367-e8}

is a product of the performance gradient ßf and the fitness gradient ßw (Arnold, 1983Go). Eq.(7) can be used to estimate ßf and ßw even if the performance surface or the fitness surface is curvilinear. In these instances of curvilinearity, however, we need a more complicated model than (7) to estimate curvature itself. We shall now turn to that more complicated model.



View larger version (9K):
[in this window]
[in a new window]
 
FIG. 1. Path diagram corresponding to the linear approximation (7) to the individual selection surface for a set of morphological traits (morphology), denoted z. The effects of morphology on fitness, w, are mediated through a set of performance variables, f. The performance gradient, ßf, is the slope of the fitted performance surface. The fitness gradient is ßw

 
If the performance surface is curved we can approximate its curvature with a quadratic function. In this quadratic approximation, performance is a function of quadratic variables, zzT, which represent the squares and products of trait values, as well as the trait values themselves,


{i1540-7063-043-03-0367-e9}

{gamma}f is a matrix of nonlinear performance gradients. We can also visualize this equation with a path diagram (Fig. 2). In the two trait case,


{i1540-7063-043-03-0367-e10}

where {gamma}11 and {gamma}22 describe curvature in character dimensions z1 and z2, respectively, and {gamma}12 describes the orientation of the surface. Quadratic surfaces can take a limited variety of shapes in two dimensions. The range of possibilities is shown in Figure 3. In three trait dimensions the range is considerably greater (Fig. 4).



View larger version (11K):
[in this window]
[in a new window]
 
FIG. 2. Path diagram corresponding to the quadratic approximation (9) to the individual selection surface for morphology. The nonlinear performance gradient {gamma}f is the curvature of the fitted surface. Other conventions as in Figure 1

 


View larger version (72K):
[in this window]
[in a new window]
 
FIG. 3. Various shapes for quadratic performance surfaces for two traits, z1 and z2. High and low points on the surface are indicated by + and –, respectively. Dotted lines indicate performance lines of least resistance (eigenvectors or principal axes). (a) A surface with a peak. (b) A surface with a saddle. (c) A surface with a ridge. (d) A surface with a rising ridge. These surfaces differ only in their coefficients for correlational performance, {gamma}f12. From Phillips and Arnold (1989)Go reproduced with permission

 


View larger version (34K):
[in this window]
[in a new window]
 
FIG. 4. Various shapes for quadratic performance surfaces for three traits (y1, y2, y3) plotted on their principal axes (dotted lines). Each figure is a contour representation of a four-dimensional surface. Performance is the fourth axis and is represented by a surface of higher performance nested within a surface of lower performance. Dotted lines represent performance lines of least resistance. Axis labels are shown in (a). The regions of highest performance are the origin in (b) and (c), a line in (d) and (f), a plane in (e) and a curved surface in (g) and (h). From Phillips and Arnold (1989)Go, reproduced with permission

 
What is the relationship of this quadratic surface to the actual performance surface, which might—after all—have a more complicated shape? Regardless of the shape of the actual performance surface, ßf will estimate its average slope and {gamma}f will estimate its average curvature (Fig. 5). ßf will also estimate its slope of the performance landscape in the vicinity of the population mean, z, but in general the curvature of the performance landscape at that point is {gamma}f–ß (Fig. 6). Note that this last formula is a the general one for the curvature of the performance landscape, because in those instances in which the permformance means lies directly under the optimum of the landscape (Fig. 5), ßf = 0 and so ß = 0 (Phillips and Arnold, 1989Go).



View larger version (13K):
[in this window]
[in a new window]
 
FIG. 5. The individual performance surface (lower panel) and the performance landscape (upper panel) for a single morpological trait. In the lower panel, the bottom solid curve represent the distribution, p(z), of the morphological variable z; the upper solid curve represents the individual performance surface, f(z); the dashed curve represents the quadratic approximation to the individual performance surface. In the upper panel, loge mean performance, ln f, is a function of mean morphology, z, in the population. The vertical dotted line marks the position of z. Note that the mean lies at the optimum of the performance landscape. After Phillips and Arnold (1989)Go

 


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 6. The individual performance surface (lower panel) and the performance landscape (upper panel) for a single morpological trait when the mean morphology does not lie at the optimum of the performance landscape. In this situation, the curvature of the quadratic approximation to the individual performance surface, {gamma}f, evaluated at the morpological mean, overestimates the curvature of the performance landscape, {gamma}f–ß. Conventions as in Figure 3. After Phillips and Arnold (1989)Go

 
It is useful to characterize the performance surface in terms of its principal axes, which represent performance lines of least resistance. The eigenvectors (or principal components) of the performance surface can be calculated from the matrix {gamma}f (Phillips and Arnold, 1989Go). The leading eigenvector of matrix {gamma}f represents the performance line of least resistance. Of all the lines that we might draw on the performance surface, movement along this line causes the smallest change in performance. The next best orthogonal line in this sense of least change in performance is the second eigenvector. The first and second eigenvectors of quadratic surfaces in Figure 3 are shown with dotted lines. To calculate the performance line of least resistance for the performance landscape we need to take the leading eigenvalue of the matrix {gamma}f–ß.

What about the total nonlinear selection gradient for a trait or a pair of traits? If we assume that the relationship between performance and fitness can be approximated with a linear fitness gradient, viz.,

(11)

then the total individual selection surface can be approximated by


{i1540-7063-043-03-0367-e12}

Notice that the total nonlinear selection gradient for the traits is a function of the linear fitness gradient, ßw, and the nonlinear performance gradient, {gamma}f. Using (3–6), the curvature of the adaptive landscape for the traits in the vicinity of the population mean is ß ({gamma}f–ß < 1, we can conclude that the curvature of the adaptive landscape will be less than the curvature of the performance landscape (Fig. 7).



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 7. The individual performance surface (lower panel) and the performance landscape (middle panel), and the adaptive landscape (upper panel) for a single morpological trait when the mean morphology lies at the optimum of the performance and adaptive landscapes. In the upper panel, loge mean fitness, ln W, is a function of mean morphology, z. Other conventions as in Figure 5

 

    SELECTION ON PERFORMANCE
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
In the preceding discussion we focused on morphology and how selection on morphology was mediated by its effects on performance. Our perspective was to view fitness from the perspective of fitness, with performance an intervening variable (Fig. 1). To consider selection on performance, we need a new perspective. We need to view fitness from the perspective of performance (Fig. 8). Turning to the issue of how selection acts on performance, we have two surfaces to consider: the individual selection surface for performance (or the selection surface for performance), in which individual fitness is a function of individual values for performance (Fig. 9, bottom), and the adaptive landscape for performance (or the performance adaptive landscape), in which the average value of fitness is a function of the population's average values for performance (Fig. 9, top). We can estimate the coefficients of the individual selection surface for performance using the quadratic approximation


{i1540-7063-043-03-0367-e13}

This model is the same as (3); we are merely considering performance, f, to be the trait of interest. To fit this model we would need data on performance and fitness for a series of individuals in a population (Arnold, 1983Go, 1986Go). ßw is the average slope of the individual selection surface for performance, as well as the slope of the adaptive landscape for performance in the vicinity of average performance,


{i1540-7063-043-03-0367-e14}

{gamma}w is the average curvature of the individual selection surface for performance,


{i1540-7063-043-03-0367-e15}

The curvature of the adaptive landscape for performance (Fig. 9, top) in the vicinity of mean performance is


{i1540-7063-043-03-0367-e16}



View larger version (15K):
[in this window]
[in a new window]
 
FIG. 8. Path diagram representing the quadratic approximation (13) to the individual selection surface for performance

 


View larger version (13K):
[in this window]
[in a new window]
 
FIG. 9. The individual selection surface (lower panel) and the adaptive landscape (upper panel) for a single performance variable. In the lower panel, the bottom solid curve represent the distribution, p( f), of the performance variable f; the upper solid curve represents the individual selection surface, w( f); the dashed curve represents the quadratic approximation to the individual selection surface. In the upper panel, loge mean fitness, ln W, is a function of mean morphology, f, in the population. The vertical dashed line marks the position of f. Note that when the mean lies at the optimum of the adaptive landscape, as in this portrayal, there is no directional selection on performance. After Phillips and Arnold (1989)Go

 

    EVOLUTION OF THE PERFORMANCE SURFACE AND THE ADAPTIVE LANDSCAPE
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
It is the adaptive landscape, in conjunction with quantitative inheritance, that governs trait evolution, eq. (1). In the short term, we can consider the adaptive landscape to be a fixed entity, with constant slope and curvature, but in the long term we need to consider the possibility that the landscape itself might evolve. That evolution might reflect change in abiotic conditions, relationships with other species or trait interactions. Before we consider landscape evolution, however, we need to consider our vision of the landscape itself.

The most tractable vision of the adaptive landscape is one that is restricted to the immediate vicinity of a population's mean (Arnold et al., 2001Go). This local vision is a practical one, because—as we have seen—we can use the tool of quadratic approximation to characterize the slope, curvature and orientation of the landscape in the vicinity of z or f. We can have local visions of multiple landscapes. To visualize multiple populations or species, we superimpose their landscapes and view them all at once in a space of average trait or performance values (Arnold et al., 2001Go, Fig. 9). The textbook vision of an adaptive landscape differs from this local vision of individual or multiple landscapes. A complex, rugged topography is often depicted in textbooks. Arnold et al. (2001)Go called it the global view of the landscape, for it portrays features near and far from the population mean. Dawkins (1996)Go has called it Mt. Improbable. The quandary of a population stuck on a local summit, down slope from the summit of Mt. Improbable, is a carry over from Wright's (1931Go, 1932Go) landscape in a space of gene frequencies. It is not clear, however, that a rugged topography, and the quandries resulting from local peaks, apply in the case of phenotypic characters. More importantly, we have no device that can tell us about landscape features far from the population's average. Mt. Improbable is, after all, Mt. Unknowable. This is the sense in which the local vision of the landscape is the most tractable. We have devices to give us a local vision. But if we adopt a local vision, how can the population ever evolve on the landscape? The answer is that the landscape's peak can move in evolutionary time (Simpson, 1944Go), causing the average trait or performance values to track that movement in response to directional selection.

A variety of models have been proposed for the evolution of the adaptive landscape (Arnold et al., 2001Go; Hansen and Martins, 1996Go). The simplest models all share the feature that the curvature and orientation of the landscape remains constant. Only the position of the landscape's peak evolves. With these provisos, we can recognize four broad categories of possibilities: (a) the peak undergoes random (e.g., Brownian) motion, (b) peak movement is selectively constrained about a fixed or moving point, (c) the peak moves at a consistent rate and in a particular, characteristic direction, (d) peak movement is episodic rather than steady. Obviously, there are other possibilities and even the ones listed could be combined in various ways. I want to focus on just one possibility because of its apparent relevance to the comparative study of performance. I want to consider the issue of whether the peak moves in ways that can be deduced by studying performance within populations. Is the movement of the adaptive landscape selectively constrained?

The position of the adaptive landscape may tend to evolve along performance lines of least resistance. The basis of this prediction is the proposition that the performance lines of least resistance arise from fundamental biomechanical relationships. The major features of the performance surface distill these fundamental relationships within populations. In a study of garter snakes, for example, Arnold and Bennett (1988)Go found that the performance surface that relates crawling performance to body and tail vertebral numbers had a major axis that promoted a positive correlation between these two parts of the vertebral column. Thus, the snakes with the fastest crawling speeds tended to lie close to the performance line of least resistance. A change in vertebral numbers that moved a hypothetical snake along the line had relatively little effect on performance; a change that moved a snake off the line lowered performance. An unknown biomechanical law that gives a performance premium to well proportioned snakes presumably underlies these results. Even though the biomechanical law is unknown in this case, it's domain of operation might transcend the limits of vertebral variation in a single population. Our proposal is then that the peak of the performance landscape tends to evolve along the performance line of least resistance. Changed ecological circumstances cause the position of the adaptive landscape to change, but that evolutionary movement is selectively constrained by the same biomechanical laws that prevail within populations. A particularly evocative prospect is presented by the model of trivariate stabilizing selection shown in Figure 4f. Movement of this surface (or its corresponding landscape) along the selective line of least resistance would correspond to concentric series of pipes sliding simultaneously along their long axes. This proposal is different from Schluter's (1996)Go proposition that evolution occurs along genetic lines of least resistance. Schluter's lines correspond to the principal axes of the G-matrix in (1–2). See Arnold et al. (2001)Go for a discussion and comparison of the two proposals.

Can we explain adaptive radiations by movement of adaptive (or performance) peaks along a selective line of least resistance? In other words, is this model capable of generating the elongate clusters of population or species means that typically characterize interspecific allometry or morphospace? One point that emerges from models of peak movement is that movement along a selective line of least resistance probably must be sustained rather than random to explain the most common kind of patterns in trait space. For example, consider a Brownian motion model in which trait evolution is constrained in the same multivariate pattern that constrains performance-based selection within populations. Because movement of the peak is about a fixed point, the long term expectation is for no correlation in species or population means. For simplicity, assume that the performance surface is the same among all descendant populations in a radiation and takes the form of a Gaussian surface so that


{i1540-7063-043-03-0367-eq1}

where {omega}f is a matrix of coefficients, analogous to variances and covariances, that describes the shape and orientation of the performance surface, and {theta}f is the peak of that surface. Suppose that the peak of the surface undergoes selectively-constrained Brownian movement. The peak moves as a normally distributed random variable with mean 0 and a variance-covariance matrix k{omega}f, where k is a scalar constant of proportionality (Arnold et al., 2001Go). Each population tends to track its moving peak, and so, after some elapsed number of generations, the population means have a dispersion pattern that reflects the shape and orientation of the performance surface. Using results from Lande (1979)Go and Hansen and Martins (1996)Go, the covariance of population means of the populations after t generations is

where Q(t) = exp[–tG({omega}f + P)–1]. Under some conditions this model predicts correlation in population means. Such patterns of covariance will be aligned with the major axes of the {omega}-matrix, which can in turn be estimated by the major axes of the {gamma}-matrix. These patterns of population covariance, however, decay very rapidly. Under plausible sets of parameter values, the covariance in population means decays to virtually nothing in less than 10,000 generations. Although the peak of the performance surface undergoes Brownian motion in a pattern that reflects the shape and orientation of the selection surface, the long term expectation is for population means to be close to the original position of the peak. In other words, Brownian motion seems unlikely to produce elongate patterns in population or species means. Our exploration of this model suggests that the position of the peak must move along selective lines of least resistance, rather than hover about a fixed point, to produce a radiation that reflects the within-population performance surface.

Comparative data can be used to test for selectively-constrained movement of the adaptive landscape. In the ideal case, we would have estimates of the performance surface from different populations or species. With those data we could test for homogeneity among populations in the {gamma}-matrix. If variation in {gamma} were not too great, we could estimate the pattern of dispersion of optima to see if that pattern coincides with the average line of selective least resistance within populations (i.e., with the average leading eigenvector of the {gamma}-matrix). Such an ideal case might require estimates of ten or more performance surfaces to accurately describe the multivariate dispersion of optima. If only one estimate of the performance surface is available, we might assume that mean of each population is close to its performance optimum, and then use the dispersion of population means as a surrogate for the dispersion of optima. Comparing the leading eigenvectors of the among-population variance-covariance matrix and the performance surface ({gamma}-matrix) provides a test for evolution along selectively-constrained lines of least resistance. The garter snake test case was of this kind and yielded remarkably strong correspondence between the leading eigenvectors of the among-population dispersion matrix (n = 9 populations) and the {gamma}-matrix for crawling speed; the vector correlation was 0.995. Thus, vertebral numbers in garter snakes evolve along performance lines of least resistance.

The perspective that the performance surface is curved and related to the adaptive landscape highlights two unresolved, empirical issues: (a) What is the time scale for the evolution of the curvature and orientation of the performance surface? The usefulness of the performance surface in understanding the evolution of morphology depends on the evolutionary stability of the surface. Because the surface springs from underlying components of the organism and their interactions, it seems likely that the performance surface and morphology evolve on the same time scale. Comparative studies of surfaces could test this proposition directly. (b) What are appropriate models for the evolution of the performance optimum? Evolution of the optimum along performance lines of least resistance is obviously just one possibility. Nevertheless, this proposal represents a simple hypothesis that cuts through the tangle of other possibilities (Emerson and Arnold, 1989Go). Here again, comparative studies of surfaces are a vital step in validating our models of surface and landscape evolution.


    ACKNOWLEDGMENTS
 
I thank J. G. Kingsolver and R. B. Huey for organizing the symposium and providing me with an opportunity to revisit a model after a lapse of 20 years. I am grateful to the Evolutionary Morphology Group at the University of Chicago and to A. F. Bennett and R. B. Huey for discussions, more than 20 years ago, that precipitated the morphology-performance-fitness idea. A. G. Jones, R. Lande, M. E. Pfrender, and P. C. Phillips made me think more precisely about selection, surfaces and landscapes. I thank Joel Kingsolver and an anonymous reviewer for comments on the manuscript. The preparation of this article was supported in part by NSF grants DEB-9903934 and IRCEB-0110666.


    FOOTNOTES
 
1 From the Symposium Selection and the Evolution of Performance in Nature presented at the Annual Meeting of the Society for Integrative and Comparative Biology, 4–8 January 2003, at Toronto, Canada. Back

2 E-mail: arnolds{at}science.oregonstate.edu Back


    References
 TOP
 SYNOPSIS
 INTRODUCTION
 THE ADAPTIVE LANDSCAPE FOR...
 SELECTION ON MORPHOLOGY VIA...
 SELECTION ON PERFORMANCE
 EVOLUTION OF THE PERFORMANCE...
 References
 
Arnold, S. J. 1983. Morphology, performance and fitness. Amer. Zool, 23:347-361.

Arnold, S. J. 1986. Laboratory and field approaches to the study of adaptation. In M. Feder and G. Lauder (eds.), Predator-prey relationships, pp. 157–179. University of Chicago Press, Chicago.

Arnold, S. J., and A. F. Bennett. 1988. Behavioural variation in natural populations. V. Morphological correlates of locomotion in the garter snake Thamnophis radix. Biol. J. of the Linnean Society, 34:175-190.

Arnold, S. J., M. E. Pfrender, and A. G. Jones. 2001. The adaptive landscape as a conceptual bridge between micro- and macroevolution. Genetica, 112:/1139-32.

Bartholomew, G. A. 1958. The role of physiology in the distribution of terrestrial vertebrates. In C. L. Hubbs (ed.), Zoogeography, pp. 81–95. Amer. Soc. Adv. Sci., Washington, DC.

Bennett, A. F. 1980. The thermal dependence of lizard behaviour. Animal Behaviour, 28:752-762.[CrossRef]

Bornberg-Bauer, E., and H. S. Chan. 1999. Modeling evolutionary landscapes: Mutational stability, topography, and superfunnels in sequence space. Proc. Natl. Acad. Sci. U.S.A, 96:10689-10694.[Abstract/Free Full Text]

Brodie, E. D., III, A. J. Moore, and F. J. Janzen. 1995. Visualizing and quantifying natural selection. Trends in Ecol. Evol, 10:313-318.[CrossRef]

Dawkins, R. 1996. Climbing mount improbable. W. W. Norton and Co., New York, New York.

Dobzhansky, Th. 1937. Genetics and the origin of species. Columbia Univ. Press, New York, New York.

Eldredge, N., and J. Cracraft. 1980. Phylogenetic patterns and the evolutionary process. Columbia Univ. Press, New York, New York.

Emerson, S. B., and S. J. Arnold. 1989. Intra- and interspecific relationships between morphology, performance, and fitness. In D. B. Wake and G. Roth (eds.), Complex organismal functions: Integration and evolution in vertebrates, pp. 295–314. Wiley, New York.

Fisher, R. A. 1930. The genetical theory of natural selection. Clarendon, Oxford, UK.

Fontana, W., and P. Schuster. 1998. Continuity in evolution: On the nature of transitions. Science, 280:1451-1455.[Abstract/Free Full Text]

Govindarajan, S., and R. A. Goldstein. 1997. The foldability landscape of model proteins. Biopolymers, 42:427-438.[CrossRef][Web of Science][Medline]

Hansen, T. F., and E. P. Martins. 1996. Translating between microevolutionary process and macroevolutionary patterns: The correlation structure of interspecific data. Evolution, 50:1404-1417.[CrossRef]

Lande, R. 1976. Natural selection and random genetic drift in phenotypic evolution. Evolution, 30:314-334.[CrossRef][Web of Science]

Lande, R. 1979. Quantitative genetic analysis of multivariate evolution, applied to brain-body size allometry. Evolution, 33:402-416.[CrossRef][Web of Science]

Lande, R. 1980. The genetic covariance between characters maintained by pleiotropic mutation. Genetics, 94:203-215.[Abstract/Free Full Text]

Lande, R., and S. J. Arnold. 1983. The measurement of selection on correlated characters. Evolution, 37:1210-1226.[CrossRef][Web of Science]

Pearson, K. 1903. Mathematical contributions to the theory of evolution. XI. On the influence of natural selection on the variability and correlations of organs. Phil. Trans. Roy. Soc. London A, 200:1-66.

Phillips, P. C., and S. J. Arnold. 1989. Visualizing multivariate selection. Evolution, 43:1209-1222.[CrossRef]

Provine, W. B. 1986. Sewall Wright and evolutionary biology. University of Chicago Press, Chicago, Illinois.

Schluter, D. 1988. Estimating the form of natural selection on a quantitative trait. Evolution, 42:849-861.[CrossRef][Web of Science]

Schluter, D. 1996. Adaptive radiation along genetic lines of least resistance. Evolution, 50:1766-1774.[CrossRef][Web of Science]

Schluter, D., and D. Nychka. 1994. Exploring fitness surfaces. Amer. Nat, 143:597-616.[CrossRef]

Schmalhausen, I. I. 1949. Factors of evolution, the theory of stabilizing selection. Univ. Chicago Press, Chicago, Illinois.

Simpson, G. G. 1944. Tempo and mode in evolution. Columbia Univ. Press, New York, New York.

Simpson, G. G. 1953. The major features of evolution. Columbia University Press, New York.

Stinchcombe, J. R., M. T. Rutter, D. S. Burdick, P. Tiffin, M. D. Rausher, and R. Mauricio. 2002. Testing for environmentally induced bias in phenotypic estimates of natural selection. Amer. Natur, 160:511-523.[CrossRef]

Wainwright, P. C., and S. M. Reilly. 1994. Ecological morphology. Univ. Chicago Press, Chicago, Illinois.

Wright, S. 1931. Evolution in Mendelian populations. Genetics, 16:97-159.[Free Full Text]

Wright, S. 1932. The roles of mutation, inbreeding, crossbreeding, and selection in evolution. Proc. VI Internat. Cong. of Genet, 1:356-366.

Wright, S. 1945. Tempo and mode in evolution: A critical review. Ecology, 26:415-419.[CrossRef]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
GeneticsHome page
J. Hunt, M. W. Blows, F. Zajitschek, M. D. Jennions, and R. Brooks
Reconciling Strong Stabilizing Selection with the Maintenance of Genetic Variation in a Natural Population of Black Field Crickets (Teleogryllus commodus)
Genetics, October 1, 2007; 177(2): 875 - 880.
[Abstract] [Full Text] [PDF]


Home page
J. Exp. Biol.Home page
M. Nikinmaa and W. Waser
Molecular and cellular studies in evolutionary physiology of natural vertebrate populations: influences of individual variation and genetic components on sampling and measurements
J. Exp. Biol., June 1, 2007; 210(11): 1847 - 1857.
[Abstract] [Full Text] [PDF]


Home page
J. Exp. Biol.Home page
J. H. Long Jr, T. J. Koob, K. Irving, K. Combie, V. Engel, N. Livingston, A. Lammert, and J. Schumacher
Biomimetic evolutionary analysis: testing the adaptive value of vertebrate tail stiffness in autonomous swimming robots
J. Exp. Biol., December 1, 2006; 209(23): 4732 - 4746.
[Abstract] [Full Text] [PDF]


Home page
Integr. Comp. Biol.Home page
J. G. Swallow, J. S. Rhodes, and T. Garland Jr.
Phenotypic and Evolutionary Plasticity of Organ Masses in Response to Voluntary Exercise in House Mice
Integr. Comp. Biol., June 1, 2005; 45(3): 426 - 437.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (22)
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Arnold, S. J.
Right arrow Search for Related Content
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?