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.     

Cross‐decades stability of an avian hybrid zone

Hybrid zones are particularly valuable for understanding the evolution of partial reproductive isolation between differentiated populations. An increasing number of hybrid zones have been inferred to move over time, but in most such cases zone movement has not been tested with long‐term genomic data. The hybrid zone between Townsend's Warblers (Setophaga townsendi) and Hermit Warblers (S. occidentalis) in the Washington Cascades was previously inferred to be moving from northern S. townsendi southwards towards S. occidentalis, based on plumage and behavioural patterns as well as a 2000‐km genetic wake of hermit mitochondrial DNA (mtDNA) in coastal Townsend's Warblers. We directly tested whether hybrid zone position has changed over 2–3 decades by tracking plumage, mtDNA and nuclear genomic variation across the hybrid zone over two sampling periods (1987–94 and 2015–16). Surprisingly, there was no significant movement in genomic or plumage cline centres between the two time periods. Plumage cline widths were narrower than expected by neutral diffusion, consistent with a ‘tension zone’ model, in which selection against hybrids is balanced by movement of parental forms into the zone. Our results indicate that this hybrid zone is either stable in its location or moving at a rate that is not detectable over 2–3 decades. Despite considerable gene flow, the stable clines in multiple phenotypic and genotypic characters over decades suggest evolutionary stability of this young pair of sister species, allowing divergence to continue. We propose a novel biogeographic scenario to explain these patterns: rather than the hybrid zone having moved thousands of kilometres to its current position, inland Townsend's met coastal Hermit Warbler populations along a broad front of the British Columbia and Alaska coast and hybridization led to replacement of the Hermit Warbler plumage with Townsend's Warbler plumage patterns along this coastline. Hence, hybrid zones along British Columbia and Alaska moved only a short distance from the inland to the coast, whereas the Hermit Warbler phenotype appears stable in Washington and further south. This case provides an example of the complex biogeographic processes that have led to the distribution of current phenotypes within and among closely related species.     

Quantifying the spatiotemporal dynamics in a chorus frog (Pseudacris) hybrid zone over 30 years

Although theory suggests that hybrid zones can move or change structure over time, studies supported by direct empirical evidence for these changes are relatively limited. We present a spatiotemporal genetic study of a hybrid zone between Pseudacris nigrita and P. fouquettei across the Pearl River between Louisiana and Mississippi. This hybrid zone was initially characterized in 1980 as a narrow and steep “tension zone,” in which hybrid populations were inferior to parentals and were maintained through a balance between selection and dispersal. We reanalyzed historical tissue samples and compared them to samples of recently collected individuals using microsatellites. Clinal analyses indicate that the cline has not shifted in roughly 30 years but has widened significantly. Anthropogenic and natural changes may have affected selective pressure or dispersal, and our results suggest that the zone may no longer best be described as a tension zone. To the best of our knowledge, this study provides the first evidence of significant widening of a hybrid cline but stasis of its center. Continued empirical study of dynamic hybrid zones will provide insight into the forces shaping their structure and the evolutionary potential they possess for the elimination or generation of species.     

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological,genetic and climatic data to examine stability

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south‐eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90‐km transect through the hybrid zone and examined variation in 26 microsatellites, the D‐loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection‐gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.     

Molecular patterns of introgression in a classic hybrid zone between the Australian tree frogs,Litoria ewingii and L. paraewingi: evidence of a tension zone

Hybrid zones provide a rare opportunity to explore the processes involved in reproductive isolation and speciation. The southern hybrid zone between the southeastern Australian tree frogs Litoria ewingii and L. paraewingi has been comprehensively studied over the last 40 years, primarily using reproductive compatibility experiments and male advertisement calls. We used mitochondrial DNA (mtDNA) and eight nuclear microsatellite markers to characterize this hybrid zone along a historically studied transect and to test various dispersal‐dependent and dispersal‐independent hybrid zone models. The species are genetically distinct and the level of hybridization within the contact zone is low, with the majority of admixed individuals representing later‐generation hybrids. Based on previous experimental genetic compatibility studies, we predicted that hybrids with L. paraewingi mtDNA would be more frequent than hybrids with L. ewingii mtDNA. Surprisingly, a greater proportion of the identified hybrids had L. ewingii mtDNA. Geographical cline analyses showed a sharp transition in allele frequencies across the transect, and both the mtDNA and microsatellite data showed concordant cline centres, but were best supported by a model that allowed width to vary. Overall, the L. ewingii–L. paraewingi hybrid zone is best characterized as a tension zone, due to the narrow cline width, concordant genetic clines and low levels of hybridization.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

