A comparison of microsatellites and genome-wide SNPs for the detection of admixture brings the first molecular evidence for hybridization between Mustela eversmanii and M. putorius (Mustelidae, Carnivora)

Red, gold, and green: comparative genomics of polymorphic leopards from South Africa

An important goal of comparative and functional genomics is to connect genetic polymorphisms to phenotypic variation. Leopards (Panthera pardus) from northern South Africa are particularly diverse, as here a unique color morph occurs, as well as two deeply diverged southern (SA) and central African (CA) mitochondrial clades, stemming from Pleistocene refugia. Here, we present the first whole genomes of a red leopard and a black (captive) leopard, and wildtypes belonging to the CA and SA mitochondrial clades, to evaluate genome-wide diversity, divergence, and high-impact mutations that may relate to their phenotype. In the black leopard, we found long runs of homozygosity (ROHs), low nucleotide diversity across the genome, and a large number of homozygous structural variants, likely resulting from inbreeding to maintain this color morph in captivity. In red leopards, runs of homozygosity were slightly longer compared to wildtype leopards, with potential deleterious mutations relating to its phenotype, including impaired vision. When assessing population structure, we found no divergence between CA and SA leopards and the rest of Africa, whether comparing single nucleotide or structural variants. This illustrates the homogenizing effect of introgression, and highlights that although leopards in northern South Africa may be phenotypically unique, they are not genetically different.

Hybridization and Introgression in Black Hakes (Merluccius polli and M. senegalensis): Evolutionary Dynamics and Conservation Implications in the Contact Zone Exploited by Multi-Species Fisheries

Hybridization is relatively common between closely related species that share part of their distribution. Understanding its dynamics is important both for conservation purposes and to determine its role as an evolutionary mechanism. Here we have studied the case of black hakes (Merluccius polli and Merluccius senegalensis) in its contact zone. The area of study is located in the FAO fishing area 34, in Mauritania and Senegal waters, where both species are exploited jointly in multi-species fisheries involving national and foreign fleets. Using a ddRADSeq approach and based on a set of 5820 SNPs and a total of 240 individuals, we identified one F1 hybrid and several backcrossed individuals among 90 M. polli samples and none in 90 M. senegalensis samples obtained in 2020, suggesting unidirectional introgression towards M. polli. Hybridization signals were not found in any of the 60 historical samples from 2000. Excluding the hybrids and developing two separate sets of SNPs (5093 SNPs for M. polli and 2794 SNPs for M. senegalensis), our results detected two distinct genetic clusters within M. polli that show different genetic diversity estimates, with one of the clusters showing a higher potential vulnerability to exploitation. This pattern was observed in both contemporary and historical samples, and both groups presented subtle depth segregation. Moreover, 109 outlier loci were identified between the two groups, that could be developed into molecular markers to further study differentiation between both clusters and contribute to improved stock assessment and management of these important demersal resources.

Genomic Introgression Between Critically Endangered and Stable Species of Darwin's Tree Finches on the Galapagos Islands

Natural hybridisation among rare or endangered species and stable congenerics is increasingly topical for the conservation of species-level diversity under anthropogenic impacts. Evidence for beneficial genes being introgressed into or selected for in hybrids raises concurrent questions about its evolutionary significance. In Darwin's tree finches on the island of Floreana (Galapagos Islands, Ecuador), the Critically Endangered medium tree finch (Camarhynchus pauper) undergoes introgression with the stable small tree finch (Camarhynchus parvulus), and hybrids regularly backcross with C. parvulus. Earlier studies in 2005-2013 documented an increase in the frequency of Camarhynchus hybridisation on Floreana using field-based and microsatellite data. With single nucleotide polymorphism (SNP) data from the same Floreana tree finches sampled in 2005 and 2013 (n = 95), we examine genome-wide divergence across parental and hybrid birds and evidence for selection in hybrids. We found that just 18% of previously assigned hybrid birds based on microsatellites could be assigned to hybrids using SNPs. Over half of the previously assigned hybrids (63%) were reassigned to C. parvulus, though parental species showed concordance with prior assignments. Of 4869 private alleles found in hybrid birds, 348 were at a high frequency (≥ 0.30) that exceeded their parental species of origin 89%-96% of the time. For private alleles detected in both years (N = 536) between 11%-76% of alleles underwent a frequency increase and 13%-61% a frequency decrease between 2005 and 2013, which was sensitive to sampling effort. We identified 28 private alleles that were candidates under selection via local PCA and outlier tests. Alleles were annotated to genes associated with inflammation, immunity, brain function and development. We provide evidence that introgression among a critically endangered and stable species of Darwin's tree finch across years may aid in the retention of adaptive alleles and genetic diversity in birds threatened with extinction.

