Butterfly eyespots evolved via cooption of an ancestral gene-regulatory network that also patterns antennae, legs, and wings.

Butterfly eyespots are beautiful novel traits with an unknown developmental origin. Here we show that eyespots likely originated via cooption of parts of an ancestral appendage gene-regulatory network (GRN) to novel locations on the wing. Using comparative transcriptome analysis, we show that eyespots cluster most closely with antennae, relative to multiple other tissues. Furthermore, three genes essential for eyespot development, Distal-less (Dll), spalt (sal), and Antennapedia (Antp), share similar regulatory connections as those observed in the antennal GRN. CRISPR knockout of cis-regulatory elements (CREs) for Dll and sal led to the loss of eyespots, antennae, legs, and also wings, demonstrating that these CREs are highly pleiotropic. We conclude that eyespots likely reused an ancient GRN for their development, a network also previously implicated in the development of antennae, legs, and wings.

The people behind the papers - Heidi Connahs, Sham Tlili, Timothy Saunders and Antónia Monteiro.

Butterfly eyespots are striking examples of animal patterning, but their developmental origins are still relatively poorly understood. A new paper in Development - the result of a collaboration between two Singapore-based labs - now combines CRISPR-Cas9 gene targeting with theoretical modelling to address the role of the Distal-less transcription factor in eyespot patterning. We caught up with co-first authors Heidi Connahs and Sham Tlili, and their respective supervisors Timothy Saunders (Assistant Professor at the Mechanobiology Institute, National University of Singapore) and Antónia Monteiro (Associate Professor at the Department of Biological Sciences, National University of Singapore and Yale-NUS College) to find out more about the story.

The pivotal role of aristaless in development and evolution of diverse antennal morphologies in moths and butterflies.

BACKGROUND: Antennae are multi-segmented appendages and main odor-sensing organs in insects. In Lepidoptera (moths and butterflies), antennal morphologies have diversified according to their ecological requirements. While diurnal butterflies have simple, rod-shaped antennae, nocturnal moths have antennae with protrusions or lateral branches on each antennal segment for high-sensitive pheromone detection. A previous study on the Bombyx mori (silk moth) antenna, forming two lateral branches per segment, during metamorphosis has revealed the dramatic change in expression of antennal patterning genes to segmentally reiterated, branch-associated pattern and abundant proliferation of cells contributing almost all the dorsal half of the lateral branch. Thus, localized cell proliferation possibly controlled by the branch-associated expression of antennal patterning genes is implicated in lateral branch formation. Yet, actual gene function in lateral branch formation in Bombyx mori and evolutionary mechanism of various antennal morphologies in Lepidoptera remain elusive. RESULTS: We investigated the function of several genes and signaling specifically in lateral branch formation in Bombyx mori by the electroporation-mediated incorporation of siRNAs or morpholino oligomers. Knock down of aristaless, a homeobox gene expressed specifically in the region of abundant cell proliferation within each antennal segment, during metamorphosis resulted in missing or substantial shortening of lateral branches, indicating its importance for lateral branch formation. aristaless expression during metamorphosis was lost by knock down of Distal-less and WNT signaling but derepressed by knock down of Notch signaling, suggesting the strict determination of the aristaless expression domain within each antennal segment by the combinatorial action of them. In addition, analyses of pupal aristaless expression in antennae with various morphologies of several lepidopteran species revealed that the aristaless expression pattern has a striking correlation with antennal shapes, whereas the segmentally reiterated expression pattern was observed irrespective of antennal morphologies. CONCLUSIONS: Our results presented here indicate the significance of aristaless function in lateral branch formation in B. mori and imply that the diversification in the aristaless expression pattern within each antennal segment during metamorphosis is one of the significant determinants of antennal morphologies. According to these findings, we propose a mechanism underlying development and evolution of lepidopteran antennae with various morphologies.

Divergent RNA Localisation Patterns of Maternal Genes Regulating Embryonic Patterning in the Butterfly Pararge aegeria.

