A genetic linkage map for cattle

Individual identification of sika deer (Cervus nippon) using short tandem repeat analysis for investigating illegal carcass disposal in Japan

The population of Japanese sika deer (Cervus nippon) are controlled by hunting to prevent damage to various crops in many areas in Japan. Hunters are subsidized by submitting the tail to the local government; animal carcasses must be properly disposed of after the hunt, and abandonment of hunted deer in the field is prohibited by law. However, there have been many carcasses of sika deer being abandoned without proper disposal. In such cases, individual identification by DNA analysis is required to match the abandoned deer and submitted tail and identify the suspect. When identifying individual wildlife by DNA analysis, it is crucial to select appropriate markers that consider both the procedure of the analysis and the animal species. To evaluate availability of Short Tandem Repeat (STR) analysis for the identification of sika deer, this study aimed to construct an STR database for sika deer in Japan and to evaluate the discrimination power of STR markers, using an identification kit for a closely related species of cattle and STR markers of the sika deer. The results showed polymorphism at six STR loci from the Bovine Genotypes Panel 3.1 Kit and two STR markers for sika deer, suggesting that these loci may be useful for sika deer identification. The coincidence rate for the three STR loci (CSSM019, TGLA53, ETH10) was 7.63 × 10-⁴, which was considered sufficient for identification of the sika deer population. This study was the first to evaluate the availability of sika deer STR analysis for individual identification in Japan and was expected to have applications in crime scene and wildlife forensics.

Mountain gazelles' (Gazella gazella) males use mutual dung middens in favorable locations

Males of the mountain gazelle deposit dung middens (different colors and shapes represent middens of different haplotypes) in preferable forest plots and countermark the same middens (two color circles) at the boundaries of their territories.

Fertility and early embryonic development in a CD46-edited Gir heifer with reduced susceptibility to BVDV†

Bovine viral diarrhea virus (BVDV) infection during pregnancy is a significant contributor to reproductive failures in cattle. The bovine receptor for BVDV (CD46) was previously edited with a six amino acid substitution (G82QVLAL to A82LPTFS) and shown to have significantly reduced BVDV susceptibility in a Gir heifer calf. Since a role for CD46 has been proposed in mammalian fertilization, our objective was to assess the edited heifer's fertilization rates, early embryonic development, and germline transmission conformation of the edit. Cumulus oocyte complexes were collected from the edited heifer and unedited females, fertilized with semen from an unedited bull and cultured until the blastocyst stage. Ultrasound examinations and serum progesterone concentration were also monitored to confirm estrous cyclicity in the CD46-edited heifer. Estrous cyclicity was normal with visualization of a corpus luteum and elevated progesterone concentrations. Fertilization rates and blastocyst development were not different in oocytes from edited and unedited controls. Genome sequence analysis of blastocysts confirmed germline transmission of either edited allele from the heifer. Subsequently, the CD46-edited heifer was artificially inseminated with semen from an unedited Gir bull and fertility status was confirmed with a diagnosed conception at Day 35 of gestation. Thus, a six amino acid substitution in CD46 did not negatively affect fertilization of edited oocytes or early embryonic development when fertilized with semen from an unedited bull. An edited bull is still needed to similarly evaluate reproductive function of sperm cells carrying this CD46 edit.

Relatedness of White-Tailed Deer from Culling Efforts Within Chronic Wasting Disease Management Zones in Minnesota

In white-tailed deer (Odocoileus virginianus), closely related females form social groups, avoiding other social groups. Consequently, females infected with chronic wasting disease (CWD) are more likely to infect social group members. Culling has been used to reduce CWD transmission in high-risk areas; however, its effectiveness in removing related individuals has not been assessed. We analyzed 11 microsatellites and a mitochondrial DNA fragment to assess: (1) the genetic structure in white-tailed deer in Minnesota, USA and (2) the effectiveness of localized culling to remove related deer. For (1), we genotyped deer culled in 2019 and 2021 in three CWD management zones, and deer collected in between zones. For (2), we only included culled deer, defining "culled groups" as deer obtained in the same township-range-section and year. We compared mean relatedness among deer from the same culled group (intra-group relatedness) and among deer from different culled groups (inter-group relatedness). We did not find evidence of genetic structure, suggesting that an outbreak in any of the management zones could naturally spread to the others. Culling removed deer that were on average more related than expected by chance (intra-group relatedness > inter-group relatedness), and most highly-related deer were culled in the same bait site.

Prion forensics: a multidisciplinary approach to investigate CWD at an illegal deer carcass disposal site

Infectious prions are resistant to degradation and remain infectious in the environment for several years. Chronic wasting disease (CWD) has been detected in cervids inhabiting North America, the Nordic countries, and South Korea. CWD-prion spread is partially attributed to carcass transport and disposal. We employed a forensic approach to investigate an illegal carcass dump site connected with a CWD-positive herd. We integrated anatomic, genetic, and prion amplification methods to discover CWD-positive remains from six white-tailed deer (Odocoileus virginianus) and, using microsatellite markers, confirmed a portion originated from the CWD-infected herd. This approach provides a foundation for future studies of carcass prion transmission risk.

Detection of putative loci affecting milk yield in Turkish Awassi sheep using microsatellite markers

On the basis of comparisons between bovine and ovine genome mapping information, the aim of the study was to analyze the genetic diversity of selected DNA microsatellites from the bovine genome and to investigate their correlation with the average daily milk yield in Awassi sheep. 18 informative microsatellite markers were selected from the significant QTL regions affecting milk yield identified in the bovine genome in previous studies. The selected microsatellite markers were then amplified by PCR as reciprocal amplifications on the genomic DNA of Awassi sheep, with standard daily milk yield records. Thus, in this study, 18 microsatellite markers associated with milk yield in the bovine genome were examined for both determination of genetic polymorphism within the flock and the effects of marker loci on average daily milk yield in Awassi sheep. Allele frequencies of markers were determined based on the results of fragment analysis. The analysis of variance showed that the 123 bp allele at the marker locus BMS1341 on BTA2 significantly influenced the average daily milk yield of Ivesi sheep (P < 0.01). On the other hand, the BMS381 locus with a 115 bp allele on BTA2, the MCM140 locus with a 185 bp allele on BTA6, the BMS2721 locus with a 155 bp allele, the BM1237 locus with 174 and 180 bp alleles on BTA7, and finally, the BMS1967 locus with a 117 bp allele, the BM4208 locus with 176 and 182 bp alleles, and the INRA locus with a185 bp allele on BTA8 showed moderately significant effects on the average daily milk yield of Ivesi ewes (P < 0.05).

