Bottlenecks, drift and differentiation: the population structure and demographic history of sika deer (Cervus nippon) in the Japanese archipelago

We assessed genetic differentiation and diversity in 14 populations of sika deer (Cervus nippon) from Japan and four populations of sika deer introduced to the UK, using nine microsatellite loci. We observed extreme levels of differentiation and significant differences in diversity between populations. Our results do not support morphological subspecies designations, but are consistent with previous mitochondrial DNA analyses which suggest the existence of two genetically distinct lineages of sika deer in Japan. The source of sika introduced to the UK was identified as Kyushu. The underlying structure of Japanese populations probably derives from drift in separate glacial refugia and male dispersal limited by distance. This structure has been perturbed by bottlenecks and habitat fragmentation, resulting from human activity from the mid-nineteenth century. Most current genetic differentiation and differences in diversity among populations probably result from recent drift. Coalescent model analysis suggests sika on each of the main Japanese islands have experienced different recent population histories. Hokkaido, which has large areas of continuous habitat, has maintained high levels of gene flow. In Honshu the population is highly fragmented and is likely to have been evolving by drift alone. In Kyushu there has been a balance between gene flow and drift but all the populations have experienced high levels of drift. Habitat fragment size was not significantly associated with genetic diversity in populations but there was a significant correlation between habitat fragment size and effective population size.     

江西桃红岭国家级自然保护区梅花鹿生境适宜性评价

华南梅花鹿(Cervus nippon)被IUCN列入濒危物种,也是我国国家Ⅰ级重点保护动物。目前种群仅分布于江西、浙江、安徽等狭窄的区域内,形成多个孤立种群,生境破碎和丧失被认为是限制梅花鹿种群增长的主要原因。于2011年3月至2013年3月采用样线法和样方法对桃红岭国家级自然保护区梅花鹿栖息地进行了野外调查,利用空间模拟方法,结合地理信息系统(GIS)技术的空间分析功能,以植被类型、坡度、坡向、海拔和人类干扰活动作为评价因子进行了生境适宜性评价。结果表明,桃红岭地区以森林为主,各类林地面积约9 488.15 hm~2,占75.90%,植被类型分为落叶阔叶林、针叶林、常绿阔叶林、针阔混交林、竹林、灌丛、草丛和芭茅丛,面积分别为1664.57、1638.63、3438.21、1247.15、87.85、1143.88、60.92 hm~2和206.94 hm~2。在不考虑人类活动影响时,梅花鹿的适宜生境和次适宜生境面积分别是2233.99 hm~2和2980.24 hm~2,分别占保护区总面积的18.61%和24.83%;而考虑人类活动影响时,梅花鹿的适宜生境和次适宜生境面积分别是1224.04 hm~2和2164.70 hm~2,分别占保护区总面积的10.20%和18.04%。由于梅花鹿的生境受到居民点、主要道路、农田耕作、森林采伐等人类活动的强烈影响,导致大量适宜和次适宜生境丧失、隔离,景观破碎度指数由0.4345增加到0.5898。以潜在可利用生境面积计算,保护区梅花鹿环境容纳量为(568±160)只,而以实际可利用生境面积计算,则只能容纳(368±105)只。适宜生境的丧失和破碎可能是限制桃红岭梅花鹿国家级自然保护区梅花鹿种群恢复的重要因素,在此基础上,通过实际调查提出了管理措施。     

Fine‐scale genetic structure in an eastern Alpine black grouse Tetrao tetrix metapopulation

