Performance and precision of double digestion RAD (ddRAD) genotyping in large multiplexed datasets of marine fish species
Maroso, F.; Hillen, J.E.J.; Pardo, B.G.; Gkagkavouzis, K.; Coscia, I.; Hermida, M.; Franch, R.; Hellemans, B.; Van Houdt, J.; Simionati, B.; Taggart, J.B.; Nielsen, E.E.; Maes, G.; Ciavaglia, S.A.; Webster, L.M.I.; Volckaert, F.A.M.; Martínez, P.; Bargelloni, L.; Ogden, R.; AquaTrace Consortium (2018). Performance and precision of double digestion RAD (ddRAD) genotyping in large multiplexed datasets of marine fish species. Marine Genomics 39: 64-72. https://dx.doi.org/10.1016/j.margen.2018.02.002
In: Marine Genomics. Elsevier: Amsterdam. ISSN 1874-7787; e-ISSN 1876-7478, more
| |
Keywords |
Dicentrarchus labrax (Linnaeus, 1758) [WoRMS]; Scophthalmus maximus (Linnaeus, 1758) [WoRMS]; Sparus aurata Linnaeus, 1758 [WoRMS] Marine/Coastal |
Author keywords |
ddRAD; European sea bass; GBS; Gilthead sea bream; Sequencing precision;Turbot |
Authors | | Top |
- Maroso, F.
- Hillen, J.E.J., more
- Pardo, B.G.
- Gkagkavouzis, K.
- Coscia, I., more
- Hermida, M.
- Franch, R.
|
- Hellemans, B., more
- Van Houdt, J., more
- Simionati, B.
- Taggart, J.B.
- Nielsen, E.E.
- Maes, G., more
- Ciavaglia, S.A.
|
- Webster, L.M.I.
- Volckaert, F.A.M., more
- Martínez, P.
- Bargelloni, L.
- Ogden, R.
- AquaTrace Consortium
|
Abstract |
The development of Genotyping-By-Sequencing (GBS) technologies enables cost-effective analysis of large numbers of Single Nucleotide Polymorphisms (SNPs), especially in “non-model” species. Nevertheless, as such technologies enter a mature phase, biases and errors inherent to GBS are becoming evident. Here, we evaluated the performance of double digest Restriction enzyme Associated DNA (ddRAD) sequencing in SNP genotyping studies including high number of samples. Datasets of sequence data were generated from three marine teleost species (>5500 samples, >2.5 × 1012 bases in total), using a standardized protocol. A common bioinformatics pipeline based on STACKS was established, with and without the use of a reference genome. We performed analyses throughout the production and analysis of ddRAD data in order to explore (i) the loss of information due to heterogeneous raw read number across samples; (ii) the discrepancy between expected and observed tag length and coverage; (iii) the performances of reference based vs. de novo approaches; (iv) the sources of potential genotyping errors of the library preparation/bioinformatics protocol, by comparing technical replicates. Our results showed use of a reference genome and a posteriori genotype correction improved genotyping precision. Individual read coverage was a key variable for reproducibility; variance in sequencing depth between loci in the same individual was also identified as an important factor and found to correlate to tag length. A comparison of downstream analysis carried out with ddRAD vs single SNP allele specific assay genotypes provided information about the levels of genotyping imprecision that can have a significant impact on allele frequency estimations and population assignment. The results and insights presented here will help to select and improve approaches to the analysis of large datasets based on RAD-like methodologies. |
|