Genetic Analysis Reveals Dispersal Patterns of Japanese Serow in Two Different Habitats of a Mountainous Region

Dispersal increases the costs of feeding and predation risk in the new environment and is reported to be biased toward habitats similar to the natal region in some mammals. The benefits and costs of dispersal often differ between sexes, and most mammals show male-biased dispersal in relation to a polygamous mating system. Japanese serow is generally a solitary and monogamous species. However, recent studies have shown that the sociality of serows on Mt. Asama differs between habitat types. In the mountain forests with low forage availability, solitary habits and social monogamy were observed, while, in alpine grasslands, female grouping and social polygyny were observed, which is probably due to abundant forage availability. We investigated the effects of habitat characteristics and sociality on the dispersal of serows using fecal and tissue samples from two different habitats on Mt. Asama. The Fst value between the two areas was significantly positive, and the mean relatedness within areas was significantly higher than that between areas, which suggests limited gene flow and natal habitat-biased dispersal. Bayesian clustering analysis showed unidirectional gene flow from forest to grassland, which was probably due to the high forage availability of the grassland. Analyses of the assignment index and mean relatedness did not show male-biased dispersal, even in the grassland, where serows were polygynous. Thus, polygyny in the grassland is not linked to male-biased dispersal. In summary, our study suggests that dispersal patterns in Japanese serows are affected by habitat rather than social differences.

Hybridization in an isolated population of blesbok and red hartebeest

Hybridization in antelope species has been widely reported in South African national parks and provincial reserves as well as on private land due to anthropogenic effects. In a closed management setting, hybridization may occur due to the crossbreeding of closely related species with unequal sex ratios, resulting in either sterile or fertile offspring. In this study, we used molecular techniques to evaluate the risk of anthropogenic hybridization between blesbok (Damaliscus pygargus phillipsi) and red hartebeest (Alcelaphus buselaphus caama) in an isolated group that purposely included the two species with unequal sex ratios (one red hartebeest male and 19 male and female blesbok). Genetic analysis based on microsatellites confirmed the presence of seven hybrid individuals. Mitochondrial analysis verified that hybridization occurred between blesbok females and the red hartebeest male. STRUCTURE and NEWHYBRIDS classified the hybrids as F1. It is suspected that the hybrid individuals were sterile as the males had undeveloped testes and only F1 hybrids were detected. Thus, the risk of hybridization between these two species may be limited in the wild. In captive settings, genetic monitoring should be included in management plans for blesbok and red hartebeest to ensure that the long-term consequences of wasted reproductive effort are limited.

Evaluation of the Impact of Population Management on the Genetic Parameters of Selected Spiral-Horned Antelopes

The rapid loss of biodiversity and the associated reduction and fragmentation of habitats means that ex situ populations have become an important part of species conservation. These populations, which are often established from a small number of founders, require careful management to avoid the negative effects of genetic drift and inbreeding. Although the inclusion of molecular data is recommended, their availability for captive breeding management remains limited. The aim of this study was to evaluate the relationship between the levels of genetic diversity in six spiral-horned antelope taxa bred under human care and their respective management strategies, conservation status, demography, and geographic origin, using 10 nuclear DNA microsatellite loci and mitochondrial control region DNA sequences. Our findings include associations between genetic diversity and management intensity but also with the diversity and contribution of wild populations to captive founders, with some populations apparently composed of animals from divergent wild lineages elevating captive genetic diversity. When population sizes are large, the potential advantages of maximizing genetic diversity in widely outcrossed populations may need careful consideration with respect to the potential disruption of adaptive diversity. Genetic data serve as a robust tool for managing captive populations, yet their interpretation necessitates a comprehensive understanding of species biology and history.

Differentiation of Bos grunniens and Bos taurus based on STR locus polymorphism

Differentiation of closely related biological species using molecular genetic analysis is important for breeding farm animals, creating hybrid lines, maintaining the genetic purity of breeds, lines and layering. Bos grunniens and Bos taurus differentiation based on STR locus polymorphism will help maintain the genetic isolation of these species and identify hybrid individuals. The aim of this study is to assess the differentiating potential of 15 microsatellite loci to distinguish between domestic yak (B. grunniens) bred in the Kalmak-Ashuu highland region (Kochkor district, Naryn region, Kyrgyz Republic) and cattle (B. taurus) of three breeds (Aberdeen-Angus, Holstein and Alatau) using molecular genetic analysis. The samples were genotyped at 15 microsatellite loci (ETH3, INRA023, TGLA227, TGLA126, TGLA122, SPS115, ETH225, TGLA53, BM2113, BM1824, ETH10, BM1818, CSSM66, ILSTS006 and CSRM60). Twelve of the loci were from the standard markers panel recommended by ISAG. Statistical analysis was performed using GenAlEx v.6.503, Structure v.2.3.4, PAST v.4.03, and POPHELPER v1.0.10. The analysis of the samples' subpopulation structure using the Structure v.2.3.4 and 15 STR locus genotyping showed that the accuracy of assigning a sample to B. taurus was 99.6 ± 0.4 %, whereas the accuracy of assigning a sample to B. grunniens was 99.2 ± 2.6 %. Of the 15 STRs, the greatest potential to differentiate B. grunniens and B. taurus was found in those with the maximal calculated FST values, including BM1818 (0.056), BM1824 (0.041), BM2113 (0.030), CSSM66 (0.034) and ILSTS006 (0.063). The classification accuracy of B. grunniens using only these five microsatellite loci was 98.8 ± 3.4 %, similar for B. taurus, 99.1 ± 1.2 %. The proposed approach, based on the molecular genetic analysis of 5 STR loci, can be used as an express test in Kyrgyzstan breeding and reproduction programs for B. grunniens.

The first Brazilian bovine breed: structure and genetic diversity of the Curraleiro Pé-duro

Background

The production of animal-based foods from native breeds have a synergistic relationship with the regional culture, the local climate, and mainly the maintenance of alternative genetic resources for a system with a lower environmental impact. Thus the efficiency of conservation and production depends on assessing the variability of these local breeds. In the case of Curraleiro Pé-duro cattle, the most adapted individuals have undergone natural selection over five hundred years in the Brazilian savannas, mating with little or no human interference. The peculiarities of these biomes, where the regional flora is the food base and cattle is raised in extensive areas, likely influenced the genetic composition of the different groups that make up the first cattle breed of Brazil.

