From Patterning Genes to Process: Unraveling the Gene Regulatory Networks That Pattern Heliconius Wings

McMillan, W. Owen; Livraghi, Luca; Concha, Carolina; Hanly, Joseph J.

doi:10.3389/fevo.2020.00221

REVIEW article

Front. Ecol. Evol., 24 July 2020

Sec. Evolutionary Developmental Biology

Volume 8 - 2020 | https://doi.org/10.3389/fevo.2020.00221

This article is part of the Research Topic Evo-Devo of Color Pattern Formation View all 20 articles

From Patterning Genes to Process: Unraveling the Gene Regulatory Networks That Pattern Heliconius Wings

$\r\nW. Owen McMillan*&#x;$ W. Owen McMillan^1*†

Luca Livraghi^1,2†

Carolina Concha^1†

Joseph J. Hanly^1,3†

¹Smithsonian Tropical Research Institute, Gamboa, Panama
²Department of Zoology, University of Cambridge, Cambridge, United Kingdom
³Department of Biological Sciences, The George Washington University, Washington, DC, United States

Butterfly wing patterns have emerged as exceptional model systems with which to link the developmental and genetic processes that generate morphological variation with the ecological and evolutionary processes that shape this variation in natural populations. Among butterflies, research on species within the genus Heliconius has provided remarkable opportunities to explore how phenotypic diversity is generated within the context of an extraordinary adaptive radiation. Wing pattern diversity among the 48 species and hundreds of intraspecific variants arose within the last 12–14 million years and includes striking pattern convergence between distantly related species, as well as marked pattern divergence between closely related populations and species. Here, we synthesize recent research aimed at gaining a mechanistic understanding of how this variation is generated. This research integrates decades of controlled crossing experiments, and the discovery of major wing patterning genes (optix, aristaless1, WntA and cortex) with recent functional genetic manipulation using CRISPR/Cas9 targeted mutagenesis. The emerging data provides a rich framework with which to explore the repeatability of evolution, particularly within the context of how natural selection acts on divergent gene regulatory networks to generate both highly convergent, as well as highly divergent phenotypes. Overall, the functional data show that the gene regulatory networks underlying pattern variation diverge rapidly in Heliconius; yet these networks retain enough flexibility so that natural selection can drive the evolution of nearly identical patterns from different developmental genetic starting points. Moreover, for the first time this research is starting to illuminate the links between the genetic changes modulating pattern variation and how they influence the larger gene networks that are ultimately responsible for patterning a butterfly wing. There are still large gaps in our understanding, but current research priorities are well laid out and experimental methodologies are in place to address them. The challenge is to synthesize diverse research strategies into a cohesive picture of morphological evolution.

“…on these expanded membranes, nature writes, as on a tablet, the story of the modifications of species.” Henry Walter Bates, in “Naturalist on the river Amazons” (1864).

Background and Overview

On the 28th of May 1848, Henry Walter Bates, then only 23 years old, stepped off the boat in Para, Brazil with his friend Alfred Russell Wallace. Bates was to spend the next 11 years exploring the Amazon Basin and his collections, particularly his butterfly collections, provided some of the most visually intuitive examples of evolution by natural selection (Bates, 1864). During his explorations, Bates was struck by the remarkable similarity in the vivid wing patterns of Heliconius and other butterfly and moth species and how these patterns would change rapidly as he traveled across the region, as if “touched by an enchanter’s wand” (Bates, 1879). His insights led to the development of mimicry theory, whereby natural selection drives the evolution of phenotypic similarity (Pasteur, 1982). In the intervening 150 years, research has started to identify the genetic and developmental processes responsible for the evolution of butterfly wing patterns generally, including the mimetic wing patterns of Heliconius. For this reason, the study of butterfly wing patterns continues to speak to research questions at the forefront of biology today, much as it did during that early dawn of evolutionary theory.

Here we review recent research into the genetic and developmental basis of wing pattern change. We focus on the vivid wing patterns of Heliconius, which are used both as aposematic warning signals and as conspecific mating cues (Benson, 1972; Mallet et al., 1990; Jiggins et al., 2001; Kapan, 2001; Chouteau et al., 2017; Merrill et al., 2019). This combination of natural and sexual selection is associated with remarkable variation in wing patterns both within and between species, including numerous cases of parallel evolution of similar wing patterns between distantly related species (Figure 1). Classical crossing work done in the 1950s and 1960s by Beebe, Crane, Sheppard, Turner, and others (Turner, 1972; Sheppard et al., 1985; Beebe, 2005) established that pattern variation within a specific Heliconius species had a relatively simple genetic basis, with most variation explained by a small number of autosomal loci. Later, linkage mapping and positional cloning demonstrated that the same suite of major patterning loci are modulating pattern variation broadly across the genus and that variation around these loci is responsible for both convergent change between species and divergent change within species (Joron et al., 2006; Kapan et al., 2006; Kronforst et al., 2006; Huber et al., 2015). The advent of more powerful genomic technologies culminated in the identification of specific genes that modulate pattern variation across the genus (Reed et al., 2011; Martin et al., 2012; Nadeau et al., 2016; Westerman et al., 2018). These loci join a growing list shown to play a general role in wing pattern variation in butterflies and moths (Table 1).

FIGURE 1

