Molecular evolution of a long wavelength-sensitive opsin in mimetic Heliconius butterflies (Lepidoptera: Nymphalidae)

A cure for the blues: opsin duplication and subfunctionalization for short-wavelength sensitivity in jewel beetles (Coleoptera: Buprestidae)

Background

Arthropods have received much attention as a model for studying opsin evolution in invertebrates. Yet, relatively few studies have investigated the diversity of opsin proteins that underlie spectral sensitivity of the visual pigments within the diverse beetles (Insecta: Coleoptera). Previous work has demonstrated that beetles appear to lack the short-wavelength-sensitive (SWS) opsin class that typically confers sensitivity to the “blue” region of the light spectrum. However, this is contrary to established physiological data in a number of Coleoptera. To explore potential adaptations at the molecular level that may compensate for the loss of the SWS opsin, we carried out an exploration of the opsin proteins within a group of beetles (Buprestidae) where short-wave sensitivity has been demonstrated. RNA-seq data were generated to identify opsin proteins from nine taxa comprising six buprestid species (including three male/female pairs) across four subfamilies. Structural analyses of recovered opsins were conducted and compared to opsin sequences in other insects across the main opsin classes—ultraviolet, short-wavelength, and long-wavelength.

Results

All nine buprestids were found to express two opsin copies in each of the ultraviolet and long-wavelength classes, contrary to the single copies recovered in all other molecular studies of adult beetle opsin expression. No SWS opsin class was recovered. Furthermore, the male Agrilus planipennis (emerald ash borer—EAB) expressed a third LWS opsin at low levels that is presumed to be a larval copy. Subsequent homology and structural analyses identified multiple amino acid substitutions in the UVS and LWS copies that could confer short-wavelength sensitivity.

Conclusions

This work is the first to compare expressed opsin genes against known electrophysiological data that demonstrate multiple peak sensitivities in Coleoptera. We report the first instance of opsin duplication in adult beetles, which occurs in both the UVS and LWS opsin classes. Through structural comparisons of known insect opsins, we suggest that opsin duplication and amino acid variation within the chromophore binding pocket explains sensitivity in the short-wavelength portion of the visible light spectrum in these species. These findings are the first to reveal molecular complexity of the color vision system within beetles.

Electronic supplementary material

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

Sexual dimorphism in the compound eye of Heliconius erato: a nymphalid butterfly with at least five spectral classes of photoreceptor

Most butterfly families expand the number of spectrally-distinct photoreceptors in their compound eye by opsin gene duplications together with lateral filter pigments, however most nymphalid genera have limited diversity, with only three or four spectral types of photoreceptor. Here we examine the spatial pattern of opsin expression and photoreceptor spectral sensitivities in Heliconius erato, a nymphalid with duplicate ultraviolet opsin genes, UVRh1 and UVRh2. We find that the H. erato compound eye is sexually dimorphic. Females express the two UV opsin proteins in separate photoreceptors, but males do not express UVRh1. Intracellular recordings confirmed that females have three short wavelength-sensitive photoreceptors (λmax=356 nm, ∼390 nm and 470 nm), while males have two (λmax=390 nm and ∼470 nm). We also found two long wavelength-sensitive photoreceptors (green, λmax ∼555 nm, and red, λmax ∼600 nm), which express the same LW opsin. The red cell's shifted sensitivity is probably due to perirhabdomal filtering pigments. Sexual dimorphism of the UV-absorbing rhodopsins may reflect the females' need to discriminate conspecifics from co-mimics. Red-green color vision may be used to detect differences in red coloration on Heliconius wings, or for host-plant identification. Among nymphalids so far investigated, only H. erato is known to possess five spectral classes of photoreceptor; sexual dimorphism of the eye via suppression of one class of opsin (here UVRh1 in males) has not—to our knowledge—been reported in any animal.

UV photoreceptors and UV-yellow wing pigments in Heliconius butterflies allow a color signal to serve both mimicry and intraspecific communication

