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.

Variation in cyanogenic compounds concentration within a Heliconius butterfly community: does mimicry explain everything?

BACKGROUND: Aposematic species advertise their unpalatability using warning signals such as striking coloration. Given that predators need to sample aposematic prey to learn that they are unprofitable, prey with similar warning signals share the cost of predator learning. This reduction in predation risk drives evolutionary convergence of warning signals among chemically defended prey (Müllerian mimicry). Whether such warning signal convergence is associated to similar defence levels among co-mimics is still an open question that has rarely been tested in wild populations. We quantified variation in cyanide-based (CN) chemical protection in wild caught individuals of eight aposematic Heliconius butterfly species belonging to four sympatric mimicry rings. We then tested for correlations between chemical protection and ecological species-specific traits. RESULTS: We report significant differences in CN concentrations both within and between sympatric species, even when accounting for the phylogeny, and within and between mimicry rings, even after considering inter-specific variation. We found significant correlations between CN concentration and both hostplant specialization and gregarious behaviour in adults and larvae. However, differences in CN concentrations were not significantly linked to mimicry ring abundance, although the two most toxic species did belong to the rarest mimicry ring. CONCLUSIONS: Our results suggest that mimicry can explain the variation in the levels of chemical defence to a certain extent, although other ecological factors are also relevant to the evolution of such variability.

Into the Andes: multiple independent colonizations drive montane diversity in the Neotropical clearwing butterflies Godyridina.

Understanding why species richness peaks along the Andes is a fundamental question in the study of Neotropical biodiversity. Several biogeographic and diversification scenarios have been proposed in the literature, but there is confusion about the processes underlying each scenario, and assessing their relative contribution is not straightforward. Here, we propose to refine these scenarios into a framework which evaluates four evolutionary mechanisms: higher speciation rate in the Andes, lower extinction rates in the Andes, older colonization times and higher colonization rates of the Andes from adjacent areas. We apply this framework to a species-rich subtribe of Neotropical butterflies whose diversity peaks in the Andes, the Godyridina (Nymphalidae: Ithomiini). We generated a time-calibrated phylogeny of the Godyridina and fitted time-dependent diversification models. Using trait-dependent diversification models and ancestral state reconstruction methods we then compared different biogeographic scenarios. We found strong evidence that the rates of colonization into the Andes were higher than the other way round. Those colonizations and the subsequent local diversification at equal rates in the Andes and in non-Andean regions mechanically increased the species richness of Andean regions compared to that of non-Andean regions ('species-attractor' hypothesis). We also found support for increasing speciation rates associated with Andean lineages. Our work highlights the importance of the Andean slopes in repeatedly attracting non-Andean lineages, most likely as a result of the diversity of habitats and/or host plants. Applying this analytical framework to other clades will bring important insights into the evolutionary mechanisms underlying the most species-rich biodiversity hotspot on the planet.

Molecular phylogenetics of the neotropical butterfly subtribe Oleriina (Nymphalidae: Danainae: Ithomiini).

The Oleriina is one of the most speciose subtribes of the neotropical nymphalid butterfly tribe Ithomiini. They are widely distributed across the Andes and Amazonian lowlands and like other ithomiines they are involved in complex mimicry rings. This subtribe is of particular interest because it contains the most diverse ithomiine genus, Oleria, as well as two genera, Megoleria and Hyposcada, that feed on hostplants not utilized elsewhere in the tribe. Here we present the first comprehensive species-level phylogeny for the Oleriina, representing 83% of recognised species in the group, and based on 6698bp from eight mitochondrial (mt) and nuclear (nc) genes. Topologies are largely congruent for ncDNA and the concatenated dataset and the genera Oleria, Hyposcada and Megoleria are recovered and well-supported, although strongly discordant genealogy between mtDNA and ncDNA suggest possible introgression among Hyposcada and Megoleria. A fourth clade containing the type species of Ollantaya is consistently recovered, and this recently synonymized name is resurrected. Clear subdivisions within Oleria separate the genus into four species groups, onega, amalda, makrena and aegle, which also correspond to differing biogeographic and elevation range characteristics. Unlike other ithomiine genera, the Oleriina show homogeneity in mimetic wing pattern, in sharp contrast to the emerging paradigm that mimetic shifts have enhanced diversification in the tribe. Our results show a potentially more important role for geographic isolation in the diversification of the Oleriina compared to other Ithomiini studied to date and provide a framework for more detailed biogeographical studies, in addition to a rare opportunity for comparative analyses with other neotropical groups.

Contrasting patterns of Andean diversification among three diverse clades of Neotropical clearwing butterflies.

The Neotropical region is the most biodiverse on Earth, in a large part due to the highly diverse tropical Andean biota. The Andes are a potentially important driver of diversification within the mountains and for neighboring regions. We compared the role of the Andes in diversification among three subtribes of Ithomiini butterflies endemic to the Neotropics, Dircennina, Oleriina, and Godyridina. The diversification patterns of Godyridina have been studied previously. Here, we generate the first time-calibrated phylogeny for the largest ithomiine subtribe, Dircennina, and we reanalyze a published phylogeny of Oleriina to test different biogeographic scenarios involving the Andes within an identical framework. We found common diversification patterns across the three subtribes, as well as major differences. In Dircennina and Oleriina, our results reveal a congruent pattern of diversification related to the Andes with an Andean origin, which contrasts with the Amazonian origin and multiple Andean colonizations of Godyridina. In each of the three subtribes, a clade diversified in the Northern Andes at a faster rate. Diversification within Amazonia occurred in Oleriina and Godyridina, while virtually no speciation occurred in Dircennina in this region. Dircennina was therefore characterized by higher diversification rates within the Andes compared to non-Andean regions, while in Oleriina and Godyridina, we found no difference between these regions. Our results and discussion highlight the importance of comparative approaches in biogeographic studies.

Diversification of clearwing butterflies with the rise of the Andes.

AIM: Despite the greatest butterfly diversity on Earth occurring in the Neotropical Andes and Amazonia, there is still keen debate about the origins of this exceptional biota. A densely sampled calibrated phylogeny for a widespread butterfly subtribe, Oleriina (Nymphalidae: Ithomiini) was used to estimate the origin, colonization history and diversification of this species-rich group. LOCATION: Neotropics. METHODS: Ancestral elevation and biogeographical ranges were reconstructed using data generated from detailed range maps and applying the dispersal-extinction-cladogenesis model using stratified palaeogeographical time slice matrices. The pattern of diversification through time was examined by comparing constant and variable rate models. We also tested the hypothesis that a change in elevation is associated with speciation. RESULTS: The Oleriina likely originated in the Andes in the Early to Middle Miocene and rapidly diversified to include four genera all of which also originated in the Andes. These clades, together with four species groups, experienced varying spatial and temporal patterns of diversification. An overall early burst and decreasing diversification rate is identified, and this pattern is reflected for most subclades. MAIN CONCLUSIONS: Changes in the palaeogeological landscape, particularly the prolonged uplift of the Andes, had a profound impact on the diversification of the subtribe. The Oleriina mostly remained within the Andes and vicariant speciation resulted in some instances. Dynamic dispersal occurred with the disappearance of geological barriers such as the Acre System and the subtribe exploited newly available habitats. Our results confirm the role of the Andean uplift in the evolution of Neotropical biodiversity.