No evidence for nuclear introgression despite complete mtDNA replacement in the Carpathian newt (Lissotriton montandoni)

Patterns of interspecific introgression may vary geographically, and the distribution of introgressed variants can yield insight into the historical dynamics of genetic interactions between hybridizing species. Urodele amphibians, often characterized by limited mobility, deep intraspecific genetic structuring and vulnerability to climatic changes, constitute suitable models for such historical inferences. Here, we combine an extensive survey of the mitochondrial (mtDNA) and nuclear (15 microsatellites) genomes in the Carpathian newt, Lissotriton montandoni (Lm) with species distribution modelling (SDM). Populations of the smooth newt, L. vulgaris (Lv) from the areas surrounding the Lm range were also sampled to test whether gene flow between these hybridizing species extends beyond the area of strict syntopy. The extent of introgression differs dramatically between the mitochondrial genome and the nuclear genome. While multiple, spatially and temporally distinct introgression events from Lv resulted in complete mtDNA replacement in Lm, there was little evidence of recent interspecific nuclear gene flow in the assayed markers. Microsatellite differentiation within Lm defines three units, probably derived from separate glacial refugia, located in the northern, eastern and southern part of the Carpathians. In situ survival and range fragmentation of Lm are supported by SDM, corroborating the role of the Carpathians as a major refugial area. Our results, in combination with previous reports of extensive introgression of the major histocompatibility complex (MHC) genes, emphasize the complexity of historical gene exchange between Lm and Lv.     

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker‐based analysis reveals cattle‐specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY_cattle = 2.66 ± 0.53% and Q_cattle = 0.69 ± 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai‐Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak‐cattle hybridization is primarily driven to produce F₁ hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make‐up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.     

Extreme mtDNA homogeneity in continental Asian populations

Mitochondrial DNA (mtDNA) variation in continental Asia has not been well-studied. Here, we report mtDNA HV1 sequences for 84 Xi'an and 82 Changsha Han Chinese, 89 Honshu Japanese, and 35 Vietnamese. Comparison of these sequences with other Asian mtDNA sequences reveals high variability within populations, but extremely low differentiation among Asian populations. Correlations between genetic distance and geographic distance, based on mtDNA and Y chromosome variation, indicate a higher migration rate in females than in males. This may reflect patrilocality, as suggested previously, but another plausible hypothesis is that the demographic expansion associated with the spread of agriculture in Asia may be responsible for the extreme genetic homogeneity in Asia.     

中国17个黄牛品种mtDNA变异特征与多态性分析

中国黄牛品种资源丰富,尚有28个地方固有品种.为了进一步深入了解这些宝贵遗传资源,本研究通过mtDNA变异特征与多态性分析揭示这些来自中国不同地域地方黄牛的母系起源与分子系统学特征.在17个品种84个体的mtDNA D-loop全序列中,一共检测到了102个核苷酸替代突变位.由此定义的53个单倍型被类聚为2个明显的单倍群：普通牛和瘤牛.mtDNA D-loop全序列变异的第一个特征是转换发生的频率远高于颠换;第二个特征是缺失与替代突变共存;第三个特征是缺失突变率比较高.所有D-loop全序列的核苷酸多样性和单倍型多样性分别为0.026 78±0.000 50和0.919±0.027.普通牛D-loop单倍型在北方牛种群中占有优势(80%～100%),而瘤牛单倍型在南方牛种群中占有优势(42.9%～100%),2种不同单倍型在中原牛种群中的分布也存在差异.2种不同单倍型在中国不同地域17个黄牛品种中的差异性分布揭示出了瘤牛mtDNA基因在中国黄牛中自南而北、由高到低的流动模式,这种基因流动模式的形成可能是由历史事件、地理隔离以及气候环境差异等造成的.     

Multiple hybridization events between Drosophila simulans and Drosophila mauritiana are supported by mtDNA introgression

