Estimating mobility using sparse data: Application to human genetic variation

Liisa Loog; Marta Mirazón Lahr; Mirna Kovacevic; Andrea Manica; Anders Eriksson; Mark G. Thomas

doi:10.1073/pnas.1703642114

Edited by Eske Willerslev, University of Copenhagen, Copenhagen, Denmark, and approved September 28, 2017 (received for review March 3, 2017)

Significance

Migratory activity is a critical factor in shaping processes of biological and cultural change through time. We introduce a method to estimate changes in underlying migratory activity that can be applied to genetic, morphological, or cultural data and is well-suited to samples that are sparsely distributed in space and through time. By applying this method to ancient genome data, we infer a number of changes in human mobility in Western Eurasia, including higher mobility in pre- than post-Last Glacial Maximum hunter–gatherers, and oscillations in Holocene mobility with peaks centering on the Neolithic transition and the beginnings of the Bronze Age and the Late Iron Age.

Abstract

Mobility is one of the most important processes shaping spatiotemporal patterns of variation in genetic, morphological, and cultural traits. However, current approaches for inferring past migration episodes in the fields of archaeology and population genetics lack either temporal resolution or formal quantification of the underlying mobility, are poorly suited to spatially and temporally sparsely sampled data, and permit only limited systematic comparison between different time periods or geographic regions. Here we present an estimator of past mobility that addresses these issues by explicitly linking trait differentiation in space and time. We demonstrate the efficacy of this estimator using spatiotemporally explicit simulations and apply it to a large set of ancient genomic data from Western Eurasia. We identify a sequence of changes in human mobility from the Late Pleistocene to the Iron Age. We find that mobility among European Holocene farmers was significantly higher than among European hunter–gatherers both pre- and postdating the Last Glacial Maximum. We also infer that this Holocene rise in mobility occurred in at least three distinct stages: the first centering on the well-known population expansion at the beginning of the Neolithic, and the second and third centering on the beginning of the Bronze Age and the late Iron Age, respectively. These findings suggest a strong link between technological change and human mobility in Holocene Western Eurasia and demonstrate the utility of this framework for exploring changes in mobility through space and time.

One of the major goals of population history inference is to assess the role played by past mobility in shaping patterns of genetic, phenotypic, and cultural variation. It is well recognized that the past movement of people shapes geographic patterns of genetic variation (1) and the subsequent ecological and evolutionary properties of populations (2). This is due to the fact that gene flow changes allele frequencies, shapes genetic drift, and can affect (3) or even mimic (4) natural selection processes. It is also recognized that migration activity can influence cultural evolutionary processes (5, 6). However, despite the general agreement that mobility has played an important role in shaping past and present patterns of genetic, phenotypic, and cultural variation among humans, relatively little is known about its temporal and geographic variation in the past (7).

Inferring past mobility is challenged by the sparseness and unevenness of sampling in time and space. As a result, studies of prehistorical mobility are typically limited to descriptive approaches, where major attested migration episodes or events are used as a proxy for general mobility. Data sources such as stable isotopes have enabled some quantification of mobility by allowing researchers to identify individuals within an archaeological community who have migrated into a region during their lifetime (e.g., ref. 8). The underlying logic behind this approach is that differences between isotope ratios—particularly strontium—within organisms reflect the isotope ratios acquired from the local environment (as a result of variation in underlying geology) (9). However, it is challenging to extrapolate within-community mobility rates to migration rates across larger geographic regions or over long time periods. Furthermore, isoscapes are still often poorly characterized, and isotope ratios can be relatively constant over large areas (9, 10) and so are not always informative.

Most standard population genetic tools used for quantifying population structure, such as ADMIXTURE analysis (11), f statistics (12), and TREEMIX (11), are poorly suited for estimating underlying mobility change through time. In classical population genetic analysis, estimators of migration rates between hypothesized subpopulations have been developed, including statistics such as F_ST (13). Some of these statistics have also been applied to large sets of quantitative trait data, such as variation in craniometric morphology (e.g., refs. 14 and 15). However, such statistics quantify differentiation among a set of contemporaneous samples and only inform on migration rates under idealized demographic scenarios—such as gene flow between discrete subpopulations—and are also influenced by other factors, such as subpopulation split times and population size fluctuations. Furthermore, these estimators reflect past migration between hypothesized subpopulations over large periods of time and therefore lack temporal resolution. Some researchers interpret the estimated ages and geographic distribution of clades on a phylogenetic tree of uniparental genetic systems (mtDNA or the Y chromosome) as proxies for the rate of spread of populations (e.g., ref. 16). However, such approaches do not permit a formal quantification of mobility and have been criticized as a tool for demographic inference (17 ⇓–19).

Thus, existing methods allow us to identify migration episodes to some extent but lack the temporal resolution and formal quantification of underlying mobility, are poorly suited to spatially and temporarily sparsely sampled data, and do not permit systematic comparison between different time periods or geographic regions. To overcome these problems, we present an estimator of past mobility that is particularly suited to sparsely distributed morphological, cultural, or genetic variation data and provide a first application to a large set of genome-wide data from ancient individuals from across Western Eurasia. We define mobility as the average distance moved by entities in a given time period.

Estimating Past Migration Rates

Under a general model of identity-by-descent with modification and isolation by distance (IBD) (20, 21), trait (genetic, morphological, or cultural) differences between any two entities (individuals or populations) increase monotonically as a function of both the temporal and spatial distance between them. We therefore expect that trait differences between entities correlate with temporal as well as spatial distances. However, the extent to which spatial and temporal differences explain observed trait variation depends on the level of spatial population structure and therefore on the level of mobility. If mobility was low (i.e., strong spatial structure), then we would expect differences between entities to be more strongly correlated with space, relative to time, whereas if mobility was high we would expect time to explain a relatively larger proportion of differences between entities (because of the homogenizing effects of high mobility across space).

Given that both spatial and temporal distances are expected to correlate with trait differences among entities, a matrix combining both spatial and temporal distance information should give a stronger correlation than either matrix alone (extra correlation, EC). However, since spatial and temporal distances are measured in different units (e.g., kilometers and years), combining them requires a scaling factor (S). Here, we show that the scaling factor value (S_max) that maximizes the correlation between a trait difference matrix and a Euclidian distance matrix combining the spatial and temporal distance matrices provides an estimator of mobility over the period and region covered by the data (see Materials and Methods and Fig. 1). For convenience, we use a geometric interpretation of the scaling factor S_max as an angle, α, in the plane defined by the spatial and temporal distances [α = atan(S_max), illustrated in Fig. 1, Inset; see Materials and Methods].

Fig. 1.

Illustration of the principle of maximum time–space correlation. The black dots show a typical dependence of the correlation between genetic and time–space distances on the scaling factor angle α (in degrees). Here space alone (α = 0) is a better predictor of genetic differences than time alone (α = 90), but the best predictor (highest correlation) is found at an intermediate angle, indicated by the vertical red line. (Inset) Geometrical interpretation of the Scaling Factor (S_max) as an angle (α).

To test the reliability and the robustness of S_max in recovering information about past mobility, we simulated data under a spatiotemporally explicit 2D model, which includes simple population dynamics with population growth, density dependence, and mobility (modeled as a Gaussian random walk), and generated variation data under different mobility parameter values (see Materials and Methods). We assessed the ability of S_max to infer simulated mobility values by correlation across simulations. We found a strong, positive linear relationship between the simulated average migration distance (i.e., mobility) and values of S_max (Fig. 2, R² = 0.8), thus demonstrating the utility of this statistic as an estimator for relative mobility. However, for this result to hold, it is important that the trait differences are generated under an approximately constant mutation rate and vary neutrally within a population.

