These aren’t the loci you’e looking for: Principles of effective SNP filtering for molecular ecologists

Abstract

Sequencing reduced-representation libraries of restriction site-associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.

A spatial genomic approach identifies time lags and historical barriers to gene flow in a rapidly fragmenting Appalachian landscape

“…The effect of missing data in landscape genomic studies has yet to be thoroughly explored, and we suggest based on our results that some spatially informed analyses may be robust to the recovery of patterns of weak population structure despite the inclusion of a high level of missing genotypes, but that parameter estimation at this geographical scale (e.g., migration rates) may be more strongly influenced. Therefore, when possible, we second the recommendation of others (O'Leary, Puritz, Willis, Hollenbeck, & Portnoy, ; Wagner et al, ) for researchers to explore the sensitivity of their results across a range of different missing data thresholds.…”

Age estimation in a long‐lived seabird (Ardenna tenuirostris) using DNA methylation‐based biomarkers

“…Reduced representation bisulphite sequencing can also be used to quantify CpG DNAm, but does require a higher quantity of initial genomic DNA (Meissner et al, ). The output of these various techniques depends upon several molecular, platform and bioinformatic factors and choices, which is discussed in detail elsewhere (O'Leary, Puritz, Willis, Hollenbeck, & Portnoy, ). Our results now show that the DREAM method can also be used to quantify global DNAm and screen for aDMPs in non‐model animals.…”

A versatile Rapture (RAD‐Capture) platform for genotyping marine turtles

“…This makes continued collaboration essential to advancing our understanding of marine turtles, particularly for researchers with access to such resources to continue efforts to increase capacity elsewhere, such as through visiting scientist training partnerships and creation of shared genetic databases. Particularly given the influence that bioinformatics parameters (e.g., filtering criteria, assembly methodology, genotyping thresholds) can have on results (O'Leary, Puritz, Willis, Hollenbeck, & Portnoy, ), it is imperative for researchers to include metadata and analysis details to ensure robust and comparable data across laboratories and over time.…”