The study of speciation has advanced considerably in the last decades because of the increased application of molecular tools. In particular, the quantification of gene flow between recently diverged species could be addressed. Drosophila simulans and Drosophila mauritiana diverged, probably allopatrically, from a common ancestor approximately 250,000 years ago. However, these species share one mitochondrial DNA (mtDNA) haplotype indicative of a recent episode of introgression. To study the extent of gene flow between these species, we took advantage of a large sample of D. mauritiana and employed a range of different markers, i.e. nuclear and mitochondrial sequences, and microsatellites. This allowed us to detect two new mtDNA haplotypes (MAU3 and MAU4). These haplotypes diverged quite recently from haplotypes of the siII group present in cosmopolitan populations of D. simulans. The mean divergence time of the most diverged haplotype (MAU4) is approximately 127,000 years, which is more than 100,000 years before the assumed speciation time. Interestingly, we also found some evidence for gene flow at the nuclear level because an excess of putatively neutral loci shows significantly reduced differentiation between D. simulans and D. mauritiana. Our results suggest that these species are exchanging genes more frequently than previously thought.     

Genetic diversity,introgression and relationships among West/Central African cattle breeds

Genetic diversity, introgression and relationships were studied in 521 individuals from 9 African Bos indicus and 3 Bos taurus cattle breeds in Cameroon and Nigeria using genotype information on 28 markers (16 microsatellite, 7 milk protein and 5 blood protein markers). The genotypes of 13 of the 16 microsatellite markers studied on three European (German Angus, German Simmental and German Yellow) and two Indian (Nelore and Ongole) breeds were used to assess the relationships between them and the African breeds. Diversity levels at microsatellite loci were higher in the zebu than in the taurine breeds and were generally similar for protein loci in the breeds in each group. Microsatellite allelic distribution displayed groups of alleles specific to the Indian zebu, African taurine and European taurine. The level of the Indian zebu genetic admixture proportions in the African zebus was higher than the African taurine and European taurine admixture proportions, and ranged from 58.1% to 74.0%. The African taurine breed, Muturu was free of Indian zebu genes while its counter Namchi was highly introgressed (30.2%). Phylogenic reconstruction and principal component analysis indicate close relationships among the zebu breeds in Cameroon and Nigeria and a large genetic divergence between the main cattle groups – African taurine, European taurine and Indian zebu, and a central position for the African zebus. The study presents the first comprehensive information on the hybrid composition of the individual cattle breeds of Cameroon and Nigeria and the genetic relationships existing among them and other breeds outside of Africa. Strong evidence supporting separate domestication events for the Bos species is also provided.     

Mitochondrial DNA variation in cattle of South China: origin and introgression

Ten restriction endonucleases were used to investigate the mitochondrial DNA restriction fragment length polymorphism (mtDNA RFLP) of 11 native cattle breeds and one cultivated cattle breed in South China. Twenty-three restriction morphs were detected, which can be sorted into five haplotypes. A phylogenetic tree of the haplotypes was constructed by using the 'upgMa' method. Our study showed that haplotype I and II are identical to the zebu (Bos indicus) and taurine (Bos taurus) haplotypes, respectively. Zebu and taurine were the two major origins of cattle populations in South China, and the zebu probably had more influence on the native cattle population than taurine did. Haplotype III is identical to haplotype I of yak (Bos grunniens), which was only detected in the Diqing cattle breed. Haplotype IV was detected for the first time. This haplotype, found only in Dehong cattle, might be from an independent domestication event, probably from another Bos indicus population. Divergence of haplotypes I and IV occurred about 268,000-535,000 years ago, much earlier than the 10,000-year history of cattle husbandry. Our results also suggest a secondary introgession of mtDNA from yak to Diqing cattle.     

Detection of mitochondrial DNA from domestic cattle in bison on Santa Catalina Island

