Butterfly Eyespot Patterns: Evidence for Specification by a Morphogen Diffusion Gradient

In this paper we describe a test for Nijhout's (1978, 1980a) hypothesis that the eyespot patterns on butterfly wings are the result of a threshold reaction of the epidermal cells to a concentration gradient of a diffusing degradable morphogen produced by focal cells at the centre of the future eyespot. The wings of the nymphalid butterfly, Bicyclus anynana, have a series of eyespots, each composed of a white pupil, a black disc and a gold outer ring. In earlier extirpation and transplantation experiments (Nijhout 1980a; French and Brakefield, 1995) it has been established that these eyespots are indeed organised around groups of signalling cells active during the first hours of pupal development. If these cells were to supply the positional information for eyespot formation in accordance with Nijhout's diffusion-degradation gradient model, then, when two foci are close together, the signals should sum, and this effect should be apparent in the detailed shape of the resulting pigment pattern. We give an equation for the form of the contours that would be obtained in this manner. We use this to test the morphogen gradient hypothesis on measurements of the outlines of fused eyespots obtained either by grafting focal cells close together, or by using a mutation (Spotty) that produces adjacent fused eyespots. The contours of the fused patterns were found to satisfy our equation, thus corroborating Nijhout's hypothesis to the extent possible with this particular type of experiment.     

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.     

A Model for Selection of Eyespots on Butterfly Wings

Unsolved Problem

The development of eyespots on the wing surface of butterflies of the family Nympalidae is one of the most studied examples of biological pattern formation.However, little is known about the mechanism that determines the number and precise locations of eyespots on the wing. Eyespots develop around signaling centers, called foci, that are located equidistant from wing veins along the midline of a wing cell (an area bounded by veins). A fundamental question that remains unsolved is, why a certain wing cell develops an eyespot, while other wing cells do not.

Key Idea and Model

We illustrate that the key to understanding focus point selection may be in the venation system of the wing disc. Our main hypothesis is that changes in morphogen concentration along the proximal boundary veins of wing cells govern focus point selection. Based on previous studies, we focus on a spatially two-dimensional reaction-diffusion system model posed in the interior of each wing cell that describes the formation of focus points. Using finite element based numerical simulations, we demonstrate that variation in the proximal boundary condition is sufficient to robustly select whether an eyespot focus point forms in otherwise identical wing cells. We also illustrate that this behavior is robust to small perturbations in the parameters and geometry and moderate levels of noise. Hence, we suggest that an anterior-posterior pattern of morphogen concentration along the proximal vein may be the main determinant of the distribution of focus points on the wing surface. In order to complete our model, we propose a two stage reaction-diffusion system model, in which an one-dimensional surface reaction-diffusion system, posed on the proximal vein, generates the morphogen concentrations that act as non-homogeneous Dirichlet (i.e., fixed) boundary conditions for the two-dimensional reaction-diffusion model posed in the wing cells. The two-stage model appears capable of generating focus point distributions observed in nature.

Result

We therefore conclude that changes in the proximal boundary conditions are sufficient to explain the empirically observed distribution of eyespot focus points on the entire wing surface. The model predicts, subject to experimental verification, that the source strength of the activator at the proximal boundary should be lower in wing cells in which focus points form than in those that lack focus points. The model suggests that the number and locations of eyespot foci on the wing disc could be largely controlled by two kinds of gradients along two different directions, that is, the first one is the gradient in spatially varying parameters such as the reaction rate along the anterior-posterior direction on the proximal boundary of the wing cells, and the second one is the gradient in source values of the activator along the veins in the proximal-distal direction of the wing cell.     

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.     

Generation of butterfly wing eyespot patterns: a model for morphological determination of eyespot and parafocal element

The determination of color patterns of butterfly wing eyespots has been explained by the morphogen concentration gradient model. The induction model has been proposed recently as a more realistic alternative, in which the eyespot-specifying signal does not depend entirely on focal activity. However, this model requires further elaboration and supporting evidence to be validated. Here, I examined various color patterns of nymphalid butterflies to propose the mechanics of the induction model. Based on cases in which an eyespot light ring is identical to the background in color, I propose that eyespots are fundamentally composed of dark rings and non-dark "background" spaces between them. In the induction model, the dark-ring-inducing signal that is released from a prospective eyespot focus (the primary organizing center) as a slow-moving wave effects both selfenhancement and peripheral induction of the dark-ring-inhibitory signal at the secondary organizing centers, resulting in an eyespot that has alternate dark and light rings. Moreover, there are cases in which an unseen "imaginary light ring" surrounds an eyespot proper and in which PFEs are integrated into the eyespot. It appears that PFEs constitute a periodic continuum of eyespot dark rings; thus, a background space between the eyespot and a PFE is mechanistically equivalent to eyespot light rings. The eyespot dark-ring-inducing signals and PFE-inducing signal are likely to be identical in quality, but released at different times from the same organizing center. Computer simulations based on the reaction-diffusion system support the feasibility of the induction model.     

