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

Wound Healing in Butterfly Pupal Wing Tissues: Real-Time In Vivo Imaging of Long-Range Cell Migration, Cluster Formation, and Calcium Oscillations

Insects can repair wounds and regenerate body parts in response to physical damage. Wound healing in butterfly pupal wing tissues is developmentally interesting because ectopic color patterns develop during healing, suggesting that normal and damage-induced color patterns may use similar mechanisms. Here we physiologically investigated wound healing and ectopic color pattern formation in butterfly pupal wing tissues using the blue pansy butterfly Junonia orithya. In response to physical puncture damage, various ectopic color patterns are formed around the damage site. After the wounding operation, we observed hemocytes migrating over long distances along the wing veins (lacunae) toward the damage site, where hemocytes and epidermal cells formed cellular clusters. Calcium oscillations were observed in cells at and near the damage site. Calcium oscillations were transiently affected by ruthenium red, an inhibitor of calcium transporters and channels, and ruthenium red caused various abnormalities in the scales of adult wings. These results suggest that cell migration, cluster formation, and calcium oscillations play important roles in wound healing and scale development at and near the damage site. Ectopic color patterns may develop in response to local calcium oscillations as a consequence of the evolutionary co-option of the healing process for normal development.

Butterfly pupal wing tissue with an eyespot organizer

Butterfly wing eyespots are developmentally determined at the early pupal stage, when prospective eyespot focal cells underneath the pupal cuticle focal spot function as eyespot organizers in the pupal wing tissue. Here, we performed light microscopy and transmission electron microscopy (TEM) to describe cellular structures of pupal wing tissue with an eyespot organizer immediately after pupation using the Blue Pansy butterfly Junonia orithya. The pupal forewing dorsal epidermis was a pseudostratified monolayer of vertically elongated epidermal cells. The apical portion of the cells adhered laterally to one another, but their medial and basal portions were thinner than the apical portion and were tilted to enclose cells at the center, forming a cellular cluster. The cellular cluster at the organizer was relatively large laterally and vertically. The apical portion of the cells and its corresponding cuticle at the organizer were thicker than those in the surroundings. The innermost cuticle layer was being synthesized, indicating high cuticle synthesis and secretion activity of the cells. At the medial and basal portions of the dorsal epidermis, there were many intracellular and extracellular vacuole-like globules, most likely containing extracellular matrix molecules. Some of the basal processes from epidermal cells extended to form protrusions of the basement membrane, which was often attended by hemocytes. These results suggest that the butterfly eyespot organizer is composed of a single or a few cellular clusters that secrete more cuticle than surrounding clusters, supporting the pupal cuticle hypothesis that cuticle formation is critical for eyespot color pattern determination in butterflies.

Socket Array Irregularities and Wing Membrane Distortions at the Eyespot Foci of Butterfly Wings Suggest Mechanical Signals for Color Pattern Determination

Eyespot foci on butterfly wings function as organizers of eyespot color patterns during development. Despite their importance, focal structures have not been examined in detail. Here, we microscopically examined scales, sockets, and the wing membrane in the butterfly eyespot foci of both expanded and unexpanded wings using the Blue Pansy butterfly Junonia orithya. Images from a high-resolution light microscope revealed that, although not always, eyespot foci had scales with disordered planar polarity. Scanning electron microscopy (SEM) images after scale removal revealed that the sockets were irregularly positioned and that the wing membrane was physically distorted as if the focal site were mechanically squeezed from the surroundings. Focal areas without eyespots also had socket array irregularities, but less frequently and less severely. Physical damage in the background area induced ectopic patterns with socket array irregularities and wing membrane distortions, similar to natural eyespot foci. These results suggest that either the process of determining an eyespot focus or the function of an eyespot organizer may be associated with wing-wide mechanics that physically disrupt socket cells, scale cells, and the wing membrane, supporting the physical distortion hypothesis of the induction model for color pattern determination in butterfly wings.

Live Detection of Intracellular Chitin in Butterfly Wing Epithelial Cells In Vivo Using Fluorescent Brightener 28: Implications for the Development of Scales and Color Patterns

Chitin is the major component of the extracellular cuticle and plays multiple roles in insects. In butterflies, chitin builds wing scales for structural colors. Here, we show that intracellular chitin in live cells can be detected in vivo with fluorescent brightener 28 (FB28), focusing on wing epithelial cells of the small lycaenid butterfly Zizeeria maha immediately after pupation. A relatively small number of cells at the apical surface of the epithelium were strongly FB28-positive in the cytosol and seemed to have extensive ER-Golgi networks, which may be specialized chitin-secreting cells. Some cells had FB28-positive tadpole-tail-like or rod-like structures relative to the nucleus. We detected FB28-positive hexagonal intracellular objects and their associated structures extending toward the apical end of the cell, which may be developing scale bases and shafts. We also observed FB28-positive fibrous intracellular structures extending toward the basal end. Many cells were FB28-negative in the cytosol, which contained FB28-positive dots or discs. The present data are crucial to understanding the differentiation of the butterfly wing epithelium, including scale formation and color pattern determination. The use of FB28 in probing intracellular chitin in live cells may be applicable to other insect systems.

