A hybrid dilemma: a molecular investigation of South African bontebok (Damaliscus pygargus pygargus) and blesbok (Damaliscus pygargus phillipsi)

Human-mediated hybridization poses a serious threat to the genetic integrity of threatened inter-fertile species or subspecies, such as the bontebok (Damaliscus pygargus pygargus) and blesbok (Damaliscus pygargus phillipsi) in South Africa. The identification of pure and admixed populations is key to sound biodiversity conservation management and practices. Where identification based on morphometric measurements are complex, molecular analysis can provide a tool to identify hybrids or introgression in populations. This study provides the first molecular analyses on pure bontebok and blesbok as well as putative hybrid populations and samples of unknown purity. Analysis of a panel of 13 microsatellite markers by applying a Bayesian model-based clustering programme provided evidence of clear partitioning between the two subspecies. Hybridization was detected in 33 % (40 of 121) of the samples with unknown purity. Currently, bontebok is listed as protected under the Convention on International Trade in Endangered Species of Wild Fauna and Flora and is listed as endangered under the United States of America’s Endangered Species Act. Our research found that the microsatellite markers used in this study could accurately identify pure and hybrid animals and can be used as a tool to conserve the genetic integrity of both subspecies.     

Patterns of hybridization among cutthroat trout and rainbow trout in northern Rocky Mountain streams

Introgressive hybridization between native and introduced species is a growing conservation concern. For native cutthroat trout and introduced rainbow trout in western North America, this process is thought to lead to the formation of hybrid swarms and the loss of monophyletic evolutionary lineages. Previous studies of this phenomenon, however, indicated that hybrid swarms were rare except when native and introduced forms of cutthroat trout co‐occurred. We used a panel of 86 diagnostic, single nucleotide polymorphisms to evaluate the genetic composition of 3865 fish captured in 188 locations on 129 streams distributed across western Montana and northern Idaho. Although introgression was common and only 37% of the sites were occupied solely by parental westslope cutthroat trout, levels of hybridization were generally low. Of the 188 sites sampled, 73% contained ≤5% rainbow trout alleles and 58% had ≤1% rainbow trout alleles. Overall, 72% of specimens were nonadmixed westslope cutthroat trout, and an additional 3.5% were nonadmixed rainbow trout. Samples from seven sites met our criteria for hybrid swarms, that is, an absence of nonadmixed individuals and a random distribution of alleles within the sample; most (6/7) were associated with introgression by Yellowstone cutthroat trout. In streams with multiple sites, upstream locations exhibited less introgression than downstream locations. We conclude that although the widespread introduction of nonnative trout within the historical range of westslope cutthroat trout has increased the incidence of introgression, sites containing nonadmixed populations of this taxon are common and broadly distributed.     

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co‐occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first‐generation hybrids. Later‐generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual‐based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.     

Evaluating the ability of Bayesian clustering methods to detect hybridization and introgression using an empirical red wolf data set

Bayesian clustering methods have emerged as a popular tool for assessing hybridization using genetic markers. Simulation studies have shown these methods perform well under certain conditions; however, these methods have not been evaluated using empirical data sets with individuals of known ancestry. We evaluated the performance of two clustering programs, baps and structure , with genetic data from a reintroduced red wolf (Canis rufus) population in North Carolina, USA. Red wolves hybridize with coyotes (C. latrans), and a single hybridization event resulted in introgression of coyote genes into the red wolf population. A detailed pedigree has been reconstructed for the wild red wolf population that includes individuals of 50–100% red wolf ancestry, providing an ideal case study for evaluating the ability of these methods to estimate admixture. Using 17 microsatellite loci, we tested the programs using different training set compositions and varying numbers of loci. structure was more likely than baps to detect an admixed genotype and correctly estimate an individual's true ancestry composition. However, structure was more likely to misclassify a pure individual as a hybrid. Both programs were outperformed by a maximum‐likelihood‐based test designed specifically for this system, which never misclassified a hybrid (50–75% red wolf) as a red wolf or vice versa. Training set composition and the number of loci both had an impact on accuracy but their relative importance varied depending on the program. Our findings demonstrate the importance of evaluating methods used for detecting admixture in the context of endangered species management.     

Genetic and morphological evidence of a geographically widespread hybrid zone between two crocodile species,Crocodylus acutus and Crocodylus moreletii

