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.     

Parallel pattern of differentiation at a genomic island shared between clinal and mosaic hybrid zones in a complex of cryptic seahorse lineages

Diverging semi‐isolated lineages either meet in narrow clinal hybrid zones, or have a mosaic distribution associated with environmental variation. Intrinsic reproductive isolation is often emphasized in the former and local adaptation in the latter, although both reduce gene flow between groups. Rarely are these two patterns of spatial distribution reported in the same study system. Here, we report that the long‐snouted seahorse Hippocampus guttulatus is subdivided into discrete panmictic entities by both types of hybrid zones. Along the European Atlantic coasts, a northern and a southern lineage meet in the southwest of France where they coexist in sympatry—i.e., in the same geographical zone—with little hybridization. In the Mediterranean Sea, two lineages have a mosaic distribution, associated with lagoon‐like and marine habitats. A fifth lineage was identified in the Black Sea. Genetic homogeneity over large spatial scales contrasts with isolation maintained in sympatry or close parapatry at a fine scale. A high variation in locus‐specific introgression rates provides additional evidence that partial reproductive isolation must be maintaining the divergence. We find that fixed differences between lagoon and marine populations in the Mediterranean Sea belong to the most differentiated SNPs between the two Atlantic lineages, against the genome‐wide pattern of structure that mostly follow geography. These parallel outlier SNPs cluster on a single chromosome‐wide island of differentiation. Since Atlantic lineages do not map to lagoon‐sea habitat variation, genetic parallelism at the genomic island suggests a shared genetic barrier contributes to reproductive isolation in contrasting contexts–i.e., spatial versus ecological. We discuss how a genomic hotspot of parallel differentiation could have evolved and become associated both with space and with a patchy environment in a single study system.     

Site-specific genetic divergence in parallel hybrid zones suggests nonallopatric evolution of reproductive barriers

The evolution of reproductive isolation in the presence of gene flow is supported by theoretical models but rarely by data. Empirical support might be gained from studies of parallel hybrid zones between interbreeding taxa. We analysed gene flow over two hybrid zones separating ecotypes of Littorina saxatilis to test the expectation that neutral genetic markers will show site-specific differences if barriers have evolved in situ. Distinct ecotypes found in contrasting shore habitats are separated by divergent selection and poor dispersal, but hybrid zones appear between them. Swedish islands formed by postglacial uplift 5000 years ago provide opportunities to assess genetic structure in a recently evolved system. Each island houses a discrete population containing subpopulations of different ecotypes. Hybrid zones between ecotypes may be a product of ecological divergence occurring on each island or a consequence of secondary overlap of ecotypes of allopatric origin that have spread among the islands. We used six microsatellite loci to assess gene flow and genetic profiles of hybrid zones on two islands. We found reduced gene flow over both hybrid zones, indicating the presence of local reproductive barriers between ecotypes. Nevertheless, subpopulations of different ecotypes from the same island were genetically more similar to each other than were subpopulations of the same ecotype from different islands. Moreover, neutral genetic traits separating the two ecotypes across hybrid zones were site-specific. This supports a scenario of in situ origin of ecotypes by ecological divergence and nonallopatric evolution of reproductive barriers.     

Speciation with gene flow: Evidence from a complex of alpine butterflies (Coenonympha,Satyridae)

Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

The maintenance of hybrid zones across a disturbance gradient