Understanding genetic consequences of habitat fragmentation is crucial for the management and conservation of wildlife populations, especially in case of species sensitive to environmental changes and landscape alteration. In central Europe, the Alps are the core area of black grouse Tetrao tetrix distribution. There, black grouse dispersal is limited by high altitude mountain ridges and recent black grouse habitats are known to show some degree of natural fragmentation. Additionally, substantial anthropogenic fragmentation has occurred within the past ninety years. Facing losses of peripheral subpopulations and ongoing range contractions, we explored genetic variability and the fine‐scale genetic structure of the Alpine black grouse metapopulation at the easternmost fringe of the species’ Alpine range. Two hundred and fifty tissue samples and non‐invasive faecal and feather samples of eleven a priori defined subpopulations were used for genetic analysis based on nine microsatellite loci. Overall, eastern Alpine black grouse show similar amounts of genetic variation (H_O = 0.65, H_E = 0.66) to those found in more continuous populations like in Scandinavia. Despite of naturally and anthropogenically fragmented landscapes, genetic structuring was weak (global F_ST < 0.05), suggesting that the actual intensity of habitat fragmentation does not completely hamper dispersal, but probably restricts it to some extent. The most peripheral subpopulations at the edge of the species range show signs of genetic differentiation. The present study gives new insights into the population genetic structure of black grouse in the eastern Alps and provides a more fine‐scale view of genetic structure than previously available. Our findings will contribute to monitor the current and future status of the population under human pressures and to support supra‐regional land use planning as well as decision making processes in responsibilities of public administration.     

Genetic variation and population structure of the Japanese sika deer (Cervus nippon) in Hokkaido Island, based on mitochondrial D-loop sequences

Mitochondrial DNA (mtDNA) D -loop region sequences (602 bp) from 141 samples of the sika deer Cervus nippon collected from Hokkaido Island of Japan were investigated to elucidate population genetic structure. All animals possessed seven repeat units (38 or 39 bp each) in the sequences. Comparison of the 602-bp sequences showed four sites of transitional mutations (A↔G or C↔T). Based on combination of the substitutions, six D -loop haplotypes (a–f types) were identified in the Hokkaido population, suggesting the occurrence of at least six maternal lineages. Distribution maps of the haplotypes constructed using the Geographic Information System showed that the distribution patterns differed from haplotype to haplotype. In particular, distribution of the major three types (a-, b-, and c-types) almost overlapped with three main areas of coniferous forests in Hokkaido. These results suggest that expansion of the sika deer population could have occurred through the habitat of coniferous forests after the historical bottleneck in Hokkaido.     

Testing hypotheses driving genetic structure in the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla

Habitat fragmentation has often been implicated in the decline of many species. For habitat specialists and/or sedentary species, loss of habitat can result in population isolation and lead to negative genetic effects. However, factors other than fragmentation can often be important and also need to be considered when assessing the genetic structure of a species. We genotyped individuals from 13 populations of the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla in Florida to test three alternative hypotheses regarding the effects that habitat fragmentation might have on genetic structure. A map of potential habitat developed from recent satellite imagery suggested that Brown‐headed Nuthatch populations in southern Florida occupied smaller and more isolated habitat patches (i.e. were more fragmented) than populations in northern Florida. We also genotyped individuals from a small, isolated Brown‐headed Nuthatch population on Grand Bahama Island. We found that populations associated with more fragmented habitat in southern Florida had lower allelic richness than populations in northern Florida (P = 0.02), although there were no differences in heterozygosity. Although pairwise estimates of F_ST were low overall, values among southern populations were generally higher than northern populations. Population assignment tests identified K = 3 clusters corresponding to a northern cluster, a southern cluster and a unique population in southeast Florida; using sampling localities as prior information revealed K = 7 clusters, with greater structure only among southern Florida populations. The Bahamas population showed moderate to high differentiation compared with Florida populations. Overall, our results suggest that fragmentation could affect gene flow in Brown‐headed Nuthatch populations and is likely to become more pronounced over time.     

Characteristics of Sasa nipponica grassland as a summer forage resource for sika deer on Mt Ohdaigahara, central Japan

