ddradseqtools: a software package for in silico simulation and testing of double-digest RADseq experiments.
Double-digested RADseq (ddRADseq) is a NGS methodology that generates reads from thousands of loci targeted by restriction enzyme cut sites, across multiple individuals. To be statistically sound and economically optimal, a ddRADseq experiment has a preliminary design stage that needs to consider issues related to the selection of enzymes, particular features of the genome of the focal species, possible modifications to the library construction protocol, coverage needed to minimize missing data, and the potential sources of error that may impact upon the coverage. We present ddradseqtools, a software package to help ddRADseq experimental design by (i) the generation of in silico double-digested fragments; (ii) the construction of modified ddRADseq libraries using adapters with either one or two indexes and degenerate base regions (DBRs) to quantify PCR duplicates; and (iii) the initial steps of the bioinformatics preprocessing of reads. ddradseqtools generates single-end (SE) or paired-end (PE) reads that may bear SNPs and/or indels. The effect of allele dropout and PCR duplicates on coverage is also simulated. The resulting output files can be submitted to pipelines of alignment and variant calling, to allow the fine-tuning of parameters. The software was validated with specific tests for the correct operability of the program. The correspondence between in silico settings and parameters from ddRADseq in vitro experiments was assessed to provide guidelines for the reliable performance of the software. ddradseqtools is cost-efficient in terms of execution time, and can be run on computers with standard CPU and RAM configuration.
Forest Genetics and Physiology Research Group, Technical University of Madrid (UPM), Ciudad Universitaria s/n, Madrid, Spain.
F
F
Mora-Márquez
Island Ecology and Evolution Research Group, IPNA-CSIC, Tenerife, Canary Islands, Spain.
V
V
García-Olivares
Island Ecology and Evolution Research Group, IPNA-CSIC, Tenerife, Canary Islands, Spain.|School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ, UK.
B C
BC
Emerson
Forest Genetics and Physiology Research Group, Technical University of Madrid (UPM), Ciudad Universitaria s/n, Madrid, Spain.
U
U
López de Heredia
D019295:Computational Biology / Q000379:methods*; D004247:DNA / Q000737:chemistry / Q000378:metabolism; D004262:DNA Restriction Enzymes / Q000378:metabolism; D059014:High-Throughput Nucleotide Sequencing / Q000379:methods*; D012984:Software
D004247:DNA; D004262:DNA Restriction Enzymes
PCR duplicates; allele dropout; coverage; double-digested RADseq; in silico simulation