A phylogenomic perspective on diversity,hybridization and evolutionary affinities in the stickleback genus Pungitius

Hybridization and convergent evolution are phenomena of broad interest in evolutionary biology, but their occurrence poses challenges for reconstructing evolutionary affinities among affected taxa. Sticklebacks in the genus Pungitius are a case in point: evolutionary relationships and taxonomic validity of different species and populations in this circumpolarly distributed species complex remain contentious due to convergent evolution of traits regarded as diagnostic in their taxonomy, and possibly also due to frequent hybridization among taxa. To clarify the evolutionary relationships among different Pungitius species and populations globally, as well as to study the prevalence and extent of introgression among recognized species, genomic data sets of both reference genome‐anchored single nucleotide polymorphisms and de novo assembled RAD‐tag loci were constructed with RAD‐seq data. Both data sets yielded topologically identical and well‐supported species trees. Incongruence between nuclear and mitochondrial DNA‐based trees was found and suggested possibly frequent hybridization and mitogenome capture during the evolution of Pungitius sticklebacks. Further analyses revealed evidence for frequent nuclear genetic introgression among Pungitius species, although the estimated proportions of autosomal introgression were low. Apart from providing evidence for frequent hybridization, the results challenge earlier mitochondrial and morphology‐based hypotheses regarding the number of species and their affinities in this genus: at least seven extant species can be recognized on the basis of genetic data. The results also shed new light on the biogeographical history of the Pungitius‐complex, including suggestion of several trans‐Arctic invasions of Europe from the Northern Pacific. The well‐resolved phylogeny should facilitate the utility of this genus as a model system for future comparative evolutionary studies.     

Population structure of Atlantic mackerel inferred from RAD‐seq‐derived SNP markers: effects of sequence clustering parameters and hierarchical SNP selection

Restriction‐site‐associated DNA sequencing (RAD‐seq) and related methods are revolutionizing the field of population genomics in nonmodel organisms as they allow generating an unprecedented number of single nucleotide polymorphisms (SNPs) even when no genomic information is available. Yet, RAD‐seq data analyses rely on assumptions on nature and number of nucleotide variants present in a single locus, the choice of which may lead to an under‐ or overestimated number of SNPs and/or to incorrectly called genotypes. Using the Atlantic mackerel (Scomber scombrus L.) and a close relative, the Atlantic chub mackerel (Scomber colias), as case study, here we explore the sensitivity of population structure inferences to two crucial aspects in RAD‐seq data analysis: the maximum number of mismatches allowed to merge reads into a locus and the relatedness of the individuals used for genotype calling and SNP selection. Our study resolves the population structure of the Atlantic mackerel, but, most importantly, provides insights into the effects of alternative RAD‐seq data analysis strategies on population structure inferences that are directly applicable to other species.     

Parallel genetic divergence among coastal–marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome‐wide data set of RAD‐derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal–marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine‐scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

PMERGE: Computational filtering of paralogous sequences from RAD‐seq data

Restriction‐site associated DNA sequencing (RAD‐seq) can identify and score thousands of genetic markers from a group of samples for population‐genetics studies. One challenge of de novo RAD‐seq analysis is to distinguish paralogous sequence variants (PSVs) from true single‐nucleotide polymorphisms (SNPs) associated with orthologous loci. In the absence of a reference genome, it is difficult to differentiate true SNPs from PSVs, and their impact on downstream analysis remains unclear. Here, we introduce a network‐based approach, PMERGE that connects fragments based on their DNA sequence similarity to identify probable PSVs. Applying our method to de novo RAD‐seq data from 150 Atlantic salmon (Salmo salar) samples collected from 15 locations across the Southern Newfoundland coast allowed the identification of 87% of total PSVs identified through alignment to the Atlantic salmon genome. Removal of these paralogs altered the inferred population structure, highlighting the potential impact of filtering in RAD‐seq analysis. PMERGE is also applied to a green crab (Carcinus maenas) data set consisting of 242 samples from 11 different locations and was successfully able to identify and remove the majority of paralogous loci (62%). The PMERGE software can be run as part of the widely used Stacks analysis package.     

Spatial and temporal genetic structure of the planktonic Sagitta setosa (Chaetognatha) in European seas as revealed by mitochondrial and nuclear DNA markers

