title: Comparing three different methods to detect selective loci using dominant markers.
abstract: We carried out a simulation study to compare the efficiency of three alternative programs (DFDIST, DETSELD and BAYESCAN) to detect loci under directional selection from genome-wide scans using dominant markers. We also evaluated the efficiency of correcting for multiple testing those methods that use a classical probability approach. Under a wide range of scenarios, we conclude that BAYESCAN appears to be more efficient than the other methods, detecting a usually high percentage of true selective loci as well as less than 1% of outliers (false positives) under a fully neutral model. In addition, the percentage of outliers detected by this software is always correlated with the true percentage of selective loci in the genome. Our results show, nevertheless, that false positives are common even with a combination of methods and multitest correction, suggesting that conclusions obtained from this approach should be taken with extreme caution.
authors.affiliation: Departamento de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain.
authors.forename: A
authors.initials: A
authors.lastname: Pérez-Figueroa
authors.affiliation: Departamento de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain.
authors.forename: M J
authors.initials: MJ
authors.lastname: García-Pereira
authors.affiliation: Departamento de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain.
authors.forename: M
authors.initials: M
authors.lastname: Saura
authors.affiliation: Departamento de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain.
authors.forename: E
authors.initials: E
authors.lastname: Rolán-Alvarez
authors.affiliation: Departamento de Bioquímica, Genética e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain.
authors.forename: A
authors.initials: A
authors.lastname: Caballero
mesh_terms: D003198:Computer Simulation; D005799:Genes, Dominant; D056426:Genetic Loci; D016678:Genome; D023281:Genomics / Q000379:methods; D012641:Selection, Genetic; D012984:Software
references.citation: Acheré, V., Favre, J.M., Besnard, G. & Jeandroz, S. 2005. Genomic organization of molecular differentiation in Norway spruce (Picea abies). Mol. Ecol. 14: 3191-3201.
references.citation: Beaumont, M.A. & Balding, D.J. 2004. Identifying adaptive genetic divergence among populations from genome scans. Mol. Ecol. 13: 969-980.
references.citation: Beaumont, M.A. & Nichols, R.A. 1996. Evaluating loci for use in the genetic analysis of population structure. Proc. Biol. Sci. 263: 1619-1626.
references.citation: Benjamini, Y. & Hochberg, Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57: 289-300.
references.citation: Bensch, S. & Akesson, M. 2005. Ten years of AFLP in ecology and evolution: why so few animals? Mol. Ecol. 14: 2899-2914.
references.citation: Bitocchi, E., Nanni, L., Rossi, M., Bellucci, E., Giardini, E., Buonamici, A., Vendramin, G.G. & Papa, R. 2009. Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy. Mol. Ecol. 18: 603-621.
references.citation: Blel, H., Panfili, J., Guinand, B., Berrebi, P., Said, K. & Durand, J.-D. 2010. Selection footprint at the first intron of the Prl gene in natural populations of the flathead mullet (Mugil cephalus, L. 1758). J. Exp. Mar. Biol. Ecol. 387: 60-67.
references.citation: Bonhomme, M., Blancher, A., Jalil, M.F. & Crouau-Roy, B. 2007. Factors shaping genetic variation in the MHC of natural non-human primate populations. Tissue Antigens 70: 398-411.
references.citation: Bonin, A., Taberlet, P., Miaud, C. & Pompanon, F. 2006. Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria). Mol. Biol. Evol. 23: 773-783.
references.citation: Bonin, A., Nicole, F., Pompanon, F., Miaud, C. & Taberlet, P. 2007. Population adaptive index: a new method to help measure intraspecific genetic diversity and prioritize populations for conservation. Conserv. Biol. 21: 697-708.
references.citation: Bryja, J., Charbonnel, N., Berthier, K., Galan, M. & Cosson, J.-F. 2007. Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles. Mol. Ecol. 16: 5084-5097.
references.citation: Caballero, A., Quesada, H. & Rolán-Álvarez, E. 2008. Impact of amplified fragment length polymorphism size homoplasy on the estimation of population genetic diversity and the detection of selective loci. Genetics 179: 539-554.
references.citation: Carvajal-Rodriguez, A., De Uña-Alvarez, J. & Rolan-Alvarez, E. 2009. A new multitest correction (SGoF) that increases its statistical power when increasing the number of tests. BMC Bioinformatics 10: 209. See http://webs.uvigo.es/acraaj/SGoF.htm
references.citation: Chen, L. & Yang, G. 2009. A genomic scanning using AFLP to detect candidate loci under selection in the finless porpoise (Neophocaena phocaenoides). Genes Genet. Syst. 84: 307-313.
references.citation: Egan, S.P., Nosil, P. & Funk, D. 2008. Selection and genomic differentiation during ecological speciation: isolating the contributions of host association via comparative genome scan of Neochlamisus bebbianae leaf beetles. Evolution 62: 1162-1181.
references.citation: Eveno, E., Collada, C., Guevara, M.A., Leger, V., Soto, A., Diaz, L., Leger, P., Gonzalez-Martinez, S.C., Cervera, M.T., Plomion, C. & Garnier-Gere, P.H. 2008. Contrasting patterns of selection at Pinus pinaster Ait. Drought stress candidate genes as revealed by genetic differentiation analyses. Mol. Biol. Evol. 25: 417-437.
references.citation: Excoffier, L., Hofer, T. & Foll, M. 2009. Detecting loci under selection in a hierarchically structured population. Heredity 103: 285-298.
references.citation: Flint, J., Bond, J., Rees, D.C., Boyce, A.J., Roberts-Thomson, J.M., Excoffier, L., Clegg, M.A., Beaumont, M.A., Nichols, R.A. & Harding, R.M. 1999. Minisatellite mutational processes reduce Fst estimates. Hum. Genet. 105: 567-576.
references.citation: Foll, M. & Gaggiotti, O. 2008. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics 180: 977-995.
references.citation: Freeland, J.R., Biss, P., Conrad, K.F. & Silvertown, J. 2010. Selection pressures have caused genome-wide population differentiation of Anthoxanthum odoratum despite the potential for high gene flow. J. Evol. Biol. 23: 776.
references.citation: Gaggiotti, O., Bekkevold, D., Jorgensen, H.B.H., Foll, M., Carvalho, G.R., Andre, C. & Ruzzante, D.E. 2009. Disentangling the effects of evolutionary, demographic, and environmental factors influencing genetic structure of natural populations: Atlantic herring as a case study. Evolution 63: 2939-2951.
references.citation: Gagnaire, P.A., Albert, V., Jónsson, B. & Bernatchez, L. 2009. Natural selection influences AFLP intraspecific genetic variability and introgression patterns in Atlantic eels. Mol. Ecol. 18: 1678-1691.
references.citation: Galindo, J., Morán, P. & Rolán-Álvarez, E. 2009. Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis. Mol. Ecol. 18: 919-930.
references.citation: Holderegger, R., Herrmann, D., Poncet, B., Gugerli, F., Thuiller, W., Taberlet, P., Gielly, L., Rioux, D., Brodbeck, S., Aubert, S. & Manel, S. 2008. Land ahead: using genome scans to identify molecular markers of adaptive relevance. Plant Ecol. Divers, 1: 273-283.
references.citation: Joost, S., Bonin, A., Bruford, M.W., Després, S.L., Conord, C., Erhardt, G. & Taberlet, P. 2007. A spatial analysis method (SAM) to detect candidate loci for selection: towards a landscape genomics approach to adaptation. Mol. Ecol. 16: 3955-3969.
references.citation: Joost, S., Kalbermatten, M. & Bonin, A. 2008. Spatial analysis method (SAM): a software tool combining molecular and environmental data to identify candidate loci for selection. Mol. Ecol. Resour. 8: 957-960.
references.citation: Jump, A.S., Hunt, J.M., Martínez-Izquierdo, J.A. & Peñuelas, J. 2006. Natural selection and climate change: temperature-linked spatial and temporal trends in gene frequency in Fagus sylvatica. Mol. Ecol. 15: 3469-3480.
references.citation: Kane, N.C. & Rieseberg, L. 2007. Selective sweeps reveal candidate genes for adaptation to drought and salt tolerance in common sunflower, Helianthus annuus. Genetics 175: 1823-1834.
references.citation: Knapen, D., De Wolf, H., Knaepkens, G., Bervoets, L., Eens, M., Blust, R. & Verheyen, E. 2009. Historical metal pollution in natural gudgeon populations: inferences from allozyme, microsatellite and condition factor analysis. Aquat. Toxicol. 95: 17-26.
references.citation: Lewontin, R. & Krakauer, J. 1973. Distribution of gene frequency as a test of theory of selective neutrality of polymorphisms. Genetics 74: 175-195.
references.citation: Manel, S., Conord, C. & Després, L. 2009. Genome scan to assess the respective role of host-plant and environmental constraints on the adaptation of a widespread insect. BMC Evol. Biol. 9: 288-297.
references.citation: Marden, J.I. 2000. Hypothesis testing: from p values to Bayes factors. J. Am. Stat. Assoc. 95: 1316-1320.
references.citation: Mealor, B.A. & Hild, A. 2006. Potential selection in native grass populations by exotic invasion. Mol. Ecol. 15: 2291-2300.
references.citation: Medugorac, I., Medugorac, A., Russ, I., Veit-Kensch, C.E., Taberlet, P., Luntz, B., Mix, H.M. & Förster, M. 2009. Genetic diversity of European cattle breeds highlights the conservation value of traditional unselected breeds with high effective population size. Mol. Ecol. 18: 3394-3410.
references.citation: Meudt, H.M. & Clarke, A.C. 2007. Almost forgotten or latest practice? AFLP applications, analyses and advances Trends Plant Sci. 12: 106-117.
references.citation: Meyer, C.-L., Vitalis, R., Saumitou-Laprade, P. & Castric, C. 2009. Genomic pattern of adaptive divergence in Arabidopsis halleri, a model species for tolerance to heavy metal. Mol. Ecol. 18: 2050-2062.
references.citation: Michalski, S., Durka, W., Jentsch, A., Kreyling, J., Pompe, S., Schweiger, O., Willner, E. & Beierkuhnlein, C. 2010. Evidence for genetic differentiation and divergent selection in an autotetraploid forage grass (Arrhenatherum elatius). Theor. Appl. Genet. 120: 1151-1162.
references.citation: Miller, N.J., Ciosi, M., Sappington, T.W., Ratcliffe, S.T., Spencer, J.L. & Giuillemaud, T. 2007. Genome scan of Diabrotica virgifera virgifera for genetic variation associated with crop rotation tolerance. J. Appl. Entomol. 131: 378-385.
references.citation: Mueller, U.G. & Wolfenbarger, L. 1999. AFLP genotyping and fingerprinting. Trends Ecol. Evol. 14: 389-394.
references.citation: Murray, M.C. & Hare, M.P. 2006. A genomic scan for divergent selection in a secondary contact zone between Atlantic and Gulf of Mexico oysters, Crassostrea virginica. Mol. Ecol. 15: 4229-4242.
references.citation: Nei, M. & Chakravarti, A. 1977. Drift variances of Fst and Gst statistics obtained from a finite number of isolated populations. Theor. Popul. Biol. 11: 307-325.
references.citation: Nielsen, R. 2005. Molecular signatures of natural selection. Annu. Rev. Genet. 39: 197-218.
references.citation: Nielsen, E.E., Hemmer-Hansen, J., Poulsen, N.A., Loeschcke, V., Moen, T., Johansen, T., Mittelholzer, C., Taranger, G.-L., Ogden, R. & Carvalho, G.R. 2009. Genomic signatures of local directional selection in a high gene flow marine organism: the Atlantic cod (Gadus morhua). BMC Evol. Biol. 9: 276-286.
references.citation: Nosil, P., Egan, S.P. & Funk, D. 2008. Heterogeneous genomic differentiation between walking-stick ecotypes: “isolation by adaptation” and multiple roles for divergent selection. Evolution 62: 316-336.
references.citation: O’Malley, K.G., Camara, M.D. & Banks, M. 2007. Candidate loci reveal genetic differentiation between temporally divergent migratory runs of Chinook salmon (Oncorhynchus tshawytscha). Mol. Ecol. 16: 4930-4941.
references.citation: Oetjen, K. & Reusch, T.B.H. 2007. Genome scans detect consistent divergent selection among subtidal vs. intertidal populations of the marine angiosperm Zostera marina. Mol. Ecol. 16: 5156-5157.
references.citation: Oetjen, K., Ferber, S., Dankert, I. & Reusch, T.B.H. 2010. New evidence for habitat-specific selection in Wadden Sea Zostera marina populations revealed by genome scanning using SNP and microsatellite markers. Mar. Biol. 157: 81-89.
references.citation: Papa, R., Bellucci, E., Rossi, M., Leonardi, S., Rau, D., Gepts, P., Nanni, L. & Attene, G. 2007. Tagging the signatures of domestication in common bean (Phaseolus vulgaris) by means of pooled DNA samples. Ann. Bot. London 100: 1039-1051.
references.citation: Paris, M., Boyer, S., Bonin, A., Collado, A., David, J. & Despres, L. 2010. Genome scan in the mosquito Aedes rusticus: population structure and detection of positive selection after insecticide treatment. Mol. Ecol. 19: 325-337.
references.citation: Parisod, C. & Joost, S. 2010. Divergent selection in trailing-versus leading-edge populations of Biscutella laevigata. Ann. Bot. London 105: 655-660.
references.citation: Raeymaekers, J.A.M., Houdt, J.K.J.V., Larmuseau, M.H.D., Geldof, S. & Volckaert, F.A.M. 2007. Divergent selection as revealed by PST and QTL-based FST in three-spined stickleback (Gasterosteus aculeatus) populations along a coastal-inland gradient. Mol. Ecol. 16: 891-905.
references.citation: Ridgway, T., Riginos, C., Davis, J. & HoeghGuldberg, O. 2008. Genetic connectivity patterns of Pocillopora verrucosa in southern African Marine Protected Areas. Mar. Ecol. Prog. Ser. 354: 161-168.
references.citation: Robertson, A. 1975. Gene frequency distributions as a test of selective neutrality. Genetics 81: 775-785.
references.citation: Santalla, M., Ron, A.M. & La Fuente, M. 2010. Integration of genome and phenotypic scanning gives evidence of genetic structure in Mesoamerican common bean (Phaseolus vulgaris L.) landraces from the southwest of Europe. Theor. Appl. Genet. 120: 1635-1651.
references.citation: Savolainen, V., Anstett, M.-C., Lexer, C., Hutton, I., Clarkson, J.J., Norup, M.V., Powell, M.P., Springate, D., Salamin, N. & Baker, W. 2006. Sympatric speciation in palms on an oceanic island. Nat. Lett. 411: 210-213.
references.citation: Scotti-Saintagne, C., Mariette, S., Porth, I., Goicoechea, P.G., Barreneche, T., Bodènés, C., Burg, K. & Kremer, A. 2004. Genome scanning for interspecific differentiation between two closely related oak species [Quercus robur L. and Q. petraea (Matt.) Liebl.]. Genetics 168: 1615-1626.
references.citation: Smith, T.B., Milá, B., Grether, G.F., Slabbekoorn, H., Sepil, I., Buermann, W., Saatchi, S. & Pollinger, J.P. 2008. Evolutionary consequences of human disturbance in a rainforest bird species from Central Africa. Mol. Ecol. 17: 58-71.
references.citation: Sokal, R.R. & Rohlf, F.J. 1995. Biometry. Freeman & Co., New York.
references.citation: Tiranti, B. & Negri, V. 2007. Selective microenvironmental effects play a role in shaping genetic diversity and structure in a Phaseolus vulgaris L. landrace: implications for on-farm conservation. Mol. Ecol. 16: 4942-4955.
references.citation: Tsumura, Y., Kado, T., Takahashi, T., Tani, N., Ujino-Ihara, T. & Iwata, H. 2007. Genome scan to detect genetic structure and adaptive genes of natural populations of Cryptomeria japonica. Genetics 176: 2393-2403.
references.citation: Vasemägi, A., Nilsson, J. & Primmer, C.R. 2005. Expressed sequence tag-linked microsatellites as a source of gene-associated polymorphisms for detecting signatures of divergent selection in Atlantic salmon (Salmo salar L.). Mol. Biol. Evol. 22: 1067-1076.
references.citation: Vitalis, R., Dawson, K. & Boursot, P. 2001. Interpretation of variation across marker loci as evidence of selection. Genetics 158: 1811-1823.
references.citation: Vitalis, R., Dawson, K., Boursot, P. & Belkhir, K. 2003. DetSel 1.0: a computer program to detect markers responding to selection. J. Hered. 94: 429-431.
references.citation: Volis, S. 2008. Detection of signatures of positive selection in naturally occurring genetic variation. In: Population Genetics Research Progress (V.T. Koven, ed), pp. 279-310. Nova Science Publishers, Hauppauge, NY, USA.
references.citation: Weir, B.S. & Cockerham, C.C. 1984. Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370.
references.citation: Wilding, C.S., Butlin, R.K. & Grahame, J. 2001. Differential gene exchange between parapatric morphs of Littorina saxatilis detected using AFLP markers. J. Evol. Biol. 14: 611-619.
references.citation: Zhivotovsky, L.A. 1999. Estimating population structure in diploids with multilocus dominant DNA markers. Mol. Ecol. 8: 907-913.