Mutants highlight the modular control of butterfly eyespot patterns

Ultrabithorax modifies a regulatory network of genes essential for butterfly eyespot development in a wing sector-specific manner

Nymphalid butterfly species often have a different number of eyespots in forewings and hindwings, but how the hindwing identity gene Ultrabithorax (Ubx) drives this asymmetry is not fully understood. We examined a three-gene regulatory network for eyespot development in the hindwings of Bicyclus anynana butterflies and compared it with the same network previously described for forewings. We also examined how Ubx interacts with each of these three eyespot-essential genes. We found similar genetic interactions between the three genes in fore- and hindwings, but we discovered three regulatory differences: Antennapedia (Antp) merely enhances spalt (sal) expression in the eyespot foci in hindwings, but is not essential for sal activation, as in forewings; Ubx upregulates Antp in all hindwing eyespot foci but represses Antp outside these wing regions; and Ubx regulates sal in a wing sector-specific manner, i.e. it activates sal expression only in the sectors that have hindwing-specific eyespots. We propose a model for how the regulatory connections between these four genes evolved to produce wing- and sector-specific variation in eyespot number.

Phylogenetic analysis of the complete mitochondrial genome of the Blomfild's Beauty butterfly Smyrna blomfildia (Fabricius 1781) (Insecta: Lepidoptera: Nymphalidae: Nymphalini)

The Blomfild's Beauty butterfly Smyrna blomfildia (Fabricius 1781) (Lepidoptera: Nymphalidae: Nymphalini) is a sexually dimorphic species found in Mexico, Central, and South America. Males are territorial and are more vibrantly colored than females. Genome skimming by Illumina sequencing allowed the assembly of a complete circular mitochondrial genome (mitogenome) of 15,149 bp from S. blomfildia consisting of 83.9% AT nucleotides, 13 protein-coding genes, 22 tRNAs, two rRNAs, and a control region in the typical butterfly gene order. The S. blomfilda COX1 gene features an atypical start codon (CGA) while ATP6, COX1, COX2, CYTB, ND1, ND3, ND4, and ND5 display partial stop codons completed by the addition of 3' A residues to the mRNA. Bayesian phylogenetic reconstruction places Smyrna as a member of the tribe Nymphalini and sister to a clade containing genera Araschnia, Vanessa, Polygonia, and Aglais, which differs from its classic taxonomic placement in tribe Coeini.

The Fractal Geometry of the Nymphalid Groundplan: Self-Similar Configuration of Color Pattern Symmetry Systems in Butterfly Wings

Simple Summary

Highly diverse color patterns of butterfly wings can be explained as modifications of an archetypical color pattern of nymphalid butterflies called the nymphalid groundplan. The nymphalid groundplan contains three major symmetry systems and a discal symmetry system, but their relationships have been elusive. Here, the morphological and spatial relationships among these symmetry systems were studied based on cross-species color-pattern comparisons of the hindwings in nymphalid butterflies. It was shown that all symmetry systems can be expressed as various structures, suggesting the equivalence (homology) of these systems in developmental potential. In some cases, the discal symmetry system is circularly surrounded by the central symmetry system, which may then be surrounded by the border and basal symmetry systems, indicating a unified supersymmetry system covering the entire wing. These results suggest that butterfly color patterns are hierarchically constructed; one system is nested within another system, which is a self-similar relationship that achieves the fractal geometry. This self-similarity is likely mediated by the serial induction of organizers during development, and a possible mechanism is proposed for symmetry breaking of the system morphology, which contributes to the diversity of butterfly wing color patterns.

Abstract

The nymphalid groundplan is an archetypical color pattern of nymphalid butterflies involving three major symmetry systems and a discal symmetry system, which share the basic morphogenesis unit. Here, the morphological and spatial relationships among these symmetry systems were studied based on cross-species comparisons of nymphalid hindwings. Based on findings in Neope and Symbrenthia, all three major symmetry systems can be expressed as bands, spots, or eyespot-like structures, suggesting equivalence (homology) of these systems in developmental potential. The discal symmetry system can also be expressed as various structures. The discal symmetry system is circularly surrounded by the central symmetry system, which may then be surrounded by the border and basal symmetry systems, based mainly on findings in Agrias, indicating a unified supersymmetry system covering the entire wing. The border symmetry system can occupy the central part of the wing when the central symmetry system is compromised, as seen in Callicore. These results suggest that butterfly color patterns are hierarchically constructed in a self-similar fashion, as the fractal geometry of the nymphalid groundplan. This self-similarity is likely mediated by the serial induction of organizers, and symmetry breaking of the system morphology may be generated by the collision of opposing signals during development.

Multiple Loci Control Eyespot Number Variation on the Hindwings of Bicyclus anynana Butterflies

The underlying genetic changes that regulate the appearance and disappearance of repeated traits, or serial homologs, remain poorly understood. One hypothesis is that variation in genomic regions flanking master regulatory genes, also known as input-output genes, controls variation in trait number, making the locus of evolution almost predictable. Another hypothesis implicates genetic variation in up- or downstream loci of master control genes. Here, we use the butterfly Bicyclus anynana, a species that exhibits natural variation in eyespot number on the dorsal hindwing, to test these two hypotheses. We first estimated the heritability of dorsal hindwing eyespot number by breeding multiple butterfly families differing in eyespot number and regressing eyespot numbers of offspring on midparent values. We then estimated the number and identity of independent genetic loci contributing to eyespot number variation by performing a genome-wide association study with restriction site-associated DNA sequencing from multiple individuals varying in number of eyespots sampled across a freely breeding laboratory population. We found that dorsal hindwing eyespot number has a moderately high heritability of ∼0.50 and is characterized by a polygenic architecture. Previously identified genomic regions involved in eyespot development, and novel ones, display high association with dorsal hindwing eyespot number, suggesting that homolog number variation is likely determined by regulatory changes at multiple loci that build the trait, and not by variation at single master regulators or input-output genes.

