INTROGRESSION OF COYOTE MITOCHONDRIAL DNA INTO SYMPATRIC NORTH AMERICAN GRAY WOLF POPULATIONS

Mitochondrial DNA (mtDNA) genotypes of gray wolves and coyotes from localities throughout North America were determined using restriction fragment length polymorphisms. Of the 13 genotypes found among the wolves, 7 are clearly of coyote origin, indicating that genetic transfer of coyote mtDNA into wolf populations has occurred through hybridization. The transfer of mtDNA appears unidirectional from coyotes into wolves because no coyotes sampled have a wolf-derived mtDNA genotype. Wolves possessing coyote-derived genotypes are confined to a contiguous geographic region in Minnesota, Ontario, and Quebec, and the frequency of coyote-type mtDNA in these wolf populations is high (>50%). The ecological history of the hybrid zone suggests that hybridization is taking place in regions where coyotes have only recently become abundant following conversion of forests to farmlands. Dispersing male wolves unable to find conspecific mates may be pairing with female coyotes in deforested areas bordering wolf territories. Our results demonstrate that closely related species of mobile terrestrial vertebrates have the potential for extensive genetic exchange when ecological conditions change suddenly.     

Patterns of population subdivision, gene flow and genetic variability in the African wild dog (Lycaon pictus)

African wild dogs are large, highly mobile carnivores that are known to disperse over considerable distances and are rare throughout much of their geographical range. Consequently, genetic variation within and differentiation between geographically separated populations is predicted to be minimal. We determined the genetic diversity of mitochondrial DNA (mtDNA) control region sequences and microsatellite loci in seven populations of African wild dogs. Analysis of mtDNA nucleotide diversity suggests that, historically, wild dog populations have been small relative to other large carnivores. However, population declines due to recent habitat loss have not caused a dramatic reduction in genetic diversity. We found one historical and eight recent mtDNA genotypes in 280 individuals that defined two highly divergent clades. In contrast to a previous, more limited, mtDNA analysis, sequences from these clades are not geographically restricted to eastern or southern African populations. Rather, we found a large admixture zone spanning populations from Botswana, Zimbabwe and south-eastern Tanzania. Mitochondrial and microsatellite differentiation between populations was significant and unique mtDNA genotypes and alleles characterized the populations. However, gene flow estimates (Nm) based on microsatellite data were generally greater than one migrant per generation. In contrast, gene flow estimates based on the mtDNA control region were lower than expected given differences in the mode of inheritance of mitochondrial and nuclear markers which suggests a male bias in long-distance dispersal.     

Y‐chromosome evidence supports asymmetric dog introgression into eastern coyotes

Hybridization has played an important role in the evolutionary history of Canis species in eastern North America. Genetic evidence of coyote–dog hybridization based on mitochondrial DNA (mtDNA) is lacking compared to that based on autosomal markers. This discordance suggests dog introgression into coyotes has potentially been male biased, but this hypothesis has not been formally tested. Therefore, we investigated biparentally, maternally, and paternally inherited genetic markers in a sample of coyotes and dogs from southeastern Ontario to assess potential asymmetric dog introgression into coyotes. Analysis of autosomal microsatellite genotypes revealed minimal historical and contemporary admixture between coyotes and dogs. We observed only mutually exclusive mtDNA haplotypes in coyotes and dogs, but we observed Y‐chromosome haplotypes (Y‐haplotypes) in both historical and contemporary coyotes that were also common in dogs. Species‐specific Zfy intron sequences of Y‐haplotypes shared between coyotes and dogs confirmed their homology and indicated a putative origin from dogs. We compared Y‐haplotypes observed in coyotes, wolves, and dogs profiled in multiple studies, and observed that the Y‐haplotypes shared between coyotes and dogs were either absent or rare in North American wolves, present in eastern coyotes, but absent in western coyotes. We suggest the eastern coyote has experienced asymmetric genetic introgression from dogs, resulting from predominantly historical hybridization with male dogs and subsequent backcrossing of hybrid offspring with coyotes. We discuss the temporal and spatial dynamics of coyote–dog hybridization and the conditions that may have facilitated the introgression of dog Y‐chromosomes into coyotes. Our findings clarify the evolutionary history of the eastern coyote.     

