Structural analysis of eyespots: dynamics of morphogenic signals that govern elemental positions in butterfly wings

Phenotypic Plasticity of the Mimetic Swallowtail Butterfly Papilio polytes: Color Pattern Modifications and Their Implications in Mimicry Evolution

Simple Summary

Diverse butterfly wing color patterns are evolutionary products in response to environmental changes in the past. Environmental stress, such as temperature shock, is known to induce color pattern modifications in various butterfly species, and this phenotypic plasticity plays an important role in color pattern evolution. However, the potential contributions of phenotypic plasticity to mimicry evolution have not been evaluated. Here, we focused on the swallowtail butterfly Papilio polytes, which has nonmimetic and mimetic forms in females, to examine its plastic phenotypes. Cold shock and heat shock treatments in the nonmimetic form induced color pattern modifications that were partly similar to those of the mimetic form, and nonmimetic females were more sensitive than males and mimetic females. These results suggest that phenotypic plasticity in nonmimetic females might have provided a basis of natural selection for mimetic color patterns during evolution.

Abstract

Butterfly wing color patterns are sensitive to environmental stress, such as temperature shock, and this phenotypic plasticity plays an important role in color pattern evolution. However, the potential contributions of phenotypic plasticity to mimicry evolution have not been evaluated. Here, we focused on the swallowtail butterfly Papilio polytes, which has nonmimetic and mimetic forms in females, to examine its plastic phenotypes. In the nonmimetic form, medial white spots and submarginal reddish spots in the ventral hindwings were enlarged by cold shock but were mostly reduced in size by heat shock. These temperature-shock-induced color pattern modifications were partly similar to mimetic color patterns, and nonmimetic females were more sensitive than males and mimetic females. Unexpectedly, injection of tungstate, a known modification inducer in nymphalid and lycaenid butterflies, did not induce any modification, but fluorescent brightener 28, another inducer discovered recently, induced unique modifications. These results suggest that phenotypic plasticity in nonmimetic females might have provided a basis of natural selection for mimetic color patterns during evolution.

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.

Scale granules and colours: Sexual dimorphism in Trichonis (Lepidoptera: Lycaenidae, Theclinae)

A large fraction of dorsal wing surface ground scales show an unusual granulated nature, composed of material apparently extruded from the scale lumen in male individuals of both Trichonis Hewitson, 1865 species in the tribe Eumaeini, a rare Guyanian-Amazonian genus. Only a few not-granulated male specimens are known, females are not granulated. The granulated scales are investigated by various microscopic (optical, scanning and transmission electron microscopy, focused ion beam lamella cutting) and spectroscopic (optical reflectance, energy-dispersive X-ray (EDS), Raman) techniques. The characteristic blue colour unique in the South American representatives of the tribe is documented and analysed. EDS spectra show that the granules contain additional calcium and oxygen as compared with the un-granulated regions of the same scale. Electron diffraction (inside the TEM) did not reveal any crystalline component in the granules. The granulated wing surfaces of the males exhibit a UV absorption band at 280 nm, characteristic for biogenic CaCO₃; therefore, the material of the granules is tentatively identified as CaCO₃. It is shown that the granules influence the optical properties of the dorsal wing surface resulting in a characteristic spectrum.

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

Simple Summary

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

Abstract

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

Speciation in North American Junonia from a genomic perspective

