Search

International audience; The biodiversity crisis has highlighted the need to assess and map biodiversity in order to prioritize conservation efforts. Clearwing butterflies (tribe Ithomiini) have been proposed as biological indicators for habitat quality in Neotropical forests, which contain the world's richest biological communities. Here, we provide maps of different facets of Ithomiini diversity across the Neotropics to identify areas of evolutionary and ecological importance for conservation and evaluate their overlap with current anthropogenic threats.Location: Neotropics.Methods: We ran species distribution models on a data set based on 28,986 georeferenced occurrences representing 388 ithomiine species to generate maps of geographic rarity, taxonomic, phylogenetic and Müllerian mimetic wing pattern diversity. We quantified and mapped the overlap of diversity hotspots with areas threatened by or providing refuge from current anthropogenic pressures.Results: The eastern slopes of the Andes formed the primary hotspot of taxonomic, phylogenetic and mimetic diversity, with secondary hotspots in Central America and the Atlantic Forest. Most diversity indices were strongly spatially correlated. Nevertheless, species-poor communities on the Pacific slopes of the Andes also sheltered some of the geographically rarest species. Overall, tropical montane forests that host high species and mimetic diversity as well as rare species and mimicry rings appeared particularly under threat.Main conclusions: Remote parts of the Upper Amazon may act as refuges against current anthropogenic pressures for a limited portion of Ithomiini diversity. Furthermore, it is likely that the current threat status may worsen with ongoing climate change and deforestation. In this context, the tropical Andes occupy a crucial position as the primary hotspot for multiple facets of biodiversity for ithomiine butterflies, as they do for angiosperms, tetrapods and other insect taxa. Our results support the role of ithomiine butterflies as a suitable flagship indicator group for Neotropical butterfly diversity and reinforce the position of the tropical Andes as a flagship region for biodiversity conservation in general, and insect and butterfly conservation in particular.

During courtship, male butterflies of many species produce androconial secretions containing male sex pheromones (MSPs) that communicate species identity and affect female choice. MSPs are thus likely candidates as reproductive barriers, yet their role in speciation remains poorly studied. Although

