The genetic basis of an adaptive radiation: warning colour in two Heliconius species

Mimetic colour pattern races of Heliconius butterflies provide a striking example of adaptive radiation and numerous crossing experiments have investigated the genetics of these racial differences. However, colour pattern differentiation between closely related Heliconius species has not been previously studied. Here we present data from crosses between H. erato cyrbia and its sister species, H. himera. The genetic architecture underlying colour pattern divergence between these species is identical to that observed between races of H. erato. As in inter-racial crosses, colour pattern differences resulted from segregation at a few major loci. Evidence from 1321 offspring in 4 F₁, 17 backcross, 7 F₂ and 21 further crosses showed that two major loci controlled most of the colour pattern differences between H. erato and H. himera. There were strong interactions between these loci in their patterns of expression and evidence for other loci with relatively minor phenotypic effects. More importantly, based on patterns of expression within broods and linkage with Aconitase, we conclude that these major loci were homologous with those known to be responsible for colour pattern differences within H. erato. Our crosses also permit a re-evaluation of the relationships between colour pattern races of H. erato. This suggests that H. e. hydara, which occurs across a major mtDNA break, is the ancestral phenotype from which other races have evolved. Based on this assumption, we find no evidence to support the recent suggestion that apparently homologous colour pattern alleles have arisen multiple times.     

Quantitative genetics of butterfly wing color patterns

Developmental processes exert their influence on the evolution of complex morphologies through the genetic correlations they engender between traits. Butterfly wing color patterns provide a model system to examine this connection between development and evolution. In butterflies, the nymphalid groundplan is a framework used to decompose complex wing patterns into their component pattern elements. The first goal of this work has been to determine whether the components of the nymphalid groundplan are the products of independent developmental processes. To test this hypothesis, the genetic correlation matrices for two species of butterflies, Precis coenia and Precis evarete, were estimated for 27 wing pattern characters. The second purpose was to test the hypothesis that the differentiation of serial homologs lowers their genetic correlations. The “eyespots” found serially repeated across the fore- and hindwing and on the dorsal and ventral wing surfaces provided an opportunity to test this hypothesis. The genetic correlation matrices of both species were very similar. The pattern of genetic correlation measured between the different types of pattern elements and between the homologous repeats of a pattern element supported the first hypothesis of developmental independence among the elements of the groundplan. The correlation pattern among the differentiated serial homologs was similarly found to support the second hypothesis: pairs of eyespots that had differentiated had lower genetic correlations than pairs that were similar in morphology. The implications of this study are twofold: First, the apparent developmental independence among the distinct elements of wing pattern has facilitated the vast diversification in morphology found in butterflies. Second, the lower genetic correlations betweendifferentiated homologs demonstrates that developmental constraints can in fact be broken. The extent to which genetic correlations readily change, however, remains unknown. © 1994 Wiley-Liss, Inc.     

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.     

Homologies in the colour patterns of the genus Heliconius (Lepidoptera: Nymphalidae)

The colour patterns of Heliconius butterflies are composed from a relatively simple set of pattern elements whose homologues are recognizable throughout the genus. Although Heliconius colour patterns look quite different from those of most nymphalids, these pattern elements are seen to derive from the generalized nymphalid groundplan. The differences arise primarily from the loss or positional shift of certain pattern elements, a high degree of fusion between individual pattern elements, and, in the forewing, asymmetries of the pattern elements relative to the wing-cell midline. The scheme of homologies we present is consistent with what is currently known about the comparative morphology and developmental physiology of colour pattern formation in Lepidoptera, and provides a framework for the interpretation of developmental, evolutionary and genetic studies in Heliconius.     

Genetic evidence for a sibling species of Heliconius charithonia (Lepidoptera; Nymphalidae)