Mimetic wing coloration evolves in butterflies in the context of predator confusion. Unless butterfly eyes have adaptations for discriminating mimetic color variation, mimicry also carries a risk of confusion for the butterflies themselves. Heliconius butterfly eyes, which express recently duplicated ultraviolet (UV) opsins, have such an adaptation. To examine bird and butterfly color vision as sources of selection on butterfly coloration, we studied yellow wing pigmentation in the tribe Heliconiini. We confirmed, using reflectance and mass spectrometry, that only Heliconius use 3-hydroxy-DL-kynurenine (3-OHK), which looks yellow to humans but reflects both UV- and long-wavelength light, whereas butterflies in related genera have chemically unknown yellow pigments mostly lacking UV reflectance. Modeling of these color signals reveals that the two UV photoreceptors of Heliconius are better suited to separating 3-OHK from non-3-OHK spectra compared with the photoreceptors of related genera or birds. The co-occurrence of potentially enhanced UV vision and a UV-reflecting yellow wing pigment could allow unpalatable Heliconius private intraspecific communication in the presence of mimics. Our results are the best available evidence for the correlated evolution of a color signal and color vision. They also suggest that predator visual systems are error prone in the context of mimicry.

Contrasting modes of evolution of the visual pigments in Heliconius butterflies

The adult compound eyes of passion-vine butterflies in the genus Heliconius contain one more UV opsin than other butterflies. Together with an 11-cis-3-hydroxyretinal chromophore, their four opsin genes UVRh1, UVRh2, BRh, and LWRh produce four rhodopsins that are UV-, blue-, or long wavelength absorbing. One of the Heliconius UV opsin genes, UVRh2, was found to have evolved under positive selection following recent gene duplication, using the branch-site test of selection. Using a more conservative test, the small-sample method, we confirm our prior finding of positive selection of UVRh2 and provide new statistical evidence of episodic evolution, that is, positive selection followed by purifying selection. We also newly note that one of the positively selected amino acid sites contains substitutions with known spectral tuning effects in avian ultraviolet- and violet-sensitive visual pigments. As this is one of a handful of described examples of positive selection of any specific gene in any butterfly where functional variation between copies has been characterized, we were interested in examining the molecular and physiological context of this adaptive event by examining the UV opsin genes in contrast to the other visual pigment genes. We cloned BRh and LWRh from 13 heliconiine species and UVRh1 and UVRh2 from Heliconius elevatus. In parallel, we performed in vivo epi-microspectrophotometric experiments to estimate the wavelength of peak absorbance, λ(max), of several rhodopsins in seven heliconiine species. In contrast to UVRh2, we found both physiological and statistical evidence consistent with purifying selection on UVRh1, BRh, and LWRh along the branch leading to the common ancestor of Heliconius. These results underscore the utility of combining molecular and physiological experiments in a comparative context for strengthening evidence for adaptive evolution at the molecular level.

Color discrimination in the red range with only one long-wavelength sensitive opsin

The basic precondition for color vision is the presence of at least two receptor types with different spectral sensitivities. The sensitivity of a receptor is mostly defined by the opsin-based visual pigment expressed in it. We show here, through behavioral experiments, that the nymphalid butterfly Heliconius erato, although it expresses short and medium wavelength opsins and only one long wavelength opsin, discriminates colors in the long-wavelength range (590 nm, 620 nm and 640 nm), whereas another nymphalid, Vanessa atalanta, despite having color vision, is unable to do so. In the eyes of H. erato we identified filtering pigments very close to the rhabdom which differ between ommatidia and produce the yellow and red ommatidial reflection seen under orthodromic illumination. The eyes of V. atalanta lack the filtering pigments, and reflect a homogeneous orange. We hypothesize that the filtering pigments found in the eyes of H. erato may shift the spectral sensitivity peak of the long wavelength receptors in some ommatidia towards longer wavelengths. The comparison of the signals between the two new receptor types makes color discrimination in the red range possible. To our knowledge, this is the first behavioral proof of color vision based on receptors expressing the same opsin.

The eyes of a patrolling butterfly: visual field and eye structure in the Orange Sulphur, Colias eurytheme (Lepidoptera, Pieridae)

Sensory information plays a critical role in determining an animal's behavior on both proximate and evolutionary timescales. Butterflies, like many other insects, use vision extensively over their lifetimes, and yet relatively little work has been published to date on their visual capabilities. We describe the visual system of a pierid butterfly, Colias eurytheme, with the ultimate goal of better understanding its role in shaping the behavior of this animal. We made several measurements: visual field dimensions, eye surface area, interommatidial angle (Deltaphi), facet diameter (D), and eye parameter (p). C. eurytheme had a large visual field and considerable regional variation in visual acuity, as inferred by Deltaphi and D. When compared to females, males had larger eye surface areas, smaller Deltaphi, and larger D in all regions except ventrally. Both sexes had proportionally large eye surface areas compared to other butterflies. Minimum p in males was small, indicating that some regions of their eyes may operate close to the diffraction limit. Finally, we found that both eye surface area and D scaled positively, but with negative allometry to body size. We discuss the relevance of these visual characteristics to the biology and behavior of C. eurytheme.

