Search

How frequent is gene flow between species? The pattern of evolution is typically portrayed as a phylogenetic tree, implying that speciation is a series of splits between lineages. Yet gene flow between good species is increasingly recognized as an important mechanism in the diversification of radiations, often spreading adaptive traits and leading to a complex pattern of phylogenetic incongruence. This process has thus far been studied in cases involving few species, or geographically restricted to spaces like islands, but not on the scale of a continental radiation. Previous studies have documented gene flow, adaptive introgression and hybrid speciation in a small subsection of the charismatic Neotropical butterflies Heliconius. Using genome-wide resequencing of 40 out of 45 species in the genus we demonstrate for the first time that admixture has played a role throughout the evolution of Heliconius and the sister genus Eueides. Modelling of phylogenetic networks based on 6848 orthologous autosomal genes (Maximum Pseudo-Likelihood Networks) or 5,483,419 high quality SNPs (Ancestral Recombination Graph) uncovers nine new cases of interspecific gene flow at up to half of the genome. However, f4 statistics of admixture show that the extent of the process has varied between subgenera. Evidence for introgression is found at all five loci controlling the colour and shape of the mimetic wing patterns, including in the putative hybrid species H. hecalesia, characterised by an unusual hindwing. Due to hybridization and incomplete coalescence during rapid speciation, individual gene trees show rampant discordance. Although reduced gene flow and faster coalescence are expected at the Z chromosome, we discover high levels of conflict between the 416 sex-linked loci. Despite this discordant pattern, both concatenation and multispecies coalescent approaches yield surprisingly consistent and fully supported genome-wide phylogenies. We conclude that the imposition of the bifurcating tree model without testing for interspecific gene flow may distort our perception of adaptive radiations and thus the ability to study trait evolution in a comparative framework.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

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.

The taxonomic identity of the butterfly Heliconius melpomene f. pyritosa var. fumigata Zikán, 1937 (Nymphalidae: Heliconiinae) is discussed based on the discovery of the specimen based on which this name was assigned. The specimen is not a variation of Eueides tales surdus Stichel, 1903, as previously stated, but is in fact a variation of H. melpomene f. pyritosa, which is a synonym of Heliconius erato amalfreda Riffarth, 1901. Keywords: Eueides, Heliconiini, Heliconius erato

A new phylogenetic hypothesis for Heliconius and related genera is presented, based on DNA sequence data from mtDNA combined with a region of the wingless gene. This study also adds eight new taxa to a previous cladistic hypothesis based on the mtDNA alone. Simultaneous phylogenetic analysis of the two gene regions together supports a topology largely in agreement with traditional views of heliconiine relationships based on morphology and suggests that the mtDNA support for the sister relationship between Eueides and H. charithonia is due to convergent evolution of homoplasious mtDNA sites.

The genus Heliconius has been revised more than a dozen times, yet relationships among many of its species groups remain obscure. A reliable phylogenetic hypothesis is desirable, because the genus has been a model system for studies of tropical community ecology, mimicry, and ecological genetics for the last 3 decades. A new cladogram for 35 species of Heliconius and the related genera Eueides, Laparus, and Neruda is presented, based on mitochondrial sequences spanning part of the COI gene and the COII gene. The data support most traditionally recognized species groups and also the monophyly of the above four genera with respect to other heliconiine outgroups. However, Heliconius is paraphyletic with respect to the other three genera. These data will allow a reinvestigation of problematical morphological, behavioral, and cytological traits in the group.

Chromosome counts in meiotic metaphase plates in the gonads of 67 of the probable 68 species of mimetic neotropical heliconian butterflies (Nymphalidae), representing 1524 individuals in 617 subspecies and geographically separate populations from southern Texas to northern Argentina, revealed a consistent haploid number of n = 21 in the genus Heliconius (except for the most advanced species with n = 33, 37, 56, and 60) and n = 31 in the more primitive genera (Eueides, Dryas, Dryadula, Agraulis, and Dione), with a transitional genus (Neruda) showing three species with n = 28-32, 21-22 + 5-10 "microchromosomes", and 20-22 + 1-5 "microchromosomes". The genus Laparus, with a single polymorphic species doris, probably an offshoot of early Heliconius, shows wide karyotypic variation (n = 20-30, 38) sometimes even within a single individual. The two most primitive genera also show much variation: Podotricha has two species with n = 9 and n = 26-29; and Philaethria shows many phenotypically similar species, two with n = 29 and a still uncertain number (at least 3) with n = 88 (most common), 67-72 (most widespread), 62 (very restricted geographically), 52, 21, and 12. Several interspecific hybrids (Heliconius cydno x H. melpomene) showed normal chromosome pairing, while deficient pairing was seen in intersubspecific hybrids in Eueides tales and Heliconius sara. The importance of these results in the evolutionary study of polytypic tropical species is discussed.

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridize 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. First, we whole-genome-sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Second, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Third, 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 N50 length of 2.1 Mb, an N50 number 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 yr.

Background While the ecological factors that drive phenotypic radiations are often well understood, less is known about the generative mechanisms that cause the emergence and subsequent diversification of novel features. Heliconius butterflies display an extraordinary diversity of wing patterns due in part to mimicry and sexual selection. Identifying the genetic drivers of this crucible of evolution is now within reach, as it was recently shown that cis-regulatory variation of the optix transcription factor explains red pattern differences in the adaptive radiations of the Heliconius melpomene and Heliconius erato species groups. Results Here, we compare the developmental expression of the Optix protein across a large phylogenetic sample of butterflies and infer that its color patterning role originated at the base of the neotropical passion-vine butterfly clade (Lepidoptera, Nymphalidae, Tribe: Heliconiini), shortly predating multiple Optix-driven wing pattern radiations in the speciose Heliconius and Eueides genera. We also characterize novel Optix and Doublesex expression in the male-specific pheromone wing scales of the basal heliconiines Dryas and Agraulis, thus illustrating that within the Heliconinii lineage, Optix has been evolutionarily redeployed in multiple contexts in association with diverse wing features. Conclusions Our findings reveal that the repeated co-option of Optix into various aspects of wing scale specification was associated with multiple evolutionary novelties over a relatively short evolutionary time scale. In particular, the recruitment of Optix expression in colored scale cell precursors was a necessary condition to the explosive diversification of passion-vine butterfly wing patterns. The novel deployment of a gene followed by spatial modulation of its expression in a given cell type could be a common mode of developmental innovation for triggering phenotypic radiations.