Heliconius charithonia is a widespread species which, unlike many Heliconius, is non-mimetic and shows little racial differentiation. Only one form, 'peruvianus', which occurs in the dry forest habitats of western Ecuador and Peru, has a distinct and clearly mimetic colour pattern. Here it was shown that H. peruvianus was distinct from H. charithonia bassleri at allozyme loci (D = 0.25 over 22 loci). This differentiation was ten times greater than that between H. charithonia sampled from Ecuador and the Caribbean (D = 0.027) and was consistent with analysis of mitochondrial sequence data (3.4-4% sequence divergence between H. peruvianus and H. charithonia). One individual with a H. peruvianus colour pattern and allozyme genotype was collected in an area where H. charithonia was known to be common, demonstrating that contact between the taxa occurs in western Ecuador. Furthermore, the allozyme genotype of another individual was heterozygous for four of five diagnostic loci and was most likely an F1 hybrid between H. charithonia and H. peruvianus. These data imply that H. charithonia and H. peruvianus are distinct species which hybridize occasionally. This species pair show many similarities with H. erato and H. himera, which are similarly differentiated genetically and also show ecological and colour pattern differences. These species fulfil some of the predictions of both allopatric refugium and parapatric adaptationist models of speciation in the neotropics, suggesting that elements of both hypotheses may be true.     

Hybridization and transgressive exploration of colour pattern and wing morphology in Heliconius butterflies

Hybridization can generate novel phenotypes distinct from those of parental lineages, a phenomenon known as transgressive trait variation. Transgressive phenotypes might negatively or positively affect hybrid fitness, and increase available variation. Closely related species of Heliconius butterflies regularly produce hybrids in nature, and hybridization is thought to play a role in the diversification of novel wing colour patterns despite strong stabilizing selection due to interspecific mimicry. Here, we studied wing phenotypes in first‐ and second‐generation hybrids produced by controlled crosses between either two co‐mimetic species of Heliconius or between two nonmimetic species. We quantified wing size, shape and colour pattern variation and asked whether hybrids displayed transgressive wing phenotypes. Discrete traits underlain by major‐effect loci, such as the presence or absence of colour patches, generate novel phenotypes. For quantitative traits, such as wing shape or subtle colour pattern characters, hybrids only exceed the parental range in specific dimensions of the morphological space. Overall, our study addresses some of the challenges in defining and measuring phenotypic transgression for multivariate traits and our data suggest that the extent to which transgressive trait variation in hybrids contributes to phenotypic diversity depends on the complexity and the genetic architecture of the traits.     

Mimicry: developmental genes that contribute to speciation

Despite renewed interest in the role of natural selection as a catalyst for the origin of species, the developmental and genetic basis of speciation remains poorly understood. Here we describe the genetics of Müllerian mimicry in Heliconius cydno and H. melpomene (Lepidoptera: Nymphalidae), sister species that recently diverged to mimic other Heliconius. This mimetic shift was a key step in their speciation, leading to pre- and postmating isolation. We identify 10 autosomal loci, half of which have major effects. At least eight appear to be homologous with genes known to control pattern differences within each species. Dominance has evolved under the influence of identifiable "modifier" loci rather than being a fixed characteristic of each locus. Epistasis is found at many levels: phenotypic interaction between specific pairs of genes, developmental canalization due to polygenic modifiers so that patterns are less sharply defined in hybrids, and overall fitness through ecological selection against nonmimetic hybrid genotypes. Most of the loci are clustered into two genomic regions or "supergenes," suggesting color pattern evolution is constrained by preexisting linked elements that may have arisen via tandem duplication rather than having been assembled by natural selection. Linkage, modifiers, and epistasis affect the strength of mimicry as a barrier to gene flow between these naturally hybridizing species and may permit introgression in genomic regions unlinked to those under disruptive selection. Müllerian mimics in Heliconius use different genetic architectures to achieve the same mimetic patterns, implying few developmental constraints. Therefore, although developmental and genomic constraints undoubtedly influence the evolutionary process, their effects are probably not strong in comparison with natural selection.     