The parapatric distribution of genetically divergent lineages in hybrid zones can be maintained by ecological differences (dispersal-independent 'ecotonal' hybrid zones), by frequency- and density-dependent interference when they intermingle and mate (dispersal-dependent 'tension' hybrid zones), or by both processes acting together. One potentially important ecological factor that has received little theoretical attention is gradients in habitat disturbance. Such gradients may be particularly important in contact zones in which the interacting lineages differ in their sexual system (e.g., self-fertile versus obligately outcrossing) because self-fertility promotes the colonization of open patches. Here we use a spatially explicit metapopulation model to examine the dynamics of a dispersal-dependent ecotonal hybrid zone across a gradient in the rate of habitat disturbance, where competing lineages differ in their sexual system. We found that self-fertility promoted the maintenance of one lineage over its outcrossing counterpart at high extinction rates, predominantly because self-fertility confers reproductive assurance. Additionally, greater seed and pollen production promoted a lineage's persistence by reducing the seed fertility of its counterpart through hybridization. Our results draw attention to the joint effects of ecological and endogenous selection in regulating the location of hybrid zones. Our study also casts new light on the maintenance of the parapatric distribution of incompatible lineages of Spanish populations of the plant Mercurialis annua. In particular, we expect the rate of movement of a contact zone in eastern Spain to increase as it moves further south, contrary to earlier predictions.     

To hybridize or not to hybridize: what separates two genetic lineages of the Chalk‐hill Blue Polyommatus coridon (Lycaenidae,Lepidoptera) along their secondary contact zone throughout eastern Central Europe?

As a consequence of postglacial range expansion, hybrid zones evolved where different genetic lineages met. In this study, we analysed the Chalk‐hill Blue Polyommatus coridon all along the contact zone of two expansive lineages. This zone stretches from the sandy areas of north‐eastern Germany, along the mountain ranges of the German–Czech border and throughout the eastern Alps. We studied allozymes (19 loci) of 38 populations (1542 individuals) and compared these data sets against 15 populations of the western and 15 populations of the eastern lineages and found different degrees of hybridization. Thus, the calcareous regions of Thuringia and Sachsen‐Anhalt were mostly colonized by the western lineage. The middle mountain ranges between Bavaria and Bohemia represented a strong barrier blocking further expansion and thus completely impeding hybridization in this region. More intense hybridization was detected in the populations of the eastern Alps, especially in the north‐eastern part, where the Danube most probably acted as an expansion corridor for both lineages followed by intensive hybridization. In the south‐eastern Alps, hybrid populations were mostly detected in the easternmost parts and along the larger river valley of Drava and Mur; pure western populations dominated in the other areas of this region. These results show that the degree of hybridization along a contact zone is correlated with the ecological demands of a species and the regional physical geographic circumstances. This finding was proved for the Chalk‐hill Blue in our study but is also the most likely scenario in most animal and plant species.     

THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE

Studies of the genetics of hybrid zones can provide insight into the genomic architecture of species boundaries. By examining patterns of introgression of multiple loci across a hybrid zone, it may be possible to identify regions of the genome that have experienced selection. Here, we present a comparison of introgression in two replicate transects through the house mouse hybrid zone through central Europe, using data from 41 single nucleotide markers. Using both genomic and geographic clines, we found many differences in patterns of introgression between the two transects, as well as some similarities. We found that many loci may have experienced the effects of selection at linked sites, including selection against hybrid genotypes, as well as positive selection in the form of genotypes introgressed into a foreign genetic background. We also found many positive associations of conspecific alleles among unlinked markers, which could be caused by epistatic interactions. Different patterns of introgression in the two transects highlight the challenge of using hybrid zones to identify genes underlying isolation and raise the possibility that the genetic basis of isolation between these species may be dependent on the local population genetic make-up or the local ecological setting.     

Rapid genetic assimilation of native wall lizard populations (Podarcis muralis) through extensive hybridization with introduced lineages

The Common Wall Lizard (Podarcis muralis) has established more than 150 non-native populations in Central Europe, stemming from eight geographically distinct evolutionary lineages. While the majority of these introduced populations are found outside the native range, some of these populations also exist at the northern range margin in southwestern Germany. To (i) infer the level of hybridization in contact zones of alien and native lineages; and (ii) compare the genetic diversity among purebred introduced, native and hybrid populations, we used a combination of maternally inherited markers (mtDNA: cytb) and Mendelian markers (microsatellites). Our results suggest a rapid genetic assimilation of native populations by strong introgression from introduced lineages. Discordant patterns of mtDNA and nDNA variation within hybrid populations may be explained by directed mate choice of females towards males of alien lineages. In contrast to previous studies, we found a nonlinear relationship between genetic diversity and admixture level. The genetic diversity of hybrid populations was substantially higher than in introduced and native populations belonging to a single lineage, but rapidly reaching a plateau of high genetic diversity at an admixture level of two. However, even introduced populations with low founder sizes and from one source population retained moderate levels of genetic diversity and no evidence for a genetic bottleneck was found. The extent of introgression and the dominance of alien haplotypes in mixed populations indicate that introductions of non-native lineages represent a serious threat to the genetic integrity of native populations due to the rapid creation of hybrid swarms.     

