Variable extent of hybridization between invasive sika (Cervus nippon) and native red deer (C. elaphus) in a small geographical area

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.

Molecular Evidence Reveals the Sympatric Distribution of Cervus nippon yakushimae and Cervus nippon taiouanus on Jeju Island, South Korea

Non-native species threaten native ecosystems and species, particularly on islands where rates of endemism and vulnerability to threats are high. Understanding species invasion will aid in providing insights into ecological and evolutionary processes. To identify the non-native sika deer (Cervus nippon) population in Jeju, South Korea, and their phylogenetic affinities, we collected tissue samples from roadkill and the World Natural Heritage Headquarters in Jeju. Mitochondrial DNA cytochrome B (CytB) gene sequences were analyzed to determine two distinct CytB haplotypes. Phylogenetic analysis using maximum likelihood tree revealed two haplotypes of CytB clustered into two different groups representing two subspecies: C. n. yakushimae, native to Japan, and C. n. taiouanus, native to Taiwan. The tentative divergence time between the two subspecies was estimated at 1.81 million years. Our study confirmed that the two subspecies of sika deer are sympatric in the natural ecosystem of Jeju Island. This study provides valuable information to help government and conservation agencies understand alien species and determine control policies for conserving native biodiversity in South Korea.

Development of 16 novel EST-SSR markers for species identification and cross-genus amplification in sambar, sika, and red deer

Deer genera around the globe are threatened by anthropogenic interference. The translocation of alien species and their subsequent genetic introgression into indigenous deer populations is particularly harmful to the species of greatest conservation concern. Products derived from deer, including venison and antler velvet, are also at risk of fraudulent labeling. The current molecular markers used to genetically identify deer species were developed from genome sequences and have limited applicability for cross-species amplification. The absence of efficacious diagnostic techniques for identifying deer species has hampered conservation and wildlife crime investigation efforts. Expressed sequence tag-simple sequence repeat (EST-SSR) markers are reliable tools for individual and species identification, especially in terms of cross-species genotyping. We conducted transcriptome sequencing of sambar (Rusa unicolor) antler velvet and acquired 11,190 EST-SSRs from 65,074 newly assembled unigenes. We identified a total of 55 unambiguous amplicons in sambar (n = 45), which were selected as markers to evaluate cross-species genotyping in sika deer (Cervus nippon, n = 30) and red deer (Cervus elaphus, n = 46), resulting in cross-species amplification rates of 94.5% and 89.1%, respectively. Based on polymorphic information content (>0.25) and genotyping fidelity, we selected 16 of these EST-SSRs for species identification. This marker set revealed significant genetic differentiation based on the fixation index and genetic distance values. Principal coordinate analysis and STRUCTURE analysis revealed distinct clusters of species and clearly identified red-sika hybrids. These markers showed applicability across different genera and proved suitable for identification and phylogenetic analyses across deer species.

Weaving place‐based knowledge for culturally significant species in the age of genomics: Looking to the past to navigate the future

Relationships with place provide critical context for characterizing biocultural diversity. Yet, genetic and genomic studies are rarely informed by Indigenous or local knowledge, processes, and practices, including the movement of culturally significant species. Here, we show how place‐based knowledge can better reveal the biocultural complexities of genetic or genomic data derived from culturally significant species. As a case study, we focus on culturally significant southern freshwater kōura (crayfish) in Aotearoa me Te Waipounamu (New Zealand, herein Aotearoa NZ). Our results, based on genotyping‐by‐sequencing markers, reveal strong population genetic structure along with signatures of population admixture in 19 genetically depauperate populations across the east coast of Te Waipounamu. Environment association and differentiation analyses for local adaptation also indicate a role for hydroclimatic variables—including temperature, precipitation, and water flow regimes—in shaping local adaptation in kōura. Through trusted partnerships between community and researchers, weaving genomic markers with place‐based knowledge has both provided invaluable context for the interpretation of data and created opportunities to reconnect people and place. We envisage such trusted partnerships guiding future genomic research for culturally significant species in Aotearoa NZ and beyond.

