Spalt expression and the development of melanic color patterns in pierid butterflies

Contrasting Size Distributions among the Wing Spots of a Pierid Butterfly Suggest Different Selective Histories

Males of the small cabbage butterfly Pieris rapae crucivora have two dark or melanic spots in the central white area of each dorsal forewing, an anterior spot (aS) and a posterior spot (pS). We used characteristics of the size distributions of these spots to infer how selection has acted on them during their evolution. Our study reveals that the aS size distribution is normal while that of pS is very right-skewed. Moreover, aS size is larger and less variable than pS size. These results suggest that the aS has been under stabilizing selection while the pS has not. The context in which this selection is acting is not yet clear, nor is it clear why pS persists as a wing marking.

The genetic basis of wing spots in Pieris canidia butterflies

Spots in pierid butterflies and eyespots in nymphalid butterflies are likely non-homologous wing colour pattern elements, yet they share a few features in common. Both develop black scales that depend on the function of the gene spalt, and both might have central signalling cells. This suggests that both pattern elements may be sharing common genetic circuitry. Hundreds of genes have already been associated with the development of nymphalid butterfly eyespot patterns, but the genetic basis of the simpler spot patterns on the wings of pierid butterflies has not been investigated. To facilitate studies of pierid wing patterns, we report a high-quality draft genome assembly for Pieris canidia, the Indian cabbage white. We then conducted transcriptomic analyses of pupal wing tissues sampled from the spot and non-spot regions of P. canidia at 3-6 h post-pupation. A total of 1352 genes were differentially regulated between wing tissues with and without the black spot, including spalt, Krüppel-like factor 10, genes from the Toll, Notch, TGF-β, and FGFR signalling pathways, and several genes involved in the melanin biosynthetic pathway. We identified 14 genes that are up-regulated in both pierid spots and nymphalid eyespots and propose that spots and eyespots share regulatory modules despite their likely independent origins.

The Genetic Mechanisms Underlying the Concerted Expression of the yellow and tan Genes in Complex Patterns on the Abdomen and Wings of Drosophila guttifera

How complex morphological patterns form is an intriguing question in developmental biology. However, the mechanisms that generate complex patterns remain largely unknown. Here, we sought to identify the genetic mechanisms that regulate the tan (t) gene in a multi-spotted pigmentation pattern on the abdomen and wings of Drosophila guttifera. Previously, we showed that yellow (y) gene expression completely prefigures the abdominal and wing pigment patterns of this species. In the current study, we demonstrate that the t gene is co-expressed with the y gene in nearly identical patterns, both transcripts foreshadowing the adult abdominal and wing melanin spot patterns. We identified cis-regulatory modules (CRMs) of t, one of which drives reporter expression in six longitudinal rows of spots on the developing pupal abdomen, while the second CRM activates the reporter gene in a spotted wing pattern. Comparing the abdominal spot CRMs of y and t, we found a similar composition of putative transcription factor binding sites that are thought to regulate the complex expression patterns of both terminal pigmentation genes y and t. In contrast, the y and t wing spots appear to be regulated by distinct upstream factors. Our results suggest that the D. guttifera abdominal and wing melanin spot patterns have been established through the co-regulation of y and t, shedding light on how complex morphological traits may be regulated through the parallel coordination of downstream target genes.

Distal-less and spalt are distal organisers of pierid wing patterns

Two genes, Distal-less (Dll) and spalt (sal), are known to be involved in establishing nymphalid butterfly wing patterns. They function in several ways: in the differentiation of the eyespot’s central signalling cells, or foci; in the differentiation of the surrounding black disc; in overall scale melanisation (Dll); and in elaborating marginal patterns, such as parafocal elements. However, little is known about the functions of these genes in the development of wing patterns in other butterfly families. Here, we study the expression and function of Dll and sal in the development of spots and other melanic wing patterns of the Indian cabbage white, Pieris canidia, a pierid butterfly. In P. canidia, both Dll and Sal proteins are expressed in the scale-building cells at the wing tips, in chevron patterns along the pupal wing margins, and in areas of future scale melanisation. Additionally, Sal alone is expressed in the future black spots. CRISPR knockouts of Dll and sal showed that each gene is required for the development of melanic wing pattern elements, and repressing pteridine granule formation, in the areas where they are expressed. We conclude that both genes likely play ancestral roles in organising distal butterfly wing patterns, across pierid and nymphalid butterflies, but are unlikely to be differentiating signalling centres in pierids black spots. The genetic and developmental mechanisms that set up the location of spots and eyespots are likely distinct in each lineage.

