Localization of Müllerian mimicry genes on a dense linkage map of Heliconius erato

Aculeate Hymenopterans as Aposematic and Mimetic Models

Aposematism and mimicry are complex phenomena which have been studied extensively; however, much of our knowledge comes from just a few focal groups, especially butterflies. Aposematic species combine a warning signal with a secondary defense that reduces their profitability as prey. Aculeate hymenopterans are an extremely diverse lineage defined by the modification of the ovipositor into a stinger which represents a potent defense against predators. Aculeates are often brightly colored and broadly mimicked by members of other arthropod groups including Diptera, Lepidoptera, Coleoptera, and Araneae. However, aculeates are surprisingly understudied as aposematic and mimetic model organisms. Recent studies have described novel pigments contributing to warning coloration in insects and identified changes in cis-regulatory elements as potential drivers of color pattern evolution. Many biotic and abiotic factors contribute to the evolution and maintenance of conspicuous color patterns. Predator distribution and diversity seem to influence the phenotypic diversity of aposematic velvet ants while studies on bumble bees underscore the importance of intermediate mimetic phenotypes in transition zones between putative mimicry rings. Aculeate hymenopterans are attractive models for studying sex-based intraspecific mimicry as male aculeates lack the defense conferred by the females’ stinger. In some species, evolution of male and female color patterns appears to be decoupled. Future studies on aposematic aculeates and their associated mimics hold great promise for unraveling outstanding questions about the evolution of conspicuous color patterns and the factors which determine the composition and distribution of mimetic communities.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation F_n individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

Association between molecular markers and behavioral phenotypes in the immatures of a butterfly

Newly hatched caterpillars of the butterfly Heliconius erato phyllis routinely cannibalize eggs. In a manifestation of kin recognition they cannibalize sibling eggs less frequently than unrelated eggs. Previous work has estimated the heritability of kin recognition in H. erato phyllis to lie between 14 and 48%. It has furthermore been shown that the inheritance of kin recognition is compatible with a quantitative model with a threshold. Here we present the results of a preliminary study, in which we tested for associations between behavioral kin recognition phenotypes and AFLP and SSR markers. We implemented two experimental approaches: (1) a cannibalism test using sibling eggs only, which allowed for only two behavioral outcomes (cannibal and non-cannibal), and (2) a cannibalism test using two sibling eggs and one unrelated egg, which allowed four outcomes [cannibal who does not recognize siblings, cannibal who recognizes siblings, "super-cannibal" (cannibal of both eggs), and "super non-cannibal" (does not cannibalize eggs at all)]. Single-marker analyses were performed using χ2 tests and logistic regression with null markers as covariates. Results of the χ2 tests identified 72 associations for experimental design 1 and 73 associations for design 2. Logistic regression analysis of the markers found to be significant in the χ2 test resulted in 20 associations for design 1 and 11 associations for design 2. Experiment 2 identified markers that were more frequently present or absent in cannibals who recognize siblings and super non-cannibals; i.e. in both phenotypes capable of kin recognition.

Complex modular architecture around a simple toolkit of wing pattern genes

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing color patterns. Here, we used full genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around color pattern genes is highly modular, with narrow genomic intervals associated with specific differences in color and pattern. This modular architecture explains the diversity of color patterns and provides a flexible mechanism for rapid morphological diversification.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.

A lack of Wolbachia-specific DNA in samples from apollo butterfly (Parnassius apollo, Lepidoptera: Papilionidae) individuals with deformed or reduced wings

Various insects contain maternally inherited endosymbiotic bacteria which can cause reproductive alterations, modulation of some physiological responses (like immunity, heat shock response, and oxidative stress response), and resistance to viral infections. In butterflies, Wolbachia sp. is the most frequent endosymbiont from this group, occurring in about 30 % of species tested to date. In this report, the presence of Wolbachia-specific DNA has been detected in apollo butterfly (Parnassius apollo). In the isolated population of this insect occurring in Pieniny National Park (Poland), malformed individuals with deformed or reduced wings appear with an exceptionally high frequency. Interestingly, while total DNA isolated from most (about 85 %) normal insects contained Wolbachia-specific sequences detected by PCR, such sequences were absent in a large fraction (70 %) of individuals with deformed wings and in all tested individuals with reduced wings. These results indicate for the first time the correlation between malformation of wings and the absence of Wolbachia sp. in insects. Although the lack of the endosymbiotic bacteria cannot be considered as the sole cause of the deformation or reduction of wings, one might suggest that Wolbachia sp. could play a protective role in the ontogenetic development of apollo butterfly.

Divergence with gene flow across a speciation continuum of Heliconius butterflies

Background

A key to understanding the origins of species is determining the evolutionary processes that drive the patterns of genomic divergence during speciation. New genomic technologies enable the study of high-resolution genomic patterns of divergence across natural speciation continua, where taxa pairs with different levels of reproductive isolation can be used as proxies for different stages of speciation. Empirical studies of these speciation continua can provide valuable insights into how genomes diverge during speciation.

Methods

We examine variation across a handful of genomic regions in parapatric and allopatric populations of Heliconius butterflies with varying levels of reproductive isolation. Genome sequences were mapped to 2.2-Mb of the H. erato genome, including 1-Mb across the red color pattern locus and multiple regions unlinked to color pattern variation.

Results

Phylogenetic analyses reveal a speciation continuum of pairs of hybridizing races and incipient species in the Heliconius erato clade. Comparisons of hybridizing pairs of divergently colored races and incipient species reveal that genomic divergence increases with ecological and reproductive isolation, not only across the locus responsible for adaptive variation in red wing coloration, but also at genomic regions unlinked to color pattern.

Discussion

