Published in Evolution & Development 5, issue 3, 269-280, 2003 1 which should be used for any reference to this work
Mimicry: developmental genes that contribute to speciation
Russell E. Naisbit,*,1 Chris D. Jiggins, and James Mallet
The Galton Laboratory, Department of Biology, University College London, London NW1 2HE, UK The Smithsonian Tropical Research Institute, Apartado 2072, Balboa, Panama
*Author for correspondence (e-mail: [email protected]) 1Present address: Laboratoire d’Ecologie Animale et Entomologie, Institut de Zoologie, Université de Neuchâtel, CH-2007 Neuchâtel, Switzerland.
SUMMARY Despite renewed interest in the role of natural overall fitness through ecological selection against nonmi- selection as a catalyst for the origin of species, the develop- metic hybrid genotypes. Most of the loci are clustered into mental and genetic basis of speciation remains poorly under- two genomic regions or “supergenes,” suggesting color pat- stood. Here we describe the genetics of Müllerian mimicry in tern evolution is constrained by preexisting linked elements Heliconius cydno and H. melpomene (Lepidoptera: Nymph- that may have arisen via tandem duplication rather than hav- alidae), sister species that recently diverged to mimic other ing been assembled by natural selection. Linkage, modifi- Heliconius. This mimetic shift was a key step in their specia- ers, and epistasis affect the strength of mimicry as a barrier tion, leading to pre- and postmating isolation. We identify 10 to gene flow between these naturally hybridizing species autosomal loci, half of which have major effects. At least eight and may permit introgression in genomic regions unlinked appear to be homologous with genes known to control pat- to those under disruptive selection. Müllerian mimics in Hel- tern differences within each species. Dominance has evolved iconius use different genetic architectures to achieve the under the influence of identifiable “modifier” loci rather than same mimetic patterns, implying few developmental con- being a fixed characteristic of each locus. Epistasis is found straints. Therefore, although developmental and genomic at many levels: phenotypic interaction between specific pairs constraints undoubtedly influence the evolutionary process, of genes, developmental canalization due to polygenic mod- their effects are probably not strong in comparison with nat- ifiers so that patterns are less sharply defined in hybrids, and ural selection.
INTRODUCTION many factors of small effect and a few of large effect (Orr 1998, 1999). The emerging field of evolutionary developmental biology Second, how do development and selection interact in spe- seeks to explain changes in ontogeny that lead from altered ciation and macroevolution, and can microevolution and genotype to altered phenotype, a process that has been called macroevolution be explained by the same developmental pro- “developmental reprogramming” (Arthur 2000). It will pro- cesses? Most work has focused on major events in body plan vide a fundamental contribution to evolutionary theory if it diversification (Holland 2000; Shankland and Seaver 2000), yields answers to two major questions. with fewer studies of intraspecific or interspecies differences First, to what extent does development constrain or drive (Stern 1998; Beldade et al. 2002). It therefore remains unclear evolution? In other words, are there emergent properties of if there are fundamentally distinct levels of diversification, or developmental reprogramming so that directionality in evo- if the higher levels can be extrapolated from microevolution lution is not the sole preserve of natural selection acting on (Leroi 2000). It has recently become apparent that speciation random variation? Such “developmental drive” would act in is often caused by normal processes of adaptive differentiation conjunction with selection rather than in opposition, leading under divergent natural selection (Filchak et al. 2000; Rundle to beneficial but suboptimal evolution, so that detectable et al. 2000; Jiggins et al. 2001; Podos 2001), so that macroevo- phylogenetic inertia should result (Arthur 2001). Does adap- lution might be simply extrapolated from microevolution. Al- tive divergence typically proceed by the substitution of many though a few studies have identified both the key traits initiat- genes of minor effect (Fisher 1930), or can genes of large ef- ing reproductive isolation and their developmental genetic fect contribute to adaptation and speciation (Orr and Coyne basis (Schemske and Bradshaw 1999; Hawthorne and Via 1992; Coyne and Orr 1998)? Recent theory suggests that 2001; Peichel et al. 2001), the field is still largely dominated when natural selection optimizes quantitative traits, an expo- by studies of hybrid sterility and inviability (Coyne and Orr nential distribution of gene effects will become fixed, with 1998; Orr and Presgraves 2000).
2
3
Genetic dissection of speciation may also provide an- early developmental pathways are redeployed in wing pat- swers to questions such as what roles do linkage, dominance, tern formation at a much later stage (Carroll et al. 1994; Bru- and epistasis play in adaptive speciation? Are these genetic netti et al. 2001; Beldade et al. 2002). constraints incidental properties of newly evolved genes, or do the constraints themselves evolve? For instance, tight Mimicry in Heliconius linkage between functionally related loci may exist because Heliconius butterflies are warningly colored and unpalatable of recent tandem duplication, or it may be adaptive, because and are often Müllerian mimics of other Heliconius or ith- clustering of loci into linked blocks can preserve associa- omiine butterflies (Turner 1984; Sheppard et al. 1985; Mal- tions between coadapted alleles (Noor et al. 2001; Rieseberg let et al. 1998). Extensive work has shown that major genes 2001). Epistasis in particular is fundamental to speciation, control color pattern differences between geographic races because it is the production of ecologically or intrinsically within species (Turner and Crane 1962; Sheppard et al. maladapted gene combinations in hybrids that causes repro- 1985; Mallet 1989; Linares 1996, 1997), but few studies in- ductive isolation (Whitlock et al. 1995; Turelli and Orr volve more than one species (Nijhout 1991; Jiggins and Mc- 2000). Millan 1997). Mimicry genes in Heliconius are developmen- Mimicry in butterflies provides a useful and visually ap- tal or regulatory rather than structural for two reasons. First, pealing example of adaptive evolution in which to answer although changing pattern they do not alter the ability to pro- these questions, where a link can be made between develop- duce pigment. For example, black (melanin), red or brown mental genetics and fitness in the field (Fisher 1930; Gold- (xanthommatins), and yellow (3-hydroxykynurenine) are all schmidt 1945; Turner 1984; Mallet and Barton 1989; Kapan usually present somewhere on the butterfly. Second, as with 2001). Mimicry is also implicated in speciation, because re- other butterfly color pattern genes (Goldschmidt 1945), there productive isolation can arise as an incidental by-product of are correlations between the microstructure of wing scales adaptive divergence in color pattern (Bates 1862; Darwin and the pigments laid down (Gilbert et al. 1988; Janssen et 1863; Vane-Wright 1978; Turner 1981; Mallet et al. 1998). al. 2001). Goldschmidt (1940), who proposed that speciation occurred The genetics of mimicry is usually studied within an eco- via major “systemic mutations” fundamentally different logical genetics (microevolutionary) framework. However, from normal quantitative variation, argued that major genes divergence in mimicry of the butterflies Heliconius cydno affecting mimicry were among the few examples of systemic and H. melpomene (Fig. 1A) has strongly affected mate pref- mutations occurring within species. In 1945, Goldschmidt erences and has led to maladaptive nonmimetic hybrids made the case for mimicry as an example of a wider phenom- (Mallet et al. 1998; Jiggins et al. 2001). The species also dis- enon where “developmental constraints” were as important play female hybrid sterility (Naisbit et al. 2002) and reduced as natural selection in evolution. Despite being morphologi- hybrid mating success (Naisbit et al. 2001), but these almost cally simple two-dimensional traits, butterfly wing color pat- certainly evolved after initial divergence, whereas traits such terns are still poorly understood at the developmental level as differences in microhabitat (Mallet and Gilbert 1995) and (Carroll et al. 1994; Koch et al. 1998; Brunetti et al. 2001; host plant use (Smiley 1978) provide only weak barriers to McMillan et al. 2002). One emerging generalization is that gene flow. Divergence in mimicry was therefore a key step
Fig. 1. (A) Heliconius melpomene, H. cydno, and their nonmimetic F1 hybrid. All pairs of wings are shown at 60% life size, with the upper surface on the right and the lower surface on the left. (B) Interaction between the N and B loci in the forewing band. A dash indicates an allele undetermined because of dominance. (C) Control of the yellow hindwing bar of H. melpomene by the Yb locus. The effect of the H. melpomene allele is shown to the left and that of the H. cydno allele to the right. The white band in the wings on the right is added by the Sb locus. (D) Control of the white hindwing submarginal band of H. cydno by the Sb locus and modification of dominance by J. The H. cydno submarginal band phenotype is shown top left and that of H. melpomene bottom right. The center row shows modification of dominance in Sb heterozygotes by the J locus. (E) Control of the brown forceps-shaped marking of H. cydno by the Br locus. The effect of the H. cydno allele is shown to the left and that of the H. melpomene allele to the right. (F) Far left, control of forewing band color by the K locus. For this and the following two loci, the effect of the H. melpomene allele is shown in the top row, above that of the H. cydno allele. Center left, control of the color of the underside of the red forewing band by the Vf locus. Center right, control of the anterior half of the forewing white hourglass by the Ac locus. Far right, variation in the width of the red portion of the forewing band.
(G) Variation in the width of the white portion of the forewing band. The wing in the center shows the phenotype of an F 1 hybrid. (H) Left, the red line (arrowed) of H. melpomene at the base of the forewing lower surface controlled by the G locus. Center and right, two hybrid phenotypes produced by recombination within the N-Sb-Vf-Yb linkage group. Both are the result of a cross-over between N and N B Sb-Vf-Yb in a backcross of F1 male to H. melpomene female. Center, genotype N N Sb1Sb1Vf2Vf2ybyb, right, genotype B B B B N N Sb1Sb3Vf1Vf2Ybcyb. Shown at 80% life size. A forewing of genotype N N Vf1Vf2 is shown in F, center left lower wings. The only B B other cross-over seen within this linkage group was between Yb and Vf-Sb-N, producing the genotype N N Sb1Sb1Vf2Vf2Ybyb. The hind- wing is shown in D, lower right, and the forewing was like that of H. melpomene.
4
in speciation of H. cydno and H. melpomene, and the devel- vines for oviposition. Eggs were collected daily, and caterpillars opmental genes controlling color pattern differences are as fed on new growth of Passiflora biflora. much “speciation genes” as those for hybrid sterility in Dro- We were able to obtain offspring only from crosses between male sophila (Orr and Presgraves 2000; Ting et al. 2000). H. melpomene and female H. cydno due to strongly asymmetrical mate preferences. Sterility of F females conformed to Haldane’s rule and Here we investigate genes that determine differences in 1 prevented F crosses (Naisbit et al. 2002), so color pattern segregation mimicry between H. cydno and H. melpomene to answer the 2 was examined in backcrosses using fertile F1 males. Single gene con- following questions: trol was inferred when a 1:1 ratio of distinct phenotypes segregated in 1. What is the genetic architecture of mimicry and does it the backcross to the parental species bearing the recessive form of the trait, and unless indicated in the results, all tests are for deviation from suggest developmental constraints? Are individual gene this 1:1 ratio. Where a continuous distribution of intermediate pheno- effects major or minor? What roles do epistasis, linkage, types was produced for any single pattern element, control was judged and dominance play in the evolution of color pattern? polygenic. Homology was inferred if gene effects and linkage were 2. Are mimicry genes that contribute to speciation ho- identical with loci previously described from interracial crosses within mologous to those used in mimetic shifts within each either species. Summaries of individual genotypes are given in Tables species? 1–3, but segregation ratios and recombination frequencies include data 3. Do partners in Müllerian mimicry use homologous ge- from 19 additional individuals (14 from backcrosses to H. cydno and 5 netic variation to achieve the same patterns? If so, de- from backcrosses to H. melpomene) that could not be scored at all loci velopmental constraints could be much more impor- due to wing damage or failure to eclose fully. tant in the evolution of mimicry (Goldschmidt 1945; We investigated the clustering of color pattern loci into tight Nijhout 1991) than classically believed. linkage groups by comparing the extent of linkage among the 10 color pattern loci with a null model assuming random distribution of loci across chromosomes to perform a test similar to that in Turner (1984, p. 158). The null distribution of loci per chromo- MATERIALS AND METHODS some is approximately but not exactly Poisson distributed, where the dispersion variance/mean 1. Clustering was tested here Crosses were performed in Gamboa, Panama between September numerically by assigning 10 loci randomly onto 21 chromosomes 1999 and March 2000 using Heliconius cydno chioneus and H. mel- one million times. We used as a test statistic the log-likelihood pomene rosina collected from nearby forest in Soberanía National ratio (G), where 10/21 loci are expected on average per chromo- Park. To obtain crosses, we isolated virgin females with older some. A more conservative test was also run using eight loci, be- males. After mating, females were kept individually in 1 1 2-m cause two pairs of putative loci that did not recombine (B and G, outdoor insectaries and supplied with pollen sources (Lantana and Sb and Vf) might each represent pleiotropic effects of a single Psiguria), artificial nectar (10% sugar solution), and Passiflora gene.