Monitoring coyote population dynamics by genotyping faeces

Reliable population estimates are necessary for effective conservation and management, and faecal genotyping has been used successfully to estimate the population size of several elusive mammalian species. Information such as changes in population size over time and survival rates, however, are often more useful for conservation biology than single population estimates. We evaluated the use of faecal genotyping as a tool for monitoring long-term population dynamics, using coyotes (Canis latrans) in the Alaska Range as a case study. We obtained 544 genotypes from 56 coyotes over 3 years (2000-2002). Tissue samples from all 15 radio-collared coyotes in our study area had > or = 1 matching faecal genotypes. We used flexible maximum-likelihood models to study coyote population dynamics, and we tested model performance against radio telemetry data. The staple prey of coyotes, snowshoe hares (Lepus americanus), dramatically declined during this study, and the coyote population declined nearly two-fold with a 1(1/2)-year time lag. Survival rates declined the year after hares crashed but recovered the following year. We conclude that long-term monitoring of elusive species using faecal genotyping is feasible and can provide data that are useful for wildlife conservation and management. We highlight some drawbacks of standard open-population models, such as low precision and the requirement of discrete sampling intervals, and we suggest that the development of open models designed for continuously collected data would enhance the utility of faecal genotyping as a monitoring tool.     

A Bayesian approach to the inference of diploid genotypes using haploid genotypes

Morris and Spieth (1978) described a method of calculating unbiased estimates of diploid genotype frequencies given information on the genotypes of haploid cells derived from diploid individuals. They concluded that three haploids per diploid would minimize sampling variance of genotype frequencies, given a fixed total number of haploids examined. If the identity of individual diploid genotypes is needed, Morris and Spieth (1978) stated that more haploids should be collected per diploid. We extend this work by showing from a Bayesian perspective that the probability of misclassification of individuals depends not only on the number of haploids sampled, but also on the genetic structure of the population since misclassification error will increase as the frequency of heterozygotes increases. Since information on the genetic structure (allele frequencies, inbreeding coefficient) of a population is rarely known prior to the initiation of an empirical study, the usefulness of our Bayesian approach is in experimental design, by revealing the magnitude of possible misclassification errors given a particular choice of number of haploids.     

Spatial genetic and morphologic structure of wolves and coyotes in relation to environmental heterogeneity in a Canis hybrid zone

Eastern wolves have hybridized extensively with coyotes and gray wolves and are listed as a ‘species of special concern’ in Canada. However, a distinct population of eastern wolves has been identified in Algonquin Provincial Park (APP) in Ontario. Previous studies of the diverse Canis hybrid zone adjacent to APP have not linked genetic analysis with field data to investigate genotype‐specific morphology or determine how resident animals of different ancestry are distributed across the landscape in relation to heterogeneous environmental conditions. Accordingly, we studied resident wolves and coyotes in and adjacent to APP to identify distinct Canis types, clarify the extent of the APP eastern wolf population beyond the park boundaries and investigate fine‐scale spatial genetic structure and landscape–genotype associations in the hybrid zone. We documented three genetically distinct Canis types within the APP region that also differed morphologically, corresponding to putative gray wolves, eastern wolves and coyotes. We also documented a substantial number of hybrid individuals (36%) that were admixed between 2 or 3 of the Canis types. Breeding eastern wolves were less common outside of APP, but occurred in some unprotected areas where they were sympatric with a diverse combination of coyotes, gray wolves and hybrids. We found significant spatial genetic structure and identified a steep cline extending west from APP where the dominant genotype shifted abruptly from eastern wolves to coyotes and hybrids. The genotypic pattern to the south and northwest was a more complex mosaic of alternating genotypes. We modelled genetic ancestry in response to prey availability and human disturbance and found that individuals with greater wolf ancestry occupied areas of higher moose density and fewer roads. Our results clarify the structure of the Canis hybrid zone adjacent to APP and provide unique insight into environmental conditions influencing hybridization dynamics between wolves and coyotes.     

