A genetic cause of male mate preference.

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

Sex-linked gene traffic underlies the acquisition of sexually dimorphic UV color vision in Heliconius butterflies.

The acquisition of novel sexually dimorphic traits poses an evolutionary puzzle: How do new traits arise and become sex-limited? Recently acquired color vision, sexually dimorphic in animals like primates and butterflies, presents a compelling model for understanding how traits become sex-biased. For example, some Heliconius butterflies uniquely possess UV (ultraviolet) color vision, which correlates with the expression of two differentially tuned UV-sensitive rhodopsins, UVRh1 and UVRh2. To discover how such traits become sexually dimorphic, we studied Heliconius charithonia, which exhibits female-specific UVRh1 expression. We demonstrate that females, but not males, discriminate different UV wavelengths. Through whole-genome shotgun sequencing and assembly of the H. charithonia genome, we discovered that UVRh1 is present on the W chromosome, making it obligately female-specific. By knocking out UVRh1, we show that UVRh1 protein expression is absent in mutant female eye tissue, as in wild-type male eyes. A PCR survey of UVRh1 sex-linkage across the genus shows that species with female-specific UVRh1 expression lack UVRh1 gDNA in males. Thus, acquisition of sex linkage is sufficient to achieve female-specific expression of UVRh1, though this does not preclude other mechanisms, like cis-regulatory evolution from also contributing. Moreover, both this event, and mutations leading to differential UV opsin sensitivity, occurred early in the history of Heliconius. These results suggest a path for acquiring sexual dimorphism distinct from existing mechanistic models. We propose a model where gene traffic to heterosomes (the W or the Y) genetically partitions a trait by sex before a phenotype shifts (spectral tuning of UV sensitivity).

Insect opsins and evo-devo: what have we learned in 25 years?

The visual pigments known as opsins are the primary molecular basis for colour vision in animals. Insects are among the most diverse of animal groups and their visual systems reflect a variety of life histories. The study of insect opsins in the fruit fly Drosophila melanogaster has led to major advances in the fields of neuroscience, development and evolution. In the last 25 years, research in D. melanogaster has improved our understanding of opsin genotype-phenotype relationships while comparative work in other insects has expanded our understanding of the evolution of insect eyes via gene duplication, coexpression and homologue switching. Even so, until recently, technology and sampling have limited our understanding of the fundamental mechanisms that evolution uses to shape the diversity of insect eyes. With the advent of genome editing and in vitro expression assays, the study of insect opsins is poised to reveal new frontiers in evolutionary biology, visual neuroscience, and animal behaviour. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

Multiple Mechanisms of Photoreceptor Spectral Tuning in Heliconius Butterflies.

The evolution of color vision is often studied through the lens of receptor gain relative to an ancestor with fewer spectral classes of photoreceptor. For instance, in Heliconius butterflies, a genus-specific UVRh opsin duplication led to the evolution of UV color discrimination in Heliconius erato females, a rare trait among butterflies. However, color vision evolution is not well understood in the context of loss. In Heliconius melpomene and Heliconius ismenius lineages, the UV2 receptor subtype has been lost, which limits female color vision in shorter wavelengths. Here, we compare the visual systems of butterflies that have either retained or lost the UV2 photoreceptor using intracellular recordings, ATAC-seq, and antibody staining. We identify several ways these butterflies modulate their color vision. In H. melpomene, chromatin reorganization has downregulated an otherwise intact UVRh2 gene, whereas in H. ismenius, pseudogenization has led to the truncation of UVRh2. In species that lack the UV2 receptor, the peak sensitivity of the remaining UV1 photoreceptor cell is shifted to longer wavelengths. Across Heliconius, we identify the widespread use of filtering pigments and co-expression of two opsins in the same photoreceptor cells. Multiple mechanisms of spectral tuning, including the molecular evolution of blue opsins, have led to the divergence of receptor sensitivities between species. The diversity of photoreceptor and ommatidial subtypes between species suggests that Heliconius visual systems are under varying selection pressures for color discrimination. Modulating the wavelengths of peak sensitivities of both the blue- and remaining UV-sensitive photoreceptor cells suggests that Heliconius species may have compensated for UV receptor loss.