Recently, the sika deer, Cervus nippon Temminck, population has increased on Mt Ohdaigahara, central Japan. The dwarf bamboo, Sasa nipponica Makino et Shibata, is a primary forage plant for sika deer in this area. To demonstrate the characteristics of S. nipponica grassland, especially as summer forage for sika deer, the habitat use intensity of sika deer was estimated by fecal densities, and biomass, growth rate, removal by deer and crude protein content were examined. Sika deer utilized the S.nipponica grassland on Mt Ohdaigahara during summer when the biomass, growth rate and crude protein content of S. nipponica were high. The recent increase in the deer population seems to be partly due to S.nipponica grassland being a favorable summer habitat.     

Molecular biodiversity in the moss <Emphasis Type="Italic">Leptodon smithii</Emphasis> (Neckeraceae) in relation to habitat disturbance and fragmentation

Bryophytes seem particularly suitable to investigate genetic diversity in relation to habitat disturbance due to their large employment as bioindicators and to the recent application of molecular markers to moss population studies. Genetic variation and structure were analysed in seven urban, extraurban and remote populations of Leptodon smithii, an epiphytic moss of Quercus ilex, a phanerogamic species of Mediterranean climax vegetation. A total of 210 individual shoots were DNA extracted and amplified with internal simple sequence repeat (ISSR) primers, and 54 haplotypes were identified. An uneven distribution of haplotype number and frequencies was observed among sites, with a higher number of haplotypes and more homogeneous haplotype frequencies in the extraurban/remote populations. Molecular diversity indices were overall higher in the extraurban sites than in the urban ones. Multilocus linkage disequilibrium values were in line with the occurrence of sexual/asexual reproduction in the seven populations. The isolation-by-distance model was not supported by Mantel test among sites; however, within-population fixation index (F_ST) highlighted a clear relation between genetic and physic distances among trees, suggesting a limited dispersal range for L. smithii’s spores. The genetic structure was mainly affected by population size, wood structure and extent, and genetic drift consequent to habitat fragmentation and human-induced disturbance.     

Population structure and genetic variation in the endangered Giant Kangaroo Rat (<Emphasis Type="Italic">Dipodomys ingens</Emphasis>)

Populations of the endangered giant kangaroo rat, Dipodomys ingens (Heteromyidae), have suffered increasing fragmentation and isolation over the recent past, and the distribution of this unique rodent has become restricted to 3% of its historical range. Such changes in population structure can significantly affect effective population size and dispersal, and ultimately increase the risk of extinction for endangered species. To assess the fine-scale population structure, gene flow, and genetic diversity of remnant populations of Dipodomys ingens, we examined variation at six microsatellite DNA loci in 95 animals from six populations. Genetic subdivision was significant for both the northern and southern part of the kangaroo rat’s range although there was considerable gene flow among southern populations. While regional gene diversity was relatively high for this endangered species, hierarchical F-statistics of northern populations in Fresno and San Benito counties suggested non-random mating and genetic drift within subpopulations. We conclude that effective dispersal, and therefore genetic distances between populations, is better predicted by ecological conditions and topography of the environment than linear geographic distance between populations. Our results are consistent with and complimentary to previous findings based on mtDNA variation of giant kangaroo rats. We suggest that management plans for this endangered rodent focus on protection of suitable habitat, maintenance of connectivity, and enhancement of effective dispersal between populations either through suitable dispersal corridors or translocations.     

东北虎豹国家公园梅花鹿的空间遗传格局及其影响因素

