Deriving genotypes from RAD-seq short-read data using Stacks

Nicolas C. Rochette, Julian M. Catchen

Research output: Contribution to journalArticlepeer-review


R estriction site-associated DNANA sequencing (RARAD-seq) allows for the genome-wide discovery and genotyping of single-nucleotide polymorphisms in hundreds of individuals at a time in model and nonmodel species alike. However, converting short-read sequencing data into reliable genotype data remains a nontrivial task, especially as RARAD-seq is used in systems that have very diverse genomic properties. Here, we present a protocol to analyze RARAD-seq data using the Stacks pipeline. This protocol will be of use in areas such as ecology and population genetics. It covers the assessment and demultiplexing of the sequencing data, read mapping, inference of RARAD loci, genotype calling, and filtering of the output data, as well as providing two simple examples of downstream biological analyses. We place special emphasis on checking the soundness of the procedure and choosing the main parameters, given the properties of the data. The procedure can be completed in 1 week, but determining definitive methodological choices will typically take up to 1 month.

Original languageEnglish (US)
Pages (from-to)2640-2659
Number of pages20
JournalNature Protocols
Issue number12
StatePublished - Dec 1 2017

ASJC Scopus subject areas

  • General Biochemistry, Genetics and Molecular Biology


Dive into the research topics of 'Deriving genotypes from RAD-seq short-read data using Stacks'. Together they form a unique fingerprint.

Cite this