Increased genetic marker density reveals high levels of admixture between red deer and introduced Japanese sika in Kintyre, Scotland

Abstract Hybridization is a natural process at species range boundaries, but increasing numbers of species are hybridizing due to direct or indirect human activities. In such cases of anthropogenic hybridization, subsequent introgression can threaten the survival of native species. To date, many such systems have been studied with too few genetic markers to assess the level of threat resulting from advanced backcrossing. Here, we use 44,999 single nucleotide polymorphisms (SNPs) and the ADMIXTURE program to study two areas of Scotland where a panel of 22 diagnostic microsatellites previously identified introgression between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon). In Kintyre, we reclassify 26% of deer from the pure species categories to the hybrid category whereas in the NW Highlands we only reclassify 2%. As expected, the reclassified individuals are mostly advanced backcrosses. We also investigate the ability of marker panels selected on different posterior allele frequency criteria to find hybrids assigned by the full marker set and show that in our data, ancestry informative markers (i.e. those that are highly differentiated between the species, but not fixed) are better than diagnostic markers (those markers that are fixed between the species) because they are more evenly distributed in the genome. Diagnostic loci are concentrated on the X chromosome to the detriment of autosomal coverage.

of anthropogenic hybridization found that most cases of admixture result in the breakdown of reproductive isolation between previously distinct species (Grabenstein & Taylor, 2018). Biologists have flagged anthropogenic hybridization as a conservation concern, and the rapidity with which hybrid swarms are formed (Grabenstein & Taylor, 2018) suggests that this is a problem that is both common and urgent. Pure species may no longer exist as such, and the evolutionary trajectory of a species may be radically changed due to the introgression of large numbers of alleles which have already proved viable in an allopatric species.
Backcrossing of hybrid individuals with parental species leads to many admixed individuals that are difficult to differentiate from parental species individuals (McFarlane & Pemberton, 2019). When hybrids mate with parental species individuals (instead of other hybrid individuals), bimodal hybrid zones can be formed, where there are many individuals with a small proportion of introgression from another species. This leads to a large proportion of individuals of hybrid ancestry in the population that are difficult to distinguish either genetically or phenotypically from parental species. In the case of genetic analysis, this is because with each generation of backcrossing a doubling of markers is needed to detect introgression with the same degree of confidence (Boecklen & Howard, 1997). However, many current studies of anthropogenic hybridization use an average of 20 markers to detect introgression (Todesco et al., 2016), which is too few to reliably detect even second-generation backcrosses in most systems (Boecklen & Howard, 1997;Vähä & Primmer, 2006).
The limited detection power of the marker panels that are typically used in studies of anthropogenic hybridization is insufficient to detect ongoing introgression and it seems likely that many thousands of markers will be necessary to accurately identify the majority of hybrid individuals (McFarlane & Pemberton, 2019).
The use of genomic markers (i.e. thousands of SNPs as compared to tens of microsatellite markers) increases the detection probability of admixed individuals, particularly in systems where backcrossing has occurred for multiple generations. Theoretically with each generation of backcrossing, the proportion of individuals that are homozygous at diagnostic markers for alleles from one parental species increases, with the power to detect backcrosses sharply decreasing with each additional generation (Boecklen & Howard, 1997). For these reasons, it has been suggested that when determining how many markers to genotype researchers consider not only divergence between hybridizing species (as suggested by (Vähä & Primmer, 2006)), but also the number of elapsed generations since hybridization has begun and the recombination rate of the taxa studied (McFarlane & Pemberton, 2019). Empirically, use of high-density markers has indeed led to an increase in detection of backcrossed individuals. For example in Italy where wolves and dogs have hybridized, 1%-5% of wolves were found to be hybrid individuals when 16-18 microsatellite markers were used (Randi, 2008), whereas 62% of sampled individuals were determined to be hybrid when 170,000 single nucleotide polymorphism (SNP) markers were used (Pilot et al., 2018). Although the individuals sampled were not the same, and the study goals differed, it is telling that there is such a discrepancy in the estimated proportion of hybrid individuals when a high-density marker panel was used. Similarly, all individuals in a westslope cutthroat trout-rainbow trout hybrid zone that were previously thought to be parental species were discovered to be hybrid individuals when re-analysed with more markers (Boyer, Muhlfeld, & Allendorf, 2008;Hohenlohe et al., 2013). Use of higher density marker panels has the potential to change the qualitative and quantitative understanding of hybrid systems, including providing information on the direction of backcrossing, the number of generations since backcrossing began and the distribution (e.g. hybrid swarm, unimodal or bimodal) of the hybrid system (McFarlane & Pemberton, 2019).
Even with the use of genomic markers, not all markers will contribute equally to hybrid detection. Diagnostic markers (DM), those that have completely fixed allele differences between populations, have more power than markers that are polymorphic in either population (Pritchard, Stephens, & Donnelly, 2000;Vähä & Primmer, 2006). Ancestry informative markers (AIM), which have substantial allele frequency differences between parental populations, are also more powerful than markers which are highly polymorphic within one or both parental populations (Shriver et al., 2003). Thus, if a situation requires the use of fewer markers (e.g. to lower the cost of genotyping individuals, or because of the need for rapid testing (Senn et al., 2019)), it would be best to focus on diagnostic or ancestry informative markers, with the caveat that more markers will always lead to improved detection of hybrid individuals (McFarlane & Pemberton, 2019). For example, Galaverni et al. (2017) compared 48, 96 and 192 AIMs taken from a set of approximately 25,000 markers that had been thinned for linkage disequilibrium (LD). They found that each of these marker sets performed similarly, although there were no DM included in this comparison. A comparison of hybrid detection that includes both diagnostic markers and ancestry informative markers would contribute to confidence when researchers use AIM, especially since there are not always DM available.
It should be noted that both diagnostic and ancestry informative markers can have a biased distribution across the genome due to the evolutionary history of a population. If there is heterogeneity in introgression across the genome, then using too few markers may miss hot spots of introgression. For example, genes that are associated with postzygotic reproductive isolation are expected to be concentrated on sex chromosomes, due to the smaller effective population size and limited recombination of sex chromosomes (Qvarnström & Bailey, 2009). Empirically, a disproportionate number of genes associated with reproductive isolation have been traced to the sex chromosomes (Coyne & Orr, 2004), including in mouse hybrids (Payseur, Krenz, & Nachman, 2004), swordtail hybrids (Schartl, 1988) and flycatcher hybrids (Saether et al., 2007). Markers in strong LD with such genes are likely to be diagnostic (DM). When selecting markers for hybrid detection, the exclusive use of diagnostic markers will ignore introgression in other parts of the genome, where DM are not present. Further, thinning for LD is likely to exclude many DM, decreasing the power of the analysis to detect hybrid individuals.
AIM may also be disproportionately distributed across the genome, but perhaps not with the same intensity as markers that have fixed between species.
Hybridization between native red deer (Cervus elaphus) and Japanese sika (C. nippon) in Scotland is a scenario in which the use of high-density markers to explore the extent of backcrossing seems very likely to provide new insights. Sika were introduced from Japan to Ireland in 1860 and from Ireland to Scotland beginning in 1870 (Ratcliffe, 1987). As the hybridization that occurred after these introductions was facilitated by human actions, this is considered anthropogenic hybridization (Allendorf et al., 2001;Grabenstein & Taylor, 2018). Based on earlier work using microsatellite markers, there appear to be two different scenarios of admixture between red deer and sika in Scotland, which we have studied here using two separate populations. Sika were introduced to the Kintyre Peninsula in 1893, and phenotypically hybrid individuals were first reported in the region in the 1974 (Ratcliffe, 1987). Given a generation time of approximately 5 years in red deer (Coulson et al., 1998), we might expect at least 6-7 generations of backcrossing in this system, as our data set includes deer sampled as early as 2006. However, we cannot discount the possibility of hybridization prior to introduction to Scotland, as hybridization was recorded in the Irish source population in the 1880s (Powerscourt, 1884). Previous research in Kintyre used 22 diagnostic microsatellites to determine that there is a range of admixture proportions and that backcrossed individuals are common Smith et al., 2018). The use of a diagnostic mitochondrial marker demonstrated the existence of backcrosses with very low levels of nuclear introgression, beyond the detection of the microsatellite markers (Smith et al., 2018).
Hybrid individuals have phenotypes that are closely correlated with their admixture score . While intermediate hybrids may be detectable by phenotype, advanced backcrosses are easily confounded with the parental species and are difficult for deer hunters ("stalkers") to identify in the field (Senn & Pemberton, 2009;Smith, Carden, Coad, Birkitt, & Pemberton, 2014).
Sika were introduced to the north of Scotland between 1889 and 1900 (Ratcliffe, 1987). Hybrid individuals have been documented in this region, but in contrast to Kintyre, there does not seem to be a hybrid swarm-the few hybrids detected by microsatellites and mtDNA are the result of extensive backcrossing (Smith et al., 2018).
In this study, we screened populations in the NW highlands, where these previously detected hybrids were sampled.
Anthropogenic hybridization between red deer and Japanese sika provides an opportunity to use SNP markers to search for cryptic backcrossed individuals that were previously identified as parental individuals. Using substantially more markers than were previously available, this study has three aims: (1) to determine how many individuals are hybrids when scored using greatly increased marker density; (2) to determine the distribution of hybridization and introgression in each study area, specifically whether the hybrid zone is uni-or bi-directional and whether hybrids mate with other hybrids or with parental species individuals; (3) to study the performance F I G U R E 1 A map of Scotland, showing the approximate sampling locations for deer from Kintyre (white) and the NW Highlands (pink) that were genotyped on a 50K cervine SNP array. The map is from Google Maps, accessed using ggmap (Kahle & Wickham, 2013)  All samples were shot as encountered, that is without bias in terms of phenotype (Smith et al., 2018). Sample collection consisted of ear tissue and has previously been described in Senn and Pemberton (2009) and Smith et al. (2018). Samples were either preserved in 95% ethanol or frozen.

