Maintaining genetic integrity of coexisting wild and domestic populations: Genetic differentiation between wild and domestic Rangifer with long traditions of intentional interbreeding

Abstract This study investigates the genetic effect of an indigenous tradition of deliberate and controlled interbreeding between wild and domestic Rangifer. The results are interpreted in the context of conservation concerns and debates on the origin of domestic animals. The study is located in Northeastern Zabaĭkal'e, Russia at approximately 57 degrees North latitude. Blood and skin samples, collected from wild and domestic Rangifer, are analyzed for their mtDNA and microsatellite signatures. Local husbandry traditions are documented ethnographically. The genetic data are analyzed with special reference to indigenous understandings of the distinctions between local domestic types and wild Rangifer. The genetic results demonstrate a strong differentiation between wild and domestic populations. Notably low levels of mtDNA haplotype sharing between wild and domestic reindeer, suggest mainly male‐mediated gene flow between the two gene pools. The nuclear microsatellite results also point to distinct differences between regional domestic clusters. Our results indicate that the Evenki herders have an effective breeding technique which, while mixing pedigrees in the short term, guards against wholesale introgression between wild and domestic populations over the long term. They support a model of domestication where wild males and domestic females are selectively interbred, without hybridizing the two populations. Our conclusions inform a debate on the origins of domestication by documenting a situation where both wild and domestic types are in constant interaction. The study further informs a debate in conservation biology by demonstrating that certain types of controlled introgression between wild and domestic types need not reduce genetic diversity.


| INTRODUCTION
Nature conservation aims to maintain biological diversity, which is often defined by the highest level of genetic diversity thought necessary to ensure the evolutionary potential for any species (Slatkin, 1987;Frankham, 2010). Genetic introgression from non-native, translocated representatives of the same species, or from domesticated versions of a species, is often represented as a key concern or even a threat in the conservation literature, as it may reduce fitness through the introduction of maladaptive traits or by disrupting positive epistatic interactions (Abbott, 2013;Mager, 2013;Wayne & Shaffer, 2016). On the other hand, bioarchaeologists often assume that domestic animals are descended from wild prototypes, and thereby imply that gene flow from wild to tame is necessary to maintain the fitness of the domestic type (Vigne, 2011;Zeder, 2012). However, little is known about the techniques that herders use to manage domestic populations in places where wild populations may pose a risk, or a benefit, to the local herds.
Understanding such management techniques would be particularly interesting, as they would necessarily influence the genetic integrity of both domesticated and wild populations. Traditionally, biologists have assumed that herdsmen enforced absolute reproductive isolation between wild and domestic forms (Price, 1984;Driscoll, 2009). This strategy would also maintain the genetic integrity of wild populations unintentionally. However, new models of domestication suggest that while spatial separation of wild and domestic types is often sought, various domestic species ranging from pigs to horses have also been cross-bred for specific characteristics by allowing controlled introgression from wild herds. (Ottoni, 2013;Jónsson, 2014;Frantz, 2015). For most domestic breeds, this idea remains a hypothesis as the original free-ranging wild forms have become largely extinct or restricted to a few isolated areas (Clutton-Brock, 1987). The exact husbandry techniques, which may have been employed in the Neolithic, can only be guessed at. However, this is not the case for reindeer.
Reindeer, Rangifer tarandus, are a species commonly considered to be in an early phase of domestication and which often coexist with wild forms (Baskin, 2000;Reimers & Colman, 2009). Today, almost 50% of the approximate 3,000,000 reindeer in the Old World are wild animals, and wild and domestic herds are managed in close coexistence in many areas of Eurasia and in Alaska (Syroechkovskii, 1995;Baskin, 2005). Therefore, this species provides a rare opportunity to link techniques of interbreeding of domestic herds and their wild relatives to their genetic signatures.
Despite their value as a proxy for studying the history of domestication, most existing theories of the origin of reindeer husbandry have not approached this problem genetically. The most authoritative theories are cultural historical models which survey the geographic diffusion of different styles of pastoralism (Wiklund, 1918;Maksimov, 1928;Vasilevich & Levin, 1951). They tend to focus on whether or not varying husbandry techniques could have been independently invented, or conversely, if they could have been simply copied and applied to local populations of wild Rangifer. They tend to assume that any Rangifer-wild or tame-can be harnessed given enough patience and if one had the right equipment.
The question of whether domestic reindeer are genetically distinct has never been tested.
Varying degrees of gene flow between wild and domestic reindeer have been reported in several recent studies (Mager, 2013;Colson, 2014;Røed, 2014). These studies primarily focus on degree of introgression of domestic lineages into the native wild gene pool, and which population and environmental factors may affect this. However, the complementary side to this model-the introgression of wild genes into domestic lineages-has been rarely studied. This process may in fact shed more light on the impact of introgression on species diversity. Northeastern Zaibaĭkal'e constitutes an area of approximately 132,000 square kilometers defined by the Stanovoe mountain range.
It is characterized by an alteration of alpine peaks of 1,200 m with year-round snow packs, interspersed with high open meadow plateaus (500-100 m) and stretches of taiga. Both wild and domestic reindeer herds coexist today in the region. It is likely that this coexistence has persisted for a long time. Zabaĭkal'e has been suggested to be one possible origin point for reindeer husbandry (Wiklund, 1918;Maksimov, 1928;Pomishin, 1990), with Okladnikov and Mazin (1976) speculating that its roots go back to the second millennium BCE. Furthermore, the region is particularly interesting because of its specific herding techniques. It has been documented that local Evenki herdsmen maintain a traditional practice of selectively cross-breeding domestic female Rangifer with wild males to produce offspring named locally as bai͡ unchikan [bai͡ unchikar -pl], which are often valued as transport reindeer (Shirokogoroff, 1929;Vodop'i͡ anov, 1970b;Davydov, 2014).
Local herdsmen claim that controlled intermixture with wild populations improves the resilience and strength of the domestic herd. One influential cultural-historical study used this practice of interbreeding as a criterion for distinguishing the style of reindeer husbandry in Zabaĭkal'e from that of neighboring regions (Vasilevich, 1964). Evenki reindeer husbandry in Zabaĭkal'e thereby presents itself as a unique case for studying how introgression of wild alleles into the domestic gene pool may have been managed during the early processes of domestication. To shed light on the potential for introgression and relevant circumstances, we combined ethnographic and genetic methods.
In this article, we test the hypothesis that Zabaĭkal Evenkis are able to enforce a strict genetic separation of wild and tame Rangifer, despite local breeding practices that tolerate or even to encourage controlled inter-breeding. We hypothesize further that their use of strict criteria to control genetic introgression through culling or castration is key to maintaining this divide. A strong genetic difference between the wild and domestic populations would therefore document their pastoral skill in maintaining and controlling specific reindeer pedigrees.
Establishing the identity and level of distinctiveness of wild and domestic herds is further important to settle an old debate on whether or not a domesticated harness reindeer constitutes a distinct type of reindeer, compared to wild reindeer.