Little is known about the spatial and temporal scales at which planktonic organisms are genetically structured. A previous study of mitochondrial DNA (mtDNA) in the holoplanktonic chaetognath Sagitta setosa revealed strong phylogeographic structuring suggesting that Northeast (NE) Atlantic, Mediterranean and Black Sea populations are genetically disjunct. The present study used a higher sampling intensity and a combination of mitochondrial and four microsatellite markers to reveal population structuring between and within basins. Between basins, both marker sets indicated significant differentiation confirming earlier results that gene flow is probably absent between the respective S. setosa populations. At the within-basin scale, we found no evidence of spatial or temporal structuring within the NE Atlantic. In the Mediterranean basin, both marker sets indicated significant structuring, but only the mtDNA data indicated a sharp genetic division between Adriatic and all other Mediterranean populations. Data were inconclusive about population structuring in the Black Sea. The levels of differentiation indicated by the two marker sets differed substantially, with far less pronounced structure detected by microsatellite than mtDNA data. This study also uncovered the presence of highly divergent mitochondrial lineages that were discordant with morphology, geography and nuclear DNA. We thus propose the hypothesis that highly divergent mitochondrial lineages may be present within interbreeding S. setosa populations.     

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M‐lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C‐lineage). Effective conservation actions require reliable molecular tools to identify pure‐bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q‐values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.     

Chromosome polymorphisms track trans‐Atlantic divergence and secondary contact in Atlantic salmon

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans‐Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans‐Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans‐Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post‐glacial trans‐Atlantic secondary contact, contrasting with low European ancestry genome‐wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large‐scale intraspecific variation in genomic architecture of northern species.     

Range‐wide population structure of European sea bass Dicentrarchus labrax

The euryhaline European sea bass Dicentrarchus labrax L., inhabiting the coasts of the eastern Atlantic Ocean and Mediterranean Sea, has had many opportunities for differentiation throughout its large natural range. However, evidence for this has been incompletely documented geographically and with an insufficient number of markers. Therefore, its full range was sampled at 22 sites and individuals were genotyped with a suite of mapped markers, including 14 microsatellite loci (N = 536) and 46 neutral or gene‐linked single nucleotide polymorphisms (SNPs; N = 644). We confirm that the Atlantic and Mediterranean basins harbour two distinct lineages. Within the Atlantic Ocean no pattern was obvious based on the microsatellite and SNP genotypes, except for a subtle difference between South‐eastern and North‐eastern Atlantic sea bass attributed to limited introgression of alleles of Mediterranean origin. SNP genotypes of the Mediterranean lineage differentiated into three groups, probably under the influence of geographical isolation. The Western Mediterranean group showed genetic homogeneity without evidence for outlier loci. The Adriatic group appeared as a distinct unit. The Eastern Mediterranean group showed a longitudinal gradient of genotypes and most interestingly an outlier locus linked to the somatolactin gene. Overall, the spatial pattern fits those observed with other taxa of between‐basin segregation and within‐basin connectivity, which concurs well with the swimming capabilities of European sea bass. Evidence from a few outlier loci in this and other studies encourages further exploration of its regional connectivity and adaptive evolution.     

Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing

A major barrier to evolutionary studies of sex determination and sex chromosomes has been a lack of information on the types of sex‐determining mechanisms that occur among different species. This is particularly problematic in groups where most species lack visually heteromorphic sex chromosomes, such as fish, amphibians and reptiles, because cytogenetic analyses will fail to identify the sex chromosomes in these species. We describe the use of restriction site‐associated DNA (RAD) sequencing, or RAD‐seq, to identify sex‐specific molecular markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the lizard Anolis carolinensis. We performed RAD‐seq on seven male and ten female A. carolinensis and found one male‐specific molecular marker. Anolis carolinensis has previously been shown to possess male heterogamety and the recently published A. carolinensis genome facilitated the characterization of the sex‐specific RAD‐seq marker. We validated the male specificity of the new marker using PCR on additional individuals and also found that it is conserved in some other Anolis species. We discuss the utility of using RAD‐seq to identify sex‐determining mechanisms in other species with cryptic or homomorphic sex chromosomes and the implications for the evolution of male heterogamety in Anolis.     

Defining conservation units in a stocking‐induced genetic melting pot: unraveling native and multiple exotic genetic imprints of recent and historical secondary contact in Adriatic grayling

The definition of conservation units is crucial for the sustainable management of endangered species, though particularly challenging when recent and past anthropogenic and natural gene flow might have played a role. The conservation of the European grayling, Thymallus thymallus, is particularly complex in its southern distribution area, where the Adriatic evolutionary lineage is endangered by a long history of anthropogenic disturbance, intensive stocking and potentially widespread genetic introgression. We provide mtDNA sequence and microsatellite data of 683 grayling from 30 sites of Adriatic as well as Danubian and Atlantic origin. We apply Bayesian clustering and Approximate Bayesian Computation (ABC) to detect microgeographic population structure and to infer the demographic history of the Adriatic populations, to define appropriate conservation units. Varying frequencies of indigenous genetic signatures of the Adriatic grayling were revealed, spanning from marginal genetic introgression to the collapse of native gene pools. Genetic introgression involved multiple exotic source populations of Danubian and Atlantic origin, thus evidencing the negative impact of few decades of stocking. Within the Adige River system, a contact zone of western Adriatic and eastern Danubian populations was detected, with ABC analyses suggesting a historical anthropogenic origin of eastern Adige populations, most likely founded by medieval translocations. Substantial river‐specific population substructure within the Adriatic grayling Evolutionary Significant Unit points to the definition of different conservation units. We finally propose a catalog of management measures, including the legal prohibition of stocking exotic grayling and the use of molecular markers in supportive‐ and captive‐breeding programs.     