Morphological Murals: The Scaling and Allometry of Butterfly Wing Patterns

The color patterns of butterflies moths are exceptionally diverse, but are very stable within a species, so that most species can be identified on the basis of their color pattern alone. The color pattern is established in the wing imaginal disc during a prolonged period of growth and differentiation, beginning during the last larval instar and ending during the first few days of the pupal stage. During this period, a variety of diffusion and reaction-diffusion signaling mechanisms determine the positions and sizes of the various elements that make up the overall color pattern. The patterning occurs while the wing is growing from a small imaginal disc to a very large pupal wing. One would therefore expect that some or all aspects of the color pattern would be sensitive to the size of the developmental field on which pattern formation takes place. To study this possibility, we analyzed the color patterns of Junonia coenia from animals whose growth patterns were altered by periodic starvation during larval growth, which produced individuals with a large range of variation in body size and wing size. Analyses of the color patterns showed that the positions and size of the pattern elements scaled perfectly isometrically with wing size. This is a puzzling finding and suggests the operation of a homeostatic or robustness mechanism that stabilizes pattern in spite of variation in the growth rate and final size of the wing.

Modularity, criticality, and evolvability of a developmental gene regulatory network

The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.

Induced mutations alter patterns of quantitative variation, phenotypic integration, and plasticity to elevated CO2 in Arabidopsis thaliana

A key step toward predicting responses to climate change is characterizing genetic variation in populations. While short-term responses will likely be shaped by currently available genetic variation, longer-term evolutionary responses will depend on the supply of novel variation by, ultimately, mutation. Studying mutational contributions to phenotypic variation can provide insights into the extent of potential variation on which selection may operate in future human-altered environments. Here we used the chemical mutagen ethyl methanesulfonate (EMS) to explore mutational contributions to phenotypic variation, integration, and plasticity to elevated carbon dioxide (eCO₂) in three accessions of Arabidopsis thaliana. We found that (1) mutagenesis increased broad-sense heritabilities and variation in plasticity to eCO₂ (genotype by environment interactions); (2) mutational effects varied among the three genetic backgrounds; (3) induced mutations had non-random (biased) effects on patterns of phenotypic integration. To our knowledge, this is the first study to address the effects of chemically induced mutations on phenotypic plasticity to eCO₂ in a model plant. We discuss our results in light of emerging insights from theoretical and empirical quantitative genetics, suggest potential avenues of research, and identify approaches that may help advance our understanding of climate-driven evolution in plants.

A new A-P compartment boundary and organizer in holometabolous insect wings

Decades of research on the highly modified wings of Drosophila melanogaster has suggested that insect wings are divided into two Anterior-Posterior (A-P) compartments separated by an axis of symmetry. This axis of symmetry is created by a developmental organizer that establishes symmetrical patterns of gene expression that in turn pattern the A-P axis of the wing. Butterflies possess more typical insect wings and butterfly wing colour patterns provide many landmarks for studies of wing structure and development. Using eyespot colour pattern variation in Vanessa butterflies, here we show an additional A-P axis of symmetry running between wing sectors 3 and 4. Boundaries of Drosophila mitotic clones suggest the existence of a previously undetected Far-Posterior (F-P) compartment boundary that coincides with this additional A-P axis. A similar compartment boundary is evident in butterfly mosaic gynandromorphs. We suggest that this additional compartment boundary and its associated developmental organizer create an axis of wing colour pattern symmetry and a gene expression-based combinatorial code, permitting each insect wing compartment to acquire a unique identity and allowing for the individuation of butterfly eyespots.

Developmental and evolutionary mechanisms shaping butterfly eyespots

Butterfly eyespots are visually compelling models to study the reciprocal interactions between evolutionary and developmental processes that shape phenotypic variation. They are evolutionarily diversified, ecologically relevant, and developmentally tractable, and have made key contributions to linking genotype, development, phenotype and fitness. Advances in the availability of analytical tools (e.g. gene editing and visualization techniques) and resources (e.g. genomic and transcriptomic data) are boosting the detailed dissection of the mechanisms underlying eyespot development and evolution. Here, we review current knowledge on the ecology, development, and evolution of butterfly eyespots, with focus on recent advances. We also highlight a number of unsolved mysteries in our understanding of the patterns and processes underlying the diversification of these structures.

Physiological Perturbation Reveals Modularity of Eyespot Development in the Painted Lady Butterfly, Vanessa cardui

Butterfly eyespots are complex morphological traits that can vary in size, shape and color composition even on the same wing surface. Homology among eyespots suggests they share a common developmental basis and function as an integrated unit in response to selection. Despite strong evidence of genetic integration, eyespots can also exhibit modularity or plasticity, indicating an underlying flexibility in pattern development. The extent to which particular eyespots or eyespot color elements exhibit modularity or integration is poorly understood, particularly following exposure to novel conditions. We used perturbation experiments to explore phenotypic correlations among different eyespots and their color elements on the ventral hindwing of V. cardui. Specifically, we identified which eyespots and eyespot features are most sensitive to perturbation by heat shock and injection of heparin-a cold shock mimic. For both treatments, the two central eyespots (3 + 4) were most affected by the experimental perturbations, whereas the outer eyespot border was more resistant to modification than the interior color elements. Overall, the individual color elements displayed a similar response to heat shock across all eyespots, but varied in their response to each other. Graphical modeling also revealed that although eyespots differ morphologically, regulation of eyespot size and colored elements appear to be largely integrated across the wing. Patterns of integration, however, were disrupted following heat shock, revealing that the strength of integration varies across the wing and is strongest between the two central eyespots. These findings support previous observations that document coupling between eyespots 3 + 4 in other nymphalid butterflies.

Nymphalid eyespots are co-opted to novel wing locations following a similar pattern in independent lineages

BACKGROUND

