Inferences from a rapidly moving hybrid zone

Genetic traces of hybrid zone movement across a fragmented habitat

Theoretical and empirical studies suggest that the structure and position of hybrid zones can change over time. Evidence for moving hybrid zones has been directly inferred by repeated sampling over time, or indirectly through the detection of genetic footprints left by the receding species and the resulting asymmetric patterns of introgression across markers. We here investigate a hybrid zone formed by two subspecies of the Iberian golden-striped salamander, Chioglossa lusitanica, using a panel of 35 nuclear loci (31 SNPs and 4 allozymes) and one mitochondrial locus in a transect in central Portugal. We found concordant and coincident clines for most of the nuclear loci (n=22, 63%), defining a narrow hybrid zone of ca. 6 km wide, with the centre positioned ca. 15 km south of the Mondego river. Asymmetric introgression was observed at another 14 loci. Their clines are displaced towards the north, with positions located either close to the Mondego river (n=6), or further northwards (n=8). We interpret these profiles as genetic traces of the southward displacement of C. lusitanica lusitanica by C. l. longipes over the wider Mondego river valley. We noted the absence of significant linkage disequilibrium and we inferred low levels of effective selection per locus against hybrids, suggesting that introgression in the area of species replacement occurred under a neutral diffusion process. A species distribution model suggests that the C. lusitanica hybrid zone coincides with a narrow corridor of fragmented habitat. From the position of the displaced clines, we infer that patches of locally suitable habitat trapped some genetic variants that became disassociated from the southward moving hybrid zone. This study highlights the influence of habitat availability on hybrid zone movement.

Weak coupling among barrier loci and waves of neutral and adaptive introgression across an expanding hybrid zone

Hybridization can serve as an evolutionary stimulus, but we have little understanding of introgression at early stages of hybrid zone formation. We analyze reproductive isolation and introgression between a range‐limited and a widespread species. Reproductive barriers are estimated based on differences in flowering time, ecogeographic distributions, and seed set from crosses. We find an asymmetrical mating barrier due to cytonuclear incompatibility that is consistent with observed clusters of coincident and concordant tension zone clines (barrier loci) for mtDNA haplotypes and nuclear SNPs. These groups of concordant clines are spread across the hybrid zone, resulting in weak coupling among barrier loci and extensive introgression. Neutral clines had nearly equal introgression into both species’ ranges, whereas putative cases of adaptive introgression had exceptionally wide clines with centers shifted toward one species. Analyses of cline shape indicate that secondary contact was initiated within the last 800 generations with the per‐generation dispersal between 200 and 400 m, and provide some of the first estimates of the strength of selection required to account for observed levels of adaptive introgression. The weak species boundary between these species appears to be in early stages of dissolution, and ultimately will precipitate genetic swamping of the range‐limited species.

Displaced clines in an avian hybrid zone (Thamnophilidae: Rhegmatorhina) within an Amazonian interfluve

Secondary contact between species often results in the formation of a hybrid zone, with the eventual fates of the hybridizing species dependent on evolutionary and ecological forces. We examine this process in the Amazon Basin by conducting the first genomic and phenotypic characterization of the hybrid zone formed after secondary contact between two obligate army-ant-followers: the White-breasted Antbird (Rhegmatorhina hoffmannsi) and the Harlequin Antbird (Rhegmatorhina berlepschi). We found a major geographic displacement (∼120 km) between the mitochondrial and nuclear clines, and we explore potential hypotheses for the displacement, including sampling error, genetic drift, and asymmetric cytonuclear incompatibilities. We cannot exclude roles for sampling error and genetic drift in contributing to the discordance; however, the data suggest expansion and unidirectional introgression of hoffmannsi into the distribution of berlepschi.

Four Species Linked by Three Hybrid Zones: Two Instances of Repeated Hybridization in One Species Group (Genus Liolaemus)