Fig. 2.

Correlation between simulated movement rate (d_mig) and estimated scaling factor (S_max). Each black circle represents a single simulation. The colors correspond to the density of circles (see the color scale bar). The black line shows the best linear fit between d_mig and S_max (R² = 0.8), demonstrating that the scaling factor captures the underlying mobility in the simulated world.

Migration Rates Among Pleistocene Hunter–Gatherers and Early Farmers

Recent advances in sequencing technologies have allowed genomic data retrieval from a large sample of past individuals (e.g., refs. 22 ⇓⇓⇓–26). Although these studies have not explicitly quantified underlying mobility in the past, they have suggested several periods of large-scale population turnover in Western Eurasia.

Given that the S_max statistic is able to recover information on past mobility in simulated data, we applied the method to a sample (n = 329) of previously published genome-wide genotype data covering a time period from the beginning of the Upper Paleolithic to the Iron Age to explore changes in past human mobility in Western Eurasia (see Materials and Methods). We also constructed nonparametric confidence intervals to account for date and sampling uncertainty and estimated P values for the S_max statistic by permutation under the null hypothesis of no IBD in space and time, which allowed us to quantify the robustness of our estimates and identify time periods during which data are too sparse for the S_max statistic to be informative (see Material and Methods). First, we explored the extent to which mobility differed between pre- and post-Last Glacial Maximum (LGM) hunter–gatherers (Fig. 3). We found the average (median) mobility rates to be higher (α = 18.1; 95% CI: 14.9–87.7; P = 0.08) among pre-LGM hunter–gatherers, temporally ranging from 37,000–26,000 y ago, compared with post-LGM hunter–gatherers (α = 9.9; 95% CI: 9.5–10.9; P = 0.03), temporally ranging from 19,000–5,000 y ago. We also estimated mobility rates for Holocene farmers, temporally ranging from 10,000–1,000 y ago, and found even higher values (α = 34.8; 95% CI: 33.9–35.3; P < 0.0001) than for both hunter–gatherer groups (see SI Appendix, Table S2 for full results).

Fig. 3.

Boxplot showing the mobility rate estimates (from jackknifing and date resampling) among pre-LGM hunter–gatherers temporally ranging from 37,000–26,000 y ago (n = 19), post-LGM hunter–gatherers temporally ranging from 19,000–5,000 y ago (n = 47), and Holocene farmers temporally ranging from 10,000–1,000 y ago (n = 263). The black solid lines are the medians of the distributions. The boxes represent the interquartile ranges, and the whiskers show the spans of the distributions.

Because Holocene western Eurasia is particularly well sampled for ancient genomic DNA, we performed a sliding window analysis to explore changes in mobility over the last 14,000 y in more detail (Fig. 4), using 4,000 y-wide windows to ensure sufficient temporal signal within each window. We inferred a reduction in mobility rate between 14,000 and 9,000 y ago, before the start of the Neolithic transition (Fig. 4A). However, throughout most of this period, the P values are not significant (Fig. 4B). Because of the small sample size in the windows covering this time period (Fig. 4C), there is no significant correlation between genetic and temporal distances, and as a result, we do not observe any EC and so lack power to estimate mobility (see Materials and Methods and SI Appendix, Fig. S3). We consequently treat the inferred decline in mobility in this time range with caution. Second, we infer a substantial increase in mobility centered on the beginning of the Neolithic, with a peak centered around 7,500 y ago (Fig. 4A). Notably, the inferred mobility rate does not remain at this level throughout the Holocene. Instead, we infer a Late Neolithic drop in mobility before a second increase centered on the beginning of the Bronze Age, around 5,000 y ago, then a decline in the Late Bronze Age and Early Iron Age, before a final increase centered on the Late Iron Age (Fig. 4A and SI Appendix, Table S3 for full results for each window).

Fig. 4.

Estimation of mobility through time from empirical data. (A) Relative mobility rate estimates in Western Eurasia over the last 14,000 y, using a 4,000-y sliding window (121 windows). The solid black line represents the mean α value from 10,000 date resampled iterations; the colored area represents the 95% confidence intervals of the jackknife distribution. (B) P values for each 4,000-y window under the null hypothesis of no EC, constructed by calculating the proportion of permuted datasets where the calculated EC value was as high or higher than the average EC value from the empirical dataset (see Material and Methods). The red dotted line represents the level above which 5% or more of the permuted datasets result in EC values as high or higher than the empirical dataset. (C) Sample size for each 4,000-y window, averaged over 10,000 date resampled iterations.

To validate the efficacy of our method to identify changes in migration rate on the time scales found in the empirical dataset (Fig. 4), we modified our simulations to represent a population experiencing two changes in migration rate, resulting in three episodes of constant migration rate. We observe a good correspondence between changes in S_max and the simulated migration rate (SI Appendix, Fig. S4), supporting our interpretation of the empirical results in Fig. 4.

Finally, we compare the performance of the S_max statistic to a simple IBD through time approach, where (the slope of) the linear relation between the genetic distances and geographic distances is used as an indicator of the level of past migratory activity: High levels of migration correspond to shallow IBD patterns. We observe a trend of decreasing spatial structure, consistent with the cumulative effects of a series of high migratory activity episodes over this period. However, this approach fails to recover the timing of those changes in migratory activity in specific periods (SI Appendix, Fig. S5). Our method overcomes this lack of power to identify changes in migratory activity by explicitly considering the temporal dimension of the data.

Discussion

Through spatiotemporally explicit simulations, we have shown that the S_max statistic can be used as a reliable proxy for the underlying relative mobility of individuals within a given time period and geographic region. Because our statistic is based on correlations, it is well suited for analyzing data from archeological and paleontological contexts, where preservation can vary significantly across different geographical areas and temporal ranges, and samples are commonly sparsely distributed across space and time. Nevertheless, in the extreme case of just a small number of sites from different geographic locations or temporal periods, spurious estimates of migratory activity may arise. The permutation procedure introduced in this study can be used to identify when the S_max estimator is uninformative. We choose only to consider relative changes in the value of the S_max estimator and do not attempt to interpret its values in absolute terms. This is because, while our intuition is that mutation rate and population size will not affect the relationship between absolute values of the S_max estimator and the true mobility rate, we admit the possibility that other factors may. Selection in response to ecological and environmental factors could also reduce the utility of the S_max statistic as a proxy for mobility because local selection can create confounding spatial or temporal population structures. However, this is a common problem for any analysis assuming neutral evolution and can be dealt with by focusing on putatively neutrally varying traits or loci.

The S_max statistic offers a robust alternative to existing methods for the quantification of IBD patterns in temporally heterogeneous datasets. In population genetics, correlations between trait differences and geographic distances are commonly used to infer past population structure and connectivity between populations (27). In such approaches, temporal structures in data are usually either ignored or mathematically controlled for by using partial least squares (e.g., ref. 28), but both of these practices have been criticized (29 ⇓–31), and we show that while such approaches can inform the cumulative effects of migration in terms of structure reduction, they are unable to recover temporal changes in migratory activity. Partial least squares analysis assumes that the effect of time on genetic differences can be decoupled from the effect of space, which is generally not the case. We avoid this problem by integrating space and time into a single distance measure. Finally, because the statistic contains information about both spatial and temporal structuring of the populations, it can be used as a potentially informative summary statistic in quantitative model-fitting frameworks such as Approximate Bayesian Computation (32).