We observe high levels of divergence between the incipient species H. erato and H. himera, suggesting that divergence may accumulate early in the speciation process. Comparisons of genomic divergence between the incipient species and allopatric races suggest that limited gene flow cannot account for the observed high levels of divergence between the incipient species.

Conclusions

Our results provide a reconstruction of the speciation continuum across the H. erato clade and provide insights into the processes that drive genomic divergence during speciation, establishing the H. erato clade as a powerful framework for the study of speciation.

Electronic supplementary material

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

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.

Unraveling the Tangled Skein: The Evolution of Transcriptional Regulatory Networks in Development

The molecular and genetic basis for the evolution of anatomical diversity is a major question that has inspired evolutionary and developmental biologists for decades. Because morphology takes form during development, a true comprehension of how anatomical structures evolve requires an understanding of the evolutionary events that alter developmental genetic programs. Vast gene regulatory networks (GRNs) that connect transcription factors to their target regulatory sequences control gene expression in time and space and therefore determine the tissue-specific genetic programs that shape morphological structures. In recent years, many new examples have greatly advanced our understanding of the genetic alterations that modify GRNs to generate newly evolved morphologies. Here, we review several aspects of GRN evolution, including their deep preservation, their mechanisms of alteration, and how they originate to generate novel developmental programs.

Molecular basis of wing coloration in a Batesian mimic butterfly, Papilio polytes

Batesian mimicry protects animals from predators through resemblance with distasteful models in shape, color pattern, or behavior. To elucidate the wing coloration mechanisms involved in the mimicry, we investigated chemical composition and gene expression of the pale yellow and red pigments of a swallowtail butterfly, Papilio polytes, whose females mimic the unpalatable butterfly Pachliopta aristolochiae. Using LC/MS, we showed that the pale yellow wing regions in non-mimetic females consist of kynurenine and N-β-alanyldopamine (NBAD). Moreover, qRT-PCR showed that kynurenine/NBAD biosynthetic genes were upregulated in these regions in non-mimetic females. However, these pigments were absent in mimetic females. RNA-sequencing showed that kynurenine/NBAD synthesis and Toll signaling genes were upregulated in the red spots specific to mimetic female wings. These results demonstrated that drastic changes in gene networks in the red and pale yellow regions can switch wing color patterns between non-mimetic and mimetic females of P. polytes.

Ultrastructure and morphogenesis of the wing scales in Heliconius erato phyllis (Lepidoptera: Nymphalidae): what silvery/brownish surfaces can tell us about the development of color patterning?

Usually the literature on Heliconius show three types of scales, classified based on the correlation between color and ultrastructure: type I - white and yellow, type II - black, and type III - orange and red. The ultrastructure of the scales located at the silvery/brownish surfaces of males/females is for the first time described in this paper. Besides, we describe the ontogeny of pigmentation, the scale morphogenesis and the maturation timing of scales fated to different colors in Heliconius erato phyllis. The silvery/brownish surfaces showed ultrastructurally similar scales to the type I, II and III. The ontogeny of pigmentation follows the sequence red, black, silvery/brownish and yellow. The maturation of yellow-fated scales, however, occurred simultaneously with the red-fated scales, before the pigmentation becomes visible. In spite of the scales at the silvery/brownish surfaces being ultrastructurally similar to the yellow, red and black scales, they mature after them; this suggests that the maturation timing does not show a relationship with the scale ultrastructure, with the deposition timing of the yellow pigment. The analysis of H. erato phyllis scale morphogenesis, as well as the scales ultrastructure and maturation timing, provided new findings into the developmental architecture of color pattern in Heliconius.

Genomic architecture of adaptive color pattern divergence and convergence in Heliconius butterflies

Identifying the genetic changes driving adaptive variation in natural populations is key to understanding the origins of biodiversity. The mosaic of mimetic wing patterns in Heliconius butterflies makes an excellent system for exploring adaptive variation using next-generation sequencing. In this study, we use a combination of techniques to annotate the genomic interval modulating red color pattern variation, identify a narrow region responsible for adaptive divergence and convergence in Heliconius wing color patterns, and explore the evolutionary history of these adaptive alleles. We use whole genome resequencing from four hybrid zones between divergent color pattern races of Heliconius erato and two hybrid zones of the co-mimic Heliconius melpomene to examine genetic variation across 2.2 Mb of a partial reference sequence. In the intergenic region near optix, the gene previously shown to be responsible for the complex red pattern variation in Heliconius, population genetic analyses identify a shared 65-kb region of divergence that includes several sites perfectly associated with phenotype within each species. This region likely contains multiple cis-regulatory elements that control discrete expression domains of optix. The parallel signatures of genetic differentiation in H. erato and H. melpomene support a shared genetic architecture between the two distantly related co-mimics; however, phylogenetic analysis suggests mimetic patterns in each species evolved independently. Using a combination of next-generation sequencing analyses, we have refined our understanding of the genetic architecture of wing pattern variation in Heliconius and gained important insights into the evolution of novel adaptive phenotypes in natural populations.

Multi-allelic major effect genes interact with minor effect QTLs to control adaptive color pattern variation in Heliconius erato

