Mathematical modeling of the eyespots in butterfly wings

Butterfly wing color patterns are a representative model system for studying biological pattern formation, due to their two-dimensional simple structural and high inter- and intra-specific variabilities. Moreover, butterfly color patterns have demonstrated roles in mate choice, thermoregulation, and predator avoidance via disruptive coloration, attack deflection, aposematism, mimicry, and masquerade. Because of the importance of color patterns to many aspects of butterfly biology and their apparent tractability for study, color patterns have been the subjects of many attempts to model their development. Early attempts focused on generalized mechanisms of pattern formation such as reaction-diffusion, diffusion gradient, lateral inhibition, and threshold responses, without reference to any specific gene products. As candidate genes with expression patterns that resembled incipient color patterns were identified, genetic regulatory networks were proposed for color pattern formation based on gene functions inferred from other insects with wings, such as Drosophila. Particularly detailed networks incorporating the gene products, Distal-less (Dll), Engrailed (En), Hedgehog (Hh), Cubitus interruptus (Ci), Transforming growth factor-β (TGF-β), and Wingless (Wg), have been proposed for butterfly border ocelli (eyespots) which helps the investigation of the formation of these patterns. Thus, in this work, we develop a mathematical model including the gene products En, Hh, Ci, TGF-β, and Wg to mimic and investigate the eyespot formation in butterflies. Our simulations show that the level of En has peaks in the inner and outer rings and the level of Ci has peaks in the inner and middle rings. The interactions among these peaks activate cells to produce white, black, and yellow pigments in the inner, middle, and outer rings, respectively, which captures the eyespot pattern of wild type Bicyclus anynana butterflies. Additionally, our simulations suggest that lack of En generates a single black spot and lack of Hh or Ci generates a single white spot, and a deficiency of TGF-β or Wg will cause the loss of the outer yellow ring. These deficient patterns are similar to those observed in the eyespots of Vanessa atalanta, Vanessa altissima, and Chlosyne nycteis. Thus, our model also provides a hypothesis to explain the mechanism of generating the deficient patterns in these species.

Playing the blues, the greens and the reds with cellulose-based structural colours

Structural vivid colours can arise from the interference of light reflected from structures exhibiting periodicity on scales in the range of visible wavelengths. This effect is observed with light reflected from cell-walls of some plants and exoskeletons of certain insects. Sometimes the colour sequence observed for these structures consists of nearly circular concentric rings that vary in colour from Red, Orange, Yellow, Green, Cyan to Blue, from the periphery to the centre, similarly to the colour scheme sequence observed for the rainbow (ROYGB). The sequence of colours has been found for solid films obtained from droplets of aqueous cellulose nanocrystals (CNCs) suspensions and attributed to a "coffee ring" effect. In this work, coloured lyotropic solutions and solid films obtained from a cellulose derivative in the presence of trifluoroacetic acid (TFA), which acts as a "reactive solvent", are revisited. The systems were investigated with spectroscopy, using circularly and linearly polarised light, coupled with a polarised optical microscope (POM) and scanning electron microscopy (SEM). The lyotropic cholesteric liquid crystalline solutions were confined in capillaries to simplify 1D molecular diffusion along the capillary where an unexpected sequence of the structural colours was observed. The development and reappearance of the sequence of vivid colours seem consistent with the reaction-diffusion of the "reactive solvent" in the presence of the cellulosic chains. The strong TFA acts as an auto-catalyst for the chemical reaction between TFA and the hydroxyl groups, existing along the cellulosic chain, and diffuses to the top and bottom along the capillaries, carrying dissolved cellulosic chains. Uncovering the precise mechanism of colour sequence and evolution over time in cellulosic lyotropic solutions has important implications for future optical/sensors applications and for the understanding of the development of cellulose-based structures in nature.

Symmetry systems on the wings of Dichromodes Guenée (Lepidoptera: Geometridae) are unconstrained by venation

The nymphalid groundplan, an idealized schematic illustrating the essential elements of butterfly wing patterns, predicts a consistent relationship between color pattern and wing venation. Moths in the family Geometridae have wing shapes and patterns that often resemble those of butterflies, and until recently, this family was believed to be among butterflies’ closest relatives. However, an examination of the geometrid genus Dichromodes Guenée, 1858 shows no consistent relationship between the central symmetry system and wing venation. Whereas the distal edge of the central symmetry system is predicted to reach the costal margin proximal to the Subcostal vein in butterflies and acronictine moths, it has no consistent relationship with the Subcostal, Radius, or Radial Sector 1 veins in Dichromodes. This finding highlights developmental diversity that was previously overlooked due to the overwhelming preference for butterflies in studies of lepidopteran wing patterns.

The genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

Activation of butterfly eyespots by Distal-less is consistent with a reaction-diffusion process

Eyespots on the wings of nymphalid butterflies represent colorful examples of pattern formation, yet the developmental origins and mechanisms underlying eyespot center differentiation are still poorly understood. Using CRISPR-Cas9 we re-examine the function of Distal-less (Dll) as an activator or repressor of eyespots, a topic that remains controversial. We show that the phenotypic outcome of CRISPR mutations depends upon which specific exon is targeted. In Bicyclus anynana, exon 2 mutations are associated with both missing and ectopic eyespots, and also exon skipping. Exon 3 mutations, which do not lead to exon skipping, produce only null phenotypes, including missing eyespots, lighter wing coloration and loss of scales. Reaction-diffusion modeling of Dll function, using Wnt and Dpp as candidate morphogens, accurately replicates these complex crispant phenotypes. These results provide new insight into the function of Dll as a potential activator of eyespot development, scale growth and melanization, and suggest that the tuning of Dll expression levels can generate a diversity of eyespot phenotypes, including their appearance on the wing.This article has an associated 'The people behind the papers' interview.

A Literal Elytral Rainbow: Tunable Structural Colors Using Single Diamond Biophotonic Crystals in Pachyrrhynchus congestus Weevils

The brilliant colors of many insects arise from the interference of incident light with complex nanostructured biomaterials that are present in their cuticle. Here, the rainbow-colored spots on the elytra of a snout weevil, Pachyrrhynchus congestus pavonius (Coleoptera: Curculionidae), are investigated using synchrotron small-angle X-ray scattering, scanning electron microscopy, microspectrophotometry, and photonic bandgap modeling. It is shown that the iridescent scales present in the rainbow-hued spots are due to a 3D photonic crystal network of chitin in air with a single diamond (Fd-3m) symmetry. In many insects, different orientations of photonic crystal domains are used to create various hues. In this weevil, however, both the chitin volume fractions as well as the lattice parameters of the biologically self-assembled single diamond nanostructure are varied to achieve the remarkable tuning of the structural colors across the visible light spectrum on a scale-by-scale basis. Uncovering the precise mechanism of color tuning employed by this weevil has important implications for further structural and developmental research on biophotonic nanostructures and may provide fresh impetus for bioinspired and biomimetic multifunctional applications, as synthesis of photonic crystals at visible length scales is currently challenging.

Waiting in the wings: what can we learn about gene co-option from the diversification of butterfly wing patterns?

A major challenge is to understand how conserved gene regulatory networks control the wonderful diversity of form that we see among animals and plants. Butterfly wing patterns are an excellent example of this diversity. Butterfly wings form as imaginal discs in the caterpillar and are constructed by a gene regulatory network, much of which is conserved across the holometabolous insects. Recent work in Heliconius butterflies takes advantage of genomic approaches and offers insights into how the diversification of wing patterns is overlaid onto this conserved network. WntA is a patterning morphogen that alters spatial information in the wing. Optix is a transcription factor that acts later in development to paint specific wing regions red. Both of these loci fit the paradigm of conserved protein-coding loci with diverse regulatory elements and developmental roles that have taken on novel derived functions in patterning wings. These discoveries offer insights into the 'Nymphalid Ground Plan', which offers a unifying hypothesis for pattern formation across nymphalid butterflies. These loci also represent 'hotspots' for morphological change that have been targeted repeatedly during evolution. Both convergent and divergent evolution of a great diversity of patterns is controlled by complex alleles at just a few genes. We suggest that evolutionary change has become focused on one or a few genetic loci for two reasons. First, pre-existing complex cis-regulatory loci that already interact with potentially relevant transcription factors are more likely to acquire novel functions in wing patterning. Second, the shape of wing regulatory networks may constrain evolutionary change to one or a few loci. Overall, genomic approaches that have identified wing patterning loci in these butterflies offer broad insight into how gene regulatory networks evolve to produce diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Focusing on butterfly eyespot focus: uncoupling of white spots from eyespot bodies in nymphalid butterflies