Ecological and genetic associations across a Heliconius hybrid zone

Differences in habitat use can bridge early and late stages of speciation by initiating assortative mating. Heliconius colour pattern races might select habitats over which each pattern confers a relative fitness advantage because signal efficacy of wing patterns can vary by environment. Thus habitat preferences could serve to promote the evolution of mimetic colour patterns for mate choice. Here I compare colour pattern genotype and phenotype frequencies to environmental variation across the H. erato hydara x H. erato erato hybrid zone in French Guiana to determine whether races exhibit habitat preferences. I found that genotype and phenotype frequencies correspond to differences in land cover moreso than to other environmental factors. Temporal shifts in colour pattern genotypes, phenotypes and land cover also were associated at individual sample sites, which further suggests that H. erato races differ in habitat use and that habitat preferences may promote speciation among Heliconius butterflies.     

Ants exhibit asymmetric hybridization in a mosaic hybrid zone

Research on hybridization between species provides unparalleled insights into the pre‐ and postzygotic isolating mechanisms that drive speciation. In social organisms, colony‐level incompatibilities may provide additional reproductive barriers not present in solitary species, and hybrid zones offer an opportunity to identify these barriers. Here, we use genotyping‐by‐sequencing to sequence hundreds of markers in a hybrid zone between two socially polymorphic ant species, Formica selysi and Formica cinerea. We characterize the zone, determine the frequency of hybrid workers, infer whether hybrid queens or males are produced and investigate whether hybridization is influenced by colony social organization. We also compare cuticular hydrocarbon profiles and aggression levels between the two species. The hybrid zone exhibits a mosaic structure. The asymmetric distribution of hybrids skewed towards F. cinerea suggests a pattern of unidirectional nuclear gene flow from F. selysi into F. cinerea. The occurrence of backcrossed individuals indicates that hybrid queens and/or males are fertile, and the presence of the F. cinerea mitochondrial haplotype in 97% of hybrids shows that successful F1 hybrids will generally have F. cinerea mothers and F. selysi fathers. We found no evidence that social organization contributes to speciation, because hybrids occur in both single‐queen and multiple‐queen colonies. Strongly differentiated cuticular hydrocarbon profiles and heightened interspecific aggression further reveal that species recognition cues are both present and perceived. The discovery of fertile hybrids and asymmetrical gene flow is unusual in ants, and this hybrid zone will therefore provide an ideal system with which to investigate speciation in social insects.     

Narrow hybrid zones in spite of very low population differentiation in neutral markers in an island bird species complex

Patterns of phenotypic and genic frequencies across hybrid zones provide insight into the origin and evolution of reproductive isolation. The Reunion grey white‐eye, Zosterops borbonicus, exhibits parapatrically distributed plumage colour forms across the lowlands of the small volcanic island of Reunion (Mascarene archipelago). These forms meet and hybridize in regions that are natural barriers to dispersal (rivers, lava fields). Here, we investigated the relationship among patterns of differentiation at neutral genetic (microsatellite) markers, phenotypic traits (morphology and plumage colour) and niche characteristics across three independent hybrid zones. Patterns of phenotypic divergence revealed that these hybrid zones are among the narrowest ever documented in birds. However, the levels of phenotypic divergence stand in stark contrast to the lack of clear population neutral genetic structure between forms. The position of the hybrid zones coincides with different natural physical barriers, yet is not associated with steep changes in vegetation and related climatic variables, and major habitat transitions are shifted from these locations by at least 18 km. This suggests that the hybrid zones are stabilized over natural dispersal barriers, independently of environmental boundaries, and are not associated with niche divergence. A striking feature of these hybrid zones is the very low levels of genetic differentiation in neutral markers between forms, suggesting that phenotypic divergence has a narrow genetic basis and may reflect recent divergence at a few linked genes under strong selection, with a possible role for assortative mating in keeping these forms apart.     

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome—that is, the genetic architecture of reproductive isolation—can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male‐transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.     

