Postzygotic isolation over multiple generations of hybrid descendents in a natural hybrid zone: how well do single-generation estimates reflect reproductive isolation?

A history of studies of reproductive isolation between Drosophila pseudoobscura and D. persimilis

Drosophila pseudoobscura and D. persimilis are a sister species pair that have been used as a model for studies of reproductive isolation and speciation for almost 100 years owing to their close evolutionary history, well characterized genetic differences, and overlapping geographic distribution. There are extensive analyses of both pre- and post-zygotic isolation, including studies of courtship divergence, conspecific sperm precedence (CSP) and how reinforcement by natural selection may or may not act to strengthen isolation in sympatry. Post-zygotic analyses explore the underlying mechanics of reproductive isolation; how inversions may give rise to initial speciation events and misexpression of key genes typically found within inversion regions render hybrid offspring unfit or inviable. We aim here to present a history of studies of reproductive isolation between this species pair, looking at how the field has developed over the last century and identifying the open questions and gaps within the literature.

A new bunting species in South China revealed by an integrative taxonomic investigation of the Emberiza godlewskii complex (Aves, Emberidae)

Species boundaries of birds across the Euro-Asian continent, especially the widely distributed passerines with multiple subspecies, are not completely resolved. The Rock Bunting complex, Emberiza cia and E. godlewskii, is subject to ongoing debate due to geographic variation in its morphology, and deep genetic splits. In this study, the phylogenetic relationships of the E. cia/godlewskii complex were evaluated based on multilocus coalescence methods combined with analyses of plumage color and morphological variation. Although the northern and southern populations of E. godlewski formed reciprocally monophyletic groups, the northern E. godlewskii clade is the sister group to E. cia, rendering E. godlewskii as currently defined paraphyletic. The significant differences in morphometry and plumage color reflectance are complementary and support the phenotypic validity of northern and southern E. godlewskii. Thus, we recommend the following taxonomic revisions: the northern and southern E. godlewskii subspecies should be recognized as different full species, such as E. godlewskii and E. yunnanensis, respectively. This resolves the issue that E. godlewskii, as currently defined, is paraphyletic. Our research provides direct evidence of cryptic species in old world buntings and highlights the underestimation of species diversity in birds in East Asia.

Development of InDel markers for interspecific hybridization between hill pigeons and feral pigeons based on whole-genome re-sequencing

Interspecific hybridization occurs among birds, and closely related sister taxa tend to hybridize at a high rate. Genomic hybridization markers are useful for understanding the patterns and processes of hybridization and for conserving endangered species in captivity and the wild. In this study, we developed genomic hybridization markers for the F1 progeny of the sister taxa feral pigeons (Columba livia var. domestica) and endangered hill pigeons (Columba rupestris) (family Columbidae). Using whole-genome re-sequencing data, we performed genome-wide analysis for insertion/deletion (InDel) polymorphisms and validated using primers. We conducted polymerase chain reaction (PCR) and agarose gel electrophoresis to identify species-specific InDels. We produced eight F1 hybrids of hill and feral pigeons, and their samples were tested by re-performing analyses and sequencing using 11 species-specific InDel polymorphisms. Eight InDel markers simultaneously amplified two DNA fragments from all F1 hybrids, and there was no abnormality in the sequencing results. The application of genomic tools to detect hybrids can play a crucial role in the assessment of hybridization frequency in the wild. Moreover, systematic captive propagation efforts with hybrids can help control the population decline of hill pigeons.

Variable Signatures of Selection Despite Conserved Recombination Landscapes Early in Speciation

Recently diverged taxa often exhibit heterogeneous landscapes of genomic differentiation, characterized by regions of elevated differentiation on an otherwise homogeneous background. While divergence peaks are generally interpreted as regions responsible for reproductive isolation, they can also arise due to background selection, selective sweeps unrelated to speciation, and variation in recombination and mutation rates. To investigate the association between patterns of recombination and landscapes of genomic differentiation during the early stages of speciation, we generated fine-scale recombination maps for six southern capuchino seedeaters (Sporophila) and two subspecies of White Wagtail (Motacilla alba), two recent avian radiations in which divergent selection on pigmentation genes has likely generated peaks of differentiation. We compared these recombination maps to those of Collared (Ficedula albicollis) and Pied Flycatchers (Ficedula hypoleuca), non-sister taxa characterized by moderate genomic divergence and a heterogenous landscape of genomic differentiation shaped in part by background selection. Although recombination landscapes were conserved within all three systems, we documented a weaker negative correlation between recombination rate and genomic differentiation in the recent radiations. All divergence peaks between capuchinos, wagtails, and flycatchers were located in regions with lower-than-average recombination rates, and most divergence peaks in capuchinos and flycatchers fell in regions of exceptionally reduced recombination. Thus, co-adapted allelic combinations in these regions may have been protected early in divergence, facilitating rapid diversification. Despite largely conserved recombination landscapes, divergence peaks are specific to each focal comparison in capuchinos, suggesting that regions of elevated differentiation have not been generated by variation in recombination rate alone.