Population genomics of free-ranging Great Plains white-tailed and mule deer reflects a long history of interspecific hybridization

Hybridization is a natural process at species-range boundaries that may variably promote the speciation process or break down species barriers but minimally will influence management outcomes of distinct populations. White-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) have broad and overlapping distributions in North America and a recognized capacity for interspecific hybridization. In response to contemporary environmental change to any of one or multiple still-unknown factors, mule deer range is contracting westward accompanied by a westward expansion of white-tailed deer, leading to increasing interactions, opportunities for gene flow, and associated conservation implications. To quantify genetic diversity, phylogenomic structure, and dynamics of hybridization in sympatric populations of white-tailed and mule deer, we used mitochondrial cytochrome b data coupled with SNP loci discovered with double-digest restriction site-associated DNA sequencing. We recovered 25,018 SNPs across 92 deer samples from both species, collected from two regions of western Kansas. Eight individuals with unambiguous external morphology representing both species were of hybrid origin (8.7%), and represented the product of multi-generational backcrossing. Mitochondrial data showed both ancient and recent directional discordance with morphological species assignments, reflecting a legacy of mule deer males mating with white-tailed deer females. Mule deer had lower genetic diversity than white-tailed deer, and both mitochondrial and nuclear data suggest contemporary mule deer effective population decline. Landscape genetic analyses show relative isolation between the two study regions for white-tailed deer, but greater connectivity among mule deer, with predominant movement from north to south. Collectively, our results suggest a long history of gene flow between these species in the Great Plains and hint at evolutionary processes that purge incompatible functional genomic elements as a result of hybridization. Surviving hybrids evidently may be reproductive, but with unknown consequences for the future integrity of these species, population trajectories, or relative susceptibility to emerging pathogens.

Admixture mapping reveals loci for carcass mass in red deer x sika hybrids in Kintyre, Scotland

We deployed admixture mapping on a sample of 386 deer from a hybrid swarm between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon) sampled in Kintyre, Scotland to search for quantitative trait loci (QTLs) underpinning phenotypic differences between the species. These two species are highly diverged genetically [Fst between pure species, based on 50K single nucleotide polymorphism (SNPs) = 0.532] and phenotypically: pure red have on average twice the carcass mass of pure sika in our sample (38.7 kg vs 19.1 kg). After controlling for sex, age, and population genetic structure, we found 10 autosomal genomic locations with QTL for carcass mass. Effect sizes ranged from 0.191 to 1.839 kg and as expected, in all cases the allele derived from sika conferred lower carcass mass. The sika population was fixed for all small carcass mass alleles, whereas the red deer population was typically polymorphic. GO term analysis of genes lying in the QTL regions are associated with oxygen transport. Although body mass is a likely target of selection, none of the SNPs marking QTL are introgressing faster or slower than expected in either direction.

Extent and direction of introgressive hybridization of mule and white-tailed deer in western Canada

Hybridization of mule deer (Odocoileus hemionus) and white-tailed deer (O. virginianus) appears to be a semi-regular occurrence in western North America. Previous studies confirmed the presence of hybrids in a variety of sympatric habitats, but their developing molecular resources limited identification to the earliest, most admixed generations. For this reason, estimates of hybrid production in wild populations often rely on anecdotal reports. As well, white-tailed deer populations' continued encroachment into historically mule deer-occupied habitats due to changes in land use, habitat homogenization, and a warming climate may increase opportunities for interspecific encounters. We sought to quantify the prevalence and extent of hybrid deer in the prairies of western Canada using a SNP assay with enhanced discriminating power. By updating the available molecular resources, we sought to identify and characterize previously cryptic introgression. We also investigated the influence of various parameters on hybridity by way of logistic regression. We observed overall hybridization rates of ~1.0%, slightly lower than that reported by previous studies, and found white-tailed-like hybrids to be more common than their mule deer-like counterparts. Here, we build upon past studies of hybridization in North American deer by increasing hybrid detection power, expanding sample sizes, demonstrating a new molecular resource applicable to future research and observing asymmetrical directionality of introgression.

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.

