Hox genes are essential for the development of eyespots in Bicyclus anynana butterflies

The eyespot patterns found on the wings of nymphalid butterflies are novel traits that originated first in hindwings and subsequently in forewings, suggesting that eyespot development might be dependent on Hox genes. Hindwings differ from forewings in the expression of Ultrabithorax (Ubx), but the function of this Hox gene in eyespot development as well as that of another Hox gene Antennapedia (Antp), expressed specifically in eyespots centers on both wings, are still unclear. We used CRISPR-Cas9 to target both genes in Bicyclus anynana butterflies. We show that Antp is essential for eyespot development on the forewings and for the differentiation of white centers and larger eyespots on hindwings, whereas Ubx is essential not only for the development of at least some hindwing eyespots but also for repressing the size of other eyespots. Additionally, Antp is essential for the development of silver scales in male wings. In summary, Antp and Ubx, in addition to their conserved roles in modifying serially homologous segments along the anterior–posterior axis of insects, have acquired a novel role in promoting the development of a new set of serial homologs, the eyespot patterns, in both forewings (Antp) and hindwings (Antp and Ubx) of B. anynana butterflies. We propose that the peculiar pattern of eyespot origins on hindwings first, followed by forewings, could be due to an initial co-option of Ubx into eyespot development followed by a later, partially redundant, co-option of Antp into the same network.     

Eyespots deflect predator attack increasing fitness and promoting the evolution of phenotypic plasticity

Some eyespots are thought to deflect attack away from the vulnerable body, yet there is limited empirical evidence for this function and its adaptive advantage. Here, we demonstrate the conspicuous ventral hindwing eyespots found on Bicyclus anynana butterflies protect against invertebrate predators, specifically praying mantids. Wet season (WS) butterflies with larger, brighter eyespots were easier for mantids to detect, but more difficult to capture compared to dry season (DS) butterflies with small, dull eyespots. Mantids attacked the wing eyespots of WS butterflies more frequently resulting in greater butterfly survival and reproductive success. With a reciprocal eyespot transplant, we demonstrated the fitness benefits of eyespots were independent of butterfly behaviour. Regardless of whether the butterfly was WS or DS, large marginal eyespots pasted on the hindwings increased butterfly survival and successful oviposition during predation encounters. In previous studies, DS B. anynana experienced delayed detection by vertebrate predators, but both forms suffered low survival once detected. Our results suggest predator abundance, identity and phenology may all be important selective forces for B. anynana. Thus, reciprocal selection between invertebrate and vertebrate predators across seasons may contribute to the evolution of the B. anynana polyphenism.     

THE EVOLUTIONARY GENETICS AND DEVELOPMENTAL BASIS OF WING PATTERN VARIATION IN THE BUTTERFLY BICYCLUS ANYNANA

We have studied interactions between developmental processes and genetic variation for the eyespot color pattern on the adult dorsal forewing of the nymphalid butterfly, Bicyclus anynana. Truncation selection was applied in both an upward and a downward direction to the size of a single eyespot consisting of rings with wing scales of differing color pigments. High heritabilities resulted in rapid responses to selection yielding divergent lines with very large or very small eyespots. Strong correlated responses occurred in most of the other eyespots on each wing surface. The cells at the center of a presumptive eyespot (the “focus”) act in the early pupal stage to establish the adult wing pattern. The developmental fate of the scale cells within an eyespot is specified by the “signaling” properties of the focus and the “response” thresholds of the epidermis. The individual eyespots can be envisaged as developmental homologues. Grafting experiments performed with the eyespot foci of the selected lines showed that additive genetic variance exists for both the response and, in particular, the signaling components of the developmental system. The results are discussed in the context of how constraints on the evolution of this wing pattern may be related to the developmental organization.     

Mutants highlight the modular control of butterfly eyespot patterns