Figure 1. Evolution’s canvas: Wing pattern mimicry among species in the melpomene/cydno/silvaniform clade (MCS) (above line) and the erato/sapho/sara clade (ESS) (below line). The two lineages diverged from each other about 12 MYA (Kozak et al., 2015), yet have recently evolved a complement of identically patterned wings. Species in the MCS clade are: (A) H. hecale zuleika, (B) H. atthis, (C) H. melpomene ecuadorensis, (D) H. pachinus, (E) H. cydno chioneus, (F) H. m. plesseni, (G) H. m. cythera, (H), H. m. rosina. Species in the ESS clade are: (J) H. sapho sapho, (K) H. erato notabilis, (L) H. e. cyrbia, (M) H. e. demophoon, (N) H. sara sara, (O) H. peruvianus, (P) H. e. etylus, (Q) H. hewitsoni, (R) H. hecalesia formosus. Heliconius doris (I) is from a third clade that falls between the ESS and MCS lineages.

TABLE 1

Table 1. Wing genes that have been identified by linkage or association as involved with development of wing color or pattern in butterflies and moths.

The discovery of “wing patterning” genes in Heliconius allows us to make broad general conclusions about how morphological diversity is generated. Foremost, patterning loci identified in Heliconius have also been shown to affect pattern variation broadly across butterflies and even moths. For example, the gene cortex appears to be a major target of evolution across the lepidoptera (Nadeau et al., 2016; van’t Hof et al., 2016). In Heliconius, the locus acts across both wing surfaces and underlies variation in white, yellow, orange, and black pattern elements (Joron et al., 2006). Cortex has also been implicated in pattern variation in other butterflies (VanKuren et al., 2019), the peppered moth Biston betularia and other geometrids (van’t Hof et al., 2016, 2019), as well as the silkmoth Bombyx mori (Ito et al., 2016). Secondly, nearly all the patterning genes thus far identified are highly conserved at the amino acid level and are used extensively throughout embryonic and larval development, and later simply redeployed on a developing wing to affect pattern variation. The transcription factor optix, for example, was first discovered in Drosophila, where ectopic expression of the gene in the developing antennae generates ectopic eyes (Seimiya and Gehring, 2000). It is also expressed in the Drosophila pupal wing and is involved with vein patterning there (Organista et al., 2015). In Heliconius, optix is expressed in the optic lobe (Martin et al., 2014), where it presumably plays a role in neural development. However, in Heliconius, optix has also evolved new expression domains and is expressed in the developing pupal wing from about 60 h after pupation, where its expression perfectly prefigures red pattern variation on Heliconius wings (Reed et al., 2011). Thirdly, pattern variation is produced mostly by regulating when and where these genes are expressed during development. The gene WntA is a member of the Wnt family of signaling ligands and has similarly evolved new patterning roles, with dramatic shifts in WntA expression fueling pattern variation in Heliconius (Martin et al., 2012; Concha et al., 2019) and other nymphalid butterflies (Gallant et al., 2014; Martin and Reed, 2014; Mazo-Vargas et al., 2017).

This review synthesizes recent research aimed at gaining a more mechanistic understanding of how wing color pattern variation in Heliconius is generated. We touch only briefly on the biology and ecology of the group that provides the context for understanding their remarkable pattern evolution. This aspect of Heliconius has been reviewed many times (see Supple et al., 2014; Kronforst and Papa, 2015), culminating in a recent book (Jiggins, 2017). Similarly, we do not review the body of research that led to the identification of the four major patterning genes. This was a sustained effort going back more than 50 years, which, coupled with ideas from evolutionary developmental genetics (Brakefield and French, 1999; Nijhout, 2001; Carroll, 2008), allowed the remarkable progress that has been made over the last 3 years. Rather, we focus on the research that is using the discovery of key patterning loci to begin to unravel the regulatory networks that lay at the heart of pattern diversity. Specifically, we outline how one can cast the origin of pattern diversity on a butterfly wing as a basic question of scale cell fate, and then discuss how recent experiments on Heliconius are beginning to illuminate the role that major patterning loci play in the determination of scale cell fate. Finally, we highlight the research challenges and opportunities moving forward.

Wing Patterns as a Model for the Evolution of Scale Cell Fate

The identification of the patterning loci that have been the targets of natural selection in Heliconius is a major research step; however, it is important not to lose sight of the fact that genes do not make phenotypes directly (Gawne et al., 2018). Rather the products of genes interact with each other, and with parts of the genome and the environment to form complex webs of interactions that unfold through development (Rudel and Sommer, 2003). In the case of butterfly wing patterns, phenotype is realized in the patterned differentiation of wing epithelial cells into specific scale cell color types (Figure 2). Understanding how butterfly wing patterns form will require us to dissect how these cell types are specified, how they become spatially arranged, and how they differentiate.

FIGURE 2