Hybrid zones represent natural laboratories to study gene flow, divergence and the nature of species boundaries between closely related taxa. We evaluated the level and extent of hybridization between Crocodylus moreletii and Crocodylus acutus using genetic and morphological data on 300 crocodiles from 65 localities. To our knowledge, this is the first genetic study that includes the entire historic range and sympatric zone of the two species. Contrary to expectations, Bayesian admixture proportions and maximum‐likelihood estimates of hybrid indexes revealed that most sampled crocodiles were admixed and that the hybrid zone is geographically extensive, extending well beyond their historical region of sympatry. We identified a few geographically isolated, nonadmixed populations of both parental species. Hybrids do not appear to be F₁s or recent backcrosses, but rather are more likely later‐generation hybrids, suggesting that hybridization has been going on for several to many generations and is mostly the result of natural processes. Crocodylus moreletii is not the sister species of C. acutus, suggesting that the hybrid zone formed from secondary contact rather than primary divergence. Nonadmixed individuals from the two species were distinguishable based on morphological characters, whereas hybrids had a complex mosaic of morphological characters that hinders identification in the wild. Very few nonadmixed C. acutus and C. moreletii populations exist in the wild. Consequently, the last nonadmixed C. moreletii populations have become critically endangered. Indeed, not only the parental species but also the naturally occurring hybrids should be considered for their potential conservation value.     

Death by sex in an Australian icon: a continent‐wide survey reveals extensive hybridization between dingoes and domestic dogs

Hybridization between domesticated animals and their wild counterparts can disrupt adaptive gene combinations, reduce genetic diversity, extinguish wild populations and change ecosystem function. The dingo is a free‐ranging dog that is an iconic apex predator and distributed throughout most of mainland Australia. Dingoes readily hybridize with domestic dogs, and in many Australian jurisdictions, distinct management strategies are dictated by hybrid status. Yet, the magnitude and spatial extent of domestic dog–dingo hybridization is poorly characterized. To address this, we performed a continent‐wide analysis of hybridization throughout Australia based on 24 locus microsatellite DNA genotypes from 3637 free‐ranging dogs. Although 46% of all free‐ranging dogs were classified as pure dingoes, all regions exhibited some hybridization, and the magnitude varied substantially. The southeast of Australia was highly admixed, with 99% of animals being hybrids or feral domestic dogs, whereas only 13% of the animals from remote central Australia were hybrids. Almost all free‐ranging dogs had some dingo ancestry, indicating that domestic dogs could have poor survivorship in nonurban Australian environments. Overall, wild pure dingoes remain the dominant predator over most of Australia, but the speed and extent to which hybridization has occurred in the approximately 220 years since the first introduction of domestic dogs indicate that the process may soon threaten the persistence of pure dingoes.     

Differentiation at mitochondrial and nuclear loci between the blesbok (Damaliscus pygargus phillipsi) and bontebok (D. p. pygargus): implications for conservation strategy

Damaliscus pygargus is an endemic species in South Africa belonging to the contemporary antelope tribe Alcelaphine. The species has been subdivided into two subspecies based upon phenotypic differences and historical geographic isolation. We found D. p. pygargus (bontebok) and D. p. phillipsi (blesbok) to be significantly differentiated based on microsatellite and mtDNA control region markers. Bontebok genetic diversity was depauperate with just one mtDNA haplotype, not found in blesbok populations, and reduced heterozygosity and polymorphism at microsatellite DNA loci compared to the blesbok. Erosion of molecular genetic variation in bontebok may have been due to anthropogenic impact causing isolation and reduced population size. Our data based on 34 bontebok and 42 blesbok indicate that the classification of alpha taxonomy should be reconsidered in this genus (in support of earlier similar suggestions), and that careful management should seek to avoid hybridization and retain the remaining diversity in the bontebok.     

Fitness of natural willow hybrids in a pioneer mosaic hybrid zone