Estimation of population bottlenecks during systemic movement of tobacco mosaic virus in tobacco plants

More often than not, analyses of virus evolution have considered that virus populations are so large that evolution can be explained by purely deterministic models. However, virus populations could have much smaller effective numbers than the huge reported census numbers, and random genetic drift could be important in virus evolution. A reason for this would be population bottlenecks during the virus life cycle. Here we report a quantitative estimate of population bottlenecks during the systemic colonization of tobacco leaves by Tobacco mosaic virus (TMV). Our analysis is based on the experimental estimation of the frequency of different genotypes of TMV in the inoculated leaf, and in systemically infected leaves, of tobacco plants coinoculated with two TMV genotypes. A simple model, based on the probability that a leaf in coinoculated plants is infected by just one genotype and on the frequency of each genotype in the source, was used to estimate the effective number of founders for the populations in each leaf. Results from the analysis of three leaves per plant in plants inoculated with different combinations of three TMV genotypes yielded highly consistent estimates. Founder numbers for each leaf were small, in the order of units. This would result in effective population numbers much smaller than the census numbers and indicates that random effects due to genetic drift should be considered for understanding virus evolution within an infected plant.     

Mitochondrial DNA polymorphism in three Antillean island populations of the fruit bat, Artibeus jamaicensis

The Neotropical fruit bat, Artibeus jamaicensis, occurs throughout Latin America and on many islands in the Caribbean. Populations from Jamaica (in the Greater Antilles) to Barbados (in the Lesser Antilles) have been classified as a subspecies (A.j. jamaicensis) separate from that on the Lesser Antillean island of St. Vincent (A.j. schwartzi). Mitochondrial DNA (mtDNA) was isolated from 54 individuals collected on these islands, analyzed by digestion with restriction endonucleases, and the restriction sites were mapped. Three different mtDNA genotypes (16,000 +/- 200 bp) were identified: J-1 (16 animals from Jamaica, one from St. Vincent, 15 from Barbados), J-2 (two animals from Jamaica), and SV-1 (18 animals from St. Vincent, two from Barbados). The J-1 and J-2 genotypes were estimated to differ from each other by only 0.4%, but the SV-1 genotype differed from J-1 and J-2 by 8.1%-10.5%. The estimated sequence divergence between SV-1 and J-1 is unusually large for mammals that are regarded as conspecific. Restriction mapping showed that the differences among the genotypes (presence or absence of particular restriction sites) were located throughout the genome. The presence of the J-1 mtDNA genotype on Jamaica and on St. Vincent and Barbados (1,400 km away) demonstrates that maternal lineages in these bats are not necessarily confined to single islands or limited geographic regions. The presence of the J-1 mtDNA genotype within the A. j. schwartzi population on St. Vincent and the presence of the SV-1 genotype in two specimens of A. j. jamaicensis from Barbados document genetic exchange between subspecific populations on these islands, which are separated by 180 km of open water.     

Phytophthora infestans isolates from Northern China show high virulence diversity but low genotypic diversity

Phenotypic and genotypic characteristics of 48 Phytophthora infestans isolates , collected in five provinces in Northern China between 1997 and 2003, were determined and compared with reference isolates. Characterisation included mating type, virulence, mitochondrial DNA (mtDNA) haplotype and DNA fingerprinting patterns based on simple sequence repeats (SSR) and amplified fragment length polymorphisms (AFLP). All isolates had the A1 mating type, mtDNA haplotype IIa and an identical SSR genotype (designated as SG-01-01) that differed from SSR genotypes found in the reference isolates, including those representing the 'old' US-1 lineage that dominated the P. infestans population worldwide prior to 1980. In contrast, the virulence spectra were highly variable and virulence to all resistance genes present in the standard differential set ( R1 to R11 ) was found. AFLP analysis revealed some diversity; eight different AFLP genotypes were found that could be grouped into two major clusters. This study shows that there is very little genotypic diversity in the P. infestans population in Northern China. The occurrence of many different races within this rather uniform population is discussed in the framework of recent insights into the molecular determinants of avirulence in potato– P. infestans 'gene-for-gene' interactions.     