Locus-specific introgression in young hybrid swarms: Drift may dominate selection

Closely related species that have previously inhabited geographically separated ranges are hybridizing at an increasing rate due to human disruptions. These human-mediated hybrid zones can be used to study reproductive isolation between species at secondary contact, including examining locus-specific rates of introgression. Introgression is expected to be heterogenous across the genome, reflecting variation in selection. Those loci that introgress especially slowly are good candidates for being involved in reproductive isolation, while those loci that introgress quickly may be involved in adaptive introgression. In the context of conservation, policy makers are especially concerned about introduced alleles moving quickly into the background of a native or endemic species, as these alleles could replace the native alleles in the population, leading to extinction via hybridization. We applied genomic cline analyses to 44,997 SNPs to identify loci introgressing more or less when compared to the genome wide expectation in a human-mediated hybridizing population of red deer and sika in Kintyre Scotland. We found 11.4% of SNPs had cline centres that were significantly different from the genome wide expectation, and 17.6% of all SNPs had excess rates of introgression. Based on simulations, we believe that many of these markers have diverged from the genome-wide average due to drift, rather than because of selection, and we suggest that these simulations can be useful as a null distribution for future studies of genomic clines. Future work on red deer and sika could determine the policy implications of allelic-replacement due to drift rather than selection, and could use replicate, geographically distinct hybrid zones to narrow down those loci that are responding to selection.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.

Breeding behaviour predicts patterns of natural hybridization in North American minnows (Cyprinidae)