Delineating species boundaries in phylogenetic groups undergoing recent radiation is a daunting challenge akin to discretizing continuity. Here, we propose a general approach exemplified by American butterflies from the genus Junonia Hübner notorious for the variety of similar phenotypes, ease of hybridization, and the lack of consensus about their classification. We obtain whole-genome shotgun sequences of about 200 specimens. We reason that discreteness emerges from continuity by means of a small number of key players, and search for the proteins that diverged markedly between sympatric populations of different species, while keeping low polymorphism within these species. Being 0.25% of the total number, these three dozen 'speciation' proteins indeed partition pairs of Junonia populations into two clusters with a prominent break in between, while all proteins taken together fail to reveal this discontinuity. Populations with larger divergence from each other, comparable to that between two sympatric species, form the first cluster and correspond to different species. The other cluster is characterized by smaller divergence, similar to that between allopatric populations of the same species and comprise conspecific pairs. Using this method, we conclude that J. genoveva (Cramer), J. litoralis Brévignon, J. evarete (Cramer), and J. divaricata C. & R. Felder are restricted to South America. We find that six species of Junonia are present in the United States, one of which is new: Junonia stemosa Grishin, sp.n. (i), found in south Texas and phenotypically closest to J. nigrosuffusa W. Barnes & McDunnough (ii) in its dark appearance. In the pale nudum of the antennal club, these two species resemble J. zonalis C. & R. Felder (iii) from Florida and the Caribbean Islands. The pair of sister species, J. grisea Austin & J. Emmel (iv) and J. coenia Hübner (v), represent the classic west/east U.S.A. split. The mangrove feeder (as caterpillar), dark nudum J. neildi Brévignon (vi) enters south Texas as a new subspecies Junonia neildi varia Grishin ssp.n. characterized by more extensive hybridization with and introgression from J. coenia, and, as a consequence, more variable wing patterns compared with the nominal J. n. neildi in Florida. Furthermore, a new mangrove-feeding species from the Pacific Coast of Mexico is described as Junonia pacoma Grishin sp.n. Finally, genomic analysis suggests that J. nigrosuffusa may be a hybrid species formed by the ancestors of J. grisea and J. stemosa sp.n.

Morphological and Spatial Diversity of the Discal Spot on the Hindwings of Nymphalid Butterflies: Revision of the Nymphalid Groundplan

Simple Summary

Butterfly wing color patterns are diverse, but they can be understood as modifications of the common scheme called the nymphalid groundplan. The discal spot is relatively small, but it is one of the important components of the nymphalid groundplan. Using many hindwing specimens of the family Nymphalidae, the morphological and spatial diversity of the discal spot was studied. The discal spot is expressed as a small or narrow spot, a pair of parallel bands, a diamond or oval structure, a large dark spot, a few fragmented spots, or a white structure. The discal spot is always located in a central portion of the wing defined by the wing veins, and this portion is sandwiched by a pair of bands of the central symmetry system, another important component of the nymphalid groundplan. On the basis of these results, the present study revises the nymphalid groundplan in minor points; the discal spot is an independent and diverse miniature symmetry system nested within the central symmetry system. Due to the involvement of wing veins to define the locations of the discal spot, the present study suggests the possible developmental dynamics of butterfly color pattern formation that produces color pattern diversity.

Abstract

Diverse butterfly wing color patterns are understood through the nymphalid groundplan, which mainly consists of central, border, and basal symmetry systems and a discal spot. However, the status of the discal spot remains unexplored. Here, the morphological and spatial diversity of the discal spot was studied in nymphalid hindwings. The discal spot is expressed as a small or narrow spot, a pair of parallel bands, a diamond or oval structure, a large dark spot, a few fragmented spots, or a white structure. In some cases, the discal spot is morphologically similar to and integrated with the central symmetry system (CSS). The discal spot is always located in a distal portion of the discal cell defined by the wing veins, which is sandwiched by the distal and proximal bands of the CSS (dBC and pBC) and is rarely occupied by border ocelli. The CSS occasionally has the central band (cBC), which differs from the discal spot. These results suggest that the discal spot is an independent and diverse miniature symmetry system nested within the CSS and that the locations of the discal spot and the CSS are determined by the wing veins at the early stage of wing development.

Butterfly eyespot color pattern formation requires physical contact of the pupal wing epithelium with extracellular materials for morphogenic signal propagation

BACKGROUND