Colour pattern homology and evolution in Vanessa butterflies (Nymphalidae: Nymphalini): eyespot characters

Ocelli are serially repeated colour patterns on the wings of many butterflies. Eyespots are elaborate ocelli that function in predator avoidance and deterrence as well as in mate choice. A phylogenetic approach was used to study ocelli and eyespot evolution in Vanessa butterflies, a genus exhibiting diverse phenotypes among these serial homologs. Forty‐four morphological characters based on eyespot number, arrangement, shape and the number of elements in each eyespot were defined and scored. Ocelli from eight wing cells on the dorsal and ventral surfaces of the forewing and hindwing were evaluated. The evolution of these characters was traced over a phylogeny of Vanessa based on 7750 DNA base pairs from 10 genes. Our reconstruction predicts that the ancestral Vanessa had 5 serially arranged ocelli on all four wing surfaces. The ancestral state on the dorsal forewing and ventral hindwing was ocelli arranged in two heterogeneous groups. On the dorsal hindwing, the ancestral state was either homogenous or ocelli arranged in two heterogeneous groups. On the ventral forewing, we determined that the ancestral state was organized into three heterogeneous groups. In Vanessa, almost all ocelli are individuated and capable of independent evolution relative to other colour patterns except for the ocelli in cells ?1 and 0 on the dorsal and ventral forewings, which appear to be constrained to evolve in parallel. The genus Vanessa is a good model system for the study of serial homology and the interaction of selective forces with developmental architecture to produce diversity in butterfly colour patterns.     

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.     

Deflective and intimidating eyespots: a comparative study of eyespot size and position in Junonia butterflies

Eyespots are conspicuous circular features found on the wings of several lepidopteran insects. Two prominent hypotheses have been put forth explaining their function in an antipredatory role. The deflection hypothesis posits that eyespots enhance survival in direct physical encounters with predators by deflecting attacks away from vital parts of the body, whereas the intimidation hypothesis posits that eyespots are advantageous by scaring away a potential predator before an attack. In the light of these two hypotheses, we investigated the evolution of eyespot size and its interaction with position and number within a phylogenetic context in a group of butterflies belonging to the genus Junonia. We found that larger eyespots tend to be found individually, rather than in serial dispositions. Larger size and conspicuousness make intimidating eyespots more effective, and thus, we suggest that our results support an intimidation function in some species of Junonia with solitary eyespots. Our results also show that smaller eyespots in Junonia are located closer to the wing margin, thus supporting predictions of the deflection hypothesis. The interplay between size, position, and arrangement of eyespots in relation to antipredation and possibly sexual selection, promises to be an interesting field of research in the future. Similarly, further comparative work on the evolution of absolute eyespot size in natural populations of other butterfly groups is needed.     

Eyespot placement and assembly in the green alga Chlamydomonas

The eyespot organelle of the green alga Chlamydomonas allows the cell to phototax toward (or away) from light to maximize the light intensity for photosynthesis and minimize photo-damage. At cytokinesis, the eyespot is resorbed at the cleavage furrow and two new eyespots form in the daughter cells 180 degrees from each other. The eyespots are positioned asymmetrically with respect to the microtubule cytoskeleton. Eyespots are assembled from all three chloroplast membranes and carotenoid-filled granules, which form a sandwich structure overlaid by the tightly apposed plasma membrane. This review describes (1) my interest in cellular asymmetry and organelle biology, (2) isolation of mutations that describe four genes governing eyespot placement and assembly, (3) the characterization of the EYE2 gene, which encodes a thioredoxin superfamily member, and (4) the characterization of the MIN1 gene, which is required for the layered organization of granules and membranes in the eyespot. BioEssays 25:410-416, 2003.     

Thioredoxin-family protein EYE2 and Ser/Thr kinase EYE3 play interdependent roles in eyespot assembly