Using the S_max statistic on ancient genomic data, we identified a sequence of changes in human mobility from the late Pleistocene to the Iron Age in western Eurasia. We find some support for reduced mobility in west Eurasian post-LGM hunter–gatherers compared with pre-LGM populations. The reasons for this result are, as yet, unclear, although possible explanations include reduced resource availability in pre-LGM Western Eurasia, requiring larger foraging ranges compared with post-LGM conditions (33, 34) and/or residual post-LGM population structure following recolonization of northern latitudes from LGM southern refugia (35). Using a sliding window analysis, we find some suggestion of a decline in post-LGM hunter–gatherer mobility leading up to the Neolithic transition. However, we caution against overinterpretation of this result as the estimated P values for the S_max statistic under the null hypothesis of no EC are mostly not significant. We find strong support for a rise in mobility during the Neolithic transition in western Eurasia, likely corresponding to a well-established demic expansion of farmers, originating in the Middle East and resulting in the spread of farming technologies throughout most of Western Eurasia (36 ⇓–38). This is followed by an inferred mobility decline toward the end of the Neolithic, possibly related to the terminal phase of the spread of farming across most of Western Eurasia, and increased sedentism (39, 40). We also find strong support for a rise in mobility centered on the onset of the Bronze Age. From previous ancient DNA studies, this period has been associated with large-scale migration of Eurasian steppe populations, particularly those related to the Yamnaya culture, into Central and Northern Europe (22, 23). However, the emergence of the first civilizations and the concomitant establishment of far-reaching trade networks, as well as technological innovations such as horse-based transport (41), could also contribute to this increase in mobility (42). Finally, our sliding window analysis indicates a mobility reduction centered on the Late Bronze Age and Early Iron Age, starting around 3,000 y ago, before a final increase centered on the Late Iron Age in Western Eurasia (Fig. 4A). One possible explanation for this pattern is a significant increase in trade and warfare during that period (43 ⇓–45). Overall, our analysis suggests a strong link between technological change and human mobility in Holocene Western Eurasia. However, it should be noted that we have used wide windows (4,000 y), which necessarily reduces chronological resolution.

A major strength of our method is its applicability to any set of neutrally evolving heritable traits where differences between individuals can be quantified and increase monotonically with geographic distance and temporal difference. This means that, in principle, the S_max statistics could allow the quantification of migratory activity in temporal and environmental contexts where obtaining ancient genetic data is not feasible, by using phenotypic data such as variation in cranial morphology, which has been shown to fit the pattern of neutral evolution and closely follow the patterns observed in analyses of neutral genetic data in humans (46, 47). Another exciting possibility is the quantification of movement based on cultural variation data, provided that appropriate near-neutral traits are used (e.g., refs. 48 ⇓–50). While it should not be assumed that the movement of artifacts always coincides with that of people, contrasting measures of movement based on genetics and cultural artifacts obtained under the same conceptual framework would allow quantification of demic vs. cultural diffusion processes. This might permit identification periods and regions where genetic, phenotypic, and cultural processes are coupled, or decoupled. Given its robustness and flexibility, we anticipate that the S_max estimator will be applicable to a wide range of genetic, phenotypic, and cultural traits, allowing the quantification of mobility in a wide variety of scenarios in which this type of analysis has previously been challenging.

Materials and Methods

The proposed migration rate estimator, S_max, is the value of a scaling factor combining spatial and temporal distance matrices into a single distance matrix that maximizes its correlation with a matrix of trait distances. To estimate that value, the geographical, temporal, and trait distance matrices are calculated as described below.

Geographic, Temporal, and Trait Distances.

The geographic distance between all sample pairs was calculated in kilometers using the Haversine Eq. (51) to account for the curvature of the Earth as follows:Gij=2r arcsin(sin((φj−φi)/2)2+cos(φi)cos(φj)sin((λi−λj)/2)2),[1]where G is the distance in kilometers between individuals i and j; φ_i and φ_j are the latitude coordinates of individuals i and j, respectively; λ_i and λ_j are the longitude coordinates of individuals i and j, respectively; and r is the radius of the earth in kilometers.

Temporal distances between samples were calculated as time in years between sample pairs. Previously reported date ranges based on stratigraphy or direct radiocarbon dating were used for all individuals (SI Appendix, Table S1). In all analyses, sample dates were randomly drawn from a uniform distribution corresponding to the upper and lower bounds of a time period for a given specimen.

Genetic distances were calculated as pairwise proportion of alleles that are not identical by state (pairwise heterozygosity), using the function ibs.dist from the Bioconductor package SNPstats v.1.18.0 (52) in the R statistical analysis environment v3.2.2 (53).

The S_max Estimator.

To consider the full range of scaling factors on a finite interval, we choose to represent S as the tangent of an angle α between 0° and 90°, where α = 0 corresponds to S = 0 (geographic variation alone explains the observed trait distances between entities) and α = π/2 corresponds to S = ∞ (temporal variation alone explains the observed trait distances between entities). Formally, the time–space product matrix (D) was calculated as follows:Dij=Gij2+(STij)2,[2]where i and j are the specimens considered, D is the time–space product matrix, G is the geographical distance matrix, T the temporal distance matrix (given by the difference in age of the samples), and S is the scaling factor [S = tan(α)].

To find the scaling factor, S_max, that maximizes the correlation between the trait distance matrix and D, the time–space product matrix, we calculated the Pearson correlation coefficient between these matrices for 200 (500 for the simulated data) scaling factor values (Fig. 1). The scaling factor value in the time–space product matrix that produced the strongest correlation with the trait distance matrix is recorded as S_max, the mobility estimator.

Simulation Tests.

The reliability and the robustness of the S_max statistic in recovering information about past mobility was explored using a spatiotemporally explicit simulation model. The simulation world consists of a grid of 8,000 by 8,000 demes. Each simulation starts with one entity placed in a randomly chosen deme and lasts 20,000 generations. The model simulated exponential population growth to a carrying capacity of 10,000 entities, followed by a stochastic birth–death process (54) mobility, and trait mutation. We generated spatiotemporal trait variation data under different mobility parameter values using the same S_max estimation protocols as described above for each dataset. A total of 10,000 independent replicates of the simulations and analyses were generated, and the utility of the S_max statistic in recovering information about mobility was assessed by correlation.

The migratory process was modeled as Gaussian random walks: In each generation, each entity moves independently in the x and y directions by distances picked randomly from a normal distribution with mean = 0 and SD = σ_mig. This corresponds to the average distance moved in a single step (d_mig) of π/2 σ_mig = 1.2533 σ_mig. Thus, d_mig is the parameter of interest. We choose 1,000 random values of d_mig from a uniform distribution with a range of 1–100. We modeled drift as a Moran-type birth–death process (54). At each generation, each entity undergoes binary fission with probability P = 0.1, creating a duplicate of itself at the same location. The two entities subsequently move and evolve independent of each other. When the number of entities reaches or exceeds the carrying capacity (10,000), excess entities are deleted at random among all entities present in that generation. Mutation was modeled as a one-dimensional Gaussian random walk for each trait (N_traits = 50). Each trait was assigned an initial value of 1,000, and new (mutated) values were picked from a random normal distribution with mean equal to the current value and SD fixed at 0.05.