Air temperature drives the evolution of mid-infrared optical properties of butterfly wings.

This study uncovers a correlation between the mid-infrared emissivity of butterfly wings and the average air temperature of their habitats across the world. Butterflies from cooler climates have a lower mid-infrared emissivity, which limits heat losses to surroundings, and butterflies from warmer climates have a higher mid-infrared emissivity, which enhances radiative cooling. The mid-infrared emissivity showed no correlation with other investigated climatic factors. Phylogenetic independent contrasts analysis indicates the microstructures of butterfly wings may have evolved in part to regulate mid-infrared emissivity as an adaptation to climate, rather than as phylogenetic inertia. Our findings offer new insights into the role of microstructures in thermoregulation and suggest both evolutionary and physical constraints to butterflies' abilities to adapt to climate change.

True UV color vision in a female butterfly with two UV opsins.

In true color vision, animals discriminate between light wavelengths, regardless of intensity, using at least two photoreceptors with different spectral sensitivity peaks. Heliconius butterflies have duplicate UV opsin genes, which encode ultraviolet and violet photoreceptors, respectively. In Heliconius erato, only females express the ultraviolet photoreceptor, suggesting females (but not males) can discriminate between UV wavelengths. We tested the ability of H. erato, and two species lacking the violet receptor, Heliconius melpomene and Eueides isabella, to discriminate between 380 and 390 nm, and between 400 and 436 nm, after being trained to associate each stimulus with a sugar reward. We found that only H. erato females have color vision in the UV range. Across species, both sexes show color vision in the blue range. Models of H. erato color vision suggest that females have an advantage over males in discriminating the inner UV-yellow corollas of Psiguria flowers from their outer orange petals. Moreover, previous models ( McCulloch et al., 2017) suggested that H. erato males have an advantage over females in discriminating Heliconius 3-hydroxykynurenine (3-OHK) yellow wing coloration from non-3-OHK yellow wing coloration found in other heliconiines. These results provide some of the first behavioral evidence for female H. erato UV color discrimination in the context of foraging, lending support to the hypothesis ( Briscoe et al., 2010) that the duplicated UV opsin genes function together in UV color vision. Taken together, the sexually dimorphic visual system of H. erato appears to have been shaped by both sexual selection and sex-specific natural selection.

A two-step method for identifying photopigment opsin and rhodopsin gene sequences underlying human color vision phenotypes.

Purpose: To present a detailed, reliable long range-PCR and sequencing (LR-PCR-Seq) procedure to identify human opsin gene sequences for variations in the long wavelength-sensitive (OPN1LW), medium wavelength-sensitive (OPN1MW), short wavelength-sensitive (OPN1SW), and rhodopsin (RHO) genes. Methods: Color vision was assessed for nine subjects using the Farnsworth-Munsell 100 hue test, Ishihara pseudoisochromatic plates, and the Rabin cone-contrast threshold procedure (ColorDX, Konan Medical). The color vision phenotypes were normal trichromacy (n = 3), potential tetrachromacy (n = 3), dichromacy (n = 2), and unexplained low color vision (n = 1). DNA was isolated from blood or saliva and LR-PCR amplified into individual products: OPN1LW (4,045 bp), OPN1MW (4,045 bp), OPN1SW (3,326 bp), and RHO (6,715 bp). Each product was sequenced using specific internal primer sets. Analysis was performed with Mutation Surveyor software. Results: The LR-PCR-Seq technique identified known single nucleotide polymorphisms (SNPs) in OPN1LW and OPN1MW gene codons (180, 230, 233, 277, and 285), as well as those for lesser studied codons (174, 178, 236, 274, 279, 298 and 309) in the OPN1LW and OPN1MW genes. Additionally, six SNP variants in the OPN1MW and OPN1LW genes not previously reported in the NCBI dbSNP database were identified. An unreported poly-T region within intron 5(c.36+126) of the rhodopsin gene was also found, and analysis showed it to be highly conserved in mammalian species. Conclusions: This LR-PCR-Seq procedure (single PCR reaction per gene followed by sequencing) can identify exonic and intronic SNP variants in OPN1LW, OPN1MW, OPN1SW, and rhodopsin genes. There is no need for restriction enzyme digestion or multiple PCR steps that can introduce errors. Future studies will combine the LR-PCR-Seq with perceptual behavior measures, allowing for accurate correlations between opsin genotypes, retinal photopigment phenotypes, and color perception behaviors.

