RESUMO
Introduced species can threaten native taxa in multiple ways, including competition and hybridization, which can reduce fitness, alter ecological niches or swamp native genomes. Encroachment and hybridization by introduced species also provide opportunities to study the dynamics of invasiveness and hybridization during early stages following contact. We used 33 microsatellites, 51 single nucleotide polymorphisms and a mtDNA marker to characterize the extent and spatial pattern of encroachment and hybridization between a native, endemic subspecies of red fox (Vulpes vulpes patwin) and an introduced red fox population composed of highly admixed, phylogenetically divergent stock, resulting from a century of domestication. Both nuclear and mtDNA markers indicated that hybridization was primarily restricted to a narrow zone where the two populations came into contact. Although a few introgressed genotypes were detected in the interior of the native range, we found no immigrant foxes or F(1) or F(2) hybrids there, suggesting native foxes excluded introduced individuals. We speculate that the observed interbreeding at the periphery was facilitated by low densities. In total, 98% of mtDNA haplotypes in the native range were native and 96% of the nuclear ancestry was estimated to be native. Although the introduced range had expanded fivefold over the past four decades, native and non-native haplotypes from museum samples collected in and near the native range three decades earlier showed a similar geographic distribution as today, suggesting that the native range and hybrid zone were relatively stable. We hypothesize that the monogamous mating system of red foxes and other wild canids may enhance their resistance to hybridization because of greater fitness consequences associated with mate discrimination.