Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Origin(s) of the diploid hybrid species Helianthus deserticola (Asteraceae)

Homoploid hybrid speciation has traditionally been considered a rare event, dependent on the establishment of both a novel, balanced genotype and reproductive isolating barriers between the new species and its progenitors. However, more recent studies have shown that synthetic hybrids converge toward the chromosomal structure of natural hybrids after only a few generations, suggesting that this phenomenon may be more frequent than previously assumed. Here, the possibility that the diploid hybrid species Helianthus deserticola arose from more than one hybrid speciation event was investigated using patterns of variation from cpDNA, 18 nuclear microsatellite loci, and population interfertility. Helianthus deserticola contains cpDNA haplotypes characteristic of both parental species, is polyphyletic with one parental species based on nine microsatellite loci, and has a high degree of interfertility among populations. The data are consistent with either a single origin followed by introgression with the parental species or multiple origins. Analysis of microsatellite variation places the origin of H. deserticola between 170?000 and 63?000 years before present, making it unlikely that anthropogenic disturbances influenced its origin. Finally, the hybrid species generally has lower levels of genetic diversity but higher levels of differentiation among populations than either parental species.     

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.     

Multilocus test for introgression between the cactophilic species Drosophila mojavensis and Drosophila arizonae

Information obtained from laboratory studies regarding the efficacy of barriers to gene flow (reproductive isolation) between species is often incomplete or misleading, so detailed genetic analyses are needed to determine whether hybridization and introgression occur in nature. Previous laboratory studies of the cactophilic species Drosophila mojavensis and Drosophila arizonae suggest that reproductive isolation is incomplete and that gene flow may occur in sympatry. We sampled 18 nuclear and one mitochondrial loci from multiple populations of D. arizonae and D. mojavensis to test for the signature of recent or historic gene flow between these two species. We located chromosomal regions that were inverted between these species and analyzed those regions independently of others. Statistical tests for introgression using all loci or only collinear loci failed to reject expectations of an isolation model. Further tests using average nucleotide differences between species and phylogenetic analyses also failed to find support for introgression between D. mojavensis and D. arizonae. Additional ecological and behavioral studies of these species in their natural habitats are required to explain why the signature of gene flow was not detected at the DNA sequence level in populations when laboratory studies suggest such gene flow should be possible.     

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

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.     

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.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Molecular and morphological patterns of introgression between two large white-headed gull species in a zone of recent secondary contact

Incomplete reproductive isolation promotes gene flow between diverging taxa. However, any gene encoding for traits involved in the reproductive barriers will be less prone to introgression than neutral markers. Comparing introgression rates among loci is thus informative of the number and functions of loci involved in the reproductive barriers. This study aimed at identifying possible mechanisms of restriction to gene flow across a zone of recent secondary contact between Larus argentatus and Larus cachinnans by comparing introgression patterns for nine microsatellite loci, a fragment of mitochondrial DNA and a set of phenotypic traits. The low linkage disequilibrium between neutral nuclear markers indicated introgression without any barrier to gene flow. However, asymmetric introgression of mitochondrial DNA suggested that interspecific crosses may be more successful in one direction. The introgression rate for phenotypic traits was variable and low compared to neutral molecular markers. This was particularly evident in colouration of bare parts: individuals with intermediate colouration were scarcer in sympatry than expected if the genomes recombined freely. We hypothesized that one of these variables, the orbital ring colour, may play a role in mate choice, acting as an incomplete premating barrier through assortative mating. This study emphasizes that multilocus approaches are useful to discriminate among possible mechanisms responsible for the maintenance of hybrid zones.     

The rates of introgression and barriers to genetic exchange between hybridizing species: sex chromosomes vs autosomes

Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.     

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.     

Human disturbance causes the formation of a hybrid swarm between two naturally sympatric fish species