Transcription factors underlying wing margin color patterns and pupal cuticle markings in butterflies

Background

The diversity of butterfly color patterns can be attributed to a relatively small number of pattern elements that are homologous across Lepidoptera. Although genes involved in patterning some of these elements have been identified, the development of several major elements remains poorly understood. To identify genes underlying wing pupal cuticle markings and wing margin color patterns, we examined expression of the candidate transcription factors Engrailed/Invected (En/Inv), Distal-less (Dll), Cubitus interruptus (Ci), and Spalt in two nymphalids: Junonia coenia and Bicyclus anynana.

Results

We found that En/Inv, Dll, and Ci mark domains on the J. coenia last-instar forewing disc that closely correspond to the position and shape of pupal cuticle markings. We also found that Spalt demarcates wing margin color patterns in both J. coenia and B. anynana, and that CRISPR/Cas9 deletions in the spalt gene result in reduction and loss of wing margin color patterns in J. coenia. These data demonstrate a role for spalt in promoting wing margin color patterning, in addition to its previously described role in eyespot patterning.

Conclusion

Our observations support the model that a core set of regulatory genes are redeployed multiple times, and in multiple roles, during butterfly wing pattern development. Of these genes, spalt is of special interest as it plays a dual role in both eyespot and margin color pattern development.

Origin of the mechanism of phenotypic plasticity in satyrid butterfly eyespots

Plasticity is often regarded as a derived adaptation to help organisms survive in variable but predictable environments, however, we currently lack a rigorous, mechanistic examination of how plasticity evolves in a large comparative framework. Here, we show that phenotypic plasticity in eyespot size in response to environmental temperature observed in Bicyclus anynana satyrid butterflies is a complex derived adaptation of this lineage. By reconstructing the evolution of known physiological and molecular components of eyespot size plasticity in a comparative framework, we showed that 20E titer plasticity in response to temperature is a pre-adaptation shared by all butterfly species examined, whereas expression of EcR in eyespot centers, and eyespot sensitivity to 20E, are both derived traits found only in a subset of species with eyespots.

A well-known family of butterflies have circular patterns on their wings that look like eyes. These eye-like markings help deflect predators away from the butterfly’s body so they attack the outer edges of their wings. However, in certain seasons, such as the dry season in Africa, the best way for this family to survive is by not drawing any attention to their bodies. Thus, butterflies born during this season shrink the size of their eyespots so they can hide among the dry leaves.

How this family of butterflies are able to change the size of these eye-like spots has only been studied in the species Bicyclus anynana. During development low temperatures, which signify the beginning of the dry season, reduce the amount of a hormone called 20E circulating in the blood of this species. This changes the behavior of hormone-sensitive cells in the eyespots making them smaller in size. But it remains unclear how B. anynana evolved this remarkable tactic and whether its relatives have similar abilities.

Now, Bhardwaj et al. show that B. anynana is the only one of its relatives that can amend the size of its eyespots in response to temperature changes. In the experiments, 13 different species of butterflies, mostly from the family that has eyespots, were developed under two different temperatures. Low temperatures caused 20E hormone levels to decrease in all 13 species. However, most of these species did not develop smaller eyespots in response to this temperature change. This includes species that are known to have larger and smaller eyespots depending on the season.

Like B. anynana, four of the species studied have receptors for the 20E hormone at the center of their eyespots. However, changing 20E hormone levels in these species did not reduce eyespot size.

These results show that although temperature changes alter hormone levels in a number of species, only B. anynana have taken advantage of this mechanism to regulate eyespot size. In addition, Bhardwaj et al. found that this unique mechanism evolved from several genetic changes over millions of years. Other species likely use other environmental cues to trigger seasonal changes in the size of their eyespots.

Dissection of Larval and Pupal Wings of Bicyclus anynana Butterflies