东北梅花鹿是东北虎豹国家公园主要的大型食草动物之一,是东北虎豹的主要猎物,对针阔混交林群落的维持有关键的作用,探究其遗传多样性及空间遗传格局对东北梅花鹿的保护以及国家公园生态系统的健康至关重要。在国家公园珲春保护区内,通过非损伤方法获得遗传样本,利用微卫星标记,研究该梅花鹿种群的空间遗传格局及其影响因素。结果表明：本研究区梅花鹿种群平均期望杂合度为0.721,遗传多样性较为丰富。有限的扩散能力常常导致种群在遗传距离上具有显著的空间自相关模式,本研究区梅花鹿种群在0-1km距离等级内在遗传距离上具有显著的空间自相关现象,据此可推测,该地区梅花鹿扩散距离为1km左右。STRUCTURE分析表明,珲春地区梅花鹿种群不存在明显的遗传分化。各种空间变量可以显著影响物种的遗传分化。本研究选取海拔、坡度、坡向、地表起伏率、人类干扰5个变量,研究其对梅花鹿种群遗传结构的影响,这5个变量多被认为与大中型哺乳动物扩散阻碍相关。依据5个变量建立了336个阻力模型,并进行偏曼特尔检验。其中,依据海拔、坡向、地表起伏率、人类干扰假设建立的246个阻力模型与遗传距离之间的关系并不显著,综合所有变量的15个生境适宜性模型阻力模型与遗传距离的关系也都不显著。在依据坡度假设建构的75个阻力模型中,只有1个模型与遗传距离有显著的正相关关系,该模型同时也是在控制空间自相关影响后,在所有模型中与遗传距离相关性最高的模型。根据该模型推测,最适宜梅花鹿扩散的坡度为10°,梅花鹿可能倾向于利用缓坡进行扩散。结果对东北虎豹国家公园梅花鹿种群的保护具有重要意义。     

Two Genetically Distinct Lineages of the Sika Deer, Cervus nippon, in Japanese Islands: Comparison of Mitochondrial D-Loop Region Sequences

To investigate genetic diversity among populations of the sika deer, Cervus nippon, nucleotide sequences (705–824 bases) of the mitochondrial D-loop regions were determined in animals from 13 localities in the Japanese islands. Phylogenetic trees constructed by the sequences indicated that the Japanese sika deer is separated into two distinct lineages: the northern Japan group (the Hokkaido island and most of the Honshu mainland) and the southern Japan group (a part of the southern Honshu mainland, the Kyushu island, and small islands around the Kyushu island). All sika deer examined in this study shared four to seven units of repetitive sequences (37 to 40 bases each) within the D-loop sequences. The number of tandem repeats was different among the populations, and it was specific to each population. Six or seven repeats occurred in populations of the northern Japan group, while four or five repeats occurred in populations of the southern Japan group. Each repeat unit included several nucleotide substitutions, compared with others, and 26 types were identified from 31 animals. Sequences of the first, second, and third units in arrays were clearly different between the northern and the southern groups. Based on these D-loop data, colonization and separation of the sika deer populations in the Japanese islands were estimated to have occurred less than 0.5 million years before present. Our results provide an invaluable insight into better understanding the evolutionary history, phylogeny, taxonomy, and population genetics of the sika deer.     

Landscape structure affects food quality of sika deer (<Emphasis Type="Italic">Cervus nippon</Emphasis>) evidenced by fecal nitrogen levels

Evaluating the quality of wildlife habitat is essential for understanding and predicting population dynamics in heterogeneous environments. We used fecal nitrogen levels as an indicator of habitat quality of sika deer (Cervus nippon) and explored important landscape elements influencing nitrogen levels, taking deer density into account. We established 92 plots differing in deer density and landscape structure on the Boso Peninsula, central Japan, and collected fecal samples along a 1-km transect at each plot. The regression models involving two independent variables, i.e., deer density and the length of forest edge within an area of 100 or 200 m from the transect, were selected based on the Akaike’s Information Criterion (AIC). Levels of fecal nitrogen were positively correlated with the length of the forest edge and negatively correlated with population density of deer. The area of 100 or 200 m from the transect most likely reflected the behavioral scale of the deer. Coverage of palatable understory vegetation increased with proximity to forest edge and decreased with deer density. Variability in the level of fecal nitrogen could thus be explained by food availability in the landscape. These results suggest that landscape alterations increase the carrying capacity of sika deer and thereby increase impacts upon the ecosystem.     

Genetic structure of,and hybridisation between,red (Cervus elaphus) and sika (Cervus nippon) deer in Ireland

