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

Transcriptome analysis and temporal expression patterns of wing development-related genes in Lymantria dispar (Lepidoptera: Erebidae)

Spongy moth, Lymantria dispar Linnaeus (Lepidoptera: Erebidae), stands as a pervasive international threat, marked by its designation as one of the "world's 100 worst invasive species" by IUCN, owing to its voracious leaf-eating habits encompassing over 500 plant species. Its strong flight ability facilitates its spread and invasion. The present study aims to uncover differential gene expression, utilizing the Illumina Novaseq6000 sequencing platform for comprehensive transcriptome sequencing and bioinformatic analysis of total RNA extracted from larvae and pupae. Results revealed pivotal processes of protein functional structure conformation, transport, and signal transduction in functional gene annotation during the 2 developmental stages of spongy moth. 18 functional genes, namely, Distal-less (Dll), Wingless (Wg), Decapentaplegic (Dpp), Hedgehog (Hh), Cubitus interruptus (Ci), Patched (Ptc), Apterous (Ap), Serrate (Ser), Fringe (Fng), Achaete (Ac), Engrailed (En), Vestigial (Vg), Scute (Sc), Invected (Inv), Scalloped (Sd), Ultrabithorax (Ubx), Serum Response Factor (SRF), and Spalt-major, associated with wing development were identified, and their expression levels were meticulously assessed through real-time quantitative PCR (RT-qPCR) in 1st-6th instar larvae and male and female pupae wing discs. The results showed that 18 genes exhibited expression. Furthermore, the relative expression values of wing development-related genes were significantly higher in the pupae stage than in the larval stage. The relative expression values of male and female pupae were also significantly different. The RT-qPCR results were in general agreement with the results of transcriptome analysis. This study establishes a foundational understanding of the developmental mechanisms governing the formation of spongy moth wings.

The transition to flying insects: lessons from evo-devo and fossils

Insects are the only arthropod group to achieve powered flight, which facilitated their explosive radiation on land. It remains a significant challenge to understand the evolutionary transition from nonflying (apterygote) to flying (pterygote) insects due to the large gap in the fossil record. Under such a situation, ontogenic information has historically been used to compensate for fossil evidence. Recent evo-devo studies support and refine a paleontology-based classical hypothesis that an ancestral exite incorporated into the body wall contributed to the origin of insect wings. The modern hypothesis locates an ancestral precoxa leg segment with an exite within the hexapod lateral tergum, reframing the long-standing debate on the insect wing origin. A current focus is on the contributions of the incorporated exite homolog and surrounding tissues, such as the pleuron and the medial bona fide tergum, to wing evolution. In parallel, recent analyses of Paleozoic fossils have confirmed thoracic and abdominal lateral body outgrowths as transitional wing precursors and suggested their possible role as respiratory organs in aquatic or semiaquatic environments. These recent studies have revised our understanding of the transition to flying insects. This review highlights recent progress in both evo-devo and paleontology, and discusses future challenges, including the evolution of metamorphic development.

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.

Genomic screening of fish-specific genes in gnathostomes and their functions in fin development

In this study, we comprehensively searched for fish-specific genes in gnathostomes that contribute to development of the fin, a fish-specific trait. Many previous reports suggested that animal group-specific genes are often important for group-specific traits. Clarifying the roles of fish-specific genes in fin development of gnathostomes, for example, can help elucidate the mechanisms underlying the formation of this trait. We first identified 91 fish-specific genes in gnathostomes by comparing the gene repertoire in 16 fish and 35 tetrapod species. RNA-seq analysis narrowed down the 91 candidates to 33 genes that were expressed in the developing pectoral fin. We analyzed the functions of approximately half of the candidate genes by loss-of-function analysis in zebrafish. We found that some of the fish-specific and fin development-related genes, including fgf24 and and1/and2, play roles in fin development. In particular, the newly identified fish-specific gene qkia is expressed in the developing fin muscle and contributes to muscle morphogenesis in the pectoral fin as well as body trunk. These results indicate that the strategy of identifying animal group-specific genes is functional and useful. The methods applied here could be used in future studies to identify trait-associated genes in other animal groups.

A genetic cause of male mate preference

A gene for mate preference has been shared between hybridizing butterfly species.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

A colourful duplication

A genetic duplication event during evolution allowed male wood tiger moths to have either yellow or white patterns on their wings.