The eyespot of the biflagellate unicellular green alga Chlamydomonas reinhardtii is a complex organelle that facilitates directional responses of the cell to environmental light stimuli. The eyespot, which assembles de novo after every cell division and is associated with the daughter four-membered (D4) microtubule rootlet, comprises an elliptical patch of rhodopsin photoreceptors on the plasma membrane and stacks of carotenoid-rich pigment granule arrays in the chloroplast. Two loci, EYE2 and EYE3, define factors involved in the formation and organization of the eyespot pigment granule arrays. Whereas EYE3, a serine/threonine kinase of the ABC1 family, localizes to pigment granules, EYE2 localization corresponds to an area of the chloroplast envelope in the eyespot. EYE2 is positioned along, and adjacent to, the D4 rootlet in the absence of pigment granules. The eyespot pigment granule array is required for maintenance of the elliptical shape of both the overlying EYE2 and channelrhodopsin-1 photoreceptor patches. We propose a model of eyespot assembly wherein rootlet and photoreceptor direct EYE2 to an area of the chloroplast envelope, where it acts to facilitate assembly of pigment granule arrays, and EYE3 plays a role in the biogenesis of the pigment granules.     

The role of eyespots as anti-predator mechanisms, principally demonstrated in the Lepidoptera

Eyespots are found in a variety of animals, in particular lepidopterans. The role of eyespots as antipredator mechanisms has been discussed since the 19th Century, with two main hypotheses invoked to explain their occurrence. The first is that large, centrally located eyespots intimidate predators by resembling the eyes of the predators' own enemies; the second, though not necessarily conflicting, hypothesis is that small, peripherally located eyespots function as markers to deflect the attacks of predators to non-vital regions of the body. A third possibility is also proposed; that eyespots intimidate predators merely because they are novel or rarely encountered salient features. These hypotheses are reviewed, with special reference given to avian predators, since these are likely to be the principal visually hunting predators of the lepidopterans considered. Also highlighted is the necessity to consider the potential influence of sexual selection on lepidopteran wing patterns, and the genetics and development of eyespot formation.     

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.     

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.     

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.     

Body size determines eyespot size and presence in coral reef fishes

Numerous organisms display conspicuous eyespots. These eye‐like patterns have been shown to effectively reduce predation by either deflecting strikes away from nonvital organs or by intimidating potential predators. While investigated extensively in terrestrial systems, determining what factors shape eyespot form in colorful coral reef fishes remains less well known. Using a broadscale approach we ask: How does the size of the eyespot relate to the actual eye, and at what size during ontogeny are eyespots acquired or lost? We utilized publicly available images to generate a dataset of 167 eyespot‐bearing reef fish species. We measured multiple features relating to the size of the fish, its eye, and the size of its eyespot. In reef fishes, the area of the eyespot closely matches that of the real eye; however, the eyespots “pupil” is nearly four times larger than the real pupil. Eyespots appear at about 20 mm standard length. However, there is a marked decrease in the presence of eyespots in fishes above 48 mm standard length; a size which is tightly correlated with significant decreases in documented mortality rates. Above 75–85 mm, the cost of eyespots appears to outweigh their benefit. Our results identify a “size window” for eyespots in coral reef fishes, which suggests that eyespot use is strictly body size‐dependent within this group.     

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.     

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.     

Conserved developmental processes and the formation of evolutionary novelties: examples from butterfly wings

The origin and diversification of evolutionary novelties-lineage-specific traits of new adaptive value-is one of the key issues in evolutionary developmental biology. However, comparative analysis of the genetic and developmental bases of such traits can be difficult when they have no obvious homologue in model organisms. The finding that the evolution of morphological novelties often involves the recruitment of pre-existing genes and/or gene networks offers the potential to overcome this challenge. Knowledge about shared developmental processes obtained from extensive studies in model organisms can then be used to understand the origin and diversification of lineage-specific structures. Here, we illustrate this approach in relation to eyespots on the wings of Bicyclus anynana butterflies. A number of spontaneous mutations isolated in the laboratory affect eyespots, lepidopteran-specific features, and also processes that are shared by most insects. We discuss how eyespot mutants with disturbed embryonic development may help elucidate the genetic pathways involved in eyespot formation, and how venation mutants with altered eyespot patterns might shed light on mechanisms of eyespot development.     

Models for cell differentiation and generation of polarity in diffusion-governed morphogenetic fields

Models based on molecular mechanisms are presented for pattern formation in developing organisms. It is assumed that there exists a diffusion governed gradient in the morphogenetic field. It is shown that cellular differentiation and the subsequent pattern formation result from the interaction of the diffusing morphogen with the genetic regulatory mechanism of cells. In a second stage it is shown that starting from a homogeneous distribution of morphogen, polarity can be generated spontaneously in the morphogenetic field giving rise to the establishment of a gradient. The stability of these gradients is demonstrated. The onset of a morphogenetic gradient and pattern formation are combined in a single coherent model. Size invariance and its biological implications are discussed.