This study investigated the levels of genetic diversity and variation exhibited by red and sika deer in Ireland, along with the extent and regional location of hybridisation between these two species. Bi-parental (microsatellites) and maternally-inherited (mitochondrial DNA) genetic markers were utilised that allowed comparisons between 85 red deer from six localities and 47 sika deer from 3 localities in Ireland. Population genetic structure was assessed using Bayesian analysis, indicating the existence of two genetic clusters in sika deer and three clusters in red deer. Levels of genetic diversity were low in both red and sika deer. These genetic data presented herein indicate a recent introduction of sika deer and subsequent translocations in agreement with historical data. The origins of the current red deer populations found in Ireland, based on genetic data presented in this study, still remain obscure. All hybrid deer (red/sika) found in this study were found in Wicklow, Galway and Mayo where the ‘red-like’ deer exhibited sika deer alleles/haplotypes, and vice versa in the case of Wicklow. Molecular methods proved invaluable in the identification of the hybrid deer because identification of hybrids based on phenotypic external appearances (pelage and body proportions) can be misleading. Areas where red and sika deer are sympatric need to be assessed for the level and extent of hybridisation occurring and thus need to be managed in order to protect the genetic integrity of ‘pure’ red deer populations.     

Genetic discontinuities in a continuously distributed and highly mobile ungulate,the Norwegian moose

Many species with currently continuously distributed populations have histories of geographic range shifts and successive shifts between decline or fragmentation, growth and spatial expansion. The moose (Alces alces) colonised Scandinavia after the last ice age. Historic records document a high abundance and a wide distribution across Norway in the middle ages, but major decline and fragmentation in the eighteenth and nineteenth centuries. After growth and expansion during the twentieth century, the Norwegian population is currently abundant and continuously distributed. We examined the distribution of genetic variation, differentiation and admixture in Norwegian moose, using 15 microsatellites. We assessed whether admixture has homogenised the population or if there are any genetic structures or discontinuities that can be related to recent or ancient shifts in demography or distribution. The Bayesian clustering algorithm STRUCTURE without any spatial information showed that there is currently a genetic dichotomy dividing the population into one southern and one northern subpopulation. Including spatial information, the Bayesian clustering algorithm TESS, which considers gradients of genetic variation and spatial autocorrelation, suggests that the population is divided into three subpopulations along a latitudinal axis, the southern one identical to the one identified with STRUCTURE. Present convergence zones of high admixture separate the identified subpopulations, which are delimited by genetic discontinuities corresponding to geographic barriers against dispersal, e.g. wide fiords and mountain ranges. The distribution of the subpopulations is supported by spatial autocorrelation analysis. However, some loci are not in Hardy–Weinberg equilibrium and the STRUCTURE analysis suggests that a lower hierarchical structure may exist within the southernmost subpopulation. No bottlenecks or founder events are indicated by the levels of genetic variation, rather a high degree of private alleles in the northern subpopulations indicates introgression. Coalescent-based Approximate Bayesian Computation estimates unambiguously suggest that the genetic structure is a result of an ancient divergence event and a more recent admixture event a few centuries ago. This indicates that the central Scandinavian subpopulation constitutes a relatively recent convergence zone of secondary contact.     

Genetic structure of red deer population in northeastern Poland in relation to the history of human interventions

We studied the genetic structure of a red deer (Cervus elaphus) population in 8 woodlands of northeastern Poland and 1 in western Belarus and compared it with the documented history of the population in the region. Red deer nearly went extinct in the region in the 18th and 19th centuries. In the mid-19th century, reintroductions began and continued until the mid-1960s. Animals were translocated from various sites in Poland and other European countries. We genotyped 303 individuals using 14 microsatellite loci and sequenced 253 individuals for a fragment of the control region (mitochondrial DNA [mtDNA]). The microsatellite analyses demonstrate that 3 genetically separate subpopulations exist, but 4 according to mtDNA. All haplotypes found in northeastern Poland are closely related to haplotypes from northern and northwestern Europe. The only individuals that could have originated from autochthonous red deer populations, rather than introductions, were found in Napiwoda Forest. The present regional genetic structure of the species is consistent with the known history of red deer translocations. Current patterns of genetic diversity in these populations are determined by the interaction of past human management and contemporary natural migrations. © 2012 The Wildlife Society.     

