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.     

Evidence of linked selection on the Z chromosome of hybridizing hummingbirds*

Levels of genetic differentiation vary widely along the genomes of recently diverged species. What processes cause this variation? Here, I analyze geographic population structure and genome-wide patterns of variation in the Rufous, Allen's, and Calliope Hummingbirds (Selasphorus rufus/Selasphorus sasin/Selasphorus calliope) and assess evidence that linked selection on the Z chromosome drives patterns of genetic differentiation in a pair of hybridizing species. Demographic models, introgression tests, and genotype clustering analyses support a reticulate evolutionary history consistent with divergence during the late Pleistocene followed by gene flow across migrant Rufous and Allen's Hummingbirds during the Holocene. Relative genetic differentiation () is elevated, and within-population diversity (π) is depressed on the Z chromosome in all interspecific comparisons. The ratio of Z to autosomal within-population diversity is much lower than that expected from population size effects alone, and Tajima's D is depressed on the Z chromosome in S. rufus and S. calliope. These results suggest that conserved structural features of the genome play a prominent role in shaping genetic differentiation through the early stages of speciation in northern Selasphorus hummingbirds, and that the Z chromosome is a likely site of genes underlying behavioral and morphological variation in the group.     

Multilocus analyses of admixture and introgression among hybridizing Heliconius butterflies

Introgressive hybridization is an important evolutionary process and new analytical methods provide substantial power to detect and quantify it. In this study we use variation in the frequency of 657 AFLP fragments and DNA sequence variation from 15 genes to measure the extent of admixture and the direction of interspecific gene flow among three Heliconius butterfly species that diverged recently as a result of natural selection for Miillerian mimicry, and which continue to hybridize. Bayesian clustering based on AFLP genotypes correctly delineated the three species and identified four H. cydno, three H. pachinus, and three H. melpomene individuals that were of mixed ancestry. Gene genealogies revealed substantial shared DNA sequence variation among all three species and coalescent simulations based on the Isolation with Migration (IM) model pointed to interspecific gene flow as its cause. The IM simulations further indicated that interspecific gene flow was significantly asymmetrical, with greater gene flow from H. pachinus into H. cydno (2Nm = 4.326) than the reverse (2Nm = 0.502), and unidirectional gene flow from H. cydno and H. pachinus into H. melpomene (2Nm = 0.294 and 0.252, respectively). These asymmetries are in the directions expected based on the genetics of wing patterning and the probability that hybrids of various phenotypes will survive and reproduce in different mimetic environments. This empirical demonstration of extensive interspecific gene flow is in contrast to a previous study which found little evidence of gene flow between another pair of hybridizing Heliconius species, H. himera and H. erato, and it highlights the critical role of natural selection in maintaining species diversity. Furthermore, these results lend support to the hypotheses that phenotypic diversification in the genus Heliconius has been fueled by introgressive hybridization and that reinforcement has driven the evolution of assortative mate preferences.     

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.     

Evidence for bimodal hybrid zones between two species of char (Pisces: Salvelinus) in northwestern North America

Dolly Varden (Salvelinus malma, Pisces: Salmonidae) and bull trout (Salvelinus confluentus) have widely overlapping, but largely parapatric ranges in watersheds in northwestern North America from Washington State to northern British Columbia. Genetic analysis of natural populations using diagnostic molecular markers revealed widespread local sympatry and hybridization with hybrids comprising 0-25% of the local samples. In a detailed analysis of hybridization using four nuclear DNA markers and mitochondrial DNA within the Thutade Lake watershed, northcentral British Columbia, hybrid genotypes constituted up to 9% of the population of juvenile char. There were significant deviations from Hardy-Weinberg, gametic, and cytonuclear equilibria, and local samples showed bimodal frequency distributions of genotypes. Pure parental and inferred backcross genotypes were most common, and F1 and F(n) hybrids were comparatively rare. Interspecific hybridization was asymmetrical, with most F1 hybrids (five of six) bearing S. confluentus mtDNA. The introgression of nuclear and mitochondrial alleles was asymmetrical, with S. confluentus mtDNA and Growth Hormone 2 introgressing into S. malma significantly more than either introgression of the three other nuclear loci, or introgression of S. malma alleles into S. confluentus. Substantial prezygotic isolation between the species likely depends on the large body size difference between them in sympatry: S. malma have small bodies and a stream resident life history (12-21 cm adult fork length at maturity), while S. confluentus are larger and adfluvial, i.e., they migrate to Thutade Lake where they grow to maturity before returning to tributary streams to spawn (40-90 cm at maturity). These traits may limit interspecific pairings because of size assortative pairing and size-dependent reproductive habitat use.     