Disentangling Population History and Character Evolution among Hybridizing Lineages.

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

Infrared optical and thermal properties of microstructures in butterfly wings.

While surface microstructures of butterfly wings have been extensively studied for their structural coloration or optical properties within the visible spectrum, their properties in infrared wavelengths with potential ties to thermoregulation are relatively unknown. The midinfrared wavelengths of 7.5 to 14 µm are particularly important for radiative heat transfer in the ambient environment, because of the overlap with the atmospheric transmission window. For instance, a high midinfrared emissivity can facilitate surface cooling, whereas a low midinfrared emissivity can minimize heat loss to surroundings. Here we find that the midinfrared emissivity of butterfly wings from warmer climates such as Archaeoprepona demophoon (Oaxaca, Mexico) and Heliconius sara (Pichincha, Ecuador) is up to 2 times higher than that of butterfly wings from cooler climates such as Celastrina echo (Colorado) and Limenitis arthemis (Florida), using Fourier-transform infrared (FTIR) spectroscopy and infrared thermography. Our optical computations using a unit cell approach reproduce the spectroscopy data and explain how periodic microstructures play a critical role in the midinfrared. The emissivity spectrum governs the temperature of butterfly wings, and we demonstrate that C. echo wings heat up to 8 °C more than A. demophoon wings under the same sunlight in the clear sky of Irvine, CA. Furthermore, our thermal computations show that butterfly wings in their respective habitats can maintain a moderate temperature range through a balance of solar absorption and infrared emission. These findings suggest that the surface microstructures of butterfly wings potentially contribute to thermoregulation and provide an insight into butterflies' survival.

Drift and Directional Selection Are the Evolutionary Forces Driving Gene Expression Divergence in Eye and Brain Tissue of Heliconius Butterflies.

Investigating gene expression evolution over micro- and macroevolutionary timescales will expand our understanding of the role of gene expression in adaptation and speciation. In this study, we characterized the evolutionary forces acting on gene expression levels in eye and brain tissue of five Heliconius butterflies with divergence times of ∼5-12 MYA. We developed and applied Brownian motion (BM) and Ornstein-Uhlenbeck (OU) models to identify genes whose expression levels are evolving through drift, stabilizing selection, or a lineage-specific shift. We found that 81% of the genes evolve under genetic drift. When testing for branch-specific shifts in gene expression, we detected 368 (16%) shift events. Genes showing a shift toward upregulation have significantly lower gene expression variance than those genes showing a shift leading toward downregulation. We hypothesize that directional selection is acting in shifts causing upregulation, since transcription is costly. We further uncovered through simulations that parameter estimation of OU models is biased when using small phylogenies and only becomes reliable with phylogenies having ≥ 50 taxa. Therefore, we developed a new statistical test based on BM to identify highly conserved genes (i.e., evolving under strong stabilizing selection), which comprised 3% of the orthoclusters. In conclusion, we found that drift is the dominant evolutionary force driving gene expression evolution in eye and brain tissue in Heliconius Nevertheless, the higher proportion of genes evolving under directional than under stabilizing selection might reflect species-specific selective pressures on vision and the brain that are necessary to fulfill species-specific requirements.

Evolution of Phototransduction Genes in Lepidoptera.

