Carlos F. Arias

Out of the Andes: patterns of diversification in clearwing butterflies

Global biodiversity peaks in the tropical forests of the Andes, a striking geological feature that has likely been instrumental in generating biodiversity by providing opportunities for both vicariant and ecological speciation. However, the role of these mountains in the diversification of insects, which dominate biodiversity, has been poorly explored using phylogenetic methods. Here we study the role of the Andes in the evolution of a diverse Neotropical insect group, the clearwing butterflies. We used dated species‐level phylogenies to investigate the time course of speciation and to infer ancestral elevation ranges for two diverse genera. We show that both genera likely originated at middle elevations in the Andes in the Middle Miocene, contrasting with most published results in vertebrates that point to a lowland origin. Although we detected a signature of vicariance caused by the uplift of the Andes at the Miocene–Pliocene boundary, most sister species were parapatric without any obvious vicariant barrier. Combined with an overall decelerating speciation rate, these results suggest an important role for ecological speciation and adaptive radiation, rather than simple vicariance.

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.

Hybrid incompatibility is consistent with a hybrid origin of Heliconius heurippa Hewitson from its close relatives, Heliconius cydno Doubleday and Heliconius melpomene Linnaeus.

Shared ancestral variation and introgression complicates the reconstruction of phylogenetic relationships among closely related taxa. Here we use overall genomic compatibility as an alternative estimate of species relationships in a group where divergence is rapid and genetic exchange is common. Heliconius heurippa, a butterfly species endemic to Colombia, has a colour pattern genetically intermediate between H. cydno and H. melpomene: its hindwing is nearly indistinguishable from that of H. melpomene and its forewing band is an intermediate phenotype between both species. This observation has lead to the suggestion that the pattern of H. heurippa arose through hybridization. We present a genetic analysis of hybrid compatibility in crosses between the three taxa. Heliconius heurippa × H. cydno and female H. melpomene × male H. heurippa yield fertile and viable F1 hybrids, but male H. melpomene × female H. heurippa crosses yield sterile F1 females. In contrast, Haldanes rule has previously been detected between H. melpomene and H cydno in both directions. Therefore, H. heurippa is most closely related to H. cydno, with some evidence for introgression of genes from H. melpomene. The results are compatible with the hypothesis of a hybrid origin for H. heurippa. In addition, backcrosses using F1 hybrid males provide evidence for a large Z(X)‐chromosome effect on sterility and for recessive autosomal sterility factors as predicted by Dominance Theory.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome‐scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene–Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome‐scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

Sharp genetic discontinuity across a unimodal Heliconius hybrid zone

Hybrid zones are powerful natural systems to study evolutionary processes to gain an understanding of adaptation and speciation. In the Cauca Valley (Colombia), two butterfly races, Heliconius cydno cydnides and Heliconius cydno weymeri, meet and hybridize. We characterized this hybrid zone using a combination of mitochondrial DNA (mtDNA) sequences, amplified fragment length polymorphisms (AFLPs), microsatellites and sequences for nuclear loci within and outside of the genomic regions that cause differences in wing colour pattern. The hybrid zone is largely composed of individuals of mixed ancestry. However, there is strong genetic discontinuity between the hybridizing races in mtDNA and, to a lesser extent, in all nuclear markers surveyed. The mtDNA clustering of H. c. cydnides with the H. cydno race from the Magdalena Valley and H. c. weymeri with the H. cydno race from the pacific coast suggests that H. c. cydnides colonized the Cauca Valley from the north, whereas H. c. weymeri did so by crossing the Andes in the southern part, implying a secondary contact origin. Colonization of the valley by H. cydno was accompanied by mimicry shift. Strong ecological isolation, driven by locally adaptive differences in mimetic wing patterns, is playing an important role in maintaining the hybrid zone. However, selection on wing pattern alone is not sufficient to explain the genetic discontinuity observed. There is evidence for differences in male mating preference, but the contribution of additional barriers needs further investigation. Overall, our results support the idea that speciation is a cumulative process, where the combination of multiple isolation barriers, combined with major phenotypic differences, facilitates population divergence in face of gene flow.

Diversity of Aplochiton Fishes (Galaxiidea) and the Taxonomic Resurrection of A. marinus

