Signatures of hybridization and speciation in genomic patterns of ancestry*

Genomes sampled from hybrid zones between nascent species provide important clues into the speciation process. With advances in genome sequencing and single nucleotide polymorphism (SNP) genotyping, it is now feasible to measure variation in gene flow with high genomic resolution. This progress motivates the development of conceptual and analytical frameworks for hybrid zones that complement well‐established cline approaches. We extend the perspective that genomic distributions of ancestry are sensitive indicators of hybridization history. We use simulations to examine the behavior of the number of ancestry junctions—a simple summary of genomic patterns—in hybrid zones under increasingly realistic scenarios. Neutral simulations revealed that ancestry junction number is shaped by population structure, migration rate, and population size. Modeling multiple genetic architectures of hybrid dysfunction, with an emphasis on epistatic hybrid incompatibilities, showed that selection reduces junction number near loci that confer reproductive barriers. The magnitude of this signature was affected by the form of selection, dominance, and genomic location (autosome vs. sex chromosome) of incompatible loci. Our results suggest that researchers can identify loci involved in reproductive isolation by scanning hybrid genomes for local reductions in junction number. We outline necessary directions for future theory and method development to realize this goal.     

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.     

Geographic variation in the structure of oak hybrid zones provides insights into the dynamics of speciation

Studying geographic variation in the rate of hybridization between closely related species could provide a useful window on the evolution of reproductive isolation. Reinforcement theory predicts greater prezygotic isolation in areas of prolonged contact between recently diverged species than in areas of recent contact, which implies that old contact zones would be dominated by parental phenotypes with few hybrids (bimodal hybrid zones), whereas recent contact zones would be characterized by hybrid swarms (unimodal hybrid zones). Here, we investigate how the hybrid zones of two closely related Chinese oaks, Quercus mongolica and Q. liaotungensis, are structured geographically using both nuclear and chloroplast markers. We found that populations of Q. liaotungensis located around the Changbai Mountains in Northeast China, an inferred glacial refugium, were introgressed by genes from Q. mongolica, suggesting historical contact between the two species in this region. However, these introgressed populations form sharp bimodal hybrid zones with Q. mongolica. In contrast, populations of Q. liaotungensis located in North China, which show no sign of ancient introgression with Q. mongolica, form unimodal hybrid zones with Q. mongolica. These results are consistent with the hypothesis that selection against hybrids has had sufficient time to reinforce the reproductive barriers between Q. liaotungensis and Q. mongolica in Northeast China but not in North China.     

Hybridization promotes speciation in Coenonympha butterflies

Hybridization has become a central element in theories of animal evolution during the last decade. New methods in population genomics and statistical model testing now allow the disentangling of the complexity that hybridization brings into key evolutionary processes such as local adaptation, colonization of new environments, species diversification and extinction. We evaluated the consequences of hybridization in a complex of three alpine butterflies in the genus Coenonympha, by combining morphological, genetic and ecological analyses. A series of approximate Bayesian computation procedures based on a large SNP data set strongly suggest that the Darwin's Heath (Coenonympha darwiniana) originated through hybridization between the Pearly Heath (Coenonympha arcania) and the Alpine Heath (Coenonympha gardetta) with different parental contributions. As a result of hybridization, the Darwin's Heath presents an intermediate morphology between the parental species, while its climatic niche seems more similar to the Alpine Heath. Our results also reveal a substantial genetic and morphologic differentiation between the two geographically disjoint Darwin's Heath lineages leading us to propose the splitting of this taxon into two different species.     

Sperm-related phenotypes implicated in both maintenance and breakdown of a natural species barrier in the house mouse

The house mouse hybrid zone (HMHZ) is a species barrier thought to be maintained by a balance between dispersal and natural selection against hybrids. While the HMHZ is characterized by frequency discontinuities for some sex chromosome markers, there is an unexpected large-scale regional introgression of a Y chromosome across the barrier, in defiance of Haldane''s rule. Recent work suggests that a major force maintaining the species barrier acts through sperm traits. Here, we test whether the Y chromosome penetration of the species barrier acts through sperm traits by assessing sperm characteristics of wild-caught males directly in a field laboratory set up in a Y introgression region of the HMHZ, later calculating the hybrid index of each male using 1401 diagnostic single nucleotide polymorphisms (SNPs). We found that both sperm count (SC) and sperm velocity were significantly reduced across the natural spectrum of hybrids. However, SC was more than rescued in the presence of the invading Y. Our results imply an asymmetric advantage for Y chromosome introgression consistent with the observed large-scale introgression. We suggest that selection on sperm-related traits probably explains a large component of patterns observed in the natural hybrid zone, including the Y chromosome penetration.     