Eyespot color pattern formation on butterfly wings is sensitive to physical damage and physical distortion as well as physical contact with materials on the surface of wing epithelial tissue at the pupal stage. Contact-mediated eyespot color pattern changes may imply a developmental role of the extracellular matrix in morphogenic signal propagation. Here, we examined eyespot responses to various contact materials, focusing on the hindwing posterior eyespots of the blue pansy butterfly, Junonia orithya.

RESULTS

Contact with various materials, including both nonbiological and biological materials, induced eyespot enlargement, reduction, or no change in eyespot size, and each material was characterized by a unique response profile. For example, silicone glassine paper almost always induced a considerable reduction, while glass plates most frequently induced enlargement, and plastic plates generally produced no change. The biological materials tested here (fibronectin, polylysine, collagen type I, and gelatin) resulted in various responses, but polylysine induced more cases of enlargement, similar to glass plates. The response profile of the materials was not readily predictable from the chemical composition of the materials but was significantly correlated with the water contact angle (water repellency) of the material surface, suggesting that the surface physical chemistry of materials is a determinant of eyespot size. When the proximal side of a prospective eyespot was covered with a size-reducing material (silicone glassine paper) and the distal side and the organizer were covered with a material that rarely induced size reduction (plastic film), the proximal side of the eyespot was reduced in size in comparison with the distal side, suggesting that signal propagation but not organizer activity was inhibited by silicone glassine paper.

CONCLUSIONS

These results suggest that physical contact with an appropriate hydrophobic surface is required for morphogenic signals from organizers to propagate normally. The binding of the apical surface of the epithelium with an opposing surface may provide mechanical support for signal propagation. In addition to conventional molecular morphogens, there is a possibility that mechanical distortion of the epithelium that is propagated mechanically serves as a nonmolecular morphogen to induce subsequent molecular changes, in accordance with the distortion hypothesis for butterfly wing color pattern formation.

Insights into eyespot color-pattern formation mechanisms from color gradients, boundary scales, and rudimentary eyespots in butterfly wings

Butterfly eyespot color patterns are traditionally explained by the gradient model, where positional information is stably provided by a morphogen gradient from a single organizer and its output is a set of non-graded (or graded) colors based on pre-determined threshold levels. An alternative model is the induction model, in which the outer black ring and the inner black core disk of an eyespot are specified by graded signals from the primary and secondary organizers that also involve lateral induction. To examine the feasibility of these models, we analyzed eyespot color gradients, boundary scales, and rudimentary eyespots in various nymphalid butterflies. Most parts of eyespots showed color gradients with gradual or fluctuating changes with sharp boundaries in many species, but some species had eyespots that were composed of a constant color within a given part. Thus, a plausible model should be flexible enough to incorporate this diversity. Some boundary scales appeared to have two kinds of pigments, and others had "misplaced" colors, suggesting an overlapping of two signals and a difficulty in assuming sharp threshold boundaries. Rudimentary eyespots of three Junonia species revealed that the outer black ring is likely determined first and the inner yellow or red ring is laterally induced. This outside-to-inside determination together with the lateral induction may favor the induction model, in which dynamic signal interactions play a major role. The implications of these results for the ploidy hypothesis and color-pattern rules are discussed.

Real-Time In Vivo Imaging of the Developing Pupal Wing Tissues in the Pale Grass Blue Butterfly Zizeeria maha: Establishing the Lycaenid System for Multiscale Bioimaging

To systematically analyze biological changes with spatiotemporal dynamics, it is important to establish a system that is amenable for real-time in vivo imaging at various size levels. Herein, we focused on the developing pupal wing tissues in the pale grass blue butterfly, Zizeeria maha, as a system of choice for a systematic multiscale approach in vivo in real time. We showed that the entire pupal wing could be monitored throughout development using a high-resolution bright-field time-lapse imaging system under the forewing-lift configuration; we recorded detailed dynamics of the dorsal and ventral epithelia that behaved independently for peripheral adjustment. We also monitored changes in the dorsal hindwing at the compartmental level and directly observed evaginating scale buds. We also employed a confocal laser microscopy system with multiple fluorescent dyes for three-dimensional observations at the tissue and cellular levels. We discovered extensive cellular clusters that may be functionally important as a unit of cellular communication and differentiation. We also identified epithelial discal and marginal dents that may function during development. Together, this lycaenid forewing system established a foundation to study the differentiation process of epithelial cells and can be used to study biophysically challenging mechanisms such as the determination of color patterns and scale nanoarchitecture at the multiscale levels.