The colorful wings of butterflies are emerging as model systems for evolutionary and developmental studies. Some of these studies focus on localizing gene transcripts and proteins in wings at the larval and pupal stages using techniques such as immunostaining and in situ hybridization. Other studies quantify mRNA expression levels or identify regions of open chromatin that are bound by proteins at different stages of wing development. All these techniques require dissection of the wings from the animal but a detailed video protocol describing this procedure has not been available until now. Here, we present a written and accompanying video protocol where we describe the tools and the method we use to remove the larval and pupal wings of the African Squinting Bush Brown butterfly Bicyclus anynana. This protocol should be easy to adapt to other species.

Molecular mechanisms underlying simplification of venation patterns in holometabolous insects

How mechanisms of pattern formation evolve has remained a central research theme in the field of evolutionary and developmental biology. The mechanism of wing vein differentiation in Drosophila is a classic text-book example of pattern formation using a system of positional-information, yet very little is known about how species with a different number of veins pattern their wings, and how insect venation patterns evolved. Here, we examine the expression pattern of genes previously implicated in vein differentiation in Drosophila in two butterfly species with more complex venation Bicyclus anynana and Pieris canidia. We also test the function of some of these genes in B. anynana. We identify both conserved as well as new domains of decapentaplegic, engrailed, invected, spalt, optix, wingless, armadillo, blistered, and rhomboid gene expression in butterflies, and propose how the simplified venation in Drosophila might have evolved via loss of decapentaplegic, spalt and optix gene expression domains, silencing of vein inducing programs at Spalt-expression boundaries, and changes in gene expression of vein maintenance genes.

Conservation and flexibility in the gene regulatory landscape of heliconiine butterfly wings

BACKGROUND

Many traits evolve by cis-regulatory modification, by which changes to noncoding sequences affect the binding affinity for available transcription factors and thus modify the expression profile of genes. Multiple examples of cis-regulatory evolution have been described at pattern switch genes responsible for butterfly wing pattern polymorphism, including in the diverse neotropical genus Heliconius, but the identities of the factors that can regulate these switch genes have not been identified.

RESULTS

We investigated the spatial transcriptomic landscape across the wings of three closely related butterfly species, two of which have a convergently evolved co-mimetic pattern and the other having a divergent pattern. We identified candidate factors for regulating the expression of wing patterning genes, including transcription factors with a conserved expression profile in all three species, and others, including both transcription factors and Wnt pathway genes, with markedly different profiles in each of the three species. We verified the conserved expression profile of the transcription factor homothorax by immunofluorescence and showed that its expression profile strongly correlates with that of the selector gene optix in butterflies with the Amazonian forewing pattern element 'dennis.'

CONCLUSION

Here we show that, in addition to factors with conserved expression profiles like homothorax, there are also a variety of transcription factors and signaling pathway components that appear to vary in their expression profiles between closely related butterfly species, highlighting the importance of genome-wide regulatory evolution between species.

Molecular logic behind the three-way stochastic choices that expand butterfly colour vision

Butterflies rely extensively on colour vision to adapt to the natural world. Most species express a broad range of colour-sensitive Rhodopsin proteins in three types of ommatidia (unit eyes), which are distributed stochastically across the retina. The retinas of Drosophila melanogaster use just two main types, in which fate is controlled by the binary stochastic decision to express the transcription factor Spineless in R7 photoreceptors. We investigated how butterflies instead generate three stochastically distributed ommatidial types, resulting in a more diverse retinal mosaic that provides the basis for additional colour comparisons and an expanded range of colour vision. We show that the Japanese yellow swallowtail (Papilio xuthus, Papilionidae) and the painted lady (Vanessa cardui, Nymphalidae) butterflies have a second R7-like photoreceptor in each ommatidium. Independent stochastic expression of Spineless in each R7-like cell results in expression of a blue-sensitive (Spineless(ON)) or an ultraviolet (UV)-sensitive (Spineless(OFF)) Rhodopsin. In P. xuthus these choices of blue/blue, blue/UV or UV/UV sensitivity in the two R7 cells are coordinated with expression of additional Rhodopsin proteins in the remaining photoreceptors, and together define the three types of ommatidia. Knocking out spineless using CRISPR/Cas9 (refs 5, 6) leads to the loss of the blue-sensitive fate in R7-like cells and transforms retinas into homogeneous fields of UV/UV-type ommatidia, with corresponding changes in other coordinated features of ommatidial type. Hence, the three possible outcomes of Spineless expression define the three ommatidial types in butterflies. This developmental strategy allowed the deployment of an additional red-sensitive Rhodopsin in P. xuthus, allowing for the evolution of expanded colour vision with a greater variety of receptors. This surprisingly simple mechanism that makes use of two binary stochastic decisions coupled with local coordination may prove to be a general means of generating an increased diversity of developmental outcomes.