Variation in the number of repeated traits, or serial homologs, has contributed greatly to animal body plan diversity. Eyespot color patterns of nymphalid butterflies, like arthropod and vertebrate limbs, are an example of serial homologs. These eyespot color patterns originated in a small number of wing sectors on the ventral hindwing surface and later appeared in novel wing sectors, novel wings, and novel wing surfaces. However, the details of how eyespots were co-opted to these novel wing locations are currently unknown.

RESULTS

We used a large data matrix of eyespot/presence absence data, previously assembled from photographs of contemporary species, to perform a phylogenetic investigation of eyespot origins in nine independent nymphalid lineages. To determine how the eyespot gene regulatory network acquired novel positional information, we used phylogenetic correlation analyses to test for non-independence in the origination of eyespots. We found consistent patterns of eyespot gene network redeployment in the nine lineages, where eyespots first redeployed from the ventral hindwing to the ventral forewing, then to new sectors within the ventral wing surface, and finally to the dorsal wing surface. Eyespots that appeared in novel wing sectors modified the positional information of their serial homolog ancestors in one of two ways: by changing the wing or surface identity while retaining sector identity, or by changing the sector identity while retaining wing and surface identity.

CONCLUSIONS

Eyespot redeployment to novel sectors, wings, and surfaces happened multiple times in different nymphalid subfamilies following a similar pattern. This indicates that parallel mutations altering expression of the eyespot gene regulatory network led to its co-option to novel wing locations over time.

Nymphalid eyespot serial homologues originate as a few individualized modules

Serial homologues are repeated traits that share similar development but occur in different parts of the body. Variation in number of repeats accounts for substantial diversity in animal form and considerable work has focused on identifying the factors accounting for this variation. Little is known, however, about how serial homologues originally become repeated, or about the relative timing of repeat individuation relative to repeat origin. Here, we show that the serially repeated eyespots on nymphalid butterfly wings most likely arose as a small cluster of units on the ventral hindwing that were later co-opted to the dorsal and anterior wing surfaces. Based on comparative analyses of over 400 species, we found support for a model of eyespot origin followed by redeployment, rather than by the conventional model, where eyespots arose as a complete row of undifferentiated units that later gained individuation. In addition, eyespots most likely evolved from simpler pattern elements, single-coloured spots, which were already individuated among different wing sectors. Finally, the late appearance of eyespots on the dorsal, hidden wing surface further suggests that these novel complex traits originally evolved for one function (thwarting predator attacks) and acquired a second function (sexual signalling) when moved to a different body location. This broad comparative analysis illustrates how serial homologues may initially evolve as a few units serving a particular function and subsequently become repeated in novel body locations with new functions.

Origin, development, and evolution of butterfly eyespots

This article reviews the latest developments in our understanding of the origin, development, and evolution of nymphalid butterfly eyespots. Recent contributions to this field include insights into the evolutionary and developmental origin of eyespots and their ancestral deployment on the wing, the evolution of eyespot number and eyespot sexual dimorphism, and the identification of genes affecting eyespot development and black pigmentation. I also compare features of old and more recently proposed models of eyespot development and propose a schematic for the genetic regulatory architecture of eyespots. Using this schematic I propose two hypotheses for why we observe limits to morphological diversity across these serially homologous traits.

System-dependent regulations of colour-pattern development: a mutagenesis study of the pale grass blue butterfly

Developmental studies on wing colour patterns have been performed in nymphalid butterflies, but efficient genetic manipulations, including mutagenesis, have not been well established. Here, we have performed mutagenesis experiments in a lycaenid butterfly, the pale grass blue Zizeeria maha, to produce colour-pattern mutants. We fed the P-generation larvae an artificial diet containing the mutagen ethyl methane sulfonate (EMS), and the F₁- and F₂-generation adults showed various aberrant colour patterns: dorsoventral transformation, anterioposterior background colouration gap, weak contrast, disarrangement of spots, reduction of the size of spots, loss of spots, fusion of spots, and ectopic spots. Among them, the disarrangement, reduction, and loss of spots were likely produced by the coordinated changes of many spots of a single wing around the discal spot in a system-dependent manner, demonstrating the existence of the central symmetry system. The present study revealed multiple genetic regulations for system-dependent and wing-wide colour-pattern determination in lycaenid butterflies.

Baculovirus-mediated gene transfer in butterfly wings in vivo: an efficient expression system with an anti-gp64 antibody

BACKGROUND

Candidate genes for color pattern formation in butterfly wings have been known based on gene expression patterns since the 1990s, but their functions remain elusive due to a lack of a functional assay. Several methods of transferring and expressing a foreign gene in butterfly wings have been reported, but they have suffered from low success rates or low expression levels. Here, we developed a simple, practical method to efficiently deliver and express a foreign gene using baculovirus-mediated gene transfer in butterfly wings in vivo.

RESULTS

A recombinant baculovirus containing a gene for green fluorescent protein (GFP) was injected into pupae of the blue pansy butterfly Junonia orithya (Nymphalidae). GFP fluorescence was detected in the pupal wings and other body parts of the injected individuals three to five days post-injection at various degrees of fluorescence. We obtained a high GFP expression rate at relatively high virus titers, but it was associated with pupal death before color pattern formation in wings. To reduce the high mortality rate caused by the baculovirus treatment, we administered an anti-gp64 antibody, which was raised against baculovirus coat protein gp64, to infected pupae after the baculovirus injection. This treatment greatly reduced the mortality rate of the infected pupae. GFP fluorescence was observed in pupal and adult wings and other body parts of the antibody-treated individuals at various degrees of fluorescence. Importantly, we obtained completely developed wings with a normal color pattern, in which fluorescent signals originated directly from scales or the basal membrane after the removal of scales. GFP fluorescence in wing tissues spatially coincided with anti-GFP antibody staining, confirming that the fluorescent signals originated from the expressed GFP molecules.

CONCLUSIONS

Our baculovirus-mediated gene transfer system with an anti-gp64 antibody is reasonably efficient, and it can be an invaluable tool to transfer, express, and functionally examine foreign genes in butterfly wings and also in other non-model insect systems.

The pleiotropic structure of the genotype-phenotype map: the evolvability of complex organisms