Methods

To evaluate the composition, diversity, variation, differentiation, and genetic structure of the populations studied, samples of hair follicles from 474 individuals of different animal categories (calves, yearlings, heifers, cows, and bulls) from three farms, defined as subpopulations "A", "B", and "C", were collected. The animals were genotyped for 17 microsatellite markers using a DNA sequencer. After verification of monomorphic alleles, alleles outside the expected size range, and for the presence of stutter bands, the results were subjected to statistical analysis.

Results

The markers used were suitable for the proposed application with a mean Polymorphism Information Content (PIC) of 0.62. On average, the effective alleles were 4.25 per marker, with mean heterozygosities of 0.74 (observed and expected), which was lower in herd A (0.70) in comparison to herds B (0.77) and C (0.74). The analysis of molecular variance (AMOVA) revealed a higher rate of variation within herds (98.5%) and lower among herds (1.5%) (FSTranging from 0.00723 and 0.03198; p-values < 0.05). However no significant differences among herds where found with the Mantel test based on geographic distances. The formation of genetic clusters of all animals sampled with the software Structure resulted in minimum cluster values, with two main genetic groups (K = 2) observed among the evaluated animals. Therefore, based on PIC and heterozygosity values, a wide genetic diversity was observed, despite little differences in population structure (AMOVA, FST, and Structure results) among sampling sites.

Genetic Diversity of Wisent Bison bonasus Based on STR Loci Analyzed in a Large Set of Samples

Wisent Bison bonasus is an example of a species saved from extinction and reintroduced into nature after a few decades of captive breeding. There were only twelve founders of the Lowland–Caucasian line (LC) and even fewer (seven out of twelve) of Lowland (LB) animals. The genetic diversity in studies based on pedigree or markers is very low. In this paper, we present a summary of the long-term genetic monitoring conducted for the worldwide population of European bison. We summarized the long-term genetic monitoring studies conducted on the worldwide population of wisents to date. We genotyped 2227 wisents from two genetic lines (LC and LB) and different populations at ten microsatellite loci. We found low polymorphism, with only 2.7 alleles per locus, and much lower values of observed heterozygosity (0.380 and 0.348 in the LC and LB lines, respectively) than expected heterozygosity. The difference between the lines is only noticeable in allele proportions, so the number of markers is not enough to distinguish the two genetic lines. We also present the genetic distance among four free-roaming populations that are geographically close to each other. We found that the genetic distance of one of them is larger than that of the others, which could be the effect of genetic drift.

A historic religious sanctuary may have preserved ancestral genetics of Japanese sika deer (Cervus nippon)

Deer have been a major resource for human populations for thousands of years. Anthropogenic activities, such as hunting, have influenced the genetic structure and distribution of deer populations. In Japan, wild Japanese sika deer (Cervus nippon) have been hunted since ancient times but have also been historically protected as sacred animals in several sanctuaries. Sika deer have been protected for over a thousand years in the religious sanctuary around the Kasuga Taisha Shrine on the Kii Peninsula, located in the center of Japan. Here, we used short sequence repeats (SSR) and mitochondrial DNA (mtDNA) to investigate the genetic diversity, population structure, and demography of Japanese sika deer inhabiting the Kii Peninsula, Japan, and discuss possible anthropogenic influences. Using SSR, three distinct genetic groups were distinguished on the Kii Peninsula: an Eastern genetic group, a Western genetic group, and an isolated genetic group with individuals in the religious sanctuary of Kasuga Taisha Shrine in Nara city. The isolated genetic sanctuary group had only the mtDNA haplotype S4. The SSR genotype data suggested a newer divergence time of the genetic groups of the religious sanctuary than would have occurred as a result of Late Quaternary climate change. This time scale coincided with the establishment of the sanctuary with Kasuga Taisha Shrine. Thus, the religious protection conserved genetic variation over a thousand years.

Population Structure of Curraleiro Pé-Duro Cattle and its Relationship With the Serological Profile Against Pathogens of Economic and Zoonotic Interest

Curraleiro Pé-Duro (CPD) cattle, a Brazilian local breed, are animals that are highly resistant to infectious and parasitic diseases. Strategies for the conservation of the breed and the genetic resistance to diseases should also consider the characteristics of the breed. The objective of this study was to analyze the diversity and population structure of the CPD breed using microsatellite markers and to correlate the serological profiles for causative agents of brucellosis, leptospirosis, neosporosis, leukosis, infectious bovine rhinotracheitis, and bovine viral diarrhea. DNA samples of 608 bovines were amplified and genotyped using 28 microsatellite markers for breed characterization. The genotypes were assigned to three clusters, indicating a substructure population related to the geographic distance. The observed heterozygosity was lower than that expected in most loci, and fixation index (Fst) in the subpopulation values ranged from 0.03073 (ETH3) to 0.402 (INRA63) on cluster 1, from 0.00 (INRA32) to 0.39359 (INRA63) on cluster 2, and from 0.00 (INRA32) to 0.28483 (TGLA53) on cluster 3. The Brucella spp., Neospora caninum, and Bovine Leukemia Virus frequencies were significantly different (p < 0.05) between clusters. The differences between the occurrences of diseases were not sufficient to indicate a genetic subpopulation with increased resistance to infections.

De novo chromosome-length assembly of the mule deer (Odocoileus hemionus) genome

The mule deer (Odocoileus hemionus) is an ungulate species that is distributed in a range from western Canada to central Mexico. Mule deer are an essential source of food for many predators, are relatively abundant, and commonly make broad migration movements. A clearer understanding of the mule deer genome can improve our knowledge of its population genetics, movements, and demographic history, aiding in conservation efforts. Their large population size, continuous distribution, and diversity of habitat make mule deer excellent candidates for population genomics studies; however, few genomic resources are currently available for this species. Here, we sequence and assemble the mule deer genome into a highly contiguous chromosome-length assembly for use in future research using long-read sequencing and Hi-C technologies. We also provide a genome annotation and compare demographic histories of the mule deer and white-tailed deer using the pairwise sequentially Markovian coalescent model. We expect this assembly to be a valuable resource in the continued study and conservation of mule deer.