Extent and direction of introgressive hybridization of mule and white-tailed deer in western Canada

Hybridization of mule deer (Odocoileus hemionus) and white-tailed deer (O. virginianus) appears to be a semi-regular occurrence in western North America. Previous studies confirmed the presence of hybrids in a variety of sympatric habitats, but their developing molecular resources limited identification to the earliest, most admixed generations. For this reason, estimates of hybrid production in wild populations often rely on anecdotal reports. As well, white-tailed deer populations' continued encroachment into historically mule deer-occupied habitats due to changes in land use, habitat homogenization, and a warming climate may increase opportunities for interspecific encounters. We sought to quantify the prevalence and extent of hybrid deer in the prairies of western Canada using a SNP assay with enhanced discriminating power. By updating the available molecular resources, we sought to identify and characterize previously cryptic introgression. We also investigated the influence of various parameters on hybridity by way of logistic regression. We observed overall hybridization rates of ~1.0%, slightly lower than that reported by previous studies, and found white-tailed-like hybrids to be more common than their mule deer-like counterparts. Here, we build upon past studies of hybridization in North American deer by increasing hybrid detection power, expanding sample sizes, demonstrating a new molecular resource applicable to future research and observing asymmetrical directionality of introgression.

Stable species boundaries despite ten million years of hybridization in tropical eels

Genomic evidence is increasingly underpinning that hybridization between taxa is commonplace, challenging our views on the mechanisms that maintain their boundaries. Here, we focus on seven catadromous eel species (genus Anguilla) and use genome-wide sequence data from more than 450 individuals sampled across the tropical Indo-Pacific, morphological information, and three newly assembled draft genomes to compare contemporary patterns of hybridization with signatures of past introgression across a time-calibrated phylogeny. We show that the seven species have remained distinct for up to 10 million years and find that the current frequencies of hybridization across species pairs contrast with genomic signatures of past introgression. Based on near-complete asymmetry in the directionality of hybridization and decreasing frequencies of later-generation hybrids, we suggest cytonuclear incompatibilities, hybrid breakdown, and purifying selection as mechanisms that can support species cohesion even when hybridization has been pervasive throughout the evolutionary history of clades.

Body condition (Wr ) and reproductive potential of bighead and silver carp hybrids: Postzygotic selection in the Mississippi River Basin

Invasive bighead (Hypophthalmichthys nobilis) and silver carp (H. molitrix) are reproductively isolated in their native range, but form a bimodal, multigenerational hybrid swarm within the Mississippi River Basin (MRB). Despite observed F1 hybrid superiority in experimental settings, effects of postzygotic selection on bighead and silver carp hybrids have not been tested in a natural system. Individual parent and hybrid genotypes were resolved at 57 species-specific loci and used to evaluate postzygotic selection for body condition (Wr ) and female reproductive potential (presence of spawning stage gonads and gonadosomatic index [GSI]) in the MRB during 2009-2011. Body condition in the Marseilles Reach, Illinois River declined with a decrease in species-specific allele frequency from 1.0 to 0.4 for each species and early generation hybrids (F1, F2, and first-generation backcross) had lower mean Wr than late generation hybrids (2nd+ generation backcrosses) and parentals. Proportions of stage IV and stage V (spawning stage) female gonads differed between bighead and silver carp, but not among parentals and their early and late generation hybrids within the MRB. Mean GSI values did not differ between parentals and hybrids. Because reproductive potential did not differ between hybrids and parentals, our results suggest that early generation hybrids occur in low frequency either as a factor of poor condition (Wr ) and postreproductive survival, infrequent reproductive encounters by parental bighead and silver carp, or selection pressures acting on juvenile or immature life stages. Our results suggest that a combination of genetic and environmental factors may contribute to the postzygotic success of bighead and silver carp hybrids in the Mississippi River Basin.

Development of a Novel Mule Deer Genomic Assembly and Species-Diagnostic SNP Panel for Assessing Introgression in Mule Deer, White-Tailed Deer, and Their Interspecific Hybrids