Hybridization of Bombina bombina and B. variegata (Anura, Discoglossidae) at a sharp ecotone in western Ukraine: comparisons across transects and over time

Abstract Bombina bombina and B. variegata are two anciently diverged toad taxa that have adapted to different breeding habitats yet hybridize freely in zones of overlap where their parapatric distributions meet. Here, we report on a joint genetic and ecological analysis of a hybrid zone in the vicinity of Stryi in western Ukraine. We used five unlinked allozyme loci, two nuclear single nucleotide polymorphisms and a mitochondrial DNA haplotype as genetic markers. Parallel allele frequency clines with a sharp central step occur across a sharp ecotone, where transitions in aquatic habitat, elevation, and terrestrial vegetation coincide. The width of the hybrid zone, estimated as the inverse of the maximum gradient in allele frequency, is 2.3 km. This is the smallest of four estimates derived from different clinal transects across Europe. We argue that the narrow cline near Stryi is mainly due to a combination of habitat distribution and habitat preference. Adult toads show a preference for either ponds ( B. bombina ) or puddles ( B. variegata ), which is known to affect the distribution of genotypes within the hybrid zones. At Stryi, it should cause a reduction of the dispersal rate across the ecotone and thus narrow the cline. A detailed comparison of all five intensively studied Bombina transects lends support to the hypothesis that habitat distribution plus habitat preference can jointly affect the structure of hybrid zones and, ultimately, the resulting barriers to gene flow between differentiated gene pools. This study also represents a resampling of an area that was last studied more than 70 years ago. Our allele-frequency clines largely coincide with those that were described then on the basis of morphological variation. However, we found asymmetrical introgression of B. variegata genes into B. bombina territory along the bank of a river.     

Species Differentiation on a Dynamic Landscape: Shifts in Metapopulation Genetic Structure Using the Chronology of the Hawaiian Archipelago

Species formation during adaptive radiation often occurs in the context of a changing environment. The establishment and arrangement of populations, in space and time, sets up ecological and genetic processes that dictate the rate and pattern of differentiation. Here, we focus on how a dynamic habitat can affect genetic structure, and ultimately, differentiation among populations. We make use of the chronology and geographical history provided by the Hawaiian archipelago to examine the initial stages of population establishment and genetic divergence. We use data from a set of 6 spider lineages that differ in habitat affinities, some preferring low elevation habitats with a longer history of connection, others being more specialized for high elevation and/or wet forest, some with more general habitat affinities. We show that habitat preferences associated with lineages are important in ecological and genetic structuring. Lineages that have more restricted habitat preferences are subject to repeated episodes of isolation and fragmentation as a result of lava flows and vegetation succession. The initial dynamic set up by the landscape translates over time into discrete lineages. Further work is needed to understand how genetic changes interact with a changing set of ecological interactions amongst a shifting mosaic of landscapes to achieve species formation.     

Inter-species hybridization among Neotropical cats of the genus Leopardus, and evidence for an introgressive hybrid zone between L. geoffroyi and L. tigrinus in southern Brazil