BACKGROUND

Developmental studies on butterfly wing color patterns often focus on eyespots. A typical eyespot (such as that of Bicyclus anynana) has a few concentric rings of dark and light colors and a white spot (called a focus) at the center. The prospective eyespot center during the early pupal stage is known to act as an organizing center. It has often been assumed, according to gradient models for positional information, that a white spot in adult wings corresponds to an organizing center and that the size of the white spot indicates how active that organizing center was. However, there is no supporting evidence for these assumptions. To evaluate the feasibility of these assumptions in nymphalid butterflies, we studied the unique color patterns of Calisto tasajera (Nymphalidae, Satyrinae), which have not been analyzed before in the literature.

RESULTS

In the anterior forewing, one white spot was located at the center of an eyespot, but another white spot associated with either no or only a small eyespot was present in the adjacent compartment. The anterior hindwing contained two adjacent white spots not associated with eyespots, one of which showed a sparse pattern. The posterior hindwing contained two adjacent pear-shaped eyespots, and the white spots were located at the proximal side or even outside the eyespot bodies. The successive white spots within a single compartment along the midline in the posterior hindwing showed a possible trajectory of a positional determination process for the white spots. Several cases of focus-less eyespots in other nymphalid butterflies were also presented.

CONCLUSIONS

These results argue for the uncoupling of white spots from eyespot bodies, suggesting that an eyespot organizing center does not necessarily differentiate into a white spot and that a prospective white spot does not necessarily signify organizing activity for an eyespot. Incorporation of these results in future models for butterfly wing color pattern formation is encouraged.

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.

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.

Evolutionary history of the recruitment of conserved developmental genes in association to the formation and diversification of a novel trait

Background

The origin and modification of novel traits are important aspects of biological diversification. Studies combining concepts and approaches of developmental genetics and evolutionary biology have uncovered many examples of the recruitment, or co-option, of genes conserved across lineages for the formation of novel, lineage-restricted traits. However, little is known about the evolutionary history of the recruitment of those genes, and of the relationship between them -for example, whether the co-option involves whole or parts of existing networks, or whether it occurs by redeployment of individual genes with de novo rewiring. We use a model novel trait, color pattern elements on butterfly wings called eyespots, to explore these questions. Eyespots have greatly diversified under natural and sexual selection, and their formation involves genetic circuitries shared across insects.

Results

We investigated the evolutionary history of the recruitment and co-recruitment of four conserved transcription regulators to the larval wing disc region where circular pattern elements develop. The co-localization of Antennapedia, Notch, Distal-less, and Spalt with presumptive (eye)spot organizers was examined in 13 butterfly species, providing the largest comparative dataset available for the system. We found variation between families, between subfamilies, and between tribes. Phylogenetic reconstructions by parsimony and maximum likelihood methods revealed an unambiguous evolutionary history only for Antennapedia, with a resolved single origin of eyespot-associated expression, and many homoplastic events for Notch, Distal-less, and Spalt. The flexibility in the (co-)recruitment of the targeted genes includes cases where different gene combinations are associated with morphologically similar eyespots, as well as cases where identical protein combinations are associated with very different phenotypes.

Conclusions

The evolutionary history of gene (co-)recruitment is consistent with both divergence from a recruited putative ancestral network, and with independent co-option of individual genes. The diversity in the combinations of genes expressed in association with eyespot formation does not parallel diversity in characteristics of the adult phenotype. We discuss these results in the context of inferring homology. Our study underscores the importance of widening the representation of phylogenetic, morphological, and genetic diversity in order to establish general principles about the mechanisms behind the evolution of novel traits.

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.

Wingless and aristaless2 define a developmental ground plan for moth and butterfly wing pattern evolution