Following a burn-in period of 10,000 generations, entities were sampled from simulations with a probability of 0.00001 at each generation. The x and y coordinates, time of sampling in generations, and the values for the 50 traits were recorded for all sampled entities.

Pairwise trait distances between all sampled entities in each of the simulated datasets were calculated using the Euclidean distance formula as follows:Mij=∑k=1n(dik−djk)2,[3]

where M_ij is the distance between the two entities i and j; d_ik and d_jk are the values of the trait k for individuals i and j respectively; and n is the number of recorded traits.

Out of 10,000 simulations, 9,866 (98.66%) resulted in an EC greater than zero. To match the simulated data with the empirical data, we filtered the simulated data based on the measured EC values and removed all simulations that produced an EC value smaller than 0.001. This resulted in 9,155 simulations being used in the correlation analysis.

To assess the reliability of the S_max statistic in recovering information about mobility, R² values were calculated for the correlation between the simulated d_mig values and their corresponding S_max values.

Human Mobility in Late Pleistocene and Holocene.

We considered genome-wide data comprising 354,199 SNPs typed in 329 West Eurasian (i.e., west of the Ural mountains) individuals (SI Appendix, Fig. S2) temporally ranging from ∼39,000–1,000 y before present (see SI Appendix, Fig. S1). We merged the overlapping SNPs typed in archaeological samples published in refs. 22 ⇓⇓⇓–26 and 55 ⇓⇓⇓⇓⇓–61 (see SI Appendix, Table S1 for list of samples and references) that met the geographic and temporal criteria described above. No additional bioinformatic processing of the data were carried out for this study.

The 329 individuals were assigned to one of following three groups based on their estimated age and subsistence strategy based on their archaeological context: pre-LGM hunter–gatherers n = 19 (temporally ranging from 39,000 y B.P. to 26,000 y B.P.), post-LGM hunter–gatherers n = 47 (temporally ranging from 19,000 y B.P. to 5,000 y B.P.), and Holocene farmers n = 263 (temporally ranging from 10,000 y B.P. to 500 y B.P.).

Sliding window analysis was performed on all individuals in the dataset postdating 16,000 y B.P. The S_max statistic was estimated for 121 overlapping 4,000-y windows, each differing by 100 y.

To take age uncertainty into account, we report the mean scaling factor angle from 10,000 replicates with sample dates randomly resampled from their age ranges. We estimated 95% confidence intervals through a jackknifing procedure in which a randomly chosen sample in each window was removed from analysis, and the 0.025 and 0.095 quantiles were calculated from the resulting distribution.

To estimate the IBD signal through time, we fitted a linear model of genetic distances as a function of geographic distances in each time window (with sample jackknifing and age resampling as before, using the lm function from the R package base version 3.2.2) (53) and reported the slope of the line.

Confidence Intervals and Robustness of S_max Estimator.

We tested the assumption that there is an IBD pattern by correlating the genetic (trait) distance matrices in all time bins and in all windows with the respective geographic distance matrices and the date-resampled temporal distance matrices and calculated the P values for these correlations. We find a positive and statistically significant IBD pattern in space in all windows (SI Appendix, Fig. S3 A and B, respectively and SI Appendix, Fig. S6). The isolation by temporal distance pattern is positive and significant for most windows, but some windows show negative correlations or are not significant. We find that these windows correspond to time periods where we observe low EC (SI Appendix, Fig. S3C) and also low P values for the EC (Fig. 4B).

To account for the uncertainty in sample ages, we calculated the scaling factor angle 10,000 times using dates resampled at random from a uniform distribution for each sample, as described above, and report the average of the scaling factor angle of the given distribution as a point estimate.

We also performed a leave-one-out analysis (10,000 replicates, combined with sample date resampling) to explore the combined effect of sampling and dating uncertainty and constructed approximate equal-tailed 95% confidence intervals for all groups and windows.

To assess the statistical significance of S_max estimates, we consider the EC—defined as the Pearson correlation coefficient between the trait difference matrix and the time–space product matrix when S = S_max, minus the Pearson correlation coefficient between the trait difference matrix and either the temporal or geographical distance matrix alone, whichever is higher.

To obtain a null-distribution of EC, we permuted trait data for individuals among the spatiotemporal sample locations 10,000 times and calculated EC for each permutation, as described above. Finally, we calculate the proportion of EC values from the permuted datasets that are equally high or higher than that obtained from the observed data. This permutation test permits assessment of how frequently the EC for the observed data are produced by chance alone or, alternatively, as the result of the method used for estimating the S_max statistic. The resultant P value is the probability of observing an equally high or higher EC value in permuted, supposedly signal-less data, and provides an indication of the information content of each dataset.

Simulated Scenario of Changing Migration Rate.

We modified our simulations to represent a population experiencing two changes in migration rate, resulting in three episodes of constant migration rate. We assumed a generation time of 25 y and chose the effective population size to be 2N_e = 10,000, standard figures in population genetic models of European populations (62). We next chose three levels of migration with relative magnitude on par with what was inferred from the empirical data: m₁ = 0.0002, m₂ = 0.01, and m₃ = 0.05. To ensure equilibrium conditions during the start of the sampling period, we discarded the first 10,000 steps of the simulation (using migration rate m₁). We then simulated a time period of 20,000 y, divided into three episodes with constant migration rate: m₁ for 25,000–15,000 y ago, m₂ for 15,000–10,000 y ago, and m₃ for the last 5,000 y of the simulation. This roughly corresponds to the time spans associated with Mesolithic hunter–gatherers, Neolithic farmers, and post-Neolithic cultures in our empirical dataset. From a population genetic point of view, whole genome data as used in the empirical estimates correspond to a large number of approximately independent replicates. Because our model does not include recombination, we accounted for this effect by increasing the sample size to 10,000 individuals. SI Appendix, Fig. S4 shows the migration rate estimation using the S_max statistic using a 4,000 y-wide sliding window.

R version 3.2.2 (53) was used for analyses throughout this manuscript. The correlations between temporal, geographic, and trait distance matrices were calculated using the mantel (method = “pearson”) function in R package Vegan version 2.3.0 (63). The permutation and bootstrap tests were performed using the function sample in the R package base version 3.2.2 (53).

Acknowledgments

We are very grateful to Robert Foley for valuable discussions during the formulation of the approach, to Mike Parker-Pearson for advice on Holocene migration processes, and to Tamsin O’Connell for advice on stable isotopes. L.L. was supported by Natural Environment Research Council, UK Grants NE/K005243/1 and NE/K003259/1 and European Research Council Grant 339941-ADAPT. M.G.T. was supported by Wellcome Trust Senior Investigator Award Grant 100719/Z/12/Z and Leverhulme Trust Grant RP2011-R-045. A.M. and A.E. were supported by European Research Council Consolidator Grant 647787-LocalAdaptation. M.M.L. was supported by European Research Council Advanced Grant 295907, In-Africa. M.K. was funded by the Engineering and Physical Sciences Research Council through the Centre for Mathematics and Physics in the Life Sciences and Experimental Biology.

Footnotes

↵¹To whom correspondence may be addressed. Email: LiisaLoog{at}gmail.com or m.thomas{at}ucl.ac.uk.