Vision is underpinned by phototransduction, a signaling cascade that converts light energy into an electrical signal. Among insects, phototransduction is best understood in Drosophila melanogaster. Comparison of D. melanogaster against three insect species found several phototransduction gene gains and losses, however, lepidopterans were not examined. Diurnal butterflies and nocturnal moths occupy different light environments and have distinct eye morphologies, which might impact the expression of their phototransduction genes. Here we investigated: 1) how phototransduction genes vary in gene gain or loss between D. melanogaster and Lepidoptera, and 2) variations in phototransduction genes between moths and butterflies. To test our prediction of phototransduction differences due to distinct visual ecologies, we used insect reference genomes, phylogenetics, and moth and butterfly head RNA-Seq and transcriptome data. As expected, most phototransduction genes were conserved between D. melanogaster and Lepidoptera, with some exceptions. Notably, we found two lepidopteran opsins lacking a D. melanogaster ortholog. Using antibodies we found that one of these opsins, a candidate retinochrome, which we refer to as unclassified opsin (UnRh), is expressed in the crystalline cone cells and the pigment cells of the butterfly, Heliconius melpomene. Our results also show that butterflies express similar amounts of trp and trpl channel mRNAs, whereas moths express ∼50× less trp, a potential adaptation to darkness. Our findings suggest that while many single-copy D. melanogaster phototransduction genes are conserved in lepidopterans, phototransduction gene expression differences exist between moths and butterflies that may be linked to their visual light environment.

Experimental field tests of Batesian mimicry in the swallowtail butterfly Papilio polytes.

The swallowtail butterfly Papilio polytes is known for its striking resemblance in wing pattern to the toxic butterfly Pachliopta aristolochiae and is a focal system for the study of mimicry evolution. Papilio polytes females are polymorphic in wing pattern, with mimetic and nonmimetic forms, while males are monomorphic and nonmimetic. Past work invokes selection for mimicry as the driving force behind wing pattern evolution in P. polytes. However, the mimetic relationship between P. polytes and P. aristolochiae is not well understood. In order to test the mimicry hypothesis, we constructed paper replicas of mimetic and nonmimetic P. polytes and P. aristolochiae, placed them in their natural habitat, and measured bird predation on replicas. In initial trials with stationary replicas and plasticine bodies, overall predation was low and we found no differences in predation between replica types. In later trials with replicas mounted on springs and with live mealworms standing in for the butterfly's body, we found less predation on mimetic P. polytes replicas compared to nonmimetic P. polytes replicas, consistent with the predator avoidance benefits of mimicry. While our results are mixed, they generally lend support to the mimicry hypothesis as well as the idea that behavioral differences between the sexes contributed to the evolution of sexually dimorphic mimicry.

Evolution of Sex-Biased Gene Expression and Dosage Compensation in the Eye and Brain of Heliconius Butterflies.

Differences in behavior and life history traits between females and males are the basis of divergent selective pressures between sexes. It has been suggested that a way for the two sexes to deal with different life history requirements is through sex-biased gene expression. In this study, we performed a comparative sex-biased gene expression analysis of the combined eye and brain transcriptome from five Heliconius species, H. charithonia, H. sara, H. erato, H. melpomene and H. doris, representing five of the main clades from the Heliconius phylogeny. We found that the degree of sexual dimorphism in gene expression is not conserved across Heliconius. Most of the sex-biased genes identified in each species are not sex-biased in any other, suggesting that sexual selection might have driven sexually dimorphic gene expression. Only three genes shared sex-biased expression across multiple species: ultraviolet opsin UVRh1 and orthologs of Drosophila Krüppel-homolog 1 and CG9492. We also observed that in some species female-biased genes have higher evolutionary rates, but in others, male-biased genes show the fastest rates when compared with unbiased genes, suggesting that selective forces driving sex-biased gene evolution in Heliconius act in a sex- and species-specific manner. Furthermore, we found dosage compensation in all the Heliconius tested, providing additional evidence for the conservation of dosage compensation across Lepidoptera. Finally, sex-biased genes are significantly enriched on the Z, a pattern that could be a result of sexually antagonistic selection.

Frequency dependence shapes the adaptive landscape of imperfect Batesian mimicry.