Transcriptome analysis of the painted lady butterfly, Vanessa cardui during wing color pattern development

BACKGROUND

Butterfly wing color patterns are an important model system for understanding the evolution and development of morphological diversity and animal pigmentation. Wing color patterns develop from a complex network composed of highly conserved patterning genes and pigmentation pathways. Patterning genes are involved in regulating pigment synthesis however the temporal expression dynamics of these interacting networks is poorly understood. Here, we employ next generation sequencing to examine expression patterns of the gene network underlying wing development in the nymphalid butterfly, Vanessa cardui.

RESULTS

We identified 9, 376 differentially expressed transcripts during wing color pattern development, including genes involved in patterning, pigmentation and gene regulation. Differential expression of these genes was highest at the pre-ommochrome stage compared to early pupal and late melanin stages. Overall, an increasing number of genes were down-regulated during the progression of wing development. We observed dynamic expression patterns of a large number of pigment genes from the ommochrome, melanin and also pteridine pathways, including contrasting patterns of expression for paralogs of the yellow gene family. Surprisingly, many patterning genes previously associated with butterfly pattern elements were not significantly up-regulated at any time during pupation, although many other transcription factors were differentially expressed. Several genes involved in Notch signaling were significantly up-regulated during the pre-ommochrome stage including slow border cells, bunched and pebbles; the function of these genes in the development of butterfly wings is currently unknown. Many genes involved in ecdysone signaling were also significantly up-regulated during early pupal and late melanin stages and exhibited opposing patterns of expression relative to the ecdysone receptor. Finally, a comparison across four butterfly transcriptomes revealed 28 transcripts common to all four species that have no known homologs in other metazoans.

CONCLUSIONS

This study provides a comprehensive list of differentially expressed transcripts during wing development, revealing potential candidate genes that may be involved in regulating butterfly wing patterns. Some differentially expressed genes have no known homologs possibly representing genes unique to butterflies. Results from this study also indicate that development of nymphalid wing patterns may arise not only from melanin and ommochrome pigments but also the pteridine pigment pathway.

Differential Expression of Ecdysone Receptor Leads to Variation in Phenotypic Plasticity across Serial Homologs

Bodies are often made of repeated units, or serial homologs, that develop using the same core gene regulatory network. Local inputs and modifications to this network allow serial homologs to evolve different morphologies, but currently we do not understand which modifications allow these repeated traits to evolve different levels of phenotypic plasticity. Here we describe variation in phenotypic plasticity across serial homologous eyespots of the butterfly Bicyclus anynana, hypothesized to be under selection for similar or different functions in the wet and dry seasonal forms. Specifically, we document the presence of eyespot size and scale brightness plasticity in hindwing eyespots hypothesized to vary in function across seasons, and reduced size plasticity and absence of brightness plasticity in forewing eyespots hypothesized to have the same function across seasons. By exploring the molecular and physiological causes of this variation in plasticity across fore and hindwing serial homologs we discover that: 1) temperature experienced during the wandering stages of larval development alters titers of an ecdysteroid hormone, 20-hydroxyecdysone (20E), in the hemolymph of wet and dry seasonal forms at that stage; 2) the 20E receptor (EcR) is differentially expressed in the forewing and hindwing eyespot centers of both seasonal forms during this critical developmental stage; and 3) manipulations of EcR signaling disproportionately affected hindwing eyespots relative to forewing eyespots. We propose that differential EcR expression across forewing and hindwing eyespots at a critical stage of development explains the variation in levels of phenotypic plasticity across these serial homologues. This finding provides a novel signaling pathway, 20E, and a novel molecular candidate, EcR, for the regulation of levels of phenotypic plasticity across body parts or serial homologs.