Figure 2. Pattern formation on a butterfly wing: (A) In the late larval stage, the presumptive wing imaginal disk is patterned with lacunae (which will form the veins) indicating that the spatial arrangement of the wing has already been patterned. Aside from this, the wing consists of two opposing sheets of undifferentiated epithelial cells, each attached to the basement membrane. Unlike in Drosophila, this wing disk can be straightforwardly mapped to the adult wing and will not undergo major morphometric rearrangement. The expression domains of several transcription factors and signaling molecules are known at this stage in butterflies (summarized in Hanly et al., 2019), and notably, despite scale cells having not been specified yet, the ligand WntA is already expressed in domains that prefigure adult wing pattern. (B) At pupation, the wing disk expands rapidly. During the first 18h of pupal development, scale precursor cells differentiate from their neighbors by notch-delta mediated lateral inhibition (Reed, 2004). Each scale precursor is surrounded by six non-scale epithelial cells, equating to a 3:1 epithelia to scale ratio in the wing. These scale precursor then undergo a series of divisions and apoptosis events that result in a scale cell, which will form the blade of the scale and a socket cell, which will form a supportive structure to hold the scale in place; this indicates homology between the butterfly scale and the D. melanogaster bristle, which is supported by the expression of achaete-scute in the scale precursor cell. This process also requires the presence of planar cell polarity in the tissue. (C) After the generation of the scale and socket cell, the scale begins to send out a large, actin-driven cytoplasmic protrusion which will form the blade of the scale. In Heliconius, the presumptive Type I scales begin this process earlier than their neighbors, which indicates that there is developmental heterochrony between different scale cell types. The expression of additional genes have by now been patterned to prefigure adult wing pattern; for example the gene optix is expressed specifically in cells that will become Type III red scales in Heliconius. (D) Selector genes such as optix and aristaless1 in turn induce the expression of pigment biosynthesis in specific regions across the developing wing. (E) In mid to late pupal development, scale cell differentiation occurs across the wing, including morphological differentiation and pigment synthesis and deposition. Pigment deposition occurs in a specific temporal manner, with red ommochrome appearing first, followed by melanin and finally deposition of 3-OHK (yellow) in Type I scales. (F) This results in the adult wing pattern at eclosion.

The wing patterns of butterflies and moths are composed of overlapping wing scales, which are homologous to insect bristles (Overton, 1966; Galant et al., 1998; Zhou et al., 2009). Each scale cell generally produces a single pigment and has an elaborate nanostructure. It is the combination of pigment and ultrastructural properties of the scale cells that creates the wing colors we observe (Nijhout, 1991). The fantastic array of complex patterns that we find in nature are generated by varying how different scale cell types are arranged across the wing surface. Importantly, the development of pattern on a butterfly wing is largely devoid of processes such as cell migration which complicate inference, and the developing wing tissue can be mapped accurately onto the corresponding adult phenotype (Nijhout, 1991). We therefore propose that the process of coloring a butterfly wing can be broken into a series of spatially programmed decisions, where progenitor cells acquire specific fates through the interplay of cell signaling and differential gene expression (Figure 2). Tweaks to these developmental decisions lead to the variation in wing patterns seen both within and between the 20,000 species of butterflies (Kristensen et al., 2007).

The process of pattern formation occurs throughout wing development and is linked to the processes underlying scale cell development. In the final instar of butterfly larvae, developing wing disks increase their size through cell proliferation and begin to acquire the shape and venation that is characteristic of the adult wing (Figure 2). At this point the wings are composed of a basal lamina straddled by a single cell epithelial layer on both the ventral and dorsal surfaces (Greenstein, 1972; Ghiradella and Radigan, 1976; Nardi and Magee-Adams, 1986; Dinwiddie et al., 2014). Developing cells at this stage are still composed of undifferentiated epithelia, but the inductive signals required for scale cell identity are being laid out through the expression of transcription factors and signaling molecules (summarized in Hanly et al., 2019). It is during the pupal stage that the morphological development of the wing scales begins. Shortly after pupation, Notch signaling induces scale cell differentiation through lateral inhibition resulting in an array of aligned, parallel cell rows (Reed and Serfas, 2004). Cytological studies have shown that these cells in turn undergo two rounds of cell division promoting a neural lineage and a scale cell lineage (Stossberg, 1938). Neural precursor cells undergo apoptosis, whereas the scale cell precursors undergo a second round of cell division forming a socket cell from which a scale cell arises through actin bundle elongation (Dinwiddie et al., 2014). During pupal development, the scale cells grow dramatically in size and after about 60 h, scale cells fated to become a specific color can often be distinguished (Reed et al., 2008). In our model, during this period all scale cells must acquire a cell identity through the patterned expression of major patterning genes, which in turn will trigger the expression of further differentiation factors or terminal effector genes. Cells must also communicate with each other so that the correct edges and boundaries to delimit pattern elements are appropriately arranged. By 100 h of pupal development, we know that all scale cells are undergoing actin-dependent laminar extension and depositing cuticle, including scale type-specific gene expression of pigment synthesis pathways. By this stage, the wing patterning process is complete.

Heliconius as a Model for Understanding the Gene Regulatory Networks (GRNs) Responsible for Pattern Formation

The extraordinary diversity in wing patterns among the 48 Heliconius species arose within the last 12–14 million years (Kozak et al., 2015) and includes remarkable pattern convergence between distantly related species. Pattern matching is perhaps best exemplified by convergent change in two of the major lineages within the genus - the H. erato/sapho/sara (ESS) and the H. melpomene/cydno/silvaniform (MCS) species groups (see Figure 1). The two lineages diverged from a common ancestor early in the formation of the genus and have been evolving separately for the past 12 million years, diverging in many aspects of their ecology and reproductive biology (Jiggins, 2017). Nonetheless, within the last 2.5–4.5 million years the groups have converged on a similar complement of wing pattern types, including a number of cases of near-perfect wing pattern mimicry. Although hybridization and introgression play a role in the evolution of novel patterns within the ESS and MCS clades (Edelman et al., 2019), there is no evidence for hybridization between the two groups.