Recent studies indicate that relatively few genomic regions are repeatedly involved in the evolution of Heliconius butterfly wing patterns. Although this work demonstrates a number of cases where homologous loci underlie both convergent and divergent wing pattern change among different Heliconius species, it is still unclear exactly how many loci underlie pattern variation across the genus. To address this question for Heliconius erato, we created fifteen independent crosses utilizing the four most distinct color pattern races and analyzed color pattern segregation across a total of 1271 F2 and backcross offspring. Additionally, we used the most variable brood, an F2 cross between H. himera and the east Ecuadorian H. erato notabilis, to perform a quantitative genetic analysis of color pattern variation and produce a detailed map of the loci likely involved in the H. erato color pattern radiation. Using AFLP and gene based markers, we show that fewer major genes than previously envisioned control the color pattern variation in H. erato. We describe for the first time the genetic architecture of H. erato wing color pattern by assessing quantitative variation in addition to traditional linkage mapping. In particular, our data suggest three genomic intervals modulate the bulk of the observed variation in color. Furthermore, we also identify several modifier loci of moderate effect size that contribute to the quantitative wing pattern variation. Our results are consistent with the two-step model for the evolution of mimetic wing patterns in Heliconius and support a growing body of empirical data demonstrating the importance of major effect loci in adaptive change.

Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n=28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n=31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.

Unraveling the thread of nature's tapestry: the genetics of diversity and convergence in animal pigmentation

Animals display incredibly diverse color patterns yet little is known about the underlying genetic basis of these phenotypes. However, emerging results are reshaping our view of how the process of phenotypic evolution occurs. Here, we outline recent research from three particularly active areas of investigation: melanin pigmentation in Drosophila, wing patterning in butterflies, and pigment variation in lizards. For each system, we highlight (i) the function and evolution of color variation, (ii) various approaches that have been used to explore the genetic basis of pigment variation, and (iii) conclusions regarding the genetic basis of convergent evolution which have emerged from comparative analyses. Results from these studies indicate that natural variation in pigmentation is a particularly powerful tool to examine the molecular basis of evolution, especially with regard to convergent or parallel evolution. Comparison of these systems also reveals that the molecular basis of convergent evolution is heterogeneous, sometimes involving conserved mechanisms and sometimes not. In the near future, additional work in other emerging systems will substantially expand the scope of available comparisons.

Diversification of complex butterfly wing patterns by repeated regulatory evolution of a Wnt ligand

Although animals display a rich variety of shapes and patterns, the genetic changes that explain how complex forms arise are still unclear. Here we take advantage of the extensive diversity of Heliconius butterflies to identify a gene that causes adaptive variation of black wing patterns within and between species. Linkage mapping in two species groups, gene-expression analysis in seven species, and pharmacological treatments all indicate that cis-regulatory evolution of the WntA ligand underpins discrete changes in color pattern features across the Heliconius genus. These results illustrate how the direct modulation of morphogen sources can generate a wide array of unique morphologies, thus providing a link between natural genetic variation, pattern formation, and adaptation.

Transcriptome analysis reveals novel patterning and pigmentation genes underlying Heliconius butterfly wing pattern variation

Background

Heliconius butterfly wing pattern diversity offers a unique opportunity to investigate how natural genetic variation can drive the evolution of complex adaptive phenotypes. Positional cloning and candidate gene studies have identified a handful of regulatory and pigmentation genes implicated in Heliconius wing pattern variation, but little is known about the greater developmental networks within which these genes interact to pattern a wing. Here we took a large-scale transcriptomic approach to identify the network of genes involved in Heliconius wing pattern development and variation. This included applying over 140 transcriptome microarrays to assay gene expression in dissected wing pattern elements across a range of developmental stages and wing pattern morphs of Heliconius erato.

Results

We identified a number of putative early prepattern genes with color-pattern related expression domains. We also identified 51 genes differentially expressed in association with natural color pattern variation. Of these, the previously identified color pattern “switch gene” optix was recovered as the first transcript to show color-specific differential expression. Most differentially expressed genes were transcribed late in pupal development and have roles in cuticle formation or pigment synthesis. These include previously undescribed transporter genes associated with ommochrome pigmentation. Furthermore, we observed upregulation of melanin-repressing genes such as ebony and Dat1 in non-melanic patterns.

Conclusions

This study identifies many new genes implicated in butterfly wing pattern development and provides a glimpse into the number and types of genes affected by variation in genes that drive color pattern evolution.

Evolution of a mimicry supergene from a multilocus architecture

The origin and evolution of supergenes have long fascinated evolutionary biologists. In the polymorphic butterfly Heliconius numata, a supergene controls the switch between multiple different forms, and results in near-perfect mimicry of model species. Here, we use a morphometric analysis to quantify the variation in wing pattern observed in two broods of H. numata with different alleles at the supergene locus, 'P'. Further, we genotype the broods to associate the variation we capture with genetic differences. This allows us to begin mapping the quantitative trait loci that have minor effects on wing pattern. In addition to finding loci on novel chromosomes, our data, to our knowledge, suggest for the first time that ancestral colour-pattern loci, known to have major effects in closely related species, may contribute to the wing patterns displayed by H. numata, despite the large transfer of effects to the supergene.

Comparative population genetics of a mimicry locus among hybridizing Heliconius butterfly species

The comimetic Heliconius butterfly species pair, H. erato and H. melpomene, appear to use a conserved Mendelian switch locus to generate their matching red wing patterns. Here we investigate whether H. cydno and H. pachinus, species closely related to H. melpomene, use this same switch locus to generate their highly divergent red and brown color pattern elements. Using an F2 intercross between H. cydno and H. pachinus, we first map the genomic positions of two novel red/brown wing pattern elements; the G locus, which controls the presence of red vs brown at the base of the ventral wings, and the Br locus, which controls the presence vs absence of a brown oval pattern on the ventral hind wing. The results reveal that the G locus is tightly linked to markers in the genomic interval that controls red wing pattern elements of H. erato and H. melpomene. Br is on the same linkage group but approximately 26 cM away. Next, we analyze fine-scale patterns of genetic differentiation and linkage disequilibrium throughout the G locus candidate interval in H. cydno, H. pachinus and H. melpomene, and find evidence for elevated differentiation between H. cydno and H. pachinus, but no localized signature of association. Overall, these results indicate that the G locus maps to the same interval as the locus controlling red patterning in H. melpomene and H. erato. This, in turn, suggests that the genes controlling red pattern elements may be homologous across Heliconius, supporting the hypothesis that Heliconius butterflies use a limited suite of conserved genetic switch loci to generate both convergent and divergent wing patterns.