Hybridization rapidly reduces fitness of a native trout in the wild

Human-mediated hybridization is a leading cause of biodiversity loss worldwide. How hybridization affects fitness and what level of hybridization is permissible pose difficult conservation questions with little empirical information to guide policy and management decisions. This is particularly true for salmonids, where widespread introgression among non-native and native taxa has often created hybrid swarms over extensive geographical areas resulting in genomic extinction. Here, we used parentage analysis with multilocus microsatellite markers to measure how varying levels of genetic introgression with non-native rainbow trout (Oncorhynchus mykiss) affect reproductive success (number of offspring per adult) of native westslope cutthroat trout (Oncorhynchus clarkii lewisi) in the wild. Small amounts of hybridization markedly reduced fitness of male and female trout, with reproductive success sharply declining by approximately 50 per cent, with only 20 per cent admixture. Despite apparent fitness costs, our data suggest that hybridization may spread due to relatively high reproductive success of first-generation hybrids and high reproductive success of a few males with high levels of admixture. This outbreeding depression suggests that even low levels of admixture may have negative effects on fitness in the wild and that policies protecting hybridized populations may need reconsideration.     

Sculpin hybrid zones: natural laboratories for the early stages of speciation

Firmly rooted as we are in the genomic era, it can seem incredible that as recently as 1974, Lewontin declared, 'we know virtually nothing about the genetic changes that occur in species formation'. To the contrary, we now know the genetic architecture of phenotypic differences and reproductive isolation between species for many diverse groups of plants, animals, and fungi. In recent years, detailed genetic analyses have produced a small but growing list of genes that cause reproductive isolation, several of which appear to have diverged by natural selection. Yet, a full accounting of the speciation process requires that we understand the reproductive and ecological properties of natural populations as they begin to diverge genetically, as well as the dynamics of newly evolved barriers to gene flow. One promising approach to this problem is the study of natural hybrid zones, where gene exchange between divergent populations can produce recombinant genotypes in situ . In such individuals, genomic variation might be shaped by introgression at universally adaptive or neutral loci, even as regions associated with local adaptation or reproductive isolation remain divergent. In Nolte et   al . (2009) , the authors take advantage of two independent, recently formed hybrid zones between sculpin species to investigate genome-wide patterns of reproductive isolation. Using a recently developed genomic clines method, the authors identify marker loci that are associated with isolation, and those that show evidence for adaptive introgression. Remarkably, Nolte et   al . (2009) find little similarity between the two hybrid zones in patterns of introgression, a fact that might reflect genetic variation within species or heterogeneous natural selection. In either case, their study system has the potential to provide insight into the early stages of speciation.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process of speciation.     

Genomics of isolation in hybrids

Hybrid zones are common in nature and can offer critical insights into the dynamics and components of reproductive isolation. Hybrids between diverged lineages are particularly informative about the genetic architecture of reproductive isolation, because introgression in an admixed population is a direct measure of isolation. In this paper, we combine simulations and a new statistical model to determine the extent to which different genetic architectures of isolation leave different signatures on genome-level patterns of introgression. We found that reproductive isolation caused by one or several loci of large effect caused greater heterogeneity in patterns of introgression than architectures involving many loci with small fitness effects, particularly when isolating factors were closely linked. The same conditions that led to heterogeneous introgression often resulted in a reasonable correspondence between outlier loci and the genetic loci that contributed to isolation. However, demographic conditions affected both of these results, highlighting potential limitations to the study of the speciation genomics. Further progress in understanding the genomics of speciation will require large-scale empirical studies of introgression in hybrid zones and model-based analyses, as well as more comprehensive modelling of the expected levels of isolation with different demographies and genetic architectures of isolation.     

Conservation Implications of Invasion by Plant Hybridization