A genetic linkage map of the mimetic butterfly Heliconius melpomene

Heliconius melpomene is a mimetic butterfly that exhibits great geographic variation in color pattern. We present here a genetic linkage map based on analysis of genetic markers in 73 individuals from a single F(2) family, offspring of a cross between H. m. cythera from western Ecuador and H. m. melpomene from French Guiana. A novel "three-step method" is described for the analysis of dominant markers in an F(2) cross, using outbred parental strains and taking advantage of the lack of crossing over in female Lepidoptera. This method is likely to prove useful for future mapping studies in outbred species with crossing over restricted to one sex, such as the Lepidoptera and Drosophila. The resulting linkage map has 21 linkage groups corresponding to the 21 chromosomes of H. melpomene and includes 219 AFLP markers, 23 microsatellites, 19 single-copy nuclear genes, and the color pattern switch genes Yb and Sb. The marker density is high, averaging >1/7 cM. The total map length is 1616 cM and the average chromosome length is 77 cM. The genome size of H. melpomene was estimated to be 292 Mb, giving a relationship of physical-to-map distance of 180 kb/cM. This map forms the basis for future comparative linkage analysis of color pattern evolution in Heliconius.

Eyeshine and spectral tuning of long wavelength-sensitive rhodopsins: no evidence for red-sensitive photoreceptors among five Nymphalini butterfly species

Spectral tuning of rhodopsins commonly refers to the effects of opsin amino acid substitutions on the wavelength for peak sensitivity of the rhodopsin absorption spectrum. Nymphalini butterflies provide an opportunity for identifying some of the amino acid substitutions responsible for insect rhodopsin spectral tuning because the majority of photoreceptor cells (R3-9)in the adult retina express only a single long wavelength-sensitive (LWS)opsin mRNA transcript. Therefore, the opsin genotype can be directly correlated with its phenotype. We determined the LWS opsin gene sequence from cDNA of the mourning cloak Nymphalis antiopa, and from genomic DNA of the malachite Siproeta stelenes and the peacock Inachis io.Using an epi-microspectrophotometer we examined each butterfly's eyeshine for photochemical evidence of multiple LWS rhodopsins and found only one. We then performed partial-bleaching experiments to obtain absorbance spectra for the LWS rhodopsins of all three species as well as from another nymphalid, the buckeye Junonia coenia. The isolated LWS opsin gene sequences varied in length from 1437-1612 bp and encode rhodopsins R522 (S. stelenes),R530 (I. io), R534 (N. antiopa) and, together with a previously published sequence, R510 (J. coenia). Comparative sequence analysis indicates that the S. stelenes rhodopsin is slightly blue-shifted compared to the typical 530 nm lepidopteran rhodopsin because of the presence of a S138A substitution at a homologous site that in mammalian MWS/LWS rhodopsins causes a 5 nm blue-shift. The difference in peak absorption between R522 of S. stelenes and R530 of Inachis io is therefore largely accounted for by this substitution. This suggests that spectral tuning mechanisms employing the S138A may have evolved in parallel in mammalian and butterfly MWS/LWS rhodopsins across 500 million years of evolution.

Single-copy nuclear genes recover cretaceous-age divergences in bees