In 1924, 14 American bison (Bison bison) were introduced to Santa Catalina Island, California and sporadically supplemented thereafter with additional animals. To reduce the herd and its impact on native vegetation, over 2000 animals have been exported during the past four decades. Today, the herd is estimated to contain around 250 individuals. Genetic analysis was performed on 98 animals removed from the island in 2004. Forty-four samples (45%) had domestic cattle mitochondrial DNA (mtDNA), 12 (12%) had previously reported bison haplotypes and 42 (43%) had a new haplotype differing by one base pair from a previously reported bison haplotype. A complement of five restriction enzymes was found to be useful in identifying bison with domestic cattle mtDNA.     

Male-mediated introgression of Bos indicus genes into Argentine and Bolivian Creole cattle breeds

The geographic distribution and frequency of Bos taurus and Bos indicus Y chromosome haplotypes amongst Argentine and Bolivian Creole cattle breeds were studied, using cytogenetic and molecular genetic techniques. A complete correspondence between Y chromosome morphology and the haplotype of the Y-linked microsatellite marker INRA 124 was found in all males examined. The taurine and indicine haplotypes were detected in 85.7 and 14.3% of the males studied, respectively, although these frequencies varied amongst the different breeds examined. The geographic distribution of this polymorphism suggests a pattern of zebu introgression in South America. The highest frequencies of the Zebu Y-chromosome are found in Brazilian populations (43-90%), in the eastern part of the continent, while it is absent in the southernmost breeds from Uruguay and Argentina. Bolivian breeds, at the centre of the continent, exhibit intermediate values (17-41%). This east/west and north/south gradient of male Zebu introgression could be explained by historical events and environmental factors.     

Evolutionary neutrality of mtDNA introgression: evidence from complete mitogenome analysis in roe deer

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein‐coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species‐specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein‐coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four‐fold higher d_N/d_S ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well‐preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.     

Geographic and phylogenetic evidence for dispersed nuclear introgression in a daphniid with sexual propagules

The role of among-species gene flow in eukaryotic evolution remains controversial. Putative hybrid lineages are common in water fleas, but their ecological success is often associated with polyploidy and the production of asexual propagules. Advanced hybrid lineages with sexual propagules are expected to be geographically restricted because their successful dispersal is contingent on overcoming fertility complications, assimilation by parent taxa, and competition with parent taxa. Here we provide evidence that a diploid lineage of Daphnia has been formed by introgression between distantly related species and attained a broad distribution (Nearctic) despite its requirement for sexual propagules. The evidence is based on geographical discordance, phylogenetic discordance, recombinant genotypes and additive genotypes of the nuclear internal transcribed spacer regions (ITS) and mitochondrial DNA. Additive genotypes also provided evidence of hybridization between introduced European Daphnia and North American Daphnia. We argue that the unique biology of Holarctic lacustrine water fleas and the spatial separation of lineages during Pleistocene glaciation have promoted hybridization and its evolutionary consequences.     

Sheep mitochondrial DNA variation in European, Caucasian, and Central Asian areas

Three distinct mitochondrial maternal lineages (haplotype Groups A, B, and C) have been found in the domestic sheep. Group B has been observed primarily in European domestic sheep. The European mouflon carries this haplotype group. This could suggest that European mouflon was independently domesticated in Europe, although archaeological evidence supports sheep domestication in the central part of the Fertile Crescent. To investigate this question, we sequenced a highly variable segment of mitochondrial DNA (mtDNA) in 406 unrelated animals from 48 breeds or local varieties. They originated from a wide area spanning northern Europe and the Balkans to the Altay Mountains in south Siberia. The sample included a representative cross-section of sheep breeds from areas close to the postulated Near Eastern domestication center and breeds from more distant northern areas. Four (A, B, C, and D) highly diverged sheep lineages were observed in Caucasus, 3 (A, B and C) in Central Asia, and 2 (A and B) in the eastern fringe of Europe, which included the area north and west of the Black Sea and the Ural Mountains. Only one example of Group D was detected. The other haplotype groups demonstrated signs of population expansion. Sequence variation within the lineages implied Group A to have expanded first. This group was the most frequent type only in Caucasian and Central Asian breeds. Expansion of Group C appeared most recently. The expansion of Group B involving Caucasian sheep took place at nearly the same time as the expansion of Group A. Group B expansion for the eastern European area started approximately 3,000 years after the earliest inferred expansion. An independent European domestication of sheep is unlikely. The distribution of Group A variation as well as other results are compatible with the Near East being the domestication site. Groups C and D may have been introgressed later into a domestic stock, but larger samples are needed to infer their geographical origin. The results suggest that some mitochondrial lineages arrived in northern Europe from the Near East across Russia.     