Butterfly wing patterns have long been a favorite system for studying the evolutionary radiation of complex morphologies. One of the key characteristics of the system is that wing patterns are based on a highly conserved ground plan of pattern homologies. In fact, the evolution of lepidopteran wing patterns is proposed to have occurred through the repeated gain, loss, and modification of only a handful of serially repeated elements. In this study, we examine the evolution and development of stripe wing pattern elements. We show that expression of the developmental morphogen wingless (wg) is associated with early determination of the major basal (B), discal (DI and DII), and marginal (EI) stripe patterns in a broad sampling of Lepidoptera, suggesting homology of these pattern elements across moths and butterflies. We describe for the first time a novel Lepidoptera-specific homeobox gene, aristaless2 (al2), which precedes wg expression during the early determination of DII stripe patterns. We show that al2 was derived from a tandem duplication of the aristaless gene, whereupon it underwent a rapid coding and cis-regulatory divergence relative to its more conserved paralog aristaless1 (al1), which retained an ancestral expression pattern. The al2 stripe expression domain evolutionarily preceded the appearance of the DII pattern elements in multiple lineages, leading us to speculate that al2 represented preexisting positional information that may have facilitated DII evolution via a developmental drive mechanism. In contrast to butterfly eyespot patterns, which are often cited as a key example of developmental co-option of preexisting developmental genes, this study provides an example where the origin of a major color pattern element is associated with the evolution of a novel lepidopteran homeobox gene.

Local pharmacological effects of tungstate on the color-pattern determination of butterfly wings: a possible relationship between the eyespot and parafocal element

Butterfly wing color patterns can be changed by the application of a temperature shock or pharmacological agents such as tungstate, producing a distinctive type of elemental modification called the TS (temperature shock) type. Heterochronic uncoupling between the signaling and reception steps during the color-pattern determination process has been proposed as a mechanism for TS-type changes. As an extension of this hypothesis, both the parafocal element (PFE) and the eyespot in the same wing compartment are considered to be determined by morphogenic signal(s) emitted from the same eyespot focus. However, these models need to be examined with additional experimental data. Furthermore, there is controversy as to whether the action of tungstate on wing color patterns is direct or Indirect. Using a species of nymphalid butterfly (Junonia orithya), we have devised a simple method for the local application of pharmacological agents directly on developing wings of pupae. Local tungstate application resulted in reduced eyespots and circular dislocated PFEs in the eyespot-less compartments only on the treated wing, demonstrating that tungstate directly induces color-pattern changes on wings. We further examined the eyespot-PFE relationship in normal and cold-shocked Individuals, showing that an eyespot can be superimposed on a PFE and vice versa, probably depending on the timing of their fate determination. Taken together, we propose a two-morphogen model for the normal color-pattern determination, in which the morphogenic signals for the eyespot and PFE are different from each other despite their Identical origin. This two-morphogen model is compatible with the heterochronic uncoupling model for TS-type changes.

The Intersection of Theory and Application in Elucidating Pattern Formation in Developmental Biology

We discuss theoretical and experimental approaches to three distinct developmental systems that illustrate how theory can influence experimental work and vice-versa. The chosen systems - Drosophila melanogaster, bacterial pattern formation, and pigmentation patterns - illustrate the fundamental physical processes of signaling, growth and cell division, and cell movement involved in pattern formation and development. These systems exemplify the current state of theoretical and experimental understanding of how these processes produce the observed patterns, and illustrate how theoretical and experimental approaches can interact to lead to a better understanding of development. As John Bonner said long ago'We have arrived at the stage where models are useful to suggest experiments, and the facts of the experiments in turn lead to new and improved models that suggest new experiments. By this rocking back and forth between the reality of experimental facts and the dream world of hypotheses, we can move slowly toward a satisfactory solution of the major problems of developmental biology.'

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.

Coloration using higher order optical interference in the wing pattern of the Madagascan sunset moth

Colour patterns of animals' bodies are usually produced by the spatial distribution of pigments with different colours. However, some animals use the spatial variation of colour-producing microstructures. We have studied one distinctive example of such structurally produced colour patterns, the wing of the Madagascan sunset moth, to clarify the physical rules that underlie the colour variation. It is known that the iridescent wing scale of the sunset moth has the alternate air-cuticle multilayer structure that causes optical interference. The microscopic and optical investigations of various parts of the wing have confirmed that the thickness of the cuticle layers within the scale largely varies to produce the colour pattern. However, it varies in very different ways between the dorsal and ventral sides of the hind wing; the thickness gradually varies on the dorsal side from scale to scale, while the abrupt changes are found on the ventral side to form distinctive borders between differently coloured areas. It is also revealed that an unusual coloration mechanism is involved in the green part of the ventral hind wing: the colour is caused by higher order optical interference of the highly non-ideal multilayer structure. The physical mechanism of the colour pattern formation is briefly discussed with the several mathematical models proposed so far.

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.