Hybridization is an evolutionary process that can generate diverse outcomes, such as reinforcing species boundaries, generating new species, or facilitating the introgression of locally-adapted alleles into new genomic backgrounds. Liolaemus is a highly diverse clade of South American lizards with ~260 species and as many as ten new species are described each year. Previous Liolaemus studies have detected gene flow and introgression among species using phylogenetic network methods and/or through comparisons of nuclear and mitochondrial DNA patterns, yet no study has systematically studied hybrid zones between Liolaemus species. Here, we compared three hybrid zones between four species in the Liolaemus fitzingerii group of lizards in Central Argentina where two species, L. melanops and L. xanthoviridis, each hybridize with two other species (L. shehuen and L. fitzingerii). We sampled three transects that were each ~120 km in length and sequenced both mitochondrial and genome-wide SNP data for 267 individuals. In our analyses of nuclear DNA, we also compared bi-allelic SNPs to phased alleles (50 bp RAD loci). Population structure analyses confirmed that boundaries separating species are sharp, and all clines are <65 km wide. Cline center estimates were consistent between SNPs and phased alleles, but cline width estimates were significantly different with the SNPs producing wider estimates. The mitochondrial clines are narrower and shifted 4–20 km southward in comparison to the nuclear clines in all three hybrid zones, indicating that either each of the species has sex-biased dispersal (males northward or females southward), the population densities are unequal, or that the hybrid zones are moving north over time. These comparisons indicate that some patterns of hybridization are similar across hybrid zones (mtDNA clines all narrower and shifted to the south), whereas cline width is variable. Hybridization in the L. fitzingerii group is common and geographically localized; further studies are needed to investigate whether hybrid zones act as hard species boundaries or promoters of speciation through processes such as reinforcement. Nonetheless, this study provides insights into both biotic and abiotic mechanisms helping to maintain species boundaries within the speciose Liolaemus system.

Mitochondrial DNA barcodes provide insight into the phylogeography and subspecies controversy in the widespread Neotropical white peacock butterfly, Anartia jatrophae (Nymphalidae: Nymphalinae)

The validity of subspecies designations in the common and wide-ranging Neotropical white peacock butterfly, Anartia jatrophae, has been debated for decades and remains an unsettled and contentious taxonomic issue among lepidopterists. Originally described by Linnaeus in the mid-18th century from specimens obtained from northern South America (Suriname), subsequent authors proposed a variety of subspecies names based on differences in adult external morphology among geographical populations. Many of these differences, however, were subsequently found to occur seasonally within populations, leading some to conclude that only a single polymorphic species should be recognized. Here, we have analysed both new and publicly available mitochondrial DNA barcodes, obtained from specimens collected from southern USA to northern Argentina, to assess whether they could provide insight into this long-standing controversy. Our molecular analyses, using a combination of character-based (nucleotide composition), population genetic and phylogenetic approaches, indicated the presence of at least four distinct genetic lineages that we suggest are distinct at the subspecies level, namely A. j. jatrophae, A. j. luteipicta, A. j. saturata and A. j. semifusca. Justification for these assignments and the proposed geographical distribution of each subspecies within the Americas are discussed.

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental-like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

Tracing the footprints of a moving hybrid zone under a demographic history of speciation with gene flow

A lack of optimal gene combinations, as well as low levels of genetic diversity, is often associated with the formation of species range margins. Conservation efforts rely on predictive modelling using abiotic variables and assessments of genetic diversity to determine target species and populations for controlled breeding, germplasm conservation and assisted migration. Biotic factors such as interspecific competition and hybridization, however, are largely ignored, despite their prevalence across diverse taxa and their role as key evolutionary forces. Hybridization between species with well-developed barriers to reproductive isolation often results in the production of offspring with lower fitness. Generation of novel allelic combinations through hybridization, however, can also generate positive fitness consequences. Despite this possibility, hybridization-mediated introgression is often considered a threat to biodiversity as it can blur species boundaries. The contribution of hybridization towards increasing genetic diversity of populations at range margins has only recently gathered attention in conservation studies. We assessed the extent to which hybridization contributes towards range dynamics by tracking spatio-temporal changes in the central location of a hybrid zone between two recently diverged species of pines: Pinus strobiformis and P. flexilis. By comparing geographic cline centre estimates for global admixture coefficient with morphological traits associated with reproductive output, we demonstrate a northward shift in the hybrid zone. Using a combination of spatially explicit, individual-based simulations and linkage disequilibrium variance partitioning, we note a significant contribution of adaptive introgression towards this northward movement, despite the potential for differences in regional population size to aid hybrid zone movement. Overall, our study demonstrates that hybridization between recently diverged species can increase genetic diversity and generate novel allelic combinations. These novel combinations may allow range margin populations to track favourable climatic conditions or facilitate adaptive evolution to ongoing and future climate change.