Spatial genetic structure and Ovis haplogroup as a tool for an adaptive conservation management of the endangered Cyprus mouflon

The mouflon (Ovis gmelini ophion) is the largest wild terrestrial mammal of Cyprus. Regarded as the local flagship species, its population range is limited to the mountainous Paphos Forest and adjacent areas including part of Troodos National Forest Park (> 700 Km²). This species, protected by both national and international legislation, is classified as Endangered by the IUCN, with livestock intrusion, poaching and habitat loss as the main threatening factors. We determined the spatial genetic structure by investigating sexes separately and shed further light on the Cyprus mouflon placement among Ovis haplogroups (HGs) to enforce its protection within an adaptive conservation framework. We genotyped 108 samples collected across the entire species' range at both14 loci of the microsatellite DNA and the mitochondrial DNA Control Region (CR). Microsatellite genotyping indicated that the wild population retained low genetic variability, which, however, was not associated with a level of inbreeding raising particular concern (F_IS = 0.12). An overall weak spatial genetic structure was disclosed, consistently with the limited extension of the mouflon range, the lack of significant physical barriers to dispersal and the intense gene flow mostly occurring along a northeast-southwest axis across the forest. Nevertheless, evidences of structure were found for females in compliance with their philopatric behaviour. The disclosure of unique features compared to Mediterranean and Near East conspecifics, such as the occurrence of a six-fold 76 bp-long repeated motif in the Cyprus mouflon CR, along with the outcome of a phylogenetic reconstruction (based on a far higher number of Ovis GenBank records and Cyprus haplotypes than in previous studies) inspired the proposal of a new haplogroup (HGC2) which included the Anatolian mouflon (O. g. anatolica, former HGX) as sister taxon. While both high habitat heterogeneity and low human disturbance were the main drivers in determining the overall spatial structure, future conservation efforts to preserve this valuable genetic resource should focus on avoiding possible introgressive hybridisation with co-grazing livestock to the edge of its range.

Vocal phenotype of male rutting roars and genetic markers delineate East European red deer (Cervus elaphus) from Central and West European populations

This study investigates a population of red deer Cervus elaphus, founded by 10 individuals introduced in the nineteenth century from Germany to the Voronezh region of the European part of Southern Russia and then developed without further introductions. We characterize for the first time the vocal phenotype of the Voronezh red deer male rutting calls in comparison with similar data on the Pannonian (native Central European) and Iberian (native West European) red deer obtained by the authors during preceding studies. In addition, we provide for the first time the genetic data on Pannonian red deer. In Voronezh stags, the number of roars per bout (2.85 ± 1.79) was lower than in Pannonian (3.18 ± 2.17) but higher than in Iberian (2.11 ± 1.71) stags. In Voronezh stags, the duration of main (the longest within bouts) roars was longer (2.46 ± 1.14 s) than in Pannonian (1.13 ± 0.50 s) or Iberian (1.90 ± 0.50 s) stags. The maximum fundamental frequency of main roars was similar between Voronezh (175 ± 60 Hz) and Pannonian (168 ± 61 Hz) but higher in Iberian stags (223 ± 35 Hz). Mitochondrial cytochrome b gene analysis of red deer from the three study populations partially supports the bioacoustical data, of closer similarity between Voronezh and Pannonian populations. In contrast, microsatellite DNA analysis delineates Voronezh red deer from either Pannonian or Iberian red deer. We discuss that population bottlenecking might affect the acoustics of the rutting roars, in addition to genotype.

A Microsatellite Genotyping-Based Genetic Study of Interspecific Hybridization between the Red and Sika Deer in the Western Czech Republic

Simple Summary

The sika deer is a very flexible invasive species, capable of living dynamically in both large forests and mixed environment characterized by a prevalence of agricultural land. The Japanese sika deer was introduced to the Czech Republic at the end of the 19th century. The success of an introduced alien species may consist in their hybridizing with closely related taxa. Where few barriers to gene flow exist, rapid introgression of genetic traits from one species into another frequently occurs. The current Czech sika populations embody the most abundant and expanding group in continental Europe. In western Bohemia, we confirmed the interspecific hybridization with the native red deer. In this context, the red deer gene pool is endangered. The animals proliferate steadily in all directions and will most probably spread all over the Czech Republic if no major, timely changes in game management are adopted.

Abstract

Although inter-species hybrids between the red and sika deer can be phenotypically determined only exceptionally, there is the eventuality of identification via molecular genetic analysis. We used bi-parentally inherited microsatellite markers and a Bayesian statistical framework to re-examine the proportion of hybrids in the Czech red and sika deer populations. In total, 123 samples were collected, and the nuclear dataset consisted of 2668 allelic values. The number of alleles per locus ranged from 10 (BM1818) to 22 (BM888 and T193), yielding the mean of 16 alleles per locus across the deer. The mean allelic diversity of the red deer markedly exceeded that of the Japanese sika deer. Interspecific hybrids were detected, enabling us to confirm the genetic introgression of the sika deer into the red deer populations and vice versa in western Bohemia. The mean hybrid score equaled 10.6%, with 14.3% of the hybrids being among red deer–like individuals and 6.7% among sika-like ones. At two western Bohemian locations, namely, Doupovské hory and Slavkovský les, the total percentages of hybrid animals equaled 18.8 and 8.9, respectively. No red deer alleles were detected in the sika populations of the subregions of Kladská, Žlutice, and Lány. The NeighborNet network clearly separated the seven red and sika deer sampling populations according to the geography. The knowledge gained from the evaluated data is applicable in hunting management to reduce hybridization with the European deer.

Development and implementation of a STR based forensic typing system for moose (Alces alces)

The moose is a highly prized game species in North America and is often the target of illegal harvesting. Forensic DNA analysis can be used to assist in investigations into wildlife crime by providing the link between the illegal incidents and suspects. In this study, we present the development and validation of a short-tandem repeat (STR) based forensic DNA test for the individualization of moose in Alberta and Yukon, Canada. We show that the markers used in the test are appropriate for forensic use and are sufficiently specific to moose. We also demonstrate the limit of detection and quantitation of the moose STR test on the ABI 3500 genetic analyzer. Lastly, we describe the population genetic structure of moose present in our forensic population database and make recommendations regarding the calculation of appropriately conservative forensic statistics.

Spatial familial networks to infer demographic structure of wild populations