Mule deer (Odocoileus hemionus) are endemic to a wide variety of habitats in western North America, many of which are shared in sympatry with their closely related sister-species white-tailed deer (Odocoileus virginianus), whom they hybridize with in wild populations. Although mule deer meet many ideal conditions for a molecular ecological research species, such as high abundance, ecological importance, and broad dispersal and gene flow, conservation genetic studies have been limited by a relative lack of existing genomic resources and inherent difficulties caused by introgression with white-tailed deer. Many molecular tools currently available for the study of cervids were designed using reference assemblies of divergent model species, specifically cattle (Bos taurus). Bovidae and Cervidae diverged approximately 28 million years ago, therefore, we sought to ameliorate the available resources by contributing the first mule deer whole genome sequence draft assembly with an average genome-wide read depth of 25X, using the white-tailed genome assembly (Ovir.te_1.0) as a reference. Comparing the two assemblies, we identified ∼33 million single nucleotide polymorphisms (SNPs) and insertion/deletion variants. We then verified fixed SNP differences between the two species and developed a 40-loci SNP assay capable of identifying pure mule deer, white-tailed deer, and interspecific hybrids. Assignment capacity of the panel, which was tested on simulated datasets, is reliable up to and including the third backcross hybrid generation. Identification of post-F1 hybrids will be necessary for hybrid zone population studies going forward, and the new mule deer assembly will be a valuable resource for genetic and comparative genomics studies.

Morphologically cryptic Amazonian bird species pairs exhibit strong postzygotic reproductive isolation

We possess limited understanding of how speciation unfolds in the most species-rich region of the planet-the Amazon basin. Hybrid zones provide valuable information on the evolution of reproductive isolation, but few studies of Amazonian vertebrate hybrid zones have rigorously examined the genome-wide underpinnings of reproductive isolation. We used genome-wide genetic datasets to show that two deeply diverged, but morphologically cryptic sister species of forest understorey birds show little evidence for prezygotic reproductive isolation, but substantial postzygotic isolation. Patterns of heterozygosity and hybrid index revealed that hybrid classes with heavily recombined genomes are rare and closely match simulations with high levels of selection against hybrids. Genomic and geographical clines exhibit a remarkable similarity across loci in cline centres, and have exceptionally narrow cline widths, suggesting that postzygotic isolation is driven by genetic incompatibilities at many loci, rather than a few loci of strong effect. We propose Amazonian understorey forest birds speciate slowly via gradual accumulation of postzygotic genetic incompatibilities, with prezygotic barriers playing a less important role. Our results suggest old, cryptic Amazonian taxa classified as subspecies could have substantial postzygotic isolation deserving species recognition and that species richness is likely to be substantially underestimated in Amazonia.

Evaluating genomic signatures of "the large X-effect" during complex speciation

The ubiquity of the "two rules of speciation"-Haldane's rule and the large X-effect-implies a general, special role for sex chromosomes in the evolution of intrinsic postzygotic reproductive isolation. The recent proliferation of genome-scale analyses has revealed two further general observations: (a) complex speciation involving some form of gene flow is not uncommon, and (b) sex chromosomes in male- and in female-heterogametic taxa tend to show elevated differentiation relative to autosomes. Together, these observations are consistent with speciation histories in which population genetic differentiation at autosomal loci is reduced by gene flow while natural selection against hybrid incompatibilities renders sex chromosomes relatively refractory to gene flow. Here, I summarize multilocus population genetic and population genomic evidence for greater differentiation on the X (or Z) vs. the autosomes and consider the possible causes. I review common population genetic circumstances involving no selection and/or no interspecific gene flow that are nevertheless expected to elevate differentiation on sex chromosomes relative to autosomes. I then review theory for why large X-effects exist for hybrid incompatibilities and, more generally, for loci mediating local adaptation. The observed levels of sex chromosome vs. autosomal differentiation, in many cases, appear consistent with simple explanations requiring neither large X-effects nor gene flow. Discerning signatures of large X-effects during complex speciation will therefore require analyses that go beyond chromosome-scale summaries of population genetic differentiation, explicitly test for differential introgression, and/or integrate experimental genetic data.

Hybridization and postzygotic isolation promote reinforcement of male mating preferences in a diverse group of fishes with traditional sex roles

Behavioral isolation is thought to arise early in speciation due to differential sexual and/or natural selection favoring different preferences and traits in different lineages. Instead, behavioral isolation can arise due to reinforcement favoring traits and preferences that prevent maladaptive hybridization. In darters, female preference for male coloration has been hypothesized to drive speciation, because behavioral isolation evolves before F1 inviability. However, as with many long-lived organisms, the fitness of second-generation hybrids has not been assessed because raising animals to adulthood in the laboratory is challenging. Of late, reinforcement of male preferences has been implicated in darters because male preference for conspecific females is high in sympatry but absent in allopatry in multiple species pairs. The hypothesis that reinforcement accounts for behavioral isolation in sympatry assumes that hybridization and postzygotic isolation are present. Here, we used genomic and morphological data to demonstrate that hybridization is ongoing between orangethroat and rainbow darters and used hybrids collected from nature to measure postzygotic barriers across two hybrid generations. We observed sex ratio distortion in adult F1s and a dramatic reduction in backcross survival. Our findings indicate that selection to avoid hybridization promotes the evolution of male-driven behavioral isolation via reinforcement in this system.