Unsupervised machine learning reveals mimicry complexes in bumblebees occur along a perceptual continuum

Müllerian mimicry theory states that frequency-dependent selection should favour geographical convergence of harmful species onto a shared colour pattern. As such, mimetic patterns are commonly circumscribed into discrete mimicry complexes, each containing a predominant phenotype. Outside a few examples in butterflies, the location of transition zones between mimicry complexes and the factors driving mimicry zones has rarely been examined. To infer the patterns and processes of Müllerian mimicry, we integrate large-scale data on the geographical distribution of colour patterns of social bumblebees across the contiguous United States and use these to quantify colour pattern mimicry using an innovative, unsupervised machine-learning approach based on computer vision. Our data suggest that bumblebees exhibit geographically clustered, but sometimes imperfect colour patterns, and that mimicry patterns gradually transition spatially rather than exhibit discrete boundaries. Additionally, examination of colour pattern transition zones of three comimicking, polymorphic species, where active selection is driving phenotype frequencies, revealed that their transition zones differ in location within a broad region of poor mimicry. Potential factors influencing mimicry transition zone dynamics are discussed.

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

Hybrid zones have long been of interest to biologists as natural laboratories where we can gain insight into the processes of adaptation and speciation. Repeated sampling of individual hybrid zones has been particularly useful in elucidating the dynamic balance between selection and dispersal that maintains most hybrid zones. Here, we revisit a hybrid zone between Heliconius erato butterflies in Panamá for a third time over more than 30 years. We combine a novel Bayesian extension of stepped-cline hybrid zone models with environmental data to understand the genetic and environmental causes of cline dynamics in this species. The cline has continued to move west, likely due to dominance drive, but has slowed and broadened. Environmental analyses suggest that widespread deforestation in Panamá could be leading to decreased avian predation and relaxed selection, causing the observed changes in cline dynamics.

Assortative Mating by an Obliquely Transmitted Local Cultural Trait Promotes Genetic Divergence: A Model

The effect of learned culture (e.g., birdsong dialects and human languages) on genetic divergence is unclear. Previous theoretical research suggests that because oblique learning allows phenotype transmission from individuals with no offspring to an unrelated individual in the next generation, the effect of sexual selection on the learned trait is masked. However, I propose that migration and spatially constrained learning can form statistical associations between cultural and genetic traits, which may allow selection on the cultural traits to indirectly affect the genetic traits. Here, I build a population genetic model that allows such statistical associations to form and find that sexual selection and divergent selection on the cultural trait can indeed help maintain genetic divergence through such statistical associations, while selection against genetic hybrids does not affect cultural trait divergence. Furthermore, I find that even when the cultural trait changes over time due to drift and mutation, it can still help maintain genetic divergence. These results suggest the role of obliquely transmitted traits in evolution may be underrated, and the lack of one-to-one associations between cultural and genetic traits may not be sufficient to disprove the role of culture in genetic divergence.

Climate-mediated hybrid zone movement revealed with genomics, museum collection, and simulation modeling

Climate-mediated changes in hybridization will dramatically alter the genetic diversity, adaptive capacity, and evolutionary trajectory of interbreeding species. Our ability to predict the consequences of such changes will be key to future conservation and management decisions. Here we tested through simulations how recent warming (over the course of a 32-y period) is affecting the geographic extent of a climate-mediated developmental threshold implicated in maintaining a butterfly hybrid zone (Papilio glaucus and Papilio canadensis; Lepidoptera: Papilionidae). These simulations predict a 68-km shift of this hybrid zone. To empirically test this prediction, we assessed genetic and phenotypic changes using contemporary and museum collections and document a 40-km northward shift of this hybrid zone. Interactions between the two species appear relatively unchanged during hybrid zone movement. We found no change in the frequency of hybridization, and regions of the genome that experience little to no introgression moved largely in concert with the shifting hybrid zone. Model predictions based on climate scenarios predict this hybrid zone will continue to move northward, but with substantial spatial heterogeneity in the velocity (55-144 km/1 °C), shape, and contiguity of movement. Our findings suggest that the presence of nonclimatic barriers (e.g., genetic incompatibilities) and/or nonlinear responses to climatic gradients may preserve species boundaries as the species shift. Further, we show that variation in the geography of hybrid zone movement could result in evolutionary responses that differ for geographically distinct populations spanning hybrid zones, and thus have implications for the conservation and management of genetic diversity.

