- Research article
- Open Access
Differences in the selection response of serially repeated color pattern characters: Standing variation, development, and evolution
BMC Evolutionary Biologyvolume 8, Article number: 94 (2008)
There is spectacular morphological diversity in nature but lineages typically display a limited range of phenotypes. Because developmental processes generate the phenotypic variation that fuels natural selection, they are a likely source of evolutionary biases, facilitating some changes and limiting others. Although shifts in developmental regulation are associated with morphological differences between taxa, it is unclear how underlying mechanisms affect the rate and direction of evolutionary change within populations under selection.
Here we focus on two ecologically relevant features of butterfly wing color patterns, eyespot size and color composition, which are similarly and strongly correlated across the serially repeated eyespots. Though these two characters show similar patterns of standing variation and covariation within a population, they differ in key features of their underlying development. We targeted pairs of eyespots with artificial selection for coordinated (concerted selection) versus independent (antagonistic selection) change in their color composition and size and compared evolutionary responses of the two color pattern characters.
The two characters respond to selection in strikingly different ways despite initially similar patterns of variation in all directions present in the starting population. Size (determined by local properties of a diffusing inductive signal) evolves flexibly in all selected directions. However, color composition (determined by a tissue-level response to the signal concentration gradient) evolves only in the direction of coordinated change. There was no independent evolutionary change in the color composition of two eyespots in response to antagonistic selection. Moreover, these differences in the directions of short-term evolutionary change in eyespot size and color composition within a single species are consistent with the observed wing pattern diversity in the genus.
Both characters respond rapidly to selection for coordinated change, but there are striking differences in their response to selection for antagonistic, independent change across eyespots. While many additional factors may contribute to both short- and long-term evolutionary response, we argue that the compartmentalization of developmental processes can influence the diversification of serial repeats such as butterfly eyespots, even under strong selection.
Despite the diversity of animal form in nature, a limited range of phenotypes is often observed [1–3]. Stasis, convergence, and limits in the direction of morphological evolution occur within lineages even when strong selection is implicated in divergence [4–7]. The propensity for phenotypes to evolve, or evolvability [8–11], is determined by the capacity of genetic and developmental systems to produce heritable variation and the action of selection on those variants [12, 13]. Consequently, rates of adaptive change and diversification are expected to depend critically on underlying developmental mechanisms [2, 5, 9, 14–16]. Several different lines of evidence provide more understanding of the role development plays in determining the direction of evolutionary change, including the documented association between morphological differences and changes in regulation of developmental-genetic pathways (reviewed in ), analysis of the variational properties of different developmental mechanisms [18–20], and analyses of the distribution of morphological characters in a given phylogeny . How development affects the rate and direction of evolutionary change in contemporary populations under selection, however, remains largely unexplored .
Diversification of serial repeats
The elaboration and differentiation of serially repeated characters provides clear examples of the association between underlying mechanisms and morphological evolution. The complexity and diversity of arthropod body segments, tetrapod limbs, vertebrate teeth, and nymphalid butterfly wing color patterns all result from the evolutionary differentiation of homologous serial repeats [23–26]. Individual elements often covary strongly, reflecting the effects of shared developmental pathways and mechanisms coordinating their development . This covariance can affect rates of evolutionary change [5, 15, 27], potentially limiting the independent evolution of serial repeats. As a consequence, development can bias evolution (sensu ), as certain morphological changes may be more readily achieved than others. However, recent empirical studies demonstrate that strong limits predicted by patterns of covariation do not always prevent the independent evolution of morphological traits under artificial selection [29–31]. One plausible explanation for the discrepancy between theoretical predictions and empirical results is the failure of current evolutionary genetic models to account for the details of development that determine how phenotypic variation is generated, and thereby affect the rate and direction of evolutionary change [32–34]. Artificial selection experiments are a powerful tool for exploring limits on the process of morphological evolution, particularly in the evolution of correlated characters [35–43]. Here we focus on two characters with similar patterns of covariation among a set of serial repeats. These characters, elements of the color pattern on butterfly wings, differ in key aspects of their development. We use artificial selection to explore effects of these differences on the rate and direction of evolutionary change in this complex morphological phenotype.
Evolution of eyespot patterns on butterfly wings
Within Bicyclus and other closely related butterfly genera, wing pattern diversification was probably driven in part by a long history of sexual and natural selection acting on the eyespots [44, 45], which are serially repeated along the wing margins. Each eyespot consists of concentric rings of black and gold color surrounding a white central spot . Color composition (the relative dimensions of the gold and black rings) and size are two important functional characteristics of eyespots [47, 48] that differ in critical aspects of their development. Properties of an inductive signal (e.g., signal concentration) produced by cells in the eyespot organizer, the central focus , largely determine eyespot size [46, 50, 51], whereas color composition is determined by the threshold response of surrounding tissue to the concentration gradient of the diffusing focal signal [52, 53]. Each of these two aspects of induction (signaling and response) appears to contribute differently to phenotypic variation in the eyespot pattern  and may have very different consequences for color pattern evolution.
When selection targets either the color composition or size of a single eyespot, all eyespots evolve in concert, even those on different wing surfaces [51, 53]. Despite differences in the underlying development of color composition and size, each character is strongly and positively correlated across serially repeated eyespots and their main axes of variation  in captive B. anynana populations lie in the direction of concerted changes in two or more eyespots. We used artificial selection to search for limits in this phenotype space, and selected for concerted (along the main axis of variation) and antagonistic (orthogonal to the main axis) change in the color composition or size of two eyespots. By applying similar selection to both characters and using a comparative analysis we test whether the characters differ in their propensity to evolve independently among serial repeats. Given initial similarities in patterns of standing variation for both characters, we argue in favor of the hypothesis that known differences in the developmental determination of eyespot color composition and eyespot size affect their ability to respond to antagonistic selection, and thus to evolve independently in a population under selection.
Results and Discussion
Differences in selection response between eyespot color composition and size
Color composition and size were similarly correlated between the serially repeated eyespots in the B. anynana base population before selection (color composition: rpearson ± SE = 0.42 ± 0.03, N = 1056, between ventral hindwing eyespots E4 and E6; size: rpearson ± SE = 0.52 ± 0.02, N = 2254, between anterior and posterior dorsal forewing eyespots). For both characters, most phenotypic variation in the base population lay in the direction of positive covariation between two eyespots but each exhibited substantial variation around this main axis (Fig. 1a, d). Previous experiments demonstrated the genetic basis for these phenotypic correlations and the presence of substantial additive genetic variation for the size and color of individual eyespots [51, 53, 55, 56]. As expected, both color composition and size responded rapidly to concerted selection imposed along the main axis of phenotypic variation (Fig. 1b, e; Fig. 2a, b). After 10 generations of concerted selection (see Materials and Methods), novel phenotypes well outside the range of variation in the base population were produced: color composition diverged ~2–3 SD (standard deviations, calculated from the base population) in the 'Black-Black' and 'Gold-Gold' directions relative to unselected control lines (Fig. 2a, c), and eyespot size diverged ~4–6 SD in the 'Large-Large' and 'Small-Small' directions relative to unselected controls (Fig. 2b, d). In contrast, responses to antagonistic selection demonstrate that the boundaries of phenotype space are very different for the two characters in this population (Fig. 1c, f; Fig. 2a, b). There is little potential for independent evolutionary change in eyespot color composition, compared with the strikingly flexible evolution of eyespot sizes.
After 10 generations of antagonistic selection, no novel 'uncoupled' eyespot phenotypes appeared in the 'Black-Gold' or 'Gold-Black' directions (Fig. 1c; Fig. 2a, c). The rate of response to antagonistic selection (scaled to the cumulative selection differential; Table 1) is clearly lower than the rate of response to selection for concerted changes in eyespot color composition. In some replicates there was no significant response to antagonistic selection for color composition (Table 1). Although selection can rapidly shift the mean color composition of the entire eyespot pattern (Fig. 1b; Fig. 2a, c), there are limits on the independent evolutionary change of individual eyespots. No comparable limits are apparent for eyespot size; antagonistic selection ('Large-Small' and 'Small-Large' directions; Fig. 1F, Fig. 2b, d) rapidly produced combinations of large and small eyespots that were not present in the starting population . Moreover, the rate of response to antagonistic selection on eyespot size was as high (or higher, in some replicates) as the rate of response to concerted selection (Table 1).
Potential causes of the difference in selection response
Our results show that under similar selection regimes, there are strong limits on the independent evolution of eyespot color composition that do not exist for eyespot size. This difference in evolutionary response occurred despite the fact that both characters exhibited substantial phenotypic variation in both concerted and antagonistic directions prior to selection (Fig. 1). Genetic correlations (measured in the outbred stock population; see Materials and Methods) across the targeted eyespots for color composition (rG ± SE = 0.68 ± 0.15) and size (rG ± SE = 0.58 ± 0.1) suggest that patterns of genetic covariation were also similar for both characters. The magnitude and direction of phenotypic and genetic correlations (and the relationship between the two; ) is generally taken to indicate the degree of limitation on the independent evolution of two or more characters [27, 59–62]. Our experimental approach involved investigating the realized potential for independent evolution among repeats of two characters with similar patterns of standing variation known to differ in important aspects of their developmental determination. We expected to observe qualitatively similar responses of color composition and size in all selected directions if the patterns of covariation alone were the major internal factors governing the independent evolution of eyespot characters. While we argue that the degree of compartmentalization of the developmental mechanisms differs between color composition and size, we discuss additional factors that may also explain the observed evolutionary responses.
Compartmentalization of focal signaling
Our results demonstrate much greater flexibility in the evolutionary response of eyespot sizes compared with color composition, despite the fact that both are characteristics of the same serially repeated color pattern elements on the wings of B. anynana. One explanation for this pattern is that the developmental pathways that regulate production of the focal signal (and determine eyespot size; [50, 51, 63]) are much more flexible and more easily decoupled across a wing surface than the pathways that regulate the threshold level of response to signal concentrations (and determine eyespot color composition; [52, 53, 64, 65]). Each eyespot-competent area along the wing margin produces its own eyespot organizer [51, 66]. This process appears to involve independent genetic control of organizer properties; single eyespots can be independently added to or eliminated from forewing or hindwing surfaces, without affecting the characteristics of eyespots in adjacent wing cells [67, 68]. Eyespots are deleted from wing surfaces when focal establishment fails in specific wing cells [67, 68], and eyespot-specific allelic effects are associated with changes in size of individual eyespots . The wing veins act to further compartmentalize the signaling foci across a wing surface . Compartmentalization (or modular organization; ) in the regulation of initially identical developmental programs is considered critical to the evolutionary divergence of serial repeats . Selection acting on heritable, independent variation in developmental regulation can promote evolutionary diversification, even when similarities in underlying developmental programs contribute to positive covariation among the characters .
Wing-wide responses to signal thresholds
In contrast, the observed tight coupling of color composition across eyespots may result from a lack of compartmentalization in the regulation of threshold levels of response to the concentration gradient of the diffusing focal signal. Cells in the region surrounding the eyespot focus exhibit a threshold response to local signal concentration, which induces the expression of regulatory genes that control pigment synthesis (Fig. 3a–f; [45, 52, 53]). Earlier experiments clearly demonstrate that the signaling and response components of eyespot determination can be selected independently, and that selection for eyespot color composition affects only the level of threshold response [51, 53]. In our experiment, concerted selection for eyespot color composition altered the response properties of the wing epidermis, even in regions of the wing that do not normally produce an eyespot and were not directly targeted by selection. The color composition of ectopic eyespots induced on wings of butterflies from concerted selection lines (Fig. 4) further suggest that the threshold response is not compartmentalized across comparable regions of a wing surface (spanning at least three wing cells) where there is clear compartmentalization of focal signaling [67, 68]. If threshold response is regulated across regions spanning several wing cells (or across an entire wing surface ), variation in the color composition of two eyespots (Fig. 1a) might result from random perturbations during the development of individual wing compartments, not from independent heritable variation in response properties. Although the developmental-genetic architecture of the processes that regulate eyespot color composition requires further characterization, we suggest that the extent to which focal signaling and response thresholds are compartmentalized is a critical factor determining the potential for evolutionary change in the wing color pattern.
Sources of covariation across wing surfaces
Selection targeted pairs of eyespots on two wing surfaces that differ in both morphology and ecological function. Eyespots on the dorsal  and ventral  wing surfaces are thought to function differently in B. anynana; however, all eyespots are determined by similar spatial and temporal patterns of gene expression [63, 72] regardless of their location. This shared development appears to play a primary role in shaping covariation and strongly integrating eyespots among all wing surfaces in Nymphalid butterflies [44, 73–76], although functional differentiation of some wing characters can also affect patterns of covariation . Single- [51, 53] and multi-trait selection on eyespot size  has demonstrated that the evolutionary potential of eyespots on all wing surfaces is strongly coupled and differs little among wing surfaces, despite differences in morphology or function. Antagonistic selection targeting the size of the two dorsal forewing eyespots also uncoupled individual eyespots on the ventral hindwing. The relative sizes of ventral hindwing eyespots E4 and E6 (the targets of selection for color composition) differed by ~2 SD in the antagonistic selection lines compared to unselected controls (; ANOVA for effect of selection direction: F4,495 = 67.7, P < 0.0001; Tukey's HSD for relative size of eyespots 4 and 6 in 'Large-Small' vs 'Small-Large' lines: t396 = 8.3, Padj < 0.001). Thus the ability to decouple eyespot sizes by selection is clearly not limited to the dorsal forewing. In contrast, across a similar spatial region of the ventral hindwing (two eyespots separated by an intervening wing cell), eyespot color composition is strongly coupled and antagonistic selection cannot produce independent evolutionary changes.
We suggest that known differences in the developmental determination of eyespot color composition and size contribute to qualitative differences in their responses to antagonistic selection. However, our experiments do not rule out the possibility that other factors contributed to these differences. In addition to developmental constraints or biases, factors that may limit the evolution of character combinations include the effects of past selection, the effects of inbreeding or genetic background, the action of maternal effects or gene-by-environment interactions, asymmetric gene frequencies or other differences in patterns of allelic effects, linkage disequilibrium or physical linkage, and correlation with other fitness-related characters (reviewed in [78, 79]). Regardless of the precise mechanism, differences in selection response are unlikely to have been caused by differences in the two experimental populations. Both starting populations were derived from the same outbred stock population, which was well adapted to the laboratory and maintained at high effective population size  before selection began, and both experiments were conducted under similar environmental conditions. Consequently, neither differences in allele frequencies between the starting populations, nor in the interactions between genetic background and environmental conditions across experiments  are likely to explain the qualitative differences in selection response between the two characters.
Differences in the history (form and/or intensity) of natural selection acting on eyespot color composition and eyespot size is one possible source of differences between the characters that could influence their responses to antagonistic selection. Selection favoring specific character combinations can build or contribute to phenotypic and genetic covariation [30, 36, 37, 82] via linkage disequilibrium. This type of covariation is expected to break down rapidly under antagonistic selection , and the breakdown of linkage disequilibrium created by past natural selection on eyespot size (for combinations of overall larger or smaller eyespots) may have contributed to the rapid independent evolution of eyespot sizes under antagonistic selection . Differences in past selection acting on color composition and size provides one possible mechanism for differences in the underlying basis of covariation between the two characters. As comparable patterns of covariation may be reached through many different developmental or genetic mechanisms [84, 85], this does not exclude our developmental hypothesis.
Our experiment explored whether limits exist on the evolution of particular types of eyespot combinations given standing variation in a population (Fig. 1a, d). Change in patterns of allelic interaction or allele frequency from those currently available in the population could affect the trait combinations or range of morphologies that can be reached via selection [34, 79, 84–86]. It is also possible that antagonistic selection acted on rare variants for eyespot size that, by chance, were not present for eyespot color composition in the stock population. Though our experiment does not rule out effects of several different underlying mechanisms, our results are consistent with theory suggesting that the mechanistic basis of character covariation is a more important predictor of evolutionary change than the correlation estimates alone [5, 32, 84, 85]. Because developmental interactions both shape patterns of phenotypic and genetic covariation among characters and evolve in response to selection acting on those characters [9, 87], genetic architecture, developmental mechanisms, and population-genetic processes interact in complex ways . Understanding how these factors interact to shape the evolution of complex traits, such as butterfly wing patterns, requires further detailed dissection of underlying mechanisms.
Diversification within the lineage
Phylogenetic patterns of morphological diversity are generally consistent with the idea that development limits the ability of species to occupy phenotype space (reviewed in ). In B. anynana, despite similar patterns of standing variation for eyespot color composition and size in the base population, responses to artificial antagonistic selection show that these two characteristics of the eyespot pattern do not evolve in the same way under similar antagonistic selection regimes. Moreover, these differences, found within a single B. anynana population, are consistent with wing pattern diversity across the genus (Fig. 3g–k; see also ). Among Bicyclus species, eyespot color composition seems to have evolved along a line of 'least resistance' [27, 28] parallel to the main axis of variation in B. anynana. In contrast, eyespot size has explored a larger area of potential phenotype space along the 'uncoupling' axis .
A number of studies suggest that natural selection plays a dominant role in morphological evolution and the distribution of lineages in phenotype space [29, 30, 36, 39, 57, 91]. However, a growing body of evidence shows that developmental mechanisms and individual components of developmental gene networks can differ in their effects on phenotypic variation [19, 20, 54], and presumably, bias or limit adaptive morphological evolution. Here we compared two characters that show similar patterns of standing variation in a population, but differ in an important aspect of their developmental determination. These characters evolve differently under artificial antagonistic selection, in a manner consistent with the developmental differences between them, though we cannot rule out the effects of other internal factors or processes besides, or in addition to, development. Butterfly color patterns and other complex forms we see in nature result from novel combinations of their individual elements. We suggest that compartmentalized development enables the production of novel phenotypes when selection favors independent changes in these elements. Under similar selection, developmental mechanisms that are less compartmentalized may limit the range of likely evolutionary outcomes.
In Bicyclus anynana, two eyespots are present on the dorsal forewing (anterior, A, and posterior, P) and seven on the ventral hindwing (eyespots E1–E7). In separate experiments, we selected for the relative size (the ratio of total eyespot diameter divided by the distance between the centers of A and P; ) or the color composition (the size of the black color ring relative to total eyespot diameter) of two eyespots simultaneously. Selection for eyespot size targeted the two dorsal forewing eyespots . The results of size selection showed that selection to uncouple A and P also uncoupled the relative sizes of anterior and posterior groups of ventral hindwing eyespots, and that the division between the two groups of ventral hindwing eyespots lay between E4 and E6 . Thus, we chose to select on the color composition of E4 and E6 because: 1) the size experiment showed that these two eyespots could be decoupled, even when they were targeted only indirectly by selection (see Results and Discussion); and 2) color composition is a function of total eyespot diameter, so it was necessary to choose two eyespots of initially similar size and color composition before selection.
The starting population for each experiment was derived from an outbred stock maintained in the laboratory for > 100 generations at high Ne . For each trait we established three types of lines from the stock population: 1) antagonistic selection lines where two eyespots were selected in opposite direction (large A and small P, 'Large-Small,' or the opposite combination 'Small-Large'; E4 gold and E6 black, 'Gold-Black,' or the opposite combination 'Black-Gold'); 2) concerted selection lines where two eyespots were selected in the same direction ('Gold-Gold,' or 'Black-Black'; 'Large-Large,' or 'Small-Small'); and 3) unselected control (UC) lines.
To estimate eyespot color composition, total diameter and the diameter of the inner black disc of E4 and E6 were measured in the proximal-distal axis, along the midline of each eyespot (parallel to the wing veins). Color composition of an eyespot was estimated from the orthogonal (Type II) regression of black disc diameter on total eyespot diameter. In this analysis, positive and negative residuals indicate relatively black and relatively gold eyespots, respectively. To identify individuals for selection, we performed an additional orthogonal regression of E4 color composition on E6 color composition (from the previous regression). The regression line is similar to the axis of the first principal component. For concerted selection, we ranked the residuals parallel to the regression line (similar to individual loadings on the first principal component); extreme positive or negative residuals indicated that both eyespots were relatively black or relatively gold, respectively. For antagonistic selection, we ranked the residuals perpendicular to the regression line (similar to individual loadings on the second, orthogonal, principal component axis); extreme positive residuals indicated individuals where E4 was relatively black compared to E6, and extreme negative residuals indicated individuals where E6 was relatively black compared to E4. Selection for eyespot size was made on the basis of the additive combination of the ranks for each individual eyespot as explained in . There were two replicates of each selection and each UC line, with one additional UC replicate in the size experiment. To estimate correlations between eyespot characters in the base population, we calculated Pearson product-moment correlations between the target eyespots in the starting population for each experiment.
Female butterflies were selected for 10 generations. For the color composition experiment, 1056 females were measured and randomly split into two groups and 40 females were randomly drawn from each group to establish the UC lines. Within each group, the remaining females were ranked and the 40 females with the most extreme phenotypes in each direction were selected, so that each selection line was replicated twice in total. For the size experiment, 2254 females were measured and randomly split into two groups after 40 females were drawn to establish the first UC line . All selected females were placed with ~50 randomly chosen males and allowed to lay eggs. Rearing conditions for the stock population and selected lines have been described previously . We maintained similar selection intensities in both experiments for generations 1–10: each line, per generation, we measured 140–240 females for color composition (mean ± SE: 209 ± 5) and 15–200 females for size (mean ± SE: 173 ± 3). We selected 40 females per line every generation in the color experiment; in the size experiment the number decreased to 35 females between generations 5–10 .
Each generation, trait values for each replicate line were calculated relative to their respective UC lines. A least-squares regression of selection response on cumulative selection differential was used to estimate the rate of response to selection for each replicate. We used nested analysis of covariance (ANCOVA) to test for differences in the rate of response between selection directions, with 'cumulative selection differential' as the covariate and 'replicate within direction' as a random effect. Tukey's multiple comparisons were used to compare responses of individual replicates when there was a significant interaction between the two main factors in the nested ANCOVA. Analyses were conducted using PROC GLM in SAS version 6.12 (SAS Institute, 1996).
Induction of ectopic eyespots
Damage to pupal wings in the first 18 hours after pupation induces formation of ectopic eyespots on the dorsal forewing  and can be used to assay the threshold response of wing tissue to eyespot-inducing signals [46, 53]. To determine the effect of selection for eyespot color composition on this threshold response, we damaged the pupal wings of females in each of the Black-Black, Gold-Gold, and UC lines and measured the color composition of the resulting ectopic eyespots. The experiment was conducted during generation 14 (lines were maintained under random mating conditions after generation 10, with 40 females and 55 males selected at random and ~250 offspring reared per line per generation).
To reduce the size of native eyespots (A and P) on the dorsal forewing and increase visibility of the induced ectopic, we first pierced the foci of the anterior and posterior eyespots on the left dorsal forewing of each individual 4.5 hours after pupation with a finely sharpened tungsten needle . We then induced ectopic eyespots by piercing each of the left forewings at a site in the fourth wing cell (immediately distal to the normal location of the eyespots) 17 hours after pupation. We returned operated pupae to 27°C, froze newly emerged adults (after their wings had fully hardened), and measured the total diameter and inner black diameter of ectopic eyespots and the interpupil distance between the reduced A and P eyespots on each manipulated (left) wing. We also measured the total diameter and black disc diameter of P on the right dorsal forewing, and the interpupil distance on this unmanipulated (right) wing.
We used ANCOVA to test for differences in color composition of ectopic eyespots among Black-Black, Gold-Gold, and UC lines. With 'total ectopic diameter' as the covariate, we examined the fixed effect of 'selection direction (control, black, or gold)' and the random effect of 'replicate line nested within direction' on the size of the ectopic black disc. None of the interaction terms (selection direction by covariate and random effect by covariate) were significant and we used Tukey's HSD to test for differences in the adjusted means between Black-Black, Gold-Gold, and UC lines.
Estimates of genetic covariation in the stock population
We used a paternal half-sib breeding design  to estimate quantitative genetic parameters in our outbred stock population. We randomly selected 100 virgin males from the stock population at adult eclosion; each male was mated sequentially to 2 virgin females. At hatching, ~30 eggs per female were transferred to mesh rearing cages and fed on young maize plants ad libitum until pupation. Full-sib offspring were reared together but densities were kept low to minimize interaction and competition between individuals. Rearing cages were moved every 4 days to randomize environmental effects within the growth chamber. Emerging adult offspring were allowed several hours for their wings to expand and fully harden before being frozen for later analysis.
Five female offspring were randomly selected from each of 174 full-sib families (representing 87 sires who successfully produced offspring by two dams each) and dorsal forewing eyespots A and P, and ventral hindwing eyespots E4 and E6 were measured as described above. We used our nested breeding design (full-sib dam families nested within sire families) to estimate sire, dam, and progeny variance and covariance components in SAS 6.12. Genetic correlations between eyespot pairs (and their standard errors) were calculated according to ; to eliminate potential inflation of the estimates by maternal effects, genetic correlations were calculated from the among-sire variance components only .
Muller GB, Wagner GP: Novelty in Evolution - Restructuring the Concept. Annual Review of Ecology and Systematics. 1991, 22: 229-256. 10.1146/annurev.es.22.110191.001305.
Alberch P: Ontogenesis and morphological diversification. American Zoologist. 1980, 20 (4): 653-667.
McGhee GR: The Geometry of Evolution: Adaptive Landscapes and Theoretical Morphospaces. 2007, Cambridge , Cambridge University Press
Alberch P: Developmental constraints: why St. Bernards often have an extra digit and poodles never do. American Naturalist. 1985, 126 (3): 430-433. 10.1086/284428.
Blows MW, Hoffmann AA: A reassessment of genetic limits to evolutionary change. Ecology. 2005, 86 (6): 1371-1384. 10.1890/04-1209.
Bonduriansky R: Convergent evolution of sexual shape dimorphism in diptera. Journal of Morphology. 2006, 267: 602-611. 10.1002/jmor.10426.
Wake DB, Larson A: Multidimensional analysis of an evolving lineage. Science. 1987, 238: 42-48. 10.1126/science.238.4823.42.
Kirschner M, Gerhart J: Evolvability. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (15): 8420-8427. 10.1073/pnas.95.15.8420.
Wagner GP, Altenberg L: Perspective: Complex adaptations and the evolution of evolvability. Evolution. 1996, 50 (3): 967-976. 10.2307/2410639.
West-Eberhard MJ: Evolution in the light of developmental and cell biology, and vice versa. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (15): 8417-8419. 10.1073/pnas.95.15.8417.
Dichtel-Danjoy ML, Felix MA: Phenotypic neighborhood and micro-evolvability. Trends in Genetics. 2004, 20 (5): 268-276. 10.1016/j.tig.2004.03.010.
Barton N, Partridge L: Limits to natural selection. Bioessays. 2000, 22 (12): 1075-1084. 10.1002/1521-1878(200012)22:12<1075::AID-BIES5>3.0.CO;2-M.
Houle D: Comparing Evolvability and Variability of Quantitative Traits. Genetics. 1992, 130 (1): 195-204.
Griswold CK: Pleiotropic mutation, modularity and evolvability. Evolution & Development. 2006, 8 (1): 81-93. 10.1111/j.1525-142X.2006.05077.x.
Wagner GP: The influence of variation and of developmental constraints on the rate of multivariate phenotypic evolution. Journal of Evolutionary Biology. 1988, 1: 45-66. 10.1046/j.1420-9101.1988.1010045.x.
Arthur W: Developmental drive: an important determinant of the direction of phenotypic evolution. Evolution & Development. 2001, 3 (4): 271-278. 10.1046/j.1525-142x.2001.003004271.x.
Prud'homme B, Gompel N, Carroll SB: Emerging principles of regulatory evolution. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (suppl. 1): 8605-8612. 10.1073/pnas.0700488104.
Hallgrímsson B, Brown JY, Ford-Hutchinson AF, Sheets HD, Zelditch ML, Jirik FR: The brachymorph mouse and the developmental-genetic basis for canalization and morphological integration. Evolution & Development. 2006, 8 (1): 61-73. 10.1111/j.1525-142X.2006.05075.x.
Milton CC, Batterham P, McKenzie JA, Hoffmann AA: Effect of E (sev) and Su (Raf) Hsp83 mutants and trans-heterozygotes on bristle trait means and variation in Drosophila melanogaster. Genetics. 2005, 171: 119-130. 10.1534/genetics.104.038463.
Salazar-Ciudad I, Jernvall J: How different types of pattern formation mechanisms affect the evolution of form and development. Evolution & Development. 2004, 6 (1): 6-16. 10.1111/j.1525-142X.2004.04002.x.
Richardson MK, Chipman AD: Developmental constraints in a comparative framework: a test case using variations in phalanx number during amniote evolution. J Exp Zoolog B Mol Dev Evol. 2003, 296 (1): 8-22. 10.1002/jez.b.13.
Arthur W: The effect of development on the direction of evolution: toward a twenty-first century consensus. Evolution & Development. 2004, 6 (4): 282-288. 10.1111/j.1525-142X.2004.04033.x.
Carroll SB: Chance and necessity: the evolution of morphological complexity and diversity. Nature. 2001, 409 (6823): 1102-1109. 10.1038/35059227.
Nijhout HF: Origin of butterfly wing patterns. The Character Concept in Evolutionary Biology. Edited by: Wagner GP. 2001, San Diego , Academic Press, 511-529.
Young N, Hallgrímsson B: Serial homology and the evolution of mammalian limb structure. Evolution. 2005, 59 (12): 2691-2704.
Jernvall J: Linking development with generation of novelty in mammalian teeth. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (6): 2641-2645. 10.1073/pnas.050586297.
Schluter D: Adaptive radiation along genetic lines of least resistance. Evolution. 1996, 50 (5): 1766-1774. 10.2307/2410734.
Arthur W: The interaction between developmental bias and natural selection: from centipede segments to a general hypothesis. Heredity. 2002, 89: 239-246. 10.1038/sj.hdy.6800139.
Emlen DJ: Artificial selection on horn length body size allometry in the horned beetle Onthophagus acuminatus (Coleoptera: Scarabaeidae). Evolution. 1996, 50 (3): 1219-1230. 10.2307/2410662.
Frankino WA, Zwaan BJ, Stern DL, Brakefield PM: Natural selection and developmental constraints in the evolution of allometries. Science. 2005, 307 (5710): 718-720. 10.1126/science.1105409.
Weber KE: Selection on wing allometry in Drosophila melanogaster. Genetics. 1990, 126: 975-989.
Rice SH: Developmental associations between traits: Covariance and beyond. Genetics. 2004, 166 (1): 513-526. 10.1534/genetics.166.1.513.
Roff DA: Evolutionary quantitative genetics. Evolutionary Genetics: Concepts and Case Studies. Edited by: Fox CW, Wolf JB. 2006, Oxford , Oxford University Press, 267-287.
Salazar-Ciudad I, Jernvall J: Graduality and innovation in the evolution of complex phenotypes: insights from development. J Exp Zoolog B Mol Dev Evol. 2005, 304 (6): 619-631. 10.1002/jez.b.21058.
Maynard Smith J, Sondhi KC: The genetics of a pattern. Genetics. 1960, 45: 1039-1050.
Davies RW, Workman PL: The genetic relationship of two quantitative characters in Drosophila melanogaster. I. Responses to selection and whole chromosome analysis. Genetics. 1971, 69: 353-361.
Guerra D, Pezzoli MC, Giorgi G, Garoia F, Cavicchi S: Developmental constraints in the Drosophila wing. Heredity. 1997, 79: 564-571.
Thompson JN: Studies on the nature and function of polygenic loci in Drosophila. I. Comparison of genomes from selection lines. Heredity. 1974, 33 (3): 373-387. 10.1038/hdy.1974.104.
Weber KE: How small are the smallest selectable domains of form?. Genetics. 1992, 130: 345-353.
Delph LF, Frey FM, Steven JC, Gehring JL: Investigating the independent evolution of the size of floral organs via G-matrix estimation and artificial selection. Evolution & Development. 2004, 6 (6): 438-448. 10.1111/j.1525-142X.2004.04052.x.
Rutledge JJ, Eisen EJ, Legates JE: Correlated response in skeletal traits and replicate variation in selected lines of mice. Theoretical and Applied Genetics. 1974, 46: 26-31. 10.1007/BF00281170.
Atchley WR, Rutledge JJ, Cowley DE: A multivariate statistical analysis of direct and correlated response to selection in the rat. Evolution. 1982, 36: 677-698. 10.2307/2407882.
Atchley WR, Cowley DE, Eisen EJ, Prasetyo H, Hawkins-Brown D: Correlated response in the developmenal choreographies of the mouse mandible to selection for body composition. Evolution. 1990, 44: 669-688. 10.2307/2409444.
Beldade P, Koops K, Brakefield PM: Modularity, individuality, and evo-devo in butterfly wings. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (22): 14262-14267. 10.1073/pnas.222236199.
Nijhout HF: Symmetry systems and compartments in lepidopteran wings - the evolution of a patterning mechanism. Development. 1994, 225-233.
Brakefield PM, French V: Eyespot development on butterfly wings: the epidermal response to damage. Developmental Biology. 1995, 168: 98-111. 10.1006/dbio.1995.1064.
Lyytinen A, Brakefield PM, Lindstrom L, Mappes J: Does predation maintain eyespot plasticity in Bicyclus anynana?. Proc Biol Sci. 2004, 271 (1536): 279-283. 10.1098/rspb.2003.2571.
Robertson KA, Monteiro A: Female Bicyclus anynana butterflies choose males on the basis of their dorsal UV-reflective eyespot pupils. Proc Biol Sci. 2005, 272 (1572): 1541-1546. 10.1098/rspb.2005.3142.
Nijhout HF: Pattern formation on lepidopteran wings: determination of an eyespot. Developmental Biology. 1980, 80 (2): 267-274. 10.1016/0012-1606(80)90403-0.
Beldade P, Brakefield PM: The genetics and evo-devo of butterfly wing patterns. Nature Reviews Genetics. 2002, 3 (6): 442-452.
Monteiro A, Brakefield PM, French V: The evolutionary genetics and developmental basis of wing pattern variation in the butterfly Bicyclus anynana. Evolution. 1994, 48 (4): 1147-1157. 10.2307/2410374.
Brunetti CR, Selegue JE, Monteiro A, French V, Brakefield PM, Carroll SB: The generation and diversification of butterfly eyespot color patterns. Current Biology. 2001, 11: 1578-1585. 10.1016/S0960-9822(01)00502-4.
Monteiro A, Brakefield PM, French V: Butterfly eyespots: the genetics and development of the color rings. Evolution. 1997, 51 (4): 1207-1216. 10.2307/2411050.
Nijhout F, Paulsen SM: Developmental models and polygenic characters. American Naturalist. 1997, 149 (2): 394-405. 10.1086/285996.
Beldade P, Brakefield PM: Concerted evolution and developmental integration in modular butterfly wing patterns. Evolution & Development. 2003, 5 (2): 169-179. 10.1046/j.1525-142X.2003.03025.x.
Holloway GJ, Brakefield PM, Kofman S: The genetics of wing pattern elements in the polyphenic butterfly, Bicyclus anynana. Heredity. 1993, 70: 179-186. 10.1038/hdy.1993.27.
Beldade P, Koops K, Brakefield PM: Developmental constraints versus flexibility in morphological evolution. Nature. 2002, 416 (6883): 844-847. 10.1038/416844a.
Klingenberg CP, Leamy LJ: Quantitative genetics of geometric shape in the mouse mandible. Evolution. 2001, 55 (11): 2342-2352.
Cheverud JM: Quantitative genetics and developmental constraints on evolution by selection. Journal of Theoretical Biology. 1984, 110: 155-172.
Cheverud JM: A comparison of genetic and phenotypic correlations. Evolution. 1988, 42: 958-968. 10.2307/2408911.
Lande R: Quantitative genetic analysis of multivariate evoluiton, applied to brain: body allometry. Evolution. 1979, 33: 402-416. 10.2307/2407630.
Arnold SJ: Constraints on phenotypic evolution. American Naturalist. 1992, 140: S85-S107. 10.1086/285398.
Nijhout HF, Rountree DB: Pattern induction across a homeotic boundary in the wings of Precis coenia (Lepidoptera: Nymphalidae). International Journal of Insect Morphology and Embryology. 1995, 24: 243-251. 10.1016/0020-7322(95)00004-N.
French V, Brakefield PM: Eyespot development on butterfly wings: the focal signal. Developmental Biology. 1995, 168: 112-123. 10.1006/dbio.1995.1065.
Nijhout HF: The Development and Evolution of Butterfly Wing Patterns. 1991, Washington, D.C. , Smithsonian Institute Press
Beldade P, Brakefield PM, Long AD: Contribution of Distal-less to quantitative variation in butterfly eyespots. Nature. 2002, 415 (6869): 315-318. 10.1038/415315a.
Monteiro A, Prijs J, Bax M, Hakkaart T, Brakefield PM: Mutants highlight the modular control of butterfly eyespot patterns. Evolution & Development. 2003, 5 (2): 180-187. 10.1046/j.1525-142X.2003.03029.x.
Reed RD, Serfas MS: Butterfly wing pattern evolution is associated with changes in a notch/distal-less temporal pattern formation process. Current Biology. 2004, 14: 1159-1166. 10.1016/j.cub.2004.06.046.
Weatherbee SD, Carroll SB: Selector genes and limb identity in arthropods and vertebrates. Cell. 1999, 97 (April 30): 283-286. 10.1016/S0092-8674(00)80737-0.
Klingenberg CP: Developmental constraints, modules, and evolvability. Variation. Edited by: Hallgrimsson B, Hall BK. 2005, San Diego , Academic Press, 219-247.
Brakefield PM: The power of evo-devo to explore evolutionary constraints: experiments with butterfly eyespots. Zoology. 2003, 106: 283-290. 10.1078/0944-2006-00124.
McMillan WO, Monteiro A, Kapan DD: Development and evolution on the wing. Trends in Ecology & Evolution. 2002, 17: 125-133. 10.1016/S0169-5347(01)02427-2.
Allen CE: The "eyespot module" and eyespots as modules: development, evolution, and integration of a complex phenotype. J Exp Zoolog B Mol Dev Evol. 2008, 310 (2): 179-190. 10.1002/jez.b.21186.
Paulsen SM: Quantitative genetics of butterfly wing color patterns. Developmental Genetics. 1994, 15 (1): 79-91. 10.1002/dvg.1020150109.
Paulsen SM, Nijhout HF: Phenotypic correlation structure among elements of the color pattern in Precis coenia (Lepidoptera: Nymphalidae). Evolution. 1993, 47: 593-618. 10.2307/2410073.
Breuker CJ, Gibbs M, Van Dyck H, Brakefield PM, Klingenberg CP, Van Dongen S: Integration of wings and their eyespots in the speckled wood butterfly Parage aegeria. J Exp Zoolog B Mol Dev Evol. 2007, 308 (4): 454-463. 10.1002/jez.b.21171.
Kingsolver JG, Wiernasz DC: Development, function, and the quantitative genetics of wing melanin pattern in Pieris butterflies. Evolution. 1991, 45 (6): 1480-1492. 10.2307/2409894.
Roff DA: A centennial celebration for quantitative genetics. Evolution. 2007, 61 (5): 1017-1032. 10.1111/j.1558-5646.2007.00100.x.
Rose MR, Passananti HB, Chippindale AK, Phelan JP, Matos M, Teotónio H, Mueller LD: The effects of evolution are local: evidence from experimental evolution in Drosophila. Integrative and Comparative Biology. 2005, 45: 486-491. 10.1093/icb/45.3.486.
Brakefield PM, El Filali E, Van der Laan R, Breuker CJ, Saccheri IJ, Zwaan B: Effective population size, reproductive success and sperm precedence in the butterfly, Bicyclus anynana, in captivity. Journal of Evolutionary Biology. 2001, 14 (1): 148-156. 10.1046/j.1420-9101.2001.00248.x.
Teotónio H, Rose MR: Variation in the reversibility of evolution. Nature. 2000, 408: 463-466. 10.1038/35044070.
Brodie ED: Genetic correlations between morphology and antipredator behaviour in natural populations of the garter snake Thamnophis ordinoides. Nature. 1989, 342 (6249): 542-543. 10.1038/342542a0.
Lande R: The genetic correlation between characters maintained by selection, linkage and inbreeding. Genetical Research. 1984, 44: 309-320.
Houle D: Genetic covariance of fitness correlates - what genetic correlations are made of and why it matters. Evolution. 1991, 45 (3): 630-648. 10.2307/2409916.
Gromko MH: Unpredictability of Correlated Response to Selection - Pleiotropy and Sampling Interact. Evolution. 1995, 49 (4): 685-693. 10.2307/2410322.
Wolf JB, Frankino WA, Agrawal AF, Brodie ED, Moore AJ: Developmental interactions and the constituents of quantitative variation. Evolution. 2001, 55: 232-245.
Raff RA: The Shape of Life. 1996, Chicago , University of Chicago Press
Phillips PC, McGuigan KL: Evolution of genetic variance-covariance structure. Evolutionary Genetics: Concepts and Case Studies. Edited by: Fox CW, Wolf JB. 2006, Oxford , Oxford University Press
Condamin M: Monographie du genre Bicyclus (Lepidoptera: Satyridae). Memoires de l'Institut Fondamental d'Afrique Noire. 1973, Dakar , IFAN
Brakefield PM, Roskam JC: Exploring evolutionary constraints. American Naturalist. 2006, 168: S4-S13. 10.1086/509049.
Conner JK: Artificial selection: a powerful tool for ecologists. Ecology. 2003, 84 (7): 1650-1660. 10.1890/0012-9658(2003)084[1650:ASAPTF]2.0.CO;2.
Lynch M, Walsh B: Genetics and analysis of quantitative traits. 1998, Sunderland, MA , Sinauer Associates, 980-
Becker WA: Manual of Quantitative Genetics. 1992, Pullman, WA , Academic Enterprises, 191-5
Monteiro A, Glaser G, Stockslager S, Glansdorp N, Ramos D: Comparative insights into questions of lepidopteran wing pattern homology. BMC Developmental Biology. 2006, 6: 52-10.1186/1471-213X-6-52.
Beldade P, French V, Brakefield PM: Developmental and genetic mechanisms for evolutionary diversification of serial repeats: eyespot size in Bicyclus anynana butterflies. J Exp Zoolog B Mol Dev Evol. 2008, 310 (2): 191-201. 10.1002/jez.b.21173.
We thank N. Wurzer and M. Lavrijsen and colleagues for rearing maize to support the experiment and K. Koops for assistance during the size selection experiment. J. Roskam provided photographs of Bicyclus species and M. Dix helped with artwork in Figure 2. We thank W. Frankino, T. van Dooren, J. Wolf, and three anonymous reviewers for helpful comments on previous versions of this manuscript. This work was supported by grant number 813.04.002 from the Netherlands Organization for Scientific Research (NWO).
CEA participated in the design of the study, carried out selection and wing damage experiments, performed the statistical analysis and drafted the manuscript. PB participated in the design of the study and carried out selection experiments. BJZ participated in the design and coordination of the study. PMB conceived of the study and participated in its design and coordination. All authors helped to draft the manuscript and all approved the final version.