Author contributions: M.G.T. devised the approach in discussion with M.M.L.; L.L. and M.G.T. developed the method with input from A.E.; M.K., A.E., and M.G.T. developed the simulation code with input from L.L.; L.L. performed the analyses with input from A.E. and M.G.T.; and L.L., A.E., and M.G.T. wrote the paper with input from M.M.L, M.K., and A.M.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
Data deposition: The R code used for analyses is available from the GitHub repository (https://github.com/LiisaLoog/SpaceTime-Framework.git).
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1703642114/-/DCSupplemental.

This is an open access article distributed under the PNAS license.

References

↵
1. Hanski I,
2. Gilpin M
(1991) Metapopulation dynamics: Brief history and conceptual domain. Biol J Linn Soc Lond 42:3–16.
OpenUrl CrossRef
↵
1. Lande R
(1988) Genetics and demography in biological conservation. Science 241:1455–1460.
OpenUrl Abstract/FREE Full Text
↵
1. Itan Y,
2. Powell A,
3. Beaumont MA,
4. Burger J,
5. Thomas MG
(2009) The origins of lactase persistence in Europe. PLOS Comput Biol 5:e1000491.
OpenUrl CrossRef PubMed
↵
1. Klopfstein S,
2. Currat M,
3. Excoffier L
(2006) The fate of mutations surfing on the wave of a range expansion. Mol Biol Evol 23:482–490.
OpenUrl CrossRef PubMed
↵
1. Powell A,
2. Shennan S,
3. Thomas MG
(2009) Late Pleistocene demography and the appearance of modern human behavior. Science 324:1298–1301.
OpenUrl Abstract/FREE Full Text
↵
1. Kline MA,
2. Boyd R
(2010) Population size predicts technological complexity in Oceania. Proc R Soc Lond B Biol Sci 277:2559–2564.
OpenUrl CrossRef PubMed
↵
1. Cox MP,
2. Hammer MF
(2010) A question of scale: Human migrations writ large and small. BMC Biol 8:98.
OpenUrl PubMed
↵
1. Gregoricka LA
(2013) Residential mobility and social identity in the periphery: Strontium isotope analysis of archaeological tooth enamel from southeastern Arabia. J Archaeol Sci 40:452–464.
OpenUrl
↵
1. Makarewicz CA,
2. Sealy J
(2015) Dietary reconstruction, mobility, and the analysis of ancient skeletal tissues: Expanding the prospects of stable isotope research in archaeology. J Archaeol Sci 56:146–158.
OpenUrl
↵
1. Bowen GJ
(2010) Isoscapes: Spatial pattern in isotopic biogeochemistry. Annu Rev Earth Planet Sci 38:161–187.
OpenUrl CrossRef
↵
1. Pickrell JK,
2. Pritchard JK
(2012) Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet 8:e1002967.
OpenUrl CrossRef PubMed
↵
1. Patterson N, et al.
(2012) Ancient admixture in human history. Genetics 192:1065–1093.
OpenUrl Abstract/FREE Full Text
↵
1. Wright S
(1990) Evolution in Mendelian populations. 1931. Bull Math Biol 52:241–295, discussion 201–207.
OpenUrl CrossRef PubMed
↵
1. Relethford JH
(1994) Craniometric variation among modern human populations. Am J Phys Anthropol 95:53–62.
OpenUrl CrossRef PubMed
↵
1. Betti L,
2. Balloux F,
3. Hanihara T,
4. Manica A
(2010) The relative role of drift and selection in shaping the human skull. Am J Phys Anthropol 141:76–82.
OpenUrl PubMed
↵
1. Underhill PA,
2. Kivisild T
(2007) Use of y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu Rev Genet 41:539–564.
OpenUrl CrossRef PubMed
↵
1. Goldstein DB,
2. Chikhi L
(2002) Human migrations and population structure: What we know and why it matters. Annu Rev Genomics Hum Genet 3:129–152.
OpenUrl CrossRef PubMed
↵
1. Nielsen R,
2. Beaumont MA
(2009) Statistical inferences in phylogeography. Mol Ecol 18:1034–1047.
OpenUrl CrossRef PubMed
↵
1. Pinhasi R,
2. Thomas MG,
3. Hofreiter M,
4. Currat M,
5. Burger J
(2012) The genetic history of Europeans. Trends Genet 28:496–505.
OpenUrl CrossRef PubMed
↵
1. Wright S
(1943) Isolation by distance. Genetics 28:114–138.
OpenUrl FREE Full Text
↵
1. Nei M
(1972) Genetic distance between populations. Am Nat 106:283–292.
OpenUrl CrossRef
↵
1. Allentoft ME, et al.
(2015) Population genomics of Bronze Age Eurasia. Nature 522:167–172.
OpenUrl CrossRef PubMed
↵
1. Haak W, et al.
(2015) Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 522:207–211.
OpenUrl CrossRef PubMed
↵
1. Mathieson I, et al.
(2015) Genome-wide patterns of selection in 230 ancient Eurasians. Nature 528:499–503.
OpenUrl CrossRef PubMed
↵
1. Fu Q, et al.
(2016) The genetic history of Ice Age Europe. Nature 534:200–205.
OpenUrl CrossRef PubMed
↵
1. Lazaridis I, et al.
(2016) Genomic insights into the origin of farming in the ancient Near East. Nature 536:419–424.
OpenUrl CrossRef PubMed
↵
1. Manel S,
2. Schwartz MK,
3. Luikart G,
4. Taberlet P
(2003) Landscape genetics: Combining landscape ecology and population genetics. Trends Ecol Evol 18:189–197.
OpenUrl CrossRef
↵
1. Reyes-Centeno H, et al.
(2014) Genomic and cranial phenotype data support multiple modern human dispersals from Africa and a southern route into Asia. Proc Natl Acad Sci USA 111:7248–7253.
OpenUrl Abstract/FREE Full Text
↵
1. Depaulis F,
2. Orlando L,
3. Hänni C
(2009) Using classical population genetics tools with heterochroneous data: Time matters! PLoS One 4:e5541.
OpenUrl CrossRef PubMed
↵
1. Guillot G,
2. Rousset F
(2013) Dismantling the Mantel tests. Methods Ecol Evol 4:336–344.
OpenUrl CrossRef
↵
1. Skoglund P,
2. Sjödin P,
3. Skoglund T,
4. Lascoux M,
5. Jakobsson M
(2014) Investigating population history using temporal genetic differentiation. Mol Biol Evol 31:2516–2527.
OpenUrl CrossRef PubMed
↵
1. Beaumont MA,
2. Zhang W,
3. Balding DJ
(2002) Approximate Bayesian computation in population genetics. Genetics 162:2025–2035.
OpenUrl Abstract/FREE Full Text
↵
1. Donnelly P,
2. Foley RA
1. Foley RA,
2. Lahr MM
(2001) The anthropological, demographic and ecological context of human evolutionary genetics. Genes, Fossils, and Behaviour: An Integrated Approach to Human Evolution, eds Donnelly P, Foley RA (IOS Press, Omaha, NE), pp 223–245.
↵
1. Collard IF,
2. Foley RA
(2002) Latitudinal patterns and environmental determinants of recent human cultural diversity: Do humans follow biogeographical rules? Evol Ecol Res 4:371–383.
OpenUrl
↵
1. Miller R
(2012) Mapping the expansion of the Northwest Magdalenian. Quat Int 272–273:209–230.
OpenUrl
↵
1. Menozzi P,
2. Piazza A,
3. Cavalli-Sforza L
(1978) Synthetic maps of human gene frequencies in Europeans. Science 201:786–792.
OpenUrl Abstract/FREE Full Text
↵
1. Skoglund P, et al.
(2012) Origins and genetic legacy of Neolithic farmers and hunter-gatherers in Europe. Science 336:466–469.
OpenUrl Abstract/FREE Full Text
↵
1. Hofmanová Z, et al.
(2016) Early farmers from across Europe directly descended from Neolithic Aegeans. Proc Natl Acad Sci USA 113:6886–6891.
OpenUrl Abstract/FREE Full Text
↵
1. Bocquet-Appel J-P,
2. Naji S,
3. Linden MV,
4. Kozlowski JK
(2009) Detection of diffusion and contact zones of early farming in Europe from the space-time distribution of 14C dates. J Archaeol Sci 36:807–820.
OpenUrl CrossRef
↵
1. Isern N,
2. Fort J
(2012) Modelling the effect of Mesolithic populations on the slowdown of the Neolithic transition. J Archaeol Sci 39:3671–3676.
OpenUrl
↵
1. Warmuth V, et al.
(2012) Reconstructing the origin and spread of horse domestication in the Eurasian steppe. Proc Natl Acad Sci USA 109:8202–8206.
OpenUrl Abstract/FREE Full Text
↵
1. Warmuth VM, et al.
(2013) Ancient trade routes shaped the genetic structure of horses in eastern Eurasia. Mol Ecol 22:5340–5351.
OpenUrl CrossRef
↵
1. Sherratt S,
2. Sherratt A
(1993) The growth of the Mediterranean economy in the early first millennium BC. World Archaeol 24:361–378.
OpenUrl
↵
1. Collis J
(2003) The Celts: Origins, Myths & Inventions (Tempus, Stroud, UK).
↵
1. Beaujard P
(2010) From three possible Iron-Age world-systems to a single Afro-Eurasian world-system. J World Hist 21:1–43.
OpenUrl
↵
1. Roseman CC,
2. Weaver TD
(2007) Molecules versus morphology? Not for the human cranium. BioEssays 29:1185–1188.
OpenUrl CrossRef PubMed
↵
1. von Cramon-Taubadel N,
2. Weaver TD
(2009) Insights from a quantitative genetic approach to human morphological evolution. Evol Anthropol Issues News Rev 18:237–240.
OpenUrl
↵
1. Shennan S
(2000) Population, culture history, and the dynamics of culture change. Curr Anthropol 41:811–835.
OpenUrl CrossRef
↵
1. Eerkens JW,
2. Lipo CP
(2007) Cultural transmission theory and the archaeological record: Providing context to understanding variation and temporal changes in material culture. J Archaeol Res 15:239–274.
OpenUrl CrossRef
↵
1. Lycett SJ,
2. Norton CJ
(2010) A demographic model for Palaeolithic technological evolution: The case of East Asia and the Movius Line. Quat Int 211:55–65.
OpenUrl
↵
1. Sinnott R
(1984) Virtues of the haversine. Sky Telescope 68:159.
OpenUrl
↵
1. Clayton D
(2014) SnpStats: SnpMatrix and XSnpMatrix Classes and Methods. Version 1.18.0. Available at https://www.bioconductor.org/packages/release/bioc/html/snpStats.html. Accessed January 1, 2017.
↵
1. R Core Team
(2015) R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna), Version 3.2.2. Available at https://www.R-project.org/. Accessed January 1, 2017.
↵
1. Moran PAP
(1958) Random processes in genetics. Math Proc Camb Phil Soc 54:60–71.
OpenUrl CrossRef
↵
1. Keller A, et al.
(2012) New insights into the Tyrolean Iceman’s origin and phenotype as inferred by whole-genome sequencing. Nat Commun 3:698.
OpenUrl CrossRef PubMed
↵
1. Gamba C, et al.
(2014) Genome flux and stasis in a five millennium transect of European prehistory. Nat Commun 5:5257.
OpenUrl CrossRef PubMed
↵
1. Lazaridis I, et al.
(2014) Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature 513:409–413.
OpenUrl CrossRef PubMed
↵
1. Olalde I, et al.
(2014) Derived immune and ancestral pigmentation alleles in a 7,000-year-old Mesolithic European. Nature 507:225–228.
OpenUrl CrossRef PubMed
↵
1. Seguin-Orlando A, et al.
(2014) Paleogenomics. Genomic structure in Europeans dating back at least 36,200 years. Science 346:1113–1118.
OpenUrl Abstract/FREE Full Text
↵
1. Skoglund P, et al.
(2014) Genomic diversity and admixture differs for Stone-Age Scandinavian foragers and farmers. Science 344:747–750.
OpenUrl Abstract/FREE Full Text
↵
1. Jones ER, et al.
(2015) Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat Commun 6:8912.
OpenUrl CrossRef PubMed
↵
1. Gronau I,
2. Hubisz MJ,
3. Gulko B,
4. Danko CG,
5. Siepel A
(2011) Bayesian inference of ancient human demography from individual genome sequences. Nat Genet 43:1031–1034.
OpenUrl CrossRef PubMed
↵
1. Oksanen J, et al.
(2015) vegan: Community Ecology Package. Version 2.3.0. Available at CRAN.R-project.org/package=vegan. Accessed January 1, 2017.