| Study populations
The reindeer in this study were sampled between 2012 and 2014 from wild and domestic regional populations around the settlements of Chapo-Ologo and Ti͡ ani͡ a and the encampment at Nomama. We have included also data from a previous study of an Evenki herd at Lake Nichatka sampled in 2001 and 2002 (Røv & Abe, 2002;Røed, 2008). The name of each population corresponds to the residential base for the constellation of families holding local herds around that site and hunt wild reindeer in that same region (Figure 1). During the Soviet period, some of these residential bases would have been the headquarters for one centrally organized state farm, which exercised a great influence on the reproduction of each regional domestic population. We understand a regional domestic population to be a larger grouping of animals, which for historic and geographic reasons probably have been in constant interchange.
It is difficult to estimate the total number of wild and domestic Rangifer as this cultural historic region crosses four political districts and official statistics are fragmentary (Table A1). Existing government statistics record approximately 7,000 domestic reindeer in the region (although local informants claim that number is greatly underestimated) (Federal'nai͡ a Sluzhba Statistiki, 2016). A different set of accounts list approximately 24,000 head of wild reindeer (although that figure is said to be exaggerated) (Pavlov, 2016). In every context where we spoke with local herders, we were told that wild reindeer were rare.
Wild reindeer were listed as an endangered species in Buri͡ atii͡ a until (Boĭkov, 2005. The Chapo-Ologo regional domestic population consists of samples obtained from local herds with biographies suggesting they were offspring of the original Charskiĭ state farm, which in turn combined the smaller herds of dozens of indigenous families living in the region before collectivization. The domestic data set represents herds, or portions of herds, kept by four extended families. The set of wild reindeer associated with the Chapo-Ologo domestic population was hunted primarily in the Verkhniĭ Sakukan valley with a few samples from the Amudisy Lakes. The Ti͡ ani͡ a region captures the territory of two clan communities Ti͡ ani͡ a and Tokko which have restored reindeer husbandry after it nearly disappeared in 1993 when the local state farm was disbanded. The contemporary domestic regional population was assembled from a small population of a few dozen local domestic reindeer with influxes of five to fifty reindeer from Lake Nichatka, Chapo-Ologo, Perevoz, and Ust'-Ni͡ uzhka and another five regional F I G U R E 1 The locations of the local encampments and sampling sites, as well as regional populations, for wild and domestic Rangifer in northeastern Zabaĭkal'e. The boundaries of the local political districts, which bisect the region, are indicated with numbers encampments between 1993 and 2010. Domestic animals born in this territory were grouped together as the Ti͡ ani͡ a regional population. It represents the herds or portions of herds kept by five extended families. The wild reindeer samples of this region were from animals hunted on the Ori͡ us-Milele and Usu Rivers. The domestic Nomama samples were from domestic reindeer held by the clan community corporation "Ulutki"-essentially one extended family-which is headquartered at the ice-covered headwaters of the Lena River. This region was once the home to the state farm "Severnyĭ" that was disbanded in 1976. The region was left empty of domestic reindeer for 15 years, albeit with a large population of feral reindeer. Domestic reindeer husbandry was reintroduced in the region first in 1992-1993, according to informants by first lassoing the local feral reindeer, and more substantially by the purchase of thirty head of reindeer from the Chapo-Ologo region. The wild reindeer hunted in the Nomama region is likely separated from all the other wild herds by significant mountain ranges over large geographic distances (more than 500 km). In contrast to all the other sites, the extended family at Lake Nichatka never had participated in collectivization. Their kinship networks for the most part looked northwards the Bodaĭbo district of Irkutsk oblast' (Anderson et al., 2014) which facilitated the exchange small groups of 5-10 head of reindeer. The wild reindeer from this setting were hunted within 25 km of the main residential camp and could be considered to be part of a wild population which moves freely between this site and sites across the border to Ti͡ ani͡ a region. The previously published set from Lake Nichatka are documented only by sex and type (wild/domestic) with some specific biographical information on individual animals available in Abe (2005). Of the 35 domestic and 15 wild samples originally gathered, 31 domestic and 13 wild samples were used in this study. In 2012, we interviewed one Evenki herder who worked with this local herd for further details on this regional population.

| DNA extraction
DNA was obtained from hair follicles, skin samples, or from FTA cards designed for forensic work (Smith & Burgoyne, 2004). The tissue and blood samples were all noninvasive taken from animals culled, or ear-clipped, during normal pastoral activities. DNA extraction of skin samples was performed using DNaeasy Blood & Tissue Kit (Qiagen) following the manufactures protocol. DNA from hair follicles was extracted using the chelex method (Walsh, 1991). DNA from FTA cards was extracted from a 0.4 × 0.4 cm piece of the FTA cards using the boiling method, as described in Kvie (2016).
For both the CR and microsatellite data, we used Arlequin v.3.5 (Excoffier & Lischer, 2010) to calculate pairwise genetic distances (F ST ) to examine genetic differentiation among populations. We also used Arlequin to perform analysis of molecular variance (AMOVA) (Excoffier, 1992) to estimate the proportion of genetic variation found among populations (F ST ), among populations within groups (F SC ), and among groups (F CT ). The AMOVA was run with all samples, imposing no hierarchical structure, and with the populations divided into two groups, separating between domestic and wild populations. Significance associated with the fixation indexes was evaluated through random permutation procedure using 10,000 permutations.
Bayesian assignment, as implemented in the program STRUCTURE 2.3.4 (Pritchard, 2000), was used to assess whether discontinuities existed in the distribution of genetic variation within the data set.
For each number of genetic clusters (K ∈ [1,8]), a model with uniform priors, admixture, correlated allele frequencies, 20,000 burn-ins, and 200,000 Markov Chain Monte Carlo (MCMC) iterations was run ten times. For each K-value, average posterior probability among runs and standard deviation (SD) was calculated. Clumpp (Jakobsson & Rosenberg, 2007) was used to choose the most representative of the ten runs. Data were graphed using STRUCTURE PLOT (Ramasamy, 2014). To identify the proportion of individuals with mainly domestic, wild, or mixed ancestry, we used the STRUCTURE analysis but with only two inferred populations, namely domestic and wild. A threshold posterior probability value of 0.2 < q < 0.8 was used to indicate individuals with admixed ancestry.

| Genetic analyses
We identified 18 unique CR haplotypes among the domestic populations and 20 unique CR haplotypes among the wild populations (Table 1). However, only two haplotypes are shared between the two gene pools, with two wild individuals, both from Ti͡ ani͡ a, having shared haplotype with local domestic reindeer ( Figure 2). Also, the F ST estimate for genetic differentiation between pooled wild and domestic was significant in the CR (F ST = 0.832 p < .0001) as well as in the microsatellites (F ST = 0.032, 0.010 < p < .05). Our results show a T A B L E 1 Amount of genetic variation in the microsatellite loci and in the mitochondrial CR in wild and domestic reindeer herds from northeastern Zabaĭkal'e. N = number of individuals. For microsatellites, Na gives number of different alleles, and H obs and H exp gives observed and expected heterozygosity. For the CR Nh gives number of different haplotypes F I G U R E 2 Median-joining network of mtDNA haplotypes in wild and domestic Rangifer in northeastern Zabaĭkal'e. Each circle represents unique haplotypes with area proportional to the number of reindeer sharing a haplotype, and with the "status" (wild or domestic) each color coded with wild in green and domestic in yellow higher significance level in the mtDNA dataset compared to the microsatellites. This may be explained by the clonal/maternal mode of inheritance of mtDNA resulting in an effective population size that is one-fourth of nuclear DNA. Smaller population size makes mtDNA more sensitive to population changes as it is more exposed to genetic drift (Birky, 1983).
High levels of pairwise genetic differences among the eight regional populations of wild and domestic Rangifer, respectively, was evident from the microsatellite markers (Table 2). Significant differences were found among all populations (p < .05), except between the regional domestic herd from Ti͡ ani͡ a and the domestic herds from Nomama and Nichatka (p = .06), and between the wild herd from Chapo-Ologo and the wild herds from Lake Nichatka (p = .09) and Ti͡ ani͡ a (p = .208). Pairwise genetic differentiation calculated from the CR data show no or low levels of significance, except for the regional population of Nichatka wild reindeer which is significantly different from all other populations (p < .001) ( Table 3). For the CR, these tests are not very informative due to relatively low sample sizes obtained when dividing the wild reindeer into regional populations. We further identified nine to ten alleles in each of the two pooled microsatellite sets, and approximately six alleles in each regional population, noting high diversity in all populations (Table 1). As many as 17 and 18 private alleles (approximately 20% of the total number of alleles) were identified in the pooled domestic and wild populations, respectively. As the sample size of wild reindeer is half the size of domestic reindeer, we can assume that there is a higher number of private alleles in the wild populations, implying higher levels of introgression from wild to the domestic genepool than visa-versa.
An AMOVA analysis was run on the entire dataset without imposing any hierarchical structure. It showed that most of the genetic variation is found within populations, rather than among populations, in both datasets (Table 4). When grouping the eight populations into two types (wild and domestic), we found that 2.89% (p = .028) of the variance can be attributed to differences between wild and domestic, while 93.60% (p < .001) of the variation is found within populations in the microsatellite data. A similar result was evident for the mtDNA data where 4.43% (p = .67, not significant) of the variance is found between wild and domestic, while 83% (p < .001) is found within populations.
We analyzed the STRUCTURE results from K = 1 to K = 5 to assess the mean likelihood of there being separate populations within the dataset. The analysis suggested an apparent main structure at K = 2 ( Figure 3a), which can be interpreted as a division between wild and domestic herds. The population structure at K = 3 is characterized by a separation of most domestic Nomama reindeer from the other domestic herds, while at K = 4 it appears that a third domestic population begins to cluster itself around the Lake Nichatka population (Figure 3b). At K = 5, although with a greater variety in possible clusters, there is some evidence that the regional population of wild reindeer at Nomama distinguishes itself from other wild populations.
All five clusters are dominated (70%-80%) by individuals from one of the eight populations (Table 5). Except for wild and domestic reindeer at Ti͡ ani͡ a, all populations had most of its membership belonging to one of the clusters. It is significant that nearly 80% of the domestic reindeer at Nomama presented traits from the domestic population at Chapo-Ologo. The second STRUCTURE analysis, which was directed only at two inferred populations (wild and domestic), suggests that the Nomama and Ti͡ ani͡ a domestic herds have relatively high proportion of individuals with mixed ancestry ( Table 6). The proportion of individuals with q > 0.8 was generally higher among the wild populations than among the domestic populations, implying that the main direction of gene flow has been from the wild to the domestic genepool.
To further test the genetic structure among the domestic reindeer, we re-ran the STRUCTURE analysis with the 101 domestic samples alone. This analysis supported a separation into three domestic genetic clusters (Figure 4a) with each of Nomama, Lake Nichatka, and Chapo-Ologo regional herds characterizing the three clusters. The Tiana population appears to be made up of a mixture of all three groups (Figure 4b).

| Traditional Evenki herding practices for controlling introgression
Zabaĭkal Evenkis traditionally hold reindeer in small local herds scattered across the region ranging in size from a dozen to two hundred head by extended families of four to twelve individuals. The herd structure varies by the needs of each family and by season, but T A B L E 2 Pairwise genetic differences (F ST ) based on 13 reindeer specific microsatellite loci, among regional populations of domestic (D) and wild (W) herds of reindeer across northeastern Zabaĭkal'e. Population numbers are coded to the abbreviated list in the first column approximately up to 40% of the herd might be breeding females. The balance would be nonbreeding juveniles or calves, with a handful of males kept as breeding bulls and another dozen males castrated and kept for use as transport reindeer ( Figure 5). They generally keep their local herds within the confines of one watershed and are typically not moved more than 50 km over the course of the year. According to our own ethnographic survey, every 5 or 10 years, a herd might be split and then relocated into a neighboring valley forming a new local population. Individual domestic reindeer from one to twenty head, typically bulls, might occasionally be sold and exchanged across regions over a distance of 500 km or more. Overviews of Zabaĭkal Evenki husbandry are available in Orlov (1858); Shirokogoroff (1929), Fondahl (1989) and Anderson (1991); Anderson et al. (2014).
The highly trained domestic reindeer are mainly used for milking and saddling (for transport) and are rarely slaughtered for food (Vasilevich, 1964;Shnirel'man, 1977). This restriction does not apply to the bai͡ unchikar. The local herdsmen therefore seek to keep their camps near groups of wild forest reindeer, which are an important source of food and skins for clothing and equipment. The wild populations are described as scattered, small, localized groups of ten to twenty animals with a range that intersects with the domestic populations. They do not migrate more than 30 km over the course of a year and alternate their pastures between high alpine meadows and taiga meadows in between mountain ridges (Vodop'i͡ anov, 1970a; Ovdin & Ovdin, 2007). An important detail in the interrelationships between these two populations is that the breeding cycle of each is offset. In this region, wild Rangifer drop their calves in May 2-3 weeks later than domestic Rangifer.
In the local Evenki dialect, wild Rangifer (sakzhoĭ) and domestic Rangifer (oron) are marked semantically as different types. The local practice of deliberately interbreeding domestic cows with wild bulls to produce bai͡ unchikar is noted in the literature but has never been described in detail. The first European to mention the practice was Leopol'd Shrenk (1883) who travelled in the region in the 1850s. One clue to the possible ancient nature of the practice is the fact that all the terminology is in Evenki and not in Russian (or Latin). In our interviews, we documented an "ideal" breeding procedure as practiced in post-Soviet conditions, which might differ from how this practice was done in Soviet or Imperial times. It was explained that during the autumn rut, the majority of the domestic female herd is sequestered to prevent them from interbreeding. In the case of a small herd, they might be kept overnight in a specially built corral. In the case of a larger herd, they would be herded into an easily monitored valley. In both cases, herdsmen would keep shifts night and day to monitor the herd.
However, one or two breeding females-that is, much less than 10% of the female domestic population-might be tethered for cross-breeding with a wild male. The wild male is not selected as such. He presents himself opportunistically. The interbreeding happens under close observation of the herders. The wild male is often shot after mating, so that he will not harm the female-or, worse, fission-off a group of T A B L E 3 Pairwise genetic differences (F ST ) in the mitochondrial CR among regional populations of wild (W) and domestic (D) reindeer across northeastern Zabaĭkal'e. Population numbers are coded to the abbreviated list in the first column made about the mixed-blood calves, it is that they often resist walking together with the herd as a whole that they keep to the edge of the forest and become unruly when enclosed. It is the process of observation and selection which controls genetic introgression. Calves, which present themselves of "wild," may be allowed to survive one or 2 years as juveniles but would be culled before they reach breeding age. It would be fair to say that the majority of the bai͡ unchikar in normal economic conditions are not allowed to breed. Female baikunchar might be slaughtered for meat before they reach reproductive age. Male crossbreeds may be slaughtered or might be prevented from breeding by castration or by trimming their antlers to prevent them from successfully competing with other males. Some male bai͡ unchikan castrates are trained to harness to be cargo-carrying reindeer. Those few mixed breeds allowed to interbreed would be monitored year by year and could always be castrated or culled over the next breeding season. Their calves in the next generation would also be scrutinized.
In our field research, the bai͡ unchikar presented to us were primarily calves. The two adults were castrated males, one of which was "onequarter wild." None of the declared mixed-breeds at the time of our research could lead to genetic introgression.
The degree to which introgression proceeds depends on local conditions. In our field research, the two adult castrates were born in 2005, one from a male reindeer who was himself a bai͡ unchikan likely born around 2000. The same dates would apply to one bai͡ unchikan adult female sampled at Nichatka. This year was a time of severe economic crisis in the region. It was economically impossible hire mechanized transport to bring breeding bulls in from afar, forcing the local herders to recruit the wild population to ensure that all of their cows were covered. In a dire situation like this, one would assume that a higher proportion of the interbred calves would be allowed to survive. This precise situation presented itself before us at Nomama. The previous autumn a wild bull was permitted to cover some of the cows. When bears killed the last breeding domestic male that summer, the herders were making ready to interbreed the remaining domestic females with another wild bull. However, economic opportunities can also encourage interbreeding. There is one published estimate that in 1970 the Soviet state farms of Zabaĭkal'e held between 750 and 800 "hybrid" reindeer out of a total domestic population of 25,000 (Vodop'i͡ anov, 1970b). That article makes clear that the peculiar economic conditions of Soviet state socialism-the use of head-counts and calf-weights to monitor productivity-encouraged herders to wildly increase the size of the herds through the use of interbreeding. The article also makes clear that the calves are rarely kept beyond a year and a half of age.
Finally, all herders will admit candidly that some bai͡ unchikar are born  spontaneously when a small group of domestic females is covered by a wild male unintentionally. The same strict observational rules would apply to those happenstance wild/tame calves.

| DISCUSSION
Our results suggest that the reindeer herders in northeastern Zabaĭkal'e have developed herding techniques that maintain the genetic integrity of coexisting and overlapping populations of wild and domestic reindeer, even while intentional interbreeding them. This implies that the marked genetic differentiation between wild and domestic reindeer previously reported for the Lake Nichatka population (Røed, 2008) reflects a general pattern applicable for other herds across the region as well.
From our genetic analyses, the AMOVA shows that the majority of genetic variation is found within the populations, not between groups, more specifically between wild and domestic. Nevertheless, pairwise genetic distances among all populations reveal relatively high levels of differentiation, especially in the microsatellite data. Further, by pooling the wild and domestic populations, we found significant differentiation-which is also evident from the relatively high number of private alleles found in both wild and domestic pooled populations.
A separation between wild and domestic is further supported by the STRUCTURE analysis showing a main structure of two clusters, separating between wild and domestic, although some level of admixture is evident. However, little admixture was shown from the network analysis based on the mitochondrial control region, where only two of a total of 36 haplotypes are shared between wild and domestic reindeer. Significant genetic differences between wild and domestic reindeer have also been reported from Alaska (Cronin, 2003;Mager, 2013;Colson, 2014), Norway (Røed, 2014) and Greenland (Jepsen, 2002). In Norway, wild and domestic herds have been present for several centuries and are usually kept separated in different mountain areas, with enhanced migration barriers due to increasing infrastructure (Røed, 2014). Therefore, classic spatial segregation likely explains the genetic differentiation (Slatkin, 1987). Alaska and Greenland have a relatively recent history of coexistence of wild and domestic herds, which also have different geographic origins. Domestic reindeer were introduced to Alaska from Russia in 1890s (Ellanna, 2005), and to Greenland from Norway in 1952 (Jepsen, 2002). Differences in ancestry, and presumably associated local/regional adaptability, likely have contributed to reproductive isolation, as indicated by their distinct morphological and behavioral differences (Jepsen, 2002;Mager, 2013). This behavioral factor in selection addresses the recent hypothesis of the existence of "islands of domestication" within the genome wherein selection for desirable behavioral traits allows breeders to maintain domestic behavior within an interbred population (Frantz, 2015 This is consistent with female philopatry and male-biased dispersal commonly seen among most mammal species (Greenwood, 1980;Pusey, 1987;Goudet, 2002). This sex-biased structure has also been documented for wild Rangifer in Alaska (Roffler, 2012). In northeastern Zabaĭkal'e, however, this appears to be greatly aggravated by the strongly sex-biased and tightly controlled breeding practice. Our ethnographic work suggests that this significant sex-bias may have been further augmented by the preference for exchanging or importing male domestic reindeer to improve a herd. These small taiga herds are often assembled through kinship. Thus, particular reindeer might be gifted to one family or another when their paths cross in the forest. In our work on reindeer biographies, we discovered that male reindeer in particular might be exchanged several times over between friends and relatives in different locations across hundreds of kilometers with the animal spending 1 or 2 years in various camps (Davydov, 2014). Female animals, on the other hand, are more likely to be kept within the family.
Female reindeer would be more likely to be translocated when a herd would grow beyond a certain threshold and relocated to a new valley.
This was the case for the local population at Amudisy (Chapo-Ologo region) which had been fissioned twice as the herd approached 300 head. There are published ethnographic accounts of female reindeer being given in exchange as bride price in the 19th Century (Vasilevich, 1969). Further, female reindeer might be exchanged within large transfers of twenty to fifty head designed to augment a herd, such as the oral accounts of the building of the regional population at Ti͡ ani͡ a. Our analysis of proportional membership in Table 5 singles out Ti͡ ani͡ a as a population presenting no strong affinity to any one of the selected populations.
This study was originally designed to test the hypothesis that careful husbandry could successfully keep wild and domestic populations genetically separated despite living in close proximity and despite deliberate attempts to let certain female reindeer to interbreed. The genetic results indicate this is indeed the case. Our results further suggest that it is possible to distinguish local domestic clusters, which may serve as signatures of particular breeding traditions. This finding may be important to bolster local claims to the pedigree status of their herds which could give them access to Russian state funding to preserve pedigree types (Zabrodin & Borozdin, 1989).
This study suggests that it is possible to distinguish regional domestic herds genetically despite complex social histories of collectivization and privatization. The status of the Nomama population is a particularly interesting case in point. Despite significant influxes of reindeer purchased from neighboring regions, and Chapo-Ologo in particular, the Nomama population continues to distinguish itself from other groups by its low level of genetic diversity within the mtDNA and a strongly different microsatellite signature. The mtDNA signature supports an interpretation that there may have been a founding by a small number of females-perhaps those that oral accounts suggested were lassoed from the free-ranging feral remnants of the defunct state farm. Their distinct microsatellite signature suggests that the founding population may have been augmented by hybridization with local feral reindeer. Further, the proportion population membership in Table 6 suggests a high probability (50%) that the domestic population is best described as wild. In the older cultural-historical literature, the North Baikal region, where Nomama is located, corresponds to a distinctly different tradition of reindeer herding technique from all other groups (Vasilevich, 1964). In this case, decollectivation and the restoration of family-based reindeer husbandry between 1976 and 2012 may have produced a small-scale case study of origin of local reindeer husbandry within a distinct regional population of wild and feral Rangifer stressing the importance of strict selection and introgression along the male line.
Up until now, the classic literature on the origin of reindeer husbandry and of reindeer pedigrees has relied upon the classification of saddling and harnessing technologies as a proxy for grouping together distinct indigenous pastoral traditions. The reindeer bodies themselves be they wild or tame were assumed to be standard. This study suggests that in this region of southeastern Siberia, perhaps one of several "hearths" of reindeer husbandry, wild and domestic types are understood to be phenotypically and behaviorally different by local herders, and that this difference extends to the distinct genetic signatures in each regional population. Unlike in other models of the origin of animal husbandry, this study suggests that a respect for difference between the populations leads to a breeding strategy which enhances the genetic diversity of the species as a whole.

ACKNOWLEDGMENTS
The funding for the fieldwork and laboratory work for this study was provided by the ERC Advanced Grant 295458 Arctic Domus

CONFLICT OF INTEREST
None declared.

DATA ACCESSIBILITY
Genetic mtDNA data generated for this study is available on GenBank.
Reference numbers are in Supplementary Table S1.