Species-specific color-pattern modifications of butterfly wings

We have previously shown that the systemic injection of sodium tungstate, a general protein-tyrosine phosphatase (PTPase) inhibitor, efficiently produces characteristic color-pattern modifications on the wings of the Painted Lady butterfly, Vanessa cardui. By using this method in the present study, we analyzed modification patterns of six species of Japanese butterflies. Whereas in Vanessa indica the black spots on the forewings reduced in size in response to the treatment, in Lycaena phlaeas the morphologically similar black spots enlarged in size. However, the metallic blue spots on the forewings of V. indica did enlarge in size, showing different behavior even within a single wing surface. The response patterns of Ypthima argus differed markedly from those of other species in that ectopic color-pattern elements were created. Colias erate showed minor modifications that coincidentally resembled the natural color-pattern of a closely related species, Colias palaeno. Through a comprehensive literature search, we confirmed the existence of naturally occurring aberrant color patterns with close similarities to the experimentally induced phenocopies in each of the modified species. Our results point out the possibility that a hypothetical transduction pathway with a PTPase for the scale-cell differentiation globally coordinates the wing-wide color-pattern development in butterflies.

Modelling butterfly wing eyespot patterns

Eyespots are concentric motifs with contrasting colours on butterfly wings. Eyespots have intra- and interspecific visual signalling functions with adaptive and selective roles. We propose a reaction-diffusion model that accounts for eyespot development. The model considers two diffusive morphogens and three non-diffusive pigment precursors. The first morphogen is produced in the focus and determines the differentiation of the first eyespot ring. A second morphogen is then produced, modifying the chromatic properties of the wing background pigment precursor, inducing the differentiation of a second ring. The model simulates the general structural organization of eyespots, their phenotypic plasticity and seasonal variability, and predicts effects from microsurgical manipulations on pupal wings as reported in the literature.

Butterfly Wing Pattern Evolution Is Associated with Changes in a Notch/Distal-less Temporal Pattern Formation Process

In butterflies there is a class of "intervein" wing patterns that have lines of symmetry halfway between wing veins. These patterns occur in a range of shapes, including eyespots, ellipses, and midlines, and were proposed to have evolved through developmental shifts along a midline-to-eyespot continuum. Here we show that Notch (N) upregulation, followed by activation of the transcription factor Distal-less (Dll), is an early event in the development of eyespot and intervein midline patterns across multiple species of butterflies. A relationship between eyespot phenotype and N and Dll expression is demonstrated in a loss-of-eyespot mutant in which N and Dll expression is reduced at missing eyespot sites. A phylogenetic comparison of expression time series from eight moth and butterfly species suggests that intervein N and Dll patterns are a derived characteristic of the butterfly lineage. Furthermore, prior to eyespot determination in eyespot-bearing butterflies, N and Dll are transiently expressed in a pattern that resembles ancestral intervein midline patterns. In this study we establish N upregulation as the earliest known event in eyespot determination, demonstrate gene expression associated with intervein midline color patterns, and provide molecular evidence that wing patterns evolved through addition to and truncation of a conserved midline-to-eyespot pattern formation sequence.

Using mathematical models to help understand biological pattern formation

One of the characteristics of biological systems is their ability to produce and sustain spatial and spatio-temporal pattern. Elucidating the underlying mechanisms responsible for this phenomenon has been the goal of much experimental and theoretical research. This paper illustrates this area of research by presenting some of the mathematical models that have been proposed to account for pattern formation in biology and considering their implications.

Pigmentation pattern formation in butterflies: experiments and models

Butterfly pigmentation patterns are one of the most spectacular and vivid examples of pattern formation in biology. They have attracted much attention from experimentalists and theoreticians, who have tried to understand the underlying genetic, chemical and physical processes that lead to patterning. In this paper, we present a brief review of this field by first considering the generation of the localised, eyespot, patterns and then the formation of more globally controlled patterns. We present some new results applied to pattern formation on the wing of the mimetic butterfly Papilio dardanus.

Development and the Genetics of Evolutionary Change Within Insect Species