| Microsatellite genotyping, mtDNA haplotyping
All individuals included in this study were previously genotyped at 22 diagnostic microsatellites Smith et al., 2018). Most individuals were also typed at the mitochondrial DNA control region, which includes a 39-bp sequence with a variable number of tandem repeats, of which red deer have a single repeat and Japanese sika have three repeats (Cook, Wang, & Sensabaugh, 1999;Smith et al., 2018).

| DNA extraction for SNP analysis
We used the DNeasy Blood and Tissue Kit (Qiagen) according to the manufacture's instructions to re-extract DNA for SNP analysis, except that to obtain DNA at a sufficiently high concentration, in the last step we eluted twice in 25 µl TE. DNA concentration was assayed using the Qubit™ dsDNA BR Assay Kit (Invitrogen™), and any samples below 50 ng/µl were vacuum-concentrated, re-extracted or omitted from SNP analysis.

| SNP Genotyping
Single nucleotide polymorphism genotyping was carried out using the cervine Illumina iSelect HD Custom BeadChip using an iScan instrument following manufacturer's instructions (as in Huisman,  (Brauning et al., 2015), in addition to two pools of 12 pure species red deer from the Isle of Rum, Scotland, and 1 pool of 12 pure species sika (based on microsatellite data) from Kintyre, at 10× coverage. The majority of SNPs on the chip (53K attempted loci) were selected to be polymorphic within red deer, specifically within the Rum red deer. However, during the development of the chip 2,250 SNPs were included that were expected to be diagnostic between red deer and sika, and 2,250 SNPs were included that were expected to be diagnostic between red deer and wapiti. Within each marker set, SNPs were selected to be evenly spaced throughout the genome according to map positions in the bovine genome, with which the deer genome has high homology (Johnston, Huisman, Ellis, & Pemberton, 2017;Slate et al., 2002). To check for consistency between batches and calculate the repeatability of each SNP, we used the same positive control twice on each 96-well plate (Huisman et al., 2016). Genotypes were scored using GenomeStudio (Illumina), using the clusters from (Huisman et al., 2016). SNPs that did not successfully cluster according to Huisman et al. (2016) were clustered manually. Quality control was carried out in PLINK (Purcell et al., 2007). We excluded individual samples with a call rate of less than 0.90. We deleted loci with minor allele frequency less than 0.001 and/or call rate of less than 0.90. We did not exclude SNPs based on the Hardy-Weinberg equilibrium, as any diagnostic markers fixed between the two species are not expected to be in HWE in an admixed population.