In social species, reproductive success and rates of dispersal vary among individuals resulting in spatially structured populations. Network analyses of familial relationships may provide insights on how these parameters influence population-level demographic patterns. These methods, however, have rarely been applied to genetically derived pedigree data from wild populations.Here, we use parent-offspring relationships to construct familial networks from polygamous boreal woodland caribou (Rangifer tarandus caribou) in Saskatchewan, Canada, to inform recovery efforts. We collected samples from 933 individuals at 15 variable microsatellite loci along with caribou-specific primers for sex identification. Using network measures, we assess the contribution of individual caribou to the population with several centrality measures and then determine which measures are best suited to inform on the population demographic structure. We investigate the centrality of individuals from eighteen different local areas, along with the entire population.We found substantial differences in centrality of individuals in different local areas, that in turn contributed differently to the full network, highlighting the importance of analyzing networks at different scales. The full network revealed that boreal caribou in Saskatchewan form a complex, interconnected familial network, as the removal of edges with high betweenness did not result in distinct subgroups. Alpha, betweenness, and eccentricity centrality were the most informative measures to characterize the population demographic structure and for spatially identifying areas of highest fitness levels and family cohesion across the range. We found varied levels of dispersal, fitness, and cohesion in family groups. Synthesis and applications: Our results demonstrate the value of different network measures in assessing genetically derived familial networks. The spatial application of the familial networks identified individuals presenting different fitness levels, short- and long-distance dispersing ability across the range in support of population monitoring and recovery efforts.

Phylogeography and population genetic structure of red muntjacs: evidence of enigmatic Himalayan red muntjac from India

BACKGROUND

Identifying factors shaping population genetic structure across continuous landscapes in the context of biogeographic boundaries for lineage diversification has been a challenging goal. The red muntjacs cover a wide range across multiple vegetation types, making the group an excellent model to study South and Southeast Asian biogeography. Therefore, we analysed mitogenomes and microsatellite loci, confirming the number of red muntjac lineages from India, gaining insights into the evolutionary history and phylogeography of red muntjacs.

RESULTS

Our results indicated the Northwestern population of red muntjac or the Himalayan red muntjac (M. aureus) in India as genetically diverse and well-structured, with significant genetic differentiation implying a low level of gene flow. The phylogenetic, population genetic structure, as well as species delimitation analyses, confirm the presence of the lineage from Western Himalayan in addition to the previously identified red muntjac lineages. Relatively low genetic diversity was observed in M. aureus compared to M. vaginalis, M. malabaricus and M. muntjak. The M. aureus and M. vaginalis lineages have split during the late Pleistocene, ~ 1.01 million years ago (Mya), making M. aureus the youngest lineage; whereas, M. malabaricus split earlier, ~ 2.2 Mya and appeared as the oldest lineage among red muntjacs.

CONCLUSIONS

Pronounced climate fluctuations during the Quaternary period were pivotal in influencing the current spatial distribution of forest-dwelling species' restriction to Northwestern India. Our finding confirms the distinct Himalayan red muntjac (M. aureus) within the red muntjac group from Northwestern India that should be managed as an Evolutionary Significant Unit (ESU). We recommend a reassessment of the conservation status of red muntjacs for effective conservation and management.

Population genetics and evolutionary history of the endangered Eld's deer (Rucervus eldii) with implications for planning species recovery

Eld's deer (Rucervus eldii) with three recognised subspecies (R. e. eldii, R. e. thamin, and R. e. siamensis) represents one of the most threatened cervids found in Southeast Asia. The species has experienced considerable range contractions and local extinctions owing to habitat loss and fragmentation, hunting, and illegal trade across its distribution range over the last century. Understanding the patterns of genetic variation is crucial for planning effective conservation strategies. This study investigated the phylogeography, divergence events and systematics of Eld's deer subspecies using the largest mtDNA dataset compiled to date. We also analysed the genetic structure and demographic history of R. e. eldii using 19 microsatellite markers. Our results showed that R. e. siamensis exhibits two divergent mtDNA lineages (mainland and Hainan Island), which diverged around 0.2 Mya (95% HPD 0.1-0.2), possibly driven by the fluctuating sea levels of the Early Holocene period. The divergence between R. e. eldii and R. e. siamensis occurred around 0.4 Mya (95% HPD 0.3-0.5), potentially associated with the adaptations to warm and humid climate with open grassland vegetation that predominated the region. Furthermore, R. e. eldii exhibits low levels of genetic diversity and small contemporary effective population size (median = 7, 4.7-10.8 at 95% CI) with widespread historical genetic bottlenecks which accentuates its vulnerability to inbreeding and extinction. Based on the observed significant evolutionary and systematic distance between Eld's deer and other species of the genus Rucervus, we propose to classify Eld's deer (Cervus eldii) in the genus Cervus, which is in congruent with previous phylogenetic studies. This study provides important conservation implications required to direct the ongoing population recovery programs and planning future conservation strategies.

Male recombination map of the autosomal genome in German Holstein

Background

Recombination is a process by which chromosomes are broken and recombine to generate new combinations of alleles, therefore playing a major role in shaping genome variation. Recombination frequencies (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ) between markers are used to construct genetic maps, which have important implications in genomic studies. Here, we report a recombination map for 44,696 autosomal single nucleotide polymorphisms (SNPs) according to the coordinates of the most recent bovine reference assembly. The recombination frequencies were estimated across 876 half-sib families with a minimum number of 39 and maximum number of 4236 progeny, comprising over 367 K genotyped German Holstein animals.

Results

Genome-wide, over 8.9 million paternal recombination events were identified by investigating adjacent markers. The recombination map spans 24.43 Morgan (M) for a chromosomal length of 2486 Mbp and an average of ~ 0.98 cM/Mbp, which concords with the available pedigree-based linkage maps. Furthermore, we identified 971 putative recombination hotspot intervals (defined as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ > 2.5 standard deviations greater than the mean). The hotspot regions were non-uniformly distributed as sharp and narrow peaks, corresponding to ~ 5.8% of the recombination that has taken place in only ~ 2.4% of the genome. We verified genetic map length by applying a likelihood-based approach for the estimation of recombination rate between all intra-chromosomal marker pairs. This resulted in a longer autosomal genetic length for male cattle (25.35 cM) and in the localization of 51 putatively misplaced SNPs in the genome assembly.

Conclusions

Given the fact that this map is built on the coordinates of the ARS-UCD1.2 assembly, our results provide the most updated genetic map yet available for the cattle genome.