SUMMARY The eyespots on butterfly wings are thought to be serially homologous pattern elements. Yet eyespots differ greatly in number, shape, color, and size, within and among species. To what extent do these serially homologues have separate developmental identities, upon which selection acts to create diversity? We examined x‐ray–induced mutations for the eyespots of the nymphalid butterfly Bicyclus anynana that highlight the modular control of these serially homologous wing pattern elements. These mutations reduce or eliminate individual eyespots, or groups of eyespots, with no further effect on the wing color pattern. The collection of mutants highlights a greater potential developmental repertoire than that observed across the genus Bicyclus. We studied in detail one such mutation, of codominant effect, that causes the elimination of two adjacent eyespots on the ventral hindwing. By analyzing the expression of genes known to be involved in eyespot formation, we found an alteration in the differentiation of the “organizing” cells at the eyespot's center. No such cells differentiate in the wing subdivisions lacking the two eyespots in the mutants. We propose several developmental models, based on wing compartmentalization in Drosophila, that provide the first framework for thinking about the molecular evolution of butterfly wing pattern modularity.     

Butterfly Wings Are Three-Dimensional: Pupal Cuticle Focal Spots and Their Associated Structures in Junonia Butterflies

Butterfly wing color patterns often contain eyespots, which are developmentally determined at the late larval and early pupal stages by organizing activities of focal cells that can later form eyespot foci. In the pupal stage, the focal position of a future eyespot is often marked by a focal spot, one of the pupal cuticle spots, on the pupal surface. Here, we examined the possible relationships of the pupal focal spots with the underneath pupal wing tissues and with the adult wing eyespots using Junonia butterflies. Large pupal focal spots were found in two species with large adult eyespots, J. orithya and J. almana, whereas only small pupal focal spots were found in a species with small adult eyespots, J. hedonia. The size of five pupal focal spots on a single wing was correlated with the size of the corresponding adult eyespots in J. orithya. A pupal focal spot was a three-dimensional bulge of cuticle surface, and the underside of the major pupal focal spot exhibited a hollowed cuticle in a pupal case. Cross sections of a pupal wing revealed that the cuticle layer shows a curvature at a focal spot, and a positional correlation was observed between the cuticle layer thickness and its corresponding cell layer thickness. Adult major eyespots of J. orithya and J. almana exhibited surface elevations and depressions that approximately correspond to the coloration within an eyespot. Our results suggest that a pupal focal spot is produced by the organizing activity of focal cells underneath the focal spot. Probably because the focal cell layer immediately underneath a focal spot is thicker than that of its surrounding areas, eyespots of adult butterfly wings are three-dimensionally constructed. The color-height relationship in adult eyespots might have an implication in the developmental signaling for determining the eyespot color patterns.     

Over-expression of Ultrabithorax alters embryonic body plan and wing patterns in the butterfly Bicyclus anynana

In insects, forewings and hindwings usually have different shapes, sizes, and color patterns. A variety of RNAi experiments across insect species have shown that the hox gene Ultrabithorax (Ubx) is necessary to promote hindwing identity. However, it remains unclear whether Ubx is sufficient to confer hindwing fate to forewings across insects. Here, we address this question by over-expressing Ubx in the butterfly Bicyclus anynana using a heat-shock promoter. Ubx whole-body over-expression during embryonic and larvae development led to body plan changes in larvae but to mere quantitative changes to adult morphology, respectively. Embryonic heat-shocks led to fused segments, loss of thoracic and abdominal limbs, and transformation of head limbs to larger appendages. Larval heat-shocks led to reduced eyespot size in the expected homeotic direction, but neither additional eyespots nor wing shape changes were observed in forewings as expected of a homeotic transformation. Interestingly, Ubx was found to be expressed in a novel, non-characteristic domain – in the hindwing eyespot centers. Furthermore, ectopic expression of Ubx on the pupal wing activated the eyespot-associated genes spalt and Distal-less, known to be directly repressed by Ubx in the fly?s haltere and leg primordia, respectively, and led to the differentiation of black wing scales. These results suggest that Ubx has been co-opted into a novel eyespot gene regulatory network, and that it is capable of activating black pigmentation in butterflies.     

