polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids

This R package is currently under development. See the list of future features to see where it is headed.

I'm always interested in new collaboration! If you find polyRAD to be helpful in your research, let me know if you'd be interested in sharing your data and results for coauthorship in the publication describing the software. I would also like to hear your feature requests. Contact: Lindsay Clark, University of Illinois, Urbana-Champaign.


Genotypes derived from genotyping-by-sequencing (GBS) and restriction site-associated DNA sequencing (RAD-seq) have inherent uncertainty associated with them due to sampling error, i.e. some alleles might not get sequenced at all, or might not be sequenced in exact proportion to their copy number in the genome. This package imports read depth in a variety of formats output by various bioinformatics pipelines and estimates the probability of each possible genotype for each taxon and locus. Unlike similar pipelines, polyRAD can account for population structure and variable inheritance modes (autopolyploid, allopolyploid, intermediate). Gentoypes and/or probability distributions can then be exported for downstream analysis such as genome-wide association, genomic selection, QTL mapping, or population structure analysis.

Formats supported

To hopefully answer the question, "Can I use polyRAD?":

polyRAD requires as input the sequence read depth at each allele for each sample. Alleles must also be grouped into loci. The bioinformatics pipeline that you used for SNP discovery did not have to assume polyploidy, as long as it faithfully reported allelic read depth. Genomic alignment information is optional. Right now there are data import functions for the following formats:

Currently there are export functions for the following software. Genotypes are exported as continuous rather than discrete variables.


polyRAD is not yet on CRAN. To install, go to the releases page and download the binary for the latest release. This is polyRAD_0.x.tar.gz if you are on Mac or Linux, or if you are on Windows. Don't unzip the file. Move it into your working directory in R, then run

install.packages("", repos = NULL)

replacing the file name with the one you just downloaded.


The tutorial document for the package is available on Github.


DOI for version 0.3: DOI

Version 0.1 was also presented in a poster:

Clark LV, Lipka AE, and Sacks EJ (2018) polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids. Plant and Animal Genome Conference XXVI, January 13-17, San Diego, California, USA. DOI:10.13140/RG.2.2.27134.08001


This material is based upon work supported by the National Science Foundation under Grant No. 1661490. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.

lvclark/polyRAD documentation built on May 25, 2018, 1:28 p.m.