A single origin for nymphalid butterfly eyespots followed by widespread loss of associated gene expression

Nymphalid eyespots are co-opted to novel wing locations following a similar pattern in independent lineages

BACKGROUND

Variation in the number of repeated traits, or serial homologs, has contributed greatly to animal body plan diversity. Eyespot color patterns of nymphalid butterflies, like arthropod and vertebrate limbs, are an example of serial homologs. These eyespot color patterns originated in a small number of wing sectors on the ventral hindwing surface and later appeared in novel wing sectors, novel wings, and novel wing surfaces. However, the details of how eyespots were co-opted to these novel wing locations are currently unknown.

RESULTS

We used a large data matrix of eyespot/presence absence data, previously assembled from photographs of contemporary species, to perform a phylogenetic investigation of eyespot origins in nine independent nymphalid lineages. To determine how the eyespot gene regulatory network acquired novel positional information, we used phylogenetic correlation analyses to test for non-independence in the origination of eyespots. We found consistent patterns of eyespot gene network redeployment in the nine lineages, where eyespots first redeployed from the ventral hindwing to the ventral forewing, then to new sectors within the ventral wing surface, and finally to the dorsal wing surface. Eyespots that appeared in novel wing sectors modified the positional information of their serial homolog ancestors in one of two ways: by changing the wing or surface identity while retaining sector identity, or by changing the sector identity while retaining wing and surface identity.

CONCLUSIONS

Eyespot redeployment to novel sectors, wings, and surfaces happened multiple times in different nymphalid subfamilies following a similar pattern. This indicates that parallel mutations altering expression of the eyespot gene regulatory network led to its co-option to novel wing locations over time.

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.

A Burst of miRNA Innovation in the Early Evolution of Butterflies and Moths

MicroRNAs (miRNAs) are involved in posttranscriptional regulation of gene expression. Because several miRNAs are known to affect the stability or translation of developmental regulatory genes, the origin of novel miRNAs may have contributed to the evolution of developmental processes and morphology. Lepidoptera (butterflies and moths) is a species-rich clade with a well-established phylogeny and abundant genomic resources, thereby representing an ideal system in which to study miRNA evolution. We sequenced small RNA libraries from developmental stages of two divergent lepidopterans, Cameraria ohridella (Horse chestnut Leafminer) and Pararge aegeria (Speckled Wood butterfly), discovering 90 and 81 conserved miRNAs, respectively, and many species-specific miRNA sequences. Mapping miRNAs onto the lepidopteran phylogeny reveals rapid miRNA turnover and an episode of miRNA fixation early in lepidopteran evolution, implying that miRNA acquisition accompanied the early radiation of the Lepidoptera. One lepidopteran-specific miRNA gene, miR-2768, is located within an intron of the homeobox gene invected, involved in insect segmental and wing patterning. We identified cubitus interruptus (ci) as a likely direct target of miR-2768, and validated this suppression using a luciferase assay system. We propose a model by which miR-2768 modulates expression of ci in the segmentation pathway and in patterning of lepidopteran wing primordia.

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.

Metamorphic labral axis patterning in the beetle Tribolium castaneum requires multiple upstream, but few downstream, genes in the appendage patterning network

The arthropod labrum is an anterior appendage-like structure that forms the dorsal side of the preoral cavity. Conflicting interpretations of fossil, nervous system, and developmental data have led to a proliferation of scenarios for labral evolution. The best supported hypothesis is that the labrum is a novel structure that shares development with appendages as a result of co-option. Here, we use RNA interference in the red flour beetle Tribolium castaneum to compare metamorphic patterning of the labrum to previously published data on ventral appendage patterning. As expected under the co-option hypothesis, depletion of several genes resulted in similar defects in the labrum and ventral appendages. These include proximal deletions and proximal-to-distal transformations resulting from depletion of the leg gap genes homothorax and extradenticle, large-scale deletions resulting from depletion of the leg gap gene Distal-less, and smaller distal deletions resulting from knockdown of the EGF ligand Keren. However, depletion of dachshund and many of the genes that function downstream of the leg gap genes in the ventral appendages had either subtle or no effects on labral axis patterning. This pattern of partial similarity suggests that upstream genes act through different downstream targets in the labrum. We also discovered that many appendage axis patterning genes have roles in patterning the epipharyngeal sensillum array, suggesting that they have become integrated into a novel regulatory network. These genes include Notch, Delta, and decapentaplegic, and the transcription factors abrupt, bric à brac, homothorax, extradenticle and the paralogs apterous a and apterous b.