BUTTERFLY EYESPOTS: THE GENETICS AND DEVELOPMENT OF THE COLOR RINGS

The butterfly Bicyclus anynana has a series of distal eyespots on its wings. Each eyespot is composed of a white pupil, a black disc, and a gold outer ring. We applied artificial selection to the large dorsal eyespot on the forewing to produce a line with the gold ring reduced or absent (BLACK) and another line with a reduced black disc and a broad gold ring (GOLD). High heritabilities, coupled with a rapid response to selection, produced two lines of butterflies with very different phenotypes. Other eyespots showed a correlated change in the proportion of their color rings. Surgical experiments were performed on pupal wings from the different lines at the time of eyespot pattern specification. They showed that the additive genetic variance for this trait was in the response of the wing epidermis to signaling from the organizing cells at the eyespot center (the focus). This response was found to vary across different regions of the wing and also between the sexes. The particular eyespot color composition found for each sex, as well as the maintenance of the high genetic variation, are discussed with reference to the ecology of the butterfly, sexual selection, and visual selection by predators.     

Temporal Gene Expression Variation Associated with Eyespot Size Plasticity in Bicyclus anynana

Seasonal polyphenism demonstrates an organism''s ability to respond to predictable environmental variation with alternative phenotypes, each presumably better suited to its respective environment. However, the molecular mechanisms linking environmental variation to alternative phenotypes via shifts in development remain relatively unknown. Here we investigate temporal gene expression variation in the seasonally polyphenic butterfly Bicyclus anynana. This species shows drastic changes in eyespot size depending on the temperature experienced during larval development. The wet season form (larvae reared over 24°C) has large ventral wing eyespots while the dry season form (larvae reared under 19°C) has much smaller eyespots. We compared the expression of three proteins, Notch, Engrailed, and Distal-less, in the future eyespot centers of the two forms to determine if eyespot size variation is associated with heterochronic shifts in the onset of their expression. For two of these proteins, Notch and Engrailed, expression in eyespot centers occurred earlier in dry season than in wet season larvae, while Distal-less showed no temporal difference between the two forms. These results suggest that differences between dry and wet season adult wings could be due to a delay in the onset of expression of these eyespot-associated genes. Early in eyespot development, Notch and Engrailed may be functioning as repressors rather than activators of the eyespot gene network. Alternatively, temporal variation in the onset of early expressed genes between forms may have no functional consequences to eyespot size regulation and may indicate the presence of an ''hourglass'' model of development in butterfly eyespots.     

Quantitative genetics of butterfly wing color patterns

Developmental processes exert their influence on the evolution of complex morphologies through the genetic correlations they engender between traits. Butterfly wing color patterns provide a model system to examine this connection between development and evolution. In butterflies, the nymphalid groundplan is a framework used to decompose complex wing patterns into their component pattern elements. The first goal of this work has been to determine whether the components of the nymphalid groundplan are the products of independent developmental processes. To test this hypothesis, the genetic correlation matrices for two species of butterflies, Precis coenia and Precis evarete, were estimated for 27 wing pattern characters. The second purpose was to test the hypothesis that the differentiation of serial homologs lowers their genetic correlations. The “eyespots” found serially repeated across the fore- and hindwing and on the dorsal and ventral wing surfaces provided an opportunity to test this hypothesis. The genetic correlation matrices of both species were very similar. The pattern of genetic correlation measured between the different types of pattern elements and between the homologous repeats of a pattern element supported the first hypothesis of developmental independence among the elements of the groundplan. The correlation pattern among the differentiated serial homologs was similarly found to support the second hypothesis: pairs of eyespots that had differentiated had lower genetic correlations than pairs that were similar in morphology. The implications of this study are twofold: First, the apparent developmental independence among the distinct elements of wing pattern has facilitated the vast diversification in morphology found in butterflies. Second, the lower genetic correlations betweendifferentiated homologs demonstrates that developmental constraints can in fact be broken. The extent to which genetic correlations readily change, however, remains unknown. © 1994 Wiley-Liss, Inc.     