Hybrid fitness is an important parameter to predict the evolutionary consequences of a hybridization event and to characterize hybrid zones. We studied fitness parameters of F₁ and later‐generation hybrids between the lowland species Salix purpurea and the alpine S. helvetica that have recently emerged during colonization of an alpine glacier forefield. Fruit production (number of capsules per catkin and fruit set) did not differ between hybrids and parents, but the number of seeds per capsule of F₁ hybrids was slightly lower than that of later‐generation hybrids and of the parents. Germination rates and seedling growth were tested on three substrates (pH 4.5, 7.0, and 8.0). Germination rates of seeds collected from F₁ hybrids were lower on acid and neutral substrates, but equal at pH 8.0 compared to all other groups, while the seeds from later‐generation hybrids performed as well as the parents on all three substrates. In seedling growth, the colonizer S. purpurea performed better than all other taxa on all three substrates, while hybrids resembled the subalpine species S. helvetica. Results suggest that endogenous selection acts against F₁ hybrids, but favors fitter genotypes in later‐generation hybrids. Exogenous selection via soil pH appears to be weak during seedling establishment. The pioneer vegetation on the glacier forefield may offer sufficient niche space for hybrid seedlings. Owing to the relatively high fitness of the hybrids and the scattered distribution of hybrids and parental individuals on the glacier forefield, this hybrid zone can be assigned to a mosaic model, probably facilitating gene flow and introgression between the parental species. As establishment of the hybrid zone appears to be linked to a colonization process, we propose to call it a pioneer mosaic hybrid zone.     

Y‐chromosome evidence supports widespread signatures of three‐species Canis hybridization in eastern North America

There has been considerable discussion on the origin of the red wolf and eastern wolf and their evolution independent of the gray wolf. We analyzed mitochondrial DNA (mtDNA) and a Y‐chromosome intron sequence in combination with Y‐chromosome microsatellites from wolves and coyotes within the range of extensive wolf–coyote hybridization, that is, eastern North America. The detection of divergent Y‐chromosome haplotypes in the historic range of the eastern wolf is concordant with earlier mtDNA findings, and the absence of these haplotypes in western coyotes supports the existence of the North American evolved eastern wolf (Canis lycaon). Having haplotypes observed exclusively in eastern North America as a result of insufficient sampling in the historic range of the coyote or that these lineages subsequently went extinct in western geographies is unlikely given that eastern‐specific mtDNA and Y‐chromosome haplotypes represent lineages divergent from those observed in extant western coyotes. By combining Y‐chromosome and mtDNA distributional patterns, we identified hybrid genomes of eastern wolf, coyote, gray wolf, and potentially dog origin in Canis populations of central and eastern North America. The natural contemporary eastern Canis populations represent an important example of widespread introgression resulting in hybrid genomes across the original C. lycaon range that appears to be facilitated by the eastern wolf acting as a conduit for hybridization. Applying conventional taxonomic nomenclature and species‐based conservation initiatives, particularly in human‐modified landscapes, may be counterproductive to the effective management of these hybrids and fails to consider their evolutionary potential.     

Interploidy hybridization in sympatric zones: the formation of Epidendrum fulgens × E. puniceoluteum hybrids (Epidendroideae,Orchidaceae)

Interspecific hybridization is a primary cause of extensive morphological and chromosomal variation and plays an important role in plant species diversification. However, the role of interploidal hybridization in the formation of hybrid swarms is less clear. Epidendrum encompasses wide variation in chromosome number and lacks strong premating barriers, making the genus a good model for clarifying the role of chromosomes in postzygotic barriers in interploidal hybrids. In this sense, hybrids from the interploidal sympatric zone between E. fulgens (2n = 2x = 24) and E. puniceoluteum (2n = 4x = 56) were analyzed using cytogenetic techniques to elucidate the formation and establishment of interploidal hybrids. Hybrids were not a uniform group: two chromosome numbers were observed, with the variation being a consequence of severe hybrid meiotic abnormalities and backcrossing with E. puniceoluteum. The hybrids were triploids (2n = 3x = 38 and 40) and despite the occurrence of enormous meiotic problems associated with triploidy, the hybrids were able to backcross, producing successful hybrid individuals with broad ecological distributions. In spite of the nonpolyploidization of the hybrid, its formation is a long‐term evolutionary process rather than a product of a recent disturbance, and considering other sympatric zones in Epidendrum, these events could be recurrent.     

Heterosis in hybrids within and between yeast species