Frizzled2 receives WntA signaling during butterfly wing pattern formation

Butterfly color patterns provide visible and biodiverse phenotypic readouts of the patterning processes. Although the secreted ligand WntA has been shown to instruct the color pattern formation in butterflies, its mode of reception remains elusive. Butterfly genomes encode four homologs of the Frizzled-family of Wnt receptors. Here, we show that CRISPR mosaic knockouts of frizzled2 (fz2) phenocopy the color pattern effects of WntA loss of function in multiple nymphalids. Whereas WntA mosaic clones result in intermediate patterns of reduced size, fz2 clones are cell-autonomous, consistent with a morphogen function. Shifts in expression of WntA and fz2 in WntA crispant pupae show that they are under positive and negative feedback, respectively. Fz1 is required for Wnt-independent planar cell polarity in the wing epithelium. Fz3 and Fz4 show phenotypes consistent with Wnt competitive-antagonist functions in vein formation (Fz3 and Fz4), wing margin specification (Fz3), and color patterning in the Discalis and Marginal Band Systems (Fz4). Overall, these data show that the WntA/Frizzled2 morphogen-receptor pair forms a signaling axis that instructs butterfly color patterning and shed light on the functional diversity of insect Frizzled receptors.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae)

BACKGROUND

Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics.

RESULTS

The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function.

CONCLUSIONS

Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

Clade-specific genes and the evolutionary origin of novelty; new tools in the toolkit

Clade-specific (a.k.a. lineage-specific) genes are very common and found at all taxonomic levels and in all clades examined. They can arise by duplication of previously existing genes, which can involve partial truncations or combinations with other protein domains or regulatory sequences. They can also evolve de novo from non-coding sequences, leading to potentially truly novel protein domains. Finally, since clade-specific genes are generally defined by lack of sequence homology with other proteins, they can also arise by sequence evolution that is rapid enough that previous sequence homology can no longer be detected. In such cases, where the rapid evolution is followed by constraint, we consider them to be ontologically non-novel but likely novel at a functional level. In general, clade-specific genes have received less attention from biologists but there are increasing numbers of fascinating examples of their roles in important traits. Here we review some selected recent examples, and argue that attention to clade-specific genes is an important corrective to the focus on the conserved developmental regulatory toolkit that has been the habit of evo-devo as a field. Finally, we discuss questions that arise about the evolution of clade-specific genes, and how these might be addressed by future studies. We highlight the hypothesis that clade-specific genes are more likely to be involved in synapomorphies that arose in the stem group where they appeared, compared to other genes.

Spatial and temporal regulation of Wnt signaling pathway members in the development of butterfly wing patterns

Wnt signaling members are involved in the differentiation of cells associated with eyespot and band color patterns on the wings of butterflies, but the identity and spatio-temporal regulation of specific Wnt pathway members remains unclear. Here, we explore the localization and function of Armadillo/β-catenin dependent (canonical) and Armadillo/β-catenin independent (noncanonical) Wnt signaling in eyespot and band development in Bicyclus anynana by localizing Armadillo (Arm), the expression of all eight Wnt ligand and four frizzled receptor transcripts present in the genome of this species and testing the function of some of the ligands and receptors using CRISPR-Cas9. We show that distinct Wnt signaling pathways are essential for eyespot and band patterning in butterflies and are likely interacting to control their active domains.

Divergent expression of aristaless1 and aristaless2 during embryonic appendage and pupal wing development in butterflies

BACKGROUND

Gene duplication events are critical for the evolution of new gene functions. Aristaless is a major regulator of distinct developmental processes. It is most known for its role during appendage development across animals. However, more recently other distinct biological functions have been described for this gene and its duplicates. Butterflies and moths have two copies of aristaless, aristaless1 (al1) and aristaless2 (al2), as a result of a gene duplication event. Previous work in Heliconius has shown that both copies appear to have novel functions related to wing color patterning. Here we expand our knowledge of the expression profiles associated with both ancestral and novel functions of Al1 across embryogenesis and wing pigmentation. Furthermore, we characterize Al2 expression, providing a comparative framework between gene copies within the same species, allowing us to understand the origin of new functions following gene duplication.

RESULTS