It was first noticed 100 years ago that mutations tend to affect more than one phenotypic characteristic, a phenomenon that was called 'pleiotropy'. Because pleiotropy was found so frequently, the notion arose that pleiotropy is 'universal'. However, quantitative estimates of pleiotropy have not been available until recently. These estimates show that pleiotropy is highly restricted and are more in line with the notion of variational modularity than with universal pleiotropy. This finding has major implications for the evolvability of complex organisms and the mapping of disease-causing mutations.

Single locus affects embryonic segment polarity and multiple aspects of an adult evolutionary novelty

Background

The characterization of the molecular changes that underlie the origin and diversification of morphological novelties is a key challenge in evolutionary developmental biology. The evolution of such traits is thought to rely largely on co-option of a toolkit of conserved developmental genes that typically perform multiple functions. Mutations that affect both a universal developmental process and the formation of a novelty might shed light onto the genetics of traits not represented in model systems. Here we describe three pleiotropic mutations with large effects on a novel trait, butterfly eyespots, and on a conserved stage of embryogenesis, segment polarity.

Results

We show that three mutations affecting eyespot size and/or colour composition in Bicyclus anynana butterflies occurred in the same locus, and that two of them are embryonic recessive lethal. Using surgical manipulations and analysis of gene expression patterns in developing wings, we demonstrate that the effects on eyespot morphology are due to changes in the epidermal response component of eyespot induction. Our analysis of morphology and of gene expression in mutant embryos shows that they have a typical segment polarity phenotype, consistent with the mutant locus encoding a negative regulator of Wingless signalling.

Conclusions

This study characterizes the segregation and developmental effects of alleles at a single locus that controls the morphology of a lineage-specific trait (butterfly eyespots) and a conserved process (embryonic segment polarity and, specifically, the regulation of Wingless signalling). Because no gene with such function was found in the orthologous, highly syntenic genomic regions of two other lepidopterans, we hypothesize that our locus is a yet undescribed, possibly lineage-specific, negative regulator of the conserved Wnt/Wg pathway. Moreover, the fact that this locus interferes with multiple aspects of eyespot morphology and maps to a genomic region containing key wing pattern loci in different other butterfly species suggests it might correspond to a 'hotspot' locus in the diversification of this novel trait.

Phylogenetics and biogeography of a spectacular Old World radiation of butterflies: the subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrini)

Background

Butterflies of the subtribe Mycalesina (Nymphalidae: Satyrinae) are important model organisms in ecology and evolution. This group has radiated spectacularly in the Old World tropics and presents an exciting opportunity to better understand processes of invertebrate rapid radiations. However, the generic-level taxonomy of the subtribe has been in a constant state of flux, and relationships among genera are unknown. There are six currently recognized genera in the group. Mycalesis, Lohora and Nirvanopsis are found in the Oriental region, the first of which is the most speciose genus among mycalesines, and extends into the Australasian region. Hallelesis and Bicyclus are found in mainland Africa, while Heteropsis is primarily Madagascan, with a few species in Africa. We infer the phylogeny of the group with data from three genes (total of 3139 bp) and use these data to reconstruct events in the biogeographic history of the group.

Results

The results indicate that the group Mycalesina radiated rapidly around the Oligocene-Miocene boundary. Basal relationships are unresolved, but we recover six well-supported clades. Some species of Mycalesis are nested within a primarily Madagascan clade of Heteropsis, while Nirvanopsis is nested within Lohora. The phylogeny suggests that the group had its origin either in Asia or Africa, and diversified through dispersals between the two regions, during the late Oligocene and early Miocene. The current dataset tentatively suggests that the Madagascan fauna comprises two independent radiations. The Australasian radiation shares a common ancestor derived from Asia. We discuss factors that are likely to have played a key role in the diversification of the group.

Conclusions

We propose a significantly revised classification scheme for Mycalesina. We conclude that the group originated and radiated from an ancestor that was found either in Asia or Africa, with dispersals between the two regions and to Australasia. Our phylogeny paves the way for further comparative studies on this group that will help us understand the processes underlying diversification in rapid radiations of invertebrates.

Eyespot evolution: phylogenetic insights from Junonia and related butterfly genera (Nymphalidae: Junoniini)

Butterfly eyespots have been the focus of a number of developmental and evolutionary studies. However, a phylogenetic component has rarely been explicitly incorporated in these studies. In this study, I utilize a phylogeny to trace the evolution of eyespot number and position on the wing in a group of nymphalid butterflies, the subtribe Junoniini. These butterflies have two kinds of eyespot arrangements which I refer to as Serial and Individual. In the Serial arrangement, eyespots are placed in a series on compartments 1-6 (counting from the anterior wing margin). In the Individual arrangement, eyespots are isolated on specific compartments, ranging from 1 to 4 in number. This can be divided into four subtypes based on the number and positions of eyespots. I map the evolution of these five arrangements over a phylogeny of Junoniini reconstructed with ca. 3000 base pairs of sequence data from three genes. The results show that almost all arrangements have evolved at least twice, with multiple shifts between them by addition and deletion of eyespots. I propose a model involving genetic or developmental coupling between eyespots in specific compartments to explain these shifts. I discuss their evolution in light of existing knowledge about their development. I also discuss potential explanations for functional significance of the eyespot patterns found in the group. Differential selection for and against eyespots, both at different times over the phylogeny and in different regions, have driven the evolution of eyespot arrangements. The study throws open many questions about the adaptive significance of eyespots and the developmental underpinnings of the various arrangements.

Metamodels and phylogenetic replication: a systematic approach to the evolution of developmental pathways

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

Accommodating natural and sexual selection in butterfly wing pattern evolution