The increasing number of invasive exotic plant species in many regions and the continuing alteration of natural ecosystems by humans promote hybridization between previously allopatric species; among both native as well as between native and introduced species. We review the ecological factors and mechanisms that promote such hybridization events and their negative consequences on biological diversity. Plant invasions through hybridization may occur in four different ways: hybridization between native species, hybridization between an exotic species and a native congener, hybridization between two exotics and by the introduction and subsequent spread of hybrids. The main harmful genetic effect of such hybrids on native species is the loss of both genetic diversity and of locally adapted populations, such as rare and threatened species. The spread of aggressive hybrid taxa can reduce the growth of, or replace, native species. The main factor promoting the formation of hybrids is species dispersal promoted by humans. However, the success and spread of hybrids is increased by disturbance and fragmentation of habitats, thus overcoming natural crossing barriers, and range expansions due to human activity. There are differences in flowering, pollination and seed dispersal patterns between parental species and hybrids. Hybrid resistance to pathogens and herbivores may also enhance the success of hybrids. To predict the mechanisms and consequences of invasions mediated by hybridization, extensive data on hybrid ecology and biology are needed, as well as carefully designed field experiments focused on the comparative ecology of parental populations and hybrids.     

Diverse paths to hybrid incompatibility in Arabidopsis

One of the most essential questions of biology is to understand how different species have evolved. Hybrid incompatibility, a phenomenon in which hybrids show reduced fitness in comparison with their parents, can result in reproductive isolation and speciation. Therefore, studying hybrid incompatibility provides an entry point in understanding speciation. Hybrid incompatibilities are known throughout taxa, and the underlying mechanisms have mystified scientists since the theory of evolution by means of natural selection was introduced. In plants, it is only in recent years that the high‐throughput genetic and molecular tools have become available for the Arabidopsis genus, thus helping to shed light on the different genes and molecular and evolutionary mechanisms that underlie hybrid incompatibilities. In this review, we highlight the current knowledge of diverse mechanisms that are known to contribute to hybrid incompatibility.     

Hybridization in headwater regions,and the role of rivers as drivers of speciation in Amazonian birds

Many understory birds and other groups form genetically differentiated subspecies or closely related species on opposite sides of major rivers of Amazonia, but are proposed to come into geographic contact in headwater regions where narrower river widths may present less of a dispersal barrier. Whether such forms hybridize in headwater regions is generally unknown, but has important implications to our understanding of the role of rivers as drivers of speciation. We used a dataset of several thousand single nucleotide polymorphisms to show that seven taxon pairs that differentiate across a major Amazonian river come into geographic contact and hybridize in headwater regions. All taxon pairs possessed hybrids with low numbers of loci in which alleles were inherited from both parental species, suggesting they are backcrossed with parentals, and indicating gene flow between parental populations. Ongoing gene flow challenges rivers as the sole cause of in situ speciation, but is compatible with the view that the wide river courses in the heart of Amazonia may have driven interfluvial divergence during episodes of wet forest retraction away from headwater regions. Taxa as old as 4 Ma in our Amazonian dataset continue to hybridize at contact zones, suggesting reproductive isolation evolves at a slow pace.     

Quantified reproductive isolation in Heliconius butterflies: Implications for introgression and hybrid speciation

Heliconius butterflies have become a model for the study of speciation with gene flow. For adaptive introgression to take place, there must be incomplete barriers to gene exchange that allow interspecific hybridization and multiple generations of backcrossing. The recent publication of estimates of individual components of reproductive isolation between several species of butterflies in the Heliconius melpomene–H. cydno clade allowed us to calculate total reproductive isolation estimates for these species. According to these estimates, the butterflies are not as promiscuous as has been implied. Differences between species are maintained by intrinsic mechanisms, while reproductive isolation of geographical races within species is mainly due to allopatry. We discuss the implications of this strong isolation for basic aspects of the hybrid speciation with introgression hypothesis.     

Hybridization in a warmer world