Exposure to anticoagulant rodenticides in steppe polecat (Mustela eversmanii) and European polecat (Mustela putorius) in central Europe

Poisoning caused by coumarin-type anticoagulant rodenticides (ARs) stands as the predominant method for controlling rodents globally. ARs, through secondary poisoning, pose a significant threat to predators due to their lethal and sublethal effects. We examined the concentration of accumulated ARs in liver samples of mostly road-killed steppe polecats (Mustela eversmanii) and European polecats (M. putorius) collected throughout Hungary between 2005 and 2021. The steppe polecat samples were found mainly from Eastern Hungary, while European polecats from Western Hungary. We measured the concentration of six residues by HPLC-FLD. Our analysis revealed the presence of one first-generation and four second-generation ARs in 53% of the steppe polecat (36) and 39% of the European polecat (26) samples. In 17 samples we detected the presence of at least two AR compounds. Although we did not find significant variance in AR accumulation between the two species, steppe polecats displayed greater prevalence and maximum concentration of ARs, whereas European polecat samples exhibited a more diverse accumulation of these compounds. Brodifacoum and bromadiolone were the most prevalent ARs; the highest concentrations were 0.57 mg/kg and 0.33 mg/kg, respectively. The accumulation of ARs was positively correlated with human population density and negatively correlated with the extent of the more natural habitats in both species. To the best of our knowledge, this is the first study to demonstrate anticoagulant rodenticide exposure in steppe polecats globally, and for European polecats in Central European region. Although the extent of AR accumulation in European polecat in Hungary appears comparatively lower than in many other European countries, the issue of secondary poisoning remains a serious problem as these ARs intrude into food webs. Reduced and more prudent usage of pesticides would provide several benefits for wildlife, included humans. However, we advocate a prioritization of ecosystem services through the complete prohibition of the toxicants.

Unraveling the genomic diversity and admixture history of captive tigers in the United States

Genomic studies of endangered species have primarily focused on describing diversity patterns and resolving phylogenetic relationships, with the overarching goal of informing conservation efforts. However, few studies have investigated genomic diversity housed in captive populations. For tigers (Panthera tigris), captive individuals vastly outnumber those in the wild, but their diversity remains largely unexplored. Privately owned captive tiger populations have remained an enigma in the conservation community, with some believing that these individuals are severely inbred, while others believe they may be a source of now-extinct diversity. Here, we present a large-scale genetic study of the private (non-zoo) captive tiger population in the United States, also known as "Generic" tigers. We find that the Generic tiger population has an admixture fingerprint comprising all six extant wild tiger subspecies. Of the 138 Generic individuals sequenced for the purpose of this study, no individual had ancestry from only one subspecies. We show that the Generic tiger population has a comparable amount of genetic diversity relative to most wild subspecies, few private variants, and fewer deleterious mutations. We observe inbreeding coefficients similar to wild populations, although there are some individuals within both the Generic and wild populations that are substantially inbred. Additionally, we develop a reference panel for tigers that can be used with imputation to accurately distinguish individuals and assign ancestry with ultralow coverage (0.25×) data. By providing a cost-effective alternative to whole-genome sequencing (WGS), the reference panel provides a resource to assist in tiger conservation efforts for both ex- and in situ populations.

Contrasting effect of hybridization on genetic differentiation in three rockfish species with similar life history

Hybridization can provide evolutionary benefits (e.g., population resilience to climate change) through the introduction of adaptive alleles and increase of genetic diversity. Nevertheless, management strategies may be designed based only on the parental species within a hybrid zone, without considering the hybrids. This can lead to ineffective spatial management of species, which can directly harm population diversity and negatively impact food webs. Three species of rockfish (Brown Rockfish (Sebastes caurinus), Copper Rockfish (S. auriculatus), and Quillback Rockfish (S. maliger)) are known to hybridize within Puget Sound, Washington, but genetic data from these species are used to infer population structure in the entire genus, including in species that do not hybridize. The goal of this project was to estimate the hybridization rates within the region and determine the effect of hybridization on geographic patterns of genetic structure. We sequenced 290 Brown, Copper, and Quillback rockfish using restriction-site associated DNA sequencing (RADseq) from four regions within and outside Puget Sound, Washington. We show that (i) hybridization within Puget Sound was asymmetrical, not recent, widespread among individuals, and relatively low level within the genome, (ii) hybridization affected population structure in Copper and Brown rockfish, but not in Quillback Rockfish and (iii) after taking hybridization into account we found limited directional dispersal in Brown and Copper rockfish, and evidence for two isolated populations in Quillback Rockfish. Our results suggest that rockfish population structure is species-specific, dependent on the extent of hybridization, and cannot be inferred from one species to another despite similar life history.

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe

BACKGROUND

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry.

METHODS AND RESULTS

MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank.

CONCLUSIONS

Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Influence of molecular marker type on estimating effective population size and other genetic parameters in a critically endangered parrot

Genetics is a fast-moving field, and for conservation practitioners or ecologists, it can be bewildering. The choice of marker used in studies is fundamental; in the literature, preference has recently shifted from microsatellites to single nucleotide polymorphism (SNP) loci. Understanding how marker type affects estimates of population genetic parameters is important in the context of conservation, especially because the accuracy of estimates has a bearing on the actions taken to protect threatened species. We compare parameter estimates between seven microsatellites, 3761 SNP loci, and a random subset of 100 SNPs for the exact same 324 individual swift parrots, Lathamus discolor, and also use 457 additional samples from subsequent years to compare SNP estimates. Both marker types estimated a lower H O than H E. We show that microsatellites and SNPs mainly indicate a lack of spatial genetic structure, except when a priori collection locations were used on the SNP data in a discriminant analysis of principal components (DAPC). The 100-SNP subset gave comparable results to when the full dataset was used. Estimates of effective population size (N e) were comparable between markers when the same individuals were considered, but SNPs had narrower confidence intervals. This is reassuring because conservation assessments that rely on population genetic estimates based on a few microsatellites are unlikely to be nullified by the general shift toward SNPs in the literature. However, estimates between markers and datasets varied considerably when only adult samples were considered; hence, including samples of all age groups is recommended to be used when available. The estimated N e was higher for the full SNP dataset (2010-2019) than the smaller comparison data (2010-2015), which might be a better reflection of the species status. The lower precision of microsatellites may not necessarily be a barrier for most conservation applications; however, SNPs will improve confidence limits, which may be useful for practitioners.

Unmasking microsatellite deceptiveness and debunking hybridization with SNPs in four marine copepod species of Calanus

Interspecific hybridization events are on the rise in natural systems due to climate change disrupting species barriers. Across taxa, microsatellites have long been the molecular markers of choice to identify admixed individuals. However, with the advent of high-throughput sequencing easing the generation of genome-wide datasets, incorrect reports of hybridization resulting from microsatellite technical artefacts have been uncovered in a growing number of taxa. In the marine zooplankton genus Calanus (Copepoda), whose species are used as climate change indicators, microsatellite markers have suggested hybridization between C. finmarchicus and C. glacialis, while other nuclear markers (InDels) never detected any admixed individuals, leaving the scientific community divided. Here, for the first time, we investigated the potential for hybridization among C. finmarchicus, C. glacialis, C. helgolandicus and C. hyperboreus using two large and independent SNP datasets. These were derived firstly from a protocol of target-capture applied to 179 individuals collected from 17 sites across the North Atlantic and Arctic Oceans, including sympatric areas, and second from published RNA sequences. All SNP-based analyses were congruent in showing that Calanus species are distinct and do not appear to hybridize. We then thoroughly re-assessed the microsatellites showing hybrids, with the support of published transcriptomes, and identified technical issues plaguing eight out of 10 microsatellites, including size homoplasy, paralogy, potential for null alleles and even two primer pairs targeting the same locus. Our study illustrates how deceptive microsatellites can be when applied to the investigation of hybridization.

Genome-wide variant analyses reveal new patterns of admixture and population structure in Australian dingoes

Admixture between species is a cause for concern in wildlife management. Canids are particularly vulnerable to interspecific hybridisation, and genetic admixture has shaped their evolutionary history. Microsatellite DNA testing, relying on a small number of genetic markers and geographically restricted reference populations, has identified extensive domestic dog admixture in Australian dingoes and driven conservation management policy. But there exists a concern that geographic variation in dingo genotypes could confound ancestry analyses that use a small number of genetic markers. Here, we apply genome-wide single-nucleotide polymorphism (SNP) genotyping to a set of 402 wild and captive dingoes collected from across Australia and then carry out comparisons to domestic dogs. We then perform ancestry modelling and biogeographic analyses to characterise population structure in dingoes and investigate the extent of admixture between dingoes and dogs in different regions of the continent. We show that there are at least five distinct dingo populations across Australia. We observed limited evidence of dog admixture in wild dingoes. Our work challenges previous reports regarding the occurrence and extent of dog admixture in dingoes, as our ancestry analyses show that previous assessments severely overestimate the degree of domestic dog admixture in dingo populations, particularly in south-eastern Australia. These findings strongly support the use of genome-wide SNP genotyping as a refined method for wildlife managers and policymakers to assess and inform dingo management policy and legislation moving forwards.