PHENOTYPIC CORRELATION STRUCTURE AMONG ELEMENTS OF THE COLOR PATTERN IN PRECIS COENIA (LEPIDOPTERA: NYMPHALIDAE)

The component parts of butterfly wing patterns are arranged in sets of serially homologous pattern elements, repeated from wing cell to wing cell. Measurements were made on the sizes and positions of these elements on two successive, independent, sets of specimens in order to elucidate the phenotypic correlation structure among pattern elements. That portion of the correlation between measures due to overall size variation was accounted for through two alternate methods: multiple regression on two vein length measures, which represent wing size, and a Wright-style factor analysis. The sizes of pattern elements belonging to a homologous series were found to be significantly correlated whereas those of non-homologous elements varied independently. The degree of correlation among homologs varied, and, in the case of eyespot sizes, appeared to be inversely related to the degree of their morphological divergence. Although not correlated in size, the positions of non-homologous elements that lie within the same wing cell are moderately correlated. The results support current developmental models for the ontogeny of butterfly color pattern.     

What can hybrid zones tell us about speciation? The case of Heliconius erato and H. himera (Lepidoptera: Nymphalidae)

To understand speciation we need to study the genetics and ecology of intermediate cases where interspecific hybridization still occurs. Two closely related species of Heliconius butterflies meet this criterion: Heliconius himera is endemic to dry forest and thorn scrub in southern Ecuador and northern Peru, while its sister species, H. erato , is ubiquitous in wet forest throughout south and central America. In three known zones of contact, the two species remain distinct, while hybrids are found at low frequency. Collections in southern Ecuador show that the contact zone is about 5 km wide, half the width of the narrowest clines between colour pattern races of H. erato. The narrowness of this dine argues that very strong selection (s ≅ 1) is maintaining the parapatric distributions of these two species. The zone is closely related with a habitat transition from wet to dry forest, which suggests that the narrow zone of parapatry is maintained primarily by ecological adaptation. Selection on colour pattern loci, assortative mating and hybrid inviability may also be important. The genetics of hybrids between the two species shows that the major gene control of pattern elements is similar to that found in previous studies of H. erato races, and some of the loci are homologous. This suggests that similar genetic processes are involved in the morphological divergence of species and races. Evidence from related Heliconius supports a hypothesis that ecological adaptation is the driving force for speciation in the group.     

A hybrid zone provides evidence for incipient ecological speciation in Heliconius butterflies

In Heliconius butterflies, it has been proposed that speciation occurs through a combination of divergence in ecological habitat preferences and mimetic colour patterns. Here we test this hypothesis by investigating a parapatric form of the widespread species Heliconius erato. Mendelian (colour patterns) and molecular genetic data permit us to address hypotheses about introgression and genetic differentiation between different populations. Combined analysis of colour pattern, microsatellite loci and mitochondrial DNA showed that Heliconius erato venus and Heliconius erato chestertonii form a bimodal hybrid zone implying partial reproductive isolation. In a sample of 121 individuals collected in sympatry, 25% were hybrids representing a significant deficit of heterozygotes compared to the Hardy-Weinberg expectation. Seven microsatellite loci, analysed for a subset of these individuals, showed marked differentiation between the parental taxa, and unambiguously identified two genotypic clusters concordant with our phenotypic classification of individuals. Mitochondrial DNA analysis showed H. erato venus as a monophyletic group well differentiated from H. erato chestertonii, implying a lack of historical introgression between the populations. Heliconius erato chestertonii is therefore an incipient species that maintains its integrity despite high levels of hybridization. Moreover, H. erato chestertonii is found at higher altitudes than other races of H. erato and has a distinct colour pattern and mimetic relationship. Hence, there are now two examples of parapatric incipient species related to H. erato, H. himera and H. erato chestertonii, both of which are associated with higher altitudes, more arid habitats and distinct mimetic relationships. This implies that parapatric habitat adaptation is a likely cause of speciation in this group.     

Color pattern analysis of nymphalid butterfly wings: revision of the nymphalid groundplan