Climate change is profoundly affecting the evolutionary trajectory of individual species and ecological communities, in part through the creation of novel species assemblages. How climate change will influence competitive interactions has been an active area of research. Far less attention, however, has been given to altered reproductive interactions. Yet, reproductive interactions between formerly isolated species are inevitable as populations shift geographically and temporally as a result of climate change, potentially resulting in introgression, speciation, or even extinction. The susceptibility of hybridization rates to anthropogenic disturbance was first recognized in the 1930s. To date, work on anthropogenically mediated hybridization has focused primarily on either physical habitat disturbance or species invasion. Here, I review recent literature on hybridization to identify how ecological responses to climate change will increase the likelihood of hybridization via the dissolution of species barriers maintained by habitat, time, or behavior. Using this literature, I identify several cases where novel hybrid zones have recently formed, likely as a result of changing climate. Future research should focus on identifying areas and taxonomic groups where reproductive species interactions are most likely to be influenced by climate change. Furthermore, a better understanding of the evolutionary consequences of climate‐mediated secondary contact is urgently needed. Paradoxically, hybridization is both a major conservation concern and an important source of novel genetic and phenotypic variation. Hybridization may therefore both contribute to increasing rates of extinction and stimulate the creation of novel phenotypes that will speed adaptation to novel climates. Predicting which result will occur following secondary contact will be an important contribution to conservation for many species.     

Increasing evidence of the role of gene flow in animal evolution: hybrid speciation in the yellow-rumped warbler complex

In this issue of Molecular Ecology, Brelsford et al. (2011) present strong evidence for a case of hybrid speciation within the yellow-rumped warbler complex. Although homoploid hybrid speciation has now been documented in many animals (Mallet 2007), it seems rare in tetrapods (Mavárez & Linares 2008) and it has barely even been mentioned in birds (Price 2008). Brelsford and colleagues thus present the first detailed molecular evidence suggesting that hybrid speciation can occur in birds. Brelsford et al. (2011) posit that Audubon's warbler (Dendroica auduboni) constitutes a hybrid species originating from the admixture of two distinct parental lineages, represented today by myrtle warbler (D. coronata) and black-fronted warbler (D. nigrifrons). The authors present three major lines of molecular evidence suggesting that this is not simply a case of a hybrid swarm or limited introgression.     

Hybridization and polyploidy as drivers of continuing evolution and speciation in Sorbus

Interspecific hybridization and polyploidy are pivotal processes in plant evolution and speciation. The fate of new hybrid and polyploid taxa is determined by their ability to reproduce either sexually or asexually. Hybrids and allopolyploids with odd chromosome numbers are frequently sterile but some establish themselves through asexual reproduction (vegetative or apomixis). This allows novel genotypes to become established by isolating them from gene flow and leads to complex patterns of variation. The genus Sorbus is a good example of taxonomic complexity arising from the combined effects of hybridization, polyploidy and apomixis. The Avon Gorge in South‐west Britain contains the greatest diversity of Sorbus in Europe, with three endemic species and four putative endemic novel hybrids among its 15 native Sorbus taxa. We used a combination of nuclear microsatellite and chloroplast DNA markers to investigate the evolutionary relationships among these Sorbus taxa within the Avon Gorge. We confirm the genetic identity of putative novel taxa and show that hybridization involving sexual diploid species, primarily S. aria and S. torminalis and polyploid facultative apomictic species from subgenus Aria, has been responsible for generating this biodiversity. Importantly our data show that this creative evolutionary process is ongoing within the Avon Gorge. Conservation strategies for the rare endemic Sorbus taxa should therefore consider all Sorbus taxa within the Gorge and must strive to preserve this evolutionary process rather than simply the individual rare taxa that it produces.     

Phylogeography of speciation: allopatric divergence and secondary contact between outcrossing and selfing Clarkia

The origins of hybrid zones between parapatric taxa have been of particular interest for understanding the evolution of reproductive isolation and the geographic context of species divergence. One challenge has been to distinguish between allopatric divergence (followed by secondary contact) versus primary intergradation (parapatric speciation) as alternative divergence histories. Here, we use complementary phylogeographic and population genetic analyses to investigate the recent divergence of two subspecies of Clarkia xantiana and the formation of a hybrid zone within the narrow region of sympatry. We tested alternative phylogeographic models of divergence using approximate Bayesian computation (ABC) and found strong support for a secondary contact model and little support for a model allowing for gene flow throughout the divergence process (i.e. primary intergradation). Two independent methods for inferring the ancestral geography of each subspecies, one based on probabilistic character state reconstructions and the other on palaeo-distribution modelling, also support a model of divergence in allopatry and range expansion leading to secondary contact. The membership of individuals to genetic clusters suggests geographic substructure within each taxon where allopatric and sympatric samples are primarily found in separate clusters. We also observed coincidence and concordance of genetic clines across three types of molecular markers, which suggests that there is a strong barrier to gene flow. Taken together, our results provide evidence for allopatric divergence followed by range expansion leading to secondary contact. The location of refugial populations and the directionality of range expansion are consistent with expectations based on climate change since the last glacial maximum. Our approach also illustrates the utility of combining phylogeographic hypothesis testing with species distribution modelling and fine-scale population genetic analyses for inferring the geography of the divergence process.     