Long-Range Effects of Wing Physical Damage and Distortion on Eyespot Color Patterns in the Hindwing of the Blue Pansy Butterfly Junonia orithya

Butterfly eyespot color patterns have been studied using several different approaches, including applications of physical damage to the forewing. Here, damage and distortion experiments were performed, focusing on the hindwing eyespots of the blue pansy butterfly Junonia orithya. Physical puncture damage with a needle at the center of the eyespot reduced the eyespot size. Damage at the eyespot outer rings not only deformed the entire eyespot, but also diminished the eyespot core disk size, despite the distance from the damage site to the core disk. When damage was inflicted near the eyespot, the eyespot was drawn toward the damage site. The induction of an ectopic eyespot-like structure and its fusion with the innate eyespots were observed when damage was inflicted in the background area. When a small stainless ball was placed in close proximity to the eyespot using the forewing-lift method, the eyespot deformed toward the ball. Taken together, physical damage and distortion elicited long-range inhibitory, drawing (attracting), and inducing effects, suggesting that the innate and induced morphogenic signals travel long distances and interact with each other. These results are consistent with the distortion hypothesis, positing that physical distortions of wing tissue contribute to color pattern determination in butterfly wings.

Developmental dynamics of butterfly wings: real-time in vivo whole-wing imaging of twelve butterfly species

Colour pattern development of butterfly wings has been studied from several different approaches. However, developmental changes in the pupal wing tissues have rarely been documented visually. In this study, we recorded real-time developmental changes of the pupal whole wings of 9 nymphalid, 2 lycaenid, and 1 pierid species in vivo, from immediately after pupation to eclosion, using the forewing-lift method. The developmental period was roughly divided into four sequential stages. At the very early stage, the wing tissue was transparent, but at the second stage, it became semi-transparent and showed dynamic peripheral adjustment and slow low-frequency contractions. At this stage, the wing peripheral portion diminished in size, but simultaneously, the ventral epithelium expanded in size. Likely because of scale growth, the wing tissue became deeply whitish at the second and third stages, followed by pigment deposition and structural colour expression at the fourth stage. Some red or yellow (light-colour) areas that emerged early were “overpainted” by expanding black areas, suggesting the coexistence of two morphogenic signals in some scale cells. The discal spot emerged first in some nymphalid species, as though it organised the entire development of colour patterns. These results indicated the dynamic wing developmental processes common in butterflies.

Butterfly eyespot organiser: in vivo imaging of the prospective focal cells in pupal wing tissues

Butterfly wing eyespot patterns are determined in pupal tissues by organisers located at the centre of the prospective eyespots. Nevertheless, organiser cells have not been examined cytochemically in vivo, partly due to technical difficulties. Here, we directly observed organiser cells in pupal forewing epithelium via an in vivo confocal fluorescent imaging technique, using 1-h post-pupation pupae of the blue pansy butterfly, Junonia orithya. The prospective eyespot centre was indented from the plane of the ventral tissue surface. Three-dimensional reconstruction images revealed that the apical portion of “focal cells” at the bottom of the eyespot indentation contained many mitochondria. The mitochondrial portion was connected with a “cell body” containing a nucleus. Most focal cells had globular nuclei and were vertically elongated, but cells in the wing basal region had flattened nuclei and were tilted toward the distal direction. Epithelial cells in any wing region had cytoneme-like horizontal processes. From 1 h to 10 h post-pupation, nuclear volume increased, suggesting DNA synthesis during this period. Morphological differences among cells in different regions may suggest that organiser cells are developmentally ahead of cells in other regions and that position-dependent heterochronic development is a general mechanism for constructing colour patterns in butterfly wings.

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

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

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.

