Bioinformatic processing of RAD-seq data dramatically impacts downstream population genetic inference

Methods in Ecology and Evolution

Volumn 8, Issue 8, 2017, Pages 907-917

  Shafer, Aaron B A ^a,b   Peart, Claire R ^a   Tusso, Sergio ^a   Maayan, Inbar ^a   Brelsford, Alan ^c   Wheat, Christopher W ^d   Wolf, Jochen B W ^a,e  

^a UPPSALA UNIVERSITY   (Sweden)

^b TRENT UNIVERSITY   (Canada)

^c UNIVERSITY OF LAUSANNE   (Switzerland)

^d STOCKHOLM UNIVERSITY   (Sweden)

^e UNIVERSITY OF MUNICH   (Germany)

Author keywords

bioinformatics; GBS; genotyping by sequencing; population genetics; RAD seq

Indexed keywords

EID: 85006240194     PISSN: None     EISSN: 2041210X     Source Type: Journal    
DOI: 10.1111/2041-210X.12700     Document Type: Article

References (68)

1. Andrews, K.R., Hohenlohe, P.A., Miller, M.R., Hand, B.K., Seeb, J.E. & Luikart, G. (2014). Trade-offs and utility of alternative RADseq methods: reply to Puritz et al. Molecular Ecology, 23, 5943–5946.
2. Andrews, K.R., Good, J.M., Miller, M.R., Luikart, G. & Hohenlohe, P.A. (2016) Harnessing the power of RADseq for ecological and evolutionary genomics. Nature Reviews Genetics, 17, 81–92.
3. Arnold, B., Corbett-Detig, R.B., Hartl, D. & Bomblies, K. (2013) RADseq underestimates diversity and introduces genealogical biases due to nonrandom haplotype sampling. Molecular Ecology, 22, 3179–3190.
4. Brelsford, A., Dufresnes, C. & Perrin, N. (2016) High-density sex-specific linkage maps of a European tree frog (Hyla arborea) identify the sex chromosome without information on offspring sex. Heredity, 116, 177–181.
5. Callicrate, T., Dikow, R., Thomas, J.W., Mullikin, J.C., Jarvis, E.D. & Fleischer, R.C. (2014) Genomic resources for the endangered Hawaiian honeycreepers. BMC Genomics, 15, 1098.
6. Catchen, J., Hohenlohe, P.A., Bassham, S., Amores, A. & Cresko, W.A. (2013) Stacks: an analysis tool set for population genomics. Molecular Ecology, 22, 3124–3140.
7. Chaisson, M.J.P., Wilson, R.K. & Eichler, E.E. (2015) Genetic variation and the de novo assembly of human genomes. Nature Reviews Genetics, 16, 627–640.
8. Chong, Z., Ruan, J. & Wu, C.-I. (2012) Rainbow: an integrated tool for efficient clustering and assembling RAD-seq reads. Bioinformatics, 28, 2732–2737.
9. DaCosta, J.M. & Sorenson, M.D. (2014) Amplification biases and consistent recovery of loci in a double-digest RAD-seq protocol. PLoS ONE, 9, e106713.
10. Danecek, P., Auton, A., Abecasis, G. et al. & 1000 Genomes Project Analysis Group. (2011). The variant call format and VCFtools. Bioinformatics (Oxford, England), 27, 2156–2158.
11. Davey, J.W., Cezard, T., Fuentes-Utrilla, P., Eland, C., Gharbi, K. & Blaxter, M.L. (2012) Special features of RAD sequencing data: implications for genotyping. Molecular Ecology, 22, 3151–3164.
12. DePristo, M.A., Banks, E., Poplin, R. & Garimella, K.V. (2011) A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nature Genetics, 43, 491–498.
13. Doležel, J., Bartoš, J., Voglmayr, H. & Greilhuber, J. (2003) Letter to the editor. Cytometry Part A, 51A, 127–128.
14. Du, B. & Wang, D. (2006) C-values of seven marine mammal species determined by flow cytometry. Zoological Science, 23, 1017–1020.
15. Eaton, D.A.R. (2014). PyRAD: assembly of de novo RADseq loci for phylogenetic analyses. Bioinformatics (Oxford, England), 30, 1844–1849.
16. Ekblom, R. & Wolf, J.B.W. (2014) A field guide to whole-genome sequencing, assembly and annotation. Evolutionary Applications, 7, 1026–1042.
17. Ellegren, H. (2014) Genome sequencing and population genomics in non-model organisms. Trends in Ecology & Evolution, 29, 51–63.
18. Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. (2012) CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics, 28, 3150–3152.
19. Garrison, E. & Marth, G. (2012). Haplotype-based variant detection from short-read sequencing. arXiv:1207.3907 [q-bio].
20. Gautier, M., Gharbi, K., Cezard, T., Foucaud, J., Kerdelhué, C., Pudlo, P., Cornuet, J.-M. & Estoup, A. (2013) The effect of RAD allele dropout on the estimation of genetic variation within and between populations. Molecular Ecology, 22, 3165–3178.
21. Goslee, S.C. & Urban, D.L. (2007) The ecodist package for dissimilarity-based Analysis of Ecological Data | Goslee | Journal of Statistical Software. Journal of Statistical Software, 22, 19.
22. Guo, Y., Li, J., Li, C.-I., Long, J., Samuels, D.C. & Shyr, Y. (2012) The effect of strand bias in Illumina short-read sequencing data. BMC Genomics, 13, 666.
23. Gutenkunst, R.N., Hernandez, R.D., Williamson, S.H. & Bustamante, C.D. (2009) Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data (ed. G. McVean). PLoS Genetics, 5, e1000695–11.
24. Henning, F., Lee, H.J., Franchini, P. & Meyer, A. (2014) Genetic mapping of horizontal stripes in Lake Victoria cichlid fishes: benefits and pitfalls of using RAD markers for dense linkage mapping. Molecular Ecology, 23, 5224–5240.
25. Hoffman, J.I., Simpson, F., David, P., Rijks, J.M., Kuiken, T., Thorne, M.A.S., Lacy, R.C. & Dasmahapatra, K.K. (2014) High-throughput sequencing reveals inbreeding depression in a natural population. Proceedings of the National Academy of Sciences of the United States of America, 111, 3775–3780.
26. Hohenlohe, P.A., Bassham, S., Etter, P.D., Stiffler, N., Johnson, E.A. & Cresko, W.A. (2010). Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags (P.A. Hohenlohe, S. Bassham, P.D. Etter, N. Stiffler, E.A. Johnson & W.A. Cresko, Eds.). PLoS Genetics, 6, e1000862.
27. Humble, E., Martinez-Barrio, A., Forcada, J., Trathan, P.N., Thorne, M.a.S., Hoffmann, M., Wolf, J.B.W. & Hoffman, J.I. (2016) A draft fur seal genome provides insights into factors affecting SNP validation and how to mitigate them. Molecular Ecology Resources, 16, 909–921.
28. Jeglinski, J.W.E., Wolf, J.B.W., Werner, C., Costa, D.P. & Trillmich, F. (2015) Differences in foraging ecology align with genetically divergent ecotypes of a highly mobile marine top predator. Oecologia, 179, 1041–1052.
29. Johnson, P.L.F. & Slatkin, M. (2008) Accounting for bias from sequencing error in population genetic estimates. Molecular Biology and Evolution, 25, 199–206.
30. Kjeldsen, S.R., Zenger, K.R., Leigh, K., Ellis, W., Tobey, J., Phalen, D., Melzer, A., FitzGibbon, S. & Raadsma, H.W. (2015) Genome-wide SNP loci reveal novel insights into koala (Phascolarctos cinereus) population variability across its range. Conservation Genetics, 17, 337–353.
31. Korneliussen, T.S., Albrechtsen, A. & Nielsen, R. (2014) ANGSD: analysis of Next Generation Sequencing Data. BMC Bioinformatics, 15, 1.
32. Li, H. (2011) A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics, 27, 2987–2993.
33. Li, H. (2014) Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics, 30, 2843–2851.
34. Li, H. & Durbin, R. (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 25, 1754–1760.
35. Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G. & Durbin, R. & 1000 Genome Project Data Processing Subgroup (2009) The sequence alignment/map format and SAMtools. Bioinformatics, 25, 2078–2079.
36. Li, R., Fan, W., Tian, G. et al. (2010) The sequence and de novo assembly of the giant panda genome. Nature, 463, 311–317.
37. Lunter, G. & Goodson, M. (2011) Stampy: a statistical algorithm for sensitive and fast mapping of Illumina sequence reads. Genome Research, 21, 936–939.
38. Lynch, M. (2009) Estimation of allele frequencies from high-coverage genome-sequencing projects. Genetics, 182, 295–301.
39. Mandeville, E.G., Parchman, T.L., McDonald, D.B. & Buerkle, C.A. (2015) Highly variable reproductive isolation among pairs of Catostomusspecies. Molecular Ecology, 24, 1856–1872.
40. Mantel, N. (1967) The detection of disease clustering and a generalized regression approach. Cancer Research, 27, 209–220.
41. Martin, M. (2011) Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal, 17, 10–12.
42. Martin, S.H., Dasmahapatra, K.K., Nadeau, N.J. et al. (2013) Genome-wide evidence for speciation with gene flow in Heliconius butterflies. Genome Research, gr.159426.113.
43. Mastretta Yanes, A., Arrigo, N., Alvarez, N., Jorgensen, T.H., Piñero, D. & Emerson, B.C. (2015) Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference. Molecular Ecology Resources, 15, 28–41.
44. McKenna, A., Hanna, M., Banks, E. et al. (2010) The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research, 20, 1297–1303.
45. Nielsen, R., Korneliussen, T., Albrechtsen, A., Li, Y. & Wang, J. (2012) SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data (ed. P. Awadalla). PLoS ONE, 7, e37558.
46. Nyman, T., Valtonen, M., Aspi, J., Ruokonen, M., Kunnasranta, M. & Palo, J.U. (2014) Demographic histories and genetic diversities of Fennoscandian marine and landlocked ringed seal subspecies. Ecology and Evolution, 4, 3420–3434.
47. Pante, E., Abdelkrim, J., Viricel, A., Gey, D., France, S.C., Boisselier, M.C. & Samadi, S. (2015) Use of RAD sequencing for delimiting species. Heredity, 114, 450–459.
48. Parchman, T.L., Gompert, Z., Mudge, J., Schilkey, F.D., Benkman, C.W. & Buerkle, C.A. (2012) Genome-wide association genetics of an adaptive trait in lodgepole pine. Molecular Ecology, 21, 2991–3005.
49. Peterson, B.K., Weber, J.N., Kay, E.H., Fisher, H.S. & Hoekstra, H.E. (2012) Double digest RADseq: an inexpensive method for De Novo SNP discovery and genotyping in model and non-model species (ed. L. Orlando). PLoS ONE, 7, e37135.
50. Pool, J.E., Hellmann, I., Jensen, J.D. & Nielsen, R. (2010) Population genetic inference from genomic sequence variation. Genome Research, 20, 291–300.
51. Pörschmann, U., Trillmich, F., Mueller, B. & Wolf, J.B.W. (2010) Male reproductive success and its behavioural correlates in a polygynous mammal, the Galápagos sea lion (Zalophus wollebaeki). Molecular Ecology, 19, 2574–2586.
52. Puritz, J.B., Hollenbeck, C.M. & Gold, J.R. (2014) dDocent?: a RADseq, variant-calling pipeline designed for population genomics of non-model organisms. PeerJ, 2, e431.
53. Puritz, J.B., Matz, M.V., Toonen, R.J., Weber, J.N., Bolnick, D.I. & Bird, C.E. (2014) Demystifying the RAD fad. Molecular Ecology, 23, 5937–5942.
54. Rousset, F. (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics, 145, 1219–1228.
55. Rousset, F. (2008) genepop'007: a complete re-implementation of the genepop software for Windows and Linux. Molecular Ecology Resources, 8, 103–106.
56. Shafer, A.B., Davis, C.S., Coltman, D.W. & Stewart, R.E. (2014) Microsatellite assessment of walrus (Odobenus rosmarus rosmarus) stocks in Canada. NAMMCO Scientific Publications, 9, 15–31.
57. Shafer, A.B.A., Gattepaille, L.M., Stewart, R.E.A. & Wolf, J.B.W. (2015a) Demographic inferences using short-read genomic data in an approximate Bayesian computation framework: in silico evaluation of power, biases and proof of concept in Atlantic walrus. Molecular Ecology, 24, 328–345.
58. Shafer, A.B.A., Wolf, J.B.W., Alves, P.C. et al. (2015b). Genomics and the challenging translation into conservation practice. Trends in Ecology & Evolution, 30, 78–87.
59. Shafer, A.B.A., Northrup, J.M., Wikelski, M. & Wittemyer, G. (2016a) Forecasting ecological genomics: high-tech animal instrumentation meets high-throughput sequencing. PLoS Biology, 14, e1002350.
60. Shafer, A.B.A., Peart, C.R., Tusso, S., Maayan, I., Brelsford, A., Wheat, C.W. & Wolf, J.B.W. (2016b) Data from: Bioinformatic processing of RAD-seq data dramatically impacts downstream population genetic inference. Dryad Digital Repository, http://dx.doi.org/10.5061/dryad.q14c1
61. Sims, D., Sudbery, I., Ilott, N.E., Heger, A. & Ponting, C.P. (2014) Sequencing depth and coverage: key considerations in genomic analyses. Nature Reviews Genetics, 15, 121–132.
62. Smouse, P.E., Long, J.C. & Sokal, R.R. (1986) Multiple regression and correlation extensions of the mantel test of matrix correspondence. Systematic Biology, 35, 627–632.
63. Sovic, M.G., Fries, A.C. & Gibbs, H.L. (2015) AftrRAD: a pipeline for accurate and efficient de novo assembly of RADseq data. Molecular Ecology Resources, 15, 1163–1171.
64. Torkamaneh, D., Laroche, J. & Belzile, F. (2016) Genome-wide SNP calling from genotyping by sequencing (GBS) data: a comparison of seven pipelines and two sequencing technologies. PLoS ONE, 11, e0161333.
65. Vijay, N., Poelstra, J.W., Künstner, A. & Wolf, J.B.W. (2012) Challenges and strategies in transcriptome assembly and differential gene expression quantification. A comprehensive in silicoassessment of RNA-seq experiments. Molecular Ecology, 22, 620–634.
66. Weir, B.S. & Cockerham, C.C. (1984) Estimating F-statistics for the analysis of population structure. Evolution, 38, 1358–1370.
67. Wolf, J.B.W., Harrod, C., Brunner, S., Salazar, S., Trillmich, F. & Tautz, D. (2008) Tracing early stages of species differentiation: ecological, morphological and genetic divergence of Galápagos sea lion populations. BMC Evolutionary Biology, 8, 150.
68. Zhang, J., Kobert, K., Flouri, T. & Stamatakis, A. (2014) PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics, 30, 614–620.