Visual patterns in animals may serve different functions, such as attracting mates and deceiving predators. If a signal is used for multiple functions, the opportunity arises for conflict among the different functions, preventing optimization for any one visual signal. Here we investigate the hypothesis that spatial separation of different visual signal functions has occurred in Bicyclus butterflies. Using phylogenetic reconstructions of character evolution and comparisons of evolutionary rates, we found dorsal surface characters to evolve at higher rates than ventral characters. Dorsal characters also displayed sex-based differences in evolutionary rates more often than did ventral characters. Thus, dorsal characters corresponded to our predictions of mate signalling while ventral characters appear to play an important role in predator avoidance. Forewing characters also fit a model of mate signalling, and displayed higher rates of evolution than hindwing characters. Our results, as well as the behavioural and developmental data from previous studies of Bicyclus species, support the hypothesis that spatial separation of visual signal functions has occurred in Bicyclus butterflies. This study is the first to demonstrate, in a phylogenetic framework, that spatial separation of signals used for mate signalling and those used for predator avoidance is a viable strategy to accommodate multiple signal functions. This signalling strategy has important ramifications on the developmental evolution of wing pattern elements and diversification of butterfly species.

Positional dependence of scale size and shape in butterfly wings: wing-wide phenotypic coordination of color-pattern elements and background

Butterfly wing color-patterns are a phenotypically coordinated array of scales whose color is determined as cellular interpretation outputs for morphogenic signals. Here we investigated distribution patterns of scale shape and size in relation to position and coloration on the hindwings of a nymphalid butterfly Junonia orithya. Most scales had a smooth edge but scales at and near the natural and ectopic eyespot foci and in the postbasal area were jagged. Scale size decreased regularly from the postbasal to distal areas, and eyespots occasionally had larger scales than the background. Reasonable correlations were obtained between the eyespot size and focal scale size in females. Histological and real-time individual observations of the color-pattern developmental sequence showed that the background brown and blue colors expanded from the postbasal to distal areas independently from the color-pattern elements such as eyespots. These data suggest that morphogenic signals for coloration directly or indirectly influence the scale shape and size and that the blue "background" is organized by a long-range signal from an unidentified organizing center in J. orithya.

A simulation study of mutations in the genetic regulatory hierarchy for butterfly eyespot focus determination

The color patterns on the wings of butterflies have been an important model system in evolutionary developmental biology. A recent computational model tested genetic regulatory hierarchies hypothesized to underlie the formation of butterfly eyespot foci [Evans, T.M., Marcus, J.M., 2006. A simulation study of the genetic regulatory hierarchy for butterfly eyespot focus determination. Evol. Dev. 8, 273-283]. The computational model demonstrated that one proposed hierarchy was incapable of reproducing the known patterns of gene expression associated with eyespot focus determination in wild-type butterflies, but that two slightly modified alternative hierarchies were capable of reproducing all of the known gene expressions patterns. Here we extend the computational models previously implemented in Delphi 2.0 to two mutants derived from the squinting bush brown butterfly (Bicyclus anynana). These two mutants, comet and Cyclops, have aberrantly shaped eyespot foci that are produced by different mechanisms. The comet mutation appears to produce a modified interaction between the wing margin and the eyespot focus that results in a series of comet-shaped eyespot foci. The Cyclops mutation causes the failure of wing vein formation between two adjacent wing-cells and the fusion of two adjacent eyespot foci to form a single large elongated focus in their place. The computational approach to modeling pattern formation in these mutants allows us to make predictions about patterns of gene expression, which are largely unstudied in butterfly mutants. It also suggests a critical experiment that will allow us to distinguish between two hypothesized genetic regulatory hierarchies that may underlie all butterfly eyespot foci.

Alternative models for the evolution of eyespots and of serial homology on lepidopteran wings

Serial homology is widespread in organismal design, but the origin and individuation of these repeated structures appears to differ with the different types of serial homologues, and remains an intriguing and exciting topic of research. Here I focus on the evolution of the serially repeated eyespots that decorate the margin of the wings of nymphalid butterflies. In this system, unresolved questions relate to the evolutionary steps that lead to the appearance of these serial homologues and how their separate identities evolved. I present and discuss two alternative hypotheses. The first proposes that eyespots first appeared as a row of undifferentiated repeated modules, one per wing compartment, that later become individuated. This individuation allowed eyespots to appear and disappear from specific wing compartments and also allowed eyespots to acquire different morphologies. The second hypothesis proposes that eyespots first appeared as individuated single units, or groups of units, that over evolutionary time were co-opted into new compartments on the wing. I discuss the merits of each of these alternate hypotheses by finding analogies to other systems and propose research avenues for addressing these issues in the future.

Developmental and genetic mechanisms for evolutionary diversification of serial repeats: eyespot size in Bicyclus anynana butterflies

Serially repeated pattern elements on butterfly wings offer the opportunity for integrating genetic, developmental, and functional aspects towards understanding morphological diversification and the evolution of individuality. We use captive populations of Bicyclus anynana butterflies, an emerging model in evolutionary developmental biology, to explore the genetic and developmental basis of compartmentalized changes in eyespot patterns. There is much variation for different aspects of eyespot morphology, and knowledge about the genetic pathways and developmental processes involved in eyespot formation. Also, despite the strong correlations across all eyespots in one butterfly, B. anynana shows great potential for independent changes in the size of individual eyespots. It is, however, unclear to what extent the genetic and developmental processes underlying eyespot formation change in a localized manner to enable such individualization. We use micromanipulations of developing wings to dissect the contribution of different components of eyespot development to quantitative differences in eyespot size on one wing surface. Reciprocal transplants of presumptive eyespot foci between artificial selection lines and controls suggest that while localized antagonistic changes in eyespot size rely mostly on localized changes in focal signal strength, concerted changes depend greatly on epidermal response sensitivities. This potentially reflects differences between the signal-response components of eyespot formation in the degrees of compartmentalization and/or the temporal pattern of selection. We also report on the phenotypic analysis of a number of mutant stocks demonstrating how single alleles can affect different eyespots in concert or independently, and thus contribute to the individualization of serially repeated traits.

Differences in the selection response of serially repeated color pattern characters: standing variation, development, and evolution

Background

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.

Results

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.