Hybrid speciation in Heliconius butterflies? A review and critique of the evidence

The evidence supporting the recent hypothesis of a homoploid hybrid origin for the butterfly species Heliconius heurippa is evaluated. Data from selective breeding experiments, mate-choice studies, and a wide variety of DNA markers are reviewed, and an alternative hypothesis for the origin of the species and its close relatives is proposed. A scenario of occasional red wing-pattern mutations in peripheral populations of Heliconius cydno with subsequent adaptive convergence towards sympatric mimicry rings involving H. melpomene and H. erato is offered as an alternative to the HHS hypothesis. Recent twists of this tale are addressed in a postscript.

Combined EST and proteomic analysis identifies rapidly evolving seminal fluid proteins in Heliconius butterflies

Seminal fluid proteins (SFPs) directly influence a wide range of reproductive processes, including fertilization, sperm storage, egg production, and immune response. Like many other reproductive proteins, the molecular evolution of SFPs is generally characterized by rapid and frequently adaptive evolution. However, the evolutionary processes underlying this often-documented pattern have not yet been confidently determined. A robust understanding of the processes governing SFP evolution will ultimately require identifying SFPs and characterizing their evolution in many different taxa, often where only limited genomic resources are available. Here, we report the first comprehensive molecular genetic and evolutionary analysis of SFPs conducted in Lepidoptera (moths and butterflies). We have identified 51 novel SFPs from two species of Heliconius butterflies (Heliconius erato and Heliconius melpomene) by combining "indirect" bioinformatic and expression analyses of expressed sequence tags from male accessory gland and wing tissues with "direct" proteomic analyses of spermatophores. Proteomic analyses identified fewer SFPs than the indirect criteria but gave consistent results. Of 51 SFPs, 40 were identified in both species but fewer than half could be functionally annotated via similarity searches (Blast, IPRscan, etc.). The majority of annotated Heliconius SFPs were predicted to be chymotrypsins. Comparisons of Heliconius SFPs with those from fruit fly, mosquito, honeybee, and cricket suggest that gene turnover is high among these proteins and that SFPs are rarely conserved across insect orders. Pairwise estimates of evolutionary rates between SFPs and nonreproductive proteins show that, on average, Heliconius SFPs are evolving rapidly. At least one of these SFPs is evolving adaptively (dN/dS > 1), implicating a role for positive selection in this rapid evolution. This work establishes a strong precedent for future research on the causes and consequences of reproductive protein evolution in the Lepidoptera. Butterflies and moths have an extremely rich history of organismal research, which will provide an informative ecological context for further molecular evolutionary investigations.

Genomic hotspots for adaptation: the population genetics of Müllerian mimicry in Heliconius erato

Wing pattern evolution in Heliconius butterflies provides some of the most striking examples of adaptation by natural selection. The genes controlling pattern variation are classic examples of Mendelian loci of large effect, where allelic variation causes large and discrete phenotypic changes and is responsible for both convergent and highly divergent wing pattern evolution across the genus. We characterize nucleotide variation, genotype-by-phenotype associations, linkage disequilibrium (LD), and candidate gene expression patterns across two unlinked genomic intervals that control yellow and red wing pattern variation among mimetic forms of Heliconius erato. Despite very strong natural selection on color pattern, we see neither a strong reduction in genetic diversity nor evidence for extended LD across either patterning interval. This observation highlights the extent that recombination can erase the signature of selection in natural populations and is consistent with the hypothesis that either the adaptive radiation or the alleles controlling it are quite old. However, across both patterning intervals we identified SNPs clustered in several coding regions that were strongly associated with color pattern phenotype. Interestingly, coding regions with associated SNPs were widely separated, suggesting that color pattern alleles may be composed of multiple functional sites, conforming to previous descriptions of these loci as "supergenes." Examination of gene expression levels of genes flanking these regions in both H. erato and its co-mimic, H. melpomene, implicate a gene with high sequence similarity to a kinesin as playing a key role in modulating pattern and provides convincing evidence for parallel changes in gene regulation across co-mimetic lineages. The complex genetic architecture at these color pattern loci stands in marked contrast to the single casual mutations often identified in genetic studies of adaptation, but may be more indicative of the type of genetic changes responsible for much of the adaptive variation found in natural populations.

EST and microarray analysis of horn development in Onthophagus beetles

Background

The origin of novel traits and their subsequent diversification represent central themes in evo-devo and evolutionary ecology. Here we explore the genetic and genomic basis of a class of traits that is both novel and highly diverse, in a group of organisms that is ecologically complex and experimentally tractable: horned beetles.

Results

We developed two high quality, normalized cDNA libraries for larval and pupal Onthophagus taurus and sequenced 3,488 ESTs that assembled into 451 contigs and 2,330 singletons. We present the annotation and a comparative analysis of the conservation of the sequences. Microarrays developed from the combined libraries were then used to contrast the transcriptome of developing primordia of head horns, prothoracic horns, and legs. Our experiments identify a first comprehensive list of candidate genes for the evolution and diversification of beetle horns. We find that developing horns and legs show many similarities as well as important differences in their transcription profiles, suggesting that the origin of horns was mediated partly, but not entirely, by the recruitment of genes involved in the formation of more traditional appendages such as legs. Furthermore, we find that horns developing from the head and prothorax differ in their transcription profiles to a degree that suggests that head and prothoracic horns are not serial homologs, but instead may have evolved independently from each other.