To better understand the developmental mechanisms of color pattern variation in butterfly wings, it is important to construct an accurate representation of pattern elements, known as the "nymphalid groundplan". However, some aspects of the current groundplan remain elusive. Here, I examined wing-wide elemental patterns of various nymphalid butterflies and confirmed that wing-wide color patterns are composed of the border, central, and basal symmetry systems. The central and basal symmetry systems can express circular patterns resembling eyespots, indicating that these systems have developmental mechanisms similar to those of the border symmetry system. The wing root band commonly occurs as a distinct symmetry system independent from the basal symmetry system. In addition, the marginal and submarginal bands are likely generated as a single system, referred to as the "marginal band system". Background spaces between two symmetry systems are sometimes light in coloration and can produce white bands, contributing significantly to color pattern diversity. When an element is enlarged with a pale central area, a visually similar (yet developmentally distinct) white band is produced. Based on the symmetric relationships of elements, I propose that both the central and border symmetry systems are comprised of "core elements" (the discal spot and the border ocelli, respectively) and a pair of "paracore elements" (the distal and proximal bands and the parafocal elements, respectively). Both core and paracore elements can be doubled, or outlined. Developmentally, this system configuration is consistent with the induction model, but not with the concentration gradient model for positional information.     

Within-population variation in the chemistry of life: the stoichiometry of sexual dimorphism in multiple dimensions

In sexual species, phenotypic divergence between males and females, or sexual dimorphism, is often the source of the most staggering examples of phenotypic variation in nature. Theory suggests that exaggerated sexual traits should drive sex-specific nutritional demands. Advances in spectrometry enable rapid quantification of the elements that make up individuals and traits, which can be used to assess patterns of intraspecific variation and the contribution of nutritionally-demanding sexual traits to these patterns. We measured dimorphism in the whole body stoichiometry of Hyalella amphipods and examined whether nutritional demands of exaggerated sexual traits differ from those of similar traits not under sexual selection. We found striking sexual dimorphism in multivariate whole body elemental composition (i.e., the ionome), including elements important for organismal growth and performance. In males, the exaggerated, sexually-selected claw-like appendage (posterior gnathopod) differed significantly in mass-specific stoichiometry from a similarly sized and serially homologous non-sexual trait (fifth pereopod), indicating that there are fundamental differences in the construction of sexual traits in relation to similar traits that are not under sexual selection. While sexually selected traits do differ from non-sexual traits in their ionomes, we found that possessing an exaggerated trait does not change organismal stoichiometry, indicating that trait exaggeration may not be directly driving ionomic sexual dimorphism. Finally, we found that larger traits are not comparatively larger resource sinks for any element, suggesting that the possession of larger traits is not a function of greater allocation of resources. Together, we discovered substantial sexual dimorphism at the lowest level of organization, chemical elements. Such information illuminates predictions about dimorphisms in foraging behavior, nutritional physiology, and sex-specific selection on the underlying loci. High throughput, multidimensional data on sexual divergence in stoichiometric composition is a powerful tool in understanding the evolutionary ecology of sexual dimorphisms.     

Mutants highlight the modular control of butterfly eyespot patterns

SUMMARY The eyespots on butterfly wings are thought to be serially homologous pattern elements. Yet eyespots differ greatly in number, shape, color, and size, within and among species. To what extent do these serially homologues have separate developmental identities, upon which selection acts to create diversity? We examined x‐ray–induced mutations for the eyespots of the nymphalid butterfly Bicyclus anynana that highlight the modular control of these serially homologous wing pattern elements. These mutations reduce or eliminate individual eyespots, or groups of eyespots, with no further effect on the wing color pattern. The collection of mutants highlights a greater potential developmental repertoire than that observed across the genus Bicyclus. We studied in detail one such mutation, of codominant effect, that causes the elimination of two adjacent eyespots on the ventral hindwing. By analyzing the expression of genes known to be involved in eyespot formation, we found an alteration in the differentiation of the “organizing” cells at the eyespot's center. No such cells differentiate in the wing subdivisions lacking the two eyespots in the mutants. We propose several developmental models, based on wing compartmentalization in Drosophila, that provide the first framework for thinking about the molecular evolution of butterfly wing pattern modularity.     