Conclusion

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.

The “Eyespot Module” and eyespots as modules: development, evolution, and integration of a complex phenotype

Organisms are inherently modular, yet modules also evolve in response to selection for functional integration or functional specialization of traits. For serially repeated homologous traits, there is a clear expectation that selection on the function of individual traits will reduce the integration between traits and subdivide a single ancestral module. The eyespots on butterfly wings are one example of serially repeated morphological traits that share a common developmental mechanism but are subject to natural and sexual selection for divergent functions. Here, I test two hypotheses about the organization of the eyespot pattern into independent dorsal-ventral and anterior-posterior modules, using a graphical modeling technique to examine patterns of eyespot covariation among and within wing surfaces in the butterfly Bicyclus anynana. Although there is a hierarchical and complex pattern of integration among eyespots, the results show a surprising mismatch between patterns of eyespot integration and the developmental and evolutionary eyespot units identified in previous empirical studies. These results are discussed in light of the relationships between developmental, functional, and evolutionary modules, and they suggest that developmental sources of independent trait variation are often masked by developmental sources of trait integration.

Integration of wings and their eyespots in the speckled wood butterfly Pararge aegeria

We investigated both the phenotypic and developmental integration of eyespots on the fore- and hindwings of speckled wood butterflies Pararge aegeria. Eyespots develop within a framework of wing veins, which may not only separate eyespots developmentally, but may at the same time also integrate them by virtue of being both signalling sources and barriers during eyespot development. We therefore specifically investigated the interaction between wing venation patterns and eyespot integration. Phenotypic covariation among eyespots was very high, but only eyespots in neighbouring wing cells and in homologous wing cells on different wing surfaces were developmentally integrated. This can be explained by the fact that the wing cells of these eyespots share one or more wing veins. The wing venation patterns of fore- and hindwings were highly integrated, both phenotypically and developmentally. This did not affect overall developmental integration of the eyespots. The adaptive significance of integration patterns is discussed and more specifically we stress the need to conduct studies on phenotypic plasticity of integration.

Reversed type of color-pattern modifications of butterfly wings: a physiological mechanism of wing-wide color-pattern determination

Application of cold shock or tungstate to butterfly pupae produces a unique color-pattern modification type on the adult wings, in which the color-pattern elements are dislocated toward the reduced focal elements. This modification-inducing activity has been primarily attributed to the putative cold-shock hormone (CSH) that is secreted into the hemolymph upon cold shock. Here, using a species of nymphalid butterfly Junonia almana, a new "reversed" type of the color-pattern modifications of butterfly wings was obtained by the application of heat shock or thapsigargin, a calcium-ATPase inhibitor, in which most elements were dislocated away from the enlarged focal elements. This result suggests that the endocrine secretion of CSH is sensitive to a wide range of temperature shocks, which then affects the cellular interpretation of the wing-wide positional information that is emitted from the focal locations. Ecdysteroid contributes to the wing-wide patterning primarily independently from CSH, but these two systems negatively interact with each other, probably in the intracellular signaling pathways.

The combined effect of two mutations that alter serially homologous color pattern elements on the fore and hindwings of a butterfly

BACKGROUND

The ability for serially homologous structures to acquire a separate identity has been primarily investigated for structures dependent on Hox gene input but is still incompletely understood in other systems. The fore and hindwings of butterflies are serially homologous structures as are the serially homologous eyespots that can decorate each of these wings. Eyespots can vary in number between fore and hindwings of the same individual and mutations of large effect can control the total number of eyespots that each of the wings displays. Here we investigate the genetics of a new spontaneous color pattern mutation, Missing, that alters eyespot number in the nymphalid butterfly, Bicyclus anynana. We further test the interaction of Missing with a previously described mutation, Spotty, describe the developmental stage affected by Missing, and test whether Missing is a mutant variant of the gene Distal-less via a linkage association study.

RESULTS

Missing removes or greatly reduces the size of two of the hindwing eyespots from the row of seven eyespots, with no detectable effect on the rest of the wing pattern. Offspring carrying a single Missing allele display intermediate sized eyespots at these positions. Spotty has the opposite effect of Missing, i.e., it introduces two extra eyespots in homologous wing positions to those affected by Missing, but on the forewing. When Missing is combined with Spotty the size of the two forewing eyespots decreases but the size of the hindwing spots stays the same, suggesting that these two mutations have a combined effect on the forewing such that Missing reduces eyespot size when in the presence of a Spotty mutant allele, but that Spotty has no effect on the hindwing. Missing prevents the complete differentiation of two of the eyespot foci on the hindwing. We found no evidence for any linkage between the Distal-less and Missing genes.

CONCLUSION

The spontaneous mutation Missing controls the differentiation of the signaling centers of a subset of the serial homologous eyespots present on both the fore and the hindwing in a dose-dependent fashion. The effect of Missing on the forewing, however, is only observed when the mutation Spotty introduces additional eyespots on this wing. Spotty, on the other hand, controls the differentiation of eyespot centers only on the forewing. Spotty, unlike Missing, may be under Ubx gene regulation, since it affects a subset of eyespots on only one of the serially homologous wings.

Exploring Evolutionary Constraints Is a Task for an Integrative Evolutionary Biology

Judging by the volume of writings about evolutionary constraints, they are an important topic in evolutionary biology. However, their involvement in shaping patterns of evolutionary change from morphological stasis to adaptive radiation remains contentious. This is at least in part because of the paucity of robust analyses of potential examples of constraints, whether of a more absolute or a relative nature. Here, we argue that what is needed to explore the type of constraints and bias on evolutionary change that may emerge from the way in which phenotypic variation is generated is an integrative approach applied to systems that can be tackled at different levels of biological organization. This is illustrated using research on the evolution of patterns in butterfly wing eyespots that has applied a combination of evolutionary genetics and evo-devo to an emerging model species with the beginnings of a comparative approach to describe patterns of variability among the extant taxa of two species-rich genera.