The marginal band system in nymphalid butterfly wings

Butterfly wing color patterns are highly complex and diverse, but they are believed to be derived from the nymphalid groundplan, which is composed of several color pattern systems. Among these pattern systems, the marginal band system, including marginal and submarginal bands, has rarely been studied. Here, we examined the color pattern diversity of the marginal band system among nymphalid butterflies. Marginal and submarginal bands are usually expressed as a pair of linear bands aligned with the wing margin. However, a submarginal band can be expressed as a broken band, an elongated oval, or a single dot. The marginal focus, usually a white dot at the middle of a wing compartment along the wing edge, corresponds to the pupal edge spot, one of the pupal cuticle spots that signify the locations of color pattern organizing centers. A marginal band can be expressed as a semicircle, an elongated oval, or a pair of eyespot-like structures, which suggest the organizing activity of the marginal focus. Physical damage at the pupal edge spot leads to distal dislocation of the submarginal band in Junonia almana and in Vanessa indica, suggesting that the marginal focus functions as an organizing center for the marginal band system. Taken together, we conclude that the marginal band system is developmentally equivalent to other symmetry systems. Additionally, the marginal band is likely a core element and the submarginal band a paracore element of the marginal band system, and both bands are primarily specified by the marginal focus organizing center.

A Single-Wing Removal Method to Assess Correspondence Between Gene Expression and Phenotype in Butterflies: The Case of Distal-less

It is often desirable but difficult to retrieve information on the mature phenotype of an immature tissue sample that has been subjected to gene expression analysis. This problem cannot be ignored when individual variation within a species is large. To circumvent this problem in the butterfly wing system, we developed a new surgical method for removing a single forewing from a pupa using Junonia orithya; the operated pupa was left to develop to an adult without eclosion. The removed right forewing was subjected to gene expression analysis, whereas the non-removed left forewing was examined for color patterns. As a test case, we focused on Distal-less (Dll), which likely plays an active role in inducing elemental patterns, including eyespots. The Dll expression level in forewings was paired with eyespot size data from the same individual. One third of the operated pupae survived and developed wing color patterns. Dll expression levels were significantly higher in males than in females, although male eyespots were smaller in size than female eyespots. Eyespot size data showed weak but significant correlations with the Dll expression level in females. These results demonstrate that a single-wing removal method was successfully applied to the butterfly wing system and suggest the weak and non-exclusive contribution of Dll to eyespot size determination in this butterfly. Our novel methodology for establishing correspondence between gene expression and phenotype can be applied to other candidate genes for color pattern development in butterflies. Conceptually similar methods may also be applicable in other developmental systems.

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.

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.

Distal-less induces elemental color patterns in Junonia butterfly wings

BACKGROUND

The border ocellus, or eyespot, is a conspicuous color pattern element in butterfly wings. For two decades, it has been hypothesized that transcription factors such as Distal-less (Dll) are responsible for eyespot pattern development in butterfly wings, based on their expression in the prospective eyespots. In particular, it has been suggested that Dll is a determinant for eyespot size. However, functional evidence for this hypothesis has remained incomplete, due to technical difficulties.

RESULTS

Here, we show that ectopically expressed Dll induces ectopic elemental color patterns in the adult wings of the blue pansy butterfly, Junonia orithya (Lepidoptera, Nymphalidae). Using baculovirus-mediated gene transfer, we misexpressed Dll protein fused with green fluorescent protein (GFP) in pupal wings, resulting in ectopic color patterns, but not the formation of intact eyespots. Induced changes included clusters of black and orange scales (a basic feature of eyespot patterns), black and gray scales, and inhibition of cover scale development. In contrast, ectopic expression of GFP alone did not induce any color pattern changes using the same baculovirus-mediated gene transfer system.

CONCLUSIONS