Gradual and contingent evolutionary emergence of leaf mimicry in butterfly wing patterns

Background

Special resemblance of animals to natural objects such as leaves provides a representative example of evolutionary adaptation. The existence of such sophisticated features challenges our understanding of how complex adaptive phenotypes evolved. Leaf mimicry typically consists of several pattern elements, the spatial arrangement of which generates the leaf venation-like appearance. However, the process by which leaf patterns evolved remains unclear.

Results

In this study we show the evolutionary origin and process for the leaf pattern in Kallima (Nymphalidae) butterflies. Using comparative morphological analyses, we reveal that the wing patterns of Kallima and 45 closely related species share the same ground plan, suggesting that the pattern elements of leaf mimicry have been inherited across species with lineage-specific changes of their character states. On the basis of these analyses, phylogenetic comparative methods estimated past states of the pattern elements and enabled reconstruction of the wing patterns of the most recent common ancestor. This analysis shows that the leaf pattern has evolved through several intermediate patterns. Further, we use Bayesian statistical methods to estimate the temporal order of character-state changes in the pattern elements by which leaf mimesis evolved, and show that the pattern elements changed their spatial arrangement (e.g., from a curved line to a straight line) in a stepwise manner and finally establish a close resemblance to a leaf venation-like appearance.

Conclusions

Our study provides the first evidence for stepwise and contingent evolution of leaf mimicry.　 Leaf mimicry patterns evolved in a gradual, rather than a sudden, manner from a non-mimetic ancestor. Through a lineage of Kallima butterflies, the leaf patterns evolutionarily originated through temporal accumulation of orchestrated changes in multiple pattern elements.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0229-5) contains supplementary material, which is available to authorized users.

Ancient homology underlies adaptive mimetic diversity across butterflies

Convergent evolution provides a rare, natural experiment with which to test the predictability of adaptation at the molecular level. Little is known about the molecular basis of convergence over macro-evolutionary timescales. Here we use a combination of positional cloning, population genomic resequencing, association mapping and developmental data to demonstrate that positionally orthologous nucleotide variants in the upstream region of the same gene, WntA, are responsible for parallel mimetic variation in two butterfly lineages that diverged >65 million years ago. Furthermore, characterization of spatial patterns of WntA expression during development suggests that alternative regulatory mechanisms underlie wing pattern variation in each system. Taken together, our results reveal a strikingly predictable molecular basis for phenotypic convergence over deep evolutionary time.

Little is known about the genetic basis of convergent evolution in deeply diverged species. Here, the authors show that variation in the WntA gene is associated with parallel wing pattern variation in two butterflies that diverged more than 65 million years ago.

Wnt signaling underlies evolution and development of the butterfly wing pattern symmetry systems

Most butterfly wing patterns are proposed to be derived from a set of conserved pattern elements known as symmetry systems. Symmetry systems are so-named because they are often associated with parallel color stripes mirrored around linear organizing centers that run between the anterior and posterior wing margins. Even though the symmetry systems are the most prominent and diverse wing pattern elements, their study has been confounded by a lack of knowledge regarding the molecular basis of their development, as well as the difficulty of drawing pattern homologies across species with highly derived wing patterns. Here we present the first molecular characterization of symmetry system development by showing that WntA expression is consistently associated with the major basal, discal, central, and external symmetry system patterns of nymphalid butterflies. Pharmacological manipulations of signaling gradients using heparin and dextran sulfate showed that pattern organizing centers correspond precisely with WntA, wingless, Wnt6, and Wnt10 expression patterns, thus suggesting a role for Wnt signaling in color pattern induction. Importantly, this model is supported by recent genetic and population genomic work identifying WntA as the causative locus underlying wing pattern variation within several butterfly species. By comparing the expression of WntA between nymphalid butterflies representing a range of prototypical symmetry systems, slightly deviated symmetry systems, and highly derived wing patterns, we were able to infer symmetry system homologies in several challenging cases. Our work illustrates how highly divergent morphologies can be derived from modifications to a common ground plan across both micro- and macro-evolutionary time scales.

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.