Influence of habitat fragmentation on population structure of red deer in Croatia

The genetic structure of red deer populations is under strong influence of human activities such as game management and habitat fragmentation. Using multilocus genotypes from 193 geo-referenced individuals, we evaluated the population genetic structure of three red deer populations in Croatia. The effect of habitat fragmentation on genetic structure was tested using Bayesian non-spatial and spatial clustering methods. Our results indicate levels of genetic diversity similar to the ones previously reported by other authors for stable and appropriately managed populations within all populations analyzed. The spatial clustering model was able to detect the effect of habitat fragmentation on population differentiation, supporting the use of spatially explicit methods in landscape genetics, and giving important guidelines for future road planning.     

Geographic differentiation of <Emphasis Type="Italic">Cobitis shikokuensis</Emphasis> inferred from mtDNA RFLP analysis

Mitochondrial DNA divergence among populations of the Japanese spinous loach Cobitis shikokuensis, endemic to Shikoku Island, was investigated by restricted fragment length polymorphism analysis. A total of 68 restriction sites on DNA fragments from the cytochrome b to D-loop regions and from the 12S rRNA to 16S rRNA regions, amplified by PCR, were analyzed. A total of 12 haplotypes (plus 6 in outgroups) were detected in 268 specimens collected from 19 localities in seven rivers (and 41 specimens from four localities in three rivers in outgroups). Three of the seven river populations of C. shikokuensis were shown to have unique haplotypes, and four of the seven river populations were monomorphic. The nested structure of the haplotype network for populations of C. shikokuensis exhibited two large clades corresponding to (1) populations from the Shimanto River and its neighbors and (2) two genetically divergent populations in the Shigenobu and Iwamatsu Rivers. The population from the Shimanto River, the largest river inhabited by C. shikokuensis, maintains great haplotype diversity as well as the allozyme diversity previously reported. On the other hand, populations from the Hiji River, the second largest river, which exhibited the highest allozyme diversity, were monomorphic in their mtDNA. The nested clade analysis (NCA) revealed that past fragmentation between the above two clades could occur in the initial distribution process of C. shikokuensis. The large genetic divergence of two river populations from the Shigenobu and Iwamatsu Rivers was inferred to be caused by a process of long distance colonization and fragmentation. MtDNA introgression into the Hiji River population from southern river populations was suggested. Taking genetic divergence into consideration, each river population of C. shikokuensis should be conserved separately as like a distinct species, and conservation programs for the small populations showing less genetic variability should be invoked as soon as possible.     

Bottleneck effects on the sika deer Cervus nippon population in Hokkaido, revealed by ancient DNA analysis

The population size of the sika deer Cervus nippon on Hokkaido Island of Japan had been remarkably reduced because of heavy hunting pressure since the beginning of Meiji Period and effects of heavy snow in 1879 and 1881. After that, the number of sika deer in Hokkaido has increased gradually due to the protection by the Hokkaido government. In the present study, in order to investigate the bottleneck effects, we analyzed ancient mitochondrial DNA (mtDNA) on sika deer bones excavated from archaeological sites just before Meiji Period. On 86 of 113 bones from 13 archaeological sites of Ainu Culture Period (17-19th centuries), 602 base-pair fragments of the mtDNA control region were successfully sequenced. Consequently, we found three new haplotypes (g-, h- and i-types) which had not been identified in modern sika deer. In addition, four haplotypes (a-, b-, c- and d-types) identified from modern sika deer were also found in the archaeological deer. The new haplotypes and previously reported hapoltypes from sika deer of Hokkaido were phylogenetically much closer to each other, compared with those of modern sika deer from Honshu, Kyushu and the Chinese continent. Geographical distribution patterns of haplotypes of the ancient population were different from those of the modern population in Hokkaido. Our findings indicated that their genetic diversity was reduced through the bottleneck and that population structures of sika deer were changed widely in Hokkaido due to genetic drift.     