Multiple genotypes of mitochondrial DNA within a horse population from a small region in Yunnan province of China

mtDNA genotypes of six domestic horses (three adult short horses whose heights are under 1 m and three common domestic horses) from a small region of 15 km² in Malipo county of Yunnan province of China were investigated by the technique of restriction fragment length polymorphism (RFLP) with 16 restriction endonucleases which recognize 6-bp sequences. An average of 56 fragments for an individual was obtained. Unlike other domestic animals, this population of horses exhibits high mtDNA genetic diversity. Each of the six horses has a specific mtDNA genotype showing a pattern of multiple maternal origins, as suggested by fossil and literature records. We think the population of horses is an amazing seed-resource pool of horses and hence deserves to be paid more attention from the view of conservation genetics. However, it is also remarkable that we did not find any typical mtDNA genetic markers which would discriminate between short horses and common domestic horses.     

Resequencing studies of nonmodel organisms using closely related reference genomes: optimal experimental designs and bioinformatics approaches for population genomics

Decreasing costs of next‐generation sequencing (NGS) experiments have made a wide range of genomic questions open for study with nonmodel organisms. However, experimental designs and analysis of NGS data from less well‐known species are challenging because of the lack of genomic resources. In this work, we investigate the performance of alternative experimental designs and bioinformatics approaches in estimating variability and neutrality tests based on the site‐frequency‐spectrum (SFS) from individual resequencing data. We pay particular attention to challenges faced in the study of nonmodel organisms, in particular the absence of a species‐specific reference genome, although phylogenetically close genomes are assumed to be available. We compare the performance of three alternative bioinformatics approaches – genotype calling, genotype–haplotype calling and direct estimation without calling genotypes. We find that relying on genotype calls provides biased estimates of population genetic statistics at low to moderate read depth (2–8×). Genotype–haplotype calling returns more accurate estimates irrespective of the divergence to the reference genome, but requires moderate depth (8–20×). Direct estimation without calling genotypes returns the most accurate estimates of variability and of most SFS tests investigated, including at low read depth (2–4×). Studies without species‐specific reference genome should thus aim for low read depth and avoid genotype calling whenever individual genotypes are not essential. Otherwise, aiming for moderate to high depth at the expense of number of individuals, and using genotype–haplotype calling, is recommended.     

mtDNA and the origin of the Icelanders: deciphering signals of recent population history

Previous attempts to investigate the origin of the Icelanders have provided estimates of ancestry ranging from a 98% British Isles contribution to an 86% Scandinavian contribution. We generated mitochondrial sequence data for 401 Icelandic individuals and compared these data with >2,500 other European sequences from published sources, to determine the probable origins of women who contributed to Iceland’s settlement. Although the mean number of base-pair differences is high in the Icelandic sequences and they are widely distributed in the overall European mtDNA phylogeny, we find a smaller number of distinct mitochondrial lineages, compared with most other European populations. The frequencies of a number of mtDNA lineages in the Icelanders deviate noticeably from those in neighboring populations, suggesting that founder effects and genetic drift may have had a considerable influence on the Icelandic gene pool. This is in accordance with available demographic evidence about Icelandic population history. A comparison with published mtDNA lineages from European populations indicates that, whereas most founding females probably originated from Scandinavia and the British Isles, lesser contributions from other populations may also have taken place. We present a highly resolved phylogenetic network for the Icelandic data, identifying a number of previously unreported mtDNA lineage clusters and providing a detailed depiction of the evolutionary relationships between European mtDNA clusters. Our findings indicate that European populations contain a large number of closely related mitochondrial lineages, many of which have not yet been sampled in the current comparative data set. Consequently, substantial increases in sample sizes that use mtDNA data will be needed to obtain valid estimates of the diverse ancestral mixtures that ultimately gave rise to contemporary populations.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Locating hybrid individuals in the red wolf (Canis rufus) experimental population area using a spatially targeted sampling strategy and faecal DNA genotyping