Geographic variation in hybridization across a reinforcement contact zone of chorus frogs (Pseudacris)

Reinforcement contact zones, which are secondary contact zones where species are diverging in reproductive behaviors due to selection against hybridization, represent natural laboratories for studying speciation‐in‐action. Here, we examined replicate localities across the entire reinforcement contact zone between North American chorus frogs Pseudacris feriarum and P. nigrita to investigate geographic variation in hybridization frequencies and to assess whether reinforcement may have contributed to increased genetic divergence within species. Previous work indicated these species have undergone reproductive character displacement (RCD) in male acoustic signals and female preferences due to reinforcement. We also examined acoustic signal variation across the contact zone to assess whether signal characteristics reliably predict hybrid index and to elucidate whether the degree of RCD predicts hybridization rate. Using microsatellites, mitochondrial sequences, and acoustic signal information from >1,000 individuals across >50 localities and ten sympatric focal regions, we demonstrate: (1) hybridization occurs and (2) varies substantially across the geographic range of the contact zone, (3) hybridization is asymmetric and in the direction predicted from observed patterns of asymmetric RCD, (4) in one species, genetic distance is higher between conspecific localities where one or both have been reinforced than between nonreinforced localities, after controlling for geographic distance, (5) acoustic signal characters strongly predict hybrid index, and (6) the degree of RCD does not strongly predict admixture levels. By showing that hybridization occurs in all sympatric localities, this study provides the fifth and final line of evidence that reproductive character displacement is due to reinforcement in the chorus frog contact zone. Furthermore, this work suggests that the dual action of cascade reinforcement and partial geographic isolation is promoting genetic diversification within one of the reinforced species.

Competition-driven build-up of habitat isolation and selection favoring modified dispersal patterns in a young avian hybrid zone

Competition-driven evolution of habitat isolation is an important mechanism of ecological speciation but empirical support for this process is often indirect. We examined how an on-going displacement of pied flycatchers from their preferred breeding habitat by collared flycatchers in a young secondary contact zone is associated with (a) access to an important food resource (caterpillar larvae), (b) immigration of pied flycatchers in relation to habitat quality, and (c) the risk of hybridization in relation to habitat quality. Over the past 12 years, the estimated access to caterpillar larvae biomass in the habitat surrounding the nests of pied flycatchers has decreased by a fifth due to shifted establishment possibilities, especially for immigrants. However, breeding in the high quality habitat has become associated with such a high risk of hybridization for pied flycatchers that overall selection currently favors pied flycatchers that were forced to immigrate into the poorer habitats (despite lower access to preferred food items). Our results show that competition-driven habitat segregation can lead to fast habitat isolation, which per se caused an opportunity for selection to act in favor of future "voluntarily" altered immigration patterns and possibly strengthened habitat isolation through reinforcement.

Resolving Evolutionary Relationships in Closely Related Species with Whole-Genome Sequencing Data

Using genetic data to resolve the evolutionary relationships of species is of major interest in evolutionary and systematic biology. However, reconstructing the sequence of speciation events, the so-called species tree, in closely related and potentially hybridizing species is very challenging. Processes such as incomplete lineage sorting and interspecific gene flow result in local gene genealogies that differ in their topology from the species tree, and analyses of few loci with a single sequence per species are likely to produce conflicting or even misleading results. To study these phenomena on a full phylogenomic scale, we use whole-genome sequence data from 200 individuals of four black-and-white flycatcher species with so far unresolved phylogenetic relationships to infer gene tree topologies and visualize genome-wide patterns of gene tree incongruence. Using phylogenetic analysis in nonoverlapping 10-kb windows, we show that gene tree topologies are extremely diverse and change on a very small physical scale. Moreover, we find strong evidence for gene flow among flycatcher species, with distinct patterns of reduced introgression on the Z chromosome. To resolve species relationships on the background of widespread gene tree incongruence, we used four complementary coalescent-based methods for species tree reconstruction, including complex modeling approaches that incorporate post-divergence gene flow among species. This allowed us to infer the most likely species tree with high confidence. Based on this finding, we show that regions of reduced effective population size, which have been suggested as particularly useful for species tree inference, can produce positively misleading species tree topologies. Our findings disclose the pitfalls of using loci potentially under selection as phylogenetic markers and highlight the potential of modeling approaches to disentangle species relationships in systems with large effective population sizes and post-divergence gene flow.

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater-resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size-assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.

Linked selection and recombination rate variation drive the evolution of the genomic landscape of differentiation across the speciation continuum of Ficedula flycatchers