Our work shows that the expression of both Al1 and Al2 is associated with the ancestral function of sensory appendage (leg, mouth, spines, and eyes) development in embryos. Interestingly, Al1 exhibits higher expression earlier in embryogenesis while the highest levels of Al2 expression are shifted to later stages of embryonic development. Furthermore, Al1 localization appears extranuclear while Al2 co-localizes tightly with nuclei earlier, and then also expands outside the nucleus later in development. Cellular expression of Al1 and Al2 in pupal wings is broadly consistent with patterns observed during embryogenesis. We also describe, for the first time, how Al1 localization appears to correlate with zones of anterior/posterior elongation of the body during embryonic growth, showcasing a possible new function related to Aristaless' previously described role in appendage extension.

CONCLUSIONS

Overall, our data suggest that while both gene copies play a role in embryogenesis and wing pigmentation, the duplicates have diverged temporally and mechanistically across those functions. Our study helps clarify principles behind sub-functionalization and gene expression evolution associated with developmental functions following gene duplication events.

aristaless1 has a dual role in appendage formation and wing color specification during butterfly development

BACKGROUND

Highly diverse butterfly wing patterns have emerged as a powerful system for understanding the genetic basis of phenotypic variation. While the genetic basis of this pattern variation is being clarified, the precise developmental pathways linking genotype to phenotype are not well understood. The gene aristaless, which plays a role in appendage patterning and extension, has been duplicated in Lepidoptera. One copy, aristaless1, has been shown to control a white/yellow color switch in the butterfly Heliconius cydno, suggesting a novel function associated with color patterning and pigmentation. Here we investigate the developmental basis of al1 in embryos, larvae, and pupae using new antibodies, CRISPR/Cas9, RNAi, qPCR assays of downstream targets, and pharmacological manipulation of an upstream activator.

RESULTS

We find that Al1 is expressed at the distal tips of developing embryonic appendages consistent with its ancestral role. In developing wings, we observe Al1 accumulation within developing scale cells of white H. cydno during early pupation while yellow scale cells exhibit little Al1 at this time point. Reduced Al1 expression is also associated with yellow scale development in al1 knockouts and knockdowns. We propose that Al1 expression in future white scales might be related to an observed downregulation of the enzyme Cinnabar and other genes that synthesize and transport the yellow pigment, 3-hydroxykynurenine (3-OHK). Finally, we provide evidence that Al1 activation is under the control of Wnt signaling.

CONCLUSIONS

We propose a model in which high levels of Al1 during early pupation, which are mediated by Wnt, are important for melanic pigmentation and specifying white portions of the wing while reduced levels of Al1 during early pupation promote upregulation of proteins needed to move and synthesize 3-OHK, promoting yellow pigmentation. In addition, we discuss how the ancestral role of aristaless in appendage extension may be relevant in understanding the cellular mechanism behind color patterning in the context of the heterochrony hypothesis.

High level of novelty under the hood of convergent evolution

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

Butterfly Wing Translucence Enables Enhanced Visual Signaling

The light reflected by the dorsal side of butterfly wings often functions as a signal for, e.g., mate choice, thermoregulation, and/or predator deterrence, while the ventral wing reflections are generally used for crypsis and camouflage. Here, we propose that transmitted light can also have an important role in visual signaling because, in many butterfly species, the dorsal and ventral wing sides are similarly patterned and locally more or less translucent. Extreme examples are the Japanese yellow swallowtail (Papilio xuthus Linnaeus, 1758) and the Yellow glassy tiger (Parantica aspasia Fabricius, 1787). Their wings exhibit a similar color pattern in reflected and transmitted light, which allows enhanced visual signaling, especially in flight. Contrasting cases in which the coloration and patterning of dorsal and ventral wings strongly differ are the papilionid Papilio nireus Linnaeus, 1758, and the pierid Delias nigrina Fabricius, 1775. The wings observed in reflected or transmitted light then show very different color patterns. Wing translucence thus will strongly affect a butterfly's visual signal.

Stage- and sex-specific transcriptome analyses reveal distinctive sensory gene expression patterns in a butterfly

Background

Animal behavior is largely driven by the information that animals are able to extract and process from their environment. However, the function and organization of sensory systems often change throughout ontogeny, particularly in animals that undergo indirect development. As an initial step toward investigating these ontogenetic changes at the molecular level, we characterized the sensory gene repertoire and examined the expression profiles of genes linked to vision and chemosensation in two life stages of an insect that goes through metamorphosis, the butterfly Bicyclus anynana.