Article Alerts

Email Article

Citation Tools

Request Permissions

Article Classifications

Biological Sciences
Evolution

Proceedings of the National Academy of Sciences: 114 (46)

Table of Contents

Submit

[1] ↵
Hanski I,
Gilpin M
(1991) Metapopulation dynamics: Brief history and conceptual domain. Biol J Linn Soc Lond 42:3–16.
OpenUrl CrossRef

[2] Hanski I,

[3] Gilpin M

[4] ↵
Lande R
(1988) Genetics and demography in biological conservation. Science 241:1455–1460.
OpenUrl Abstract/FREE Full Text

[5] Lande R

[6] ↵
Itan Y,
Powell A,
Beaumont MA,
Burger J,
Thomas MG
(2009) The origins of lactase persistence in Europe. PLOS Comput Biol 5:e1000491.
OpenUrl CrossRef PubMed

[7] Itan Y,

[8] Powell A,

[9] Beaumont MA,

[10] Burger J,

[11] Thomas MG

[12] ↵
Klopfstein S,
Currat M,
Excoffier L
(2006) The fate of mutations surfing on the wave of a range expansion. Mol Biol Evol 23:482–490.
OpenUrl CrossRef PubMed

[13] Klopfstein S,

[14] Currat M,

[15] Excoffier L

[16] ↵
Powell A,
Shennan S,
Thomas MG
(2009) Late Pleistocene demography and the appearance of modern human behavior. Science 324:1298–1301.
OpenUrl Abstract/FREE Full Text

[17] Powell A,

[18] Shennan S,

[19] Thomas MG

[20] ↵
Kline MA,
Boyd R
(2010) Population size predicts technological complexity in Oceania. Proc R Soc Lond B Biol Sci 277:2559–2564.
OpenUrl CrossRef PubMed

[21] Kline MA,

[22] Boyd R

[23] ↵
Cox MP,
Hammer MF
(2010) A question of scale: Human migrations writ large and small. BMC Biol 8:98.
OpenUrl PubMed

[24] Cox MP,

[25] Hammer MF

[26] ↵
Gregoricka LA
(2013) Residential mobility and social identity in the periphery: Strontium isotope analysis of archaeological tooth enamel from southeastern Arabia. J Archaeol Sci 40:452–464.
OpenUrl