Hybridization with coyotes (Canis latrans) continues to threaten the recovery of endangered red wolves (Canis rufus) in North Carolina and requires the development of new strategies to detect and remove coyotes and hybrids. Here, we combine a spatially targeted faecal collection strategy with a previously published reference genotype data filtering method and a genetic test for coyote ancestry to screen portions of the red wolf experimental population area for the presence of nonred wolf canids. We also test the accuracy of our maximum-likelihood assignment test for identifying hybrid individuals using eight microsatellite loci instead of the original 18 loci and compare its performance to the Bayesian approach implemented in newhybrids. We obtained faecal DNA genotypes for 89 samples, 73 of which were matched to 23 known individuals. The performance of two sampling strategies - comprehensive sweep and opportunistic spot-check was evaluated. The opportunistic spot-check sampling strategy required less effort than the comprehensive sweep sampling strategy but identified fewer individuals. Six hybrids or coyotes were detected and five of these individuals were subsequently captured and removed from the population. The accuracy and power of the genetic test for coyote ancestry is decreased when using eight loci; however, nonred wolf canids are identified with high frequency. This combination of molecular and traditional field-based approaches has great potential for addressing the challenge of hybridization in other species and ecosystems.     

Efficiency of haplotype frequency estimation when nuclear family information is included

In genetic studies the haplotype structure of the regarded population is expected to carry important information. Experimental methods to derive haplotypes, however, are expensive and none of them has yet become standard methodology. On the other hand, maximum likelihood haplotype estimation from unphased individual genotypes may incur inaccuracies. We therefore investigated the relative efficiency of haplotype frequency estimation when nuclear family information is included compared to estimation from experimentally derived haplotypes. Efficiency was measured in terms of variance ratios of the estimates. The variances were derived from the binomial distribution for experimentally derived haplotypes, and from the Fisher information matrix corresponding to the general likelihood function of the haplotype frequency parameters, including family information. We subsequently compared these variance ratios to the variance ratios for the case of estimation from individual genotypes. We found that the information gained from a single child compensates missing phase information to a high degree, resulting in estimates almost as reliable as those derived from observed haplotypes. Thus, if children have already been genotyped for other reasons, it is highly recommendable to include them into the estimation. If child information is not already present, it depends on the number of loci and the haplotype diversity if it is useful to genotype a single child just to reduce phase ambiguity. In general, if the number of loci is less than or equal to three or if the number of haplotypes with a frequency >5% is less than or equal to four, haplotype estimation from individuals is quite good already and the improvement gained from a single child can not compensate the genotyping effort for it. On the other hand, under scenarios with many loci and high haplotype diversity, haplotype frequency estimation from trios can be more efficient than haplotype frequency estimation from individuals also on a per genotype base.     

Reconciling nuclear microsatellite and mitochondrial marker estimates of population structure: breeding population structure of Chesapeake Bay striped bass (Morone saxatilis)

Comparative analyses of nuclear and organelle genetic markers may help delineate evolutionarily significant units or management units, although population differentiation estimates from multiple genomes can also conflict. Striped bass (Morone saxatilis) are long-lived, highly migratory anadromous fish recently recovered from a severe decline in population size. Previous studies with protein, nuclear DNA and mitochondrial DNA (mtDNA) markers produced discordant results, and it remains uncertain if the multiple tributaries within Chesapeake Bay constitute distinct management units. Here, 196 young-of-the-year (YOY) striped bass were sampled from Maryland's Choptank, Potomac and Nanticoke Rivers and the north end of Chesapeake Bay in 1999 and from Virginia's Mataponi and Rappahannock Rivers in 2001. A total of 10 microsatellite loci exhibited between two and 27 alleles per locus with observed heterozygosities between 0.255 and 0.893. The 10-locus estimate of R(ST) among the six tributaries was -0.0065 (95% confidence interval -0.0198 to 0.0018). All R(ST) and all but one theta estimates for pairs of populations were not significantly different from zero. Reanalysis of Chesapeake Bay striped bass mtDNA data from two previous studies estimated population differentiation between theta=-0.002 and 0.160, values generally similar to mtDNA population differentiation predicted from microsatellite R(ST) after adjusting for reduced effective population size and uniparental inheritance in organelle genomes. Based on mtDNA differentiation, breeding sex ratios or gene flow may have been slightly male biased in some years. The results reconcile conflicting past studies based on different types of genetic markers, supporting a single Chesapeake Bay management unit encompassing a panmictic striped bass breeding population.     