Natural hybrid zones between distinct species have been reported for many taxa, but so far, few examples involve carnivores or Neotropical mammals in general. In this study, we employed mitochondrial DNA (mtDNA) sequences and nine microsatellite loci to identify and characterize a hybrid zone between two Neotropical felids, Leopardus geoffroyi and L. tigrinus, both of which are well-established species having diverged from each other c. 1 million years ago. These two felids are mostly allopatric throughout their ranges in South America, with a narrow contact zone that includes southern Brazil. We present strong evidence for the occurrence of hybridization between these species and identify at least 14 individuals (most of them originating from the geographical contact zone) exhibiting signs of interspecific genomic introgression. The genetic structure of Brazilian L. tigrinus populations seems to be affected by this introgression process, showing a gradient of differentiation from L. geoffroyi correlated with distance from the contact zone. We also corroborate and extend previous findings of hybridization between L. tigrinus and a third related felid, L. colocolo, leading to an unusual situation for a mammal, in which the former species contains introgressed mtDNA lineages from two distinct taxa in addition to its own.     

Strong premating divergence in a unimodal hybrid zone between two subspecies of the house mouse

Although selection against hybridization is expected to generate prezygotic divergence in unimodal hybrid zones, such a pattern has been seldom described. This study aims to better understand how prezygotic mechanisms may evolve in such zones. We investigated prezygotic divergence between populations of two subspecies of mice (Mus musculus musculus and M. m. domesticus) located at the edges of their unimodal hybrid zone in Denmark, and we developed an original multiple-population choice-test design, which allows assessment of within and between subspecies variation. Our study demonstrates that a strong assortative preference characterises one of the two subspecies (musculus) and that urinary signals are involved in this subspecies recognition. Taking into account the specific genetic and geographical characteristics of the Danish hybrid zone, we discuss the influence of the above pattern on its fate and the mechanisms that could have favoured this prezygotic divergence, among which the role of recombined populations constituting the core of the zone.     

The importance (or lack thereof) of niche divergence to the maintenance of a northern species complex: the case of the long‐toed salamander (Ambystoma macrodactylum Baird)

The relative importance of ecological vs. nonecological factors for the origin and maintenance of species is an open question in evolutionary biology. Young lineages – such as the distinct genetic groups that make up the ranges of many northern species – represent an opportunity to study the importance of ecological divergence during the early stages of diversification. Yet, few studies have examined the extent of niche divergence between lineages in previously glaciated regions and the role of ecology in maintaining the contact zones between them. In this study, we used tests of niche overlap in combination with ecological niche models to explore the extent of niche divergence between lineages of the long‐toed salamander (Ambystoma macrodactylum Baird) species complex and to determine whether contact zones correspond to (divergent) niche limits. We found limited evidence for niche divergence between the different long‐toed salamander lineages, substantial overlap in the predicted distribution of suitable climatic space for all lineages and range limits that are independent of niche limits. These results raise questions as to the importance of ecological divergence to the development of this widespread species complex and highlight the potential for non‐ecological factors to play a more important role in the maintenance of northern taxa.     

Quaternary phylogeography: the roots of hybrid zones

The older history of hybrid zones is explored through consideration of recent advances in climatology, paleontology and phylogeography in the Late Cenozoic, particularly the Quaternary Period with its major climatic cycles. The fossil record shows that these ice ages and their nested millennial oscillations caused substantial changes in species distributions and with genetic evidence allows deduction of refugia and colonization routes in arctic, temperate, desert and tropical regions. The age of divergence between hybridizing lineages varies from the Late Pleistocene to the Late Miocene, implying much range change and varying selection on sister lineages. Hybridizing lineages in the Tropical and Temperate regions range in age from young to old, but those studied in the Arctic are no more than a few ice ages old and their refugial roots are not clear. Mid to low latitude regions often show parapatric patchworks of lineages and multiple refugia stable through many climatic oscillations. Particular hybrid zones may have formed more than once; while some expansions were not the same, producing reticulation and introgression in previous glacial cycles. Hybrid-zone roots are complex and deep, and considerations of their complexity can reveal evolutionary pathways of species. They are indeed windows on evolution.     

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.     

Testing for intraspecific postzygotic isolation between cryptic lineages of Pseudacris crucifer