Morphological and molecular evidence for hybridization and introgression in a willow (Salix) hybrid zone

Hybrid zones provide biologists with the opportunity to examine genetic and ecological interactions between differentiated populations. Accurate identification of hybrid genealogies is considered a necessary prerequisite to understanding observed patterns of hybridization-related phenomena. We analysed molecular and morphological data from individuals in a hybrid zone between two species of willows (Salix sericea Marshall and S. eriocephala Michaux) and report the use of randomly amplified polymorphic DNA (RAPD), chloroplast DNA (cpDNA), and ribosomal DNA (rDNA) markers, as well as vegetative morphology and foliar chemistry data to identify individuals in terms of hybrid genealogy and to infer the direction and extent of backcrossing and introgression within the hybrid zone. A novel version of a maximum likelihood estimate approach (developed for this study) was used to calculate hybrid index scores from RAPD marker data; this method produced results similar to those obtained using traditional arithmetic methods. Distribution of rDNA, cpDNA, and chemistry data were examined within the graphical context of RAPD-based hybrid index score histograms and principal component analyses (PCA) on RAPD and morphology data. Seven of the 21 plants classified as S. eriocephala in the field were possible introgressants. Another plant presented an unequivocal example of backcrossed S. sericea chemistry and RAPD markers. Inter- and intraspecific chloroplast diversity found within the hybrid zone suggests both historic introgression (perhaps in a glacial refugium), and contemporary hybridization. Patterns of inheritance and expression within the hybrid zone suggest that morphological characters are often not expressed in a simple additive fashion, and problems associated with both morphological and molecular data are considered.     

Morphological and phylogenetic evidence for hybridization and introgression in a sea star secondary contact zone

Abstract. Glacial cycles and other climatic events have been widely invoked as factors promoting divergence, secondary contact, and hybridization between populations of terrestrial organisms, but the origin and fate of secondary contact in the sea is much less well understood. We studied the distribution of morphological and genetic variation in a northwest Atlantic zone of secondary contact between congeneric sea stars of Asterias that probably separated after the Pliocene as part of the trans-Arctic interchange. These species have similar reproductive biology and can hybridize in the laboratory. However, multivariate analysis of morphological traits scored from sea stars inside and outside the zone of secondary contact clearly indicated two clusters of phenotypes that corresponded to the two taxonomic species. A quantitative analysis of this clustering pattern did not support the hypothesis of a third grouping that might correspond to intermediate hybrid phenotypes. Known F₁ hybrids from laboratory matings grouped with one of the two taxonomic species. However, a survey of mtDNA sequence variation clearly indicated that ∼13% of individuals of one species ( Asterias forbesi ) are descendants of hybridization events that resulted in introgression of haplotypes of Asterias rubens into populations of A. forbesi . We conclude that morphological phenotypes are inadequate to identify hybrids of Asterias and their descendants, and that hybridization and introgression might be common in this secondary contact zone.     

Genetic evidence for the introgression of Western NR6A1 haplotype into Chinese Licha breed associated with increased vertebral number