Genetic variation in two conserved local Romanian pig breeds using type 1 DNA markers

Analysis of the genetic variation of an endangered population is an important component for the success of conservation. Animals from two local Romanian pig breeds, the Mangalitsa and Bazna, were analysed for variation at a number of genetic loci using PCR-based DNA tests. Polymorphism was assessed at loci which 1) are known to cause phenotypic variation, 2) are potentially involved in trait differences or 3) are putative candidate genes. The traits considered are disease resistance, growth, coat colour, meat quality and prolificacy. Even though the populations are small and the markers are limited to specific genes, we found significant differences in five of the ten characterised loci. In some cases the observed allele frequencies were interesting in relation to gene function and the phenotype of the breed. These breeds are part of a conservation programme in Romania and marker information may be useful in preserving a representative gene pool in the populations. The use of polymorphisms in type 1 (gene) markers may be a useful complement to analysis based on anonymous markers.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

Genetic constraints on wing pattern variation in Lycaeides butterflies: A case study on mapping complex,multifaceted traits in structured populations

Patterns of phenotypic variation within and among species can be shaped and constrained by trait genetic architecture. This is particularly true for complex traits, such as butterfly wing patterns, that consist of multiple elements. Understanding the genetics of complex trait variation across species boundaries is difficult, as it necessitates mapping in structured populations and can involve many loci with small or variable phenotypic effects. Here, we investigate the genetic architecture of complex wing pattern variation in Lycaeides butterflies as a case study of mapping multivariate traits in wild populations that include multiple nominal species or groups. We identify conserved modules of integrated wing pattern elements within populations and species. We show that trait covariances within modules have a genetic basis and thus represent genetic constraints that can channel evolution. Consistent with this, we find evidence that evolutionary changes in wing patterns among populations and species occur in the directions of genetic covariances within these groups. Thus, we show that genetic constraints affect patterns of biological diversity (wing pattern) in Lycaeides, and we provide an analytical template for similar work in other systems.     

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced‐complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over‐represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini ‐ the largest radiation of mimetic butterflies.     

Cytogenetic Analysis of the Second Chromosome Heterochromatin of Drosophila Melanogaster

This paper reports the cytogenetic characterization of the second chromosome heterochromatin of Drosophila melanogaster. High resolution cytological analysis of a sample of translocations, inversions, deficiencies and free duplications involving the pericentric regions of the second chromosome was achieved by applying sequential Hoechst 33258 and N-chromosome banding techniques to larval neuroblast prometaphase chromosomes. Heterochromatic rearrangements were employed in a series of complementation assays and the genetic elements previously reported to be within or near the second chromosome heterochromatin were thus precisely assigned to specific heterochromatic bands. The results of this analysis reveal a nonhomogeneous distribution of loci along the second chromosome heterochromatin. The l(2)41Aa, l(2)41Ab, rolled (l(2)41Ac) and l(2)41Ad loci are located within the proximal heterochromatin of 2R, while the nine remaining loci in the left arm and two (l(2)41Ae and l(2)41Ah) in the right arm map to h35 and to h46, respectively, the most distal heterochromatic regions. In addition, a common feature of these loci revealed by the cytogenetic analysis is that they map to specific heterochromatic blocks but do not correspond to the blocks themselves, suggesting that they are not as large as the Y fertility factors or the Rsp locus. Mutations of the proximal most heterochromatic loci, l(2)41Aa and rolled, were also examined for their phenotypic effects. Extensive cell death during imaginal disc development was observed in individuals hemizygous for either the EMS 31 and rolled mutations, leading to a pattern of phenotypic defects of adult structures.