Interspecies mating in sympatric species of Syngnathus pipefish

Hybridization is thought to be an important source of novel genetic variation, and interspecific hybridization may increase the adaptive potential of wild populations. While hybridization has not been previously reported in syngnathid fishes (seahorses and pipefish), the sympatric occurrence of closely related species at high densities increases the probability of interspecies mating in this group. Southern California is home to five species of Syngnathus pipefish, and these species frequently co-occur in near-shore eelgrass beds along the California coast. Recent work has identified exceptionally high levels of genetic diversity in southern populations of Syngnathus leptorhynchus, a widespread species which ranges from Mexico to Alaska. Microsatellite genotyping and mitochondrial sequence data are used here to study the population genetics of S. leptorhynchus and S. auliscus at a site in San Diego Bay where they are found to co-occur at high densities. While no adult hybrids were detected in the study population, analysis of male broods indicates that interspecies mating is occurring between the two species. The lack of premating isolating mechanisms between these two relatives suggests that hybridization may be common in sympatric species of Syngnathus.     

Gene flow across a climatic barrier between hybridizing avian species, California and Gambel's quail (Callipepla californica and C. gambelii)

Allopatric species commonly interbreed in a restricted margin between their ranges. The particular factors that permit interbreeding between species determine the extent of hybridization and its significance for evolution and conservation. Using California quail and Gambel's quail (Callipepla californica and C. gambelii) that naturally hybridize in a narrow region between relatively mesic and xeric environments, I assessed the exchange of genetic and phenotypic traits in relation to vegetative and climatic features (temperature and precipitation) that characterize the area of range overlap, and I examined genetic and phenotypic traits within the hybrid zone over a five-year period in relation to variation in precipitation. Using microsatellite markers, this study reveals that genetic, plumage, and morphometric traits are tightly associated with vegetation, rainfall, and temperature profiles through the abrupt transition from one parental species to the other across the hybrid zone. Results show that the hybrid zone has remained clinal, stationary, and bounded over the five-year study period. There was no evidence of introgression outside the narrow hybrid zone. Interannual climatic fluctuations are associated with internal hybrid zone dynamics but did not alter the shape and position of the zone. A transect through the hybrid zone revealed rapid and episodic genetic mixing within the zone. Possible long-term consequences of this restricted hybridization for the evolution of the two parental species are discussed in the light of changing environments.     

Genetic differentiation without mimicry shift in a pair of hybridizing Heliconius species (Lepidoptera: Nymphalidae)

Butterflies in the genus Heliconius have undergone rapid adaptive radiation for warning patterns and mimicry, and are excellent models to study the mechanisms underlying diversification. In Heliconius, mimicry rings typically involve distantly related species, whereas closely related species often join different mimicry rings. Genetic and behavioural studies have n how reproductive isolation in many pairs of Heliconius taxa is largely mediated by natural and sexual selection on wing colour patterns. However, recent studies have uncovered new cases in which pairs of closely related species are near‐perfect mimics of each other. Here, we provide morphometric and genetic evidence for the coexistence of two closely related, hybridizing co‐mimetic species on the eastern slopes of the Andes, H. melpomene amaryllis and H. timareta ssp. nov. , which is described here as H. timareta thelxinoe . A joint analysis of multilocus genotyping and geometric morphometrics of wing shape shows a high level of differentiation between the two species, with only limited gene flow and mixing. Some degree of genetic mixing can be detected, but putative hybrids were rare, only one of 175 specimens being a clear hybrid. In contrast, we found phenotypic differentiation between populations of H. timareta thelxinoe , possibly indicative of strong selection for local mimicry in different communities. In this pair of species, the absence of breakdown of genetic isolation despite near‐identical wing patterns implies that factors other than wing patterns keep the two taxa apart, such as chemical or behavioural signals, or ecological adaptation along a strong altitudinal gradient. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 830–847.     

Hybridization among sympatric species of Rhododendron (Ericaceae) in Turkey: morphological and molecular evidence