The performance of hybrids relative to their parents is an important factor in speciation research. We measured the growth of 46 Saccharomyces yeast F1 interspecific and intraspecific hybrids, relative to the growth of each of their parents, in pairwise competition assays. We found that the growth of a hybrid relative to the average of its parents, a measure of mid‐parent heterosis, correlated with the difference in parental growth relative to their hybrid, a measure of phenotypic divergence, which is consistent with simple complementation of low fitness alleles in one parent by high fitness alleles in the other. Interspecific hybrids showed stronger heterosis than intraspecific hybrids. To manipulate parental phenotypic divergence independently of genotype, we also measured the competitive growth of a single interspecific hybrid relative to its parents in 12 different environments. In these assays, we not only identified a strong relationship between parental phenotypic divergence and mid‐parent heterosis as before, but, more tentatively, a weak relationship between phenotypic divergence and best‐parent heterosis, suggesting that complementation of deleterious mutations was not the sole cause of interspecific heterosis. Our results show that mating between different species can be beneficial, and demonstrate that competition assays between parents and offspring are a useful way to study the evolutionary consequences of hybridization.     

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.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Morphological variation of genetically confirmed Alouatta Pigra × A. palliata hybrids from a natural hybrid zone in Tabasco,Mexico

While hybridization has been reported for a large number of primate taxa, there is a general lack of data on hybrid morphology for wild individuals with known genetic ancestry. A confirmed hybrid zone for the closely related Neotropical primates Alouatta palliata and A. pigra has provided a unique opportunity to study primate hybrid morphological variation. Here we used molecular evidence based on mitochondrial, Y‐chromosome, and autosomal data to assess hybrid ancestry. We conducted univariate and multivariate statistical comparisons of morphometric data collected from individuals both outside and within the hybrid zone in Tabasco, Mexico. Our results show that of all the hybrids detected (N = 128), only 12% of them were approximately genetically intermediate, and none of them were first generation hybrids. Univariate pairwise comparisons among parental individuals, multigenerational backcrossed hybrids, and intermediate hybrids showed that overall, multigenerational backcrossed hybrids resemble the parental species with which they share most of their alleles. Conversely, intermediates were highly variable. Similarly, principal component analysis depicts an overlap between the parental species and their backcrosses when considering overall morphological differences. Finally, discriminant function analysis of the morphological variables was overall unreliable for classifying individuals into their assigned genotypic classes. Taken together, our results suggest that primate natural hybridization studies should incorporate molecular methods for determining ancestry, because morphology may not always be a reliable indicator of hybrid status. Hybrid zones could comprise a large number of multigenerational backcrossed hybrids that are indistinguishable from the parental species. The implications for studying hybridization in the primate fossil record are discussed. Am J Phys Anthropol, 2013. © 2012 Wiley Periodicals, Inc.     

Hybridization between two sympatrically breeding species of fur seal at Iles Crozet revealed by genetic analysis

Introgressive hybridization occurs when closely related taxa overlap in distribution and is often associated with historically isolated populations coming into contact as a result of anthropogenic disturbance. There is evolutionary and conservation interest in detecting hybridization to determine its implications on future species composition, especially for threatened and recovering taxa such as subantarctic (Arctocephalus tropicalis) and Antarctic (A. gazella) fur seals, which were driven to the brink of extinction by human exploitation. Hybridization between these species has been reported at two locations and they breed sympatrically at a third site, Iles Crozet. While hybrid individuals have previously been identified based on phenotype, individuals can be difficult to classify based on these characteristics alone. Genotypic hybrid identification has been successful in several species, including fur seals. In this study we conducted an assignment test using microsatellite data to identify hybrids and to measure the frequency of hybridization at Iles Crozet. Samples were collected from 372 individuals and screened with 6 polymorphic microsatellite markers. MtDNA genotypes were also determined for individuals identified as hybrids or backcrosses based on microsatellite genotype. Phenotype, microsatellite and mtDNA genotype were then compared in order to identify hybrids. The results indicate that 1% of the population have hybrid genotypes and at a minimum, 2.4% of the population are backcrossed to parental species. We found that the two species are genetically distinct from one another and given the low rate of hybridization it is unlikely that they will fuse. These results suggest that there is a mechanism for species recognition that acts as a barrier to hybridization. It therefore seems unlikely that fur seals are threatened by significant introgression. Further investigation of fur seal mating systems would provide valuable insight into the mechanisms governing hybridization and species recognition in mate choice.     

Widespread hybridization in the introduced hog deer population of Victoria,Australia, and its implications for conservation