We analyzed the higher level phylogeny of the bee family Halictidae based on the coding regions of three single-copy nuclear genes (long-wavelength [LW] opsin, wingless, and elongation factor 1-alpha [EF-1 alpha]). Our combined data set consisted of 2,234 aligned nucleotide sites (702 base pairs [bp] for LW opsin, 405 bp for wingless, and 1,127 bp for EF-1 alpha) and 779 parsimony-informative sites. We included 58 species of halictid bees from 33 genera, representing all subfamilies and tribes, and rooted the trees using seven outgroups from other bee families: Colletidae, Andrenidae, Melittidae, and Apidae. We analyzed the separate and combined data sets by a variety of methods, including equal weights parsimony, maximum likelihood, and Bayesian methods. Analysis of the combined data set produced a strong phylogenetic signal with high bootstrap and Bremer support and high posterior probability well into the base of the tree. The phylogeny recovered the monophyly of the Halictidae and of all four subfamilies and both tribes, recovered relationships among the subfamilies and tribes congruent with morphology, and provided robust support for the relationships among the numerous genera in the tribe Halictini, sensu Michener (2000). Using our combined nucleotide data set, several recently described halictid fossils from the Oligocene and Eocene, and recently developed Bayesian methods, we estimated the antiquity of major clades within the family. Our results indicate that each of the four subfamilies arose well before the Cretaceous-Tertiary boundary and suggest that the early radiation of halictid bees involved substantial African-South American interchange roughly coincident with the separation of these two continents in the late Cretaceous. This combination of single-copy nuclear genes is capable of recovering Cretaceous-age divergences in bees with high levels of support. We propose that LW opsin, wingless, and EF-1 alpha(F2 copy) may be useful in resolving relationships among bee families and other Cretaceous-age insect lineages.

How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets

We analyzed 12 combined mitochondrial and nuclear gene datasets in seven orders of insects using both equal weights parsimony (to evaluate phylogenetic utility) and Bayesian methods (to investigate substitution patterns). For the Bayesian analyses we used relatively complex models (e.g., general time reversible models with rate variation) that allowed us to quantitatively compare relative rates among genes and codon positions, patterns of rate variation among genes, and substitution patterns within genes. Our analyses indicate that nuclear and mitochondrial genes differ in a number of important ways, some of which are correlated with phylogenetic utility. First and most obviously, nuclear genes generally evolve more slowly than mitochondrial genes (except in one case), making them better markers for deep divergences. Second, nuclear genes showed universally high values of CI and (generally) contribute more to overall tree resolution than mitochondrial genes (as measured by partitioned Bremer support). Third, nuclear genes show more homogeneous patterns of among-site rate variation (higher values of alpha than mitochondrial genes). Finally, nuclear genes show more symmetrical transformation rate matrices than mitochondrial genes. The combination of low values of alpha and highly asymmetrical transformation rate matrices may explain the overall poor performance of mitochondrial genes when compared to nuclear genes in the same analysis. Our analyses indicate that some parameters are highly correlated. For example, A/T bias was positively and significantly associated with relative rate and CI was positively and significantly associated with alpha (the shape of the gamma distribution). These results provide important insights into the substitution patterns that might characterized high quality genes for phylogenetic analysis: high values of alpha, unbiased base composition, and symmetrical transformation rate matrices. We argue that insect molecular systematists should increasingly focus on nuclear rather than mitochondrial gene datasets because nuclear genes do not suffer from the same substitutional biases that characterize mitochondrial genes.

Male-killing Wolbachia and mitochondrial DNA: selective sweeps, hybrid introgression and parasite population dynamics

Mitochondrial DNA (mtDNA) sequences are widely used as neutral genetic markers in insects. However, patterns of mtDNA variability are confounded by the spread of maternally transmitted parasites, which are genetically linked to the mitochondria. We have investigated these effects in the butterflies Acraea encedon (which is host to two strains of male-killing Wolbachia bacteria) and A. encedana (which is host to one strain). Within a population, the mitochondria are in linkage disequilibrium with the different male-killers. Furthermore, there has been a recent selective sweep of the mtDNA, which has led to the loss of mitochondrial variation within populations and erased any geographical structure. We also found that one of the male-killers, together with the associated mtDNA, has introgressed from A. encedana into A. encedon within the last 16,000 years. Interestingly, because butterflies are female heterogametic, this will presumably have also led to the introgression of genes on the W sex chromosome. Finally, in A. encedon the mitochondria in uninfected females are unaltered by the spread of the male-killer and have diverse, geographically structured mtDNA. This means we can reject the hypothesis that the male-killer is at a stable equilibrium maintained by imperfect transmission of the bacterium. Instead, some other form of balancing selection may be maintaining uninfected females in the population and preventing the species from going extinct due to a shortage of males.

Not all butterfly eyes are created equal: rhodopsin absorption spectra, molecular identification, and localization of ultraviolet-, blue-, and green-sensitive rhodopsin-encoding mRNAs in the retina of Vanessa cardui