Origin and cross-century dynamics of an avian hybrid zone

Background

Characterizations of the dynamics of hybrid zones in space and time can give insights about traits and processes important in population divergence and speciation. We characterized a hybrid zone between tanagers in the genus Ramphocelus (Aves, Thraupidae) located in southwestern Colombia. We evaluated whether this hybrid zone originated as a result of secondary contact or of primary differentiation, and described its dynamics across time using spatial analyses of molecular, morphological, and coloration data in combination with paleodistribution modeling.

Results

Models of potential historical distributions based on climatic data and genetic signatures of demographic expansion suggested that the hybrid zone likely originated following secondary contact between populations that expanded their ranges out of isolated areas in the Quaternary. Concordant patterns of variation in phenotypic characters across the hybrid zone and its narrow extent are suggestive of a tension zone, maintained by a balance between dispersal and selection against hybrids. Estimates of phenotypic cline parameters obtained using specimens collected over nearly a century revealed that, in recent decades, the zone appears to have moved to the east and to higher elevations, and may have become narrower. Genetic variation was not clearly structured along the hybrid zone, but comparisons between historical and contemporary specimens suggested that temporal changes in its genetic makeup may also have occurred.

Conclusions

Our data suggest that the hybrid zone likey resulted from secondary contact between populations. The observed changes in the hybrid zone may be a result of sexual selection, asymmetric gene flow, or environmental change.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12862-017-1096-7) contains supplementary material, which is available to authorized users.

A genomic footprint of hybrid zone movement in crested newts

Speciation typically involves a stage in which species can still exchange genetic material. Interspecific gene flow is facilitated by the hybrid zones that such species establish upon secondary contact. If one member of a hybridizing species pair displaces the other, their hybrid zone would move across the landscape. Although theory predicts that moving hybrid zones quickly stagnate, hybrid zones tracked over one or a few decades do not always follow such a limitation. This suggests that hybrid zones have the potential to traverse considerable distances over extended periods of time. When hybrid zones move, introgression is predicted to result in biased gene flow of selectively neutral alleles, from the receding species into the advancing species. We test for such a genomic footprint of hybrid zone movement in a pair of crested newt species (genus Triturus) for which we have a priori support for westward hybrid zone movement. We perform a multilocus phylogeographical survey and conduct Bayesian clustering analysis, estimation of ancestry and heterozygosity, and geographical cline analysis. In a 600 km wide area east of the present day hybrid zone a genomic footprint constitutes empirical evidence consistent with westward hybrid zone movement. The crested newt case suggests that hybrid zone movement can occur over an extensive span of time and space. Inferring hybrid zone movement provides fundamental insight into historical biogeography and the speciation process, and we anticipate that hybrid zones will prove to be far more mobile than currently appreciated.

Changing climate mediates sapsucker (Aves: Sphyrapicus) hybrid zone movement

Hybrid zones, where two divergent taxa meet and interbreed, offer unique opportunities to investigate how climate contributes to reproductive isolation between closely related taxa and how these taxa may respond to climatic changes. Red-naped (Sphyrapicus nuchalis) and Red-breasted (Sphyrapicus ruber) sapsuckers (Aves: Picidae) hybridize along a narrow contact zone that stretches from northern California to British Columbia. The hybrid zone between these species has been studied extensively for more than 100 years and represents an excellent system for investigations of the evolution of reproductive isolation. Shifts in the proportions of phenotypes at hybrid localities since 1910 that were inferred using specimens from museum collections were confirmed using species distribution models. We predicted the historical, current, and future distributions of parental and hybrid sapsuckers using Random Forests models to quantify how climate change is affecting hybrid zone movement in the Pacific Northwest. We found observed distribution shifts of parental sapsuckers were likely the result of climate change over the past 100 years, with these shifts predicted to continue for both sapsuckers over the next 80 years. We found Red-breasted Sapsuckers are predicted to continue to expand, while Red-naped Sapsuckers are predicted to contract substantially under future climate scenarios. As a result of the predicted changes, the amount of overlap in the distribution of these sapsuckers may decrease. Using hybrid phenotypes, we found the climate niche occupied by the hybrid zone is predicted to disappear under future conditions. The disappearance of this climate niche where the two parental species come into contact and hybridize may lead to a substantial reduction in genetic introgression. Understanding the impacts of global climate change on hybrid zones may help us to better understand how speciation has been shaped by climate in the past, as well as how evolution may respond to climate change in the future.

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.

