Movement of a Heliconius hybrid zone over 30 years: A Bayesian approach

Trait Variation and Spatiotemporal Dynamics across Avian Secondary Contact Zones

A secondary contact zone (SCZ) is an area where incipient species or divergent populations may meet, mate, and hybridize. Due to the diverse patterns of interspecific hybridization, SCZs function as field labs for illuminating the on-going evolutionary processes of speciation and the establishment of reproductive isolation. Interspecific hybridization is widely present in avian populations, making them an ideal system for SCZ studies. This review exhaustively summarizes the variations in unique traits within avian SCZs (vocalization, plumage, beak, and migratory traits) and the various movement patterns of SCZs observed in previous publications. It also highlights several potential future research directions in the genomic era, such as the relationship between phenotypic and genomic differentiation in SCZs, the genomic basis of trait differentiation, SCZs shared by multiple species, and accurate predictive models for forecasting future movements under climate change and human disturbances. This review aims to provide a more comprehensive understanding of speciation processes and offers a theoretical foundation for species conservation.

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Limited movement of an avian hybrid zone in relation to regional variation in magnitude of climate change

Studies of natural hybrid zones can provide documentation of range shifts in response to climate change and identify loci important to reproductive isolation. Using a temporal (36-38 years) comparison of the black-capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee hybrid zone, we investigated movement of the western portion of the zone (western Missouri) and assessed whether loci and pathways underpinning reproductive isolation were similar to those in the eastern portion of the hybrid zone. Using 92 birds sampled along the hybrid zone transect in 2016 and 68 birds sampled between 1978 and 1980, we generated 11,669 SNPs via ddRADseq. These SNPs were used to assess movement of the hybrid zone through time and to evaluate variation in introgression among loci. We demonstrate that the interface has moved ~5 km to the northwest over the last 36-38 years, that is, at only one-fifth the rate at which the eastern portion (e.g., Pennsylvania, Ohio) of the hybrid zone has moved. Temperature trends over the last 38 years reveal that eastern areas have warmed 50% more than western areas in terms of annual mean temperature, possibly providing an explanation for the slower movement of the hybrid zone in Missouri. Our results suggest hybrid zone movement in broadly distributed species, such as chickadees, will vary between areas in response to local differences in the impacts of climate change.

Scale-dependent environmental effects on phenotypic distributions in Heliconius butterflies

Identifying the relative importance of different mechanisms responsible for the emergence and maintenance of phenotypic diversity can be challenging, as multiple selective pressures and stochastic events are involved in these processes. Therefore, testing how environmental conditions shape the distribution of phenotypes can offer important insights on local adaptation, divergence, and speciation. The red-yellow Müllerian mimicry ring of Heliconius butterflies exhibits a wide diversity of color patterns across the Neotropics and is involved in multiple hybrid zones, making it a powerful system to investigate environmental drivers of phenotypic distributions. Using the distantly related Heliconius erato and Heliconius melpomene co-mimics and a multiscale distribution approach, we investigated whether distinct phenotypes of these species are associated with different environmental conditions. We show that Heliconius red-yellow phenotypic distribution is strongly driven by environmental gradients (especially thermal and precipitation variables), but that phenotype and environment associations vary with spatial scale. While co-mimics are usually predicted to occur in similar environments at large spatial scales, patterns at local scales are not always consistent (i.e., different variables are best predictors of phenotypic occurrence in different locations) or congruent (i.e., co-mimics show distinct associations with environment). We suggest that large-scale analyses are important for identifying how environmental factors shape broad mimetic phenotypic distributions, but that local studies are essential to understand the context-dependent biotic, abiotic, and historical mechanisms driving finer-scale phenotypic transitions.

Mito-nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern "Rocky Mountain'' color form with ferruginous mid-abdominal segments (B. m. melanopygus) and a southern "Pacific'' form with black mid-abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color-controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito-nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.

Genetic, morphological and ecological variation across a sharp hybrid zone between two alpine butterfly species

Identifying the mechanisms involved in the formation and maintenance of species is a central question in evolutionary biology, and distinguishing the selective drivers of populations' divergence from demographic processes is of particular interest to better understand the speciation process. Hybrid zones are recognized to provide ideal places to investigate the genetic architecture of speciation and to identify the mechanisms allowing diverging species to maintain their integrity in the face of gene flow. Here, we studied two alpine butterfly species, Coenonympha macromma and C. gardetta, which can be found flying together and hybridizing in narrow contact zones in the southern French Alps. We characterized the genomic composition of individuals, their morphology and their local habitat requirements, within and around a hybrid zone. Genetic diversity analysis at 794 SNPs revealed that all individuals within the hybrid zone were highly admixed, which was not the case outside the hybrid zone. Cline analysis showed that, despite ongoing hybridization, 56 out of 122 loci differentially fixed or nearly so between the two species were impermeable to introgression across the sharp hybrid zone (9 km wide). We also found concordance in cline position and width among genetic, morphological and environmental variation, suggesting a coupling of different reproductive barriers. Habitat characteristics such as the presence of trees and shrubs and the start of the growing season were strongly associated with the genetic variation, and we found evidence of divergence at genetic markers associated with morphology and physiology, putatively involved in visual or environmental reproductive isolation. We discuss the various behavioural and ecological factors that might interplay to maintain current levels of divergence and gene flow between this species pair.

Genome-wide patterns of differentiation over space and time in the Queensland fruit fly

The Queensland fruit fly, Bactrocera tryoni, is a major pest of Australian horticulture which has expanded its range in association with the spread of horticulture over the last ~ 150 years. Its distribution in northern Australia overlaps that of another fruit fly pest to which some authors accord full species status, Bactrocera aquilonis. We have used reduced representation genome-wide sequencing to genotype 359 individuals taken from 35 populations from across the current range of the two taxa, plus a further 73 individuals from six of those populations collected 15–22 years earlier. We find significant population differentiation along an east–west transect across northern Australia which likely reflects limited but bidirectional gene flow between the two taxa. The southward expansion of B. tryoni has led to relatively little genetic differentiation, and most of it is associated with a move into previously marginal inland habitats. Two disjunct populations elsewhere in Australia and three on Melanesian islands are each clearly differentiated from all others, with data strongly supporting establishment from relatively few founders and significant isolation subsequently. Resequencing of historical samples from one of the disjunct Australian populations shows that its genetic profile has changed little over a 15-year period, while the Melanesian data suggest a succession of ‘island hopping’ events with progressive reductions in genetic diversity. We discuss our results in relation to the control of B. tryoni and as a model for understanding the genetics of invasion and hybridisation processes.

Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation

Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation.