Quantitative evaluation of hybridization and the impact on biodiversity conservation

Anthropogenic hybridization is an increasing conservation threat worldwide. In South Africa, recent hybridization is threatening numerous ungulate taxa. For example, the genetic integrity of the near‐threatened bontebok (Damaliscus pygargus pygargus) is threatened by hybridization with the more common blesbok (D. p. phillipsi). Identifying nonadmixed parental and admixed individuals is challenging based on the morphological traits alone; however, molecular analyses may allow for accurate detection. Once hybrids are identified, population simulation software may assist in determining the optimal conservation management strategy, although quantitative evaluation of hybrid management is rarely performed. In this study, our objectives were to describe species‐wide and localized rates of hybridization in nearly 3,000 individuals based on 12 microsatellite loci, quantify the accuracy of hybrid assignment software (STRUCTURE and NEWHYBRIDS), and determine an optimal threshold of bontebok ancestry for management purposes. According to multiple methods, we identified 2,051 bontebok, 657 hybrids, and 29 blesbok. More than two‐thirds of locations contained at least some hybrid individuals, with populations varying in the degree of introgression. HYBRIDLAB was used to simulate four generations of coexistence between bontebok and blesbok, and to optimize a threshold of ancestry, where most hybrids will be detected and removed, and the fewest nonadmixed bontebok individuals misclassified as hybrids. Overall, a threshold Q‐value (admixture coefficient) of 0.90 would remove 94% of hybrid animals, while a threshold of 0.95 would remove 98% of hybrid animals but also 8% of nonadmixed bontebok. To this end, a threshold of 0.90 was identified as optimal and has since been implemented in formal policy by a provincial nature conservation agency. Due to widespread hybridization, effective conservation plans should be established and enforced to conserve native populations that are genetically unique.     

Discrimination of hybrid classes using cross-species amplification of microsatellite loci: methodological challenges and solutions in Daphnia

Microsatellite markers are important tools in population, conservation and forensic studies and are frequently used for species delineation, the detection of hybridization and introgression. Therefore, marker sets that amplify variable DNA regions in two species are required; however, cross-species amplification is often difficult, as genotyping errors such as null alleles may occur. To estimate the level of potential misidentifications based on genotyping errors, we compared the occurrence of parental alleles in laboratory and natural Daphnia hybrids (Daphnia longispina group). We tested a set of 12 microsatellite loci with regard to their suitability for unambiguous species and hybrid class identification using F(1) hybrids bred in the laboratory. Further, a large set of 44 natural populations of D. cucullata, D. galeata and D. longispina (1715 individuals) as well as their interspecific hybrids were genotyped to validate the discriminatory power of different marker combinations. Species delineation using microsatellite multilocus genotypes produced reliable results for all three studied species using assignment tests. Daphnia galeata × cucullata hybrid detection was limited due to three loci exhibiting D. cucullata-specific null alleles, which most likely are caused by differences in primer-binding sites of parental species. Overall, discriminatory power in hybrid detection was improved when a subset of markers was identified that amplifies equally well in both species.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

The effects of locus number,genetic divergence,and genotyping error on the utility of dominant markers for hybrid identification

In surveys of hybrid zones, dominant genetic markers are often used to identify individuals of hybrid origin and assign these individuals to one of several potential hybrid classes. Quantitative analyses that address the statistical power of dominant markers in such inference are scarce. In this study, dominant genotype data were simulated to evaluate the effects of, first, the number of loci analyzed, second, the magnitude of differentiation between the markers scored in the groups that are hybridizing, and third, the level of genotyping error associated with the data when assigning individuals to various parental and hybrid categories. The overall performance of the assignment methods was relatively modest at the lowest level of divergence examined (F_st ~ 0.4), but improved substantially at higher levels of differentiation (F_st ~ 0.67 or 0.8). The effect of genotyping error was dependent on the level of divergence between parental taxa, with larger divergences tempering the effects of genotyping error. These results highlight the importance of considering the effects of each of the variables when assigning individuals to various parental and hybrid categories, and can help guide decisions regarding the number of loci employed in future hybridization studies to achieve the power and level of resolution desired.     

Chromosomal basis of viability differences in Tigriopus californicus interpopulation hybrids

Crosses between populations of Tigriopus californicus result in backcross and F2 hybrid breakdown for a variety of fitness related measures. The magnitude of this hybrid breakdown is correlated with evolutionary divergence. We assessed the chromosomal basis of viability differences in nonrecombinant backcross hybrids using markers mapped to individual chromosomes. To assess effects of evolutionary divergence we crossed one population to three different populations: two distantly related (approximately 18% mitochondrial COI sequence divergence) and one closely related (approximately 1% mitochondrial COI sequence divergence). We found that all three interpopulation crosses resulted in significant deviations from expected Mendelian ratios at a majority of the loci studied. In all but one case, deviations were due to a deficit of parental homozygotes. This pattern implies that populations of T. californicus carry a significant genetic load, and that a combination of beneficial dominance and deleterious homozygote-heterozygote interactions significantly affects hybrid viability. Pairwise tests of linkage disequilibrium detected relatively few significant interactions. For the two divergent crosses, effects of individual chromosomes were highly concordant. These two crosses also showed higher heterozygote excess in females than males across the vast majority of chromosomes.     

Epistatic and cytonuclear interactions govern outbreeding depression in the autotetraploid Campanulastrum americanum

The consequences of combining divergent genomes among populations of a diploid species often involve F1 hybrid vigor followed by hybrid breakdown in later recombinant generations. As many as 70% of plant species are thought to have polyploid origins; yet little is known about the genetic architecture of divergence in polyploids and how it may differ from diploid species. We investigated the genetic architecture of population divergence using controlled crosses among five populations of the autotetraploid herb, Campanulastrum americanum. Plants were reciprocally hybridized to produce F1, F2, and F1-backcross generations that were grown with parental types in a greenhouse and measured for performance. In contrast to diploid expectations, most F1 hybrids lacked heterosis and instead showed strong outbreeding depression for early life traits. Recombinant hybrid generations often showed a recovery of performance to levels approximating, or at times even exceeding, the parental values. This pattern was also evident for an index of cumulative fitness. Analyses of line means indicated nonadditive gene action, especially forms of digenic epistasis, often influenced hybrid performance. However, standard diploid genetic models were not adequate for describing the underlying genetic architecture in a number of cases. Differences between reciprocal hybrids indicated that cytoplasmic and/or cytonuclear interactions also contributed to divergence. An enhanced role of epistasis in population differentiation may be the norm in polyploids, which have more gene copies. This study, the first of its kind on a natural autotetraploid, suggests that gene duplication may cause polyploid populations to diverge in a fundamentally different way than diploids.     

A spatially explicit individual-based model of reinforcement in hybrid zones

Abstract An individual-based model consisting of two dioecious populations in a two-dimensional environmental grid was constructed. Each population began with, and never exceeded, 1000 individuals; extinction was allowed. Genomes consisting of 30 biallelic loci for male sexual advertisement call, female mate preference, and population origin were constructed, and lineages of each individual in the starting populations were followed for 2000 generations. Type and level of hybrid disadvantage, initial population distribution, patchiness of environmental resources, and level of mate choice were varied. Persistence of bimodal hybrid zones was nonexistent at low levels of hybrid disadvantage and universal at high levels of hybrid disadvantage, with a narrow threshold in which persistence was unpredictable. Persistence occurred at lower levels of hybrid disadvantage when populations were initially parapatric rather than sympatric, and environments were patchy rather than homogeneous. Increased divergence in mating systems occurred when hybrid disadvantage was high, hybrids were infertile, populations were initially parapatric, and increased female choice was allowed. Mating system divergence was much higher in interacting populations compared with noninteracting populations, indicating that reinforcement caused most of the observed divergence. When hybrids were infertile, reinforcement contributed to speciation, because under hybrid infertility the probability of persistence at low levels of hybrid disadvantage was positively related to mate choice. The results agree with previous one-dimensional spatial models in finding that population persistence is more likely in parapatric and patchy population distributions. In addition, the results show that hybrid infertility may facilitate the process of reinforcement and speciation.     

SNP genotyping for detecting the ‘rare allele phenomenon’ in hybrid zones

Hybrid zones are regions where genetically distinct populations meet, mate and produce offspring. In such zones, genetically less compatible gene combinations are usually generated, resulting in reduced fitness, and hybrid zones are often maintained because of continuous removal of unfit genotypes, balanced by gene flow into the zone from the parental populations (and are then referred to as ‘tension zones’). Tension zones often display unexpectedly high frequencies of gene variants that are rare outside the zone. Previous work has shown that this ‘rare allele phenomenon’ is not the result of intragenic recombination or increased mutation rates. Further understanding of the population genetics of the phenomenon requires an approach in which both the numbers of individuals and the numbers of loci is increased. Here, we report an approach using a combination of Illumina next‐generation sequencing and mass spectrophotometer genotyping to identify markers that may be used for genome‐wide investigations of the rare allele phenomenon. We test this approach on a hybrid zone in the land snail Albinaria hippolyti from Greece.     

Hybrid bridges to gene flow: a case study in milkweeds (Asclepias)

Natural hybridization occurs throughout areas of sympatry for the North American milkweeds Asclepias exaltata and A. syriaca (Asclepiadaceae), even though the formation of F1 hybrid seed is a rare event. For introgressive hybridization to proceed, F1 and advanced hybrids must be released from reproductive barriers and successfully mate with one or both parental species. I investigated the mating system of natural hybrids between A. exaltata and A. syriaca in three populations in Shenandoah National Park, Virginia. Allozyme data and a maximum-likelihood procedure were used to estimate the frequency of six genotypic classes (parentals, F1, F2, and backcrosses) of the hybridizing populations, the pollinia received by hybrid plants, and the paternal parents of seeds produced by hybrids. F1 hybrids, backcross A. syriaca, and parental A. syriaca individuals were common in three hybrid populations. Even though self-pollinations and interhybrid pollinations were common, F2 seed production and the occurrence of F2 individuals were rare in hybrid populations. Hybrid plants received more pollen from A. syriaca than A. exaltata, which resulted in the production of more backcross-A. syriaca seed than backcross-A. exaltata seed. Asclepias exaltata was rare in the hybrid populations, but A. exaltata pollinia were received by hybrids and this species sired between 15% and 36% of the seeds produced on hybrids. The potential for introgression with A. exaltata populations is lower because this species is unsuccessful as the maternal parent in interspecific and backcross hand-pollinations. The asymetry of hybridization with A. syriaca as the maternal parent is further supported by the incorporation of maternally inherited chloroplast DNA markers in hybrids. Hybrid milkweeds frequently backcross with both parental species and may be released from the reproductive barriers that limit the formation of F1 hybrids in natural populations. The direction of interspecific gene flow and introgression in milkweeds is influenced by the reproductive biology of hybrids, the constituency of the surrounding population, and failure of some crosses to produce seeds. Finally, introgressive hybridization remains an important evolutionary force even when the initial formation of F1 hybrids in natural populations is rare.     

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.     

Classifying genealogical origins in hybrid populations using dominant markers

In hybrid studies, potential for error is high when classifying genealogical origins of individuals (e.g., parental, F1, F2) based on their genotypic arrays. For codominant markers, previous researchers have considered the probability of misclassification by genotypic inspection and proposed alternative maximum-likelihood approaches to estimating genealogical class frequencies. Recently developed dominant marker systems may significantly increase the number of diagnostic loci available for hybrid studies. I examine probabilities of classification error based on the number of dominant loci. As in earlier studies, I assume that only parental and first- and second-generation hybrid crosses between two taxa potentially exist. Thirteen loci with dominant expression from each parental taxon (i.e., 26 total loci) are needed to reduce classification error below 5% for F2 individuals, compared to 13 codominant loci for the same error rate. Use of loci in similar numbers from both taxa most efficiently increases power to characterize all genealogical classes. In contrast, classification of backcrosses to one parental taxon is wholly dependent on loci from the other taxon. Use of dominant diagnostic markers may increase the power and expand the use of maximum-likelihood methods for evaluating hybrid mixtures.     

Clonal diversity,clonal persistence and rapid taxon replacement in natural populations of species and hybrids of the Daphnia longispina complex

Hybridization is common among cyclical parthenogens, especially in zooplankton species assemblages of the genus Daphnia. To explore hybridization dynamics and the extent of clonal diversity in the Daphnia longispina complex, we analysed population structure in eight permanent lakes. Based on 15 microsatellite loci, three major taxonomic units emerged: two species, D. galeata and D. longispina and their F1 hybrids, supported by factorial correspondence analysis and two Bayesian methods. At the same time, the detection of backcross classes differed between methods. Mean clonal diversity was lowest in the F1 hybrids, as expected from the high rate of asexual reproduction. Within taxa, replicated genotypes were of clonal origin, but clonal lineages persisted in subsequent years in only one of three resampled lakes. In another lake, the taxon composition changed from being dominated by hybrids to complete dominance by one parental taxon. Such a year‐to‐year taxon replacement has not been reported for the D. longispina complex before. Our data on this hybrid complex illustrate that high‐resolution genotyping is essential for the understanding of ecological and evolutionary outcomes of hybridization in partially clonal taxa.     

Effectiveness of selective genotyping for detection of quantitative trait loci: an analysis of grain and malt quality traits in three barley populations.

Marker genotype data and grain and malt quality phenotype data from three barley (Hordeum vulgare L.) mapping populations were used to investigate the feasibility of selective genotyping for detection of quantitative trait loci (QTLs). With selective genotyping, only individuals with high and low phenotypic values for the trait of interest are genotyped. Here, genotyping of 10 to 70% of each population (i.e., 5 to 35% in each tail of the phenotypic distribution) was considered. Genomic positions detected by selective genotyping were compared to QTL position estimates from interval mapping analysis using marker genotype data from the entire population. Selective genotyping reliably detected almost all of the mapped QTLs, often with only 10% of the population genotyped. Selective genotyping also detected spurious QTLs in regions of the genome where no significant QTL had been mapped. Even with additional genotyping to verify putative QTLs, the total genotyping effort for detection of QTLs for a single trait by selective genotyping was usually less than 30% of that required for conventional interval mapping. Simultaneous investigation of two or more traits by selective genotyping would require additional genotyping effort, but could still be worthwhile.     

A 34K SNP genotyping array for Populus trichocarpa: Design,application to the study of natural populations and transferability to other Populus species

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost‐effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre‐ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.     

GENETIC ARCHITECTURE AND POSTZYGOTIC REPRODUCTIVE ISOLATION: EVOLUTION OF BATESON–DOBZHANSKY–MULLER INCOMPATIBILITIES IN A POLYGENIC MODEL

The Bateson–Dobzhansky–Muller model predicts that postzygotic isolation evolves due to the accumulation of incompatible epistatic interactions, but few studies have quantified the relationship between genetic architecture and patterns of reproductive divergence. We examined how the direction and magnitude of epistatic interactions in a polygenic trait under stabilizing selection influenced the evolution of hybrid incompatibilities. We found that populations evolving independently under stabilizing selection experienced suites of compensatory allelic changes that resulted in genetic divergence between populations despite the maintenance of a stable, high‐fitness phenotype. A small number of loci were then incompatible with multiple alleles in the genetic background of the hybrid and the identity of these incompatibility loci changed over the evolution of the populations. For F₁ hybrids, reduced fitness evolved in a window of intermediate strengths of epistatic interactions, but F₂ and backcross hybrids evolved reduced fitness across weak and moderate strengths of epistasis due to segregation variance. Strong epistatic interactions constrained the allelic divergence of parental populations and prevented the development of reproductive isolation. Because many traits with varying genetic architectures must be under stabilizing selection, our results indicate that polygenetic drift is a plausible hypothesis for the evolution of postzygotic reproductive isolation.     

不同作图群体对显隐性分子标记位点的作图效率

根据位点组合和位点的有效性,发展了一种对使用３种不同的作用图群体作图显隐性分子标记的作图效率评价方法,应用该方法所评价的结果是,双单倍体（ＤＨ）群体的作图效率最高,自交群体（Ｆ２群体）与回交群体的作图效率相同,因此使用双单倍体群体作图不仅所用费用低,而且作图速度快,但只有在极少数植物中能获得双单倍体群体,对于那些不能获得双单倍体的动植物物种而言,可使用自交群体或回交群体作图。如果作高密度的分子标记     

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.     

hybridlab (version 1.0): a program for generating simulated hybrids from population samples

We present the computer program hybridlab 1.0 for simulating intraspecific hybrids from population samples of nuclear genetic markers such as microsatellites, allozymes or SNPs (single nucleotide polymorphisms). The program generates a user‐specified number of multilocus F1 hybrid genotypes between any pair of potentially hybridizing populations included in a standard input‐file of multilocus genotypes for population genetic analysis. This simple, user‐friendly program has a wide range of applications for studying natural and artificial hybridization; in particular, for evaluating the statistical power for individual assignment of parental and hybrid individuals. An example of application for Atlantic cod populations is given.     

Spatial and temporal genetic structure in a hybrid cordgrass invasion

Invasive hybrids and their spread dynamics pose unique opportunities to study evolutionary processes. Invasive hybrids of native Spartina foliosa and introduced S. alterniflora have expanded throughout San Francisco Bay intertidal habitats within the past 35 years by deliberate plantation and seeds floating on the tide. Our goals were to assess spatial and temporal scales of genetic structure in Spartina hybrid populations within the context of colonization history. We genotyped adult and seedling Spartina using 17 microsatellite loci and mapped their locations in three populations. All sampled seedlings were hybrids. Bayesian ordination analysis distinguished hybrid populations from parent species, clearly separated the population that originated by plantation from populations that originated naturally by seed and aligned most seedlings within each population. Population genetic structure estimated by analysis of molecular variance was substantial (F_ST=0.21). Temporal genetic structure among age classes varied highly between populations. At one population, the divergence between adults and 2004 seedlings was low (F_ST=0.02) whereas at another population this divergence was high (F_ST=0.26). This latter result was consistent with local recruitment of self-fertilized seed produced by only a few parental plants. We found fine-scale spatial genetic structure at distances less than ∼200 m, further supporting local seed and/or pollen dispersal. We posit a few self-fertile plants dominating local recruitment created substantial spatial genetic structure despite initial long-distance, human dispersal of hybrid Spartina through San Francisco Bay. Fine-scale genetic structure may more strongly develop when local recruits are dominated by the offspring of a few self-fertile plants.