In Australia, many species have been introduced that have since undergone drastic declines in their native range. One species of note is the hog deer (Axis porcinus) which was introduced in the 1860s to Victoria, Australia, and has since become endangered in its native range throughout South‐East Asia. There is increased interest in using non‐native populations as a source for genetic rescue; however, considerations need to be made of the genetic suitability of the non‐native population. Three mitochondrial markers and two nuclear markers were sequenced to assess the genetic variation of the Victorian population of hog deer, which identified that the Victorian population has hybrid origins with the closely related chital (Axis axis), a species that is no longer present in the wild in Victoria. In addition, the mitochondrial D‐loop region within the Victorian hog deer is monomorphic, demonstrating that mitochondrial genetic diversity is very low within this population. This study is the first to report of long‐term persistence of hog deer and chital hybrids in a wild setting, and the continual survival of this population suggests that hybrids of these two species are fertile. Despite the newly discovered hybrid status in Victorian hog deer, this population may still be beneficial for future translocations within the native range. However, more in‐depth analysis of genetic diversity within the Victorian hog deer population and investigation of hybridization rates within the native range are necessary before translocations are attempted.     

Hybridization in closely related Rhododendron species: half of all species‐differentiating markers experience serious transmission ratio distortion

An increasing number of studies of hybridization in recent years have revealed that complete reproductive isolation between species is frequently not finalized in more or less closely related organisms. Most of these species do, however, seem to retain their phenotypical characteristics despite the implication of gene flow, highlighting the remaining gap in our knowledge of how much of an organism's genome is permeable to gene flow, and which factors promote or prevent hybridization. We used AFLP markers to investigate the genetic composition of three populations involving two interfertile Rhododendron species: two sympatric populations, of which only one contained hybrids, and a further hybrid‐dominated population. No fixed differences between the species were found, and only 5.8% of the markers showed some degree of species differentiation. Additionally, 45.5% of highly species‐differentiating markers experienced significant transmission distortion in the hybrids, which was most pronounced in F1 hybrids, suggesting that factors conveying incompatibilities are still segregating within the species. Furthermore, the two hybrid populations showed stark contrasting composition of hybrids; one was an asymmetrically backcrossing hybrid swarm, while in the other, backcrosses were absent, thus preventing gene flow.     

Genomic consequences of a recent three‐way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human‐mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make‐up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double‐digested restriction site‐associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple‐way admixture with population genomic data.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Genetic differentiation and inferred dynamics of a hybrid zone between Northern Spotted Owls (Strix occidentalis caurina) and California Spotted Owls (S. o. occidentalis) in northern California

Genetic differentiation among Spotted Owl (Strix occidentalis) subspecies has been established in prior studies. These investigations also provided evidence for introgression and hybridization among taxa but were limited by a lack of samples from geographic regions where subspecies came into close contact. We analyzed new sets of samples from Northern Spotted Owls (NSO: S. o. caurina) and California Spotted Owls (CSO: S. o. occidentalis) in northern California using mitochondrial DNA sequences (mtDNA) and 10 nuclear microsatellite loci to obtain a clearer depiction of genetic differentiation and hybridization in the region. Our analyses revealed that a NSO population close to the northern edge of the CSO range in northern California (the NSO Contact Zone population) is highly differentiated relative to other NSO populations throughout the remainder of their range. Phylogenetic analyses identified a unique lineage of mtDNA in the NSO Contact Zone, and Bayesian clustering analyses of the microsatellite data identified the Contact Zone as a third distinct population that is differentiated from CSO and NSO found in the remainder of the subspecies' range. Hybridization between NSO and CSO was readily detected in the NSO Contact Zone, with over 50% of individuals showing evidence of hybrid ancestry. Hybridization was also identified among 14% of CSO samples, which were dispersed across the subspecies' range in the Sierra Nevada Mountains. The asymmetry of hybridization suggested that the hybrid zone may be dynamic and moving. Although evidence of hybridization existed, we identified no F1 generation hybrid individuals. We instead found evidence for F2 or backcrossed individuals among our samples. The absence of F1 hybrids may indicate that (1) our 10 microsatellites were unable to distinguish hybrid types, (2) primary interactions between subspecies are occurring elsewhere on the landscape, or (3) dispersal between the subspecies' ranges is reduced relative to historical levels, potentially as a consequence of recent regional fires.