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

Exploring interspecific hybridization dynamics in artificial forests of Pinus brutia and P. halepensis: Implications for sustainable afforestation

Interspecific hybridization increases genetic diversity, which is essential for coping with changing environments. Hybrid zones, occurring naturally in overlapping habitats of closely related species, can be artificially established during afforestation. The resulting interspecific hybridization may promote sustainability in artificial forests, particularly in regions facing degradation due to climate change. Currently, there is limited evidence of hybridization during regeneration of artificial forests. Here, we studied the frequency of Pinus brutia Ten. × P. halepensis Mill. hybridization in five planted forests in Israel in three stages of forest regeneration: seeds before dispersal, emerged seedlings and recruited seedlings at the end of the dry season. We found hybrids on P. brutia, but not on P. halepensis trees due to asynchronous cone production phenology. Using 94 single-nucleotide polymorphism (SNP) markers, we found hybrids at all stages, most of which were hybrids of advanced generations. The hybrid proportions increased from 4.7 ± 2.1 to 8.2 ± 1.4 and 21.6 ± 6.4 per cent, from seeds to emerged seedlings and to recruited seedlings stages, respectively. The increased hybrid ratio implies an advantage of hybrids over P. brutia during forest regeneration. To test this hypothesis, we measured seedling growth rate and morphological traits under controlled conditions and found that the hybrid seedlings exhibited selected traits of the two parental species, which likely contributed to the fitness and survival of the hybrids during the dry season. This study highlights the potential contribution of hybrids to sustainable-planted forests and contributes to the understanding of genetic changes that occur during the regeneration of artificial forests.

Development of a tetra-primer ARMS-PCR for identification of sika and red deer and their hybrids

Accurate identification of deer-derived components is significant in food and drug authenticity. Over the years, several methods have been developed to authenticate these products; however, identifying whether female deer products are hybrids is challenging. In this study, the zinc finger protein X-linked (ZFX) gene sequences of sika deer (Cervus nippon), red deer (Cervus elaphus) and their hybrid offspring were amplified and sequenced, the X221 and X428 species-specific single nucleotide polymorphisms (SNP) loci were verified, and a tetra-primer amplification refractory mutation system (T-ARMS-PCR) assay was developed to identify the parent-of-origin of female sika deer, red deer, and their hybrid deer. The T-ARMS-PCR developed based on the X221 locus could identify sika deer, red deer, and their hybrid offspring according to the presence or absence of PCR product sizes of 486 bp, 352 bp, and 179 bp, respectively, just as X428 locus could identify sika deer, red deer, and their hybrid offspring according to the presence or absence of PCR product sizes of 549 bp, 213 bp, and 383 bp, respectively. Forty products labeled deer-derived ingredients randomly purchased were tested using this assay, and the results showed that the identification results based on the two SNP loci were utterly consistent with the actual sources. In addition, this method was found to be accurate, simple, convenient, and with high specificity, thus providing an essential technical reference for deer product species identification. It is also an important supplement to the identification methods of the original ingredients of existing deer products.

Speciation without gene-flow in hybridizing deer

Under the ecological speciation model, divergent selection acts on ecological differences between populations, gradually creating barriers to gene flow and ultimately leading to reproductive isolation. Hybridisation is part of this continuum and can both promote and inhibit the speciation process. Here, we used white-tailed (Odocoileus virginianus) and mule deer (O. hemionus) to investigate patterns of speciation in hybridizing sister species. We quantified genome-wide historical introgression and performed genome scans to look for signatures of four different selection scenarios. Despite ample modern evidence of hybridisation, we found negligible patterns of ancestral introgression and no signatures of divergence with gene flow, rather localized patterns of allopatric and balancing selection were detected across the genome. Genes under balancing selection were related to immunity, MHC and sensory perception of smell, the latter of which is consistent with deer biology. The deficiency of historical gene-flow suggests that white-tailed and mule deer were spatially separated during the glaciation cycles of the Pleistocene and genome wide differentiation accrued via genetic drift. Dobzhansky-Muller incompatibilities and selection against hybrids are hypothesised to be acting, and diversity correlations to recombination rates suggests these sister species are far along the speciation continuum.