Comparative insights into questions of lepidopteran wing pattern homology

Background

Butterfly and moth eyespots can share a similar appearance, involving multiple concentric rings of colored scales, but usually occuring in non-homologous positions on the wing. Within the butterflies, on the other hand, spots that share the same homologous position may not share the concentric ring structure; and, in butterfly species that have eyespots with concentric rings, ectopic eyespots with a similar ring structure can be induced by means of a simple epidermal wound. The extent to which all these eyespots, natural or induced, share similar genes and developmental mechanisms is investigated here by means of protein in-situ localizations in selected butterfly and moth species. In addition to looking at some of the transcription factors previously identified as being involved in eyespot formation, we also tested the involvement of candidate genes from the Wingless and TGF-β signaling pathways as putative morphogens for eyespot development.

Results

Saturniid moth and nymphalid butterfly eyespots with concentric rings of color express at least two transcription factors, Distal-less and Engrailed, in the center of the future pattern. Nymphalid eyespots centers also express the ligand Wingless and an activated signal transducer, a phosphorylated Smad protein, but neither these proteins nor the previous two proteins are found in pierid spot centers, which consist of a single patch of color. Both butterfly wing patterns, however, express a third transcription factor, Spalt, a portion of whose expression domain maps to the black scales on the adult wing. Wounding a nymphalid wing, on the other hand, leads to upregulation of Distal-less, engrailed and spalt in subsets of cells around the wounding site, mimicking concentric eyespot development.

Conclusion

Wingless and TGF-β ligands are both candidate morphogens involved in nymphalid butterfly eyespot formation. These eyespots, as well as saturniid moth eyespots with concentric circles, share two genes that are associated with the differentiation of the signaling cells in nymphalid eyespots. This commonality suggests that they may be produced via the same developmental mechanism despite their non-homologous location. By contrast, pierid butterfly spots of a single color share some of the same genes but appear to be produced by a different mechanism. Eyespots with concentric rings may have co-opted a wound healing genetic network during their evolution.

Heliconius wing patterns: an evo-devo model for understanding phenotypic diversity

Evolutionary Developmental Biology aims for a mechanistic understanding of phenotypic diversity, and present knowledge is largely based on gene expression and interaction patterns from a small number of well-known model organisms. However, our understanding of biological diversification depends on our ability to pinpoint the causes of natural variation at a micro-evolutionary level, and therefore requires the isolation of genetic and developmental variation in a controlled genetic background. The colour patterns of Heliconius butterflies (Nymphalidae: Heliconiinae) provide a rich suite of naturally occurring variants with striking phenotypic diversity and multiple taxonomic levels of variation. Diversification in the genus is well known for its dramatic colour-pattern divergence between races or closely related species, and for Müllerian mimicry convergence between distantly related species, providing a unique system to study the development basis of colour-pattern evolution. A long history of genetic studies has showed that pattern variation is based on allelic combinations at a surprisingly small number of loci, and recent developmental evidence suggests that pattern development in Heliconius is different from the eyespot determination of other butterflies. Fine-scale genetic mapping studies have shown that a shared toolkit of genes is used to produce both convergent and divergent phenotypes. These exciting results and the development of new genomic resources make Heliconius a very promising evo-devo model for the study of adaptive change.

Functional evo-devo

Functional factors such as optimal design and adaptive value have been the central concern of evolutionary biology since the advent of the New Synthesis. By contrast, evolutionary developmental biology (evo-devo) has concentrated primarily on structural factors such as the ways in which body parts can be built. These different emphases have stood in the way of an integrated understanding of the role of development in evolution. Here, we try to bridge this gap by outlining the relevance of functional factors in evo-devo. We use modularity and the view of development as a flexible evolutionary system to outline a unified perspective that includes both structural and functional aspects.

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

A simulation study of the genetic regulatory hierarchy for butterfly eyespot focus determination

The color patterns on the wings of butterflies have been an important model system in evolutionary developmental biology. Two types of models have been used to study these patterns. The first type of model employs computational techniques and generalized mechanisms of pattern formation to make predictions about how color patterns will vary as parameters of the model are changed. These generalized mechanisms include diffusion gradient, reaction-diffusion, lateral inhibition, and threshold responses. The second type of model uses known genetic interactions from Drosophila melanogaster and patterns of candidate gene expression in one of several butterfly species (most often Junonia (Precis) coenia or Bicyclus anynana) to propose specific genetic regulatory hierarchies that appear to be involved in color pattern formation. This study combines these two approaches using computational techniques to test proposed genetic regulatory hierarchies for the determination of butterfly eyespot foci (also known as border ocelli foci). Two computer programs, STELLA 8.1 and Delphi 2.0, were used to simulate the determination of eyespot foci. Both programs revealed weaknesses in a genetic model previously proposed for eyespot focus determination. On the basis of these simulations, we propose two revised models for eyespot focus determination and identify components of the genetic regulatory hierarchy that are particularly sensitive to changes in model parameter values. These components may play a key role in the evolution of butterfly eyespots. Simulations like these may be useful tools for the study of other evolutionary developmental model systems and reveal similar sensitive components of the relevant genetic regulatory hierarchies.

Pharmacologic approaches to butterfly wing patterning: sulfated polysaccharides mimic or antagonize cold shock and alter the interpretation of gradients of positional information

Butterflies produce complex and diverse wing patterns by mechanisms that are generally unknown. We have employed a pharmacological approach to explore the molecular mechanisms of pattern formation. In a screen of over 200 compounds injected into developing Junonia coenia pupae, we identified several specific sulfated polysaccharides that caused widespread, dose-dependent effects on adult wing patterns. These compounds were well tolerated and permitted butterflies to eclose normally and take flight at moderate levels of effect. Heparin and closely related chondroitin sulfates caused stage-specific expansion of distal and proximal band systems and reduction and repatterning of eyespots. Dextran sulfate and fucoidan, whose structures are widely divergent from heparin and one another, caused contraction of distal and proximal systems, but had no effect on eyespots. Nonsulfated or nonpolymeric saccharides were without effect. Pattern alterations were indistinguishable from those reported for extreme cold shock and exposure to sodium tungstate and "molsin". When administered after cold shock or coinjected with heparin, dextran sulfate reversed all patterning effects. We suggest that the primary effect of polysaccharide treatments is to alter the interpretation of gradients of positional information along the proximodistal axis of the pupal wing.