An assessment of sampling designs using SCR analyses to estimate abundance of boreal caribou

Accurately estimating abundance is a critical component of monitoring and recovery of rare and elusive species. Spatial capture-recapture (SCR) models are an increasingly popular method for robust estimation of ecological parameters. We provide an analytical framework to assess results from empirical studies to inform SCR sampling design, using both simulated and empirical data from noninvasive genetic sampling of seven boreal caribou populations (Rangifer tarandus caribou), which varied in range size and estimated population density. We use simulated population data with varying levels of clustered distributions to quantify the impact of nonindependence of detections on density estimates, and empirical datasets to explore the influence of varied sampling intensity on the relative bias and precision of density estimates. Simulations revealed that clustered distributions of detections did not significantly impact relative bias or precision of density estimates. The genotyping success rate of our empirical dataset (n = 7,210 samples) was 95.1%, and 1,755 unique individuals were identified. Analysis of the empirical data indicated that reduced sampling intensity had a greater impact on density estimates in smaller ranges. The number of captures and spatial recaptures was strongly correlated with precision, but not absolute relative bias. The best sampling designs did not differ with estimated population density but differed between large and small ranges. We provide an efficient framework implemented in R to estimate the detection parameters required when designing SCR studies. The framework can be used when designing a monitoring program to minimize effort and cost while maximizing effectiveness, which is critical for informing wildlife management and conservation.

Genetic Diversity of Historical and Modern Populations of Russian Cattle Breeds Revealed by Microsatellite Analysis

Analysis of ancient and historical DNA has great potential to trace the genetic diversity of local cattle populations during their centuries-long development. Forty-nine specimens representing five cattle breeds (Kholmogor, Yaroslavl, Great Russian, Novgorod, and Holland), dated from the end of the 19th century to the first half of the 20th century, were genotyped for nine polymorphic microsatellite loci. Using a multiple-tube approach, we determined the consensus genotypes of all samples/loci analysed. Amplification errors, including allelic drop-out (ADO) and false alleles (FA), occurred with an average frequency of 2.35% and 0.79%, respectively. A significant effect of allelic length on ADO rate (r² = 0.620, p = 0.05) was shown. We did not observe significant differences in genetic diversity among historical samples and modern representatives of Kholmogor and Yaroslavl breeds. The unbiased expected heterozygosity values were 0.726–0.774 and 0.708–0.739; the allelic richness values were 2.716–2.893 and 2.661–2.758 for the historical and modern samples, respectively. Analyses of F_ST and Jost’s D genetic distances, and the results of STRUCTURE clustering, showed the maintenance of a part of historical components in the modern populations of Kholmogor and Yaroslavl cattle. Our study contributes to the conservation of biodiversity in the local Russian genetic resources of cattle.

Utilization of microsatellite markers in genotyping of Saudi Arabian camels for productivity and conservation

Camels are considered as integral and notable components of the heritage of the Kingdom of Saudi Arabia. Genetic variabilities within and among four camel populations in Saudi Arabia were evaluated using 21 simple sequence repeat (SSR) loci of 122 unrelated individuals, including three indigenous breeds [Humur (HA), Zurg (ZR), Shuguh (SG)] and one exotic breed [Sudanese (SN)]. Nineteen SSR markers generated multilocus fingerprints with a total of 225 alleles, a range of 4–23 alleles per locus, and an average of 9, 7, 7, and 6 alleles per locus in HA, ZR, SG, and SN populations, respectively. The mean multilocus FSTvalue (0.034 ± 0.005) showed non-significant population differentiation. Mean observed heterozygosity values were 0.908 for HA, 0.860 for ZR, 0.919 for SG, and 0.887 for SN, which were higher than the expected heterozygosity. An excess of heterozygotes was observed, suggesting the presence of overdominant selection or the occurrence of outbreeding. Pairwise genetic distances indicated that the three indigenous camel breeds were genetically close to each other and genetically distant to the SN population. This genetic variability assessment by microsatellite analysis is important and useful for the conservation of local camel genetic resources as well as the future development of breeding programs.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

Symposium review: Advances in sequencing technology herald a new frontier in cattle genomics and genome-enabled selection

The cattle reference genome assembly has underpinned major innovations in beef and dairy genetics through genome-enabled selection, including removal of deleterious recessive variants and selection for favorable alleles affecting quantitative production traits. The initial reference assemblies, up to and including UMD3.1 and Btau4.1, were based on a combination of clone-by-clone sequencing of bacterial artificial chromosome clones generated from blood DNA of a Hereford bull and whole-genome shotgun sequencing of blood DNA from his inbred daughter/granddaughter named L1 Dominette 01449 (Dominette). The approach introduced assembly gaps, misassemblies, and errors, and it limited the ability to assemble regions that undergo rearrangement in blood cells, such as immune gene clusters. Nonetheless, the reference supported the creation of genotyping tools and provided a basis for many studies of gene expression. Recently, long-read sequencing technologies have emerged that facilitated a re-assembly of the reference genome, using lung tissue from Dominette to resolve many of the problems and providing a bridge to place historical studies in common context. The new reference, ARS-UCD1.2, successfully assembled germline immune gene clusters and improved overall continuity (i.e., reduction of gaps and inversions) by over 250-fold. This reference properly places nearly all of the legacy genetic markers used for over a decade in the industry. In this review, we discuss the improvements made to the cattle reference; remaining issues present in the assembly; tools developed to support genome-based studies in beef and dairy cattle; and the emergence of newer genome assembly methods that are producing even higher-quality assemblies for other breeds of cattle at a fraction of the cost. The new frontier for cattle genomics research will likely include a transition from the individual Hereford reference genome, to a "pan-genome" reference, representing all the DNA segments existing in commonly used cattle breeds, bringing the cattle reference into line with the current direction of human genome research.

Genetic diversity and population genetic structure of six dromedary camel (camelus dromedarius) populations in Saudi Arabia