Aplochiton is a small genus of galaxiid fishes endemic to Patagonia and the Falkland Islands whose taxonomy is insufficiently resolved. Recent genetic analyses confirmed the existence of only two closely related species, Aplochiton taeniatus and Aplochiton zebra, while a third controversial species, Aplochiton marinus, remained lost to synonymy with A. taeniatus. Using an integrative taxonomy framework, we studied original samples and published sequences from a broad range in western Patagonia and the Falkland Islands, and generated robust species hypotheses based on single-locus (Cytochrome Oxidase subunit I; COI) species-delineation methods and known diagnostic morphological characters analyzed in a multivariate context. Results revealed three distinct evolutionary lineages that morphologically resemble, in important respects, existing nominal species descriptions. Interestingly, the lineage associated with A. marinus was unambiguously identifiable (100% accuracy) both from the genetic and morphological viewpoints. In contrast, the morphology of A. taeniatus and A. zebra overlapped substantially, mainly due to the high variability of A. taeniatus. Discriminant function analysis aided the identification of these species with 83.9% accuracy. Hence, for their unambiguous identification, genetic screening is needed. A. marinus has seldom been documented, and when recorded, it has always been found in sites with clear marine influence. It is possible that only A. marinus preserves a life cycle related to the sea akin to the hypothesized ancestral galaxiid. We did not find evidence of claimed diadromy in A. taeniatus or A. zebra, and, therefore, these should be regarded as freshwater species. Finally, a lack of phylogeographic patterns and overrepresentation of uncommon haplotypes suggested demographic expansions in recent evolutionary time, especially of A. zebra, in line with the hypothesis of large-scale range expansion and lineage spread in western Patagonia.

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole‐genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.

Genomics at the evolving species boundary

Molecular studies on hybridization date back to Dobzhansky who compared chromosomal banding patterns to determine if interspecific gene flow occurred in nature [1]. Now, the advent of high-throughput sequencing provides increasingly fine insights into genomic differentiation between incipient taxa that are changing our view of adaptation and speciation and the links between the two. Empirical data from hybridizing taxa demonstrate highly heterogeneous patterns of genomic differentiation. Although underlining reasons for this heterogeneity are complex, studies of hybridizing taxa offers some of the best insights into the regions of the genome under divergent selection and the role these regions play in species boundaries. The challenge moving forward is to develop a better theoretical framework that fully leverages these powerful natural experiments.

Renewed diversification following Miocene landscape turnover in a Neotropical butterfly radiation

The Neotropical region has experienced a dynamic landscape evolution throughout the Miocene, with the large wetland Pebas occupying western Amazonia until 11-8 my ago and continuous uplift of the Andes mountains along the western edge of South America. Although the complex dynamics between the Andes and Amazonia may have strongly affected the trajectory of Neotropical biodiversity, there is little evidence for such an influence from time-calibrated phylogenies of groups that diversified during this period. Here, we generate one of the most comprehensive time-calibrated molecular phylogenies of a group of Neotropical insects: the butterfly tribe Ithomiini. Our tree includes 340 species (87% of extant species), spanning 26 million years of diversification in the Neotropics. We investigate temporal and spatial patterns of diversification, focusing on the influence of Miocene landscape tranformations on the dynamics of speciation, extinction and biotic interchanges at the Amazonia/Andes interface. We find that Ithomiini likely began diversifying at the interface between the Andes and the Amazonia around 26.4 my ago. Five subtribes with a very low extant diversity started diversifying early in western Amazonia, but a rapid decrease in diversification rate due to increased extinction rate between 20 and 10 my ago suggests a negative impact of the Pebas wetland system on these early lineages. By contrast, the clade containing the five most species-rich subtribes (85% of extant species) was characterized by a high, positive net diversification rate. This clade diversified exclusively in the Central Andes from 20 to 10 my ago. After the demise of the Pebas system (11-8 my ago), we found a sudden increase of interchanges with the Northern Andes and Amazonia, followed by local diversification, which led to a substantial renewal of diversification. In general, ecological turnovers throughout the Miocene strongly determined the dynamics of speciation, and extinction and interchanges, and appear as a key driving force shaping the current extraordinary biodiversity of the Neotropical region.

Inheritance, distribution and genetic differentiation of a color polymorphism in Panamanian populations of the tortoise beetle, Chelymorpha alternans (Coleoptera: Chrysomelidae)

Intraspecific variation maintained in natural populations has long intrigued scientists and naturalists. One form of this variation, color polymorphisms, provide a rich opportunity to connect genotypic and phenotypic diversity within an ecological and evolutionary context. The existence of color polymorphisms in Panamanian populations of the Neotropical tortoise beetle, Chelymorpha alternans, has been suspected but never systematically explored. To characterize geographic distribution and underlying genetics we sampled a total of 3819 beetles from 28 sites across Panama, quantifying five distinct phenotypes. Two phenotypes, the “metallic” and “rufipennis” are the most widely distributed phenotypes, occurring in nearly all collecting sites. The “veraguensis” phenotype was found to be restricted to the Western end of the Isthmus and the “militaris” phenotypes restricted to sites east of the canal. Controlled matings between phenotypes and reared offspring revealed no indications of reproductive barriers, even among phenotypes which do not co-occur in nature. Color pattern phenotype is largely controlled by Mendelian assortment of four alleles competing at a single locus. A clear dominance hierarchy exists among alleles, with two being co-dominant. Genomic scans from 32 individuals revealed low levels of genetic differentiation, with a small fraction of the genome showing a high degree of divergence. The easily observed variation among populations, simple genetic architecture, and rearing capabilities, make this a promising system for investigating proximate and ultimate factors of phenotypic variation.