Genotype‐environment associations support a mosaic hybrid zone between two tidal marsh birds

Local environmental features can shape hybrid zone dynamics when hybrids are bounded by ecotones or when patchily distributed habitat types lead to a corresponding mosaic of genotypes. We investigated the role of marsh‐level characteristics in shaping a hybrid zone between two recently diverged avian taxa – Saltmarsh (Ammodramus caudacutus) and Nelson's (A. nelsoni) sparrows. These species occupy different niches where allopatric, with caudacutus restricted to coastal marshes and nelsoni found in a broader array of wetland and grassland habitats and co‐occur in tidal marshes in sympatry. We determined the influence of habitat types on the distribution of pure and hybrid sparrows and assessed the degree of overlap in the ecological niche of each taxon. To do this, we sampled and genotyped 305 sparrows from 34 marshes across the hybrid zone and from adjacent regions. We used linear regression to test for associations between marsh characteristics and the distribution of pure and admixed sparrows. We found a positive correlation between genotype and environmental variables with a patchy distribution of genotypes and habitats across the hybrid zone. Ecological niche models suggest that the hybrid niche was more similar to that of A. nelsoni and habitat suitability was influenced strongly by distance from coastline. Our results support a mosaic model of hybrid zone maintenance, suggesting a role for local environmental features in shaping the distribution and frequency of pure species and hybrids across space.     

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.     

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome‐scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well‐defined categories of locus types (e.g., neutral vs. highly selected). Our multi‐locus analysis reveals cross‐genomic variation in selective constraints on gene flow and locus‐specific flexibility in the permeability of the interspecies membrane.     

Facultative pupal mating in Heliconius erato: Implications for mate choice,female preference,and speciation

Mating systems have broad impacts on how sexual selection and mate choice operate within a species, but studies of mating behavior in the laboratory may not reflect how these processes occur in the wild. Here, we examined the mating behavior of the neotropical butterfly Heliconius erato in the field by releasing larvae and virgin females and observing how they mated. H. erato is considered a pupal‐mating species (i.e., males mate with females as they emerge from the pupal case). However, we observed only two teneral mating events, and experimentally released virgins were almost all mated upon recapture. Our study confirms the presence of some pupal‐mating behavior in H. erato, but suggests that adult mating is likely the prevalent mating strategy in this species. These findings have important implications for the role of color pattern and female mate choice in the generation of reproductive isolation in this diverse genus.     

Sharp acoustic boundaries across an altitudinal avian hybrid zone despite asymmetric introgression

Birdsong is a sexually selected trait that could play an important evolutionary role when related taxa come into secondary contact. Many songbird species, however, learn their songs through copying one or more tutors, which complicates the evolutionary outcome of such contact. Two subspecies of a presumed vocal learner, the grey‐breasted wood‐wren (Henicorhina leucophrys), replace each other altitudinally across the western slope of the Ecuadorian Andes. These subspecies are morphologically very similar, but show striking differences in their song. We examined variation in acoustic traits and genetic composition across the altitudinal range covered by both subspecies and between two allopatric populations. The acoustic boundary between the subspecies was found to be highly abrupt across a narrow elevational range with virtually no evidence of song convergence. Mixed singing and use of hetero‐subspecific song occurred in the contact zone and was biased towards the use of leucophrys song types. Hetero‐subspecific song copying by hilaris and not by leucophrys reflected a previously found asymmetric pattern of response to song playback. Using amplified fragment length polymorphisms (AFLP) markers, we detected hybridization in the contact zone and asymmetric introgression in parapatric populations, with more leucophrys alleles present in hilaris populations than vice versa. This pattern may be a trail of introgression due to upslope displacement of leucophrys by hilaris. Our data suggest that song learning may impact speciation and hybridization in contrasting ways at different spatial scales: although learning may speed up population divergence in songs, thereby enhancing assortative mating and reducing gene flow, it may at a local level also lead to the copying of heterospecific songs, therefore allowing some level of hybridization and introgression.     

Emergence patterns of novelty in European vegetation assemblages over the past 15 000 years