Population genetic structure and the effect of founder events on the genetic variability of introduced sika deer, Cervus nippon, in Germany and Austria

The relationships among 214 wild-living sika deer from five locations in Germany and two in Lower Austria were examined using mitochondrial DNA (mtDNA) control region sequence. A total of 18 haplotypes are grouped consistently into two major divergent clades, A and B, which differ by a mean of 8.4% sequence divergence. Recently introduced sika deer showed a complex pattern of population structuring, which probably results from historical vicariance in at least two unknown source populations from southeastern Asia (as previously described by morphological and mtDNA findings), and subsequent population admixture as a result of human-mediated restocking. A strong genetic differentiation among populations was indicated by a global _ST value of 0.78 reflecting mainly the differential distribution of clades A and B haplotypes. There was no association between related haplotypes and their distribution among local populations. These indicate that genealogy is a better predictor of the genetic affinity among most sika deer populations than their present-day locations. The abundant mitochondrial divergence we observed, may reflect a subspecies differentiation and could be associated with phenotypic differences among the introduced sika deer.An erratum to this article can be found at     

Evaluation of introgressive hybridization among Cervidae in Japan's Kinki District via two novel genetic markers developed from public NGS data

Hybridization and backcrossing of native populations with introduced species can lead to introgression and genetic alteration. In this study, we evaluated introgression in 43 deer from a potential hybrid zone around Okinoshima Island, Kinki District, Japan. This region witnessed the migration of a hybrid population (cross between the Formosan sika deer [Cervus nippon taiouanus] and other deer species) that could potentially breed with the native Japanese sika deer (C. n. centralis). We used an existing genetic marker for the mitochondrial cytochrome b gene and two novel markers for nuclear DNA, developed using publicly available next‐generation sequencing data. We identified one mainland deer with a mitochondrial haplotype identical to that of the Formosan sika deer as well as nuclear heterozygous sequences identical to those of Formosan and Japanese sika deer. This suggests that the mainland deer is a hybrid offspring of the Okinoshima population and native deer. However, only Japanese sika deer sequences were found in the other 42 samples, indicating limited introgression. Nevertheless, hybridization pre‐ and postintroduction in the Okinoshima population could cause multispecies introgression among Japanese sika deer, negatively affecting genetic integrity. We developed a simple test based on polymerase chain reaction–restriction fragment length polymorphism to detect introgression in natural populations. Our method can accelerate genetic monitoring of Japanese sika deer in Kinki District. In conclusion, to prevent further introgression and maintain genetic integrity of Japanese sika deer, we recommend establishing fences around Okinoshima Island to limit migration, besides a continued genetic monitoring of the native deer.     

Use of microsatellite markers to assess the spatial genetic structure of a population of sika deerCervus nippon on Kinkazan Island, Japan

To understand the relationship between social behaviour and gene distribution, we used microsatellite markers to resolve the spatial genetic structure of the sika deerCervus nippon Temminck, 1838 population on Kinkazan Island, a small island (9.6 km²) in northern Japan. We obtained 177 samples of deer which correspond to about 30% of the total population on the island. 126 were from a local population where each deer was individually identified, while 51 from other area on the island. Although there were no apparent geographical barriers on the island, the sika deer population showed local differentiation in its genetic composition. By comparing allele-sharing rates between gender and social categories, we demonstrated a higher genetic relatedness within males in a local group. The Assignment Index (AI) of each sex within a local group showed similar distributions. None of our analyses indicated a lower genetic relatedness among males than females, which contradicted our predictions based on the dispersal behaviour of males. Considering other factors, the results suggest that the range of male dispersal is quite limited on this island.