Conclusion

We have laid the foundation for a systematic analysis of the genetic basis of horned beetle development and diversification with the potential to contribute significantly to several major frontiers in evolutionary developmental biology.

Microsatellite markers associated with genes expressed in developing wings of Bicyclus anynana butterflies

Deriving useful microsatellite markers in lepidopterans has been challenging when relying on scans of genomic DNA libraries, presumably due to repetitiveness in their genomes. We assayed 96 of 320 microsatellites identified in silico from a collection of Bicyclus anynana ESTs, in 11 independent individuals from a laboratory population. From the 68 successful assays, we identified 40 polymorphic markers including 22 with BLAST-based annotation. Nine of 12 selected polymorphic markers tested in a panel of 24 wild-caught individuals converted to successful assays and were all polymorphic. We discuss how microsatellite discovery in ESTs is an efficient strategy with important attendant advantages.

A gene-based linkage map for Bicyclus anynana butterflies allows for a comprehensive analysis of synteny with the lepidopteran reference genome

Lepidopterans (butterflies and moths) are a rich and diverse order of insects, which, despite their economic impact and unusual biological properties, are relatively underrepresented in terms of genomic resources. The genome of the silkworm Bombyx mori has been fully sequenced, but comparative lepidopteran genomics has been hampered by the scarcity of information for other species. This is especially striking for butterflies, even though they have diverse and derived phenotypes (such as color vision and wing color patterns) and are considered prime models for the evolutionary and developmental analysis of ecologically relevant, complex traits. We focus on Bicyclus anynana butterflies, a laboratory system for studying the diversification of novelties and serially repeated traits. With a panel of 12 small families and a biphasic mapping approach, we first assigned 508 expressed genes to segregation groups and then ordered 297 of them within individual linkage groups. We also coarsely mapped seven color pattern loci. This is the richest gene-based map available for any butterfly species and allowed for a broad-coverage analysis of synteny with the lepidopteran reference genome. Based on 462 pairs of mapped orthologous markers in Bi. anynana and Bo. mori, we observed strong conservation of gene assignment to chromosomes, but also evidence for numerous large- and small-scale chromosomal rearrangements. With gene collections growing for a variety of target organisms, the ability to place those genes in their proper genomic context is paramount. Methods to map expressed genes and to compare maps with relevant model systems are crucial to extend genomic-level analysis outside classical model species. Maps with gene-based markers are useful for comparative genomics and to resolve mapped genomic regions to a tractable number of candidate genes, especially if there is synteny with related model species. This is discussed in relation to the identification of the loci contributing to color pattern evolution in butterflies.

Evolution of coloration patterns

There is an amazing amount of diversity in coloration patterns in nature. The ease of observing this diversity and the recent application of genetic and molecular techniques to model and nonmodel animals are allowing us to investigate the genetic basis and evolution of coloration in an ever-increasing variety of animals. It is now possible to ask questions about how many genes are responsible for any given pattern, what types of genetic changes have occurred to generate the diversity, and if the same underlying genetic changes occur repeatedly when coloration phenotypes arise through convergent evolution or parallel evolution.

Cytogenetic characterization and AFLP-based genetic linkage mapping for the butterfly Bicyclus anynana, covering all 28 karyotyped chromosomes

Background

The chromosome characteristics of the butterfly Bicyclus anynana, have received little attention, despite the scientific importance of this species. This study presents the characterization of chromosomes in this species by means of cytogenetic analysis and linkage mapping.

Methodology/Principal Findings

Physical genomic features in the butterfly B. anynana were examined by karyotype analysis and construction of a linkage map. Lepidoptera possess a female heterogametic W-Z sex chromosome system. The WZ-bivalent in pachytene oocytes of B. anynana consists of an abnormally small, heterochromatic W-chromosome with the Z-chromosome wrapped around it. Accordingly, the W-body in interphase nuclei is much smaller than usual in Lepidoptera. This suggests an intermediate stage in the process of secondary loss of the W-chromosome to a ZZ/Z sex determination system. Two nucleoli are present in the pachytene stage associated with an autosome and the WZ-bivalent respectively. Chromosome counts confirmed a haploid number of n = 28. Linkage mapping had to take account of absence of crossing-over in females, and of our use of a full-sib crossing design. We developed a new method to determine and exclude the non-recombinant uninformative female inherited component in offspring. The linkage map was constructed using a novel approach that uses exclusively JOINMAP-software for Lepidoptera linkage mapping. This approach simplifies the mapping procedure, avoids over-estimation of mapping distance and increases the reliability of relative marker positions. A total of 347 AFLP markers, 9 microsatellites and one single-copy nuclear gene covered all 28 chromosomes, with a mapping distance of 1354 cM. Conserved synteny of Tpi on the Z-chromosome in Lepidoptera was confirmed for B. anynana. The results are discussed in relation to other mapping studies in Lepidoptera.

Conclusions/Significance

This study adds to the knowledge of chromosome structure and evolution of an intensively studied organism. On a broader scale it provides an insight in Lepidoptera sex chromosome evolution and it proposes a simpler and more reliable method of linkage mapping than used for Lepidoptera to date.

EST analysis of male accessory glands from Heliconius butterflies with divergent mating systems

Background

Heliconius butterflies possess a remarkable diversity of phenotypes, physiologies, and behaviors that has long distinguished this genus as a focal taxon in ecological and evolutionary research. Recently Heliconius has also emerged as a model system for using genomic methods to investigate the causes and consequences of biological diversity. One notable aspect of Heliconius diversity is a dichotomy in mating systems which provides an unusual opportunity to investigate the relationship between sexual selection and the evolution of reproductive proteins. As a first step in pursuing this research, we report the generation and analysis of expressed sequence tags (ESTs) from the male accessory gland of H. erato and H. melpomene, species representative of the two mating systems present in the genus Heliconius.