Rapid movement and instability of an invasive hybrid swarm

Unstable hybrid swarms that arise following the introduction of non‐native species can overwhelm native congeners, yet the stability of invasive hybrid swarms has not been well documented over time. Here, we examine genetic variation and clinal stability across a recently formed hybrid swarm involving native blacktail shiner (Cyprinella venusta) and non‐native red shiner (C. lutrensis) in the Upper Coosa River basin, which is widely considered to be a global hot spot of aquatic biodiversity. Examination of phenotypic, multilocus genotypic, and mitochondrial haplotype variability between 2005 and 2011 revealed that the proportion of hybrids has increased over time, with more than a third of all sampled individuals exhibiting admixture in the final year of sampling. Comparisons of clines over time indicated that the hybrid swarm has been rapidly progressing upstream, but at a declining and slower pace than rates estimated from historical collection records. Clinal comparisons also showed that the hybrid swarm has been expanding and contracting over time. Additionally, we documented the presence of red shiner and hybrids farther downstream than prior studies have detected, which suggests that congeners in the Coosa River basin, including all remaining populations of the threatened blue shiner (Cyprinella caerulea), are at greater risk than previously thought.

Hybrid zones: windows on climate change

Defining the impacts of anthropogenic climate change on biodiversity and species distributions is currently a high priority. Niche models focus primarily on predicted changes in abiotic factors; however, species interactions and adaptive evolution will impact the ability of species to persist in the face of changing climate. Our review focuses on the use of hybrid zones to monitor responses of species to contemporary climate change. Monitoring hybrid zones provides insight into how range boundaries shift in response to climate change by illuminating the combined effects of species interactions and physiological sensitivity. At the same time, the semipermeable nature of species boundaries allows us to document adaptive introgression of alleles associated with response to climate change.

Stable Heliconius butterfly hybrid zones are correlated with a local rainfall peak at the edge of the Amazon basin

Multilocus clines between Müllerian mimetic races of Heliconius butterflies provide a classic example of the maintenance of hybrid zones and their importance in speciation. Concordant hybrid zones in the mimics Heliconius erato and H. melpomene in northern Peru were carefully documented in the 1980s, and this prior work now permits a historical analysis of the movement or stasis of the zones. Previous work predicted that these zones might be moving toward the Andes due to selective asymmetry. Extensive deforestation and climate change might also be expected to affect the positions and widths of the hybrid zones. We show that the positions and shapes of these hybrid zones have instead remained remarkably stable between 1985 and 2012. The stability of this interaction strongly implicates continued selection, rather than neutral mixing following secondary contact. The stability of cline widths and strong linkage disequilibria (gametic correlation coefficients Rmax = 0.35-0.56 among unlinked loci) over 25 years suggest that mimetic selection pressures on each color pattern locus have remained approximately constant (s ≈ 0.13-0.40 per locus in both species). Exceptionally high levels of precipitation at the edge of the easternmost Andes may act as a population density trough for butterflies, trapping the hybrid zones at the foot of the mountains, and preventing movement. As such, our results falsify one prediction of the Pleistocene Refugium theory: That the ranges of divergent species or subspecies should be centered on regions characterized by maxima of rainfall, with hybrid zones falling in more arid regions between them.

A new barrier to dispersal trapped old genetic clines that escaped the Easter Microplate tension zone of the Pacific vent mussels

Comparative phylogeography of deep-sea hydrothermal vent species has uncovered several genetic breaks between populations inhabiting northern and southern latitudes of the East Pacific Rise. However, the geographic width and position of genetic clines are variable among species. In this report, we further characterize the position and strength of barriers to gene flow between populations of the deep-sea vent mussel Bathymodiolus thermophilus. Eight allozyme loci and DNA sequences of four nuclear genes were added to previously published sequences of the cytochrome c oxidase subunit I gene. Our data confirm the presence of two barriers to gene flow, one located at the Easter Microplate (between 21°33′S and 31°S) recently described as a hybrid zone, and the second positioned between 7°25′S and 14°S with each affecting different loci. Coalescence analysis indicates a single vicariant event at the origin of divergence between clades for all nuclear loci, although the clines are now spatially discordant. We thus hypothesize that the Easter Microplate barrier has recently been relaxed after a long period of isolation and that some genetic clines have escaped the barrier and moved northward where they have subsequently been trapped by a reinforcing barrier to gene flow between 7°25′S and 14°S.