Results

Using RNA-seq, we compared gene expression in the heads of late fifth instar larvae and newly eclosed adults that were reared under identical conditions. Over 50 % of all expressed genes were differentially expressed between the two developmental stages, with 4,036 genes upregulated in larval heads and 4,348 genes upregulated in adult heads. In larvae, upregulated vision-related genes were biased toward those involved with eye development, while phototransduction genes dominated the vision genes that were upregulated in adults. Moreover, the majority of the chemosensory genes we identified in the B. anynana genome were differentially expressed between larvae and adults, several of which share homology with genes linked to pheromone detection, host plant recognition, and foraging in other species of Lepidoptera.

Conclusions

These results revealed promising candidates for furthering our understanding of sensory processing and behavior in the disparate developmental stages of butterflies and other animals that undergo metamorphosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07819-4.

Porcupine/Wntless-dependent trafficking of the conserved WntA ligand in butterflies

Wnt ligands are key signaling molecules in animals, but little is known about the evolutionary dynamics and mode of action of the WntA orthologs, which are not present in the vertebrates or in Drosophila. Here we show that the WntA subfamily evolved at the base of the Bilateria + Cnidaria clade, and conserved the thumb region and Ser209 acylation site present in most other Wnts, suggesting WntA requires the core Wnt secretory pathway. WntA proteins are distinguishable from other Wnts by a synapomorphic Iso/Val/Ala216 amino-acid residue that replaces the otherwise ubiquitous Thr216 position. WntA embryonic expression is conserved between beetles and butterflies, suggesting functionality, but the WntA gene was lost three times within arthropods, in podoplean copepods, in the cyclorrhaphan fly radiation, and in ensiferan crickets and katydids. Finally, CRISPR mosaic knockouts (KOs) of porcupine and wntless phenocopied the pattern-specific effects of WntA KOs in the wings of Vanessa cardui butterflies. These results highlight the molecular conservation of the WntA protein across invertebrates, and imply it functions as a typical Wnt ligand that is acylated and secreted through the Porcupine/Wntless secretory pathway.

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.

R2R3-MYB genes control petal pigmentation patterning in Clarkia gracilis ssp. sonomensis (Onagraceae)

Petal pigmentation patterning is widespread in flowering plants. The genetics of these pattern elements has been of great interest for understanding the evolution of phenotypic diversification. Here, we investigate the genetic changes responsible for the evolution of an unpigmented petal element on a colored background. We used transcriptome analysis, gene expression assays, cosegregation in F₂ plants and functional tests to identify the gene(s) involved in petal coloration in Clarkia gracilis ssp. sonomensis. We identified an R2R3-MYB transcription factor (CgsMYB12) responsible for anthocyanin pigmentation of the basal region ('cup') in the petal of C. gracilis ssp. sonomensis. A functional mutation in CgsMYB12 creates a white cup on a pink petal background. Additionally, we found that two R2R3-MYB genes (CgsMYB6 and CgsMYB11) are also involved in petal background pigmentation. Each of these three R2R3-MYB genes exhibits a different spatiotemporal expression pattern. The functionality of these R2R3-MYB genes was confirmed through stable transformation of Arabidopsis. Distinct spatial patterns of R2R3-MYB expression have created the possibility that pigmentation in different sections of the petal can evolve independently. This finding suggests that recent gene duplication has been central to the evolution of petal pigmentation patterning in C. gracilis ssp. sonomensis.

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.

Multicomponent structures in camouflage and mimicry in butterfly wing patterns