Within this evolutionary backdrop, much of the wing pattern variation in Heliconius is created by differences in the distribution of only three major scale cell types. Type I yellow scales contain the ommochrome precursor 3-hydroxykynurenine (3-OHK) (Brown, 1967); whereas Type I white scales lack pigment and the color is the result of the scale cell morphology (i.e., structural) (Gilbert et al., 1988). Type II scale cells are pigmented with melanin and Type III scale cells contain the ommochrome pigments xanthommatin and dihydroxanthommatin (Gilbert et al., 1988). The timing of scale cell enlargement and maturation differs between scale types, with Type I scale cells maturing first, followed by Type III scale cells. Scale cells fated to become melanic (Type II) mature last (Aymone et al., 2013). Pigment development largely follows differences in scale cell maturation. Scales fated to become yellow or white are evident as colorless patches beginning at about 6 days after pupation. At the same developmental time point, red-fated scale cells are orange and progressively darken to red over the next 24 h. After the red pigments have matured, melanic scales develop in a ‘wave’ emanating from the center of the wing; whereas yellow-fated scales remain colorless until a few hours before adult emergence. These changes in scale cell color are associated with upregulation of genes involved in the uptake and synthesis of pigment, some of which show strong pattern specific expression, while others do not (Reed et al., 2008; Ferguson and Jiggins, 2009; Hines et al., 2012).

The process of scale cell differentiation provides a useful context to begin to understand the function of the major patterning loci that have been identified in Heliconius (Figure 3). The patterning genes WntA and cortex are expressed during larval wing development well before scale cells have begun to form (Martin et al., 2012; Nadeau et al., 2016). These genes provide key aspects of the positional information signals to downstream pigment and scale cell maturation pathways. Later in pupal development (about 40 h after pupation), the Type I scale cells begin to undergo laminar extension before their Type II and III neighbors (Gilbert et al., 1988; Reed et al., 2008; Aymone et al., 2013). Type III scales mature next and begin by expressing the transcription factor optix, signaling cells to produce the red ommochrome pigments (Reed et al., 2011). Type II cells contain melanins, which are laid down later in pupal development. The white/yellow difference between Type I scale cells is due to the presence or absence of the ommochrome precursor 3-OHK. This color switch appears to be regulated by the differential expression of aristaless1 (Westerman et al., 2018), with the uptake of 3-OHK from the pupal hemolymph occurring right before adult emergence (Reed et al., 2008).

FIGURE 3

Figure 3. Major patterning switch genes in Heliconius butterflies: Most of the color pattern variation in Heliconius is the result of 4 major effect loci. Loci are unlinked and allelic variation around them is associated with both convergent and divergent change among species. Over the past decade, loci have been positionally cloned and genotype × phenotype, expression and, most recently, KOs with CRISPR/Cas9 have identified specific patterning genes. The gene cortex is located on chromosome 15 and variation around the gene is associated with differences in the yellow hindwing bar, forewing band color, and the white margin and ventral color phenotypes among different geographic morphs of H. erato and H. melpomene. In H. numata, variation around the region that contains cortex controls all variation in wing patterns and is responsible for differences in yellow, orange, and black pattern elements among sympatric color pattern morphs (Nadeau et al., 2016; Saenko et al., 2019). cortex KO results in the appearance of Type I scales across the wing surface of Heliconius butterflies (see Livraghi et al., 2020). The signaling ligand WntA is on chromosome 10 and variation around the locus is associated with differences in forewing band size, shape and position in H. erato and other Heliconius species. KO results in the expansions of color pattern elements in proximal wing regions in H. erato, but generates a wide variety of phenotypic effects across the Heliconius radiation (Concha et al., 2019). The gene optix is on chromosome 18 and variation around this gene is responsible for variation in red color pattern elements. optix KOs result in the transformation of red Type III scales to melanic Type II scales. Finally, variation around the gene aristaless1 on chromosome 1 modulates the presence or absence of yellow pigment in Type I scales. Aristaless1 KOs result in the development of yellow scales suggesting that the gene represses synthesis/deposition of the yellow pigment, 3-hydroxy-L-kynurenine.

From Patterns to Process – a Model of Wing Pattern Evolution

Substantial insights into the function that major switch genes play in pattern variation has been gained with the recent development of CRISPR/Cas9 technology in Heliconius (Mazo-Vargas et al., 2017; Zhang et al., 2017b; Westerman et al., 2018; Concha et al., 2019; Lewis et al., 2019) (Figure 3). For the first time, researchers can perform functional assays by inducing targeted mutations in coding regions as well as putative regulatory elements in order to observe the resulting mutant wing pattern phenotypes (reviewed in Livraghi et al., 2018). This research is illuminating the roles that the four major patterning genes play in pattern variation, how they interact during development to influence scale cell fate, and how gene patterning networks evolve within the context of adaptive phenotypic change.

