A single‐nucleotide polymorphism‐based approach for rapid and cost‐effective genetic wolf monitoring in Europe based on noninvasively collected samples

Noninvasive genetics based on microsatellite markers has become an indispensable tool for wildlife monitoring and conservation research over the past decades. However, microsatellites have several drawbacks, such as the lack of standardisation between laboratories and high error rates. Here, we propose an alternative single‐nucleotide polymorphism (SNP)‐based marker system for noninvasively collected samples, which promises to solve these problems. Using nanofluidic SNP genotyping technology (Fluidigm), we genotyped 158 wolf samples (tissue, scats, hairs, urine) for 192 SNP loci selected from the Affymetrix v2 Canine SNP Array. We carefully selected an optimised final set of 96 SNPs (and discarded the worse half), based on assay performance and reliability. We found rates of missing data in this SNP set of <10% and genotyping error of ~1%, which improves genotyping accuracy by nearly an order of magnitude when compared to published data for other marker types. Our approach provides a tool for rapid and cost‐effective genotyping of noninvasively collected wildlife samples. The ability to standardise genotype scoring combined with low error rates promises to constitute a major technological advancement and could establish SNPs as a standard marker for future wildlife monitoring.     

Estimation of breeding population size using DNA‐based pedigree reconstruction in brown bears

Robust estimates of demographic parameters are critical for effective wildlife conservation and management but are difficult to obtain for elusive species. We estimated the breeding and adult population sizes, as well as the minimum population size, in a high‐density brown bear population on the Shiretoko Peninsula, in Hokkaido, Japan, using DNA‐based pedigree reconstruction. A total of 1288 individuals, collected in and around the Shiretoko Peninsula between 1998 and 2020, were genotyped at 21 microsatellite loci. Among them, 499 individuals were identified by intensive genetic sampling conducted in two consecutive years (2019 and 2020) mainly by noninvasive methods (e.g., hair and fecal DNA). Among them, both parents were assigned for 330 bears, and either maternity or paternity was assigned to 47 and 76 individuals, respectively. The subsequent pedigree reconstruction indicated a range of breeding and adult (≥4 years old) population sizes: 128–173 for female breeders and 66–91 male breeders, and 155–200 for female adults and 84–109 male adults. The minimum population size was estimated to be 449 (252 females and 197 males) in 2019. Long‐term continuous genetic sampling prior to a short‐term intensive survey would enable parentage to be identified in a population with a high probability, thus enabling reliable estimates of breeding population size for elusive species.     

Single nucleotide polymorphism‐based dispersal estimates using noninvasive sampling