Most evidence for hybrid swarm formation stemming from anthropogenic habitat disturbance comes from the breakdown of reproductive isolation between incipient species, or introgression between allopatric species following secondary contact. Human impacts on hybridization between divergent species that naturally occur in sympatry have received considerably less attention. Theory predicts that reinforcement should act to preserve reproductive isolation under such circumstances, potentially making reproductive barriers resistant to human habitat alteration. Using 15 microsatellites, we examined hybridization between sympatric populations of alewife (Alosa pseudoharengus) and blueback herring (A. aestivalis) to test whether the frequency of hybridization and pattern of introgression have been impacted by the construction of a dam that isolated formerly anadromous populations of both species in a landlocked freshwater reservoir. The frequency of hybridization and pattern of introgression differed markedly between anadromous and landlocked populations. The rangewide frequency of hybridization among anadromous populations was generally 0–8%, whereas all landlocked individuals were hybrids. Although neutral introgression was observed among anadromous hybrids, directional introgression leading to increased prevalence of alewife genotypes was detected among landlocked hybrids. We demonstrate that habitat alteration can lead to hybrid swarm formation between divergent species that naturally occur sympatrically, and provide empirical evidence that reinforcement does not always sustain reproductive isolation under such circumstances.     

GENEALOGICAL DISCORDANCE AND PATTERNS OF INTROGRESSION AND SELECTION ACROSS A CRICKET HYBRID ZONE

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus , offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.     

Limited Introgression between Rock-Wallabies with Extensive Chromosomal Rearrangements

Chromosome rearrangements can result in the rapid evolution of hybrid incompatibilities. Robertsonian fusions, particularly those with monobrachial homology, can drive reproductive isolation amongst recently diverged taxa. The recent radiation of rock-wallabies (genus Petrogale) is an important model to explore the role of Robertsonian fusions in speciation. Here, we pursue that goal using an extensive sampling of populations and genomes of Petrogale from north-eastern Australia. In contrast to previous assessments using mitochondrial DNA or nuclear microsatellite loci, genomic data are able to separate the most closely related species and to resolve their divergence histories. Both phylogenetic and population genetic analyses indicate introgression between two species that differ by a single Robertsonian fusion. Based on the available data, there is also evidence for introgression between two species which share complex chromosomal rearrangements. However, the remaining results show no consistent signature of introgression amongst species pairs and where evident, indicate generally low introgression overall. X-linked loci have elevated divergence compared with autosomal loci indicating a potential role for genic evolution to produce reproductive isolation in concert with chromosome change. Our results highlight the value of genome scale data in evaluating the role of Robertsonian fusions and structural variation in divergence, speciation, and patterns of molecular evolution.     

Using differential introgression in hybrid zones to identify genomic regions involved in speciation

Hybrids between species provide information about the evolutionary processes involved in divergence. In addition to creating hybrids in the laboratory, biologists can take advantage of natural hybrid zones to understand the factors that shape gene flow between divergent lineages. In the early stages of speciation, most regions of the genome continue to flow freely between populations. Alternatively, the subset of the genome that confers reproductive barriers between nascent species is expected to reject introgression. Now enabled by advances in genomics, this perspective is motivating detailed comparisons of gene flow across genomic regions in hybrid zones. Here, I review methods for measuring and interpreting introgression at multiple loci in hybrid zones, focusing on the problem of identifying loci that contribute to reproductive isolation. Emerging patterns from multi-locus studies of hybrid zones are highlighted, including remarkable variance in introgression across the genome. Although existing methods have been useful, there is scope for development of new analytical approaches that better connect differential patterns of gene flow in hybrid zones with current knowledge of speciation mechanisms. I outline future prospects for differential introgression studies on a genomic scale.     

Long Distance Linkage Disequilibrium and Limited Hybridization Suggest Cryptic Speciation in Atlantic Cod

Hybrid zones provide unprecedented opportunity for the study of the evolution of reproductive isolation, and the extent of hybridization across individuals and genomes can illuminate the degree of isolation. We examine patterns of interchromosomal linkage disequilibrium (ILD) and the presence of hybridization in Atlantic cod, Gadus morhua, in previously identified hybrid zones in the North Atlantic. Here, previously identified clinal loci were mapped to the cod genome with most (∼70%) occurring in or associated with (<5 kb) coding regions representing a diverse array of possible functions and pathways. Despite the observation that clinal loci were distributed across three linkage groups, elevated ILD was observed among all groups of clinal loci and strongest in comparisons involving a region of low recombination along linkage group 7. Evidence of ILD supports a hypothesis of divergence hitchhiking transitioning to genome hitchhiking consistent with reproductive isolation. This hypothesis is supported by Bayesian characterization of hybrid classes present and we find evidence of common F1 hybrids in several regions consistent with frequent interbreeding, yet little evidence of F2 or backcrossed individuals. This work suggests that significant barriers to hybridization and introgression exist among these co-occurring groups of cod either through strong selection against hybrid individuals, or genetic incompatibility and intrinsic barriers to hybridization. In either case, the presence of strong clinal trends, and little gene flow despite extensive hybridization supports a hypothesis of reproductive isolation and cryptic speciation in Atlantic cod. Further work is required to test the degree and nature of reproductive isolation in this species.     