There is evidence that NR6A1 is a strong candidate for being a causal gene underlying vertebral number in pigs. The Licha Black is one of the leanest Chinese indigenous pig breeds, having an average vertebral number of 21.5. The introgression of Western germplasm into Licha Black, resulting in increased vertebral number, has been assumed but is not confirmed. This study detected allele frequencies of the NR6A1 causative mutation (c.575T>C) in 519 pigs from three Western and seven Chinese breeds including Licha Black, and evaluated the genetic variation in a 650-kb region containing NR6A1 in the 10 breeds. Allele T for increased vertebral number was fixed in Western breeds. In contrast, this allele was very rare in most of the Chinese native breeds. Notably, the T allele was present in the Licha Black at a rather higher frequency (0.585) and in the Laiwu at lower frequency (0.250). As expected, selection pressure has wiped out the genetic variability in the 650 kb region in Western breeds. Conversely, Chinese indigenous breeds showed a high degree of genetic variability in this region. However, the Licha Black displayed dramatically reduced heterozygosity at the loci proximal to the causative mutation. Moreover, a high proportion (45.9%) of Licha Black pigs and a small number (21%) of Laiwu pigs had the Western NR6A1 haplotype, and the two breeds showed closer relationships with Western commercial breeds than other Chinese breeds in the phylogenic tree. When the results are taken together, this study supports the assumption that the Western NR6A1 haplotypes were introduced into Licha Black and possibly Laiwu and are associated with increased vertebral number.     

Genetic variation within and among domesticated Atlantic salmon broodstocks in British Columbia, Canada

Atlantic salmon have been reared in the British Columbia, Canada aquaculture industry since the early 1980s. No breeding programmes spanned the entire production period and pedigree records were not kept for broodstocks prior to or since importation. Of the three recognized industry strains, two are of European ancestry ('Mowi' from Norway and 'McConnell' from Scotland) and one is of North American heritage ('Cascade' from Gaspe, Quebec). We evaluated the amount and distribution of genetic variation within industry broodstocks by surveying microsatellite variation at 11 loci in 20 broodstock groups sampled from major production facilities. Allelic richness averaged 10.9 (range 5.8-13.8), compared with a value of 20.3 obtained for a North American wild population. Pairwise genetic distances (D(S)) between samples within strains were generally less than those between strains, with samples attributed to the same strain clustering together in a neighbour-joining dendrogram. Nevertheless, average distances between samples within the European strains were high (0.41 for Mowi; 0.71 for McConnell) but lower (0.06) for the Cascade strain. The reduced intra-sample and increased intra-strain genetic variation observed for the BC domesticated samples compared with wild populations was similar to observations for European domesticated Atlantic salmon. Evidence of introgression of the Cascade strain into European broodstocks was provided by the presence of large Ssa202 alleles (confined to North America in wild populations) in some Mowi and McConnell samples. Introgression likely also contributed to the decreased intercontinental genetic distance for the domesticated samples of this study compared with that observed for wild populations.     

An evaluation of introgression of Atlantic coast striped bass mitochondrial DNA in a Gulf of Mexico population using formalin-preserved museum collections

Striped bass Morone saxatilis populations in drainages along the Gulf of Mexico coast (Gulf) were depleted in the 1950s and 1960s, probably because of anthropogenic influences. It is believed that only the Apalachicola-Chattahoochee-Flint (A-C-F) river system continually supported a naturally reproducing population of Gulf lineage. Striped bass juveniles of Atlantic coast (Atlantic) ancestry were introduced to restore population abundances in the A-C-F from the late 1960s to the mid 1970s and in many other Gulf rivers from the 1960s to the present. We previously identified mtDNA polymorphisms that were unique to ? 60% of striped bass from the A-C-F and which confirmed the continued successful natural reproduction of striped bass of Gulf maternal ancestry within the system. However, the genetic relatedness of the extant A-C-F population to ‘pure’ Gulf striped bass was not addressed. In this study, we determined the frequency of a diagnostic mtDNA XbaI polymorphism in samples of ‘pure’ Gulf striped bass that were collected from the A-C-F prior to the introduction of Atlantic fish, that were obtained from museum collections, and that were originally preserved in formalin. PCR primers were developed that allowed for amplification of a 191-bp mtDNA fragment that contained the diagnostic XbaI restriction site. Using RFLP and direct sequence analyses of the PCR amplicons, we found no significant differences in mtDNA XbaI genotype frequencies between the archived samples and extant A-C-F samples collected over a 15-year period. This indicates that significant maternally mediated introgression of Atlantic mtDNA genomes into the A-C-F gene pool has not occurred. Additionally, we found no evidence of the unique Gulf mtDNA genotype in striped bass from extant populations in Texas, Louisiana and the Mississippi River. These results highlight the importance of the A-C-F as a repository of striped bass to restore extirpated Gulf populations and the potential use of museum collections in retrospective population studies.     