Phenotypically cryptic lineages appear common in nature, yet little is known about the mechanisms that initiate and/or maintain barriers to gene flow, or how secondary contact between them might influence evolutionary trajectories. The consequences of such contact between diverging lineages depend on hybrid fitness, highlighting the potential for postzygotic isolating barriers to play a role in the origins of biological species. Previous research shows that two cryptic, deeply diverged intraspecific mitochondrial lineages of a North American chorus frog, the spring peeper (Pseudacris crucifer), meet in secondary contact in Southwestern Ontario, Canada. Our study quantified hatching success, tadpole survival, size at metamorphosis, and development time for experimentally generated pure lineage and hybrid tadpoles. Results suggest that lineages differ in tadpole survival and that F₁ hybrids may have equal fitness and higher than average mass at metamorphosis compared with pure parental crosses. These findings imply hybrid early life viability may not be the pivotal reproductive isolation barrier helping to maintain lineage boundaries. However, we observed instances of tadpole gigantism, failure to metamorphose, and bent tails in some tadpoles from hybrid families. We also speculate and provide some evidence that apparent advantages or similarities of hybrids compared with pure lineage tadpoles may disappear when tadpoles are raised with competitors of different genetic makeup. This pilot study implies that ecological context and consideration of extrinsic factors may be a key to revealing mechanisms causing negative hybrid fitness during early life stages, a provocative avenue for future investigations on barriers to gene flow among these intraspecific lineages.     

Genetic diversity and the origins of parthenogenesis in the teiid lizard Aspidoscelis laredoensis

Unisexual vertebrates typically form through hybridization events between sexual species in which reproductive mode transitions occur in the hybrid offspring. This evolutionary history is thought to have important consequences for the ecology of unisexual lineages and their interactions with congeners in natural communities. However, these consequences have proven challenging to study owing to uncertainty about patterns of population genetic diversity in unisexual lineages. Of particular interest is resolving the contribution of historical hybridization events versus post formational mutation to patterns of genetic diversity in nature. Here we use restriction site associated DNA genotyping to evaluate genetic diversity and demographic history in Aspidoscelis laredoensis, a diploid unisexual lizard species from the vicinity of the Rio Grande River in southern Texas and northern Mexico. The sexual progenitor species from which one or more lineages are derived also occur in the Rio Grande Valley region, although patterns of distribution across individual sites are quite variable. Results from population genetic and phylogenetic analyses resolved the major axes of genetic variation in this species and highlight how these match predictions based on historical patterns of hybridization. We also found discordance between results of demographic modelling using different statistical approaches with the genomic data. We discuss these insights within the context of the ecological and evolutionary mechanisms that generate and maintain lineage diversity in unisexual species. As one of the most dynamic, intriguing, and geographically well investigated groups of whiptail lizards, these species hold substantial promise for future studies on the constraints of diversification in unisexual vertebrates.     

Hybridization,introgression, and linkage evolution

Genetic mapping methods provide a unique opportunity to study the interactions of differentiated genes and genomes in a hybrid genetic background. After a brief discussion of theoretical and analytical concerns, we review the application of these methods to a wide range of evolutionary issues. Map-based studies of experimental hybrids indicate that most postzygotic reproductive barriers in plants are polygenic and that the expression of extreme or novel traits in segregating hybrids (transgressive segregation) results from the complementary action of divergent parental alleles. However, genetic studies of hybrid vigor do not concur in their interpretations of the relative roles of dominance, overdominance, and epistasis. Map-based studies of natural hybrids are much rarer, but the few existing studies confirm the polygenic basis of postzygotic barriers and demonstrate the utility of genetic linkage for detecting cryptic introgression. In addition, studies of experimental and natural hybrid lineages provide compelling evidence that homoploid hybrid speciation has occurred in nature, and that it represents a rapid and repeatable mode of speciation. Data further indicate that this mode is facilitated by strong fertility selection and high chromosomal mutation rates. We recommend that future studies of hybrid genomes focus on natural hybrids, not only because of the paucity of data in this area, but also because of the availability of highly recombinant hybrid genotypes in hybrid zones. Of particular value will be studies of long-lived or difficult-to-propagate organisms, which previously have not been amenable to genetic study.