A Spatial Framework for Understanding Population Structure and Admixture.

Abstract

Geographic patterns of genetic variation within modern populations, produced by complex histories of migration, can be difficult to infer and visually summarize. A general consequence of geographically limited dispersal is that samples from nearby locations tend to be more closely related than samples from distant locations, and so genetic covariance often recapitulates geographic proximity. We use genome-wide polymorphism data to build "geogenetic maps," which, when applied to stationary populations, produces a map of the geographic positions of the populations, but with distances distorted to reflect historical rates of gene flow. In the underlying model, allele frequency covariance is a decreasing function of geogenetic distance, and nonlocal gene flow such as admixture can be identified as anomalously strong covariance over long distances. This admixture is explicitly co-estimated and depicted as arrows, from the source of admixture to the recipient, on the geogenetic map. We demonstrate the utility of this method on a circum-Tibetan sampling of the greenish warbler (Phylloscopus trochiloides), in which we find evidence for gene flow between the adjacent, terminal populations of the ring species. We also analyze a global sampling of human populations, for which we largely recover the geography of the sampling, with support for significant histories of admixture in many samples. This new tool for understanding and visualizing patterns of population structure is implemented in a Bayesian framework in the program SpaceMix. Author SummaryIn this paper, we introduce a statistical method for inferring, for a set of sequenced samples, a map in which the distances between population locations reflect genetic, rather than geographic, proximity. Two populations that are sampled at distant locations but that are genetically similar (perhaps one was recently founded by a colonization event from the other) may have inferred locations that are nearby, while two populations that are sampled close together, but that are genetically dissimilar (e.g., are separated by a barrier), may have inferred locations that are farther apart. The result is a "geogenetic" map in which the distances between populations are effective distances, indicative of the way that PLOS Genetics |

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

“…Whole-genome relationships, as summarized using principal components analysis (PCA; Fig. S1, Supporting information), show the pattern expected based on previous research (Alcaide et al 2014; see also Bradburd et al 2016), of two highly distinct Siberian forms (viridanus and plumbeitarsus) and a progression of genomic signatures through the ring of populations to the south. Note that Alcaide et al (2014) summarized variation in the same GBS reads, but used an entirely distinct bioinformatics pipeline and based their PCA on only 2334 SNPs due to very restrictive filtering; the fact that the current study recovers similar patterns using much less restrictive filtering and approximately 250 times the number of SNPs gives strong confidence in the inferred relationships.…”

Ancestral and neo‐sex chromosomes contribute to population divergence in a dioecious plant

“…As our sampling design did not include the purported cytotype contact zone (shaded area in Fig. 1C; Smith 1969), and thus created a discontinuity in sampling, clustering may be a product of the interaction between incomplete sampling and isolation-bydistance rather than representing a discrete barrier to gene flow (see Balkenhol et al 2015;Bradburd et al 2016;Perez et al 2018). To evaluate this possibility, we examined the estimated effective migration surface (EEMS; Petkova et al 2016), which is a Bayesian estimate of effective migration rates based on SNP frequencies and a prior hypothesis of isolation-by-distance.…”

Phylogenomics at the tips: inferring lineages and their demographic history in a tropical lizard, Carlia amax

“…We applied an alternate method using the program SPACEMIX (Bradburd et al 2015) which estimates allele frequency covariances associated with geography to explore spatial patterns of genetic variation and explore the extent of admixture between lineages. SPACEMIX provides a better summary of the data than PCoA because it accounts for geographic proximity in genetic covariance estimates (Bradburd et al 2015). High covariances over large distances are indicative of nonlocal gene flow (admixture).…”