Speciation is a continuous process during which genetic changes gradually accumulate in the genomes of diverging species. Recent studies have documented highly heterogeneous differentiation landscapes, with distinct regions of elevated differentiation (“differentiation islands”) widespread across genomes. However, it remains unclear which processes drive the evolution of differentiation islands; how the differentiation landscape evolves as speciation advances; and ultimately, how differentiation islands are related to speciation. Here, we addressed these questions based on population genetic analyses of 200 resequenced genomes from 10 populations of four Ficedula flycatcher sister species. We show that a heterogeneous differentiation landscape starts emerging among populations within species, and differentiation islands evolve recurrently in the very same genomic regions among independent lineages. Contrary to expectations from models that interpret differentiation islands as genomic regions involved in reproductive isolation that are shielded from gene flow, patterns of sequence divergence (d_xy and relative node depth) do not support a major role of gene flow in the evolution of the differentiation landscape in these species. Instead, as predicted by models of linked selection, genome-wide variation in diversity and differentiation can be explained by variation in recombination rate and the density of targets for selection. We thus conclude that the heterogeneous landscape of differentiation in Ficedula flycatchers evolves mainly as the result of background selection and selective sweeps in genomic regions of low recombination. Our results emphasize the necessity of incorporating linked selection as a null model to identify genome regions involved in adaptation and speciation.

Climate adaptation and speciation: particular focus on reproductive barriers in Ficedula flycatchers

Climate adaptation is surprisingly rarely reported as a cause for the build‐up of reproductive isolation between diverging populations. In this review, we summarize evidence for effects of climate adaptation on pre‐ and postzygotic isolation between emerging species with a particular focus on pied (Ficedula hypoleuca) and collared (Ficedula albicollis) flycatchers as a model for research on speciation. Effects of climate adaptation on prezygotic isolation or extrinsic selection against hybrids have been documented in several taxa, but the combined action of climate adaptation and sexual selection is particularly well explored in Ficedula flycatchers. There is a general lack of evidence for divergent climate adaptation causing intrinsic postzygotic isolation. However, we argue that the profound effects of divergence in climate adaptation on the whole biochemical machinery of organisms and hence many underlying genes should increase the likelihood of genetic incompatibilities arising as side effects. Fast temperature‐dependent co‐evolution between mitochondrial and nuclear genomes may be particularly likely to lead to hybrid sterility. Thus, how climate adaptation relates to reproductive isolation is best explored in relation to fast‐evolving barriers to gene flow, while more research on later stages of divergence is needed to achieve a complete understanding of climate‐driven speciation.

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single-nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole-genome resequencing data from population samples, molecular ecology studies will be able to take truly genome-wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single-nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10-fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later-generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.

Genetic, morphological, and acoustic evidence reveals lack of diversification in the colonization process in an island bird

Songbirds with recently (i.e., early Holocene) founded populations are suitable models for studying incipient differentiation in oceanic islands. On such systems each colonization event represents a different evolutionary episode that can be studied by addressing sets of diverging phenotypic and genetic traits. We investigate the process of early differentiation in the spectacled warbler (Sylvia conspicillata) in 14 populations separated by sea barriers from three Atlantic archipelagos and from continental regions spanning from tropical to temperate latitudes. Our approach involved the study of sexual acoustic signals, morphology, and genetic data. Mitochondrial DNA did not provide clear population structure. However, microsatellites analyses consistently identified two genetic groups, albeit without correspondence to subspecies classification and little correspondence to geography. Coalescent analyses showed significant evidence for gene flow between the two genetic groups. Discriminant analyses could not correctly assign morphological or acoustic traits to source populations. Therefore, although theory predicting that in isolated populations genetic, morphological, or acoustic traits can lead to radiation, we have strikingly failed to document differentiation on these attributes in a resident passerine throughout three oceanic archipelagos.

Haldane's rule revisited: do hybrid females have a shorter lifespan? Survival of hybrids in a recent contact zone between two large gull species

Haldane's rule predicts that particularly high fitness reduction should affect the heterogametic sex of interspecific hybrids. Despite the fact that hybridization is widespread in birds, survival of hybrid individuals is rarely addressed in studies of avian hybrid zones, possibly because of methodological constraints. Here, having applied capture-mark-recapture models to an extensive, 19-year-long data set on individually marked birds, we estimate annual survival rates of hybrid individuals in the hybrid zone between herring (Larus argentatus) and Caspian (Larus cachinnans) gulls. In both parental species, males have a slightly higher survival rate than females (model-weighted mean ± SE: herring gull males 0.88 ± 0.01, females 0.87 ± 0.01, Caspian gull males 0.88 ± 0.01, females 0.87 ± 0.01). Hybrid males do not survive for a shorter time than nonhybrid ones (0.88 ± 0.01), whereas hybrid females have the lowest survival rate among all groups of individuals (0.83 ± 0.03). This translates to a shorter adult (reproductive) lifespan (on average by 1.7-1.8 years, i.e. ca 25%) compared with nonhybrid females. We conclude that, in line with Haldane's rule, the lower survival rate of female hybrids may contribute to selection against hybrids in this hybrid zone.