Surveys of spectral sensitivities, visual pigment spectra, and opsin gene sequences have indicated that all butterfly eyes contain ultraviolet-, blue-, and green-sensitive rhodopsins. Some species also contain a fourth or fifth type, related in amino acid sequence to green-sensitive insect rhodopsins, but red shifted in absorbance. By combining electron microscopy, epi-microspectrophotometry, and polymerase chain reaction cloning, we found that the compound eye of Vanessa cardui has the typical ultrastructural features of the butterfly retina but contains only the three common insect rhodopsins. We estimated lambda-max values and relative densities of the rhodopsins in the Vanessa retina (0.72, P530; 0.12, P470; and 0.15, P360) from microspectrophotometric measurements and calculations based on a computational model of reflectance spectra. We isolated three opsin-encoding cDNA fragments that were identified with P530, P470, and P360 by homology to the well-characterized insect rhodopsin families. The retinal mosaic was mapped by opsin mRNA in situ hybridization and found to contain three kinds of ommatidia with respect to their patterns of short wavelength rhodopsin expression. In some ommatidia, P360 or P470 was expressed in R1 and R2 opposed receptor cells; in others, one cell expressed P360, whereas its complement expressed P470. P530 was expressed in the other seven cells of all ommatidia. P470-expressing cells were abundant in the ventral retina but nearly absent dorsally. Our results indicated that there are major differences between the color vision systems of nymphalid and papilionid butterflies: the nymphalid Vanessa has a simpler, trichromatic, system than do the tetrachromatic papilionids that have been studied.

Opsin cDNA sequences of a UV and green rhodopsin of the satyrine butterfly Bicyclus anynana

The cDNAs of an ultraviolet (UV) and long-wavelength (LW) (green) absorbing rhodopsin of the bush brown Bicyclus anynana were partially identified. The UV sequence, encoding 377 amino acids, is 76-79% identical to the UV sequences of the papilionids Papilio glaucus and Papilio xuthus and the moth Manduca sexta. A dendrogram derived from aligning the amino acid sequences reveals an equidistant position of Bicyclus between Papilio and Manduca. The sequence of the green opsin cDNA fragment, which encodes 242 amino acids, represents six of the seven transmembrane regions. At the amino acid level, this fragment is more than 80% identical to the corresponding LW opsin sequences of Dryas, Heliconius, Papilio (rhodopsin 2) and Manduca. Whereas three LW absorbing rhodopsins were identified in the papilionid butterflies, only one green opsin was found in B. anynana.

Functional diversification of lepidopteran opsins following gene duplication

A comparative approach was taken for identifying amino acid substitutions that may be under positive Darwinian selection and are correlated with spectral shifts among orthologous and paralogous lepidopteran long wavelength-sensitive (LW) opsins. Four novel LW opsin fragments were isolated, cloned, and sequenced from eye-specific cDNAs from two butterflies, Vanessa cardui (Nymphalidae) and Precis coenia (Nymphalidae), and two moths, Spodoptera exigua (Noctuidae) and Galleria mellonella (Pyralidae). These opsins were sampled because they encode visual pigments having a naturally occurring range of lambda(max) values (510-530 nm), which in combination with previously characterized lepidopteran opsins, provide a complete range of known spectral sensitivities (510-575 nm) among lepidopteran LW opsins. Two recent opsin gene duplication events were found within the papilionid but not within the nymphalid butterfly families through neighbor-joining, maximum parsimony, and maximum likelihood phylogenetic analyses of 13 lepidopteran opsin sequences. An elevated rate of evolution was detected in the red-shifted Papilio Rh3 branch following gene duplication, because of an increase in the amino acid substitution rate in the transmembrane domain of the protein, a region that forms the chromophore-binding pocket of the visual pigment. A maximum likelihood approach was used to estimate omega, the ratio of nonsynonymous to synonymous substitutions per site. Branch-specific tests of selection (free-ratio) identified one branch with omega = 2.1044, but the small number of substitutions involved was not significantly different from the expected number of changes under the neutral expectation of omega = 1. Ancestral sequences were reconstructed with a high degree of certainty from these data. Reconstructed ancestral sequences revealed several instances of convergence to the same amino acid between butterfly and vertebrate cone pigments, and between independent branches of the butterfly opsin tree that are correlated with spectral shifts.