Range dynamics, rather than convergent selection, explain the mosaic distribution of red-winged blackbird phenotypes

Geographic distributions of genetic and phenotypic characters can illuminate historical evolutionary processes. In particular, mosaic distributions of phenotypically similar populations can arise from parallel evolution or from irregular patterns of dispersal and colonization by divergent forms. Two phenotypically divergent forms of the red-winged blackbird (Agelaius phoeniceus) show a mosaic phenotypic distribution, with a "bicolored" form occurring disjunctly in California and Mexico. We analyzed the relationships among these bicolored populations and neighboring typical populations, using ∼600 bp of mitochondrial DNA sequence data and 10 nuclear short tandem repeat loci. We find that bicolored populations, although separated by ∼3000 km, are genetically more similar to one other than they are to typical populations separated by ∼400 km. We also find evidence of ongoing gene flow among populations, including some evidence of asymmetric gene flow. We conclude that the current distribution of bicolored forms represents incomplete speciation, where recent asymmetric hybridization with typical A. phoeniceus is dividing the range of a formerly widespread bicolored form. This hypothesis predicts that bicolored forms may suffer extinction by hybridization. Future work will use fine-scaled geographical sampling and nuclear sequence data to test for hybrid origins of currently typical populations and to more precisely quantify the directionality of gene flow.

Parental ploidy strongly affects offspring fitness in heteroploid crosses among three cytotypes of autopolyploid Jacobaea carniolica (Asteraceae)

Reproductive interactions among cytotypes in their contact zones determine whether these cytotypes can co-exist and form stable contact zones or not. In autopolyploids, heteroploid cross-compatibilities might depend on parental ploidy, but tests of this hypothesis in autopolyploid systems with more than two ploidies are lacking. Here, we study Jacobaea carniolica, which comprises diploid, tetraploid, and hexaploid individuals regularly forming contact zones. Seeds obtained from in situ cross-pollinations within and among cytotypes were subjected to DNA flow cytometry and greenhouse germination experiments. Hybrid fitness and parental effects on hybrid fitness were tested with regression models comparing fitness parameters of early life stages. Irrespective of the direction of crosses, seed viability and seedling survival in diploid-polyploid crosses were substantially lower than in tetraploid-hexaploid crosses. In contrast, seedling growth traits indicated neither transgressive character expression nor any selection against hybrid offspring. Congruent with a model of genome dosage effects, these traits differed between reciprocal crosses, especially of diploids and tetraploids, where trait values resembled those of the maternal parent. The strong effect of parental ploidy on offspring fitness in heteroploid crosses may cause contact zones involving exclusively polyploid cytotypes to be less stable over longer terms than those involving diploids and polyploids.

When history repeats itself: exploring the genetic architecture of host-plant adaptation in two closely related lepidopteran species

The genus Ostrinia includes two allopatric maize pests across Eurasia, namely the European corn borer (ECB, O. nubilalis) and the Asian corn borer (ACB, O. furnacalis). A third species, the Adzuki bean borer (ABB, O. scapulalis), occurs in sympatry with both the ECB and the ACB. The ABB mostly feeds on native dicots, which probably correspond to the ancestral host plant type for the genus Ostrinia. This situation offers the opportunity to characterize the two presumably independent adaptations or preadaptations to maize that occurred in the ECB and ACB. In the present study, we aimed at deciphering the genetic architecture of these two adaptations to maize, a monocot host plant recently introduced into Eurasia. To this end, we performed a genome scan analysis based on 684 AFLP markers in 12 populations of ECB, ACB and ABB. We detected 2 outlier AFLP loci when comparing French populations of the ECB and ABB, and 9 outliers when comparing Chinese populations of the ACB and ABB. These outliers were different in both countries, and we found no evidence of linkage disequilibrium between any two of them. These results suggest that adaptation or preadaptation to maize relies on a different genetic architecture in the ECB and ACB. However, this conclusion must be considered in light of the constraints inherent to genome scan approaches and of the intricate evolution of adaptation and reproductive isolation in the Ostrinia spp. complex.

Oscillation of an anuran hybrid zone: morphological evidence spanning 50 years

Background