Morphological comparison of pupal wing cuticle patterns in butterflies

Butterfly wing color-patterns are determined in the prospective wing tissues during the late larval and early pupal stages. To study the cellular differentiation process of wings, morphological knowledge on pupal wings is prerequisite. Here we systematically examined morphological patterns of the pupal wing cuticular surface in a wide variety of nymphalid butterflies in relation to adult color-patterns. Several kinds of pupal wing patterns corresponding to particular adult color-pattern elements were widely observed in many species. Especially noteworthy were the pupal "focal" spots corresponding to the adult border ocelli system, which were detected in many species of Nymphalinae, Apaturinae, Argynninae, Satyrinae, and Danainae. Striped patterns on the pupal wing cuticle seen in some species of Limenitinae, Ariadnae, and Marpesiinae directly corresponded to those of the adult wings. In Vanessa cardui, eyespot-like pattern elements were tentatively produced during development in the wing tissue underneath the pupal spots and subsequently erased, suggesting a mechanism for producing novel color-patterns in the course of development and evolution. The pupal focal spots reasonably correlated with the adult eyespots in size in Precis orithya and Ypthima argus. We physically damaged the pupal focal spots and their corresponding cells underneath in these species, which abolished or inhibited the formation of the adult eyespots. Taken together, our results clarified that pupal cuticle patterns were often indicative of the adult color-patterns and apparently reflect molecular activity of organizing centers for the adult color-pattern formation at least in nymphalid butterflies.

Jumping genes and AFLP maps: transforming lepidopteran color pattern genetics

The color patterns on the wings of lepidopterans are among the most striking patterns in nature and have inspired diverse biological hypotheses such as the ecological role of aposomatic coloration, the evolution of mimicry, the role of human activities in industrial melanism, and the developmental basis of phenotypic plasticity. Yet, the developmental mechanisms underlying color pattern development are not well understood for three reasons. First, few mutations that alter color patterns have been characterized at the molecular level, so there is little mechanistic understanding of how mutant phenotypes are produced. Second, although gene expression patterns resembling adult color patterns are suggestive, there are few data available showing that gene products have a functional role in color pattern formation. Finally, because with few exceptions (notably Bombyx), genetic maps for most species of Lepidoptera are rudimentary or nonexistent, it is very difficult to characterize spontaneous mutants or to determine whether mutations with similar phenotypes are because of lesions in the same gene or different genes. Discussed here are two strategies for overcoming these difficulties: germ-line transformation of lepidopteran species using transposon vectors and amplified frequency length polymorphism-based genetic mapping using variation between divergent strains within a species or between closely related and interfertile species. These advances, taken together, will create new opportunities for the characterization of existing genetic variants, the creation of new sequence-tagged mutants, and the testing of proposed functional genetic relationships between gene products, and will greatly facilitate our understanding of the evolution and development of lepidopteran color patterns.

Germline transformation of the butterfly Bicyclus anynana

Ecological and evolutionary theory has frequently been inspired by the diversity of colour patterns on the wings of butterflies. More recently, these varied patterns have also become model systems for studying the evolution of developmental mechanisms. A technique that will facilitate our understanding of butterfly colour-pattern development is germline transformation. Germline transformation permits functional tests of candidate gene products and of cis-regulatory regions, and provides a means of generating new colour-pattern mutants by insertional mutagenesis. We report the successful transformation of the African satyrid butterfly Bicyclus anynana with two different transposable element vectors, Hermes and piggyBac, each carrying EGFP coding sequences driven by the 3XP3 synthetic enhancer that drives gene expression in the eyes. Candidate lines identified by screening for EGFP in adult eyes were later confirmed by PCR amplification of a fragment of the EGFP coding sequence from genomic DNA. Flanking DNA surrounding the insertions was amplified by inverse PCR and sequenced. Transformation rates were 5% for piggyBac and 10.2% for Hermes. Ultimately, the new data generated by these techniques may permit an integrated understanding of the developmental genetics of colour-pattern formation and of the ecological and evolutionary processes in which these patterns play a role.

Color-pattern Modifications and Speciation in Butterflies of the Genus Vanessa and its Related Genera Cynthia and Bassaris

We have previously shown that the systemic injection of sodium tungstate, a protein-tyrosine phosphatase (PTPase) inhibitor, to pupae immediately after pupation efficiently produces characteristic color-pattern modifications on the wings of many species of butterflies including Vanessa indica and Cynthia cardui. In these species, the black spots reduced in size in response to the treatment. Similar modifications are occasionally seen in the field-caught aberrant individuals. Exceptionally, however, a C. cardui individual with enlarged black spots ("reversed" modification pattern) has been reported. Here we show that these modified patterns of V. indica and C. cardui are quite similar to the normal color-patterns of other Vanessa species. V. indica with tungstate-induced modifications resembled V. tameamea, V. samani, and Bassaris itea, whereas V. dilecta, V. atalanta, and V. dejeanii are similar to the "reversed" individual. Most features seen in the experimentally-modified V. indica were observed throughout the fore- and hindwings of V. samani. In contrast, the experimentally-induced color-patterns of C. cardui did not parallel variation of Cynthia butterflies. Since it has been proposed that a hypothetical transduction pathway with a PTPase for the scale-cell differentiation globally coordinates the wing-wide color-patterns, our findings suggest that spontaneous mutations in genes in this hypothetical pathway might have played a major role in creating new color-patterns and species in the Vanessa genus but not in the Cynthia genus. This evolutionary mechanism may probably be shared more widely in Lepidoptera, although this would not be a sole determinant for the color-pattern development and evolution.