Screening for epistatic selection signatures: A simulation study

Detecting combinations of alleles that diverged between subpopulations via selection signature statistics can contribute to decipher the phenomenon of epistasis. This research focused on the simulation of genomic data from subpopulations under divergent epistatic selection (ES). We used D’IS2 and FST statistics in pairs of loci to scan the whole-genome. The results showed the ability to identify loci under additive-by-additive ES (ESaa) by reporting large statistical departures between subpopulations with a high level of divergence, while it did not show the same advantage in the other types of ES. Despite this, limitations such as the difficulty to distinguish between the quasi-complete fixation of one locus by ESaa from other events were observed. However, D’IS2 can detect loci under ESaa by defining a minimum boundary for the minor allele frequency on a multiple subpopulation analysis where ES only takes place in one subset. Even so, the major limitation was distinguishing between ES and single-locus selection (SS); therefore, we can conclude that divergent locus can be also a result of ES. The test conditions with D-statistics of both Ohta (1982a, 1982b) and Black and Krafsur (1985) did not provide evidence to differentiate ES in our simulation framework of isolated subpopulations.