The aggregate site frequency spectrum for comparative population genomic inference

Author(s)

Xue, Alexander T. ^{[1, 2]} ; Hickerson, Michael J. ^[3]

Affiliation

^[1] CUNY City Coll, Dept Biol Subprogram Ecol Evolutionary Biol & Beh, New York, NY 10031 - USA

^[2] CUNY, Grad Ctr, New York, NY 10031 - USA

^[3] Amer Museum Nat Hist, Div Invertebrate Zool, New York, NY 10024 - USA

• Document type: Journal article
• Source: Molecular Ecology; v. 24, n. 24, p. 6223-6240, DEC 2015.
• Web of Science Citations: 23

Abstract

Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co-distributed species for assemblage-level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large-scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next-generation comparative population genomic inference. (AU)

• FAPESP's process: 13/50297-0 - Dimensions US-BIOTA São Paulo: a multidisciplinary framework for biodiversity prediction in the Brazilian Atlantic forest hotspot
• Grantee: Cristina Yumi Miyaki
• Support type: BIOTA-FAPESP Program - Thematic Grants