Two major patterning genes, WntA and cortex, are expressed early in wing development and act to establish positional information for the downstream expression of scale cell selector genes. Under a simple model (Figure 4), cortex acts as a “master” regulator of scale cell identity, where scales switch between Type I and Type II/III scales depending on the expression of cortex. A certain threshold of cortex expression may initiate a differentiation cue that sets up a permissive environment for either melanization (optix negative) or red pigment synthesis (optix positive). In contrast, cells below that threshold become either yellow or white depending on the expression of aristaless1, which represses the pathway leading to the synthesis and deposition of 3-OHK pigment (Westerman et al., 2018). Allelic variation in the expression of cortex is responsible for variation in white, yellow, and red/orange patches across the radiation (Figure 3). However, the expression of the gene appears to be more complex than initially envisioned. Although previous expression studies using in situ hybridization have shown a strong correlation between the position of cortex mRNA in wing disks of fifth instar larvae and distal black pattern elements in the adult wings of H. numata and H. melpomene (Nadeau et al., 2016; Saenko et al., 2019), the color switches induced by cortex knockouts (KOs) do not seem to act in a pattern specific way in H. erato and H. melpomene (Livraghi et al., 2020). Instead, cortex KOs result in ectopic switches from melanic (Type II) and red/brown/orange (Type III) into white/yellow (Type I) scale cells across the whole surface of both wings in H. erato and H. melpomene. The switch is associated with both a shift in scale cell color and nanostructure. Thus, the loss of function of cortex causes a wing scale cell to change its identity and fully acquire the identity of a Type I scale cell. The broad effects of cortex across the wings of Heliconius underscore the potential for this gene to generate much wider pattern variation than other major effects genes. This wide range of effects is exemplified by the polymorphic H. numata, where variation in yellow, white, black and orange elements all map to this locus (Joron et al., 2006; Saenko et al., 2010; Nadeau et al., 2016).

FIGURE 4

Figure 4. Bifurcating model of interactions among color pattern loci: During wing development, scales express key effector genes that lead to alternate scale cell types. Early presumptive scale cells (PSCs) express cell specification genes such as cortex, which initiate differentiation into Type II (optix−) or Type III (optix +) scales. In the absence of cortex, scale cells differentiate into Type I scales, which differ in pigmentation state based on 3-OHK synthesis controlled by aristaless1 expression. Under this model, WntA acts as a landscape modifier, whereby Wnt signaling modifies the trans environment experienced by differentiating scale cells, in turn delineating the boundaries in which specific differentiation factors can act. Model inspired by Larry Gilbert and collaborators (Gilbert et al., 1988).

If cortex is acting as a “on-off” switch among scale cell types, the signaling ligand WntA might be best conceptualized as modifying the overall regulatory landscape used to provide the positional information for the establishment of boundaries between scale cell types. In contrast to cortex, WntA effects are limited to specific wing pattern elements and the extent of pattern induction across the wing varies both within and between species in the ESS/MCS radiation (Figure 5). The key role of WntA in orchestrating scale cell identity is illustrated by the observation that the shift in scale cell fate happens in every possible direction–from black to red, yellow, and white; as well as from any given color to black, and also from red to yellow scales (Figure 5A). Thus, WntA appears to regulate scale cell fate in a more flexible manner than cortex. Furthermore, unlike cortex where scale cells revert to Type I in cortex KO cells, the exact scale cell phenotype induced by a WntA signal is highly context-dependent (Concha et al., 2019). The complexity of these interactions is best demonstrated in the different responses of optix to WntA KOs in divergent color pattern morphs of H. erato (Figure 5B). Within a given genetic background, WntA KOs result in an upregulation of optix expression and concomitant development of red scales. However, the opposite effect is observed in yellow-banded morphs of H. erato, where the existing yellow patterns expand, even across red scales which are known to express optix in wild-type butterflies.

FIGURE 5

Figure 5. Chance and contingency on a Heliconius wing: CRISPR/Cas9 experiments demonstrate the pervasive role that WntA plays in orchestrating pattern variation and underscores how rapidly the underlying GRNs shaping wing patterns can diverge. (A) WntA seems to be broadly involved in establishing pattern boundaries across the wing and causes shifts in the ultimate scale cell fate in every possible direction. (B) WntA KO effects are highly context specific within H. erato color pattern races. While the same proximal elements are affected, the color switch varies depending on genetic background. (C) The underlying GRNs controlling mimicry have diverged substantially between the MCS/ESS lineages. This is exemplified by the co-mimetic species H. pachinus and H. hewitsoni. WntA KOs in H. hewitsoni affect pattern across most of the forewing, except the most distal regions. By contrast, H. pachinus KO effects extend more distally, and forewing band regions remain unchanged. For each Heliconius wild type and mutant image, we show a drawing of the wings, where green represents the wing color pattern domains that are changed by knocking out the WntA gene. In (A), light blue represents changes in distal patterns and orange represents an extension of black patterns.

The other two major patterning genes identified in Heliconius, optix and aristaless1, are both transcription factors that are expressed later in wing scale cell development in response to signals laid down either directly or indirectly by WntA and cortex. The transcription factor optix is expressed at about 72 h after pupation in Type III scale cells that are fated to become red. Functional knockout of this gene in H. erato results in the loss of red color patterns and is associated with changes in scale cell morphology (Zhang et al., 2017b). Consistent with a role in switching between ommochrome and melanin synthesis, an RNAseq experiment comparing optix mutant and wild type wings in two species of butterflies, Vanessa cardui and Junonia coenia, showed that optix, either directly or indirectly, regulates the expression of many known pigmentation genes involved in ommochrome and melanin synthesis (Zhang et al., 2017b). Lastly, the homeodomain transcription factor aristaless1 causes a change in Type I scale cells and promotes the development of structural white scales by repressing the synthesis or uptake of yellow 3-OHK (Westerman et al., 2018). The aristaless gene is duplicated in Lepidoptera (Martin and Reed, 2010) and in divergent color pattern races of H. cydno there is a strong peak of differentiation associated with the presence/absence of white coloration in a putative cis-regulatory region midway between the two aristaless genes, aristaless1 and aristaless2. Both genes are differentially expressed in the forewings of white and yellow H. cydno morphs, with aristaless1 strongly upregulated late in wing development of the white morph. Consistent with a role in patterning, functional knockout of aristaless1 in white morphs results in the transformation of white scale cells into yellow scale cells (Figure 3 and Westerman et al., 2018).