Results

We successfully sequenced 933 ESTs clustering into 371 unigenes from H. erato and 1033 ESTs clustering into 340 unigenes from H. melpomene. Results from the two species were very similar. Approximately a third of the unigenes showed no significant BLAST similarity (E-value <10^-5) to sequences in GenBank's non-redundant databases, indicating that a large proportion of novel genes are expressed in Heliconius male accessory glands. In both species only a third of accessory gland unigenes were also found among genes expressed in wing tissue. About 25% of unigenes from both species encoded secreted proteins. This includes three groups of highly abundant unigenes encoding repetitive proteins considered to be candidate seminal fluid proteins; proteins encoded by one of these groups were detected in H. erato spermatophores.

Conclusion

This collection of ESTs will serve as the foundation for the future identification and evolutionary analysis of male reproductive proteins in Heliconius butterflies. These data also represent a significant advance in the rapidly growing collection of genomic resources available in Heliconius butterflies. As such, they substantially enhance this taxon as a model system for investigating questions of ecological, phenotypic, and genomic diversity.

Convergent evolution in the genetic basis of Müllerian mimicry in heliconius butterflies

The neotropical butterflies Heliconius melpomene and H. erato are Müllerian mimics that display the same warningly colored wing patterns in local populations, yet pattern diversity between geographic regions. Linkage mapping has previously shown convergent red wing phenotypes in these species are controlled by loci on homologous chromosomes. Here, AFLP bulk segregant analysis using H. melpomene crosses identified genetic markers tightly linked to two red wing-patterning loci. These markers were used to screen a H. melpomene BAC library and a tile path was assembled spanning one locus completely and part of the second. Concurrently, a similar strategy was used to identify a BAC clone tightly linked to the locus controlling the mimetic red wing phenotypes in H. erato. A methionine rich storage protein (MRSP) gene was identified within this BAC clone, and comparative genetic mapping shows red wing color loci are in homologous regions of the genome of H. erato and H. melpomene. Subtle differences in these convergent phenotypes imply they evolved independently using somewhat different developmental routes, but are nonetheless regulated by the same switch locus. Genetic mapping of MRSP in a third related species, the "tiger" patterned H. numata, has no association with wing patterning and shows no evidence for genomic translocation of wing-patterning loci.

Colour pattern specification in the Mocker swallowtail Papilio dardanus: the transcription factor invected is a candidate for the mimicry locus H

The swallowtail butterfly, Papilio dardanus, is an iconic example of a polymorphic Batesian mimic. The expression of various female-limited colour forms is thought to be controlled by a single autosomal locus, termed H, whose function in determining the wing pattern remains elusive. As a step towards the physical mapping of H, we established a set of 272 polymorphic amplified fragment length polymorphism (AFLP) markers (EcoRI-MseI). Segregation patterns in a 'female-informative' brood (exploiting the absence of crossing over in female Lepidoptera) mapped these AFLPs to 30 linkage groups (putative chromosomes). The difference between the hippocoon and cenea female forms segregating in this family resides on a single one of these linkage groups, defined by 14 AFLPs. In a 'male-informative' cross (markers segregating within a linkage group), a pair of AFLPs co-segregated closely with the two female forms, except in four recombinants out of 19 female offspring. Linkage with these AFLP markers using four further female-informative families demonstrated that the genetic factor determining other morphs (poultoni, lamborni and trimeni) also maps to this same linkage group. The candidate gene invected, obtained in a screen for co-segregation of developmental genes with the colour forms, resides in a 13.9 cM interval flanked by the two AFLP markers. In the male-informative family invected co-segregated perfectly with the hippocoon/cenea factor, despite the four crossovers with the AFLPs. These findings make invected, and possibly its closely linked paralogue engrailed, strong candidates for H. This is supported by their known role in eyespot specification in nymphalid butterfly wings.

Highly conserved gene order and numerous novel repetitive elements in genomic regions linked to wing pattern variation in Heliconius butterflies

BACKGROUND

With over 20 parapatric races differing in their warningly colored wing patterns, the butterfly Heliconius erato provides a fascinating example of an adaptive radiation. Together with matching races of its co-mimic Heliconius melpomene, H. erato also represents a textbook case of Müllerian mimicry, a phenomenon where common warning signals are shared amongst noxious organisms. It is of great interest to identify the specific genes that control the mimetic wing patterns of H. erato and H. melpomene. To this end we have undertaken comparative mapping and targeted genomic sequencing in both species. This paper reports on a comparative analysis of genomic sequences linked to color pattern mimicry genes in Heliconius.

RESULTS

Scoring AFLP polymorphisms in H. erato broods allowed us to survey loci at approximately 362 kb intervals across the genome. With this strategy we were able to identify markers tightly linked to two color pattern genes: D and Cr, which were then used to screen H. erato BAC libraries in order to identify clones for sequencing. Gene density across 600 kb of BAC sequences appeared relatively low, although the number of predicted open reading frames was typical for an insect. We focused analyses on the D- and Cr-linked H. erato BAC sequences and on the Yb-linked H. melpomene BAC sequence. A comparative analysis between homologous regions of H. erato (Cr-linked BAC) and H. melpomene (Yb-linked BAC) revealed high levels of sequence conservation and microsynteny between the two species. We found that repeated elements constitute 26% and 20% of BAC sequences from H. erato and H. melpomene respectively. The majority of these repetitive sequences appear to be novel, as they showed no significant similarity to any other available insect sequences. We also observed signs of fine scale conservation of gene order between Heliconius and the moth Bombyx mori, suggesting that lepidopteran genome architecture may be conserved over very long evolutionary time scales.