| Estimating global F ST
To estimate global F st between red deer and sika in our two study areas, we considered only individuals that had been assessed as parental species from the phenotype by the ranger concerned. We used phenotype rather than the genotype from our ADMIXTURE analyses (see below) to exclude admixed individuals for this analysis to avoid a circular argument regarding the power we had to detect hybrid individuals. We used the R package "assigner" (Gosselin, Anderson, & Bradbury, 2016) to estimate Weir and Cockerham's F st (Weir & Cockerham, 1984), with 95% confidence intervals (CI) based on bootstrapping.

| Estimating admixture
For the microsatellite data, we used the STRUCTURE Q scores from Smith et al., 2018 to compare species assignment with the Q scores estimated using SNP markers (described below). However, unlike Smith et al. (2018), we used credibility intervals rather than a threshold cut-off to assign individuals to species categories (see description below).
For the SNP data, we used ADMIXTURE (Alexander & Lange, 2011) to estimate the admixture (Q) score for each individual.
ADMIXTURE uses maximum likelihood to estimate individual ancestry proportions (Alexander, Novembre, & Lange, 2009). We ran ADMIXTURE unsupervised separately on the deer from Kintyre and those from the NW highlands, as these are two separate populations and we did not want to inadvertently introduce population structure while examining hybridization (Gompert & Buerkle, 2016). We assumed two ancestral populations. We used the bootstrapping function in ADMIXTURE to estimate standard errors, from which we calculated 95% confidence intervals as 1.96*SE. We also used fastSTRUCTURE (Raj, Stephens, & Pritchard, 2014)  When we excluded markers based on LD (r 2 < .3, 12K SNPs remaining), we found a very high correlation between admixture scores with and without this LD pruning (using ADMIXTURE, Pearson's correlation = 0.9915, t = 172.1, df = 511, p < 2.2e-16) and that only 6 individuals were assigned to a difference species in the LD pruned analysis. Thus, all downstream analyses were done without accounting for LD.
To categorize animals into pure red deer, pure sika deer or hybrid, we determined whether the 95% confidence intervals for the Q scores overlapped 0.99999 for pure red deer or 0.00001 for pure sika deer. These are the highest and lowest possible Q scores produced by ADMIXTURE, and, with many markers, confidence intervals were extremely narrow. All other individuals were classified as "hybrids." These stringent criteria for assigning an individual to a parental species means that we are more likely to wrongly assign parental species individuals as hybrids than wrongly assign hybrid individuals as parental (i.e. low accuracy or a type II error with regard to wrongly assigning individuals as parental species versus low efficiency or a type I error as defined by (Vähä & Primmer, 2006)).

| Effect of varying SNP selection
We wanted to compare the efficiency of using our entire marker panel against using either only diagnostic markers, those fixed for different alleles in each parental population, ancestry informative markers, those with distinct allele frequency differences in the two parental populations (Shriver et al., 2003) or markers indiscriminately sampled from the SNP chip that were neither diagnostic nor ancestry informative (hereafter "random"). Diagnostic markers To compare the power of each marker type, without confounding by the number of markers used, we first determined the number of DM that were available, which was 629, and then randomly sampled this number of markers from each of the AIM (excluding DM) and the entire marker panel (excluding DM and AIM). We also included a marker set that was ten times as many randomly sampled SNPs We refer to these data sets as 45K, 629 DM, 629 AIM, 629 random and 6,290 random later in the text. We then reran the ADMIXTURE analysis with each of the five marker sets for the Kintyre sample and re-estimated the number of individuals in each class (red, sika, hybrid). We focused on the deer from Kintyre because the sample size was larger and there was prior evidence of a large number of hybrid individuals in this region , a finding confirmed below. To understand the uncertainty around the number of individuals in each class generated by the different marker sets, we ran each analysis 100 times and varied the "set seed" in ADMIXTURE (-s time) to allow for slight differences between replicates. Note that the 629 AIM and the 629 random and 6,290 random panels were a new random draw of markers for each bootstrap in the analysis. We present the mean estimated number of red deer, hybrid and sika called, as well as the standard errors around each estimate from the replicated analyses.
For each category of marker, we compared the number of deer per iteration that were identified to species class differently when compared to the 45K marker set. We refer to these deer as "mismatches." We used a linear model with the number of mismatches as the response variable explained by the marker category to ask which marker category gave results most similar to those from the 45K panel. We report the mean number of mismatches ± standard error, as well as an ANOVA of the linear model. Since we found large differences between marker groups, we explored the location of markers using the red deer linkage map (Johnston et al., 2017) as a possible explanation of the observed differences.

| RE SULTS
After quality control, there were 513 individual deer from ten sam-

| NW Highlands
In the NW Highlands, we found 45 red deer (mean Q = 0.9999 ± 0001), 59 sika (mean Q = 0.00001 ± 0.0000) and 4 hybrids (Q = 0.2333 ± 0.4223). This is similar to the 45 red deer, 61 sika and 2 hybrids that were previously reported in this group of individuals ( Figure 2b). The Q scores for three of these hybrid individuals ranged from 0.0199 ± 0.006 to 0.0251 ± 0.007, that is they were very sika like hybrids, while the 4th individual had a Q score of 0.8667 ± 0.133 (red like hybrid). Two of these hybrid individuals had red deer mtDNA; the other two had sika mtDNA. We also found 4 individuals that we assessed as pure species sika using the nuclear genome but that had red deer mitochondrial markers. Given the scarcity of hybrids in this region, we did not conduct further analysis on this data set.

| Diagnostic and Ancestry Informative Marker Panels in Kintyre
We found substantial differences in the number of individuals assigned to each species by each of five marker panels (all 45K, 629 DM, 629 AIM, 629 random sample and 6,290 random sample).
Generally, we found that the 629 diagnostic markers and the sample of 629 random markers estimated fewer hybrids in the population than the other marker groups (Figure 3). We found that the number of mismatches (deer identified differently by one of the

| D ISCUSS I ON
We found many more hybrid individuals on Kintyre using 45K SNPs than had previously been reported using 22 microsatellites. The majority of the newly discovered hybrids were those with either a very high or very low Q score, suggesting that we mainly have higher efficiency (i.e. we are more likely to assign true hybrids as hybrid individuals; Vähä & Primmer, 2006) than previous studies of this system Smith et al., 2018). However, it should be noted that we estimated a very similar F st as that reported in previous studies (0.53 compared with 0.58, using microsatellite markers selected to be diagnostic (Senn & Pemberton, 2009), suggesting that while increasing the number of markers allowed us to detect more hybrid individuals, more markers did not greatly affect the estimate of F st .
A possible concern is that AIMs do not represent hybridization, but rather incomplete lineage sorting (Sang & Zhong, 2000), and as a result, we are identifying some pure individuals as hybrids.
However, as the use of AIMs alone would result in relatively greater uncertainty, and we have assigned hybrids based on credible intervals that do not overlap with 0 or 1, an individual with both an extreme point estimate and high uncertainty would be assessed as a parental species. As the uncertainty due to AIMs is thus built-in to our analyses, we believe that we have been conservative in our assignment of hybrid individuals. Furthermore, the existence of many highly introgressed individuals was expected due to the detection in Kintyre of "mitochondrial hybrids" with microsatellite genotypes suggesting pure sika but with mitochondrial haplotypes of red deer (Smith et al., 2018) and from considerations about the power of the small number of microsatellites used in that study to detect hybrids (see McFarlane & Pemberton, 2019). For these reasons, incomplete lineage sorting is unlikely to account for the increased number of hybrid individuals we have assigned here.
Even with 45K markers, we are unlikely to detect all of the backcrosses that could be present in the Kintyre population. The most extreme Q score that we found in a hybrid (i.e. an individual with credible intervals that excluded 0 or 1) was 0.0078 ± 0.003, which would be the approximate average Q score for a sixth-generation backcross. Given that hybridization was already apparent 6-7 generations ago (Ratcliffe, 1987), there could be more extreme hybrid individuals that we could not detect. Indeed, there were thirteen individuals that we assessed as pure species sika, but that have red deer mitochondrial markers. These individuals appear to be backcrossed beyond our power to detect them using nuclear markers. be homozygous for parental alleles from one species for all informative markers that we have used here, making them impossible to detect (Boecklen & Howard, 1997).
A surprise in our analysis was that the diagnostic marker panel of 629 markers performed worse, in terms of finding hybrids, compared with equivalent-sized panels of ancestry informative markers ( Figure 3). Diagnostic markers are the gold standard in the field (Pritchard et al., 2000), so this observation requires explanation. markers," generated more mismatches with the full marker set than the 629 AIM panels, and the smaller sample size of 629 such markers was worse in this regard than the larger panel of 6,920 such markers. This is as expected if the markers are relatively uninformative but power is gained by having more of them.
The question arises as to why so many of the diagnostic markers were on the X chromosome. Although the selection criteria during the construction of the cervine SNP array included even spacing across the chromosomes (against the bovine map), both within the selected red deer polymorphic and the sika diagnostic loci, we suspect that the categorization into diagnostic versus ancestry informative markers has led to a concentration of X-linked markers into the diagnostic category. Since introgression is impeded on the sex chromosomes because of limited recombination, there may be long runs of fixed markers on the sex chromosomes in any hybridizing system (Baack & Rieseberg, 2007;Barton, 1979). While thinning markers for linkage will alleviate this problem to some extent, it will also limit the number of informative markers in an analysis and thus might not be the best solution when more markers will lead to higher power.
Many researchers in hybrid systems use anonymous and unmapped markers (e.g. from RADseq SNP discovery) and select diagnostic markers for analysis. Given the demonstrable downside to using diagnostic markers in our system (concentration on the X at the cost of autosome coverage, leading to fewer identified hybrids), we suggest that future analyses of anthropogenic hybrid zones should include both DM and AIM, and probably all possible markers, in their analyses. We anticipate that this will lead to increased discovery of hybrid individuals, as it has done in the red deer-sika bimodal hybrid zone on Kintyre.  . The newly detected hybrids tend towards a more southerly distribution than those hybrids that had previously been detected, although there were also newly detected hybrids in the north of Kintyre. These backcrossed individuals in the north and south regions of Kintyre provide an opportunity for sika gene flow into red deer in the rest of Scotland and thus the possibility that pure red deer could eventually become extinct in mainland Scotland.
Whether such a scenario occurs depends on the detailed evolutionary dynamics of the system, which at the moment we can only speculate about. In the Scottish landscape, sika are proving extremely successful in commercial timber forests, and stalkers find them challenging to manage by shooting due to their secretive behaviour, so introgressing sika alleles could confer these traits on red deer. Introgressed individuals (as detected by microsatellites) have intermediate phenotypes in terms of appearance and size , though differences in fecundity or female pregnancy rate were not detectable in the latter study. However, we are currently ignorant about how selection is operating on introgressed alleles and phenotypes. Future work could use higher density markers to re-examine the cline dynamics (both spatial and temporal) in Kintyre to determine whether there are temporal changes over longer time spans than hitherto studied, for example if the tension zone of hybridization has moved over time (Barton & Hewitt, 1985). In addition, it should be possible to study whether particular genome regions are introgressing faster than others and potentially what traits they underpin.
Finally, it is a debatable matter whether there is conservation management value in detecting such advanced backcrosses. If an individual is less than 0.08% (i.e. a 6-7th generation backcross as we have detected in the present study), should this individual be considered a hybrid when a management decision is made? This is a difficult question to answer (Allendorf et al., 2001), and we believe it should be answered based on policy and biological considerations, rather than determined based on the statistical power a given set of markers allows. A gene-based theory of conservation has been proposed (Petit, 2004), in which individuals of conservation concern are determined to be hybrids based on the alleles carried for specific genes of interest for the native parental population (i.e. rare alleles, or alleles associated with distinct phenotypes). However, such "gene-targeted conservation" is probably out of reach for this and many other study systems of conservation concern (Kardos & Shafer, 2018). More practically, a mixed strategy of preserving those individuals that are genetically parental species (i.e. narrow Q score CIs that overlap with 0 or 1) and score highly on emblematic red deer phenotypes (e.g. summer coat without spots, long pointed ears, large antlers) could reduce the threat to red deer from sika introgression. A similar strategy is currently in use to preserve both the genetics and phenotypes of Scottish wildcats, where individuals with both high phenotypic and genetic scores are preferentially protected and used in a captive breeding programme (Senn et al., 2019). Given the presence of so many advanced backcrosses in this system, any management considered should account for this introgression.

ACK N OWLED G EM ENTS
We thank the Forestry Commission Scotland rangers, particularly

DATA AVA I L A B I L I T Y S TAT E M E N T
All data and code are available at https ://figsh are.com/proje cts/ Incre ased_genet ic_marker_densi ty_revea ls_high_levels_of_admix ture_betwe en_red_deer_and_intro duced_sika_in_Kinty re_Scotl and/66860 .