These results suggest that Dll plays an instructive role in the development of color pattern elements in butterfly wings, although Dll alone may not be sufficient to induce a complete eyespot. This study thus experimentally supports the hypothesis of Dll function in eyespot development.

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.

The Lycaenid Central Symmetry System: Color Pattern Analysis of the Pale Grass Blue Butterfly Zizeeria maha

The nymphalid groundplan has been proposed to explain diverse butterfly wing color patterns. In this model, each symmetry system is composed of a core element and a pair of paracore elements. The development of this elemental configuration has been explained by the induction model for positional information. However, the diversity of color patterns in other butterfly families in relation to the nymphalid groundplan has not been thoroughly examined. Here, we examined aberrant color pattern phenotypes of a lycaenid butterfly, Zizeeria maha, from mutagenesis and plasticity studies as well as from field surveys. In several mutants, the third and fourth spot arrays were coordinately positioned much closer to the discal spot in comparison to the normal phenotype. In temperature-shock types, the third and fourth array spots were elongated inwardly or outwardly from their normal positions. In field-caught spontaneous mutants, small black spots were located adjacent to normal black spots. Analysis of these aberrant phenotypes indicated that the spots belonging to the third and fourth arrays are synchronously changeable in position and shape around the discal spot. Thus, these arrays constitute paracore elements of the central symmetry system of the lycaenid butterflies, and the discal spot comprises the core element. These aberrant phenotypes can be explained by the black-inducing signals that propagate from the prospective discal spot, as predicted by the induction model. These results suggest the existence of long-range developmental signals that cover a large area of a wing not only in nymphalid butterflies, but also in lycaenid butterflies.

Live Cell Imaging of Butterfly Pupal and Larval Wings In Vivo

Butterfly wing color patterns are determined during the late larval and early pupal stages. Characterization of wing epithelial cells at these stages is thus critical to understand how wing structures, including color patterns, are determined. Previously, we successfully recorded real-time in vivo images of developing butterfly wings over time at the tissue level. In this study, we employed similar in vivo fluorescent imaging techniques to visualize developing wing epithelial cells in the late larval and early pupal stages 1 hour post-pupation. Both larval and pupal epithelial cells were rich in mitochondria and intracellular networks of endoplasmic reticulum, suggesting high metabolic activities, likely in preparation for cellular division, polyploidization, and differentiation. Larval epithelial cells in the wing imaginal disk were relatively large horizontally and tightly packed, whereas pupal epithelial cells were smaller and relatively loosely packed. Furthermore, larval cells were flat, whereas pupal cells were vertically elongated as deep as 130 μm. In pupal cells, many endosome-like or autophagosome-like structures were present in the cellular periphery down to approximately 10 μm in depth, and extensive epidermal feet or filopodia-like processes were observed a few micrometers deep from the cellular surface. Cells were clustered or bundled from approximately 50 μm in depth to deeper levels. From 60 μm to 80 μm in depth, horizontal connections between these clusters were observed. The prospective eyespot and marginal focus areas were resistant to fluorescent dyes, likely because of their non-flat cone-like structures with a relatively thick cuticle. These in vivo images provide important information with which to understand processes of epithelial cell differentiation and color pattern determination in butterfly wings.

Color pattern evolution in Vanessa butterflies (Nymphalidae: Nymphalini): non-eyespot characters

A phylogenetic approach was used to study color pattern evolution in Vanessa butterflies. Twenty-four color pattern elements from the Nymphalid ground plan were identified on the dorsal and ventral surfaces of the fore- and hind wings. Eyespot characters were excluded and will be examined elsewhere. The evolution of each character was traced over a Bayesian phylogeny of Vanessa reconstructed from 7750 DNA base pairs from 10 genes. Generally, the correspondence between character states on the same surface of the two wings is stronger on the ventral side compared to the dorsal side. The evolution of character states on both sides of a wing correspond with each other in most extant species, but the correspondence between dorsal and ventral character states is much stronger in the forewing than in the hindwing. The dorsal hindwing of many species of Vanessa is covered with an extended Basal Symmetry System and the Discalis I pattern element is highly variable between species, making this wing surface dissimilar to the other wing surfaces. The Basal Symmetry System and Discalis I may contribute to behavioral thermoregulation in Vanessa. Overall, interspecific directional character state evolution of non-eyespot color patterns is relatively rare in Vanessa, with a majority of color pattern elements showing non-variable, non-directional, or ambiguous character state evolution. The ease with which the development of color patterns can be modified, including character state reversals, has likely made important contributions to the production of color pattern diversity in Vanessa and other butterfly groups.