CONCLUSION

Here we have demonstrated the tractability of progressing from a genetic linkage map to genomic sequence data in Heliconius butterflies. We have also shown that fine-scale gene order is highly conserved between distantly related Heliconius species, and also between Heliconius and B. mori. Together, these findings suggest that genome structure in macrolepidoptera might be very conserved, and show that mapping and positional cloning efforts in different lepidopteran species can be reciprocally informative.

Gene expression underlying adaptive variation in Heliconius wing patterns: non-modular regulation of overlapping cinnabar and vermilion prepatterns

Geographical variation in the mimetic wing patterns of the butterfly Heliconius erato is a textbook example of adaptive polymorphism; however, little is known about how this variation is controlled developmentally. Using microarrays and qPCR, we identified and compared expression of candidate genes potentially involved with a red/yellow forewing band polymorphism in H. erato. We found that transcripts encoding the pigment synthesis enzymes cinnabar and vermilion showed pattern- and polymorphism-related expression patterns, respectively. cinnabar expression was associated with the forewing band regardless of pigment colour, providing the first gene expression pattern known to be correlated with a major Heliconius colour pattern. In contrast, vermilion expression changed spatially over time in red-banded butterflies, but was not expressed at detectable levels in yellow-banded butterflies, suggesting that regulation of this gene may be involved with the red/yellow polymorphism. Furthermore, we found that the yellow pigment, 3-hydroxykynurenine, is incorporated into wing scales from the haemolymph rather than being synthesized in situ. We propose that some aspects of Heliconius colour patterns are determined by spatio-temporal overlap of pigment gene transcription prepatterns and speculate that evolutionary changes in vermilion regulation may in part underlie an adaptive colour pattern polymorphism.

Adaptive variation in beach mice produced by two interacting pigmentation genes

Little is known about the genetic basis of ecologically important morphological variation such as the diverse color patterns of mammals. Here we identify genetic changes contributing to an adaptive difference in color pattern between two subspecies of oldfield mice (Peromyscus polionotus). One mainland subspecies has a cryptic dark brown dorsal coat, while a younger beach-dwelling subspecies has a lighter coat produced by natural selection for camouflage on pale coastal sand dunes. Using genome-wide linkage mapping, we identified three chromosomal regions (two of major and one of minor effect) associated with differences in pigmentation traits. Two candidate genes, the melanocortin-1 receptor (Mc1r) and its antagonist, the Agouti signaling protein (Agouti), map to independent regions that together are responsible for most of the difference in pigmentation between subspecies. A derived mutation in the coding region of Mc1r, rather than change in its expression level, contributes to light pigmentation. Conversely, beach mice have a derived increase in Agouti mRNA expression but no changes in protein sequence. These two genes also interact epistatically: the phenotypic effects of Mc1r are visible only in genetic backgrounds containing the derived Agouti allele. These results demonstrate that cryptic coloration can be based largely on a few interacting genes of major effect.

The tremendous amount of variation in color patterns among organisms helps individuals survive and reproduce in the wild, yet we know surprisingly little about the genes that produce these adaptive patterns. Here we used a genomic analysis to uncover the molecular basis of a pale color pattern that camouflages beach mice inhabiting the sandy dunes of Florida's coast from predators. We identified two pigmentation genes, the melanocortin-1 receptor (Mc1r) and its ligand, the agouti signaling protein (Agouti), which together produce a light color pattern. We show that this light pigmentation results partly from a single amino acid mutation in Mc1r, which reduces the activity of the receptor but does not affect the gene's expression level, and partly from the derived Agouti allele, which shows no change in protein sequence but does exhibit an increase in mRNA expression. We also show that these two genes do not act additively to produce pale color; rather, the derived Agouti allele must be present to see any effect of Mc1r on pigmentation. Thus, the light color pattern of beach mice largely results from the physical interaction between a structural change in a receptor (reducing Mc1r activity) and a regulatory change in the receptor's antagonist (increasing Agouti expression).

Species of oldfield mice have coat colors adapted for their environment. By using genome-wide linkage mapping, the authors show that three chromosomal regions are associated with differences in pigmentation traits.

Synteny and chromosome evolution in the lepidoptera: evidence from mapping in Heliconius melpomene

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.

Patterns on the insect wing

The evolution of wings and the adaptive advantage they provide have allowed insects to become one of the most evolutionarily successful groups on earth. The incredible diversity of their shape, size, and color patterns is a direct reflection of the important role wings have played in the radiation of insects. In this review, we highlight recent studies on both butterflies and Drosophila that have begun to uncover the types of genetic variations and developmental mechanisms that control diversity in wing color patterns. In butterflies, these analyses are now possible because of the recent development of a suite of genomic and functional tools, such as detailed linkage maps and transgenesis. In one such study, extensive linkage mapping in Heliconius butterflies has shown that surprisingly few, and potentially homologous, loci are responsible for several major pattern variations on the wings of these butterflies. Parallel work on a clade of Drosophila has uncovered how cis-regulatory changes of the same gene correlate with the repeated gain and loss of pigmented wing spots. Collectively, our understanding of formation and evolution of color pattern in insect wings is rapidly advancing because of these recent breakthroughs in several different fields.

Natural hybridization in heliconiine butterflies: the species boundary as a continuum

BACKGROUND

To understand speciation and the maintenance of taxa as separate entities, we need information about natural hybridization and gene flow among species.

RESULTS

