Predicting RAD-seq marker numbers across the eukaryotic tree of life
Reyes-Herrera, Paula H.
Shank, Timothy M.
MetadataShow full item record
KeywordRAD-seq; Reduced representation sequencing; PredRAD; Experimental design; Genome size prediction; Restriction recognition sequence probability
High-throughput sequencing of reduced representation libraries obtained through digestion with restriction enzymes—generically known as restriction site associated DNA sequencing (RAD-seq)—is a common strategy to generate genome-wide genotypic and sequence data from eukaryotes. A critical design element of any RAD-seq study is knowledge of the approximate number of genetic markers that can be obtained for a taxon using different restriction enzymes, as this number determines the scope of a project, and ultimately defines its success. This number can only be directly determined if a reference genome sequence is available, or it can be estimated if the genome size and restriction recognition sequence probabilities are known. However, both scenarios are uncommon for nonmodel species. Here, we performed systematic in silico surveys of recognition sequences, for diverse and commonly used type II restriction enzymes across the eukaryotic tree of life. Our observations reveal that recognition sequence frequencies for a given restriction enzyme are strikingly variable among broad eukaryotic taxonomic groups, being largely determined by phylogenetic relatedness. We demonstrate that genome sizes can be predicted from cleavage frequency data obtained with restriction enzymes targeting “neutral” elements. Models based on genomic compositions are also effective tools to accurately calculate probabilities of recognition sequences across taxa, and can be applied to species for which reduced representation data are available (including transcriptomes and neutral RAD-seq data sets). The analytical pipeline developed in this study, PredRAD (https://github.com/phrh/PredRAD), and the resulting databases constitute valuable resources that will help guide the design of any study using RAD-seq or related methods.
© The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Genome Biology and Evolution 7 (2015): 3207-3225, doi:10.1093/gbe/evv210.
The following license files are associated with this item:
Except where otherwise noted, this item's license is described as Attribution-NonCommercial 4.0 International
Showing items related by title, author, creator and subject.
Herrera, Marcela; Nanninga, Gerrit; Planes, Serge; Jones, Geoffrey P.; Thorrold, Simon R.; Saenz-Agudelo, Pablo; Almany, Glenn R.; Berumen, Michael L. (The Royal Society, 2016-08-09)The persistence and resilience of many coral reef species are dependent on rates of connectivity among sub-populations. However, despite increasing research efforts, the spatial scale of larval dispersal remains unpredictable ...
Amino acid sequence of the ligand-binding domain of the aryl hydrocarbon receptor 1 predicts sensitivity of wild birds to effects of dioxin-like compounds Farmahin, Reza; Manning, Gillian E.; Crump, Doug; Wu, Dongmei; Mundy, Lukas J.; Jones, Stephanie P.; Hahn, Mark E.; Karchner, Sibel I.; Giesy, John P.; Bursian, Steven J.; Zwiernik, Matthew J.; Fredricks, Timothy B.; Kennedy, Sean W. (2012-07-17)The sensitivity of avian species to the toxic effects of dioxin-like compounds (DLCs) varies up to 1000-fold among species and this variability has been associated with inter-species differences in aryl hydrocarbon receptor ...
Moberg, Emily A.; Shyu, Esther; Herrera, Guillermo E.; Lenhart, Suzanne; Lou, Yuan; Neubert, Michael G. (2015-08)Marine reserves are an increasingly used and potentially contentious tool in fisheries management. Depending upon the way that individuals move, no-take marine reserves can be necessary for maximizing equilibrium rent in ...