The hybrid zone between the primarily forest-dwelling American toad, Anaxyrus americanus, and the prairie-adapted Canadian toad, A. hemiophrys, in southeastern Manitoba is known to have shifted its position during the past 50 years. Hybrid zones are areas of interbreeding between species and their movement across a landscape should reflect their underlying dynamics and environmental change. However, empirical demonstrations of hybrid zone movements over long periods of time are rare. This hybrid zone is dominated by individuals of intermediate morphology and genetic composition. We sought to determine if it had continued to move and if that movement was associated with shifts in habitat, as predicted.

Methodology/Principle Findings

We used variation in the toads’ most diagnostic morphological feature, the separation between their interorbital cranial crests, to determine the geographic position of the hybrid zone center at four times between 1960 and 2009 using maximum likelihood methods. The hybrid zone center moved west by 38 km over 19 years and then east again by 10 km over the succeeding 29 years. The position of the hybrid zone did not track either the direction or the magnitude of movement of the forest-prairie habitat transition over the same time period.

Conclusions/Significance

This is the first reported evidence of oscillation in the position of a hybrid zone. The back and forth movement indicates that neither species maintains a selective advantage over the other in the long term. However, the movement of the hybrid zone was not bounded by the breadth of the habitat transition. Its oscillation suggests that the hybrid zone is better described as being elastically tethered to the habitat transition.

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.

Discordant introgression in a rapidly expanding hybrid swarm

The erosion of species boundaries can involve rapid evolutionary change. Consequently, many aspects of the process remain poorly understood, including the formation, expansion, and evolution of hybrid swarms. Biological invasions involving hybridization present exceptional opportunities to study the erosion of species boundaries because timelines of interactions and outcomes are frequently well known. Here, we examined clinal variation across codominant and maternally inherited genetic markers as well as phenotypic traits to characterize the expansion and evolution of a hybrid swarm between native Cyprinella venusta and invasive Cyprinella lutrensis minnows. Discordant introgression of phenotype, microsatellite multilocus genotype, and mtDNA haplotype indicates that the observable expansion of the C. venusta × C. lutrensis hybrid swarm is a false invasion front. Both parental and hybrid individuals closely resembling C. lutrensis are numerically dominant in the expansion wake, indicating that the non-native parental phenotype may be selectively favored. These findings show that cryptic introgression can extend beyond the phenotypic boundaries of hybrid swarms and that hybrid swarms likely expand more rapidly than can be documented from phenotypic variation alone. Similarly, dominance of a single parental phenotype following an introduction event may lead to instances of species erosion being mistaken for species displacement without hybridization.

Relationships between spatio-temporal environmental and genetic variation reveal an important influence of exogenous selection in a pupfish hybrid zone

The importance of exogenous selection in a natural hybrid zone between the pupfishes Cyprinodon atrorus and Cyprinodon bifasciatus was tested via spatio-temporal analyses of environmental and genetic change over winter, spring and summer for three consecutive years. A critical influence of exogenous selection on hybrid zone regulation was demonstrated by a significant relationship between environmental (salinity and temperature) and genetic (three diagnostic nuDNA loci) variation over space and time (seasons) in the Rio Churince system, Cuatro Ciénegas, Mexico. At sites environmentally more similar to parental habitats, the genetic composition of hybrids was stable and similar to the resident parental species, whereas complex admixtures of parental and hybrid genotypic classes characterized intermediate environments, as did the greatest change in allelic and genotypic frequencies across seasons. Within hybrids across the entire Rio Churince system, seasonal changes in allelic and genotypic frequencies were consistent with results from previous reciprocal transplant experiments, which showed C. bifasciatus to suffer high mortality (75%) when exposed to the habitat of C. atrorus in winter (extreme temperature lows and variability) and summer (abrupt salinity change and extreme temperature highs and variability). Although unconfirmed, the distributional limits of C. atrorus and C. atrorus-like hybrids appear to be governed by similar constraints (predation or competition). The argument favouring evolutionary significance of hybridization in animals is bolstered by the results of this study, which links the importance of exogenous selection in a contemporary hybrid zone between C. atrorus and C. bifasciatus to previous demonstration of the long-term evolutionary significance of environmental variation and introgression on the phenotypic diversification Cuatro Ciénegas Cyprinodon.

Strong selection against hybrids maintains a narrow contact zone between morphologically cryptic lineages in a rainforest lizard