Plant communities are not stable over time and biological novelty is predicted to emerge due to climate change, the introduction of exotic species and land‐use change. However, the rate at which this novelty may arise over longer time periods has so far received little attention. We reconstruct the emergence of novelty in Europe for a set of baseline conditions over the past 15 000 years to assess past rates of emergence and investigate underlying causes. The emergence of novelty is baseline specific and, during the early‐Holocene, was mitigated by the rapid spread of plant taxa. Although novelty generally increases as a function of time, climate and human‐induced landscape changes contributed to a non‐linear post‐glacial trajectory of novelty with jumps corresponding to periods of rapid changes. Emergence of novelty accelerated during the past 1000 years. Historical cultural landscapes experienced a faster novelty development due to the contribution from anthropogenic land‐cover changes.     

Embryogenesis of Heliconius erato (Lepidoptera,Nymphalidae): A contribution to the anatomical development of an evo‐devo model organism

This study reports on the embryogenesis of Heliconius erato phyllis between blastoderm formation and the prehatching larval stage. Syncytial blastoderm formation occurred approximately 2 h after egg laying (AEL) and at about 4 h, the cellular blastoderm was formed. The germ band arose from the entire length of the blastoderm, and rapidly became compacted occupying approximately two‐thirds of the egg length. At about 7 h AEL, protocephalon and protocorm differentiation occurred. Continued proliferation of the germ band was followed by penetration into the yolk mass, forming a C‐shaped embryo at about 10 h. Approximately 12 h AEL, the gnathal, thoracic and abdominal segments became visible. The primordium of the mouthparts and thoracic legs formed as paired evaginations, while the prolegs formed as paired lobes. At about 30 h, the embryo reversed dorsoventrally. Approximately 32 h AEL, the protocephalon and gnathal segments fused, shifting the relative position of the rudimentary appendages in this region. At about 52 h, the embryo was U‐shaped in lateral view and at approximately 56 h, the bristles began evagination from the larval cuticle. Larvae hatched at about 72 h. We found that H. erato phyllis followed an embryonic pattern consistent with long‐germ embryogenesis. Thus, we believe that H. erato phyllis should be classified as a long‐germ lepidopteran. The study of H. erato phyllis embryogenesis provided a structural glimpse into the morphogenetic events that occur in the Heliconius egg period. This study could help future molecular approaches to understanding the evolution of Heliconius development.     

hzar: hybrid zone analysis using an R software package

We present a new software package (hzar ) that provides functions for fitting molecular genetic and morphological data from hybrid zones to classic equilibrium cline models using the Metropolis–Hastings Markov chain Monte Carlo (MCMC) algorithm. The software applies likelihood functions appropriate for different types of data, including diploid and haploid genetic markers and quantitative morphological traits. The modular design allows flexibility in fitting cline models of varying complexity. To facilitate hypothesis testing, an autofit function is included that allows automated model selection from a set of nested cline models. Cline parameter values, such as cline centre and cline width, are estimated and may be compared statistically across clines. The package is written in the R language and is available through the Comprehensive R Archive Network (CRAN; http://cran.r-project.org/ ). Here, we describe hzar and demonstrate its use with a sample data set from a well‐studied hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus). Comparisons of our results with previously published results for this hybrid zone validate the hzar software. We extend analysis of this hybrid zone by fitting additional models to molecular data where appropriate.     

Stabilization of a salamander moving hybrid zone

When related species meet upon postglacial range expansion, hybrid zones are frequently formed. Theory predicts that such zones may move over the landscape until equilibrium conditions are reached. One hybrid zone observed to be moving in historical times (1950–1979) is that of the pond‐breeding salamanders Triturus cristatus and Triturus marmoratus in western France. We identified the ecological correlates of the species hybrid zone as elevation, forestation, and hedgerows favoring the more terrestrial T. marmoratus and pond density favoring the more aquatic T. cristatus. The past movement of the zone of ca. 30 km over three decades has probably been driven by the drastic postwar reduction of the “bocage” hedgerow landscape, favoring T. cristatus over T. marmoratus. No further hybrid zone movement was observed from 1979 to the present. To explain the changing dynamics of the hybrid zone, we propose that it stalled, either because an equilibrium was found at an altitude of ca. 140 m a.s.l. or due to pond loss and decreased population densities. While we cannot rule out the former explanation, we found support for the latter. Under agricultural intensification, ponds in the study area are lost at an unprecedented rate of 5.5% per year, so that remaining Triturus populations are increasingly isolated, hampering dispersal and further hybrid zone movement.