Protein Delivery to Insect Epithelial Cells In Vivo: Potential Application to Functional Molecular Analysis of Proteins in Butterfly Wing Development

Protein delivery to cells in vivo has great potential for the functional analysis of proteins in nonmodel organisms. In this study, using the butterfly wing system, we investigated a method of protein delivery to insect epithelial cells that allows for easy access, treatment, and observation in real time in vivo. Topical and systemic applications (called the sandwich and injection methods, respectively) were tested. In both methods, green/orange fluorescent proteins (GFP/OFP) were naturally incorporated into intracellular vesicles and occasionally into the cytosol from the apical surface without any delivery reagent. However, the antibodies were not delivered by the sandwich method at all, and were delivered only into vesicles by the injection method. A membrane-lytic peptide, L17E, appeared to slightly improve the delivery of GFP/OFP and antibodies. A novel peptide reagent, ProteoCarry, successfully promoted the delivery of both GFP/OFP and antibodies into the cytosol via both the sandwich and injection methods. These protein delivery results will provide opportunities for the functional molecular analysis of proteins in butterfly wing development, and may offer a new way to deliver proteins into target cells in vivo in nonmodel organisms.

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.

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.

The cytoneme connection: direct long-distance signal transfer during development

During development, specialized cells produce signals that distribute among receiving cells to induce a variety of cellular behaviors and organize tissues. Recent studies have highlighted cytonemes, a type of specialized signaling filopodia that carry ligands and/or receptor complexes, as having a role in signal dispersion. In this Primer, we discuss how the dynamic regulation of cytonemes facilitates signal transfer in complex environments. We assess recent evidence for the mechanisms for cytoneme formation, function and regulation, and postulate that contact between cytoneme membranes promotes signal transfer as a new type of synapse (morphogenetic synapsis). Finally, we reflect on the fundamental unanswered questions related to understanding cytoneme biology.

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.

Emperors, admirals and giants, zebras, tigers and woolly bears: casting a broader net in exploring heparin effects on Lepidoptera wing patterns

Background: Studies of heparin effects on Lepidoptera wing patterns have been restricted to a small number of species. I report observations from experiments on a broader range of taxa, including first results from swallowtails, tiger moths and microlepidoptera. Methods: Heparin injections were made in prepupae and pupae of Junonia coenia (common buckeyes), Agraulis vanillae (gulf fritillaries), Heliconius charithonia (zebra longwings), Asterocampa clyton (tawny emperors) , Danaus plexippus (monarchs), Vanessa atalanta (red admirals); Heraclides cresphontes (giant swallowtails), Pterourus troilus (spicebush swallowtails), Protographium marcellus (zebra swallowtails), Battus polydamas (polydamas swallowtails); Hypercompe scribonia (giant leopard moths), Estigmene acrea (acrea moths), Hyphantria cunea (fall webworm moths) , Utetheisa ornatrix (ornate bella moths); Glyphodes sibillalis (mulberry leaftier). Results: Heparin sometimes altered the entire pattern in a dramatic way, sometimes caused changes locally. In buckeyes, the previous heparin study conducted on pupae was compared to injections made at a prepupal stage. In gulf fritillaries, zebra longwings and tawny emperors, the dramatic changes occurred throughout their wings, while in monarchs, changes were restricted to wing margins. Changes achieved in red admirals, show that heparin action is unrelated to the original color. In swallowtails, transformations were restricted to border system, indicating higher levels of stability and compartmentalization of wing patterns. In mulberry leaftier, changes were restricted to the marginal bands. In tiger moths, elongation of black markings led to merging of spots; in the ornate bella moth, it was accompanied by an expansion of the surrounding white bands, and results were compared to the effects of colder temperatures. Conclusions: Using pharmaceutical intervention demonstrates that there are many similarities and some very significant differences in the ways wing patterns are formed in different Lepidoptera lineages. By creating a range of variation one can demonstrate how one pattern can easily evolve into another, aiding in understanding of speciation and adaptation processes.

Leopards and giants, tigers and woolly bears: casting a broader net in exploring heparin effects on Lepidoptera wing patterns