Despite more than a century of biological research on the evolution and maintenance of mimetic signals, the relative frequencies of models and mimics necessary to establish and maintain Batesian mimicry in natural populations remain understudied. Here we investigate the frequency-dependent dynamics of imperfect Batesian mimicry, using predation experiments involving artificial butterfly models. We use two geographically distinct populations of Adelpha butterflies that vary in their relative frequencies of a putatively defended model (Adelpha iphiclus) and Batesian mimic (Adelpha serpa). We found that in Costa Rica, where both species share similar abundances, Batesian mimicry breaks down, and predators more readily attack artificial butterfly models of the presumed mimic, A. serpa By contrast, in Ecuador, where A. iphiclus (model) is significantly more abundant than A. serpa (mimic), both species are equally protected from predation. Our results provide compelling experimental evidence that imperfect Batesian mimicry is frequency-dependent on the relative abundance of models and mimics in natural populations, and contribute to the growing body of evidence that complex dynamics, such as seasonality or the availability of alternative prey, influence the evolution of mimetic traits.

Longwing (Heliconius) butterflies combine a restricted set of pigmentary and structural coloration mechanisms.

BACKGROUND: Longwing butterflies, Heliconius sp., also called heliconians, are striking examples of diversity and mimicry in butterflies. Heliconians feature strongly colored patterns on their wings, arising from wing scales colored by pigments and/or nanostructures, which serve as an aposematic signal. RESULTS: Here, we investigate the coloration mechanisms among several species of Heliconius by applying scanning electron microscopy, (micro)spectrophotometry, and imaging scatterometry. We identify seven kinds of colored scales within Heliconius whose coloration is derived from pigments, nanostructures or both. In yellow-, orange- and red-colored wing patches, both cover and ground scales contain wavelength-selective absorbing pigments, 3-OH-kynurenine, xanthommatin and/or dihydroxanthommatin. In blue wing patches, the cover scales are blue either due to interference of light in the thin-film lower lamina (e.g., H. doris) or in the multilayered lamellae in the scale ridges (so-called ridge reflectors, e.g., H. sara and H. erato); the underlying ground scales are black. In the white wing patches, both cover and ground scales are blue due to their thin-film lower lamina, but because they are stacked upon each other and at the wing substrate, a faint bluish to white color results. Lastly, green wing patches (H. doris) have cover scales with blue-reflecting thin films and short-wavelength absorbing 3-OH-kynurenine, together causing a green color. CONCLUSIONS: The pigmentary and structural traits are discussed in relation to their phylogenetic distribution and the evolution of vision in this highly interesting clade of butterflies.

Copy Number Variation and Expression Analysis Reveals a Nonorthologous Pinta Gene Family Member Involved in Butterfly Vision.

Vertebrate (cellular retinaldehyde-binding protein) and Drosophila (prolonged depolarization afterpotential is not apparent [PINTA]) proteins with a CRAL-TRIO domain transport retinal-based chromophores that bind to opsin proteins and are necessary for phototransduction. The CRAL-TRIO domain gene family is composed of genes that encode proteins with a common N-terminal structural domain. Although there is an expansion of this gene family in Lepidoptera, there is no lepidopteran ortholog of pinta. Further, the function of these genes in lepidopterans has not yet been established. Here, we explored the molecular evolution and expression of CRAL-TRIO domain genes in the butterfly Heliconius melpomene in order to identify a member of this gene family as a candidate chromophore transporter. We generated and searched a four tissue transcriptome and searched a reference genome for CRAL-TRIO domain genes. We expanded an insect CRAL-TRIO domain gene phylogeny to include H. melpomene and used 18 genomes from 4 subspecies to assess copy number variation. A transcriptome-wide differential expression analysis comparing four tissue types identified a CRAL-TRIO domain gene, Hme CTD31, upregulated in heads suggesting a potential role in vision for this CRAL-TRIO domain gene. RT-PCR and immunohistochemistry confirmed that Hme CTD31 and its protein product are expressed in the retina, specifically in primary and secondary pigment cells and in tracheal cells. Sequencing of eye protein extracts that fluoresce in the ultraviolet identified Hme CTD31 as a possible chromophore binding protein. Although we found several recent duplications and numerous copy number variants in CRAL-TRIO domain genes, we identified a single copy pinta paralog that likely binds the chromophore in butterflies.

Sexual Dimorphism and Retinal Mosaic Diversification following the Evolution of a Violet Receptor in Butterflies.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Ultraviolet and yellow reflectance but not fluorescence is important for visual discrimination of conspecifics by Heliconius erato.