[27] Gregoricka LA

[28] ↵
Makarewicz CA,
Sealy J
(2015) Dietary reconstruction, mobility, and the analysis of ancient skeletal tissues: Expanding the prospects of stable isotope research in archaeology. J Archaeol Sci 56:146–158.
OpenUrl

[29] Makarewicz CA,

[30] Sealy J

[31] ↵
Bowen GJ
(2010) Isoscapes: Spatial pattern in isotopic biogeochemistry. Annu Rev Earth Planet Sci 38:161–187.
OpenUrl CrossRef

[32] Bowen GJ

[33] ↵
Pickrell JK,
Pritchard JK
(2012) Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet 8:e1002967.
OpenUrl CrossRef PubMed

[34] Pickrell JK,

[35] Pritchard JK

[36] ↵
Patterson N, et al.
(2012) Ancient admixture in human history. Genetics 192:1065–1093.
OpenUrl Abstract/FREE Full Text

[37] Patterson N, et al.

[38] ↵
Wright S
(1990) Evolution in Mendelian populations. 1931. Bull Math Biol 52:241–295, discussion 201–207.
OpenUrl CrossRef PubMed

[39] Wright S

[40] ↵
Relethford JH
(1994) Craniometric variation among modern human populations. Am J Phys Anthropol 95:53–62.
OpenUrl CrossRef PubMed

[41] Relethford JH

[42] ↵
Betti L,
Balloux F,
Hanihara T,
Manica A
(2010) The relative role of drift and selection in shaping the human skull. Am J Phys Anthropol 141:76–82.
OpenUrl PubMed

[43] Betti L,

[44] Balloux F,

[45] Hanihara T,

[46] Manica A

[47] ↵
Underhill PA,
Kivisild T
(2007) Use of y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu Rev Genet 41:539–564.
OpenUrl CrossRef PubMed

[48] Underhill PA,

[49] Kivisild T

[50] ↵
Goldstein DB,
Chikhi L
(2002) Human migrations and population structure: What we know and why it matters. Annu Rev Genomics Hum Genet 3:129–152.
OpenUrl CrossRef PubMed

[51] Goldstein DB,

[52] Chikhi L

[53] ↵
Nielsen R,
Beaumont MA
(2009) Statistical inferences in phylogeography. Mol Ecol 18:1034–1047.
OpenUrl CrossRef PubMed

[54] Nielsen R,

[55] Beaumont MA

[56] ↵
Pinhasi R,
Thomas MG,
Hofreiter M,
Currat M,
Burger J
(2012) The genetic history of Europeans. Trends Genet 28:496–505.
OpenUrl CrossRef PubMed

[57] Pinhasi R,

[58] Thomas MG,

[59] Hofreiter M,

[60] Currat M,

[61] Burger J

[62] ↵
Wright S
(1943) Isolation by distance. Genetics 28:114–138.
OpenUrl FREE Full Text

[63] Wright S

[64] ↵
Nei M
(1972) Genetic distance between populations. Am Nat 106:283–292.
OpenUrl CrossRef

[65] Nei M

[66] ↵
Allentoft ME, et al.
(2015) Population genomics of Bronze Age Eurasia. Nature 522:167–172.
OpenUrl CrossRef PubMed

[67] Allentoft ME, et al.

[68] ↵
Haak W, et al.
(2015) Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 522:207–211.
OpenUrl CrossRef PubMed

[69] Haak W, et al.

[70] ↵
Mathieson I, et al.
(2015) Genome-wide patterns of selection in 230 ancient Eurasians. Nature 528:499–503.
OpenUrl CrossRef PubMed

[71] Mathieson I, et al.

[72] ↵
Fu Q, et al.
(2016) The genetic history of Ice Age Europe. Nature 534:200–205.
OpenUrl CrossRef PubMed

[73] Fu Q, et al.

[74] ↵
Lazaridis I, et al.
(2016) Genomic insights into the origin of farming in the ancient Near East. Nature 536:419–424.
OpenUrl CrossRef PubMed

[75] Lazaridis I, et al.

[76] ↵
Manel S,
Schwartz MK,
Luikart G,
Taberlet P
(2003) Landscape genetics: Combining landscape ecology and population genetics. Trends Ecol Evol 18:189–197.
OpenUrl CrossRef

[77] Manel S,

[78] Schwartz MK,

[79] Luikart G,

[80] Taberlet P

[81] ↵
Reyes-Centeno H, et al.
(2014) Genomic and cranial phenotype data support multiple modern human dispersals from Africa and a southern route into Asia. Proc Natl Acad Sci USA 111:7248–7253.
OpenUrl Abstract/FREE Full Text

[82] Reyes-Centeno H, et al.

[83] ↵
Depaulis F,
Orlando L,
Hänni C
(2009) Using classical population genetics tools with heterochroneous data: Time matters! PLoS One 4:e5541.
OpenUrl CrossRef PubMed

[84] Depaulis F,

[85] Orlando L,

[86] Hänni C

[87] ↵
Guillot G,
Rousset F
(2013) Dismantling the Mantel tests. Methods Ecol Evol 4:336–344.
OpenUrl CrossRef

[88] Guillot G,

[89] Rousset F

[90] ↵
Skoglund P,
Sjödin P,
Skoglund T,
Lascoux M,
Jakobsson M
(2014) Investigating population history using temporal genetic differentiation. Mol Biol Evol 31:2516–2527.
OpenUrl CrossRef PubMed

[91] Skoglund P,

[92] Sjödin P,

[93] Skoglund T,

[94] Lascoux M,

[95] Jakobsson M

[96] ↵
Beaumont MA,
Zhang W,
Balding DJ
(2002) Approximate Bayesian computation in population genetics. Genetics 162:2025–2035.
OpenUrl Abstract/FREE Full Text

[97] Beaumont MA,

[98] Zhang W,

[99] Balding DJ

[100] ↵
Donnelly P,
Foley RA
Foley RA,
Lahr MM
(2001) The anthropological, demographic and ecological context of human evolutionary genetics. Genes, Fossils, and Behaviour: An Integrated Approach to Human Evolution, eds Donnelly P, Foley RA (IOS Press, Omaha, NE), pp 223–245.

[101] Donnelly P,

[102] Foley RA

[103] Foley RA,

[104] Lahr MM

[105] ↵
Collard IF,
Foley RA
(2002) Latitudinal patterns and environmental determinants of recent human cultural diversity: Do humans follow biogeographical rules? Evol Ecol Res 4:371–383.
OpenUrl

[106] Collard IF,

[107] Foley RA

[108] ↵
Miller R
(2012) Mapping the expansion of the Northwest Magdalenian. Quat Int 272–273:209–230.
OpenUrl

[109] Miller R

[110] ↵
Menozzi P,
Piazza A,
Cavalli-Sforza L
(1978) Synthetic maps of human gene frequencies in Europeans. Science 201:786–792.
OpenUrl Abstract/FREE Full Text

[111] Menozzi P,

[112] Piazza A,

[113] Cavalli-Sforza L

[114] ↵
Skoglund P, et al.
(2012) Origins and genetic legacy of Neolithic farmers and hunter-gatherers in Europe. Science 336:466–469.
OpenUrl Abstract/FREE Full Text

[115] Skoglund P, et al.

[116] ↵
Hofmanová Z, et al.
(2016) Early farmers from across Europe directly descended from Neolithic Aegeans. Proc Natl Acad Sci USA 113:6886–6891.
OpenUrl Abstract/FREE Full Text