Speciation in Passerina buntings: introgression patterns of sex-linked loci identify a candidate gene region for reproductive isolation

Sex-chromosomes are thought to play an important role in speciation, but few studies of non-model organisms have investigated the relative influence of multiple sex-linked markers on reproductive isolation. We collected 222 individuals along a geographical transect spanning the hybrid zone between Passerina amoena and P. cyanea (Aves: Cardinalidae). Using maximum-likelihood cline fitting methods, we estimated locus-specific introgression rates for 10 z-linked markers. Although the cline width estimates ranged from 2.8 to 584 km, eight of 10 loci had cline widths between 224 and 271 km. We also used coalescent-based estimates of locus-specific divergence times between P. amoena and P. cyanea to test a recently proposed hypothesis of an inverse relationship between divergence time and cline width but did not find a significant association. The narrow width (2.8 km) of the cline estimated from the VLDLR9 locus indicates strong selection retarding introgression of alleles at this locus across the hybrid zone. Interestingly, a mutation in the very low density lipoprotein receptor ( VLDLR ) gene, in which VLDLR9 is an intron, is known to reduce the egg-laying ability of some chickens, suggesting a possible link between this gene region and reproductive isolation between P. amoena and P. cyanea . These results underscore the importance of sampling multiple loci to investigate introgression patterns across a chromosome or genome and support previous findings of the importance of sex-linked genes in speciation.     

Strong selection on male plumage in a hybrid zone between a hybrid bird species and one of its parents

Homoploid hybrid speciation (HHS) requires reproductive barriers between hybrid and parent species, despite incomplete reproductive isolation (RI) between the parents. Novel secondary sexual trait values in hybrids may cause prezygotic isolation from both parents, whereas signals inherited by the hybrid from one parent species may cause prezygotic isolation with the other. Here we investigate whether differences in male plumage function as a premating barrier between the hybrid Italian sparrow and one of its parent species, the house sparrow, in a narrow Alpine hybrid zone. Italian sparrow male plumage is a composite mosaic of the parental traits, with its head plumage most similar to its other parent, the Spanish sparrow. We use geographical cline analysis to examine selection on three plumage traits, 75 nuclear single nucleotide polymorphisms (SNPs) and hybrid indices based on these SNPs. Several SNPs showed evidence of restricted introgression in the Alps, supporting earlier findings. Crown colour exhibited the narrowest plumage cline, representing a 37% (range 4–65%) drop in fitness. The cline was too narrow to be due to neutral introgression. Only crown colour was significantly bimodal in the hybrid zone. Bimodality may be due to RI or a major QTL, although fitness estimates suggest that selection contributes to the pattern. We discuss the implications with respect to HHS and the species status of the Italian sparrow.     

Formation of a recent hybrid zone offers insight into the geographic puzzle and maintenance of species boundaries in musk turtles

Speciation is the result of an accumulation of reproductive barriers between populations, but pinpointing the factors involved is often difficult. However, hybrid zones can form when these barriers are not complete, especially when lineages come into contact in intermediate or modified habitats. We examine a hybrid zone between two closely related riverine turtle species, Sternotherus depressus and S. peltifer, and use dual‐digest RAD sequencing to understand how this hybrid zone formed and elucidate genomic patterns of reproductive isolation. First, the geographical extent and timing of formation of the hybrid zone is established to provide context for understanding the role of extrinsic and intrinsic reproductive isolating mechanisms in this system. The strength of selection on taxon‐specific contributions to maintenance of the hybrid zone is then inferred using a Bayesian genomic cline model. These analyses identify a role for selection inhibiting introgression in some genomic regions at one end of the hybrid zone and promoting introgression in many loci at the other. When selective pressures necessary to generate outliers to the genomic cline are considered with the geographical and temporal context of this hybrid zone, we conclude that habitat‐specific selection probably limits introgression from S. depressus to S. peltifer in the direction of river flow. However, selection is mediating rapid, unidirectional introgression from S. peltifer to S. depressus, which is probably facilitated by anthropogenic habitat alteration. These findings indicate a potentially imminent threat of population‐level genomic extinction for an already imperiled species due to ongoing human‐caused habitat alteration.