Camels are an integral and essential component of the Saudi Arabian heritage. The genetic diversity and population genetic structure of dromedary camels are poorly documented in Saudi Arabia so this study was carried out to investigate the genetic diversity of both local and exotic camel breeds. The genetic diversity was evaluated within and among camel populations using 21 microsatellite loci. Hair and blood samples were collected from 296 unrelated animals representing 4 different local breeds, namely Majaheem (MG), Maghateer (MJ), Sofr (SO), and Shaul (SH), and two exotic breeds namely Sawahli (SL) and Somali (SU). Nineteen out of 21 microsatellite loci generated multi-locus fingerprints for the studied camel individuals, with an average of 13.3 alleles per locus. Based on the genetic analyses, the camels were divided into two groups: one contained the Saudi indigenous populations (MG, MJ, SH and SO) and the other contained the non-Saudi ones (SU and SL). There was very little gene flow occurring between the two groups. The African origin of SU and SL breeds may explain their close genetic relationship. It is anticipated that the genetic diversity assessment is important to preserve local camel genetic resources and develop future breeding programs to improve camel productivity.

Phylogeography of moose in western North America

Subspecies designations within temperate species’ ranges often reflect populations that were isolated by past continental glaciation, and glacial vicariance is believed to be a primary mechanism behind the diversification of several subspecies of North American cervids. We used genetics and the fossil record to study the phylogeography of three moose subspecies (Alces alces andersoni, A. a. gigas, and A. a. shirasi) in western North America. We sequenced the complete mitochondrial genome (16,341 base pairs; n = 60 moose) and genotyped 13 nuclear microsatellites (n = 253) to evaluate genetic variation among moose samples. We also reviewed the fossil record for detections of all North American cervids to comparatively assess the evidence for the existence of a southern refugial population of moose corresponding to A. a. shirasi during the last glacial maximum of the Pleistocene. Analysis of mtDNA molecular variance did not support distinct clades of moose corresponding to currently recognized subspecies, and mitogenomic haplotype phylogenies did not consistently distinguish individuals according to subspecies groupings. Analysis of population structure using microsatellite loci showed support for two to five clusters of moose, including the consistent distinction of a southern group of moose within the range of A. a. shirasi. We hypothesize that these microsatellite results reflect recent, not deep, divergence and may be confounded by a significant effect of geographic distance on gene flow across the region. Review of the fossil record showed no evidence of moose south of the Wisconsin ice age glaciers ≥ 15,000 years ago. We encourage the integration of our results with complementary analyses of phenotype data, such as morphometrics, originally used to delineate moose subspecies, for further evaluation of subspecies designations for North American moose.

Evidence for adaptive introgression of exons across a hybrid swarm in deer

Background

Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes.

Results

We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat.

Conclusions

Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1497-x) contains supplementary material, which is available to authorized users.

Assessing introgressive hybridization in roan antelope (Hippotragus equinus): Lessons from South Africa

Biological diversity is being lost at unprecedented rates, with genetic admixture and introgression presenting major threats to biodiversity. Our ability to accurately identify introgression is critical to manage species, obtain insights into evolutionary processes, and ultimately contribute to the Aichi Targets developed under the Convention on Biological Diversity. The current study concerns roan antelope, the second largest antelope in Africa. Despite their large size, these antelope are sensitive to habitat disturbance and interspecific competition, leading to the species being listed as Least Concern but with decreasing population trends, and as extinct over parts of its range. Molecular research identified the presence of two evolutionary significant units across their sub-Saharan range, corresponding to a West African lineage and a second larger group which includes animals from East, Central and Southern Africa. Within South Africa, one of the remaining bastions with increasing population sizes, there are a number of West African roan antelope populations on private farms, and concerns are that these animals hybridize with roan that naturally occur in the southern African region. We used a suite of 27 microsatellite markers to conduct admixture analysis. Our results indicate evidence of hybridization, with our developed tests using a simulated dataset being able to accurately identify F1, F2 and non-admixed individuals at threshold values of qi > 0.80 and qi > 0.85. However, further backcrosses were not always detectable with backcrossed-Western roan individuals (46.7–60%), backcrossed-East, Central and Southern African roan individuals (28.3–45%) and double backcrossed (83.3–98.3%) being incorrectly classified as non-admixed. Our study is the first to confirm ongoing hybridization in this within this iconic African antelope, and we provide recommendations for the future conservation and management of this species.

Identification of muscle-specific candidate genes in Simmental beef cattle using imputed next generation sequencing

Genome-wide association studies (GWAS) have commonly been used to identify candidate genes that control economically important traits in livestock. Our objective was to detect potential candidate genes associated mainly with muscle development traits related to dimension of hindquarter in cattle. A next generation sequencing (NGS) dataset to imputed to 12 million single nucleotide polymorphisms (SNPs) (from 1252 Simmental beef cattle) were used to search for genes affecting hindquarter traits using a linear, mixed model approach. We also used haplotype and linkage disequilibrium blocks to further support our identifications. We identified 202 significant SNPs in the bovine BTA4 chromosome region associated with width of hind leg, based on a stringent statistical threshold (p = 0.05/ effective number of SNPs identified). After exploring the region around these SNPs, we found candidate genes that were potentially related to the associated markers. More importantly, we identified a region of approximately 280 Kb on the BTA4 chromosome that harbored several muscle specific candidate genes, genes to be in a potential region for muscle development. However, we also found candidate gene SLC13A1 on BTA4, which seems to be associated with bone disorders (such as chondrodysplasia) in Simmental beef cattle.

Genetic diversity and demographic history of introduced sika deer on the Delmarva Peninsula

The introduction of non-native species can have long-term effects on native plant and animal communities. Introduced populations are occasionally not well understood and offer opportunities to evaluate changes in genetic structure through time and major population changes such as bottleneck and or founder events. Invasive species can often evolve rapidly in new and novel environments, which could be essential to their long-term success. Sika deer are native to East Asia, and their introduction and establishment to the Delmarva Peninsula, USA, is poorly documented, but probably involved ≥1 founder and/or bottleneck events. We quantified neutral genetic diversity in the introduced population and compared genetic differentiation and diversity to the presumed source population from Yakushima Island, Japan, and a captive population of sika deer in Harrington, Delaware, USA. Based on the data from 10 microsatellite DNA loci, we observed reduced genetic variation attributable to founder events, support for historic hybridization events, and evidence that the population did originate from Yakushima Island stocks. Estimates of population structure through Bayesian clustering and demographic history derived from approximate Bayesian computation (ABC), were consistent with the hypothesized founder history of the introduced population in both timing and effective population size (approximately five effective breeding individuals, an estimated 36 generations ago). Our ABC results further supported a single introduction into the wild happening before sika deer spread throughout the Delmarva. We conclude that free-ranging sika deer on Delmarva are descended from ca. five individuals introduced about 100 years ago from captive stocks of deer maintained in the United Kingdom. Free-ranging sika deer on Delmarva have lost neutral diversity due to founder and bottleneck events, yet populations have expanded in recent decades and show no evidence of abnormalities associated with inbreeding. We suggest management practices including increasing harvest areas and specifically managing sika deer outside of Maryland.