▪ Abstract  Changes in genes and in developmental processes generate the phenotypic variation that is sorted by natural selection in adaptive evolution. We review several case studies in which artificial selection experiments in insects have led to divergent morphologies, and where further work has revealed information about the underlying changes at both the genetic and developmental levels. In addition, we examine several studies of phenotypic plasticity where multidisciplinary approaches are also beginning to reveal more about how developmental processes are modulated. Such integrated research will lead to a richer understanding of the changes in development that occur during evolutionary responses to natural selection, and it will also more rigorously examine how developmental processes can influence the tempo and direction of evolutionary change.

Spatial patterns in ant colonies

The origins of large-scale spatial patterns in biology have been an important source of theoretical speculation since the pioneering work by Turing (1952) on the chemical basis of morphogenesis. Knowing how these patterns emerge and their functional role is important to our understanding of the evolution of biocomplexity and the role played by self organization. However, so far, conclusive evidence for local activation-long-range inhibition mechanisms in real biological systems has been elusive. Here a well-defined experimental and theoretical analysis of the pattern formation dynamics exhibited by clustering behavior in ant colonies is presented. These experiments and a simple mathematical model show that these colonies do indeed use this type of mechanism. All microscopic variables have been measured and provide the first evidence, to our knowledge, for this type of self-organized behavior in complex biological systems, supporting early conjectures about its role in the organization of insect societies.

Turing patterns with pentagonal symmetry

We explore numerically the formation of Turing patterns in a confined circular domain with small aspect ratio. Our results show that stable fivefold patterns are formed over a well defined range of disk sizes, offering a possible mechanism for inducing the fivefold symmetry observed in early development of regular echinoids. Using this pattern as a seed, more complex biological structures can be mimicked, such as the pigmentation pattern of sea urchins and the plate arrangements of the calyxes of primitive camerate crinoids.

A model for colour pattern formation in the butterfly wing of Papilio dardanus

The butterfly Papilio dardanus is well known for the spectacular phenotypic polymorphism in the female of the species. We show that numerical simulations of a reaction diffusion model on a geometrically accurate wing domain produce spatial patterns that are consistent with many of those observed on the butterfly. Our results suggest that the wing coloration is due to a simple underlying stripe-like pattern of some pigment-inducing morphogen. We focus on the effect of key factors such as parameter values for mode selection, threshold values which determine colour, wing shape and boundary conditions. The generality of our approach should allow us to investigate other butterfly species. The relationship between these key factors and gene activities is discussed in the context of recent biological advances.

Incorporation of bromodeoxyuridine in regenerating fin tissue of the goldfish Carassius auratus

We have investigated the pattern of incorporation of 5-bromo-2'-deoxyuridine-5'-monophosphate (BrdU) by proliferating cells during regeneration of the tail fin of Carassius auratus. Fifteen days after amputation, intraperitoneal injection of a single dose of 0.25 mg/g wet weight of BrdU and subsequent immunocytochemical detection on sections revealed groups of replicating cells in the blastema and epidermis at different proximodistal levels. Proliferating blastemal cells were confined to a crowded, compact distal area that lost its replicative capacity laterally, causing the differentiation of scleroblasts, which synthesize the lepidotrichia hemisegments. Proximally, but centrally located, the blastemal cells did not incorporate BrdU and they differentiated giving rise to the mature intraray connective tissue. An independent cell-proliferation process was noted in the epidermis. The distal cap did not proliferate; the lateral faces of the epidermis showed high rates of cell replication in the central layer at every level of the regenerate rays; quiescent cells remained in the superficial layers. The basal epidermal cells did not incorporate BrdU when actinotrichia were present. The possible role of basal epidermal cells in the synthesis of actinotrichia, the contribution of these collagen macrofibrils to the morphogenetic process, and the different pathways of cell differentiation during fin regeneration are discussed.

Metaphors and the role of genes in development

In describing the flawless regularity of developmental processes and the correlation between changes at certain genetic loci and changes in morphology, biologists frequently employ two metaphors: that genes 'control' development, and that genomes embody 'programs' for development. Although these metaphors have an admirable sharpness and punch, they lead, when taken literally, to highly distorted pictures of developmental processes. A more balanced, and useful, view of the role of genes in development is that they act as suppliers of the material needs of development and, in some instances, as context-dependent catalysts of cellular changes, rather than as 'controllers' of developmental progress and direction. The consequences of adopting this alternative view of development are discussed.