The maternal effect genes responsible for patterning the embryo along the antero-posterior (AP) axis are broadly conserved in insects. The precise function of these maternal effect genes is the result of the localisation of their mRNA in the oocyte. The main developmental mechanisms involved have been elucidated in Drosophila melanogaster, but recent studies have shown that other insect orders often diverge in RNA localisation patterns. A recent study has shown that in the butterfly Pararge aegeria the distinction between blastodermal embryonic (i.e. germ band) and extra-embryonic tissue (i.e. serosa) is already specified in the oocyte during oogenesis in the ovariole, long before blastoderm cellularisation. To examine the extent by which a female butterfly specifies and patterns the AP axis within the region fated to be the germ band, and whether she specifies a germ plasm, we performed in situ hybridisation experiments on oocytes in P. aegeria ovarioles and on early embryos. RNA localisation of the following key maternal effect genes were investigated: caudal (cad), orthodenticle (otd), hunchback (hb) and four nanos (nos) paralogs, as well as TDRD7 a gene containing a key functional domain (OST-HTH/LOTUS) shared with oskar. TDRD7 was mainly confined to the follicle cells, whilst hb was exclusively zygotically transcribed. RNA of some of the nos paralogs, otd and cad revealed complex localisation patterns within the cortical region prefiguring the germ band (i.e. germ cortex). Rather interestingly, otd was localised within and outside the anterior of the germ cortex. Transcripts of nos-O formed a distinct granular ring in the middle of the germ cortex possibly prefiguring the region where germline stem cells form. These butterfly RNA localisation patterns are highly divergent with respect to other insects, highlighting the diverse ways in which different insect orders maternally regulate early embryogenesis of their offspring.

Double trouble: combined action of meiotic drive and Wolbachia feminization in Eurema butterflies.

Arthropod sex ratios can be manipulated by a diverse range of selfish genetic elements, including maternally inherited Wolbachia bacteria. Feminization by Wolbachia is rare but has been described for Eurema mandarina butterflies. In this species, some phenotypic and functional females, thought to be ZZ genetic males, are infected with a feminizing Wolbachia strain, wFem. Meanwhile, heterogametic WZ females are not infected with wFem. Here, we establish a quantitative PCR assay allowing reliable sexing in three Eurema species. Against expectation, all E. mandarina females, including wFem females, had only one Z chromosome that was paternally inherited. Observation of somatic interphase nuclei confirmed that W chromatin was absent in wFem females, but present in females without wFem. We conclude that the sex bias in wFem lines is due to meiotic drive (MD) that excludes the maternal Z and thus prevents formation of ZZ males. Furthermore, wFem lines may have lost the W chromosome or harbour a dysfunctional version, yet rely on wFem for female development; removal of wFem results in all-male offspring. This is the first study that demonstrates an interaction between MD and Wolbachia feminization, and it highlights endosymbionts as potentially confounding factors in MD of sex chromosomes.

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.

Single locus affects embryonic segment polarity and multiple aspects of an adult evolutionary novelty.

BACKGROUND: The characterization of the molecular changes that underlie the origin and diversification of morphological novelties is a key challenge in evolutionary developmental biology. The evolution of such traits is thought to rely largely on co-option of a toolkit of conserved developmental genes that typically perform multiple functions. Mutations that affect both a universal developmental process and the formation of a novelty might shed light onto the genetics of traits not represented in model systems. Here we describe three pleiotropic mutations with large effects on a novel trait, butterfly eyespots, and on a conserved stage of embryogenesis, segment polarity. RESULTS: We show that three mutations affecting eyespot size and/or colour composition in Bicyclus anynana butterflies occurred in the same locus, and that two of them are embryonic recessive lethal. Using surgical manipulations and analysis of gene expression patterns in developing wings, we demonstrate that the effects on eyespot morphology are due to changes in the epidermal response component of eyespot induction. Our analysis of morphology and of gene expression in mutant embryos shows that they have a typical segment polarity phenotype, consistent with the mutant locus encoding a negative regulator of Wingless signalling. CONCLUSIONS: This study characterizes the segregation and developmental effects of alleles at a single locus that controls the morphology of a lineage-specific trait (butterfly eyespots) and a conserved process (embryonic segment polarity and, specifically, the regulation of Wingless signalling). Because no gene with such function was found in the orthologous, highly syntenic genomic regions of two other lepidopterans, we hypothesize that our locus is a yet undescribed, possibly lineage-specific, negative regulator of the conserved Wnt/Wg pathway. Moreover, the fact that this locus interferes with multiple aspects of eyespot morphology and maps to a genomic region containing key wing pattern loci in different other butterfly species suggests it might correspond to a 'hotspot' locus in the diversification of this novel trait.