Genetic Characteristics of the Japanese Serow Capricornis crispus in the Kii Mountain Range, Central Japan

The Japanese serow, Capricornis crispus, is an indigenous bovid species exclusively inhabiting mountain regions in the main Japanese islands, excepting Hokkaido. It had decreased in abundance to its lowest level due to overhunting and deforestation, with its distribution severely fragmented from the middle of the 20th century, many populations of C. crispus currently facing the risk of extinction. The Kii Mountain Range (KM) on Honshu is one such location that has seen a drastic population decline of C. crispus. In this study, we examined genetic characteristics of C. crispus in KM and neighboring regions of the Chubu district, using mtDNA and microsatellite markers, in order to devise strategies for its conservation. Results for mtDNA were characterized by low nucleotide diversity with five endemic and two dominant haplotypes shared by individuals in neighboring regions. A Bayesian skyline plot indicated a gradual increase after the last glacial maximum. For microsatellites, the genetic diversity of C. crispus in KM was comparable to Shizuoka and higher than Shikoku. Recent genetic bottlenecks were strongly suggested in C. crispus in KM. Bayesian clustering showed a genetic cline between KM and neighboring regions, where multivariate analysis suggested three local populations. A Mantel test indicated male-biased dispersal. These results indicate that C. crispus in KM and neighboring regions constitute multiple local populations, connected through restricted gene flow. For the conservation of C. crispus, it is important to define small-scale conservation units, among which genetic connectivity should be facilitated to prevent further loss of genetic diversity.

Population Structure, Admixture, and Migration Patterns of Japanese Sika Deer (Cervus nippon) Inhabiting Toyama Prefecture in Japan

Rapid expansion of sika deer, in both number and distribution, in the Japanese Archipelago has resulted in serious ecological disturbance. In the present study, the population structure and migration patterns of sika deer (Cervus nippon) among Toyama and adjacent Prefectures were investigated using 11 polymorphic microsatellite loci. Deviation from Hardy-Weinberg equilibrium was detected in both total and individual regional sika deer samples from Toyama Prefecture. Results of pairwise F_ST results, factorial correspondence analysis, and STRUCTURE analysis indicated that sika deer in Toyama are not genetically distinct from those in adjacent Prefectures. Bayesian STRUCTURE results suggested the existence of two distinct clusters. However, multiple lines of genetic structure and high admixture were detected across the populations located in the central region of Toyama Prefecture. Both contemporary and historical migration analyses showed that dispersal into Toyama Prefecture from neighboring prefectures was high, especially migration from the prefecture on the east into Toyama Prefecture, and bidirectional dispersion between Toyama Prefecture and the prefecture to the south. Knowledge of such genetic structures and population dynamics is required for appropriate management and conservation of sika deer populations in the Japanese Archipelago.

Development of a Novel Mule Deer Genomic Assembly and Species-Diagnostic SNP Panel for Assessing Introgression in Mule Deer, White-Tailed Deer, and Their Interspecific Hybrids

Mule deer (Odocoileus hemionus) are endemic to a wide variety of habitats in western North America, many of which are shared in sympatry with their closely related sister-species white-tailed deer (Odocoileus virginianus), whom they hybridize with in wild populations. Although mule deer meet many ideal conditions for a molecular ecological research species, such as high abundance, ecological importance, and broad dispersal and gene flow, conservation genetic studies have been limited by a relative lack of existing genomic resources and inherent difficulties caused by introgression with white-tailed deer. Many molecular tools currently available for the study of cervids were designed using reference assemblies of divergent model species, specifically cattle (Bos taurus). Bovidae and Cervidae diverged approximately 28 million years ago, therefore, we sought to ameliorate the available resources by contributing the first mule deer whole genome sequence draft assembly with an average genome-wide read depth of 25X, using the white-tailed genome assembly (Ovir.te_1.0) as a reference. Comparing the two assemblies, we identified ∼33 million single nucleotide polymorphisms (SNPs) and insertion/deletion variants. We then verified fixed SNP differences between the two species and developed a 40-loci SNP assay capable of identifying pure mule deer, white-tailed deer, and interspecific hybrids. Assignment capacity of the panel, which was tested on simulated datasets, is reliable up to and including the third backcross hybrid generation. Identification of post-F1 hybrids will be necessary for hybrid zone population studies going forward, and the new mule deer assembly will be a valuable resource for genetic and comparative genomics studies.

Partitioning drivers of spatial genetic variation for a continuously distributed population of boreal caribou: Implications for management unit delineation

Isolation by distance (IBD) is a natural pattern not readily incorporated into theoretical models nor traditional metrics for differentiating populations, although clinal genetic differentiation can be characteristic of many wildlife species. Landscape features can also drive population structure additive to baseline IBD resulting in differentiation through isolation-by-resistance (IBR). We assessed the population genetic structure of boreal caribou across western Canada using nonspatial (STRUCTURE) and spatial (MEMGENE) clustering methods and investigated the relative contribution of IBD and IBR on genetic variation of 1,221 boreal caribou multilocus genotypes across western Canada. We further introduced a novel approach to compare the partitioning of individuals into management units (MU) and assessed levels of genetic connectivity under different MU scenarios. STRUCTURE delineated five genetic clusters while MEMGENE identified finer-scale differentiation across the study area. IBD was significant and did not differ for males and females both across and among detected genetic clusters. MEMGENE landscape analysis further quantified the proportion of genetic variation contributed by IBD and IBR patterns, allowing for the relative importance of spatial drivers, including roads, water bodies, and wildfires, to be assessed and incorporated into the characterization of population structure for the delineation of MUs. Local population units, as currently delineated in the boreal caribou recovery strategy, do not capture the genetic variation and connectivity of the ecotype across the study area. Here, we provide the tools to assess fine-scale spatial patterns of genetic variation, partition drivers of genetic variation, and evaluate the best management options for maintaining genetic connectivity. Our approach is highly relevant to vagile wildlife species that are of management and conservation concern and demonstrate varying degrees of IBD and IBR with clinal spatial genetic structure that challenges the delineation of discrete population boundaries.