Phenotypically cryptic lineages comprise an important yet understudied part of biodiversity; in particular, we have much to learn about how these lineages are formed and maintained. To better understand the evolutionary significance of such lineages, we studied a hybrid zone between two morphologically cryptic phylogeographic lineages in the rainforest lizard, Lampropholis coggeri. Analyzing a multilocus genetic dataset through cline inference, individual-based methods and population measures of disequilibrium and using simulations to explore our genetic results in context of theoretical expectations, we inferred the processes maintaining this hybrid zone. We find that these lineages meet in a hybrid zone that is narrow (≈400 m) relative to inferred dispersal rate. Further, the hybrid zone exhibits substantial genetic disequilibrium and sharply coincident and largely concordant clines. Based on our knowledge about the region's biogeography, the species' natural history, and our simulation results, we suggest that strong selection against hybrids structures this system. As all clines show a relatively narrow range of introgression, we posit that this hybrid zone might not yet be in equilibrium. Nonetheless, our results clearly show that phylogeographic lineages can evolve substantial reproductive isolation without concomitant morphological diversification, suggesting that such lineages can constitute a significant component of evolutionary diversity.

Hybridization of mouse lemurs: different patterns under different ecological conditions

Background

Several mechanistic models aim to explain the diversification of the multitude of endemic species on Madagascar. The island's biogeographic history probably offered numerous opportunities for secondary contact and subsequent hybridization. Existing diversification models do not consider a possible role of these processes. One key question for a better understanding of their potential importance is how they are influenced by different environmental settings. Here, we characterized a contact zone between two species of mouse lemurs, Microcebus griseorufus and M. murinus, in dry spiny bush and mesic gallery forest that border each other sharply without intermediate habitats between them. We performed population genetic analyses based on mtDNA sequences and nine nuclear microsatellites and compared the results to a known hybrid zone of the same species in a nearby wide gradient from dry spiny bush over transitional forest to humid littoral forest.

Results

In the spiny-gallery system, Microcebus griseorufus is restricted to the spiny bush; Microcebus murinus occurs in gallery forest and locally invades the dryer habitat of its congener. We found evidence for bidirectional introgressive hybridization, which is closely linked to increased spatial overlap within the spiny bush. Within 159 individuals, we observed 18 hybrids with mitochondrial haplotypes of both species. Analyses of simulated microsatellite data indicate that we identified hybrids with great accuracy and that we probably underestimated their true number. We discuss short-term climatic fluctuations as potential trigger for the dynamic of invasion and subsequent hybridization. In the gradient hybrid zone in turn, long-term aridification could have favored unidirectional nuclear introgression from Microcebus griseorufus into M. murinus in transitional forest.

Conclusions

Madagascar's southeastern transitional zone harbors two very different hybrid zones of mouse lemurs in different environmental settings. This sheds light on the multitude of opportunities for the formation of hybrid zones and indicates an important influence of environmental factors on secondary contact and hybridization. Our findings suggest that hybridization could enhance the adaptability of mouse lemurs without necessarily leading to a loss of distinctiveness. They point to a potential role of hybridization in Madagascar's diversification history that requires further investigation.

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography-dependent asymmetrical dispersal or selection favouring one subspecies over the other.

Scanning the European corn borer (Ostrinia spp.) genome for adaptive divergence between host-affiliated sibling species

It has recently been shown that the European corn borer, a major pest of maize crops, is actually composed of two genetically differentiated and reproductively isolated taxa, which are found in sympatry over a wide geographical range in Eurasia. Each taxon is adapted to specific host plants: Ostrinia nubilalis feeds mainly on maize, while O. scapulalis feeds mainly on hop or mugwort. Here, we present a genome scan approach as a first step towards an integrated molecular analysis of the adaptive genomic divergence between O. nubilalis and O. scapulalis. We analysed 609 AFLP marker loci in replicate samples of sympatric populations of Ostrinia spp. collected on maize, hop and mugwort, in France. Using two genome scan methods based on the analysis of population differentiation, we found a set of 28 outlier loci that departed from the neutral expectation in one or the other method (of which a subset of 14 loci were common to both methods), which showed a significantly increased differentiation between O. nubilalis and O. scapulalis, when compared to the rest of the genome. A subset of 12 outlier loci were sequenced, of which 7 were successfully re-amplified as target candidate loci. Three of these showed homology with annotated lepidopteran sequences from public nucleotide databases.