[117] Hofmanová Z, et al.

[118] ↵
Bocquet-Appel J-P,
Naji S,
Linden MV,
Kozlowski JK
(2009) Detection of diffusion and contact zones of early farming in Europe from the space-time distribution of 14C dates. J Archaeol Sci 36:807–820.
OpenUrl CrossRef

[119] Bocquet-Appel J-P,

[120] Naji S,

[121] Linden MV,

[122] Kozlowski JK

[123] ↵
Isern N,
Fort J
(2012) Modelling the effect of Mesolithic populations on the slowdown of the Neolithic transition. J Archaeol Sci 39:3671–3676.
OpenUrl

[124] Isern N,

[125] Fort J

[126] ↵
Warmuth V, et al.
(2012) Reconstructing the origin and spread of horse domestication in the Eurasian steppe. Proc Natl Acad Sci USA 109:8202–8206.
OpenUrl Abstract/FREE Full Text

[127] Warmuth V, et al.

[128] ↵
Warmuth VM, et al.
(2013) Ancient trade routes shaped the genetic structure of horses in eastern Eurasia. Mol Ecol 22:5340–5351.
OpenUrl CrossRef

[129] Warmuth VM, et al.

[130] ↵
Sherratt S,
Sherratt A
(1993) The growth of the Mediterranean economy in the early first millennium BC. World Archaeol 24:361–378.
OpenUrl

[131] Sherratt S,

[132] Sherratt A

[133] ↵
Collis J
(2003) The Celts: Origins, Myths & Inventions (Tempus, Stroud, UK).

[134] Collis J

[135] ↵
Beaujard P
(2010) From three possible Iron-Age world-systems to a single Afro-Eurasian world-system. J World Hist 21:1–43.
OpenUrl

[136] Beaujard P

[137] ↵
Roseman CC,
Weaver TD
(2007) Molecules versus morphology? Not for the human cranium. BioEssays 29:1185–1188.
OpenUrl CrossRef PubMed

[138] Roseman CC,

[139] Weaver TD

[140] ↵
von Cramon-Taubadel N,
Weaver TD
(2009) Insights from a quantitative genetic approach to human morphological evolution. Evol Anthropol Issues News Rev 18:237–240.
OpenUrl

[141] von Cramon-Taubadel N,

[142] Weaver TD

[143] ↵
Shennan S
(2000) Population, culture history, and the dynamics of culture change. Curr Anthropol 41:811–835.
OpenUrl CrossRef

[144] Shennan S

[145] ↵
Eerkens JW,
Lipo CP
(2007) Cultural transmission theory and the archaeological record: Providing context to understanding variation and temporal changes in material culture. J Archaeol Res 15:239–274.
OpenUrl CrossRef

[146] Eerkens JW,

[147] Lipo CP

[148] ↵
Lycett SJ,
Norton CJ
(2010) A demographic model for Palaeolithic technological evolution: The case of East Asia and the Movius Line. Quat Int 211:55–65.
OpenUrl

[149] Lycett SJ,

[150] Norton CJ

[151] ↵
Sinnott R
(1984) Virtues of the haversine. Sky Telescope 68:159.
OpenUrl

[152] Sinnott R

[153] ↵
Clayton D
(2014) SnpStats: SnpMatrix and XSnpMatrix Classes and Methods. Version 1.18.0. Available at https://www.bioconductor.org/packages/release/bioc/html/snpStats.html. Accessed January 1, 2017.

[154] Clayton D

[155] ↵
R Core Team
(2015) R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna), Version 3.2.2. Available at https://www.R-project.org/. Accessed January 1, 2017.

[156] R Core Team

[157] ↵
Moran PAP
(1958) Random processes in genetics. Math Proc Camb Phil Soc 54:60–71.
OpenUrl CrossRef

[158] Moran PAP

[159] ↵
Keller A, et al.
(2012) New insights into the Tyrolean Iceman’s origin and phenotype as inferred by whole-genome sequencing. Nat Commun 3:698.
OpenUrl CrossRef PubMed

[160] Keller A, et al.

[161] ↵
Gamba C, et al.
(2014) Genome flux and stasis in a five millennium transect of European prehistory. Nat Commun 5:5257.
OpenUrl CrossRef PubMed

[162] Gamba C, et al.

[163] ↵
Lazaridis I, et al.
(2014) Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature 513:409–413.
OpenUrl CrossRef PubMed

[164] Lazaridis I, et al.

[165] ↵
Olalde I, et al.
(2014) Derived immune and ancestral pigmentation alleles in a 7,000-year-old Mesolithic European. Nature 507:225–228.
OpenUrl CrossRef PubMed

[166] Olalde I, et al.

[167] ↵
Seguin-Orlando A, et al.
(2014) Paleogenomics. Genomic structure in Europeans dating back at least 36,200 years. Science 346:1113–1118.
OpenUrl Abstract/FREE Full Text

[168] Seguin-Orlando A, et al.

[169] ↵
Skoglund P, et al.
(2014) Genomic diversity and admixture differs for Stone-Age Scandinavian foragers and farmers. Science 344:747–750.
OpenUrl Abstract/FREE Full Text

[170] Skoglund P, et al.

[171] ↵
Jones ER, et al.
(2015) Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat Commun 6:8912.
OpenUrl CrossRef PubMed

[172] Jones ER, et al.

[173] ↵
Gronau I,
Hubisz MJ,
Gulko B,
Danko CG,
Siepel A
(2011) Bayesian inference of ancient human demography from individual genome sequences. Nat Genet 43:1031–1034.
OpenUrl CrossRef PubMed

[174] Gronau I,

[175] Hubisz MJ,

[176] Gulko B,

[177] Danko CG,

[178] Siepel A

[179] ↵
Oksanen J, et al.
(2015) vegan: Community Ecology Package. Version 2.3.0. Available at CRAN.R-project.org/package=vegan. Accessed January 1, 2017.

[180] Oksanen J, et al.

Main menu

User menu

Search

Estimating mobility using sparse data: Application to human genetic variation

Significance

Abstract

Estimating Past Migration Rates

Migration Rates Among Pleistocene Hunter–Gatherers and Early Farmers

Discussion

Materials and Methods

Geographic, Temporal, and Trait Distances.

The S_max Estimator.

Simulation Tests.

Human Mobility in Late Pleistocene and Holocene.

Confidence Intervals and Robustness of S_max Estimator.

Simulated Scenario of Changing Migration Rate.

Acknowledgments

Footnotes

References

Citation Manager Formats

Article Classifications

You May Also be Interested in

Similar Articles

Articles

PNAS Portals

Information

Main menu

User menu

Search

Significance

Abstract

Estimating Past Migration Rates

Migration Rates Among Pleistocene Hunter–Gatherers and Early Farmers

Discussion

Materials and Methods

Geographic, Temporal, and Trait Distances.

The Smax Estimator.

Simulation Tests.

Human Mobility in Late Pleistocene and Holocene.

Confidence Intervals and Robustness of Smax Estimator.

Simulated Scenario of Changing Migration Rate.

Acknowledgments

Footnotes

References

Citation Manager Formats

Article Classifications

Sign up for Article Alerts

Jump to section

You May Also be Interested in

Similar Articles

Articles

PNAS Portals

Information

The S_max Estimator.

Confidence Intervals and Robustness of S_max Estimator.