Impact of climate change and human‐mediated introgression on southern European Atlantic salmon populations

This study focuses on temporal changes in Atlantic salmon (Salmo salar) populations from the vulnerable periphery of the species range (northern Spain). Using microsatellite markers to assess population structuring and introgression of exogenous genes in four different temporal samples collected across 20 years, we have determined the relative weights of climate and stocking practices in shaping contemporary regional population genetic patterns. Climate, represented by the North Atlantic Oscillation Index, was identified as the main factor for determining the level of population genetic differentiation. Populations within the region have become homogenized through gene flow enhanced by straying of adult salmon from natal rivers and subsequent interchange of genes among rivers due to warmer temperatures. At the same time, and in line with documented changes in stock transfer strategies, evidence of genetic introgression from past stock transfers has decreased throughout the study period, becoming a secondary factor in erasing population structuring. The ability to disentangle the effects of climatic changes and anthropogenic factors (fisheries management practices) is essential for effective long‐term conservation of this iconic species. We emphasize the importance of evaluating all factors which may be linked to stocking practices in vulnerable species, particularly those sensitive to climate change.     

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site‐associated DNA sequencing (RAD‐Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F_ST 0.008; range 0.00–0.018) were concordant with previous microsatellite‐based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age‐1 and age‐2 families of full and half‐siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non‐native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.     

Attack of the PCR clones: Rates of clonality have little effect on RAD‐seq genotype calls

Interpretation of high‐throughput sequence data requires an understanding of how decisions made during bioinformatic data processing can influence results. One source of bias that is often cited is PCR clones (or PCR duplicates). PCR clones are common in restriction site‐associated sequencing (RAD‐seq) data sets, which are increasingly being used for molecular ecology. To determine the influence PCR clones and the bioinformatic handling of clones have on genotyping, we evaluate four RAD‐seq data sets. Data sets were compared before and after clones were removed to estimate the number of clones present in RAD‐seq data, quantify how often the presence of clones in a data set causes genotype calls to change compared to when clones were removed, investigate the mechanisms that lead to genotype call changes and test whether clones bias heterozygosity estimates. Our RAD‐seq data sets contained 30%–60% PCR clones, but 95% of RAD‐tags had five or fewer clones. Relatively few genotypes changed once clones were removed (5%–10%), and the vast majority of these changes (98%) were associated with genotypes switching from a called to no‐call state or vice versa. PCR clones had a larger influence on genotype calls in individuals with low read depth but appeared to influence genotype calls at all loci similarly. Removal of PCR clones reduced the number of called genotypes by 2% but had almost no influence on estimates of heterozygosity. As such, while steps should be taken to limit PCR clones during library preparation, PCR clones are likely not a substantial source of bias for most RAD‐seq studies.     

Is RAD‐seq suitable for phylogenetic inference? An in silico assessment and optimization

Inferring phylogenetic relationships between closely related taxa can be hindered by three factors: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD‐seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within‐species level will necessarily reduce the number of conserved and non‐duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD‐seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD‐seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD‐seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.     

Evidence for broadscale introgressive hybridization between two redfish (genus Sebastes) in the North-west Atlantic: a rare marine example

The evolutionary importance of introgressive hybridization has long been recognized by plant evolutionists, and there is now a growing recognition for its potential role in animals as well. Detailed empirical investigations of this evolutionary process, however, are still lacking in many animal groups, particularly in the marine environment. Using integrated microsatellite DNA data (eight loci analysed over 803 individuals representing 17 sampling locations) and multivariate statistical procedures (principal component, factorial correspondence and admixture proportion analyses), we: (i) provide a detailed dissection of the dynamics of introgressive hybridization between Sebastes fasciatus and S. mentella, two economically important redfishes from the North-west Atlantic; and (ii) infer the factors potentially involved in the maintenance of the hybrid zone observed in the gulf of St. Lawrence and south of Newfoundland. This study provided one of the rare examples of extensive introgressive hybridization in the ocean, and highlighted the predominant role of this process in shaping the extent of genetic diversity, interspecific differences and population structuring among redfishes from the North-west Atlantic. The extensive (average rate of introgression = 15%) but geographically circumscribed and asymmetrical pattern of introgressive hybridization, the sympatric persistence of two reproductively isolated introgressed groups, the differential patterns of linkage disequilibrium among samples, and the maintenance of genetic integrity of both species outside the defined zone of introgression despite high potential for gene flow, all implicated selection in promoting and maintaining the observed pattern of introgression.     