Premating barriers such as variation in reproductive behaviour can evolve quickly, but because gametic and postzygotic incompatibilities often evolve more slowly, circumstances that bring gametes into contact can breach the boundaries of premating isolation. In aquatic environments, the gametes of organisms with external fertilization are released into a constantly moving environment and may come into contact with heterospecific gametes. In fishes, nest association (spawning in another species' nest) is a behaviour that brings gametes from different species into close spatiotemporal proximity. These interactions might increase chances of hybridization, especially when multiple species associate with a single nest builder. This study addresses these interactions in the largest clade of North American freshwater fishes, the minnows (Cyprinidae). We compiled a list of over 17,000 hybrid specimens in conjunction with species distribution data, breeding behaviours, and an inferred phylogeny to test if breeding behaviour, in addition to evolutionary history, is an important predictor of hybridization. We find that breeding behaviour is a significant predictor of hybridization, even when phylogenetic relatedness and divergence time are accounted for. Specifically, nest associates are more likely to hybridize with other nest associates whereas non-nesting species had relatively low rates of hybridization.

Repeated hybridization of two closely related gazelle species (Gazella bennettii and Gazella subgutturosa) in central Iran

Interspecific hybridization increasingly occurs in the course of anthropogenic actions, such as species translocations and introductions, and habitat modifications or occurs in sympatric species due to the shortage of conspecific mates. Compared with anthropogenically caused hybridization, natural hybridization is more difficult to prove, but both play an important role in conservation. In this study, we detected hybridization of two gazelle sister species, Gazella bennettii (adapted to dry areas) and Gazella subgutturosa (adapted to open plains), in five habitat areas, where G. bennettii naturally occur in central Iran. The hybrids have a nuclear genomic identity (based on two introns), habitat preference, and phenotype of G. bennettii, but the mitochondrial identity (based on cyt b) of G. subgutturosa. We suggest that natural hybridization of female G. subgutturosa and male G. bennettii happened twice in central Iran in prehistoric times, based on the haplotype pattern that we found. However, we found indications of recent hybridization between both species under special circumstances, for example, in breeding centers, due to translocations, or in areas of sympatry due to the shortage of conspecific mates. Therefore, these two species must be kept separately in the breeding centers, and introduction of one of them into the habitat of the other must be strictly avoided.

Genetic population structure of invasive raccoons (Procyon lotor) in Hokkaido, Japan: Unique phenomenon caused by pet escape or abandonment

Phylogeographic studies can resolve relationships between genetic population structure of organisms and geographical distributions. Raccoons have become feral in Japan, and in Hokkaido island, they have been rapidly increasing in number and spreading since the 1970s. We analyzed mitochondrial (mtDNA) and microsatellite DNA to understand the current phylogenetic distribution and invasive founder events. Overall, Hokkaido raccoons maintained high genetic diversity (i.e., the level of heterozygosity was comparable to the original habitat, North America). Based on mtDNA distribution and microsatellite diversity, Hokkaido raccoons were divided into six management units. However, mtDNA haplotype distributions and genetic structures based on microsatellites did not always correspond to each other (e.g., two geographically and genetically separated populations showed similar mtDNA distributions). In addition, a high degree of genetic admixture was observed in every unit, and the degree of genetic differentiation was low even between regions separated by long distances. Compared with other countries in Europe where genetic distribution of introduced raccoons is more clearly structured, the current results represent a unique and complex phenomenon of pet escape/abandonment in Hokkaido: i.e., genetically related colonies were introduced into multiple regions as founder events, resulting in the current state in which raccoons are not clearly genetically differentiated even 40 years after introduction.

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

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

Microsatellite based assignment reveals history of extirpated mountain ungulate

During the early 1900s, Northern chamois (Rupicapra rupicapra) populations in the northern Dinaric Mountains were extirpated. During the 1960s and 1970s there were several reintroductions of individuals from two Northern chamois subspecies (Alpine chamois, R. r. rupicapra and Balkan chamois, R. r. balcanica) from neighbouring areas in the attempt to re-establish the population. Accurate taxonomic classification, at subspecies level, of the autochthonous extirpated population was not known. To clarify which subspecies was present before reintroduction, we genotyped four male chamois skulls originating from Velebit Mountain, collected around 25 years before the population local extinction. DNA was successfully extracted from middle layer and outer sheath of horns. Assignment based on microsatellite loci, using both Bayesian clustering in STRUCTURE (with q values between 0.55 and 0.73) and DAPC (with individual membership probabilities of 0.99 and 1.00) indicated higher assessed likelihood for the Alpine subspecies.

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.

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.

Widespread hybridization in the introduced hog deer population of Victoria, Australia, and its implications for conservation

In Australia, many species have been introduced that have since undergone drastic declines in their native range. One species of note is the hog deer (Axis porcinus) which was introduced in the 1860s to Victoria, Australia, and has since become endangered in its native range throughout South-East Asia. There is increased interest in using non-native populations as a source for genetic rescue; however, considerations need to be made of the genetic suitability of the non-native population. Three mitochondrial markers and two nuclear markers were sequenced to assess the genetic variation of the Victorian population of hog deer, which identified that the Victorian population has hybrid origins with the closely related chital (Axis axis), a species that is no longer present in the wild in Victoria. In addition, the mitochondrial D-loop region within the Victorian hog deer is monomorphic, demonstrating that mitochondrial genetic diversity is very low within this population. This study is the first to report of long-term persistence of hog deer and chital hybrids in a wild setting, and the continual survival of this population suggests that hybrids of these two species are fertile. Despite the newly discovered hybrid status in Victorian hog deer, this population may still be beneficial for future translocations within the native range. However, more in-depth analysis of genetic diversity within the Victorian hog deer population and investigation of hybridization rates within the native range are necessary before translocations are attempted.

Variation in the prion protein gene (PRNP) sequence of wild deer in Great Britain and mainland Europe

Susceptibility to prion diseases is largely determined by the sequence of the prion protein gene (PRNP), which encodes the prion protein (PrP). The recent emergence of chronic wasting disease (CWD) in Europe has highlighted the need to investigate PRNP gene diversity in European deer species, to better predict their susceptibility to CWD. Here we report a large genetic survey of six British deer species, including red (Cervus elaphus), sika (Cervus nippon), roe (Capreolus capreolus), fallow (Dama dama), muntjac (Muntiacus reevesii), and Chinese water deer (Hydropotes inermis), which establishes PRNP haplotype and genotype frequencies. Two smaller data sets from red deer in Norway and the Czech Republic are also included for comparison. Overall red deer show the most PRNP variation, with non-synonymous/coding polymorphisms at codons 98, 168, 226 and 247, which vary markedly in frequency between different regions. Polymorphisms P168S and I247L were only found in Scottish and Czech populations, respectively. T98A was found in all populations except Norway and the south of England. Significant regional differences in genotype frequencies were observed within both British and European red deer populations. Other deer species showed less variation, particularly roe and fallow deer, in which identical PRNP gene sequences were found in all individuals analysed. Based on comparison with PRNP sequences of North American cervids affected by CWD and limited experimental challenge data, these results suggest that a high proportion of wild deer in Great Britain may be susceptible to CWD.

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.

Advances in population ecology and species interactions in mammals

The study of mammals has promoted the development and testing of many ideas in contemporary ecology. Here we address recent developments in foraging and habitat selection, source–sink dynamics, competition (both within and between species), population cycles, predation (including apparent competition), mutualism, and biological invasions. Because mammals are appealing to the public, ecological insight gleaned from the study of mammals has disproportionate potential in educating the public about ecological principles and their application to wise management. Mammals have been central to many computational and statistical developments in recent years, including refinements to traditional approaches and metrics (e.g., capture-recapture) as well as advancements of novel and developing fields (e.g., spatial capture-recapture, occupancy modeling, integrated population models). The study of mammals also poses challenges in terms of fully characterizing dynamics in natural conditions. Ongoing climate change threatens to affect global ecosystems, and mammals provide visible and charismatic subjects for research on local and regional effects of such change as well as predictive modeling of the long-term effects on ecosystem function and stability. Although much remains to be done, the population ecology of mammals continues to be a vibrant and rapidly developing field. We anticipate that the next quarter century will prove as exciting and productive for the study of mammals as has the recent one.

Distinguishing the victim from the threat: SNP-based methods reveal the extent of introgressive hybridization between wildcats and domestic cats in Scotland and inform future in situ and ex situ management options for species restoration

The degree of introgressive hybridization between the Scottish wildcat and domestic cat has long been suspected to be advanced. Here, we use a 35-SNP-marker test, designed to assess hybridization between wildcat and domestic cat populations in Scotland, to assess a database of 295 wild-living and captive cat samples, and test the assumptions of the test using 3,097 SNP markers generated independently in a subset of the data using ddRAD. We discovered that despite increased genetic resolution provided by these methods, wild-living cats in Scotland show a complete genetic continuum or hybrid swarm structure when judged against reference data. The historical population of wildcats, although hybridized, clearly groups at one end of this continuum, as does the captive population of wildcats. The interpretation of pelage scores against nuclear genetic data continues to be problematic. This is probably because of a breakdown in linkage equilibrium between wildcat pelage genes as the two populations have become increasingly mixed, meaning that pelage score or SNP score alone is poor diagnostic predictors of hybrid status. Until better tools become available, both should be used jointly, where possible, when making management decisions about individual cats. We recommend that the conservation community in Scotland must now define clearly what measures are to be used to diagnose a wildcat in the wild in Scotland, if future conservation action is to be effective.

Development of Diagnostic SNP Markers To Monitor Hybridization Between Sika Deer (Cervus nippon) and Wapiti (Cervus elaphus)

Sika deer (Cervus Nippon) and wapiti (Cervus elaphus) are closely related species and their hybridization can result in significant allele-shift of their gene pool. Additive genetic effects and putative heterotic effects of their hybridization on growth performance could confer considerable economic advantage in deer farming. Here, we used double-digest restriction site-associated DNA sequencing technology (ddRAD-seq) and detected ∼320,000 genome-wide SNPs from 30 captive individuals: 7 sika deer, 6 wapiti and 17 F1 hybrids (reciprocal cross). By screening observed heterozygosity of each SNP across four taxonomic groups, we report for the first time a resource of 2,015 putative diagnostic SNP markers (species-specific SNPs for sika deer and wapiti), which can be used to design tools for assessing or monitoring the degree of hybridization between sika deer and wapiti. These ddRAD-seq data and SNP datasets are also valuable resources for genome-wide studies, including trait discovery for breeders of domestic deer.

Transmission of Helminths between Species of Ruminants in Austria Appears More Likely to Occur than Generally Assumed

Helminth infections of the gastrointestinal tract and lungs can lead to devastating economical losses to the pastoral based animal production. Farm animals can suffer from malnutrition, tissue damage, and blood loss resulting in impaired production traits and reproduction parameters. In Austria, pastures grazed by sheep, goats, and cattle overlap with the habitats of several species of wild cervids (roe deer, red deer, sika deer, and fallow deer) and bovids (mouflon, chamois, and ibex), and transmission of parasites between different ruminant species seems likely. A complete and updated overview on the occurrence of helminths of domestic and wild ruminants in Austria is presented. Based on these data, intersections of the host spectrum of the determined parasites were depicted. The “liability index” was applied to identify the ruminant species, which most likely transmit parasites between each other. A degree for host specificity was calculated for each parasite species based on the average taxonomic distance of their host species. Of the 73 identified helminth species 42 were identified as generalists, and 14 transmission experiments supported the assumed broad host specificity for 14 generalists and 1 specialist helminth species. Overall, 61 helminths were found to infect more than one host species, and 4 were found in all 10 ruminant species investigated. From these analyses, it can be concluded that a number of helminth parasites of the gastrointestinal tract and the lungs are potentially transmitted between domestic and wild ruminants in Austria. For some parasites and host species, experimental evidence is in support for possible transmission, while for other such studies are lacking. Host preference of different genotypes of the same parasite species may have a confounding effect on the evaluation of cross-transmission, but so far this has not been evaluated systematically in helminths in Austria. Further studies focusing on experimental cross-transmission and genetic characterization are needed to define the potential consequences for the epidemiology of those parasites, animal welfare, and economic impact.

Introgression of exotic Cervus (nippon and canadensis) into red deer (Cervus elaphus) populations in Scotland and the English Lake District

Since the mid-19th century, multiple introductions of Japanese sika deer (Cervus nippon nippon) and North American wapiti (C. canadensis) have taken place in the British Isles. While wapiti have generally been unsuccessful, sika have been very successful, especially in Scotland where they now overlap at least 40% of the range of native red deer (C. elaphus). Hybridization between these two species and red deer has been demonstrated in captivity and in the wild. Using a panel of 22 microsatellite loci that are highly diagnostic between red deer and sika, and moderately diagnostic between red deer and wapiti, we investigated the extent of introgression between these species in 2,943 deer sampled from around Scotland and from the English Lake District using the Bayesian clustering software STRUCTURE. We also used a diagnostic mitochondrial marker for red deer and sika. Our survey extends previous studies indicating little introgression of wapiti nuclear alleles into red deer, in particular in Northern Scotland, Kintyre, and the Lake District. We found a new area of extensive sika introgression in South Kintyre. In the North Highlands, we show for the first time geographically scattered evidence of past hybridization followed by extensive backcrossing, including one red-like individual with sika introgression, two sika-like individuals with red deer introgression, and six individuals that were apparently pure sika at the nuclear markers assessed but which carried red deer mitochondria. However, there has not been a collapse of assortative mating in this region. Similarly, in the English Lake District red deer, we found only traces of past sika introgression. No sika alleles were detected in the Central Highlands or the Hebridean red deer refugia. We make suggestions for management to prevent further spread of sika alleles into red deer and vice versa.

Human-mediated introduction of introgressed deer across Wallace's line: Historical biogeography of Rusa unicolor and R. timorensis

In this study we compared the phylogeographic patterns of two Rusa species, Rusa unicolor and Rusa timorensis, in order to understand what drove and maintained differentiation between these two geographically and genetically close species and investigated the route of introduction of individuals to the islands outside of the Sunda Shelf. We analyzed full mitogenomes from 56 archival samples from the distribution areas of the two species and 18 microsatellite loci in a subset of 16 individuals to generate the phylogeographic patterns of both species. Bayesian inference with fossil calibration was used to estimate the age of each species and major divergence events. Our results indicated that the split between the two species took place during the Pleistocene, ~1.8 Mya, possibly driven by adaptations of R. timorensis to the drier climate found on Java compared to the other islands of Sundaland. Although both markers identified two well-differentiated clades, there was a largely discrepant pattern between mitochondrial and nuclear markers. While nDNA separated the individuals into the two species, largely in agreement with their museum label, mtDNA revealed that all R. timorensis sampled to the east of the Sunda shelf carried haplotypes from R. unicolor and one Rusa unicolor from South Sumatra carried a R. timorensis haplotype. Our results show that hybridization occurred between these two sister species in Sundaland during the Late Pleistocene and resulted in human-mediated introduction of hybrid descendants in all islands outside Sundaland.

A High-Density Linkage Map Reveals Sexual Dimorphism in Recombination Landscapes in Red Deer (Cervus elaphus)

High-density linkage maps are an important tool to gain insight into the genetic architecture of traits of evolutionary and economic interest, and provide a resource to characterize variation in recombination landscapes. Here, we used information from the cattle genome and the 50 K Cervine Illumina BeadChip to inform and refine a high-density linkage map in a wild population of red deer (Cervus elaphus). We constructed a predicted linkage map of 38,038 SNPs and a skeleton map of 10,835 SNPs across 34 linkage groups. We identified several chromosomal rearrangements in the deer lineage relative to sheep and cattle, including six chromosome fissions, one fusion, and two large inversions. Otherwise, our findings showed strong concordance with map orders in the cattle genome. The sex-averaged linkage map length was 2739.7 cM and the genome-wide autosomal recombination rate was 1.04 cM/Mb. The female autosomal map length was 1.21 longer than that of males (2767.4 cM vs. 2280.8 cM, respectively). Sex differences in map length were driven by high female recombination rates in peri-centromeric regions, a pattern that is unusual relative to other mammal species. This effect was more pronounced in fission chromosomes that would have had to produce new centromeres. We propose two hypotheses to explain this effect: (1) that this mechanism may have evolved to counteract centromeric drive associated with meiotic asymmetry in oocyte production; and/or (2) that sequence and structural characteristics suppressing recombination in close proximity to the centromere may not have evolved at neo-centromeres. Our study provides insight into how recombination landscapes vary and evolve in mammals, and will provide a valuable resource for studies of evolution, genetic improvement, and population management in red deer and related species.

Male behaviour drives assortative reproduction during the initial stage of secondary contact

Phenotypic divergence in allopatry can facilitate speciation by reducing the likelihood that individuals of different lineages hybridize during secondary contact. However, few studies have established the causes of reproductive isolation in the crucial early stages of secondary contact. Here, we establish behavioural causes of assortative reproduction between two phenotypically divergent lineages of the European wall lizard (Podarcis muralis), which have recently come into secondary contact. Parentage was highly assortative in experimental contact zones. However, despite pronounced divergence in male phenotypes, including chemical and visual sexual signals, there was no evidence that females discriminated between males of the two lineages in staged interactions or under naturalistic free-ranging conditions. Instead, assortative reproduction was driven by male mate preferences and, to a lesser extent, male-male competition. The effects were more pronounced when the habitat structure promoted high lizard densities. These results emphasize that assortative reproduction can occur in the absence of female choice and that male behaviour may play an important role in limiting hybridization during the initial stages of secondary contact.

Female Sexual Preferences Toward Conspecific and Hybrid Male Mating Calls in Two Species of Polygynous Deer, Cervus elaphus and C. nippon

The behavioral processes at the basis of hybridization and introgression are understudied in terrestrial mammals. We use a unique model to test the role of sexual signals as a reproductive barrier to introgression by investigating behavioral responses to male sexual calls in estrous females of two naturally allopatric but reproductively compatible deer species, red deer and sika deer. Previous studies demonstrated asymmetries in acoustic species discrimination between these species: most but not all female red deer prefer conspecific over sika deer male calls while female sika deer exhibit no preference differences. Here, we extend this examination of acoustic species discrimination to the role of male sexual calls in introgression between parent species and hybrids. Using two-speaker playback experiments, we compared the preference responses of estrous female red and sika deer to male sexual calls from conspecifics versus red × sika hybrids. These playbacks simulate early secondary contact between previously allopatric species after hybridization has occurred. Based on previous conspecific versus heterospecific playbacks, we predicted that most female red deer would prefer conspecific calls while female sika deer would show no difference in their preference behaviors toward conspecific and hybrid calls. However, results show that previous asymmetries did not persist as neither species exhibited more preferences for conspecific over hybrid calls. Thus, vocal behavior is not likely to deter introgression between these species during the early stages of sympatry. On a wider scale, weak discrimination against hybrid sexual signals could substantially contribute to this important evolutionary process in mammals and other taxa.

Identification of deer species (Cervidae, Cetartiodactyla) in China using mitochondrial cytochrome c oxidase subunit I (mtDNA COI)

In China, many deer species are threatened and their wild populations are decreasing. In this study, a segment of mitochondrial cytochrome oxidase subunit I (COI) was used as a DNA barcode to identify Cervidae species. COI sequences of 30 individuals from nine species were determined. Together with 148 sequences from BOLD and Genbank, a total of 178 sequences from 21 species of the family Cervidae were analyzed. The results showed that all species had unique COI sequences, and there was no barcode sharing among them. The mean K2P distances within species, genus, and family were 1.3%, 3.4%, and 8.9%, respectively. The neighbor-joining (NJ) tree was in most cases concordant with modern deer classification. Three species showed maximum intraspecific divergences higher than their minimum interspecific divergences. However, all species could be discriminated by their diagnostic characters in BLOG analysis. The present study confirmed that COI barcodes can effectively distinguish Cervidae species in China.

Genetic diversity and population structure of a Sichuan sika deer (Cervus sichuanicus) population in Tiebu Nature Reserve based on microsatellite variation

Cervus sichuanicus is a species of sika deer (Cervus nippon Group). To date, research has mainly focused on quantity surveying and behavior studies, with genetic information on this species currently deficient. To provide scientific evidence to assist in the protection of this species, we collected Sichuan sika deer fecal samples from the Sichuan Tiebu Nature Reserve (TNR) and extracted DNA from those samples. Microsatellite loci of bovine were used for PCR amplification. After GeneScan, the genotype data were used to analyze the genetic diversity and population structure of the Sichuan sika deer in TNR. Results showed that the average expected heterozygosity of the Sichuan sika deer population in TNR was 0.562, equivalent to the average expected heterozygosity of endangered animals, such as Procapra przewalskii. Furthermore, 8 of 9 microsatellite loci significantly deviated from the Hardy-Weinberg equilibrium and two groups existed within the Sichuan sika deer TNR population. This genetic structure may be caused by a group of Manchurian sika deer (Cervus hortulorum) released in TNR.

Hybridization of two major termite invaders as a consequence of human activity

While hybridization of an invasive species with a native species is a common occurrence, hybridization between two invasive species is rare. Formosan subterranean termites (Coptotermes formosanus) and Asian subterranean termites (C. gestroi) are both ecologically successful and are the two most economically important termite pests in the world. Both species have spread throughout many areas of the world due to human activity; however, their distributions overlap in only three narrow areas because of distinct ecological requirements. In south Florida, where C. formosanus and C. gestroi are both invasive, the dispersal flight seasons of both species overlapped for the first time on record in 2013 and 2014. Pairings of heterospecific individuals were readily observed in the field and C. gestroi males preferentially engaged in mating behavior with C. formosanus females rather than females from their own species. In the laboratory, heterospecific and conspecific pairings had an equal colony establishment rate, but heterospecific incipient colonies had twice the growth rate of conspecific incipient colonies, suggesting a potential case of hybrid vigor. As all pre-zygotic barriers were lifted between the two species in the field, the apparent absence of post-zygotic barriers in the laboratory raises the possibility for introgressive hybridization in south Florida. While laboratory observations remain to be confirmed in the field, and the alate hybrid fertility is currently unknown, our results raise a tangible concern about the hybridization of two major destructive pest species. Such hybridization would likely be associated with a new economic impact.