Multilocus analyses indicate a mosaic distribution of hybrid populations in ground squirrels (genus Ictidomys)

DNA sequence data from mitochondrial cytochrome-b (Cytb) and Y-linked structural maintenance of chromosomes (SmcY) genes were combined with 478 nuclear loci obtained from amplified fragment length polymorphisms (AFLP) to assess the extent of hybridization and genetic spatial structure of populations in two hybridizing species of ground squirrel (Ictidomys parvidens and Ictidomys tridecemlineatus). Based on AFLP analyses of 134 individuals from 28 populations, 10 populations were identified that possessed hybrid individuals. Overall estimates of F ST values revealed strong support for population structure in the Cytb data set; however, analyses of the SmcY gene and the AFLP data indicated ongoing gene flow between species. Pairwise F ST comparisons of populations were not significant for the SmcY gene; although they were significant for the Cytb gene, indicating that these populations were structured and that gene flow was minimal. Therefore, gene flow between I. parvidens and I. tridecemlineatus appeared to be restricted to populations that exhibited hybridization. In addition, the fragmented nature of the geographic landscape suggested limited gene flow between populations. As a result, the distributional pattern of interspersed parental and hybrid populations were compatible with a mosaic hybrid zone model. Because ground squirrels display female philopatry and male-biased dispersal, the ecology of these species is compatible with this hypothesis.

Multifaceted, cross-generational costs of hybridization in sibling Drosophila species

Maladaptive hybridization, as determined by the pattern and intensity of selection against hybrid individuals, is an important factor contributing to the evolution of prezygotic reproductive isolation. To identify the consequences of hybridization between Drosophila pseudoobscura and D. persimilis, we estimated multiple fitness components for F1 hybrids and backcross progeny and used these to compare the relative fitness of parental species and their hybrids across two generations. We document many sources of intrinsic (developmental) and extrinsic (ecological) selection that dramatically increase the fitness costs of hybridization beyond the well-documented F1 male sterility in this model system. Our results indicate that the cost of hybridization accrues over multiple generations and reinforcement in this system is driven by selection against hybridization above and beyond the cost of hybrid male sterility; we estimate a fitness loss of >95% relative to the parental species across two generations of hybridization. Our findings demonstrate the importance of estimating hybridization costs using multiple fitness measures from multiple generations in an ecologically relevant context; so doing can reveal intense postzygotic selection against hybridization and thus, an enhanced role for reinforcement in the evolution of populations and diversification of species.

Demographic divergence history of pied flycatcher and collared flycatcher inferred from whole-genome re-sequencing data

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000-80,000) and census sizes (5-50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.

Restriction site-associated DNA sequencing generates high-quality single nucleotide polymorphisms for assessing hybridization between bighead and silver carp in the United States and China

Bighead carp (Hypophthalmichthys nobilis) and silver carp (H. molitrix) are invasive species and listed as US federally injurious species under the Lacy Act. They have established populations in much of the Mississippi River Basin (MRB; Mississippi, Illinois, and Missouri rivers) and are capable of producing fertile hybrids and complex introgression. Characterizing the composition of this admixture requires a large set of high-quality, evolutionarily conserved, diagnostic genetic markers to aid in the identification and management of these species in the midst of morphological ambiguity. Restriction site-associated DNA (RAD) sequencing of 45 barcoded bighead and silver carp from the United States and China produced reads that were aligned to the silver carp transcriptome yielded 261 candidate single nucleotide polymorphisms (SNPs) with fixed allelic differences between the two species. We selected the highest quality 112 SNP loci for validation using 194 putative pure-species and F1 hybrids from the MRB and putative bighead carp and silver carp pure species from China (Amur, Pearl and Yangtze rivers). Fifty SNPs were omitted due to design/amplification failure or lack of diagnostic utility. A total of 57 species-diagnostic SNPs conserved between carp species in US and Chinese rivers were identified; 32 were annotated to functional gene loci. Twenty-seven of the 181 (15%) putative pure species were identified as hybrid backcrosses after validation, including three backcrosses from the Amur River, where hybridization has not been documented previously. The 57 SNPs identified through RAD sequencing provide a diagnostic tool to detect population admixture and to identify hybrid and pure-species Asian carps in the United States and China.

Speciation: clash of the genomes

Complete genomes of hybridizing bird species demonstrate the importance of the sex chromosomes, telomeres and centromeres to the initial stages of speciation.

Low fertility of wild hybrid male flycatchers despite recent divergence

Postzygotic isolation may be important for maintaining species boundaries, particularly when premating barriers are incomplete. Little is known about the course of events leading from minor environmental mismatches affecting hybrid fitness to severe genetic incompatibilities causing sterility or inviability. We investigated whether reduced reproductive success of hybrid males was caused by suboptimal sperm traits or by more severe genetic incompatibilities in a hybrid zone of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) on the island of Öland, Sweden. About 4 per cent hybridization is observed in this population and all female hybrids are sterile. We found no sperm in the ejaculates of most sampled hybrid males, and sperm with abnormal morphology in two hybrids. Furthermore, none of the hybrids sired any offspring because of high levels of hatching failure and extra-pair paternity in their nests. These results from a natural hybrid zone suggest that the spermatogenesis of hybrid males may become disrupted despite little genetic divergence between the parental species.