Understanding how morphological structures are built is essential for appreciating the morphological complexity and divergence of organisms. One representative case of morphological structures is the camouflage and mimicry of butterfly wing patterns. Some previous studies have questioned whether camouflage and mimicry are truly structures, considering that they rely on coloration. Nevertheless, our recent study revealed that the leaf pattern of Kallima inachus butterfly wings evolved through the combination of changes in several pigment components in a block-wise manner; it remains unclear whether such block-wise structures are common in other cases of camouflage and mimicry in butterflies and how they come about. Previous studies focused solely on a set of homologous components, termed the nymphalid ground plan. In the present study, we extended the scope of the description of components by including not only the nymphalid ground plan but also other common components (i.e., ripple patterns, dependent patterns, and color fields). This extension allowed us to analyze the combinatorial building logic of structures and examine multicomponent structures of camouflage and mimicry in butterfly wing patterns. We investigated various patterns of camouflage and mimicry (e.g., masquerade, crypsis, Müllerian mimicry, Batesian mimicry) in nine species and decomposed them into an assembly of multiple components. These structural component analyses suggested that camouflage and mimicry in butterfly wing patterns are built up by combining multiple types of components. We also investigated associations between components and the kinds of camouflage and mimicry. Several components are statistically more often used to produce specific types of camouflage or mimicry. Thus, our work provides empirical evidence that camouflage and mimicry patterns of butterfly wings are mosaic structures, opening up a new avenue of studying camouflage, and mimicry from a structural perspective.

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.

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

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

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.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

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.

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

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

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.

apterous A specifies dorsal wing patterns and sexual traits in butterflies

Butterflies have evolved different colour patterns on their dorsal and ventral wing surfaces to serve different signalling functions, yet the developmental mechanisms controlling surface-specific patterning are still unknown. Here, we mutate both copies of the transcription factor apterous in Bicyclus anynana butterflies using CRISPR/Cas9 and show that apterous A, expressed dorsally, functions both as a repressor and modifier of ventral wing colour patterns, as well as a promoter of dorsal sexual ornaments in males. We propose that the surface-specific diversification of wing patterns in butterflies proceeded via the co-option of apterous A or its downstream effectors into various gene regulatory networks involved in the differentiation of discrete wing traits. Further, interactions between apterous and sex-specific factors such as doublesex may have contributed to the origin of sexually dimorphic surface-specific patterns. Finally, we discuss the evolution of eyespot number diversity in the family Nymphalidae within the context of developmental constraints due to apterous regulation.

Aristaless Controls Butterfly Wing Color Variation Used in Mimicry and Mate Choice

Neotropical Heliconius butterflies display a diversity of warningly colored wing patterns, which serve roles in both Müllerian mimicry and mate choice behavior. Wing pattern diversity in Heliconius is controlled by a small number of unlinked, Mendelian "switch" loci [1]. One of these, termed the K locus, switches between yellow and white color patterns, important mimicry signals as well as mating cues [2-4]. Furthermore, mate preference behavior is tightly linked to this locus [4]. K controls the distribution of white versus yellow scales on the wing, with a dominant white allele and a recessive yellow allele. Here, we combine fine-scale genetic mapping, genome-wide association studies, gene expression analyses, population and comparative genomics, and genome editing with CRISPR/Cas9 to characterize the molecular basis of the K locus in Heliconius and to infer its evolutionary history. We show that white versus yellow color variation in Heliconius cydno is due to alternate haplotypes at a putative cis-regulatory element (CRE) downstream of a tandem duplication of the homeodomain transcription factor aristaless. Aristaless1 (al1) and aristaless2 (al2) are differentially regulated between white and yellow wings throughout development with elevated expression of al1 in developing white wings, suggesting a role in repressing pigmentation. Consistent with this, knockout of al1 causes white wings to become yellow. The evolution of wing color in this group has been marked by retention of the ancestral yellow color in many lineages, a single origin of white coloration in H. cydno, and subsequent introgression of white color from H. cydno into H. melpomene.

A new A-P compartment boundary and organizer in holometabolous insect wings

Decades of research on the highly modified wings of Drosophila melanogaster has suggested that insect wings are divided into two Anterior-Posterior (A-P) compartments separated by an axis of symmetry. This axis of symmetry is created by a developmental organizer that establishes symmetrical patterns of gene expression that in turn pattern the A-P axis of the wing. Butterflies possess more typical insect wings and butterfly wing colour patterns provide many landmarks for studies of wing structure and development. Using eyespot colour pattern variation in Vanessa butterflies, here we show an additional A-P axis of symmetry running between wing sectors 3 and 4. Boundaries of Drosophila mitotic clones suggest the existence of a previously undetected Far-Posterior (F-P) compartment boundary that coincides with this additional A-P axis. A similar compartment boundary is evident in butterfly mosaic gynandromorphs. We suggest that this additional compartment boundary and its associated developmental organizer create an axis of wing colour pattern symmetry and a gene expression-based combinatorial code, permitting each insect wing compartment to acquire a unique identity and allowing for the individuation of butterfly eyespots.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity

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

