Alternative forms for genomic clines

Understanding factors regulating hybrid fitness and gene exchange is a major research challenge for evolutionary biology. Genomic cline analysis has been used to evaluate alternative patterns of introgression, but only two models have been used widely and the approach has generally lacked a hypothesis testing framework for distinguishing effects of selection and drift. I propose two alternative cline models, implement multivariate outlier detection to identify markers associated with hybrid fitness, and simulate hybrid zone dynamics to evaluate the signatures of different modes of selection. Analysis of simulated data shows that previous approaches are prone to false positives (multinomial regression) or relatively insensitive to outlier loci affected by selection (Barton's concordance). The new, theory‐based logit‐logistic cline model is generally best at detecting loci affecting hybrid fitness. Although some generalizations can be made about different modes of selection, there is no one‐to‐one correspondence between pattern and process. These new methods will enhance our ability to extract important information about the genetics of reproductive isolation and hybrid fitness. However, much remains to be done to relate statistical patterns to particular evolutionary processes. The methods described here are implemented in a freely available package “HIest” for the R statistical software (CRAN; http://cran.r-project.org/ ).     

bgc: Software for Bayesian estimation of genomic clines

Introgression in admixed populations can be used to identify candidate loci that might underlie adaptation or reproductive isolation. The Bayesian genomic cline model provides a framework for quantifying variable introgression in admixed populations and identifying regions of the genome with extreme introgression that are potentially associated with variation in fitness. Here we describe the bgc software, which uses Markov chain Monte Carlo to estimate the joint posterior probability distribution of the parameters in the Bayesian genomic cline model and designate outlier loci. This software can be used with next‐generation sequence data, accounts for uncertainty in genotypic state, and can incorporate information from linked loci on a genetic map. Output from the analysis is written to an HDF5 file for efficient storage and manipulation. This software is written in C++ . The source code, software manual, compilation instructions and example data sets are available under the GNU Public License at http://sites.google.com/site/bgcsoftware/ .     

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.     

The genetics of reproductive isolation and the potential for gene exchange between Drosophila pseudoobscura and D. persimilis via backcross hybrid males

Hybrid male sterility, hybrid inviability, sexual isolation, and a hybrid male courtship dysfunction reproductively isolate Drosophila pseudoobscura and D. persimilis. Previous studies of the genetic bases of these isolating mechanisms have yielded only limited information about how much and what areas of the genome are susceptible to interspecies introgression. We have examined the genetic basis of these barriers to gene exchange in several thousand backcross hybrid male progeny of these species using 14 codominant molecular genetic markers spanning the five chromosomes of these species, focusing particularly on the autosomes. Hybrid male sterility, hybrid inviability, and the hybrid male courtship dysfunction were all associated with X-autosome interactions involving primarily the inverted regions on the left arm of the X-chromosome and the center of the second chromosome. Sexual isolation from D. pseudoobscura females was primarily associated with the left arm of the X-chromosome, although both the right arm and the center of the second chromosome also contributed to it. Sexual isolation from D. persimilis females was primarily associated with the second chromosome. The absence of isolating mechanisms being associated with many autosomal regions, including some large inverted regions that separate the strains, suggests that these phenotypes may not be caused by genes spread throughout the genome. We suggest that gene flow between these species via hybrid males may be possible at loci spread across much of the autosomes.     

Low intraspecific variation for genomic isolation between hybridizing sunflower species

Barriers to gene flow between species result from selection against foreign linkage blocks in hybrids. When the geographic ranges of taxa meet at multiple locations, the opportunity exists for variation in the genetic architecture of isolating barriers. Hybrid zones between two sunflower species (Helianthus annuus and H. petiolaris) in Nebraska and California exhibited remarkably similar patterns of introgression of mapped molecular markers. Congruence among hybrid zones may result from limited intraspecific variation at loci contributing to isolation and from similar selective effects of alleles in the heterospecific genetic background. The observed consistency of introgression patterns across distantly separated hybrid zones suggests that intrinsic forces predominate in determining hybrid zone dynamics and boundaries between these sunflower species.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

The genomic and ecological context of hybridization affects the probability that symmetrical incompatibilities drive hybrid speciation

Despite examples of homoploid hybrid species, theoretical work describing when, where, and how we expect homoploid hybrid speciation to occur remains relatively rare. Here, I explore the probability of homoploid hybrid speciation due to “symmetrical incompatibilities” under different selective and genetic scenarios. Through simulation, I test how genetic architecture and selection acting on traits that do not themselves generate incompatibilities interact to affect the probability that hybrids evolve symmetrical incompatibilities with their parent species. Unsurprisingly, selection against admixture at “adaptive” loci that are linked to loci that generate incompatibilities tends to reduce the probability of evolving symmetrical incompatibilities. By contrast, selection that favors admixed genotypes at adaptive loci can promote the evolution of symmetrical incompatibilities. The magnitude of these outcomes is affected by the strength of selection, aspects of genetic architecture such as linkage relationships and the linear arrangement of loci along a chromosome, and the amount of hybridization following the formation of a hybrid zone. These results highlight how understanding the nature of selection, aspects of the genetics of traits affecting fitness, and the strength of reproductive isolation between hybridizing taxa can all be used to inform when we expect to observe homoploid hybrid speciation due to symmetrical incompatibilities.     

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.     

Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Hybrid zones allow the measurement of gene flow across the genome, producing insight into the genomic architecture of speciation. Such analysis is particularly powerful when applied to multiple pairs of hybridizing species, as patterns of genomic differentiation can then be related to age of the hybridizing species, providing a view into the build‐up of differentiation over time. We examined 33 809 single nucleotide polymorphisms (SNPs) in three hybridizing woodpecker species: Red‐breasted, Red‐naped and Yellow‐bellied sapsuckers (Sphyrapicus ruber, Sphyrapicus nuchalis and Sphyrapicus varius), two of which (ruber and nuchalis) are much more closely related than each is to the third (varius). To identify positions of SNPs on chromosomes, we developed a localization method based on comparative genomics. We found narrow clines, bimodal distributions of hybrid indices and genomic regions with decreased rates of introgression. These results suggest moderately strong reproductive isolation among species and selection against specific hybrid genotypes. We found 19 small regions of strong differentiation between species, partly shared among species pairs, but no large regions of differentiation. An association analysis revealed a single strong‐effect candidate locus associated with plumage, possibly explaining mismatch among the three species in genomic relatedness and plumage similarity. Our comparative analysis of species pairs of different age and their hybrid zones showed that moderately strong reproductive isolation can occur with little genomic differentiation, but that reproductive isolation is incomplete even with much greater genomic differentiation, implying there are long periods of time when hybridization is possible if diverging populations are in geographic contact.     

introgress: a software package for mapping components of isolation in hybrids

A new software package (introgress) provides functions for analysing introgression of genotypes between divergent, hybridizing lineages, including estimating genomic clines from multi-locus genotype data and testing for deviations from neutral expectations. The software works with co-dominant, dominant and haploid marker data, and does not require fixed allelic differences between parental populations for the sampled genetic markers. Permutation and parametric procedures generate neutral expectations for introgression and provide a basis for significance tests of observed genomic clines. The software also implements maximum likelihood estimates of hybrid index from genotypic data and a number of graphical analyses. The package is an extension of the R statistical software, is written in the R language and is freely available through the Comprehensive R Archive Network (CRAN; http://cran.r-project.org/). In this study, we describe introgress and demonstrate its use with a sample data set.     

Strong postmating reproductive isolation in Mimulus section Eunanus

Postmating reproductive isolation can help maintain species boundaries when premating barriers to reproduction are incomplete. The strength and identity of postmating reproductive barriers are highly variable among diverging species, leading to questions about their genetic basis and evolutionary drivers. These questions have been tackled in model systems but are less often addressed with broader phylogenetic resolution. In this study we analyse patterns of genetic divergence alongside direct measures of postmating reproductive barriers in an overlooked group of sympatric species within the model monkeyflower genus, Mimulus. Within this Mimulus brevipes species group, we find substantial divergence among species, including a cryptic genetic lineage. However, rampant gene discordance and ancient signals of introgression suggest a complex history of divergence. In addition, we find multiple strong postmating barriers, including postmating prezygotic isolation, hybrid seed inviability and hybrid male sterility. M. brevipes and M. fremontii have substantial but incomplete postmating isolation. For all other tested species pairs, we find essentially complete postmating isolation. Hybrid seed inviability appears linked to differences in seed size, providing a window into possible developmental mechanisms underlying this reproductive barrier. While geographic proximity and incomplete mating isolation may have allowed gene flow within this group in the distant past, strong postmating reproductive barriers today have likely played a key role in preventing ongoing introgression. By producing foundational information about reproductive isolation and genomic divergence in this understudied group, we add new diversity and phylogenetic resolution to our understanding of the mechanisms of plant speciation.     

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.     

Structure of a mosaic hybrid zone between the field crickets Gryllus firmus and G. pennsylvanicus

Hybrid zones provide insight into the nature of species boundaries and the evolution of barriers to gene exchange. Characterizing multiple regions within hybrid zones is essential for understanding both their history and current dynamics. Here, we describe a previously uncharacterized region of a well‐studied hybrid zone between two species of field crickets, Gryllus pennsylvanicus and G. firmus. We use a combination of mitochondrial DNA sequencing, morphological data, and modeling of environmental variables to identify the ecological factors structuring the hybrid zone and define patterns of hybridization and introgression. We find an association between species distribution and natural habitat; Gryllus pennsylvanicus occupies natural habitat along forest edges and natural clearings, whereas G. firmus occupies more disturbed areas in agricultural and suburban environments. Hybridization and introgression occur across patch boundaries; there is evidence of substantial admixture both in morphological characters and mtDNA, over a broad geographic area. Nonetheless, the distribution of morphological types is bimodal. Given that F₁ hybrids are viable and fertile in the lab, this suggests that strong pre‐zygotic barriers are operating in this portion of the hybrid zone.     

Maximum‐likelihood estimation of a hybrid index based on molecular markers

This paper elaborates on a hybrid index that utilizes information from genetic markers to quantify the genetic contribution of hybridizing species to individuals of unknown ancestry. Dominant markers will only lead to reliable and accurate estimates of hybrid index in later generation hybrids. In contrast, codominant markers can be fully resolved and their use is unproblematic. For both types of markers and allele frequencies that differ substantially between parental species (F_ST ≥ 0.17), a hybrid index based on 35–45 loci will have a nearly minimal confidence interval. Estimates of hybrid index are robust to modest errors in estimates of parental allele frequencies.     

Hybridization and restricted gene flow between native and introduced stocks of Alpine whitefish (Coregonus sp.) across multiple environments

Translocations of Baltic whitefish (Coregonus sp.) into Austrian Alpine lakes have created 'artificial hybrid zones', threatening the genetic integrity of native lineages. We evaluate the genetic structure of Coregonus in Austrian lakes and characterize hybridization and introgression between native and introduced lineages. Fifteen populations (N=747) were assessed for allelic variation at eight microsatellite loci and a reduced set (N=253) for variation across two mtDNA genes (cyt b and NADH-3). Bayesian approaches were used to estimate individual admixture proportions (q-values) and classify genotypes as native, introduced or hybrids. q-value distributions varied among populations highlighting differential hybridization and introgression histories. Many lakes revealed a clear distinction between native and introduced genotypes despite hybridization, whereas some locations revealed hybrid swarms. Genetic structure among lakes was congruent with morphological divergence and novelty raising speculation of multiple taxa, including a population south of the Alps, outside the putative native range of Coregonus. Although statistically congruent with inferences based on nuclear markers, mitochondrial haplotype data was not diagnostic with respect to native and non-native lineages, supporting that the Alpine region was colonized post-glacially by an admixture of mtDNA lineages, which coalesce >1 Ma. Mechanisms promoting or eroding lineage isolation are discussed, as well as a high potential to conserve native Alpine lineages despite the extensive historical use of introduced Baltic stocks.     

Discordant divergence times among Z-chromosome regions between two ecologically distinct swallowtail butterfly species

We investigate multilocus patterns of differentiation between parental populations of two swallowtail butterfly species that differ at a number of ecologically important sex-linked traits. Using a new coalescent-based approach, we show that there is significant heterogeneity in estimated divergence times among five Z-linked markers, rejecting a purely allopatric speciation model. We infer that the Z chromosome is a mosaic of regions that differ in the extent of historical gene flow, potentially due to isolating barriers that prevent the introgression of species-specific traits that result in hybrid incompatibilities. Surprisingly, a candidate region for a strong barrier to introgression, Ldh, does not show a significantly deeper divergence time than other markers on the Z chromosome. Our approach can be used to test alternative models of speciation and can potentially assign chronological order to the appearance of factors contributing to reproductive isolation between 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.     

A powerful regression-based method for admixture mapping of isolation across the genome of hybrids

We propose a novel method that uses natural admixture between divergent lineages (hybridization) to investigate the genetic architecture of reproductive isolation and adaptive introgression. Our method employs multinomial regression to estimate genomic clines and to quantify introgression for individual loci relative to the genomic background (clines in genotype frequency along a genomic admixture gradient). Loci with patterns of introgression that deviate significantly from null expectations based on the remainder of the genome are potentially subject to selection and thus of interest to understanding adaptation and the evolution of reproductive isolation. Using simulations, we show that different forms of selection modify these genomic clines in predictable ways and that our method has good power to detect moderate to strong selection for multiple forms of selection. Using individual-based simulations, we demonstrate that our method generally has a low false positive rate, except when genetic drift is particularly pronounced (e.g. low population size, low migration rates from parental populations, and substantial time since initial admixture). Additional individual-based simulations reveal that moderate selection against heterozygotes can be detected as much as 50 c m away from the focal locus directly experiencing selection, but is not detected at unlinked loci. Finally, we apply our analytical method to previously published data sets from a mouse ( Mus musculus and M. domesticus ) and two sunflower ( Helianthus petiolaris and H. annuus ) hybrid zones. This method should be applicable to numerous species that are currently the focus of research in evolution and ecology and should help bring about new insights regarding the processes underlying the origin and maintenance of biological diversity.     

Hybridization following recent secondary contact results in asymmetric genotypic and phenotypic introgression between island species of Myzomela honeyeaters

Hybridization and introgression can have important evolutionary consequences for speciation, especially during early stages of secondary contact when reproductive barriers may be weak. Few studies, however, have quantified dynamics of hybridization and introgression in systems in which recent natural dispersal across a geographic barrier resulted in secondary contact. We investigated patterns of hybridization and introgression between two Myzomela honeyeaters (M. tristrami and M. cardinalis) that recently achieved secondary contact on Makira in the Solomon Islands. Hybridization in this system was hypothesized to be a byproduct of conspecific mate scarcity during early stages of colonization. Our research, however, provides evidence of ongoing hybridization more than a century after secondary contact. Mitochondrial sequencing revealed strongly asymmetric reproductive isolation that is most likely driven by postzygotic incompatibilities rather than prezygotic behavioral barriers. Nuclear introgression was observed from the native species (M. tristrami) to the colonizing species (M. cardinalis). Nuclear introgression in the reverse direction is almost exclusively limited to birds that are phenotypically M. tristrami but possess M. cardinalis mitochondrial haplotypes, consistent with introgression of plumage‐related alleles into the genomic background of M. cardinalis. These results provide unique insight into the dynamics and consequences of hybridization and introgression during early stages of secondary contact.