Reproductive plasticity, ovarian dynamics and maternal effects in response to temperature and flight in Pararge aegeria.

In nature, ovipositing females may be subjected to multiple extrinsic and intrinsic environmental factors simultaneously. To adequately assess a species response to environmental conditions during oviposition it may therefore be necessary to consider the interaction between multiple intrinsic and extrinsic factors simultaneously. Using the butterfly, Pararge aegeria, this study examined the combined effects of extrinsic (temperature and flight) and intrinsic (body mass and age) factors on ovarian dynamics, egg provisioning and reproductive output, and explored how these effects subsequently influenced offspring fitness when egg-stage development occurred in a low humidity environment. Both temperature- and flight-mediated plasticity in female reproductive output was observed, and there were strong temperature by flight interaction effects for the traits oocyte size and egg mass. As females aged, mean daily fecundity differed across temperature treatments, but not across flight treatments. Overall, temperature had more pronounced effects on ovarian dynamics than flight. Flight mainly influenced egg mass via changes in relative water content. A mismatch between the physiological response of females to high temperature and the requirements of their offspring had a negative impact on offspring fitness via effects on egg hatching success.

Energetics of embryonic development: effects of temperature on egg and hatchling composition in a butterfly.

Phenotypic plasticity may allow an organism to adjust its phenotype to environmental needs. However, little is known about environmental effects on offspring biochemical composition and turnover rates, including energy budgets and developmental costs. Using the tropical butterfly Bicyclus anynana and employing a full-factorial design with two oviposition and two developmental temperatures, we explore the consequences of temperature variation on egg and hatchling composition, and the associated use and turnover of energy and egg compounds. At the lower temperature, larger but fewer eggs were produced. Larger egg sizes were achieved by provisioning these eggs with larger quantities of all compounds investigated (and thus more energy), whilst relative egg composition was rather similar to that of smaller eggs laid at the higher temperature. Turnover rates during embryonic development differed across developmental temperatures, suggesting an emphasis on hatchling quality (i.e. protein content) at the more stressful lower temperature, but on storage reserves (i.e. lipids) at the higher temperature. These differences may represent adaptive maternal effects. Embryonic development was much more efficient at the lower temperature, providing a possible mechanism underlying the temperature-size rule.

Butterfly genomics eclosing.

Technological and conceptual advances of the last decade have led to an explosion of genomic data and the emergence of new research avenues. Evolutionary and ecological functional genomics, with its focus on the genes that affect ecological success and adaptation in natural populations, benefits immensely from a phylogenetically widespread sampling of biological patterns and processes. Among those organisms outside established model systems, butterflies offer exceptional opportunities for multidisciplinary research on the processes generating and maintaining variation in ecologically relevant traits. Here we highlight research on wing color pattern variation in two groups of Nymphalid butterflies, the African species Bicyclus anynana (subfamily Satyrinae) and species of the South American genus Heliconius (subfamily Heliconiinae), which are emerging as important systems for studying the nature and origins of functional diversity. Growing genomic resources including genomic and cDNA libraries, dense genetic maps, high-density gene arrays, and genetic transformation techniques are extending current gene mapping and expression profiling analysis and enabling the next generation of research questions linking genes, development, form, and fitness. Efforts to develop such resources in Bicyclus and Heliconius underscore the general challenges facing the larger research community and highlight the need for a community-wide effort to extend ongoing functional genomic research on butterflies.