Relative performance of hybrid nestlings in Ficedula flycatchers: a translocation experiment

Ecological speciation predicts that hybrids should experience relatively low fitness in the local environments of their parental species. In this study, we performed a translocation experiment of nestling hybrids between collared and pied flycatchers into the nests of conspecific pairs of their parental species. Our aim was to compare the performance of hybrids with purebred nestlings. Nestling collared flycatchers are known to beg and grow faster than nestling pied flycatchers under favorable conditions, but to experience higher mortality than nestling pied flycatchers under food limitation. The experiment was performed relatively late in the breeding season when food is limited. If hybrid nestlings have an intermediate growth potential and begging intensity, we expected them to beg and grow faster, but also to experience lower survival than pied flycatchers. In comparison with nestling collared flycatchers, we expected them to beg and grow slower, but to survive better. We found that nestling collared flycatchers indeed begged significantly faster and experienced higher mortality than nestling hybrids. Moreover, nestling hybrids had higher weight and tended to beg faster than nestling pied flycatchers, but we did not detect a difference in survival between the latter two groups of nestlings. We conclude that hybrid Ficedula nestlings appear to have a better intrinsic adaptation to food limitation late in the breeding season compared with nestling collared flycatchers. We discuss possible implications for gene flow between the two species.

Inferring the demographic history of European Ficedula flycatcher populations

Background

Inference of population and species histories and population stratification using genetic data is important for discriminating between different speciation scenarios and for correct interpretation of genome scans for signs of adaptive evolution and trait association. Here we use data from 24 intronic loci re-sequenced in population samples of two closely related species, the pied flycatcher and the collared flycatcher.

Results

We applied Isolation-Migration models, assignment analyses and estimated the genetic differentiation and diversity between species and between populations within species. The data indicate a divergence time between the species of <1 million years, significantly shorter than previous estimates using mtDNA, point to a scenario with unidirectional gene-flow from the pied flycatcher into the collared flycatcher and imply that barriers to hybridisation are still permeable in a recently established hybrid zone. Furthermore, we detect significant population stratification, predominantly between the Spanish population and other pied flycatcher populations.

Conclusions

Our results provide further evidence for a divergence process where different genomic regions may be at different stages of speciation. We also conclude that forthcoming analyses of genotype-phenotype relations in these ecological model species should be designed to take population stratification into account.

The genomic landscape of species divergence in Ficedula flycatchers

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.

Increased divergence but reduced variation on the Z chromosome relative to autosomes in Ficedula flycatchers: differential introgression or the faster-Z effect?

Recent multilocus studies of congeneric birds have shown a pattern of elevated interspecific divergence on the Z chromosome compared to the autosomes. In contrast, intraspecifically, birds exhibit less polymorphism on the Z chromosome relative to the autosomes. We show that the four black-and-white Ficedula flycatcher species show greater genetic divergence on the Z chromosome than on the autosomes, and that the ratios of intraspecific polymorphism at Z-linked versus autosomal markers are below the neutral expectation of 75%. In all species pairs, we found more fixed substitutions and fewer shared polymorphisms on the Z chromosome than on the autosomes. Finally, using isolation with migration (IMa) models we estimated gene flow among the four closely related flycatcher species. The results suggest that different pattern of evolution of Z chromosomes and autosomes is best explained by the faster-Z hypothesis, since the estimated long-term gene flow parameters were close to zero in all comparisons.

A test for environmental effects on behavioral isolation in two species of killifish