Interspecific hybrids occur regularly in Heliconius and Eueides (Lepidoptera: Nymphalidae) in the wild: 26-29% of the species of Heliconiina are involved, depending on species concept employed. Hybridization is, however, rare on a per-individual basis. For one well-studied case of species hybridizing in parapatric contact (Heliconius erato and H. himera), phenotypically detectable hybrids form around 10% of the population, but for species in sympatry hybrids usually form less than 0.05% of individuals. There is a roughly exponential decline with genetic distance in the numbers of natural hybrids in collections, both between and within species, suggesting a simple "exponential failure law" of compatibility as found in some prokaryotes.

CONCLUSION

Hybridization between species of Heliconius appears to be a natural phenomenon; there is no evidence that it has been enhanced by recent human habitat disturbance. In some well-studied cases, backcrossing occurs in the field and fertile backcrosses have been verified in insectaries, which indicates that introgression is likely, and recent molecular work shows that alleles at some but not all loci are exchanged between pairs of sympatric, hybridizing species. Molecular clock dating suggests that gene exchange may continue for more than 3 million years after speciation. In addition, one species, H. heurippa, appears to have formed as a result of hybrid speciation. Introgression may often contribute to adaptive evolution as well as sometimes to speciation itself, via hybrid speciation. Geographic races and species that coexist in sympatry therefore form part of a continuum in terms of hybridization rates or probability of gene flow. This finding concurs with the view that processes leading to speciation are continuous, rather than sudden, and that they are the same as those operating within species, rather than requiring special punctuated effects or complete allopatry. Although not qualitatively distinct from geographic races, nor "real" in terms of phylogenetic species concepts or the biological species concept, hybridizing species of Heliconius are stably distinct in sympatry, and remain useful groups for predicting morphological, ecological, behavioural and genetic characteristics.

Butterfly genomics eclosing

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

Repeating patterns of mimicry

Axel Meyer discusses how evolution frequently converges on similar mimetic patterns in different species.

Transgenic approaches to study wing color pattern development in Lepidoptera

The extremely diverse lepidopteran wing patterns make useful models to study the evolution of development and the molecular changes that enable it. Until now, the implication of candidate genes in the differentiation of color patterns has relied primarily on correlational evidence, i.e., gene expression patterns in a developing wing mapping closely to the adult color pattern. The use of transgenic techniques in the Lepidoptera, including the manipulation of gene expression, will finally allow researchers to test hypotheses of gene function at various levels of the patterning hierarchy, from signaling ligands and transcription factors to pigment enzymes. Here we present an overview of transgenic techniques employed in lepidopteran systems and highlight areas where current and future research will provide exciting opportunities to deepen our understanding of the mechanisms of morphological evolution.

A conserved supergene locus controls colour pattern diversity in Heliconius butterflies

We studied whether similar developmental genetic mechanisms are involved in both convergent and divergent evolution. Mimetic insects are known for their diversity of patterns as well as their remarkable evolutionary convergence, and they have played an important role in controversies over the respective roles of selection and constraints in adaptive evolution. Here we contrast three butterfly species, all classic examples of Müllerian mimicry. We used a genetic linkage map to show that a locus, Yb, which controls the presence of a yellow band in geographic races of Heliconius melpomene, maps precisely to the same location as the locus Cr, which has very similar phenotypic effects in its co-mimic H. erato. Furthermore, the same genomic location acts as a "supergene", determining multiple sympatric morphs in a third species, H. numata. H. numata is a species with a very different phenotypic appearance, whose many forms mimic different unrelated ithomiine butterflies in the genus Melinaea. Other unlinked colour pattern loci map to a homologous linkage group in the co-mimics H. melpomene and H. erato, but they are not involved in mimetic polymorphism in H. numata. Hence, a single region from the multilocus colour pattern architecture of H. melpomene and H. erato appears to have gained control of the entire wing-pattern variability in H. numata, presumably as a result of selection for mimetic "supergene" polymorphism without intermediates. Although we cannot at this stage confirm the homology of the loci segregating in the three species, our results imply that a conserved yet relatively unconstrained mechanism underlying pattern switching can affect mimicry in radically different ways. We also show that adaptive evolution, both convergent and diversifying, can occur by the repeated involvement of the same genomic regions.

Heliconius wing patterns: an evo-devo model for understanding phenotypic diversity

Evolutionary Developmental Biology aims for a mechanistic understanding of phenotypic diversity, and present knowledge is largely based on gene expression and interaction patterns from a small number of well-known model organisms. However, our understanding of biological diversification depends on our ability to pinpoint the causes of natural variation at a micro-evolutionary level, and therefore requires the isolation of genetic and developmental variation in a controlled genetic background. The colour patterns of Heliconius butterflies (Nymphalidae: Heliconiinae) provide a rich suite of naturally occurring variants with striking phenotypic diversity and multiple taxonomic levels of variation. Diversification in the genus is well known for its dramatic colour-pattern divergence between races or closely related species, and for Müllerian mimicry convergence between distantly related species, providing a unique system to study the development basis of colour-pattern evolution. A long history of genetic studies has showed that pattern variation is based on allelic combinations at a surprisingly small number of loci, and recent developmental evidence suggests that pattern development in Heliconius is different from the eyespot determination of other butterflies. Fine-scale genetic mapping studies have shown that a shared toolkit of genes is used to produce both convergent and divergent phenotypes. These exciting results and the development of new genomic resources make Heliconius a very promising evo-devo model for the study of adaptive change.

Parallel genetic architecture of parallel adaptive radiations in mimetic Heliconius butterflies

It is unknown whether homologous loci underlie the independent and parallel wing pattern radiations of Heliconius butterflies. By comparing the locations of color patterning genes on linkage maps we show that three loci that act similarly in the two radiations are in similar positions on homologous chromosomes.