Mimicry in butterflies: co-option and a bag of magnificent developmental genetic tricks

Butterfly wing patterns are key adaptations that are controlled by remarkable developmental and genetic mechanisms that facilitate rapid evolutionary change. With swift advancements in the fields of genomics and genetic manipulations, identifying the regulators of wing development and mimetic wing patterns has become feasible even in nonmodel organisms such as butterflies. Recent mapping and gene expression studies have identified single switch loci of major effects such as transcription factors and supergenes as the main drivers of adaptive evolution of mimetic and polymorphic butterfly wing patterns. We highlight several of these examples, with emphasis on doublesex, optix, WntA and other dynamic, yet essential, master regulators that control critical color variation and sex-specific traits. Co-option emerges as a predominant theme, where typically embryonic and other early-stage developmental genes and networks have been rewired to regulate polymorphic and sex-limited mimetic wing patterns in iconic butterfly adaptations. Drawing comparisons from our knowledge of wing development in Drosophila, we illustrate the functional space of genes that have been recruited to regulate butterfly wing patterns. We also propose a developmental pathway that potentially results in dorsoventral mismatch in butterfly wing patterns. Such dorsoventrally mismatched color patterns modulate signal components of butterfly wings that are used in intra- and inter-specific communication. Recent advances-fuelled by RNAi-mediated knockdowns and CRISPR/Cas9-based genomic edits-in the developmental genetics of butterfly wing patterns, and the underlying biological diversity and complexity of wing coloration, are pushing butterflies as an emerging model system in ecological genetics and evolutionary developmental biology. WIREs Dev Biol 2018, 7:e291. doi: 10.1002/wdev.291 This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Comparative Development and Evolution > Regulation of Organ Diversity Comparative Development and Evolution > Evolutionary Novelties.

Giving eyespots a shiner: Pharmacologic manipulation of the Io moth wing pattern

Our knowledge of wing pattern formation in Lepidoptera has advanced significantly in recent years due to the careful examination of several groups of butterflies. The eyespot is a prominent feature of Lepidoptera wing pattern, especially in the family Saturniidae. The present study examined how sulfated polysaccharides affected the wing pattern formation of the Io moth, Automeris io (Saturniidae).  Prepupae and pupae of this species were subjected to injections of heparin and cold shock. While the cold shock had little to no effect on wing pattern, the aberrations resulting from heparin injections were moderate to profound and depended on the dose and the stage at which injection was made. The changes consisted of expansion of the black ring around the dorsal hindwing eyespots and distortion of discal spots on both dorsal and ventral sides of forewings, suggesting a possible link between genetic controls of these elements. Several different types of scales form the normal color pattern of Automeris io, and heparin-induced changes correspond to changes in shape of scales. The resulting aberrations are dubbed ‘Black Eye’ and ‘Comet Eye.’ Other known aberrations of Automeris io eyespots are summarized, illustrated, and named.

Wingless is a positive regulator of eyespot color patterns in Bicyclus anynana butterflies

Eyespot patterns of nymphalid butterflies are an example of a novel trait yet, the developmental origin of eyespots is still not well understood. Several genes have been associated with eyespot development but few have been tested for function. One of these genes is the signaling ligand, wingless, which is expressed in the eyespot centers during early pupation and may function in eyespot signaling and color ring differentiation. Here we tested the function of wingless in wing and eyespot development by down-regulating it in transgenic Bicyclus anynana butterflies via RNAi driven by an inducible heat-shock promoter. Heat-shocks applied during larval and early pupal development led to significant decreases in wingless mRNA levels and to decreases in eyespot size and wing size in adult butterflies. We conclude that wingless is a positive regulator of eyespot and wing development in B. anynana butterflies.

The wing pattern of Moerarchis Durrant, 1914 (Lepidoptera: Tineidae) clarifies transitions between predictive models