Quantifying dispersal within wild populations is an important but challenging task. Here we present a method to estimate contemporary, individual‐based dispersal distance from noninvasively collected samples using a specialized panel of 96 SNPs (single nucleotide polymorphisms). One main issue in conducting dispersal studies is the requirement for a high sampling resolution at a geographic scale appropriate for capturing the majority of dispersal events. In this study, fecal samples of brown bear (Ursus arctos) were collected by volunteer citizens, resulting in a high sampling resolution spanning over 45,000 km² in Gävleborg and Dalarna counties in Sweden. SNP genotypes were obtained for unique individuals sampled (n = 433) and subsequently used to reconstruct pedigrees. A Mantel test for isolation by distance suggests that the sampling scale was appropriate for females but not for males, which are known to disperse long distances. Euclidean distance was estimated between mother and offspring pairs identified through the reconstructed pedigrees. The mean dispersal distance was 12.9 km (SE 3.2) and 33.8 km (SE 6.8) for females and males, respectively. These results were significantly different (Wilcoxon's rank‐sum test: P‐value = 0.02) and are in agreement with the previously identified pattern of male‐biased dispersal. Our results illustrate the potential of using a combination of noninvasively collected samples at high resolution and specialized SNPs for pedigree‐based dispersal models.     

Using pedigree reconstruction to estimate population size: genotypes are more than individually unique marks

Estimates of population size are critical for conservation and management, but accurate estimates are difficult to obtain for many species. Noninvasive genetic methods are increasingly used to estimate population size, particularly in elusive species such as large carnivores, which are difficult to count by most other methods. In most such studies, genotypes are treated simply as unique individual identifiers. Here, we develop a new estimator of population size based on pedigree reconstruction. The estimator accounts for individuals that were directly sampled, individuals that were not sampled but whose genotype could be inferred by pedigree reconstruction, and individuals that were not detected by either of these methods. Monte Carlo simulations show that the population estimate is unbiased and precise if sampling is of sufficient intensity and duration. Simulations also identified sampling conditions that can cause the method to overestimate or underestimate true population size; we present and discuss methods to correct these potential biases. The method detected 2–21% more individuals than were directly sampled across a broad range of simulated sampling schemes. Genotypes are more than unique identifiers, and the information about relationships in a set of genotypes can improve estimates of population size.     

Phylogeography and adaptation genetics of stickleback from the Haida Gwaii archipelago revealed using genome‐wide single nucleotide polymorphism genotyping

Threespine stickleback populations are model systems for studying adaptive evolution and the underlying genetics. In lakes on the Haida Gwaii archipelago (off western Canada), stickleback have undergone a remarkable local radiation and show phenotypic diversity matching that seen throughout the species distribution. To provide a historical context for this radiation, we surveyed genetic variation at >1000 single nucleotide polymorphism (SNP) loci in stickleback from over 100 populations. SNPs included markers evenly distributed throughout genome and candidate SNPs tagging adaptive genomic regions. Based on evenly distributed SNPs, the phylogeographic pattern differs substantially from the disjunct pattern previously observed between two highly divergent mtDNA lineages. The SNP tree instead shows extensive within watershed population clustering and different watersheds separated by short branches deep in the tree. These data are consistent with separate colonizations of most watersheds, despite underlying genetic connections between some independent drainages. This supports previous suppositions that morphological diversity observed between watersheds has been shaped independently, with populations exhibiting complete loss of lateral plates and giant size each occurring in several distinct clades. Throughout the archipelago, we see repeated selection of SNPs tagging candidate freshwater adaptive variants at several genomic regions differentiated between marine–freshwater populations on a global scale (e.g. EDA, Na/K ATPase). In estuarine sites, both marine and freshwater allelic variants were commonly detected. We also found typically marine alleles present in a few freshwater lakes, especially those with completely plated morphology. These results provide a general model for postglacial colonization of freshwater habitat by sticklebacks and illustrate the tremendous potential of genome‐wide SNP data sets hold for resolving patterns and processes underlying recent adaptive divergences.     

Canadian polar bear population structure using genome‐wide markers

Predicting the consequences of environmental changes, including human‐mediated climate change on species, requires that we quantify range‐wide patterns of genetic diversity and identify the ecological, environmental, and historical factors that have contributed to it. Here, we generate baseline data on polar bear population structure across most Canadian subpopulations (n = 358) using 13,488 genome‐wide single nucleotide polymorphisms (SNPs) identified with double‐digest restriction site‐associated DNA sequencing (ddRAD). Our ddRAD dataset showed three genetic clusters in the sampled Canadian range, congruent with previous studies based on microsatellites across the same regions; however, due to a lack of sampling in Norwegian Bay, we were unable to confirm the existence of a unique cluster in that subpopulation. These data on the genetic structure of polar bears using SNPs provide a detailed baseline against which future shifts in population structure can be assessed, and opportunities to develop new noninvasive tools for monitoring polar bears across their range.     

A single‐nucleotide polymorphism of miR‐196a‐2 and vitiligo: an association study and functional analysis in a Han Chinese population

Recent evidence indicates that oxidative stress and genetic factors play an important role in the pathogenesis of vitiligo. SNPs in miRNAs involved in oxidative stress could potentially influence the development of vitiligo. In this case–control study, we investigated the association of a functional SNP of rs11614913 in miR‐196a‐2 with risk of vitiligo. A significantly lower risk of vitiligo was associated with the rs11614913 miR‐196a‐2 CC genotype (adjusted OR, 0.77; CI, 0.60–0.98). In addition, TYRP1 gene expression was considerably down‐regulated by the rs11614913 C allele in miR‐196a‐2, which lowered the levels of intracellular reactive oxygen species (ROS) and reduced the proportion of early apoptosis in human melanocytes in response to H₂O₂ treatment. Our data suggest that the rs11614913 C allele in miR‐196a‐2 confers potential protection against oxidative effects on human melanocytes through the modulation of the target gene, TYRP1, which may account for the decreased risk of vitiligo in this study population.     

A single‐nucleotide polymorphism of miR‐146a and psoriasis: an association and functional study

Epidermal growth factor receptor (EGFR), which is overexpressed in psoriatic lesions, has been proven to contribute to the hyperproliferation of keratinocytes in psoriasis. Single nucleotide polymorphisms (SNPs) involved in miRNAs that can regulate the expression of EGFR could potentially influence the development of psoriasis. The present study investigated the association between a functional SNP of rs2910164 in miR‐146a and the risk of psoriasis in the Chinese Han population. A total of 521 Han Chinese patients with psoriasis and 582 healthy controls were recruited in this study. The miR‐146a rs2910164 SNP was genotyped by polymerase chain reaction‐restriction fragment length polymorphism. Overall, a significantly increased risk of psoriasis was associated with the rs2910164 miR‐146a CG and GG genotypes (adjusted OR, 1.38; 95% CI, 1.06–1.80). Furthermore, the rs2910164G allele in miR‐146a attenuated its inhibitory regulation on the expression of EGFR as well as the proliferation of human keratinocytes, and lowered the level of miR‐146a in the psoriatic lesions. These findings indicate that the rs2910164G allele in miR‐146a weakens its suppression on the proliferation of keratinocytes probably through the decreased inhibition of the target gene, EGFR, which may account for the increased risk of psoriasis in this study population.     

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single‐nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole‐genome resequencing data from population samples, molecular ecology studies will be able to take truly genome‐wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single‐nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10‐fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later‐generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.     

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over‐exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded F_ST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between‐ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17‐SNP panel was developed in Atlantic BFT by cross‐species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (F_ST = 0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.     

Assessing performance of single‐sample molecular genetic methods to estimate effective population size: empirical evidence from the endangered Gochu Asturcelta pig breed

Estimating effective population size (N_e) using linkage disequilibrium (LD) information (N_e(LD)) has the operational advantage of using a single sample. However, N_e(LD) estimates assume discrete generations and its performance are constrained by demographic issues. However, such concerns have received little empirical attention so far. The pedigree of the endangered Gochu Asturcelta pig breed includes individuals classified into discrete filial generations and individuals with generations overlap. Up to 780 individuals were typed with a set of 17 microsatellites. Performance of N_e(LD) was compared with N_e estimates obtained using genealogical information, molecular coancestry (N_e(M)) and a temporal (two‐sample) method (N_e(JR)). Molecular‐based estimates of N_e exceeded those obtained using pedigree data. Estimates of N_e(LD) for filial generations F3 and F4 (17.0 and 17.3, respectively) were lower and steadier than those obtained using yearly or biannual samplings. N_e(LD) estimated for samples including generations overlap could only be compared with those obtained for the discrete filial generations when sampling span approached a generation interval and demographic correction for bias was applied. Single‐sample N_e(M) estimates were lower than their N_e(LD) counterparts. N_e(M) estimates are likely to partially reflect the number of founders rather than population size. In any case, estimates of LD and molecular coancestry tend to covary and, therefore, N_e(M) and N_e(LD) can hardly be considered independent. Demographically adjusted estimates of N_e(JR) and N_e(LD) took comparable values when: (1) the two samples used for the former were separated by one equivalent to discrete generations in the pedigree and (2) sampling span used for the latter approached a generation interval. Overall, the empirical evidence given in this study suggested that the advantage of using single‐sample methods to obtain molecular‐based estimates of N_e is not clear in operational terms. Estimates of N_e obtained using methods based in molecular information should be interpreted with caution.     

Assessing polar bear (Ursus maritimus) population structure in the Hudson Bay region using SNPs

Defining subpopulations using genetics has traditionally used data from microsatellite markers to investigate population structure; however, single‐nucleotide polymorphisms (SNPs) have emerged as a tool for detection of fine‐scale structure. In Hudson Bay, Canada, three polar bear (Ursus maritimus) subpopulations (Foxe Basin (FB), Southern Hudson Bay (SH), and Western Hudson Bay (WH)) have been delineated based on mark–recapture studies, radiotelemetry and satellite telemetry, return of marked animals in the subsistence harvest, and population genetics using microsatellites. We used SNPs to detect fine‐scale population structure in polar bears from the Hudson Bay region and compared our results to the current designations using 414 individuals genotyped at 2,603 SNPs. Analyses based on discriminant analysis of principal components (DAPC) and STRUCTURE support the presence of four genetic clusters: (i) Western—including individuals sampled in WH, SH (excluding Akimiski Island in James Bay), and southern FB (south of Southampton Island); (ii) Northern—individuals sampled in northern FB (Baffin Island) and Davis Strait (DS) (Labrador coast); (iii) Southeast—individuals from SH (Akimiski Island in James Bay); and (iv) Northeast—individuals from DS (Baffin Island). Population structure differed from microsatellite studies and current management designations demonstrating the value of using SNPs for fine‐scale population delineation in polar bears.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Geographic structure in the Southern Ocean circumpolar brittle star Ophionotus victoriae (Ophiuridae) revealed from mtDNA and single‐nucleotide polymorphism data

Marine systems have traditionally been thought of as “open” with few barriers to gene flow. In particular, many marine organisms in the Southern Ocean purportedly possess circumpolar distributions that have rarely been well verified. Here, we use the highly abundant and endemic Southern Ocean brittle star Ophionotus victoriae to examine genetic structure and determine whether barriers to gene flow have existed around the Antarctic continent. Ophionotus victoriae possesses feeding planktotrophic larvae with presumed high dispersal capability, but a previous study revealed genetic structure along the Antarctic Peninsula. To test the extent of genetic differentiation within O. victoriae, we sampled from the Ross Sea through the eastern Weddell Sea. Whereas two mitochondrial DNA markers (16S rDNA and COI) were employed to allow comparison to earlier work, a 2b‐RAD single‐nucleotide polymorphism (SNP) approach allowed sampling of loci across the genome. Mitochondrial data from 414 individuals suggested three major lineages, but 2b‐RAD data generated 1,999 biallelic loci that identified four geographically distinct groups from 89 samples. Given the greater resolution by SNP data, O. victoriae can be divided into geographically distinct populations likely representing multiple species. Specific historical scenarios that explain current population structure were examined with approximate Bayesian computation (ABC) analyses. Although the Bransfield Strait region shows high diversity possibly due to mixing, our results suggest that within the recent past, dispersal processes due to strong currents such as the Antarctic Circumpolar Current have not overcome genetic subdivision presumably due to historical isolation, questioning the idea of large open circumpolar populations in the Southern Ocean.     

Single‐nucleotide polymorphism discovery and diversity in the model legume Medicago truncatula

Extensive genomic resources are available in the model legume Medicago truncatula. Here, we present the discovery and design of the first array of single‐nucleotide polymorphism (SNP) markers in M. truncatula through large‐scale Sanger resequencing of genomic fragments spanning the genome, in a diverse panel of 16 M. truncatula accessions. Both anonymous fragments and fragments targeting candidate genes for flowering phenology and symbiosis were surveyed for nucleotide variation in almost 230 kb of unique genomic regions. A set of 384 SNP markers was designed for an Illumina's GoldenGate assay, genotyped on a collection of 192 inbred lines (CC192) representing the geographical range of the species and used to survey the diversity of two natural populations. Finally, 86% of the tested SNPs were of high quality and exhibited polymorphism in the CC192 collection. Even at the population level, we detected polymorphism for more than 50% of the selected SNPs. Analysis of the allele frequency spectrum in the CC192 showed a reduced ascertainment bias, mostly limited to very rare alleles (frequency <0.01). The substantial polymorphism detected at the species and population levels, the high marker quality and the potential to survey large samples of individuals make this set of SNP markers a valuable tool to improve our understanding of the effect of demographic and selective factors that shape the natural genetic diversity within the selfing species Medicago truncatula.     

Exploring a Pool‐seq‐only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool‐seq data to generate a de novo genome assembly for mining exons, upon which Pool‐seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual‐based single nucleotide polymorphisms [SNPs]) and from mapping the Pool‐seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F_ST) between the two introduced populations exceeds that of the naturally sympatric populations (F_ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high‐quality reference assembly from a divergent species. We conclude that the Pool‐seq‐only approach can be suitable for detecting and quantifying genome‐wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.     

Assessment of coyote–wolf–dog admixture using ancestry‐informative diagnostic SNPs

The evolutionary importance of hybridization as a source of new adaptive genetic variation is rapidly gaining recognition. Hybridization between coyotes and wolves may have introduced adaptive alleles into the coyote gene pool that facilitated an expansion in their geographic range and dietary niche. Furthermore, hybridization between coyotes and domestic dogs may facilitate adaptation to human‐dominated environments. We genotyped 63 ancestry‐informative single‐nucleotide polymorphisms in 427 canids to examine the prevalence, spatial distribution and the ecology of admixture in eastern coyotes. Using multivariate methods and Bayesian clustering analyses, we estimated the relative contributions of western coyotes, western and eastern wolves, and domestic dogs to the admixed ancestry of Ohio and eastern coyotes. We found that eastern coyotes form an extensive hybrid swarm, with all our samples having varying levels of admixture. Ohio coyotes, previously thought to be free of admixture, are also highly admixed with wolves and dogs. Coyotes in areas of high deer density are genetically more wolf‐like, suggesting that natural selection for wolf‐like traits may result in local adaptation at a fine geographic scale. Our results, in light of other previously published studies of admixture in Canis, revealed a pattern of sex‐biased hybridization, presumably generated by male wolves and dogs mating with female coyotes. This study is the most comprehensive genetic survey of admixture in eastern coyotes and demonstrates that the frequency and scope of hybridization can be quantified with relatively few ancestry‐informative markers.     

A method for estimating population sex ratio for sage‐grouse using noninvasive genetic samples

Population sex ratio is an important metric for wildlife management and conservation, but estimates can be difficult to obtain, particularly for sexually monomorphic species or for species that differ in detection probability between the sexes. Noninvasive genetic sampling (NGS) using polymerase chain reaction (PCR) has become a common method for identifying sex from sources such as hair, feathers or faeces, and is a potential source for estimating sex ratio. If, however, PCR success is sex‐biased, naively using NGS could lead to a biased sex ratio estimator. We measured PCR success rates and error rates for amplifying the W and Z chromosomes from greater sage‐grouse (Centrocercus urophasianus) faecal samples, examined how success and error rates for sex identification changed in response to faecal sample exposure time, and used simulation models to evaluate precision and bias of three sex assignment criteria for estimating population sex ratio with variable sample sizes and levels of PCR replication. We found PCR success rates were higher for females than males and that choice of sex assignment criteria influenced the bias and precision of corresponding sex ratio estimates. Our simulations demonstrate the importance of considering the interplay between the sex bias of PCR success, number of genotyping replicates, sample size, true population sex ratio and accuracy of assignment rules for designing future studies. Our results suggest that using faecal DNA for estimating the sex ratio of sage‐grouse populations has great potential and, with minor adaptations and similar marker evaluations, should be applicable to numerous species.     

A genomewide catalogue of single nucleotide polymorphisms in white‐beaked and Atlantic white‐sided dolphins

The field of population genetics is rapidly moving into population genomics as the quantity of data generated by high‐throughput sequencing platforms increases. In this study, we used restriction‐site‐associated DNA sequencing (RADSeq) to recover genomewide genotypes from 70 white‐beaked (Lagenorhynchus albirostris) and 43 Atlantic white‐sided dolphins (L. acutus) gathered throughout their north‐east Atlantic distribution range. Both species are at a high risk of being negatively affected by climate change. Here, we provide a resource of 38 240 RAD‐tags and 52 981 nuclear SNPs shared between both species. We have estimated overall higher levels of nucleotide diversity in white‐sided (π = 0.0492 ± 0.0006%) than in white‐beaked dolphins (π = 0.0300 ± 0.0004%). White‐sided dolphins sampled in the Faroe Islands, belonging to two pods (N = 7 and N = 11), showed similar levels of diversity (π = 0.0317 ± 0.0007% and 0.0267 ± 0.0006%, respectively) compared to unrelated individuals of the same species sampled elsewhere (e.g. π = 0.0285 ± 0.0007% for 11 Scottish individuals). No evidence of higher levels of kinship within pods can be derived from our analyses. When identifying the most likely number of genetic clusters among our sample set, we obtained an estimate of two to four clusters, corresponding to both species and possibly, two further clusters within each species. A higher diversity and lower population structuring was encountered in white‐sided dolphins from the north‐east Atlantic, in line with their preference for pelagic waters, as opposed to white‐beaked dolphins that have a more patchy distribution, mainly across continental shelves.