Bayesian inference of recent migration rates using multilocus genotypes

A new Bayesian method that uses individual multilocus genotypes to estimate rates of recent immigration (over the last several generations) among populations is presented. The method also estimates the posterior probability distributions of individual immigrant ancestries, population allele frequencies, population inbreeding coefficients, and other parameters of potential interest. The method is implemented in a computer program that relies on Markov chain Monte Carlo techniques to carry out the estimation of posterior probabilities. The program can be used with allozyme, microsatellite, RFLP, SNP, and other kinds of genotype data. We relax several assumptions of early methods for detecting recent immigrants, using genotype data; most significantly, we allow genotype frequencies to deviate from Hardy-Weinberg equilibrium proportions within populations. The program is demonstrated by applying it to two recently published microsatellite data sets for populations of the plant species Centaurea corymbosa and the gray wolf species Canis lupus. A computer simulation study suggests that the program can provide highly accurate estimates of migration rates and individual migrant ancestries, given sufficient genetic differentiation among populations and sufficient numbers of marker loci.     

Estimating major gene effects with partial information using Gibbs sampling

A method for estimating major gene effects using Gibbs sampling to infer genotype of individuals with unknown values, was compared with a standard mixed-model analysis. The purpose of this study was to evaluate the effect of including information of individuals with unknown genotypes on the estimates and their error variances (Ve) of the single-gene effects. When genotypes were known for all the individuals, results using the Gibbs method (GS) were similar to those obtained with the mixed model (MM). In the absence of selection, when information from individuals with unknown genotypes was included, GS yielded unbiased estimates of the major gene effects while reducing the Ve associated with them. This reduction in Ve depended on the gene frequency and mode of action of the major locus. For the additive effect, the reduction in Ve ranged from 29 to 69% of the total reduction which would have been obtained if all individuals had had a known genotype. Similarly the reduction in Ve found for the dominance effect ranged from 12 to 58%. Estimates using GS generally had small detectable biases when the polygenic heritability used in the analysis was inflated or estimated simultaneously. However, the benefit of using information from individuals with unknown genotypes was still maintained when comparing the mean square error of the estimates using either GS or MM when genotypes are only known for a subset of the population. When the population has been under selection, the use of Gibbs sampling to incorporate information of individuals without genotypes reduced substantially the bias and mean square error found for MM analysis on partial data. Nevertheless, there was some bias detected using Gibbs sampling. The gene frequency of the major gene in the base population was also well estimated despite its change over generations due to selection.     

Intense harvesting of eastern wolves facilitated hybridization with coyotes

Despite ethical arguments against lethal control of wildlife populations, culling is routinely used for the management of predators, invasive or pest species, and infectious diseases. Here, we demonstrate that culling of wildlife can have unforeseen impacts that can be detrimental to future conservation efforts. Specifically, we analyzed genetic data from eastern wolves (Canis lycaon) sampled in Algonquin Provincial Park (APP), Ontario, Canada from 1964 to 2007. Research culls in 1964 and 1965 killed the majority of wolves within a study region of APP, accounting for approximately 36% of the park's wolf population at a time when coyotes were colonizing the region. The culls were followed by a significant decrease in an eastern wolf mitochondrial DNA (mtDNA) haplotype (C1) in the Park's wolf population, as well as an increase in coyote mitochondrial and nuclear DNA. The introgression of nuclear DNA from coyotes, however, appears to have been curtailed by legislation that extended wolf protection outside park boundaries in 2001, although eastern wolf mtDNA haplotype C1 continued to decline and is now rare within the park population. We conclude that the wolf culls transformed the genetic composition of this unique eastern wolf population by facilitating coyote introgression. These results demonstrate that intense localized harvest of a seemingly abundant species can lead to unexpected hybridization events that encumber future conservation efforts. Ultimately, researchers need to contemplate not only the ethics of research methods, but also that future implications may be obscured by gaps in our current scientific understanding.