One of the most exquisite types of organismal adaptations in nature occurs when organisms are able to change the way they develop in anticipation of the different selective environments they will experience as adults. This leads to variation in adult morphologies that are adaptive. Environmental cues experienced during development often lead to variation in hormonal titers that can have profound effect on the way genes are regulated and on the adult morphology. Here we examine the hormonal and molecular mechanisms that allow specific traits that are repeated in an organism (butterfly eyespots) to either be sensitive to environmental cues–and develop different morphologies—or insensitive to these cues and develop the same morphology. We discover that a specific gene, a hormone receptor, that is expressed in the sensitive eyespots but absent in the insensitive eyespots, is responsible for regulating the level of sensitivity of each of the eyespots to an environmental cue. We identify a molecule that is regulating levels of environmental sensitivity, or phenotypic plasticity, across repeated traits in an organism.

RNAi screening of developmental toolkit genes: a search for novel wing genes in the red flour beetle, Tribolium castaneum

The amazing array of diversity among insect wings offers a powerful opportunity to study the mechanisms guiding morphological evolution. Studies in Drosophila (the fruit fly) have identified dozens of genes important for wing development. These genes are often called candidate genes, serving as an ideal starting point to study wing development in other insects. However, we also need to explore beyond the candidate genes to gain a more comprehensive view of insect wing evolution. As a first step away from the traditional candidate genes, we utilized Tribolium (the red flour beetle) as a model and assessed the potential involvement of a group of developmental toolkit genes (embryonic patterning genes) in beetle wing development. We hypothesized that the highly pleiotropic nature of these developmental genes would increase the likelihood of finding novel wing genes in Tribolium. Through the RNA interference screening, we found that Tc-cactus has a less characterized (but potentially evolutionarily conserved) role in wing development. We also found that the odd-skipped family genes are essential for the formation of the thoracic pleural plates, including the recently discovered wing serial homologs in Tribolium. In addition, we obtained several novel insights into the function of these developmental genes, such as the involvement of mille-pattes and Tc-odd-paired in metamorphosis. Despite these findings, no gene we examined was found to have novel wing-related roles unique in Tribolium. These results suggest a relatively conserved nature of developmental toolkit genes and highlight the limited degree to which these genes are co-opted during insect wing evolution.

Nymphalid eyespot serial homologues originate as a few individualized modules

Serial homologues are repeated traits that share similar development but occur in different parts of the body. Variation in number of repeats accounts for substantial diversity in animal form and considerable work has focused on identifying the factors accounting for this variation. Little is known, however, about how serial homologues originally become repeated, or about the relative timing of repeat individuation relative to repeat origin. Here, we show that the serially repeated eyespots on nymphalid butterfly wings most likely arose as a small cluster of units on the ventral hindwing that were later co-opted to the dorsal and anterior wing surfaces. Based on comparative analyses of over 400 species, we found support for a model of eyespot origin followed by redeployment, rather than by the conventional model, where eyespots arose as a complete row of undifferentiated units that later gained individuation. In addition, eyespots most likely evolved from simpler pattern elements, single-coloured spots, which were already individuated among different wing sectors. Finally, the late appearance of eyespots on the dorsal, hidden wing surface further suggests that these novel complex traits originally evolved for one function (thwarting predator attacks) and acquired a second function (sexual signalling) when moved to a different body location. This broad comparative analysis illustrates how serial homologues may initially evolve as a few units serving a particular function and subsequently become repeated in novel body locations with new functions.

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.

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

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

Dynamics of F-actin prefigure the structure of butterfly wing scales

The wings of butterflies and moths consist of dorsal and ventral epidermal surfaces that give rise to overlapping layers of scales and hairs (Lepidoptera, "scale wing"). Wing scales (average length ~200 µm) are homologous to insect bristles (macrochaetes), and their colors create the patterns that characterize lepidopteran wings. The topology and surface sculpture of wing scales vary widely, and this architectural complexity arises from variations in the developmental program of the individual scale cells of the wing epithelium. One of the more striking features of lepidopteran wing scales are the longitudinal ridges that run the length of the mature (dead) cell, gathering the cuticularized scale cell surface into pleats on the sides of each scale. While also present around the periphery of other insect bristles and hairs, longitudinal ridges in lepidopteran wing scales gain new significance for their creation of iridescent color through microribs and lamellae. Here we show the dynamics of the highly organized F-actin filaments during scale cell development, and present experimental manipulations of actin polymerization that reveal the essential role of this cytoskeletal component in wing scale elongation and the positioning of longitudinal ribs.