The evolution of wing pattern in Lepidoptera is a popular area of inquiry but few studies have examined microlepidoptera, with fewer still focusing on intraspecific variation. The tineid genus Moerarchis Durrant, 1914 includes two species with high intraspecific variation of wing pattern. A subset of the specimens examined here provide, to my knowledge, the first examples of wing patterns that follow both the 'alternating wing-margin' and 'uniform wing-margin' models in different regions along the costa. These models can also be evaluated along the dorsum of Moerarchis, where a similar transition between the two models can be seen. Fusion of veins is shown not to effect wing pattern, in agreement with previous inferences that the plesiomorphic location of wing veins constrains the development of colour pattern. The significant correlation between wing length and number of wing pattern elements in Moerarchis australasiella shows that wing size can act as a major determinant of wing pattern complexity. Lastly, some M. australasiella specimens have wing patterns that conform entirely to the 'uniform wing-margin' model and contain more than six bands, providing new empirical insight into the century-old question of how wing venation constrains wing patterns with seven or more bands.

Comparative transcriptomics of anal fin pigmentation patterns in cichlid fishes

Background

Understanding the genetic basis of novel traits is a central topic in evolutionary biology. Two novel pigmentation phenotypes, egg-spots and blotches, emerged during the rapid diversification of East African cichlid fishes. Egg-spots are circular pigmentation markings on the anal fins of hundreds of derived haplochromine cichlids species, whereas blotches are patches of conspicuous anal fin pigmentation with ill-defined boundaries that occur in few species that belong to basal cichlid lineages. Both traits play an important role in the breeding behavior of this group of fishes. Knowledge about the origin, homology and underlying genetics of these pigmentation traits is sparse.

Results

Here, we present a comparative transcriptomic and differential gene expression analysis of egg-spots and blotches. We first conducted an RNA sequencing experiment where we compared egg-spot tissue with the remaining portion of egg-spot-free fin tissue using six individuals of Astatotilapia burtoni. We identified 1229 differentially expressed genes between the two tissue types. We then showed that rates of evolution of these genes are higher than average estimated on whole transcriptome data. Using quantitative real-time PCR, we found that 29 out of a subset of 46 differentially expressed genes showed an analogous expression pattern in another haplochromine species’ egg-spots, Cynotilapia pulpican, strongly suggesting that these genes are involved in the egg-spot phenotype. Among these are the previously identified egg-spot gene fhl2a, two known patterning genes (hoxC12a and bmp3) as well as other pigmentation related genes such as asip. Finally, we analyzed the expression patterns of the same gene subset in two species that feature blotches instead of egg-spots, one haplochromine species (Pseudocrenilabrus philander) and one ectodine species (Callochromis macrops), revealing that the expression patterns in blotches and egg-spots are rather distinct.

Conclusions

We identified several candidate genes that will serve as an important and useful resource for future research on the emergence and diversification of cichlid fishes’ egg-spots. Only a limited degree of conservation of gene expression patterns was detected between the egg-spots of the derived haplochromines and blotches from ancestral haplochromines, as well as between the two types of blotches, suggesting an independent origin of these traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3046-y) contains supplementary material, which is available to authorized users.

The phenotypic plasticity of developmental modules

Background

Organisms develop and evolve in a modular fashion, but how individual modules interact with the environment remains poorly understood. Phenotypically plastic traits are often under selection, and studies are needed to address how traits respond to the environment in a modular fashion. In this study, tissue-specific plasticity of melanic spots was examined in the large milkweed bug, Oncopeltus fasciatus.

Results

Although the size of the abdominal melanic bands varied according to rearing temperatures, wing melanic bands were more robust. To explore the regulation of abdominal pigmentation plasticity, candidate genes involved in abdominal melanic spot patterning and biosynthesis of melanin were analyzed. While the knockdown of dopa decarboxylase (Ddc) led to lighter pigmentation in both the wings and the abdomen, the shape of the melanic elements remained unaffected. Although the knockdown of Abdominal-B (Abd-B) partially phenocopied the low-temperature phenotype, the abdominal bands were still sensitive to temperature shifts. These observations suggest that regulators downstream of Abd-B but upstream of DDC are responsible for the temperature response of the abdomen. Ablation of wings led to the regeneration of a smaller wing with reduced melanic bands that were shifted proximally. In addition, the knockdown of the Wnt signaling nuclear effector genes, armadillo1 and armadillo 2, altered both the melanic bands and the wing shape. Thus, the pleiotropic effects of Wnt signaling may constrain the amount of plasticity in wing melanic bands.

Conclusions