Phylogeny, systematics and biogeography of the genus panolis (lepidoptera: noctuidae) based on morphological and molecular evidence

The genus Panolis is a small group of noctuid moths with six recognized species distributed from Europe to East Asia, and best known for containing the widespread Palearctic pest species P. flammea, the pine beauty moth. However, a reliable classification and robust phylogenetic framework for this group of potentially economic importance are currently lacking. Here, we use morphological and molecular data (mitochondrial genes cytochrome c oxidase subunit I and 16S ribosomal RNA, nuclear gene elongation factor-1 alpha) to reconstruct the phylogeny of this genus, with a comprehensive systematic revision of all recognized species and a new one, P. ningshan sp. nov. The analysis results of maximum parsimony, maximum likelihood and Bayesian inferring methods for the combined morphological and molecular data sets are highly congruent, resulting in a robust phylogeny and identification of two clear species groups, i.e., the P. flammea species group and the P. exquisita species group. We also estimate the divergence times of Panolis moths using two conventional mutation rates for the arthropod mitochondrial COI gene with a comparison of two molecular clock models, as well as reconstruct their ancestral areas. Our results suggest that 1) Panolis is a young clade, originating from the Oriental region in China in the Late Miocene (6-10Mya), with an ancestral species in the P. flammea group extending northward to the Palearctic region some 3-6 Mya; 2) there is a clear possibility for a representative of the Palearctic clade to become established as an invasive species in the Nearctic taiga.

A Survey of Eyespot Sexual Dimorphism across Nymphalid Butterflies

Differences between sexes of the same species are widespread and are variable in nature. While it is often assumed that males are more ornamented than females, in the nymphalid butterfly genus Bicyclus, females have, on average, more eyespot wing color patterns than males. Here we extend these studies by surveying eyespot pattern sexual dimorphism across the Nymphalidae family of butterflies. Eyespot presence or absence was scored from a total of 38 wing compartments for two males and two females of each of 450 nymphalid species belonging to 399 different genera. Differences in eyespot number between sexes of each species were tallied for each wing surface (e.g., dorsal and ventral) of forewings and hindwings. In roughly 44% of the species with eyespots, females had more eyespots than males, in 34%, males had more eyespots than females, and, in the remaining 22% of the species, there was monomorphism in eyespot number. Dorsal and forewing surfaces were less patterned, but proportionally more dimorphic, than ventral and hindwing surfaces, respectively. In addition, wing compartments that frequently displayed eyespots were among the least sexually dimorphic. This survey suggests that dimorphism arises predominantly in "hidden" or "private" surfaces of a butterfly's wing, as previously demonstrated for the genus Bicyclus.

Temporal gene expression variation associated with eyespot size plasticity in Bicyclus anynana

Seasonal polyphenism demonstrates an organism's ability to respond to predictable environmental variation with alternative phenotypes, each presumably better suited to its respective environment. However, the molecular mechanisms linking environmental variation to alternative phenotypes via shifts in development remain relatively unknown. Here we investigate temporal gene expression variation in the seasonally polyphenic butterfly Bicyclus anynana. This species shows drastic changes in eyespot size depending on the temperature experienced during larval development. The wet season form (larvae reared over 24°C) has large ventral wing eyespots while the dry season form (larvae reared under 19°C) has much smaller eyespots. We compared the expression of three proteins, Notch, Engrailed, and Distal-less, in the future eyespot centers of the two forms to determine if eyespot size variation is associated with heterochronic shifts in the onset of their expression. For two of these proteins, Notch and Engrailed, expression in eyespot centers occurred earlier in dry season than in wet season larvae, while Distal-less showed no temporal difference between the two forms. These results suggest that differences between dry and wet season adult wings could be due to a delay in the onset of expression of these eyespot-associated genes. Early in eyespot development, Notch and Engrailed may be functioning as repressors rather than activators of the eyespot gene network. Alternatively, temporal variation in the onset of early expressed genes between forms may have no functional consequences to eyespot size regulation and may indicate the presence of an 'hourglass' model of development in butterfly eyespots.

Baculovirus-mediated gene transfer in butterfly wings in vivo: an efficient expression system with an anti-gp64 antibody

BACKGROUND