Spatial patterns of correlated scale size and scale color in relation to color pattern elements in butterfly wings

Complex butterfly wing color patterns are coordinated throughout a wing by unknown mechanisms that provide undifferentiated immature scale cells with positional information for scale color. Because there is a reasonable level of correspondence between the color pattern element and scale size at least in Junonia orithya and Junonia oenone, a single morphogenic signal may contain positional information for both color and size. However, this color–size relationship has not been demonstrated in other species of the family Nymphalidae. Here, we investigated the distribution patterns of scale size in relation to color pattern elements on the hindwings of the peacock pansy butterfly Junonia almana, together with other nymphalid butterflies, Vanessa indica and Danaus chrysippus. In these species, we observed a general decrease in scale size from the basal to the distal areas, although the size gradient was small in D. chrysippus. Scales of dark color in color pattern elements, including eyespot black rings, parafocal elements, and submarginal bands, were larger than those of their surroundings. Within an eyespot, the largest scales were found at the focal white area, although there were exceptional cases. Similarly, ectopic eyespots that were induced by physical damage on the J. almana background area had larger scales than in the surrounding area. These results are consistent with the previous finding that scale color and size coordinate to form color pattern elements. We propose a ploidy hypothesis to explain the color–size relationship in which the putative morphogenic signal induces the polyploidization (genome amplification) of immature scale cells and that the degrees of ploidy (gene dosage) determine scale color and scale size simultaneously in butterfly wings.     

The relationship between eyespot shape and wing shape in the butterfly Bicyclus anynana: A genetic and morphometrical approach

The African butterfly, Bicyclus anynana, normally possesses circular eyespots on its wings. Artificial selection lines, which express ellipsoidal eyespots on the dorsal surface of the forewing, were used to investigate correlated changes in wing shape. Morphometric analysis of linear wing measurements and wing scale counts provided evidence that eyespot shape was correlated with localised shape changes in the corresponding wing-cell, with overall shape changes in the wing, and with the density/arrangement of scales around the eyespot area.     

Color-pattern analysis of eyespots in butterfly wings: a critical examination of morphogen gradient models

Butterfly wing color patterns consist of many color-pattern elements such as eyespots. It is believed that eyespot patterns are determined by a concentration gradient of a single morphogen species released by diffusion from the prospective eyespot focus in conjunction with multiple thresholds in signal-receiving cells. As alternatives to this single-morphogen model, more flexible multiple-morphogen model and induction model can be proposed. However, the relevance of these conceptual models to actual eyespots has not been examined systematically. Here, representative eyespots from nymphalid butterflies were analyzed morphologically to determine if they are consistent with these models. Measurement of ring widths of serial eyespots from a single wing surface showed that the proportion of each ring in an eyespot is quite different among homologous rings of serial eyespots of different sizes. In asymmetric eyespots, each ring is distorted to varying degrees. In extreme cases, only a portion of rings is expressed remotely from the focus. Similarly, there are many eyespots where only certain rings are deleted, added, or expanded. In an unusual case, the central area of an eyespot is composed of multiple "miniature eyespots," but the overall macroscopic eyespot structure is maintained. These results indicate that each eyespot ring has independence and flexibility to a certain degree, which is less consistent with the single-morphogen model. Considering a "periodic eyespot", which has repeats of a set of rings, damage-induced eyespots in mutants, and a scale-size distribution pattern in an eyespot, the induction model is the least incompatible with the actual eyespot diversity.     

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.     

Physiologically induced color-pattern changes in butterfly wings: mechanistic and evolutionary implications