Genomic and machine learning-based screening of aquaculture-associated introgression into at-risk wild North American Atlantic salmon (Salmo salar) populations

The negative genetic impacts of gene flow from domestic to wild populations can be dependent on the degree of domestication and exacerbated by the magnitude of pre-existing genetic differences between wild populations and the domestication source. Recent evidence of European ancestry within North American aquaculture Atlantic salmon (Salmo salar) has elevated the potential impact of escaped farmed salmon on often at-risk wild North American salmon populations. Here, we compare the ability of single nucleotide polymorphism (SNP) and microsatellite (SSR) marker panels of different sizes (7-SSR, 100-SSR and 220K-SNP) to detect introgression of European genetic information into North American wild and aquaculture populations. Linear regression comparing admixture predictions for a set of individuals common to the three datasets showed that the 100-SSR panel and 7-SSR panels replicated the full 220K-SNP-based admixture estimates with low accuracy (r² of .64 and .49, respectively). Additional tests explored the effects of individual sample size and marker number, which revealed that ~300 randomly selected SNPs could replicate the 220K-SNP admixture predictions with greater than 95% fidelity. We designed a custom SNP panel (301-SNP) for European admixture detection in future monitoring work and then developed and tested a python package, salmoneuadmix (https://github.com/CNuge/SalmonEuAdmix), which uses a deep neural network to make de novo estimates of individuals' European admixture proportion without the need to conduct complete admixture analysis utilizing baseline samples. The results demonstrate the mobilization of targeted SNP panels and machine learning in support of at-risk species conservation and management.

Detection and sequence analysis of Canine morbillivirus in multiple species of the Mustelidae family

BACKGROUND

Canine morbillivirus (canine distemper virus, CDV) is a member of the Paramyxoviridae family. Canine distemper is a serious viral disease that affects many mammalian species, including members of the Mustelidae family. These animals have an elusive nature, which makes related virological studies extremely challenging. There is a significant knowledge gap about the evolution of their viruses and about the possible effects of these viruses to the population dynamics of the host animals. Spleen and lung tissue samples of 170 road-killed mustelids belonging to six species were collected between 1997 and 2022 throughout Hungary and tested for CDV with real-time RT-PCR.

RESULTS

Three species were positive for viral RNA, 2 out of 64 Steppe polecats (Mustela eversmanii), 1 out of 36 European polecats (Mustela putorius) and 2 out of 36 stone martens (Martes foina); all 18 pine martens (Martes martes), 10 least weasels (Mustela nivalis) and 6 stoats (Mustela erminea) tested negative. The complete CDV genome was sequenced in five samples using pan-genotype CDV-specific, amplicon-based Nanopore sequencing. Based on the phylogenetic analysis, all five viral sequences were grouped to the Europe/South America 1 lineage and the distribution of one sequence among trees indicated recombination of the Hemagglutinin gene. We verified the recombination with SimPlot analysis.

CONCLUSIONS

This paper provides the first CDV genome sequences from Steppe polecats and additional complete genomes from European polecats and stone martens. The infected specimens of various species originated from distinct parts of the country over a long time, indicating a wide circulation of CDV among mustelids throughout Hungary. Considering the high virulence of CDV and the presence of the virus in these animals, we highlight the importance of conservation efforts for wild mustelids. In addition, we emphasize the importance of full genomic data acquisition and analysis to better understand the evolution of the virus. Since CDV is prone to recombination, specific genomic segment analyses may provide less representative evolutionary traits than using complete genome sequences.

Extensive genome introgression between domestic ferret and European polecat during population recovery in Great Britain

The European polecat (Mustela putorius) is a mammalian predator which occurs across much of Europe east to the Ural Mountains. In Great Britain, following years of persecution the range of the European polecat contracted and by the early 1900s was restricted to unmanaged forests of central Wales. The European polecat has recently undergone a population increase due to legal protection and its range now overlaps that of feral domestic ferrets (Mustela putorius furo). During this range expansion, European polecats hybridized with feral domestic ferrets producing viable offspring. Here, we carry out population-level whole-genome sequencing on 8 domestic ferrets, 19 British European polecats, and 15 European polecats from the European mainland. We used a range of population genomics methods to examine the data, including phylogenetics, phylogenetic graphs, model-based clustering, phylogenetic invariants, ABBA-BABA tests, topology weighting, and Fst. We found high degrees of genome introgression in British polecats outside their previous stronghold, even in those individuals phenotyped as “pure” polecats. These polecats ranged from presumed F1 hybrids (gamma = 0.53) to individuals that were much less introgressed (gamma = 0.2). We quantify this introgression and find introgressed genes containing Fst outliers associated with cognitive function and sight.