Genomic correlates for migratory direction in a free-ranging cervid

Animal migrations are some of the most ubiquitous and one of the most threatened ecological processes globally. A wide range of migratory behaviours occur in nature, and this behaviour is not uniform among and within species, where even individuals in the same population can exhibit differences. While the environment largely drives migratory behaviour, it is necessary to understand the genetic mechanisms influencing migration to elucidate the potential of migratory species to cope with novel conditions and adapt to environmental change. In this study, we identified genes associated with a migratory trait by undertaking pooled genome-wide scans on a natural population of migrating mule deer. We identified genomic regions associated with variation in migratory direction, including FITM1, a gene linked to the formation of lipids, and DPPA3, a gene linked to epigenetic modifications of the maternal line. Such a genetic basis for a migratory trait contributes to the adaptive potential of the species and might affect the flexibility of individuals to change their behaviour in the face of changes in their environment.

Genotyping-in-Thousands by sequencing panel development and application for high-resolution monitoring of introgressive hybridization within sockeye salmon

Stocking programs have been widely implemented to re-establish extirpated fish species to their historical ranges; when employed in species with complex life histories, such management activities should include careful consideration of resulting hybridization dynamics with resident stocks and corresponding outcomes on recovery initiatives. Genetic monitoring can be instrumental for quantifying the extent of introgression over time, however conventional markers typically have limited power for the identification of advanced hybrid classes, especially at the intra-specific level. Here, we demonstrate a workflow for developing, evaluating and deploying a Genotyping-in-Thousands by Sequencing (GT-seq) SNP panel with the power to detect advanced hybrid classes to assess the extent and trajectory of intra-specific hybridization, using the sockeye salmon (Oncorhynchus nerka) stocking program in Skaha Lake, British Columbia as a case study. Previous analyses detected significant levels of hybridization between the anadromous (sockeye) and freshwater resident (kokanee) forms of O. nerka, but were restricted to assigning individuals to pure-stock or “hybrid”. Simulation analyses indicated our GT-seq panel had high accuracy, efficiency and power (> 94.5%) of assignment to pure-stock sockeye salmon/kokanee, F₁, F₂, and B₂ backcross-sockeye/kokanee. Re-analysis of 2016/2017 spawners previously analyzed using TaqMan® assays and otolith microchemistry revealed shifts in assignment of some hybrids to adjacent pure-stock or B₂ backcross classes, while new assignment of 2019 spawners revealed hybrids comprised 31% of the population, ~ 74% of which were B₂ backcross or F₂. Overall, the GT-seq panel development workflow presented here could be applied to virtually any system where genetic stock identification and intra-specific hybridization are important management parameters.

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.

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.

Chromosome-level genome assembly of Tarim red deer, Cervus elaphus yarkandensis

Tarim red deer (Cervus elaphus yarkandensis) is the only subspecies of red deer (of 22 subspecies) from Central Asia. This species is a desert dweller of the Tarim Basin of southern Xinjiang, China, and exhibits some unique adaptations to the dry and extreme hot climate. We report here the assembly of a Tarim red deer genome employing a 10X Genomics library, termed CEY_v1. Our genome consisted of 2.6 Gb with contig N50 and scaffold N50 of 275.5 Kb and 31.7 Mb, respectively. Around 96% of the assembled sequences were anchored onto 34 chromosomes based on the published high-quality red deer genetic linkage map. More than 94% BUSCOs complete genes (including 90.5% single and 3.6% duplicated ones) were detected in the CEY_v1 and 20,653 genes were annotated. The CEY_v1 is expected to contribute to comparative analysis of genome biology, to evolutionary studies within Cervidae, and to facilitating investigation of mechanisms underlying adaptation of this species to the extreme dry and hot climate.