We propose that when traits are regulated by distinct pre-patterning mechanisms, they can respond to the environment in a modular fashion, whereas when the environment impacts developmental regulators that are shared between different modules, phenotypic plasticity can manifest as a developmentally integrated system.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0053-7) contains supplementary material, which is available to authorized users.

Forewing color pattern in Micropterigidae (Insecta: Lepidoptera): homologies between contrast boundaries, and a revised hypothesis for the origin of symmetry systems

Background

Despite the great importance of lepidopteran wing patterns in various biological disciplines, homologies between wing pattern elements in different moth and butterfly lineages are still not understood. Among other reasons, this may be due to an incomplete understanding of the relationship between color pattern and wing venation; many individual wing pattern elements have a known relationship with venation, but a framework to unite all wing pattern elements with venation is lacking. Though plesiomorphic wing veins are known to influence color patterning even when not expressed in the adult wing, most studies of wing pattern evolution have focused on derived taxa with a reduced suite of wing veins.

Results

The present study aims to address this gap through an examination of Micropterigidae, a very early-diverged moth family in which all known plesiomorphic lepidopteran veins are expressed in the adult wing. The relationship between wing pattern and venation was examined in 66 species belonging to 9 genera. The relationship between venation and pattern element location, predicted based on moths in the family Tortricidae, holds for Sabatinca just as it does for Micropterix. However, the pattern elements that are lightly colored in Micropterix are dark in Sabatinca, and vice-versa. When plotted onto a hypothetical nymphalid wing in accordance with the relationship between pattern and venation discussed here, the wing pattern of Sabatinca doroxena very closely resembles the nymphalid groundplan.

Conclusions

The color difference in pattern elements between Micropterix and Sabatinca indicates that homologies exist among the contrast boundaries that divide wing pattern elements, and that color itself is not a reliable indicator of homology. The similarity between the wing pattern of Sabatinca doroxena and the nymphalid groundplan suggests that the nymphalid groundplan may have originated from a Sabatinca-like wing pattern subjected to changes in wing shape and reduced expression of venation.

Electronic supplementary material

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

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.

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.

Color Pattern on the Forewing of Micropterix (Lepidoptera: Micropterigidae): Insights into the Evolution of Wing Pattern and Wing Venation in Moths

Wing patterns are key taxonomic characters that have long been used in descriptions of Lepidoptera; however, wing pattern homologies are not understood among different moth lineages. Here, we examine the relationship between wing venation and wing pattern in the genus Micropterix, among the most basal extant Lepidoptera, in order to evaluate the two existing predictive models that have the potential to establish wing pattern element homologies for the order. The location of wing pattern elements along the costal margin of the wing in Micropterix is consistent with the predictions of the model proposed for Tortricidae by Brown and Powell in 1991, later modified by Baixeras in 2002. The predictive power of this model for such distantly related taxa suggests that the model may hold across various superfamilies within Lepidoptera, and supports the long-held notion that fasciae, not spots, are the most likely primitive wing pattern elements for the order. In addition, the location of wing pattern elements suggests that the wing vein commonly termed Sc1 may in fact be a different vein, which Comstock identified in Trichoptera and referred to as "a."

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.

Integrating evo-devo with ecology for a better understanding of phenotypic evolution

Evolutionary developmental biology (evo-devo) has provided invaluable contributions to our understanding of the mechanistic relationship between genotypic and phenotypic change. Similarly, evolutionary ecology has greatly advanced our understanding of the relationship between the phenotype and the environment. To fully understand the evolution of organismal diversity, a thorough integration of these two fields is required. This integration remains highly challenging because model systems offering a rich ecological and evolutionary background, together with the availability of developmental genetic tools and genomic resources, are scarce. In this review, we introduce the semi-aquatic bugs (Gerromorpha, Heteroptera) as original models well suited to study why and how organisms diversify. The Gerromorpha invaded water surfaces over 200 mya and diversified into a range of remarkable new forms within this new ecological habitat. We summarize the biology and evolutionary history of this group of insects and highlight a set of characters associated with the habitat change and the diversification that followed. We further discuss the morphological, behavioral, molecular and genomic tools available that together make semi-aquatic bugs a prime model for integration across disciplines. We present case studies showing how the implementation and combination of these approaches can advance our understanding of how the interaction between genotypes, phenotypes and the environment drives the evolution of distinct morphologies. Finally, we explain how the same set of experimental designs can be applied in other systems to address similar biological questions.

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.