Rhododendron (Ericaceae) is a large genus in which barriers to hybridization are especially weak, but many species are maintained in sympatry. Hybridization among four species of Rhododendron subsect. Pontica, which occur in sympatry in Turkey, was investigated. Material of R. ponticum, R. smirnovii, R. ungernii, and R. caucasicum and their putative hybrids was collected from the wild. Based on morphology, chloroplast DNA and nuclear ribosomal DNA restriction fragment length polymorphism (RFLP) profiles, each accession was identified as a species or hybrid combination. Five of the six possible hybrid combinations among the four species were detected. Rhododendron ponticum × R. smirnovii was represented by a single individual and R. caucasicum × R. smirnovii by one small group of hybrid plants. The combinations R. ponticum × R. ungernii and R. ungernii × R. smirnovii showed evidence of frequent backcrossing, while R. ponticum × R. caucasicum appeared unusual in that an intermediate hybrid type was abundant, whereas hybrids with phenotypes approaching either parent were rare. Possible explanations of this latter situation are discussed. The results suggest that natural hybridization among Rhododendron species is common and that ecological factors are important in maintaining integrity when species occur in sympatry.     

Genomic regions with a history of divergent selection affect fitness of hybrids between two butterfly species

Speciation is the process by which reproductively isolated lineages arise, and is one of the fundamental means by which the diversity of life increases. Whereas numerous studies have documented an association between ecological divergence and reproductive isolation, relatively little is known about the role of natural selection in genome divergence during the process of speciation. Here, we use genome-wide DNA sequences and Bayesian models to test the hypothesis that loci under divergent selection between two butterfly species (Lycaeides idas and L. melissa) also affect fitness in an admixed population. Locus-specific measures of genetic differentiation between L. idas and L. melissa and genomic introgression in hybrids varied across the genome. The most differentiated genetic regions were characterized by elevated L. idas ancestry in the admixed population, which occurs in L. idas-like habitat, consistent with the hypothesis that local adaptation contributes to speciation. Moreover, locus-specific measures of genetic differentiation (a metric of divergent selection) were positively associated with extreme genomic introgression (a metric of hybrid fitness). Interestingly, concordance of differentiation and introgression was only partial. We discuss multiple, complementary explanations for this partial concordance.     

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.     

Rejecting strictly allopatric speciation on a continental island: prolonged postdivergence gene flow between Taiwan (Leucodioptron taewanus,Passeriformes Timaliidae) and Chinese (L. canorum canorum) hwameis

Allopatry is conventionally considered the geographical mode of speciation for continental island organisms. However, strictly allopatric speciation models that assume the lack of postdivergence gene flow seem oversimplified given the recurrence of land bridges during glacial periods since the late Pliocene. Here, to evaluate whether a continental island endemic, the Taiwan hwamei (Leucodioptron taewanus, Passeriformes Timaliidae) speciated in strict allopatry, we used weighted‐regression‐based approximate Bayesian computation (ABC) to analyse the genetic polymorphism of 18 neutral nuclear loci (total length: 8500 bp) in Taiwan hwamei and its continental sister species, the Chinese hwamei (L. canorum canorum). The nonallopatry model was found to fit better with observed genetic polymorphism of the two hwamei species (posterior possibility = 0.82). We also recovered unambiguous signals of nontrivial bidirectional postdivergence gene flow (N_em » 1) between Chinese hwamei and Taiwan hwamei until 0.5 Ma. Divergence time was estimated to be 3.5 to 2 million years earlier than that estimated from mitochondrial cytochrome b sequences. Finally, using the inferred nonallopatry model to simulate genetic variation at 24 nuclear genes examined showed that the adiponectin receptor 1 gene may be under divergent adaptation. Our findings imply that the role of geographical barrier may be less prominent for the speciation of continental island endemics, and suggest a shift in speciation studies from simply correlating geographical barrier and genetic divergence to examining factors that facilitate and maintain divergence, e.g. differential selection and sexual selection, especially in the face of interpopulation gene flow.     

Evidence for sexual isolation as a prezygotic barrier to gene flow between morphologically divergent species of Rhagoletis fruit flies