Toxic Heliconius butterflies have yellow hindwing bars that - unlike those of their closest relatives - reflect ultraviolet (UV) and long wavelength light, and also fluoresce. The pigment in the yellow scales is 3-hydroxy-dl-kynurenine (3-OHK), which is found in the hair and scales of a variety of animals. In other butterflies like pierids with color schemes characterized by independent sources of variation in UV and human-visible yellow/orange, behavioral experiments have generally implicated the UV component as most relevant to mate choice. This has not been addressed in Heliconius butterflies, where variation exists in analogous color components, but moreover where fluorescence due to 3-OHK could also contribute to yellow wing coloration. In addition, the potential cost due to predator visibility is largely unknown for the analogous well-studied pierid butterfly species. In field studies with butterfly paper models, we show that both UV and 3-OHK yellow act as signals for H. erato when compared with models lacking UV or resembling ancestral Eueides yellow, respectively, but attack rates by birds do not differ significantly between the models. Furthermore, measurement of the quantum yield and reflectance spectra of 3-OHK indicates that fluorescence does not contribute to the visual signal under broad-spectrum illumination. Our results suggest that the use of 3-OHK pigmentation instead of ancestral yellow was driven by sexual selection rather than predation.

Complex dynamics underlie the evolution of imperfect wing pattern convergence in butterflies.

Adaptive radiation is characterized by rapid diversification that is strongly associated with ecological specialization. However, understanding the evolutionary mechanisms fueling adaptive diversification requires a detailed knowledge of how natural selection acts at multiple life-history stages. Butterflies within the genus Adelpha represent one of the largest and most diverse butterfly lineages in the Neotropics. Although Adelpha species feed on an extraordinary diversity of larval hosts, convergent evolution is widespread in this group, suggesting that selection for mimicry may contribute to adaptive divergence among species. To investigate this hypothesis, we conducted predation studies in Costa Rica using artificial butterfly facsimiles. Specifically, we predicted that nontoxic, palatable Adelpha species that do not feed on host plants in the family Rubiaceae would benefit from sharing a locally convergent wing pattern with the presumably toxic Rubiaceae-feeding species via reduced predation. Contrary to expectations, we found that the presumed mimic was attacked significantly more than its locally convergent model at a frequency paralleling attack rates on both novel and palatable prey. Although these data reveal the first evidence for protection from avian predators by the supposed toxic, Rubiaceae-feeding Adelpha species, we conclude that imprecise mimetic patterns have high costs for Batesian mimics in the tropics.

Gene Duplication and Gene Expression Changes Play a Role in the Evolution of Candidate Pollen Feeding Genes in Heliconius Butterflies.

Heliconius possess a unique ability among butterflies to feed on pollen. Pollen feeding significantly extends their lifespan, and is thought to have been important to the diversification of the genus. We used RNA sequencing to examine feeding-related gene expression in the mouthparts of four species of Heliconius and one nonpollen feeding species, Eueides isabella We hypothesized that genes involved in morphology and protein metabolism might be upregulated in Heliconius because they have longer proboscides than Eueides, and because pollen contains more protein than nectar. Using de novo transcriptome assemblies, we tested these hypotheses by comparing gene expression in mouthparts against antennae and legs. We first looked for genes upregulated in mouthparts across all five species and discovered several hundred genes, many of which had functional annotations involving metabolism of proteins (cocoonase), lipids, and carbohydrates. We then looked specifically within Heliconius where we found eleven common upregulated genes with roles in morphology (CPR cuticle proteins), behavior (takeout-like), and metabolism (luciferase-like). Closer examination of these candidates revealed that cocoonase underwent several duplications along the lineage leading to heliconiine butterflies, including two Heliconius-specific duplications. Luciferase-like genes also underwent duplication within lepidopterans, and upregulation in Heliconius mouthparts. Reverse-transcription PCR confirmed that three cocoonases, a peptidase, and one luciferase-like gene are expressed in the proboscis with little to no expression in labial palps and salivary glands. Our results suggest pollen feeding, like other dietary specializations, was likely facilitated by adaptive expansions of preexisting genes-and that the butterfly proboscis is involved in digestive enzyme production.