Chance and Contingency on the Heliconius Wing

In Heliconius, there are many replicated cases of mimetic convergence between species within the ESS and MCS clades (see Figure 1 for some examples), providing a powerful context to understand how patterning networks co-evolve to produce variation. One unresolved question is the extent to which similar or distinct mechanisms have generated convergent phenotypes – in other words, are phenotypic repeats based on parallel reiterations of the same process, or do distinct developmental pathways lead to the same outcome? One of the more intriguing conclusions to emerge from comparative studies between species in the two Heliconius clades is how labile the GRNs responsible for patterning a Heliconius wing are (Concha et al., 2019). WntA knockouts in co-mimetic species of Heliconius result in very different phenotypes, demonstrating that the underlying GRNs controlling mimicry have diverged substantially between the two major clades (Concha et al., 2019; Figure 5C). Therefore, while WntA has been the repeated target of selection, its role in establishing pattern elements varies greatly between co-mimics. This fact is perhaps best demonstrated by comparing the KO phenotype between the two co-mimics H. hewitsoni (ESS clade) and H. pachinus (MCS clade). The two co-occur in forests of Panama and Costa Rica and have independently evolved identical wing patterns. However, the WntA KO forewing phenotypes are very different from each other, with WntA affecting pattern from the base of the forewing to nearly the tip in H. hewitsoni, whereas, the KO phenotype in H. pachinus causes a switch from melanic (Type III) to yellow Type I scale cells much more distally. In addition, knocking out WntA in H. pachinus does not affect the black melanic patch across the midsection of the forewing (Figure 5C).

The WntA knockout effects across the radiation highlight rapid divergence in the developmental landscape responsible for pattering a Heliconius wing; including changes in cis-regulatory variation around the locus, divergence in the interpretation of WntA signals by downstream developmental pathways, and variable epistatic interactions between WntA and other patterning loci (Concha et al., 2019). This finding is consistent with a recent study showing considerable differences in the expression patterns of transcription factors and signaling molecules between the two co-mimetic species H. e. demophoon and H. m. rosina (Hanly et al., 2019). Together these findings reinforce a growing appreciation that gene regulatory networks can diverge rapidly (True and Haag, 2001) and create a complex set of historical contingences that we are only beginning to appreciate (see Kittelmann et al., 2018; Sorrells et al., 2018; Ellison and Bachtrog, 2019). On a Heliconius wing, this divergence creates extraordinary flexibility for natural selection to drive the rapid evolution of highly convergent wing patterns in unpredictable ways.

Current Challenges

While the functional genomic work on major wing patterning genes is beginning to provide a picture of the molecular events underlying the process of scale cell differentiation, there are many unresolved questions regarding the details of this process.

How Are Patterning Signals Propagated to Alter Scale Cell Fate?

The major wing patterning loci in Heliconius are “evolutionary hotspots” and have been repeatedly targeted at both micro and macro-evolutionary scales to produce pattern variation. However, wing pattern development is likely to involve the coordinated expression of hundreds of genes. At the most basic level, we currently have little understanding of how WntA or cortex signals are propagated and interpreted by developing scale cells. Indeed, cortex does not conform to the typical paradigm of a wing patterning locus (sensu; Martin and Courtier-Orgogozo, 2017). Unlike the other wing patterning genes (Table 1), cortex is neither a conserved transcription factor nor a component of a known cell signaling pathway. Rather, the gene shows the closest homology to the cdc20/fizzy family of cell-cycle regulators (Nadeau et al., 2016). In Drosophila, Fizzy proteins are known to accumulate at kinetochores and centrosomes and regulate APC/C activity through the degradation of cyclins, which has been shown to prevent mitotic entry (Raff et al., 2002; Meghini et al., 2016). Specifically, expression of fizzy-related in Drosophila follicle cells induces cell endocycling and results in polyploidy, and cortex has been shown to interact with the APC/C and downregulate cyclins during female meiosis (Schaeffer et al., 2004; Shcherbata et al., 2004; Pesin and Orr-Weaver, 2007; Swan and Schüpbach, 2007). The Heliconius ortholog also localizes to the nucleus in developing pupal wings, suggesting it may also be interacting with the APC/C in developing butterfly scales (Livraghi et al., 2020). Other work has shown a correlation between scale cell ploidy and pigmentation state where increased ploidy coincides with darker scales (Cho and Nijhout, 2013; Iwata and Otaki, 2016) and cortex may thus be modulating pigmentation state by varying endocycle rates across the wing. Alternatively, cortex may be acting by modulating the developmental rate of scale cells during wing development. Melanic scales develop slower than the other scale cell types with melanin being synthesized and deposited later in wing development (Gilbert et al., 1988; Reed et al., 2008; Aymone et al., 2013). By controlling the faster or slower development of the scales, cortex may be modulating the final pigmentation state of a scale cell. This hypothesis of “developmental heterochrony” to explain the mechanism of butterfly wing color patterning was proposed 20 years ago (Koch et al., 2000) without gaining much supporting evidence over the years.