Ecological speciation

Ecological processes are central to the formation of new species when barriers to gene flow (reproductive isolation) evolve between populations as a result of ecologically‐based divergent selection. Although laboratory and field studies provide evidence that ‘ecological speciation’ can occur, our understanding of the details of the process is incomplete. Here we review ecological speciation by considering its constituent components: an ecological source of divergent selection, a form of reproductive isolation, and a genetic mechanism linking the two. Sources of divergent selection include differences in environment or niche, certain forms of sexual selection, and the ecological interaction of populations. We explore the evidence for the contribution of each to ecological speciation. Forms of reproductive isolation are diverse and we discuss the likelihood that each may be involved in ecological speciation. Divergent selection on genes affecting ecological traits can be transmitted directly (via pleiotropy) or indirectly (via linkage disequilibrium) to genes causing reproductive isolation and we explore the consequences of both. Along with these components, we also discuss the geography and the genetic basis of ecological speciation. Throughout, we provide examples from nature, critically evaluate their quality, and highlight areas where more work is required.     

Cytonuclear incompatibility contributes to the early stages of speciation

Genetic incompatibility is a hallmark of speciation. Cytonuclear incompatibilities are proposed to be among the first genetic barriers to arise during speciation. Accordingly, reproductive isolation (RI) within species should be heavily influenced by interactions between the organelle and nuclear genomes. However, there are few clear examples of cytonuclear incompatibility within a species. Here, we show substantial postzygotic RI in first‐generation hybrids between differentiated populations of an herbaceous plant (up to 92% reduction in fitness). RI was primarily due to germination and survival, with moderate RI for pollen viability. RI for survival was asymmetric and caused by cytonuclear incompatibility, with the strength of incompatibility linearly related to chloroplast genetic distance. This cytonuclear incompatibility may be the result of a rapidly evolving plastid genome. Substantial asymmetric RI was also found for germination, but was not associated with cytonuclear incompatibility, indicating endosperm or maternal‐zygote incompatibilities. These results demonstrate that cytonuclear incompatibility contributes to RI within species, suggesting that initial rates of speciation could be influenced by rates of organelle evolution. However, other genetic incompatibilities are equally important, indicating that even at early stages, speciation can be a complex process involving multiple genes and incompatibilities.     

Interspecific crossing and genetic mapping reveal intrinsic genomic incompatibility between two Senecio species that form a hybrid zone on Mount Etna,Sicily

Studies of hybridizing species can reveal much about the genetic basis and maintenance of species divergence in the face of gene flow. Here we report a genetic segregation and linkage analysis conducted on F₂ progeny of a reciprocal cross between Senecio aethnensis and S. chrysanthemifolius that form a hybrid zone on Mount Etna, Sicily, aimed at determining the genetic basis of intrinsic hybrid barriers between them. Significant transmission ratio distortion (TRD) was detected at 34 (∼27%) of 127 marker loci located in nine distinct clusters across seven of the ten linkage groups detected, indicating genomic incompatibility between the species. TRD at these loci could not be attributed entirely to post-zygotic selective loss of F₂ individuals that failed to germinate or flower (16.7%). At four loci tests indicated that pre-zygotic events, such as meiotic drive in F₁ parents or gametophytic selection, contributed to TRD. Additional tests revealed that cytonuclear incompatibility contributed to TRD at five loci, Bateson–Dobzhansky–Muller (BDM) incompatibilities involving epistatic interactions between loci contributed to TRD at four loci, and underdominance (heterozygote disadvantage) was a possible cause of TRD at one locus. Major chromosomal rearrangements were probably not a cause of interspecific incompatibility at the scale that could be examined with current map marker density. Intrinsic genomic incompatibility between S. aethnensis and S. chrysanthemifolius revealed by TRD across multiple genomic regions in early-generation hybrids is likely to impact the genetic structure of the natural hybrid zone on Mount Etna by limiting introgression and promoting divergence across the genome.