Spontaneous long-range calcium waves in developing butterfly wings

Background

Butterfly wing color patterns emerge as the result of a regular arrangement of scales produced by epithelial scale cells at the pupal stage. These color patterns and scale arrangements are coordinated throughout the wing. However, the mechanism by which the development of scale cells is controlled across the entire wing remains elusive. In the present study, we used pupal wings of the blue pansy butterfly, Junonia orithya, which has distinct eyespots, to examine the possible involvement of Ca²⁺ waves in wing development.

Results

Here, we demonstrate that the developing pupal wing tissue of the blue pansy butterfly displayed spontaneous low-frequency Ca²⁺ waves in vivo that propagated slowly over long distances. Some waves appeared to be released from the immediate peripheries of the prospective eyespot and discal spot, though it was often difficult to identify the specific origins of these waves. Physical damage, which is known to induce ectopic eyespots, led to the radiation of Ca²⁺ waves from the immediate periphery of the damaged site. Thapsigargin, which is a specific inhibitor of Ca²⁺-ATPases in the endoplasmic reticulum, induced an acute increase in cytoplasmic Ca²⁺ levels and halted the spontaneous Ca²⁺ waves. Additionally, thapsigargin-treated wings showed incomplete scale development as well as other scale and color pattern abnormalities.

Conclusions

We identified a novel form of Ca²⁺ waves, spontaneous low-frequency slow waves, which travel over exceptionally long distances. Our results suggest that spontaneous Ca²⁺ waves play a critical role in the coordinated development of scale arrangements and possibly in color pattern formation in butterflies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-015-0067-8) contains supplementary material, which is available to authorized users.

Real-time in vivo imaging of butterfly wing development: revealing the cellular dynamics of the pupal wing tissue

Butterfly wings are covered with regularly arranged single-colored scales that are formed at the pupal stage. Understanding pupal wing development is therefore crucial to understand wing color pattern formation. Here, we successfully employed real-time in vivo imaging techniques to observe pupal hindwing development over time in the blue pansy butterfly, Junonia orithya. A transparent sheet of epithelial cells that were not yet regularly arranged was observed immediately after pupation. Bright-field imaging and autofluorescent imaging revealed free-moving hemocytes and tracheal branches of a crinoid-like structure underneath the epithelium. The wing tissue gradually became gray-white, epithelial cells were arranged regularly, and hemocytes disappeared, except in the bordering lacuna, after which scales grew. The dynamics of the epithelial cells and scale growth were also confirmed by fluorescent imaging. Fluorescent in vivo staining further revealed that these cells harbored many mitochondria at the surface of the epithelium. Organizing centers for the border symmetry system were apparent immediately after pupation, exhibiting a relatively dark optical character following treatment with fluorescent dyes, as well as in autofluorescent images. The wing tissue exhibited slow and low-frequency contraction pulses with a cycle of approximately 10 to 20 minutes, mainly occurring at 2 to 3 days postpupation. The pulses gradually became slower and weaker and eventually stopped. The wing tissue area became larger after contraction, which also coincided with an increase in the autofluorescence intensity that might have been caused by scale growth. Examination of the pattern of color development revealed that the black pigment was first deposited in patches in the central areas of an eyespot black ring and a parafocal element. These results of live in vivo imaging that covered wide wing area for a long time can serve as a foundation for studying the cellular dynamics of living wing tissues in butterflies.

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

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