Double‐digest RAD sequencing outperforms microsatellite loci at assigning paternity and estimating relatedness: A proof of concept in a highly promiscuous bird

Information on genetic relationships among individuals is essential to many studies of the behaviour and ecology of wild organisms. Parentage and relatedness assays based on large numbers of single nucleotide polymorphism (SNP) loci hold substantial advantages over the microsatellite markers traditionally used for these purposes. We present a double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) analysis pipeline that, as such, simultaneously achieves the SNP discovery and genotyping steps and which is optimized to return a statistically powerful set of SNP markers (typically 150–600 after stringent filtering) from large numbers of individuals (up to 240 per run). We explore the trade‐offs inherent in this approach through a set of experiments in a species with a complex social system, the variegated fairy‐wren (Malurus lamberti) and further validate it in a phylogenetically broad set of other bird species. Through direct comparisons with a parallel data set from a robust panel of highly variable microsatellite markers, we show that this ddRAD‐seq approach results in substantially improved power to discriminate among potential relatives and considerably more precise estimates of relatedness coefficients. The pipeline is designed to be universally applicable to all bird species (and with minor modifications to many other taxa), to be cost‐ and time‐efficient, and to be replicable across independent runs such that genotype data from different study periods can be combined and analysed as field samples are accumulated.     

Conservation genetics of redside dace (<Emphasis Type="Italic">Clinostomus elongatus</Emphasis>): phylogeography and contemporary spatial structure

Redside dace Clinostomus elongatus (Teleostei: Cyprinidae) is a species of conservation concern that is declining throughout its range as a result of habitat fragmentation, degradation and loss. We characterized the genetic structure and diversity of redside dace populations across the species range using mitochondrial and microsatellite data to inform conservation efforts and assess how historical and recent events have shaped genetic structure and diversity within and among populations. Phylogeographic structure among 28 redside dace populations throughout southern Ontario (Canada) and the United States was assessed by sequence analysis of the mitochondrial cytochrome b and ATPase 6 and 8 genes. Populations were also genotyped using ten microsatellite loci to examine genetic diversity within and among populations as well as contemporary spatial structuring. Mitochondrial DNA sequence data revealed three geographically distinct lineages, which were highly concordant with groupings identified by microsatellite analysis. The combined genetic data refute published glacial refugia hypotheses of a single Mississippian refugium or of two lineages associated with Mississippian and Atlantic refugia. Secondary contact between the two eastern groups was documented in the Allegheny River drainage and tributaries to Lake Ontario. With the exception of several allopatric populations within the Allegheny watershed, high genetic structuring among populations suggests their isolation, indicating that recovery efforts should be population-based.     

Substantial differences in bias between single‐digest and double‐digest RAD‐seq libraries: A case study

The trade‐offs of using single‐digest vs. double‐digest restriction site‐associated DNA sequencing (RAD‐seq) protocols have been widely discussed. However, no direct empirical comparisons of the two methods have been conducted. Here, we sampled a single population of Gulf pipefish (Syngnathus scovelli) and genotyped 444 individuals using RAD‐seq. Sixty individuals were subjected to single‐digest RAD‐seq (sdRAD‐seq), and the remaining 384 individuals were genotyped using a double‐digest RAD‐seq (ddRAD‐seq) protocol. We analysed the resulting Illumina sequencing data and compared the two genotyping methods when reads were analysed either together or separately. Coverage statistics, observed heterozygosity, and allele frequencies differed significantly between the two protocols, as did the results of selection components analysis. We also performed an in silico digestion of the Gulf pipefish genome and modelled five major sources of bias: PCR duplicates, polymorphic restriction sites, shearing bias, asymmetric sampling (i.e., genotyping fewer individuals with sdRAD‐seq than with ddRAD‐seq) and higher major allele frequencies. This combination of approaches allowed us to determine that polymorphic restriction sites, an asymmetric sampling scheme, mean allele frequencies and to some extent PCR duplicates all contribute to different estimates of allele frequencies between samples genotyped using sdRAD‐seq versus ddRAD‐seq. Our finding that sdRAD‐seq and ddRAD‐seq can result in different allele frequencies has implications for comparisons across studies and techniques that endeavour to identify genomewide signatures of evolutionary processes in natural populations.