Candidate genes for color pattern formation in butterfly wings have been known based on gene expression patterns since the 1990s, but their functions remain elusive due to a lack of a functional assay. Several methods of transferring and expressing a foreign gene in butterfly wings have been reported, but they have suffered from low success rates or low expression levels. Here, we developed a simple, practical method to efficiently deliver and express a foreign gene using baculovirus-mediated gene transfer in butterfly wings in vivo.

RESULTS

A recombinant baculovirus containing a gene for green fluorescent protein (GFP) was injected into pupae of the blue pansy butterfly Junonia orithya (Nymphalidae). GFP fluorescence was detected in the pupal wings and other body parts of the injected individuals three to five days post-injection at various degrees of fluorescence. We obtained a high GFP expression rate at relatively high virus titers, but it was associated with pupal death before color pattern formation in wings. To reduce the high mortality rate caused by the baculovirus treatment, we administered an anti-gp64 antibody, which was raised against baculovirus coat protein gp64, to infected pupae after the baculovirus injection. This treatment greatly reduced the mortality rate of the infected pupae. GFP fluorescence was observed in pupal and adult wings and other body parts of the antibody-treated individuals at various degrees of fluorescence. Importantly, we obtained completely developed wings with a normal color pattern, in which fluorescent signals originated directly from scales or the basal membrane after the removal of scales. GFP fluorescence in wing tissues spatially coincided with anti-GFP antibody staining, confirming that the fluorescent signals originated from the expressed GFP molecules.

CONCLUSIONS

Our baculovirus-mediated gene transfer system with an anti-gp64 antibody is reasonably efficient, and it can be an invaluable tool to transfer, express, and functionally examine foreign genes in butterfly wings and also in other non-model insect systems.

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

Background

Little is currently known about wing pattern development in the butterfly family Pieridae, which consists mostly of black melanized elements on white or yellow/orange backgrounds. A single transcription factor, Spalt (Sal), has been previously associated with the development of some pattern elements in Pieris rapae, but it is unclear to what extent Sal is associated with patterns in other pierid species.

Results

We use immunohistochemistry targeting Sal proteins across several pierids and show that Sal is associated with dense patches of melanization across species but is not associated with vein-melanization or diffuse melanization on the wing. In addition, Sal is expressed along cross-veins and wing compartment midlines that do not develop melanization. Male and female P. rapae spots are sexually dimorphic in size and this dimorphism is also present in the domains of Sal expression. Finally, by disrupting cells positioned in the center of the anterior black spots of P. rapae, before and during the time of Sal expression, spot size was reduced.

Conclusions

Our results suggest, but do not conclusively show, that pierid spots may develop in a manner similar to that of nymphalid eyespots, that is, containing a group of signaling cells at the center of the pattern responsible for the differentiation of the complete spot, and that spots and eyespots share at least one signal-response gene in common, the transcription factor Sal. We propose that focal differentiation and focal signaling mechanisms evolved prior to the split of the nymphalid and pierid lineages.

Differential involvement of Hedgehog signaling in butterfly wing and eyespot development

Butterfly eyespots may have evolved from the recruitment of pre-existent gene circuits or regulatory networks into novel locations on the wing. Gene expression data suggests one such circuit, the Hedgehog (Hh) signaling pathway and its target gene engrailed (en), was recruited from a role in patterning the anterior-posterior insect wing axis to a role patterning butterfly eyespots. However, while Junonia coenia expresses hh and en both in the posterior compartment of the wing and in eyespot centers, Bicyclus anynana lacks hh eyespot-specific expression. This suggests that Hh signaling may not be functioning in eyespot development in either species or that it functions in J. coenia but not in B. anynana. In order to test these hypotheses, we performed functional tests of Hh signaling in these species. We investigated the effects of Hh protein sequestration during the larval stage on en expression levels, and on wing size and eyespot size in adults. Hh sequestration led to significantly reduced en expression and to significantly smaller wings and eyespots in both species. But while eyespot size in B. anynana was reduced proportionately to wing size, in J. coenia, eyespots were reduced disproportionately, indicating an independent role of Hh signaling in eyespot development in J. coenia. We conclude that while Hh signaling retains a conserved role in promoting wing growth across nymphalid butterflies, it plays an additional role in eyespot development in some, but not all, lineages of nymphalid butterflies. We discuss our findings in the context of alternative evolutionary scenarios that led to the differential expression of hh and other Hh pathway signaling members across nymphalid species.