The mechanism of how the locus containing cortex affects pattern variation is probably complex. Although cortex is clearly involved in determining scale cell identity, there are other genes that are physically linked to cortex that might also play a role in patterning. Most notably, the genes domeless and washout show pattern specific expression in both H. melpomene and H. numata (Nadeau et al., 2016; Saenko et al., 2019). Domeless encodes a transmembrane receptor of the JAK/STAT signaling pathway (Brown et al., 2001) while washout is involved in actin polymerization through the activation of the Arp2/3 complex (Liu et al., 2009; Verboon et al., 2015), raising the interesting possibility that it may be involved in the formation of adult scales through the control of actin bundle elongation (Dinwiddie et al., 2014). It is therefore possible that multiple linked genes around cortex contribute to wing patterning evolution in Heliconius. Theory predicts that in some cases the linking of such genes might be favored by evolution (Charlesworth and Charlesworth, 1976), leading to the emergence of “supergenes”: tightly linked functional genetic elements that allow a switch between phenotypes to be maintained in a stable local polymorphism (Joron et al., 2011; Le Poul et al., 2014). Indeed, there is some indication that such a system is acting to control polymorphism in H. numata, where all wing pattern variation is controlled by inversions surrounding cortex, domeless and washout (Joron et al., 2011; Nadeau et al., 2016; Saenko et al., 2019). However, we currently have no functional validation for other genes at this locus, nor evidence for interactions between them.

WntA, in contrast, is a classic cell signaling molecule belonging to the family of Wnt signaling ligands. However, we also currently lack an understanding of how WntA signals are transduced into a coordinated scale differentiation program. Wnt signaling may act through the canonical pathway, by which β-catenin is permitted to translocate to the nucleus and act as a transcriptional cofactor (reviewed in Collu et al., 2014), or alternatively it may signal through a non-canonical modality such as the planar cell polarity pathway (PCP) (Lawrence and Casal, 2018), or by a third mechanism such as Wnt-calcium signaling (Agrawal and Hasan, 2015). Understanding the mechanism by which WntA causes patterned scale cell differentiation will be an important first step in unraveling the larger gene regulatory networks that pattern the wings.

How Does New Regulatory Variation Evolve?

One idea to emerge from developmental biology is that loci that have become the nexus of morphological change, so called hotspot genes, are often associated with complex cis-regulatory variation (Stern and Orgogozo, 2009). This is certainly true of the major patterning loci identified in Heliconius. Expression, association, and CRISPR studies have all revealed that cis-regulatory variation is important in driving wing patterning differences. Nonetheless, we still don’t understand how new regulatory regions evolve and what role, if any, ancient regulatory architecture plays in wing pattern variation. Association studies in Heliconius suggest that there is a high degree of modularity with regards to the cis-regulatory elements controlling color pattern variation (Supple et al., 2013; Wallbank et al., 2016; Enciso-Romero et al., 2017; Van Belleghem et al., 2017). For example, different red color pattern elements in H. melpomene and H. erato are associated with discrete genomic intervals in non-coding regions around optix (Wallbank et al., 2016; Van Belleghem et al., 2017). A similar modular architecture is also observed for variation in forewing band shape (WntA) and the presence or absence of a yellow hindwing bar (cortex) within H. erato (Van Belleghem et al., 2017). The idea is that new wing pattern phenotypes can evolve rapidly by simply reshuffling these discrete regulatory regions. In the MCS lineage, for example, the phylogenetic patterns are consistent with at least four instances of between species hybridization that allowed these elements to cross species boundaries, caused pattern change, and probably promoted speciation in the group (Wallbank et al., 2016).

A model where “enhancer shuffling” drives morphological differentiation, however, was not strongly supported by recent work that has more directly probed the relationship between regulatory and phenotypic variation (Lewis et al., 2019). Here, the authors identified and excised a number of putative cis-regulatory regions in genomic intervals associated with red pattern variation in H. erato using a combination of epigenetic profiling and genotype × phenotype analysis. Surprisingly, they found that deletions of a single cis-regulatory element were highly pleiotropic and affected red pattern elements in both the forewing and hindwing. The level of pleiotropy observed suggests that more constraints are acting at color pattern loci than previously hypothesized and this is difficult to reconcile under a model where a single regulatory region affects a single pattern feature. A possible resolution to this apparent conflict between the association and functional data may be explained by the idiosyncratic properties of phenotypes segregating in Heliconius hybrid zones, whereby a single locus can act in a modular fashion in localized population structures; however, because of an ever shifting trans-regulatory landscape, excising a specific enhancer after patterns have become established and genetic backgrounds have diverged is not sufficient to induce the predicted phenotype (see Lewis and Van Belleghem, unpublished, for a discussion of this topic). At this point, we clearly need more functional studies that assess the evidence for and against modularity in driving the evolution of Heliconius wing patterns.

How Labile Is Pattern Variation on the Wing Across the Lepidoptera?