Historical vicariance and male-mediated gene flow in the toad-headed lizards Phrynocephalus przewalskii

Using mitochondrial and microsatellite DNA data and a population genetic approach, we tested male‐mediated gene flow in the toad‐headed lizards Phrynocephalus przewalskii. The mitochondrial DNA (ND2 gene), on the one hand, revealed two major lineages and a strong population genetic structure (F_ST = 0.692; F_ST′ = 0.995). The pairwise differences between the two lineages ranged from 2.1% to 6.4% and the geographical division of the two lineages coincided with a mountain chain consisting of the Helan and Yin Mountains, suggesting a historical vicariant pattern. On the other hand, the nuclear microsatellite DNA revealed a significant but small population genetic structure (F_ST = 0.017; F_ST′ = 0.372). The pairwise F_ST among the nine populations examined with seven microsatellite DNA loci ranged from 0.0062 to 0.0266; the assignment test failed to detect any naturally occurring population clusters. Furthermore, the populations demonstrated a weak isolation by distance and a northeast to southwest clinal variation, rather than a vicariant pattern. A historical vicariant event followed by male‐mediated gene flow appears to be the best explanation for the data. Approximately 2–5 Ma, climatic change may have created an uninhabitable zone along the Helan‐Yin mountain chain and initiated the divergence between the two mitochondrial lineages. With further climatic changes, males were able to disperse across the mountain chain, causing sufficient gene flow that eventually erased the vicariant pattern and drastically reduced the population genetic structure, while females remained philopatric and maintained the mitochondrial DNA (mtDNA) divergence. Although polygyny mating system and female philopatry may partially contribute to the reduced movement of females, other hypotheses, such as female intrasexual aggression, should also be explored.     

Unexpected population genetic structure of European roe deer in Poland: an invasion of the mtDNA genome from Siberian roe deer

Introgressive hybridization is a widespread evolutionary phenomenon which may lead to increased allelic variation at selective neutral loci and to transfer of fitness‐related traits to introgressed lineages. We inferred the population genetic structure of the European roe deer (Capreolus capreolus) in Poland from mitochondrial (CR and cyt b) and sex‐linked markers (ZFX, SRY, DBY4 and DBY8). Analyses of CR mtDNA sequences from 452 individuals indicated widespread introgression of Siberian roe deer (C. pygargus) mtDNA in the European roe deer genome, 2000 km from the current distribution range of C. pygargus. Introgressed individuals constituted 16.6% of the deer studied. Nearly 75% of them possessed haplotypes belonging to the group which arose 23 kyr ago and have not been detected within the natural range of Siberian roe deer, indicating that majority of present introgression has ancient origin. Unlike the mtDNA results, sex‐specific markers did not show signs of introgression. Species distribution modelling analyses suggested that C. pygargus could have extended its range as far west as Central Europe after last glacial maximum. The main hybridization event was probably associated with range expansion of the most abundant European roe deer lineage from western refugia and took place in Central Europe after the Younger Dryas (10.8–10.0 ka BP). Initially, introgressed mtDNA variants could have spread out on the wave of expansion through the mechanism of gene surfing, reaching high frequencies in European roe deer populations and leading to observed asymmetrical gene flow. Human‐mediated introductions of C. pygargus had minimal effect on the extent of mtDNA introgression.