Background: Studies of heparin effects on Lepidoptera wing patterns have been restricted to a small number of species. I report observations from experiments on a broader range of taxa, including first results from swallowtails, tiger moths and microlepidoptera. Methods: Heparin injections were made in prepupal and early pupal stages of the following species representing several Lepidoptera lineages: Junonia coenia, Agraulis vanillae, Asterocampa clyton (Nymphalidae); Heraclides cresphontes, Pterourus troilus, Eurytides marcellus (Papilionidae); Hypercompe scribonia, Estigmene acrea, Hyphantria cunea (Erebidae); and Glyphodes sibillalis (Crambidae). Heparin-induced changes in wing pattern are illustrated, and advantages of using prepupal vs. pupal stages for this type of pharmacological manipulation of wing patterns are discussed. Results: In buckeyes, heparin-induced changes consisted of loss of orange parafocal elements as marginal and submarginal bands shifted basally. In gulf fritillaries, changes in black and silver markings were similar to those found in wild aberrant individuals. In tawny emperor, intermediate and extreme levels of transformation were achieved, pointing to homology of this species' unusual eyespots to those of other nymphalids. In swallowtails, heparin-induced changes were much more restricted and proved harder to achieve, possibly indicating higher levels of stability and compartmentalization of wing patterns in this butterfly family. In tiger moths, elongation of black markings occurred so that normally distinct spots sometimes merged; in leopard moth, these changes were restricted to areas adjacent to discal spot. In pyraloid moth, changes were mostly restricted to expansion of wing marginal bands and hindwing border. Conclusions: Variation in wing pattern response to heparin found between different species and families in this study warrants further taxonomic widening of exploration of wing pattern formation mechanisms in Lepidoptera. While there are many similarities, there also seem to be very significant differences in the ways wing patterns are formed in different families of butterflies and moths.

Emperors, admirals and giants, zebras, tigers and woolly bears: casting a broader net in exploring heparin effects on Lepidoptera wing patterns

Background: Studies of heparin effects on Lepidoptera wing patterns have been restricted to a small number of species. I report observations from experiments on a broader range of taxa, including first results from swallowtails, tiger moths and microlepidoptera. Methods: Heparin injections were made in prepupae and pupae of Junonia coenia (common buckeyes), Agraulis vanillae (gulf fritillaries), Heliconius charithonia (zebra longwings), Asterocampa clyton (tawny emperors) , Danaus plexippus (monarchs), Vanessa atalanta (red admirals); Heraclides cresphontes (giant swallowtails), Pterourus troilus (spicebush swallowtails), Protographium marcellus (zebra swallowtails), Battus polydamas (polydamas swallowtails); Hypercompe scribonia (giant leopard moths), Estigmene acrea (acrea moths), Hyphantria cunea (fall webworm moths) , Utetheisa ornatrix (ornate bella moths); Glyphodes sibillalis (mulberry leaftier). Results: Heparin sometimes altered the entire pattern in a dramatic way, sometimes caused changes locally. In buckeyes, the previous heparin study conducted on pupae was compared to injections made at a prepupal stage. In gulf fritillaries, zebra longwings and tawny emperors, the dramatic changes occurred throughout their wings, while in monarchs, changes were restricted to wing margins. Changes achieved in red admirals, show that heparin action is unrelated to the original color. In swallowtails, transformations were restricted to border system, indicating higher levels of stability and compartmentalization of wing patterns. In mulberry leaftier, changes were restricted to the marginal bands. In tiger moths, elongation of black markings led to merging of spots; in the ornate bella moth, it was accompanied by an expansion of the surrounding white bands, and results were compared to the effects of colder temperatures. Conclusions: Using pharmaceutical intervention demonstrates that there are many similarities and some very significant differences in the ways wing patterns are formed in different Lepidoptera lineages. By creating a range of variation one can demonstrate how one pattern can easily evolve into another, aiding in understanding of speciation and adaptation processes.

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.

Wound healing, calcium signaling, and other novel pathways are associated with the formation of butterfly eyespots

Background

One hypothesis surrounding the origin of novel traits is that they originate from the co-option of pre-existing genes or larger gene regulatory networks into novel developmental contexts. Insights into a trait’s evolutionary origins can, thus, be gained via identification of the genes underlying trait development, and exploring whether those genes also function in other developmental contexts. Here we investigate the set of genes associated with the development of eyespot color patterns, a trait that originated once within the Nymphalid family of butterflies. Although several genes associated with eyespot development have been identified, the eyespot gene regulatory network remains largely unknown.

Results

In this study, next-generation sequencing and transcriptome analyses were used to identify a large set of genes associated with eyespot development of Bicyclus anynana butterflies, at 3-6 h after pupation, prior to the differentiation of the color rings. Eyespot-associated genes were identified by comparing the transcriptomes of homologous micro-dissected wing tissues that either develop or do not develop eyespots in wild-type and a mutant line of butterflies, Spotty, with extra eyespots. Overall, 186 genes were significantly up and down-regulated in wing tissues that develop eyespots compared to wing tissues that do not. Many of the differentially expressed genes have yet to be annotated. New signaling pathways, including the Toll, Fibroblast Growth Factor (FGF), extracellular signal–regulated kinase (ERK) and/or Jun N-terminal kinase (JNK) signaling pathways are associated for the first time with eyespot development. In addition, several genes involved in wound healing and calcium signaling were also found to be associated with eyespots.

Conclusions

Overall, this study provides the identity of many new genes and signaling pathways associated with eyespots, and suggests that the ancient wound healing gene regulatory network may have been co-opted to cells at the center of the pattern to aid in eyespot origins. New transcription factors that may be providing different identities to distinct wing sectors, and genes with sexually dimorphic expression in the eyespots were also identified.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4175-7) contains supplementary material, which is available to authorized users.