In Heliconius, much of the pattern variation across the genus is driven by variation at homologous genes. However, even in this relatively recent radiation, developmental systems drift appears to be shifting the wiring of the underlying GRNs (Concha et al., 2019). Expanding functional studies across the lepidoptera will be necessary to start addressing questions relating to the evolutionary constraints acting across GRNs as a whole. Moreover, while mapping studies across Heliconius consistently implicate four major patterning genes, variation in other loci across the genome also have quantitative effects on pattern variation (Papa et al., 2013). Indeed, more recent mapping experiments show that structural coloration in Heliconius is not achieved using any of the known patterning loci and that different loci are used to produce structural colors in the two distantly related co-mimics H. erato and H. melpomene (Parnell et al., 2018; Brien et al., 2019). Also, while cortex has been repeatedly mapped in within-species pattern variation, recent studies have begun to identify other genes, not previously identified in Heliconius, as responsible for pattern variation (Table 1). In this light, it will be very interesting to determine if these other conserved components of GRNs underlying wing patterning in other lepidoptera contribute to pattern variation in Heliconius.

While individual components have been identified through evolutionary approaches, very few regulatory connections have been resolved, and we still don’t know the broad applicability of discoveries made in Heliconius to other butterflies or moths. From an evolutionary perspective, Heliconius wings represent a highly derived type of a more ancestral patterning system called the Nymphalid Ground Plan (NGP) (Nijhout and Wray, 1988). In this sense, Heliconius wings may not represent the best system with which to address questions relating to the conservation of GRNs in other lepidoptera, as it is difficult to draw homologies between color pattern elements at larger evolutionary scales. For example, while WntA has been implicated in color pattern evolution across both Heliconius and other nymphalids (Martin et al., 2012; Gallant et al., 2014; Van Belleghem et al., 2017), the pattern elements it affects vary substantially. In species displaying NGP-like patterns, WntA functions as an organizing factor of Central and External Symmetry Systems in at least three distantly related nymphalids (Mazo-Vargas et al., 2017), suggesting a conserved role across large evolutionary scales. On the other hand, it is noteworthy that all of the four major patterning loci discovered in Heliconius affect pattern in other butterflies and moths. In this light, the apparent flexibility in the GRNs that we are observing in Heliconius may reflect the nature of wing patterning networks generally. Thus, while a particular set of conserved loci might have a function in generating pattern variation, how they generate that variation may differ substantially among lineages. For example, the eyespot patterns that are found in many nymphalid butterflies appear to have evolved once; but the gene regulatory networks underlying these patterns have diverged significantly over the past 65 million years such that the eyespots of different species are controlled by different sets of genes (Oliver et al., 2012).

Linking Genotype to Phenotype

In multicellular organisms, phenotype is ultimately realized in the patterned differentiation of cells into specific functional types. In this respect, butterfly wings have a unique combination of features which make them an ideal system to uncover how changes in a DNA sequence are translated and interpreted through development to generate variation: a relatively simple developmental system, large phenotypic diversity, and experimental tractability.

There is clearly much to be discovered, but the last few years have witnessed remarkable progress in our understanding of pattern formation on butterfly wings. Research priorities are clearly laid out and technical improvements, including techniques for characterizing the developmental landscape of protein–DNA interactions (Park, 2009; Rao et al., 2014; Buenrostro et al., 2015; Lewis et al., 2016, 2019) and the expression profiles of individual cells (Kulkarni et al., 2019), coupled with more efficient functional manipulation provide us with the necessary instruments to further uncover the gene regulatory networks underlying wing pattern variation. In this respect, one of the key technical limitations facing Heliconius is the lack of transgenic tools. Our current use of CRISPR as a functional tool has been limited to the analysis of mosaic loss-of-function mutants, where only a subset of cells carry the desired mutation on the adult wing. While this can be a strength, allowing the analysis of the effects of highly pleiotropic genes and permitting the observation of effects at clone boundaries, in many cases this makes it difficult to determine the extent to which a given induced mutation affects pattern across the wing. These shortcomings could be bypassed by developing technologies that allow researchers to routinely produce gain-of-function mutations; as well as those that would allow better spatiotemporal control over genetic manipulations. For example, by placing reporter constructs under specific inducible elements, the ability of a putative cis-regulatory element to drive expression in the tissue of interest could be directly visualized (Tong et al., 2014). Nonetheless, if we stay on pace we will have come a long way toward fulfilling the prophecy made by Bates (1864) and “the study of butterflies—creatures selected as the types of airiness and frivolity—instead of being despised, will someday be valued as one of the most important branches of Biological science.”

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

WM was supported by the Smithsonian Tropical Research Institute and the Smithsonian Scholarly Studies Program. LL was supported by BBSRC BB/R007500/1. JH was supported by the NSF award IOS-1656553. CC was supported by the Smithsonian Institute Biodiversity Genomics Fellowship and NSF IOS-1656389.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank the Butterfly Ecology and Evolution Research group at STRI, as well as the Butterfly Genetics group at the University of Cambridge, the Martin Lab at GWU and the Papa Lab at UPR-RP for constructive conversations that helped to formulate this review. In addition, we thank Marcus Kronforst, Bob Reed, and Steven Van Belleghem for providing images of butterflies, and two reviewers for their helpful feedback on this manuscript.

References

Agrawal, T., and Hasan, G. (2015). Maturation of a central brain flight circuit in Drosophila requires Fz2/Ca2+ signaling. eLife 4:e07046. doi: 10.7554/eLife.07046.001