PyRAD: assembly of de novo RADseq loci for phylogenetic analyses

By Deren A. R. Eaton. bioRxiv, 2013

732 citations recorded by the source corpus.

Restriction-site associated genomic markers are a powerful tool for investigating evolutionary questions at the population level, but are limited in their utility at deeper phylogenetic scales where fewer orthologous loci are typically recovered across disparate taxa. While this limitation stems in part from mutations to restriction recognition sites that disrupt data generation, an alternative source of data loss comes from the failure to identify homology during bioinformatic analyses. Clustering methods that allow for lower similarity thresholds and the inclusion of indel variation will perform better at assembling RADseq loci at the phylogenetic scale. PyRAD is a pipeline to assemble de novo RADseq loci with the aim of optimizing coverage across phylogenetic data sets. It utilizes a wrapper around an alignment-clustering algorithm which allows for indel variation within and between samples, as well as for incomplete overlap among reads (e.g., paired-end). Here I compare PyRAD with the program Stacks in their performance analyzing a simulated RADseq data set that includes indel variation. Indels disrupt clustering of homologous loci in Stacks but not in PyRAD, such that the latter recovers more shared loci across disparate taxa. I show through re-analysis of an empirical RADseq data set that indels are a common feature of such data, even at shallow phylogenetic scales. PyRAD utilizes parallel processing as well as an optional hierarchical clustering method which allow it to rapidly assemble phylogenetic data sets with hundreds of sampled individuals. Availability Software is written in Python and freely available at http://www.dereneaton.com/software/ Supplement Scripts to completely reproduce all simulated and empirical analyses are available in the Supplementary Materials.

Subfields: Molecular Biology, Plant Science

View Paper: https://www.semanticscholar.org/paper/a0a2de722f134347dd97ffbf8c755ec5e997ee44

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

In Molecular Biology

In Plant Science

Citations in the corpus

ipyrad: Interactive assembly and analysis of RADseq datasets (2020, methodology): “ipyrad provides a command line interface that is easy to use and which inherits interaction logic from its predecessor (Eaton, 2014).”
Restriction site‐associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference (2015, methodology): “…assembly and genotyping tools for RADseq data have recently been released, such as RaPiD (Willing et al. 2011), RADtools (Baxter et al. 2011), graph-based distance clustering approaches (Peterson et al. 2012), Stacks (Catchen et al. 2011, 2013), Rain-bow (Chong et al. 2012) and pyRAD (Eaton 2014).”