Identification of Ancestry Informative Marker (AIM) Panels to Assess Hybridisation between Feral and Domestic Sheep

Simple Summary

Once present in the entirety of Europe, mouflon (wild sheep) became extinct due to intense hunting, but remnant populations survived and became feral on the Mediterranean islands of Corsica and Sardinia. Although now protected by regional laws, Sardinian mouflon is threatened by crossbreeding with domestic sheep causing genetic hybridisation. The spread of domestic genes can be detrimental for wild populations as it dilutes the genetic features that characterise them. This work aimed to identify diagnostic tools that could be applied to monitor the level of hybridisation between mouflon and domestic sheep. Tens of thousands of genetic markers known as single nucleotide polymorphisms (SNPs) were screened and we identified the smallest number of SNPs necessary to discriminate between pure mouflon and sheep. We produced four SNP panels of different sizes which were able to assess the hybridisation level of a mouflon and we verified that the SNP panels efficacy is independent of the domestic sheep breed involved in the hybrid. The implementation of these results into actual diagnostic tools will help the conservation of this unique and irreplaceable mouflon population, and the methodology applied can easily be transferred to other case studies of interest.

Abstract

Hybridisation of wild populations with their domestic counterparts can lead to the loss of wildtype genetic integrity, outbreeding depression, and loss of adaptive features. The Mediterranean island of Sardinia hosts one of the last extant autochthonous European mouflon (Ovis aries musimon) populations. Although conservation policies, including reintroduction plans, have been enforced to preserve Sardinian mouflon, crossbreeding with domestic sheep has been documented. We identified panels of single nucleotide polymorphisms (SNPs) that could act as ancestry informative markers able to assess admixture in feral x domestic sheep hybrids. The medium-density SNP array genotyping data of Sardinian mouflon and domestic sheep (O. aries aries) showing pure ancestry were used as references. We applied a two-step selection algorithm to this data consisting of preselection via Principal Component Analysis followed by a supervised machine learning classification method based on random forest to develop SNP panels of various sizes. We generated ancestry informative marker (AIM) panels and tested their ability to assess admixture in mouflon x domestic sheep hybrids both in simulated and real populations of known ancestry proportions. All the AIM panels recorded high correlations with the ancestry proportion computed using the full medium-density SNP array. The AIM panels proposed here may be used by conservation practitioners as diagnostic tools to exclude hybrids from reintroduction plans and improve conservation strategies for mouflon populations.

Genomic insights into ruminant evolution: from past to future prospects

Ruminants (Ruminantia) are among the most successful herbivorous mammals, exhibiting wide-ranging morphological and ecological characteristics (such as headgear and multichambered stomach) and including various key livestock species (e.g., cattle, buffalo, yak, sheep, and goat). Understanding their evolution is of great significance not only in scientific research but also in applications potential for human society. The rapid growth of genomic resources provides unprecedented opportunities to dissect the evolutionary histories and molecular mechanisms underlying the distinct characteristics of ruminants. Here we summarize our current understanding of the genetic, morphological, and ecological diversity of ruminants and provide prospects for future studies.

The Value of Reference Genomes in the Conservation of Threatened Species

Conservation initiatives are now more crucial than ever-over a million plant and animal species are at risk of extinction over the coming decades. The genetic management of threatened species held in insurance programs is recommended; however, few are taking advantage of the full range of genomic technologies available today. Less than 1% of the 13505 species currently listed as threated by the International Union for Conservation of Nature (IUCN) have a published genome. While there has been much discussion in the literature about the importance of genomics for conservation, there are limited examples of how having a reference genome has changed conservation management practice. The Tasmanian devil (Sarcophilus harrisii), is an endangered Australian marsupial, threatened by an infectious clonal cancer devil facial tumor disease (DFTD). Populations have declined by 80% since the disease was first recorded in 1996. A reference genome for this species was published in 2012 and has been crucial for understanding DFTD and the management of the species in the wild. Here we use the Tasmanian devil as an example of how a reference genome has influenced management actions in the conservation of a species.