A mechanistic understanding of the butterfly wing color-pattern determination can be facilitated by experimental pattern changes. Here I review physiologically induced color-pattern changes in nymphalid butterflies and their mechanistic and evolutionary implications. A type of color-pattern change can be elicited by elemental changes in size and position throughout the wing, as suggested by the nymphalid groundplan. These changes of pattern elements are bi-directional and bi-sided dislocation toward or away from eyespot foci and in both proximal and distal sides of the foci. The peripheral elements are dislocated even in the eyespot-less compartments. Anterior spots are more severely modified, suggesting the existence of an anterior-posterior gradient. In one species, eyespots are transformed into white spots with remnant-like orange scales, and such patterns emerge even at the eyespot-less "imaginary" foci. A series of these color-pattern modifications probably reveal "snap-shots" of a dynamic morphogenic signal due to heterochronic uncoupling between the signaling and reception steps. The conventional gradient model can be revised to account for these observed color-pattern changes.     

A single origin for nymphalid butterfly eyespots followed by widespread loss of associated gene expression

Understanding how novel complex traits originate involves investigating the time of origin of the trait, as well as the origin of its underlying gene regulatory network in a broad comparative phylogenetic framework. The eyespot of nymphalid butterflies has served as an example of a novel complex trait, as multiple genes are expressed during eyespot development. Yet the origins of eyespots remain unknown. Using a dataset of more than 400 images of butterflies with a known phylogeny and gene expression data for five eyespot-associated genes from over twenty species, we tested origin hypotheses for both eyespots and eyespot-associated genes. We show that eyespots evolved once within the family Nymphalidae, approximately 90 million years ago, concurrent with expression of at least three genes associated with early eyespot development. We also show multiple losses of expression of most genes from this early three-gene cluster, without corresponding losses of eyespots. We propose that complex traits, such as eyespots, may have originated via co-option of a large pre-existing complex gene regulatory network that was subsequently streamlined of genes not required to fulfill its novel developmental function.     

Automatic recognition and measurement of butterfly eyespot patterns

A favorite wing pattern element in butterflies that has been the focus of intense study in evolutionary and developmental biology, as well as in behavioral ecology, is the eyespot. Because the pace of research on these bull's eye patterns is accelerating we sought to develop a tool to automatically detect and measure butterfly eyespot patterns in digital images of the wings. We used a machine learning algorithm with features based on circularity and symmetry to detect eyespots on the images. The algorithm is first trained with examples from a database of images with two different labels (eyespot and non-eyespot), and subsequently is able to provide classification for a new image. After an eyespot is detected the radius measurements of its color rings are performed by a 1D Hough Transform which corresponds to histogramming. We trained software to recognize eyespot patterns of the nymphalid butterfly Bicyclus anynana but eyespots of other butterfly species were also successfully detected by the software.     

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

SUMMARY 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.     

Wnt signaling underlies evolution and development of the butterfly wing pattern symmetry systems

Most butterfly wing patterns are proposed to be derived from a set of conserved pattern elements known as symmetry systems. Symmetry systems are so-named because they are often associated with parallel color stripes mirrored around linear organizing centers that run between the anterior and posterior wing margins. Even though the symmetry systems are the most prominent and diverse wing pattern elements, their study has been confounded by a lack of knowledge regarding the molecular basis of their development, as well as the difficulty of drawing pattern homologies across species with highly derived wing patterns. Here we present the first molecular characterization of symmetry system development by showing that WntA expression is consistently associated with the major basal, discal, central, and external symmetry system patterns of nymphalid butterflies. Pharmacological manipulations of signaling gradients using heparin and dextran sulfate showed that pattern organizing centers correspond precisely with WntA, wingless, Wnt6, and Wnt10 expression patterns, thus suggesting a role for Wnt signaling in color pattern induction. Importantly, this model is supported by recent genetic and population genomic work identifying WntA as the causative locus underlying wing pattern variation within several butterfly species. By comparing the expression of WntA between nymphalid butterflies representing a range of prototypical symmetry systems, slightly deviated symmetry systems, and highly derived wing patterns, we were able to infer symmetry system homologies in several challenging cases. Our work illustrates how highly divergent morphologies can be derived from modifications to a common ground plan across both micro- and macro-evolutionary time scales.