Aim Climatic changes throughout the Pleistocene have strongly modified species distributions. We examine how these range shifts have affected the genetic diversity of a montane butterfly species and whether the genetic diversity in the extant populations is threatened by future climate change. Location Europe. Taxon Erebia epiphron Lepidoptera: Nymphalidae. Methods We analyzed mtDNA to map current genetic diversity and differentiation of E. epiphron across Europe to identify population refugia and postglacial range shifts. We used species distribution modeling (SDM) to hindcast distributions over the last 21,000 years to identify source locations of extant populations and to project distributions into the future (2070) to predict potential losses in genetic diversity. Results We found substantial genetic diversity unique to specific regions within Europe (total number of haplotypes = 31, number of unique haplotypes = 27, H d = 0.9). Genetic data and SDM hindcasting suggest long‐term separation and survival of discrete populations. Particularly, high rates of unique diversity in postglacially colonized sites in England (H d = 0.64) suggest this population was colonized from a now extinct cryptic refugium. Under future climate change, SDMs predict loss of climate suitability for E. epiphron, particularly at lower elevations (, The genetic diversification of cold‐adapted mountain species, such as E. epiphron, has been shaped by Pleistocene glaciations, resulting in unique genetic diversity in isolated populations. The unique genetic diversity in mountain and cold‐adapted species is at under future climate warming, and we predict E. epiphron will lose 38%–64% of its range in the future, resulting in the loss of genetic diversity, reducing its ability to adapt.

Much debate surrounds the importance of top-down and bottom-up effects in the Southern Ocean, where the harvesting of over two million whales in the mid twentieth century is thought to have produced a massive surplus of Antarctic krill. This excess of krill may have allowed populations of other predators, such as seals and penguins, to increase, a top-down hypothesis known as the ‘krill surplus hypothesis’. However, a lack of pre-whaling population baselines has made it challenging to investigate historical changes in the abundance of the major krill predators in relation to whaling. Therefore, we used reduced representation sequencing and a coalescent-based maximum composite likelihood approach to reconstruct the recent demographic history of the Antarctic fur seal, a pinniped that was hunted to the brink of extinction by 18th and 19th century sealers. In line with the known history of this species, we found support for a demographic model that included a substantial reduction in population size around the time period of sealing. Furthermore, maximum likelihood estimates from this model suggest that the recovered, post-sealing population at South Georgia may have been around two times larger than the pre-sealing population. Our findings lend support to the krill surplus hypothesis and illustrate the potential of genomic approaches to shed light on long-standing questions in population biology.

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

Hybrids between diverging populations are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex – a pattern known as Haldane’s rule. The genetics of Haldane’s Rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies ofHeliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females, and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary, and as a function of oocyte development. The Z chromosome has a major effect, similar to the “large X effect” inDrosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. Among loci differentially expressed between females with arrested vs. non-arrested ovary development, we identified six candidate genes known also fromDrosophila melanogasterandParage aegeriaoogenesis. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera. We demonstrate that sterility is produced by multiple complex epistastic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane’s Rule.

Phenotypic convergence between taxa can be caused by divergent genetic evolution (different genetic pathways), parallel genetic evolution (convergent mutations), or collateral evolution (shared ancestry). Heliconius butterflies have bright mimetic color patterns shared between multiple species, making an excellent ....., Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex. While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato. Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

In a series of dawn-to-dusk studies, we examined the nature and accessibility of nectar rewards for pollinating insects by monitoring insect visits and the secretion rate and standing crop of nectar in the British native plant species Salvia pratensis , Stachys palustris , S. officinalis , Lythrum salicaria , Linaria vulgaris , the non-native Calendula officinalis , Petunia × hybrida , Salvia splendens , and the possibly introduced Saponaria officinalis . We also compared single with double variants of Lotus corniculatus , Saponaria officinalis , Petunia × hybrida and Calendula officinalis . All the British species studied are nectar-rich and are recommended for pollinator-friendly gardens. They showed maximal secretion rates of about 10–90 μg sugar per flower h −1 , and most had mean standing crops of about 5–60 μg sugar per flower. In all British species studied, the corolla was deep enough for the relatively long-tongued bumblebee Bombus pascuorum , but the shallower flowers of Lythrum salicaria were also much visited by shorter-tongued bees and hoverflies, as well as by butterflies. The exotic Salvia splendens , presumably coevolved with hummingbirds in the Neotropics, has such deep flowers that British bees cannot reach the nectar except by crawling down the corolla. With a secretion rate approaching 300 μg sugar per flower h −1 and little depletion by insects, S. splendens accumulated high standing crops of nectar. S. splendens , and single and double flowers of the two probably moth-pollinated species Petunia × hybrida and Saponaria officinalis , received few daytime visits despite abundant nectar but Calendula was well visited by hoverflies and bees. We compared single and double variants of Lotus corniculatus , Petunia × hybrida and Calendula officinalis , and also Saponaria officinalis , the last being probably introduced in Britain (Stace, 1997 New flora of the British Isles. 2nd edn. Cambridge: Cambridge University Press). In Petunia , Saponaria and Lotus , double flowers secreted little or no nectar. In Calendula , where doubling involved a change in the proportion of disc and ray florets rather than modification of individual flower structure, double and single capitula had similar standing crops of nectar. Except in Calendula , exotic or double flowers were little exploited by insect visitors. In the exotics, this was probably due to the absence or scarcity of coevolved pollinators, coupled, in double flowers, with the absence of nectar. Copyright 2001 Annals of Botany Company

Phylogeny confirming the Alcolapia sequences discussed are indeed CPS III

Table of primers for sequencing CPS III from Alcolapia species

Geographic isolation is suggested to be among the most important processes in the generation of cichlid fish diversity in East Africa's Great Lakes, both through isolation by distance and fluctuating connectivity caused by changing lake levels. However, even broad scale phylogeographic patterns are currently unknown in many non‐cichlid littoral taxa from these systems. To begin to address this, we generated restriction‐site‐associated DNA sequence (RADseq) data to investigate phylogeographic structure throughout Lake Tanganyika (LT) in two broadly sympatric rocky shore catfish species from independent evolutionary radiations with differing behaviors: the mouthbrooding claroteine, Lophiobagrus cyclurus, and the brood‐parasite mochokid, Synodontis multipunctatus. Our results indicated contrasting patterns between these species, with strong lake‐wide phylogeographic signal in L. cyclurus including a deep divergence between the northern and southern lake basins. Further structuring of these populations was observed across a heterogeneous habitat over much smaller distances. Strong population growth was observed in L. cyclurus sampled from shallow shorelines, suggesting population growth associated with the colonization of new habitats following lake‐level rises. Conversely, S. multipunctatus, which occupies a broader depth range, showed little phylogeographic structure and lower rates of population growth. Our findings suggest that isolation by distance and/or habitat barriers may play a role in the divergence of non‐cichlid fishes in LT, but this effect varies by species.

We here pioneer a low-cost assembly strategy for 20 Heliconiini genomes to characterize the evolutionary history of the rapidly radiating genusHeliconius. A bifurcating tree provides a poor fit to the data, and we therefore explore a reticulate phylogeny forHeliconius. We probe the genomic architecture of gene flow, and develop a new method to distinguish incomplete lineage sorting from introgression. We find that most loci with non-canonical histories arose through introgression, and are strongly underrepresented in regions of low recombination and high gene density. This is expected if introgressed alleles are more likely to be purged in such regions due to tighter linkage with incompatibility loci. Finally, we identify a hitherto unrecognized inversion, and show it is a convergent structural rearrangement that captures a known color pattern switch locus within the genus. Our multi-genome assembly approach enables an improved understanding of adaptive radiation.

Identifying the traits causing reproductive isolation and the order in which they evolve is fundamental to understanding speciation. Here, we quantify prezygotic and intrinsic postzygotic isolation among allopatric, parapatric, and sympatric populations of the butterflies Heliconius elevatus and Heliconius pardalinus. Sympatric populations from the Amazon (H. elevatus and H. p. butleri) exhibit strong prezygotic isolation and rarely mate in captivity; however, hybrids are fertile. Allopatric populations from the Amazon (H. p. butleri) and Andes (H. p. sergestus) mate freely when brought together in captivity, but the female F1 hybrids are sterile. Parapatric populations (H. elevatus and H. p. sergestus) exhibit both assortative mating and sterility of female F1s. Assortative mating in sympatric populations is consistent with reinforcement in the face of gene flow, where the driving force, selection against hybrids, is due to disruption of mimicry and other ecological traits rather than hybrid sterility. In contrast, the lack of assortative mating and hybrid sterility observed in allopatric populations suggests that geographic isolation enables the evolution of intrinsic postzygotic reproductive isolation. Our results show how the types of reproductive barriers that evolve between species may depend on geography.

Tetrapods and fish have adapted distinct carbamoyl-phosphate synthase (CPS) enzymes to initiate the ornithine urea cycle during the detoxification of nitrogenous wastes. We report evidence that in the ureotelic subgenus of extremophile fishOreochromis Alcolapia, CPS III has undergone convergent evolution and adapted its substrate affinity to ammonia, which is typical of terrestrial vertebrate CPS I. Unusually, unlike in other vertebrates, the expression of CPS III inAlcolapiais localized to the skeletal muscle and is activated in the myogenic lineage during early embryonic development with expression remaining in mature fish. We propose that adaptation inAlcolapiaincluded both convergent evolution of CPS function to that of terrestrial vertebrates, as well as changes in development mechanisms redirectingCPS IIIgene expression to the skeletal muscle.

Recent advances in high throughput sequencing have transformed the study of wild organisms by facilitating the generation of high quality genome assemblies and dense genetic marker datasets. These resources have the potential to significantly advance our understanding of diverse phenomena at the level of species, populations and individuals, ranging from patterns of synteny through rates of linkage disequilibrium (LD) decay and population structure to individual inbreeding. Consequently, we used PacBio sequencing to refine an existing Antarctic fur seal (Arctocephalus gazella) genome assembly and genotyped 83 individuals from six populations using restriction site associated DNA (RAD) sequencing. The resulting hybrid genome comprised 6,169 scaffolds with an N50 of 6.21 Mb and provided clear evidence for the conservation of large chromosomal segments between the fur seal and dog (Canis lupus familiaris). Focusing on the most extensively sampled population of South Georgia, we found that LD decayed rapidly, reaching the background level ofr2= 0.09 by around 26 kb, consistent with other vertebrates but at odds with the notion that fur seals experienced a strong historical bottleneck. We also found evidence for population structuring, with four main Antarctic island groups being resolved. Finally, appreciable variance in individual inbreeding could be detected, reflecting the strong polygyny and site fidelity of the species. Overall, our study contributes important resources for future genomic studies of fur seals and other pinnipeds while also providing a clear example of how high throughput sequencing can generate diverse biological insights at multiple levels of organisation.

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.

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (FST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.

Neotropical Heliconius butterflies are members of various mimicry rings characterized by diverse colour patterns. In the present study we investigated whether a similar diversity is observed in the chemistry of volatile compounds present in male wing androconia. Recent research has shown that these androconia are used during courting of females. Three to five wild-caught male Heliconius individuals of 17 species and subspecies were analyzed by GC/MS. Most of the identified compounds originate from common fatty acids precursors, including aldehydes, alcohols, acetates or esters preferentially with a C18 and C20 chain, together with some alkanes. The compounds occurred in species-specific mixtures or signatures. For example, octadecanal is characteristic for H. melpomene, but variation in composition between the individuals was observed. Cluster analysis of compound occurrence in individual bouquets and analyses based on biosynthetic motifs such as functional group, chain length, or basic carbon-backbone modification were used to reveal structural patterns. Mimetic pairs contain different scent bouquets, but also some compounds in common, whereas sympatric species, both mimetic and non-mimetic, have more distinct compound compositions. The compounds identified here may play a role in mate choice thus helping maintain species integrity in a butterfly genus characterized by pervasive interspecific gene flow. © 2017, Springer Science+Business Media, LLC.

Multilocus clines between Mullerian mimetic races of Heliconius butterflies provide a classic example of the maintenance of hybrid zones and their importance in speciation. Concordant hybrid zones in the mimics Heliconius erato and H. melpomene in northern Peru were carefully documented in the 1980s, and this prior work now permits a historical analysis of the movement or stasis of the zones. Previous work predicted that these zones might be moving toward the Andes due to selective asymmetry. Extensive deforestation and climate change might also be expected to affect the positions and widths of the hybrid zones. We show that the positions and shapes of these hybrid zones have instead remained remarkably stable between 1985 and 2012. The stability of this interaction strongly implicates continued selection, rather than neutral mixing following secondary contact. The stability of cline widths and strong linkage disequilibria (gametic correlation coefficients Rmax = 0.35–0.56 among unlinked loci) over 25 years suggest that mimetic selection pressures on each color pattern locus have remained approximately constant (s ≈ 0.13–0.40 per locus in both species). Exceptionally high levels of precipitation at the edge of the easternmost Andes may act as a population density trough for butterflies, trapping the hybrid zones at the foot of the mountains, and preventing movement. As such, our results falsify one prediction of the Pleistocene Refugium theory: That the ranges of divergent species or subspecies should be centered on regions characterized by maxima of rainfall, with hybrid zones falling in more arid regions between them.

Most speciation events probably occur gradually, without complete and immediate reproductive isolation, but the full extent of gene flow between diverging species has rarely been characterized on a genome-wide scale. Documenting the extent and timing of admixture between diverging species can clarify the role of geographic isolation in speciation. Here we use new methodology to quantify admixture at different stages of divergence in Heliconius butterflies, based on whole-genome sequences of 31 individuals. Comparisons between sympatric and allopatric populations of H. melpomene, H. cydno, and H. timareta revealed a genome-wide trend of increased shared variation in sympatry, indicative of pervasive interspecific gene flow. Up to 40% of 100-kb genomic windows clustered by geography rather than by species, demonstrating that a very substantial fraction of the genome has been shared between sympatric species. Analyses of genetic variation shared over different time intervals suggested that admixture between these species has continued since early in speciation. Alleles shared between species during recent time intervals displayed higher levels of linkage disequilibrium than those shared over longer time intervals, suggesting that this admixture took place at multiple points during divergence and is probably ongoing. The signal of admixture was significantly reduced around loci controlling divergent wing patterns, as well as throughout the Z chromosome, consistent with strong selection for Müllerian mimicry and with known Z-linked hybrid incompatibility. Overall these results show that species divergence can occur in the face of persistent and genome-wide admixture over long periods of time.

Butterflies in the genus Heliconius have undergone rapid adaptive radiation for warning patterns and mimicry, and are excellent models to study the mechanisms underlying diversification. In Heliconius, mimicry rings typically involve distantly related species, whereas closely related species often join different mimicry rings. Genetic and behavioural studies have n how reproductive isolation in many pairs of Heliconius taxa is largely mediated by natural and sexual selection on wing colour patterns. However, recent studies have uncovered new cases in which pairs of closely related species are near-perfect mimics of each other. Here, we provide morphometric and genetic evidence for the coexistence of two closely related, hybridizing co-mimetic species on the eastern slopes of the Andes, H. melpomene amaryllis and H. timareta ssp. nov., which is described here as H. timareta thelxinoe. A joint analysis of multilocus genotyping and geometric morphometrics of wing shape shows a high level of differentiation between the two species, with only limited gene flow and mixing. Some degree of genetic mixing can be detected, but putative hybrids were rare, only one of 175 specimens being a clear hybrid. In contrast, we found phenotypic differentiation between populations of H. timareta thelxinoe, possibly indicative of strong selection for local mimicry in different communities. In this pair of species, the absence of breakdown of genetic isolation despite near-identical wing patterns implies that factors other than wing patterns keep the two taxa apart, such as chemical or behavioural signals, or ecological adaptation along a strong altitudinal gradient. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109, 830–847.

As the generation and use of genomic datasets is becoming increasingly common in all areas of biology, the need for resources to collate, analyse and present data from independent (Tier 1) species-level genome projects into well supported clade-oriented (Tier 2) databases and provide a mechanism for these data to be propagated to pan-taxonomic (Tier 3) databases is becoming more pressing. Lepbase is a Tier 2 genomic resource for the Lepidoptera, supporting a research community using genomic approaches to understand evolution, speciation, olfaction, behaviour and pesticide resistance in a wide range of target species. Lepbase offers a core set of tools to make genomic data widely accessible including an Ensembl genome browser, text and sequence homology searches and bulk downloads of consistently presented and formatted datasets. As a part of the taxonomic community that we serve, we are working directly with Lepidoptera researchers to prioritise analyses and add tools that will be of most value to current research questions.

Ecomorphological differentiation is a key feature of adaptive radiations, with a general trend for specialization and niche expansion following divergence. Ecological opportunity afforded by invasion of a new habitat is thought to act as an ecological release, facilitating divergence, and speciation. Here, we investigate trophic adaptive morphology and ecology of an endemic clade of oreochromine cichlid fishes (Alcolapia) that radiated along a herbivorous trophic axis following colonization of an isolated lacustrine environment, and demonstrate phenotype‐environment correlation. Ecological and morphological divergence of the Alcolapia species flock are examined in a phylogenomic context, to infer ecological niche occupation within the radiation. Species divergence is observed in both ecology and morphology, supporting the importance of ecological speciation within the radiation. Comparison with an outgroup taxon reveals large‐scale ecomorphological divergence but shallow genomic differentiation within the Alcolapia adaptive radiation. Ancestral morphological reconstruction suggests lake colonization by a generalist oreochromine phenotype that diverged in Lake Natron to varied herbivorous morphologies akin to specialist herbivores in Lakes Tanganyika and Malawi.

An important goal in evolutionary biology is to understand the genetic changes underlying novel morphological structures. We investigated the origins of a complex wing pattern found among Amazonian Heliconius butterflies. Genome sequence data from 142 individuals across 17 species identified narrow regions associated with two distinct red colour pattern elements, dennis and ray. We hypothesise that these modules in non-coding sequence represent distinct cis-regulatory loci that control expression of the transcription factor optix, which in turn controls red pattern variation across Heliconius. Phylogenetic analysis of the two elements demonstrated that they have distinct evolutionary histories and that novel adaptive morphological variation was created by shuffling these cis-regulatory modules through recombination between divergent lineages. In addition, recombination of modules into different combinations within species further contributes to diversity. Analysis of the timing of diversification in these two regions supports the hypothesis of introgression moving regulatory modules between species, rather than shared ancestral variation. The dennis phenotype introgressed into Heliconius melpomene at about the same time that ray originated in this group, while ray introgressed back into H. elevatus much more recently. We show that shuffling of existing enhancer elements both within and between species provides a mechanism for rapid diversification and generation of novel morphological combinations during adaptive radiation., Novel butterfly wing patterns have arisen via the sharing of gene regulatory elements between lineages through introgression and recombination; new combinations of regulatory elements offer a rapid route to the evolution of new patterns., Author Summary Butterflies show an amazing diversity of patterns on their wings. In fact, most of the 18,000 species of butterfly can be distinguished on the basis of their wing pattern. Much of this diversity is thought to arise through novel switches in the genome that turn genes on in new contexts during wing development, thereby producing new patterns. Here we study a set of switches that control the expression of optix, a gene that places red patches onto the wings of Heliconius butterflies. We show that two patterning switches—one that produces red rays on the hindwing and the other a red patch on the base of the forewing—are located adjacent to one another in the genome. These switches have each evolved just once among a group of 16 species but have then been repeatedly shared between species by hybridisation and introgression. Despite the fact that they are now part of a common pattern in the Amazon basin, these two pattern components actually arose in completely different species before being brought together through hybridisation. In addition, recombination among these switches has produced new combinations of patterns within species. Such sharing of genetic variation is one way in which mimicry can evolve, whereby patterns are shared between species to send a common signal to predators. Our work suggests a new mechanism for generating evolutionary novelty, by shuffling these genetic switches among lineages and within species.

The Heliconius butterflies are a diverse recent radiation comprising multiple levels of divergence with ongoing gene flow between species. The recently sequenced genome of Heliconius melpomene allowed us to investigate the genomic evolution of this group using dense RAD marker sequencing. Phylogenetic analysis of 54 individuals robustly supported reciprocal monophyly of H. melpomene and Heliconius cydno and refuted previous phylogenetic hypotheses that H. melpomene may be paraphylectic with respect to H. cydno. Heliconius timareta also formed a monophyletic clade closely related but distinct from H. cydno with Heliconius heurippa falling within this clade. We find evidence for genetic admixture between sympatric populations of the sister clades H. melpomene and H. cydno/timareta, particularly between H. cydno and H. melpomene from Central America and between H. timareta and H. melpomene from the eastern slopes of the Andes. Between races, divergence is primarily explained by isolation by distance and there is no detectable genetic population structure between parapatric races, suggesting that hybrid zones between races are not zones of secondary contact. Our results also support previous findings that colour pattern loci are shared between populations and species with similar colour pattern elements. Furthermore, this pattern is almost unique to these genomic regions, with only a very small number of other loci showing significant similarity between populations and species with similar colour patterns.

The understanding of mimicry has relied on a strong biosystematic framework ever since early naturalists first recognized this textbook example of natural selection. We follow in this tradition, applying new biosystematics information to resolve problems in an especially difficult genus of tropical butterflies. Mechanitis species are important components of Neotropical mimetic communities. However, their colour pattern variability has presented challenges for systematists, and has made it difficult to study the very mimicry they so nicely illustrate. The South American Mechanitis mazaeus and relatives have remained particularly intractable. Recent systematists have recognized one highly polytypic species, whereas earlier work recognized the melanic Andean foothill races as a distinct species: Mechanitis messenoides. Recent molecular evidence suggests M. mazaeus and M. messenoides are genetically well differentiated, but evidence of morphological and ecological differences indicative of separate species was still lacking. Thus, it remains to be conclusively demonstrated whether this is an extreme case of a polymorphic mimetic species, or whether distinct co-mimetic lineages are involved. Here we provide evidence that M. mazaeus and M. messenoides are ecologically distinct and identify consistent morphological differences in both adult and immature stages. These ecological and morphological differences are correlated with mitochondrial sequence data. In spite of some overlap in almost all traits, wing shape, adult colour pattern, and larval colour pattern differ between the two species, in addition to clutch size and larval host use in local sympatry. Although three well-differentiated mitochondrial DNA (mtDNA) haplogroups were identified within these two species, one for M. mazaeus and two within M. messenoides, no morphological or ecological differences were found between two mtDNA haplogroups, both of which appear to belong to M. messenoides. We conclude that M. mazaeus and M. messenoides are distinct although highly polymorphic species, each with multiple sympatric co-mimetic forms, and suggest that further work is needed to clarify the identity of other phenotypes and subspecies of Mechanitis. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106, 540–560.

Attempts by biogeographers to understand biotic diversification in the Amazon have often employed contemporary species distribution patterns to support particular theories, such as Pleistocene rainforest refugia, rather than to test among alternative hypotheses. Suture zones, narrow regions where multiple contact zones and hybrid zones between taxa cluster, have been seen as evidence for past expansion of whole biotas that have undergone allopatric divergence in vicariant refuges. We used coalescent analysis of mutilocus sequence data to examine population split times in 22 pairs of geminate taxa in ithomiine and heliconiine butterflies. We test a hypothesis of simultaneous divergence across a suture zone in NE Peru. Our results reveal a scattered time course of diversification in this suture zone, rather than a tight cluster of split times. Additionally, we find rapid diversification within some lineages such as Melinaea contrasting with older divergence within lineages such as the Oleriina (Hyposcada and Oleria). These results strongly reject simple vicariance as a cause of the suture zone. At the same time, observed lineage effects are incompatible with a series of geographically coincident vicariant events which should affect all lineages similarly. Our results suggest that Pleistocene climatic forcing cannot readily explain this Peruvian suture zone. Lineage-specific biological traits, such as characteristic distances of gene flow or varying rates of parapatric divergence, may be of greater importance.

Supergenes are groups of tightly linked loci whose variation is inherited as a single Mendelian locus and are a common genetic architecture for complex traits under balancing selection1. Supergene alleles are long-range haplotypes with numerous mutations underlying distinct adaptive strategies, often maintained in linkage disequilibrium through the suppression of recombination by chromosomal rearrangements2–5. However, the mechanism governing the formation of supergenes is not well understood, and poses the paradox of establishing divergent functional haplotypes in face of recombination1,6. Here, we show that the formation of the supergene alleles encoding mimicry polymorphism in the butterfly Heliconius numata is associated with the introgression of a divergent, inverted chromosomal segment. Haplotype divergence and linkage disequilibrium indicate that supergene alleles, each allowing precise wing-pattern resemblance to distinct butterfly models, originate from over a million years of independent chromosomal evolution in separate lineages. These “superalleles” have evolved from a chromosomal inversion captured by introgression and maintained in balanced polymorphism, triggering supergene inheritance. This mode of evolution is likely to be a common feature of complex structural polymorphisms associated with the coexistence of distinct adaptive syndromes, and shows that the reticulation of genealogies may have a powerful influence on the evolution of genetic architectures in nature.

Mimicry and extensive geographical subspecies polymorphism combine to make species in the ithomiine butterfly genus Mechanitis (Lepidoptera; Nymphalidae) difficult to determine. We use mitochondrial DNA (mtDNA) barcoding, nuclear sequences and amplified fragment length polymorphism (AFLP) genotyping to investigate species limits in this genus. Although earlier biosystematic studies based on morphology described only four species, mtDNA barcoding revealed eight well-differentiated haplogroups, suggesting the presence of four new putative 'cryptic species'. However, AFLP markers supported only one of these four new 'cryptic species' as biologically meaningful. We demonstrate that in this genus, deep genetic divisions expected on the basis of mtDNA barcoding are not always reflected in the nuclear genome, and advocate the use of AFLP markers as a check when mtDNA barcoding gives unexpected results.

Amplified fragment length polymorphisms (AFLPs) are widely used for phylogenetic reconstruction in plants but their use in animal taxa has been limited, and in mammals rare. In addition, their use has been largely limited to shallow relationships amongst species or subspecies. Here, we genotype 23 pinniped species for 310 AFLP markers and find a strong phylogenetic signal, with individuals coclustering within species, and overall a good agreement between our phylogeny and those constructed using mitochondrial DNA and nuclear sequences even at nodes approximately 15 million years old. Although supporting the existing ideas about pinniped relationships, our data shed light on relationships within the hitherto relatively unresolved Phocine species group, and provide further supporting evidence for raising two subspecies of Zalophus californianus, Z. c. californianus and Z. c. wollebaeki, to species level. Plotting AFLP divergence time estimates against those based on both mtDNA and nuclear sequences we find strong linear relationships, suggesting that the different markers are evolving in a clocklike fashion. These data further emphasize the utility of AFLP markers as general tools for phylogenetic reconstruction.

In the absence of detailed pedigree records, researchers have attempted to estimate individuals' levels of inbreeding using molecular markers, generally making use of heterozygosity measures based on microsatellite markers. Here we report and validate a method for estimating an individual's inbreeding coefficient, f, using amplified fragment length polymorphism (AFLP) markers. We use simulations to confirm that our measure scales appropriately with f when allele frequencies can be estimated from a subset of outbred individuals. We also present an approach for obtaining satisfactory estimates even in the absence of an independent set of known outbred individuals from which to estimate allele frequencies. We then test our method against empirical data from 179 wild and captive-bred old-field mice, Peromyscus polionotus subgriseus, comprising pedigree-based estimates of f, along with genetic data from 94 AFLP markers and 12 microsatellites. Inbreeding estimates based on both AFLP and microsatellite markers were found to correlate strongly with pedigree-based inbreeding coefficients. Owing to their ease of amplification in any species, AFLP markers may prove to be a valuable new tool for estimating f in natural populations and for examining correlations between heterozygosity and fitness.

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridise in the wild. Here, we present a substantially improved assembly of the Heliconius melpomene genome, developed using novel methods that should be applicable to improving other genome assemblies produced using short read sequencing. Firstly, we whole genome sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Secondly, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Thirdly, we incorporated ~20x coverage of Pacific Biosciences sequencing and scaffolded the haploid genome using an assembly of this long read sequence. These improvements result in a genome of 795 scaffolds, 275 Mb in length, with an L50 of 2.1 Mb, an N50 of 34 and with 99% of the genome placed and 84% anchored on chromosomes. We use the new genome assembly to confirm that the Heliconius genome underwent 10 chromosome fusions since the split with its sister genus Eueides, over a period of about 6 million years.

Background Although hybridization is thought to be relatively rare in animals, the raw genetic material introduced via introgression may play an important role in fueling adaptation and adaptive radiation. The butterfly genus Heliconius is an excellent system to study hybridization and introgression but most studies have focused on closely related species such as H. cydno and H. melpomene. Here we characterize genome-wide patterns of introgression between H. besckei, the only species with a red and yellow banded ‘postman’ wing pattern in the tiger-striped silvaniform clade, and co-mimetic H. melpomene nanna. Results We find a pronounced signature of putative introgression from H. melpomene into H. besckei in the genomic region upstream of the gene optix, known to control red wing patterning, suggesting adaptive introgression of wing pattern mimicry between these two distantly related species. At least 39 additional genomic regions show signals of introgression as strong or stronger than this mimicry locus. Gene flow has been on-going, with evidence of gene exchange at multiple time points, and bidirectional, moving from the melpomene to the silvaniform clade and vice versa. The history of gene exchange has also been complex, with contributions from multiple silvaniform species in addition to H. besckei. We also detect a signature of ancient introgression of the entire Z chromosome between the silvaniform and melpomene/cydno clades. Conclusions Our study provides a genome-wide portrait of introgression between distantly related butterfly species. We further propose a comprehensive and efficient workflow for gene flow identification in genomic data sets. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-0889-0) contains supplementary material, which is available to authorized users.

Müllerian mimicry among Neotropical Heliconiini butterflies is an excellent example of natural selection, associated with the diversification of a large continental-scale radiation. Some of the processes driving the evolution of mimicry rings are likely to generate incongruent phylogenetic signals across the assemblage, and thus pose a challenge for systematics. We use a data set of 22 mitochondrial and nuclear markers from 92% of species in the tribe, obtained by Sanger sequencing and de novo assembly of short read data, to re-examine the phylogeny of Heliconiini with both supermatrix and multispecies coalescent approaches, characterize the patterns of conflicting signal, and compare the performance of various methodological approaches to reflect the heterogeneity across the data. Despite the large extent of reticulate signal and strong conflict between markers, nearly identical topologies are consistently recovered by most of the analyses, although the supermatrix approach failed to reflect the underlying variation in the history of individual loci. However, the supermatrix represents a useful approximation where multiple rare species represented by short sequences can be incorporated easily. The first comprehensive, time-calibrated phylogeny of this group is used to test the hypotheses of a diversification rate increase driven by the dramatic environmental changes in the Neotropics over the past 23 myr, or changes caused by diversity-dependent effects on the rate of diversification. We find that the rate of diversification has increased on the branch leading to the presently most species-rich genus Heliconius, but the change occurred gradually and cannot be unequivocally attributed to a specific environmental driver. Our study provides comprehensive comparison of philosophically distinct species tree reconstruction methods and provides insights into the diversification of an important insect radiation in the most biodiverse region of the planet.

Males with elaborate secondary sexual traits can enhance their mating success both by choosing to display at a site that increases the efficacy of their signal and by displaying at a site where the population is structured to provide inclusive fitness benefits via kin selection. Theoretical models and field studies reveal that, across a range of species and especially those with nonresource-based mating systems, males do display with kin. Bowerbirds use a nonresource-based mating system where males show spatial site fidelity and females visit the sites for the purpose of mate choice. However, neighbours can interfere with bower display by marauding (destroying and decoration stealing). Models suggest that local nonaggression pacts provide the best conditional strategy for males, a situation that may be facilitated by kin-mediated inclusive fitness benefits if males displayed together. We used a novel method of calculating relatedness, based on amplified fragment length polymorphism (AFLP), in a population of spotted bowerbirds, Chlamydera maculata, and found no evidence of fine-scale genetic structuring of display sites. Furthermore, males did not display next to kin, but instead an owner's neighbours were those most likely to maraud his bower. Marauders did not discriminate between kin and nonkin as their victims. Such theft and destruction may affect the mating success of the targeted bower owner.

Recent work has shown that elaborate secondary sexual traits and the corresponding preferences for them may be transmitted culturally rather than by genetic inheritance. Evidence for such cultural transmission commonly invokes spatial patterns of local similarity, with neighbouring individuals or populations appearing similar to each other. Alternative explanations for local similarity include ecological similarity of neighbouring environments and confounding genetic effects caused by aggregations of kin. We found that bowers built by male spotted bowerbirds, Chlamydera maculata, within a single population showed fine-scale similarities between neighbours in the decorations displayed on them. Such similarities did not covary with local decoration availability, local display environment or kinship and could not be explained by stealing behaviour by neighbours. Instead, we suggest that these similarities are products of local tradition, either culturally transmitted by neighbouring males who regularly inspect neighbours' bowers, or as localized responses to variable individual female preferences.

While Anartia butterflies have served as model organisms for research on the genetics of speciation, no phylogeny has been published to describe interspecific relationships. Here, we present a molecular phylogenetic analysis of Anartia species relationships, using both mitochondrial and nuclear genes. Analyses of both data sets confirm earlier predictions of sister species pairings based primarily on genital morphology. Yet both the mitochondrial and nuclear gene phylogenies demonstrate that Anartia jatrophae is not sister to all other Anartia species, but rather that it is sister to the Anartia fatima-Anartia amathea lineage. Traditional biogeographic explanations for speciation across the genus relied on A. jatrophae being sister to its congeners. These explanations invoked allopatric divergence of sister species pairs and multiple sympatric speciation events to explain why A. jatrophae flies alongside all its congeners. The molecular phylogenies are more consistent with lineage divergence due to vicariance, and range expansion of A. jatrophae to explain its sympatry with congeners. Further interpretations of the tree topologies also suggest how morphological evolution and eco-geographic adaptation may have set species range boundaries.

Anartia fatima and Anartia amathea (Lepidoptera: Nymphalidae) are sister taxa whose ranges abut in a narrow hybrid zone in eastern Panama. At the center of the zone, hybrids are abundant, although deviations from Hardy-Weinberg and linkage disequilibria are strong, due in part to assortative mating. We measured differences across the zone in four wing color-pattern characters, three allozyme loci, and mitochondrial haplotype. Wing pattern, allozyme, and mitochondrial clines were coincident (i.e., had the same positions) and concordant (i.e., all markers had similar cline shapes, about 28 km wide). Repeated samples demonstrated that the hybrid zone has been moving eastwards at an average rate of 2.5 km/year over the past 20 years, accompanied by an equivalent movement of the mtDNA cline. No introgression of mtDNA haplotypes were found in the ''wake'' of the moving cline, as might be expected for a neutral marker. The concordance of morphological and mtDNA clines between 1994 and 2000, in spite of hybrid zone movement, suggests strong epistasis between the mitochondrial genome and nuclear loci. Cline move- ment is achieved mainly by pure fatima immigrating into amathea populations; hybrids had little effect, and were presumably outcompeted by fitter purefatima genotypes. This movement can be explained if random dispersal of 7- 19 km.gen 21/2 is coupled with a competitive advantage to A. fatima genomes of 2-5%. Hybrid zone motion is equivalent to Phase III of Wright's shifting balance. Hybrid zone movement has rarely been considered likely in the past, but our results show that it may be more important in biogeography and evolution than generally realized.

Ten new dinucleotide microsatellite loci were isolated from the Antarctic fur seal Arctocephalus gazella. These markers should prove useful for studying the reproductive ecology of Antarctic fur seals and other related pinniped species.

Multilocus clines between Müllerian mimetic races of Heliconius butterflies provide a classic example of the maintenance of hybrid zones and their importance in speciation. Concordant hybrid zones in the mimics Heliconius erato and H. melpomene in northern Peru were carefully documented in the 1980s, and this prior work now permits a historical analysis of the movement or stasis of the zones. Previous work predicted that these zones might be moving toward the Andes due to selective asymmetry. Extensive deforestation and climate change might also be expected to affect the positions and widths of the hybrid zones. We show that the positions and shapes of these hybrid zones have instead remained remarkably stable between 1985 and 2012. The stability of this interaction strongly implicates continued selection, rather than neutral mixing following secondary contact. The stability of cline widths and strong linkage disequilibria (gametic correlation coefficients Rmax = 0.35-0.56 among unlinked loci) over 25 years suggest that mimetic selection pressures on each color pattern locus have remained approximately constant (s ≈ 0.13-0.40 per locus in both species). Exceptionally high levels of precipitation at the edge of the easternmost Andes may act as a population density trough for butterflies, trapping the hybrid zones at the foot of the mountains, and preventing movement. As such, our results falsify one prediction of the Pleistocene Refugium theory: That the ranges of divergent species or subspecies should be centered on regions characterized by maxima of rainfall, with hybrid zones falling in more arid regions between them.

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.

Tropical butterflies in the genus Heliconius have long been models in the study of the stages of speciation. Heliconius are unpalatable to predators, and many species are notable for multiple geographic populations with striking warning colour pattern differences associated with Mullerian mimicry. A speciation continuum is evident in Heliconius hybrid zones. Examples range from hybrid zones across which (a) there is little genetic differentiation other than at mimicry loci, but where hybrids are common, (b) to 'bimodal' hybrid zones with strong genetic divergence and few hybrids, (c) through to 'good' sympatric species, with hybridization extremely rare or absent. Now, in this issue of Molecular Ecology, Arias et al. (2012) have found an intermediate case in Colombian Heliconius cydno showing evidence for assortative mating and molecular differences, but where hybrids are abundant.