1. Certain groups of fruit flies in the genus Rhagoletis (Diptera: Tephritidae) are exemplars for sympatric speciation via host plant shifting. Flies in these species groups are morphologically similar and overlap in their geographic ranges, yet attack different, non‐overlapping sets of host plants. Ecological adaptations related to differences in host choice and preference have been shown to be important prezygotic barriers to gene flow between these taxa, as Rhagoletis flies mate on or near the fruit of their respective host plants. Non‐host‐related assortative mating is generally absent or present at low levels between these sympatrically diverging fly populations. 2. However, some Rhagoletis taxa occasionally migrate to ‘non‐natal’ plants that are the primary hosts of other, morphologically differentiated fly species in the genus. These observations raise the question of whether sexual isolation may reduce courtship and copulation between morphologically divergent species of Rhagoletis flies, contributing to their prezygotic isolation along with host‐specific mating. 3. Using reciprocal multiple‐choice mating trials, we measured sexual isolation among nine species pairs of morphologically differentiated Rhagoletis flies. Complete sexual isolation was observed in eight of the nine comparisons, while partial sexual isolation was observed in the remaining comparison. 4. We conclude that sexual isolation can be an effective prezygotic barrier to gene flow contributing to substantial reproductive isolation between many morphologically distinct Rhagoletis species, even in the absence of differential host plant choice and host‐associated mating.     

Genetic,morphological, and niche variation in the widely hybridizing Rhus integrifolia-Rhus ovata species complex

Hybridization and introgression are common processes among numerous plant species that present both challenges and opportunities for studies of species delimitation, phylogenetics, taxonomy and adaptation. Rhus integrifolia and R. ovata are two ecologically important shrubs native to the southwestern USA and Mexico, and are known to hybridize frequently, but the morphological, genetic and ecological implications of hybridization in these species are poorly studied on a broad geographic scale. Analyses were conducted using leaf morphology, genetic variation of plastid and nuclear loci, and species distribution models for both species and their putative hybrid introgressants across 19 localities in California and Arizona, USA. These analyses revealed evidence for morphological and genetic distinction among localities comprising putative parental species, but a high degree of morpho-genetic intermediacy among localities with putative hybrids. Comparison of morphological and genetic population structure among localities revealed evidence for putative local adaptation or widespread phenotypic plasticity. Multiple regression models identified a weak but statistically significant negative association between leaf area and precipitation. Finally, species distribution modeling inferred northward range shifts over time, with both species predicted to occupy more coastal regions in the future, possibly increasing the frequency of hybridization among them. These findings underscore the importance of integrative assessment of multiple data sources in the study of hybridizing species and highlight the R. integrifolia-ovata complex as a powerful model for investigating the adaptive implications of hybridization.     

Differences in timing of migration and response to sexual signalling drive asymmetric hybridization across a migratory divide

Ecological traits and sexual signals may both contribute to the process of ecological speciation. Here we investigate the roles of an ecological trait, seasonal migratory behaviour and a sexual trait, song, in restricting or directing gene flow across a migratory divide in the Swainson's thrush (Catharus ustulatus). We show that short-distance migratory ecotypes wintering in Central America arrive earlier at the breeding grounds than long-distance migratory ecotypes wintering primarily in South America, providing the potential for some premating isolation. Playback experiments suggest that early- and late-arriving forms recognize each other as competitors, but that the early-arriving form responds more aggressively to a broader spectrum of stimuli. Genetic analysis suggests that hybridization occurs more often between males of the early-arriving ecotype and females of the late-arriving ecotype. Together our results suggest that differences in arrival times may reduce the temporal coincidence of mate choice, but asymmetry in response to heterotypic song may hinder complete divergence. These data provide further insight into the roles of ecological traits and sexual signals during the incipient stages of speciation.     

Extensive Introgression of Mitochondrial DNA Found between two Genetically Divergent Forms of Threespine Stickleback, Gasterosteus Aculeatus, around Japan

The genetic structure and phylogenetic relationships of Gasterosteus aculeatus around Japan were investigated for 15 populations of the Japan Sea and Pacific Ocean forms previously detected by allozyme analyses. PCR-RFLP analysis of mtDNA ND5/6 gene indicated two major divergent clades: (1) Alaskan population and (2) Russian and Japanese populations. Japanese populations were clearly subdivided into G. aculeatus populations and hariyo (G. aculeatus leiurus). However, the relationships among Japanese G. aculeatus populations did not support the genetic divergence between Japan Sea and Pacific Ocean forms shown by allozyme analyses. This suggests that mtDNA genetype similarity between the two forms could have resulted from complete replacement of mtDNA from one form to the other. MtDNA introgression is possible in two direction: from Japan Sea form to Pacific Ocean form and the reverse. The haplotype composition of each population observed in the present study suggests that the second direction is more likely. The migrant number of females per generation (Nm) between populations of the two forms was remarkablely large (9.083–infinite), suggesting that gene flow is taking place between the populations at the present time or might have occurred until recent years.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.