Behavioral isolation is a common and potent mechanism of reproductive isolation. Determining the extent to which behavioral isolation varies with environmental conditions is critical to understanding speciation and the maintenance of species boundaries. Here, we tested the effect of salinity on behavioral isolation (female species recognition, male-male competition, and male species recognition) between two closely related killifish (Lucania goodei and L. parva) that differ in salinity tolerance. We performed no-choice assays and behavioral trials where males could compete and court females in fresh water (0 ppt) and brackish water (15 ppt). We found high levels of behavioral isolation that did not vary as a function of salinity. In behavioral trials, male species recognition of females was strong and asymmetric between the two species. Lucania goodei males preferred conspecifics and rarely courted or mated with L. parva females. Lucania parva males preferred conspecifics but readily courted and mated with L. goodei females. This asymmetry matches previously documented asymmetries in hybrid offspring fitness. Crosses between L. parva males and L. goodei females produce fully viable/fertile hybrids, but crosses between L. goodei males and L. parva females produce males with reduced fertility. Hence, behavioral isolation may have evolved in part due to reinforcement.

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.

Positive feedback between ecological and reproductive character displacement in a young avian hybrid zone

Character displacement can reduce costly interspecific interactions between young species. We investigated the mechanisms behind divergence in three key traits-breeding habitat choice, timing of breeding, and plumage coloration-in Ficedula flycatchers. We found that male pied flycatchers became expelled from the preferred deciduous habitat into mixed forest as the superior competitor, collared flycatchers, increased in numbers. The peak in food abundance differs between habitats, and the spatial segregation was paralleled by an increased divergence in timing of breeding between the two species. Male pied flycatchers vary from brown to black with brown coloration being more frequent in sympatry with collared flycatchers, a pattern often proposed to result from selection against hybridization, that is, reinforcement. In contrast to this view, we show that brown male pied flycatchers more often hybridize than black males. Male pied flycatcher plumage coloration influenced the territory obtained in areas of co-occurrence with collared flycatchers, and brown male pied flycatchers experienced higher relative fitness than black males when faced with heterospecific competition. We suggest that allopatric divergence in resource defense ability causes a feedback loop at secondary contact where male pied flycatchers with the most divergent strategy compared to collared flycatchers are favored by selection.

What do we need to know about speciation?

Speciation has been a major focus of evolutionary biology research in recent years, with many important advances. However, some of the traditional organising principles of the subject area no longer provide a satisfactory framework, such as the classification of speciation mechanisms by geographical context into allopatric, parapatric and sympatry classes. Therefore, we have asked where speciation research should be directed in the coming years. Here, we present a distillation of questions about the mechanisms of speciation, the genetic basis of speciation and the relationship between speciation and diversity. Our list of topics is not exhaustive; rather we aim to promote discussion on research priorities and on the common themes that underlie disparate speciation processes.

No evidence for Z-chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene-based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post-zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z-chromosome. Could rearrangements of the Z-chromosome between the species explain their reproductive isolation? We developed high coverage Z-chromosome linkage maps for both species, using gene-based markers and large-scale SNP genotyping. Best order maps contained 57-62 gene markers with an estimated average density of one every 1-1.5 Mb. We estimated the recombination rates in flycatcher Z-chromosomes to 1.1-1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co-linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post-zygotic reproductive isolation between these two Ficedula species.

Speciation in Ficedula flycatchers

Speciation in animals often requires that population divergence goes through three major evolutionary stages, i.e. ecological divergence, development of sexual isolation and the build-up of genetic incompatibility. There is theoretical consensus regarding favourable conditions required for speciation to reach its final and irreversible stage, but empirical tests remain rare. Here, we review recent research on processes of speciation, based on studies in hybrid zones between collared (Ficedula albicollis) and pied flycatchers (Ficedula hypoleuca). A major advantage of this study system is that questions concerning all three major sources of reproductive isolation and their interconnections can be addressed. We conclude that (i) ecological divergence is caused by divergence in life-history traits, (ii) females prefer mates of their own species based on differences in both plumage and song characteristics, (iii) male plumage characteristics have diverged but their song has converged in sympatry, (iv) there is genetic incompatibility in accordance with Haldane's rule, and (v) the Z-chromosome appears to be a hotspot for genes involved in sexual isolation and genetic incompatibility. We discuss how identification of the genes underlying the three major sources of reproductive isolation can be used to draw conclusions about links between the processes driving their evolution.

Introgressive hybridization in North American hakes after secondary contact

High levels of bidirectional introgressive hybridization were found between the two Atlantic North American hakes Merluccius albidus and M. bilinearis in their overlapping distribution area between the 34 and the 42 degrees N parallels, employing mitochondrial and nuclear DNA markers. Absence of F(1) hybrids, and varied levels of bidirectional introgression, indicate long-time hybridization and backcrossing. Based on the evolutionary history of the genus Merluccius, originated in this area from the ancestor of the present M. bilinearis by the rise of the Panama Isthmus, secondary contact between the two species has probably been promoted by northwards displacement of M. albidus. Higher introgression rates in southern areas of M. albidus could be explained by restricted gene flow in that area which may allow long-term accumulation of introgressed genes.