Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Jan 9;16(1):e2003703.
doi: 10.1371/journal.pbio.2003703. eCollection 2018 Jan.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation

Affiliations

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation

Torsten Günther et al. PLoS Biol. .

Abstract

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500-6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east-west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Mesolithic samples and their genetic affinities.
(A) Map of the Mesolithic European samples used in this study. The pie charts show the model-based [18,19] estimates of genetic ancestry for each SHG individual. The map also displays the ice sheet covering Scandinavia 10,000 cal BP (most credible [solid line] and maximum extend [dashed line] following [10]). Newly sequenced individuals are shown with bold and italic site names. SF11 is excluded from this map due to its low coverage (0.1×). Additional European EHG and WHG individuals used in this study derive from sites outside this map. The map was plotted using the R package rworldmap [28]. (B) Magnified section of genetic similarity among ancient and modern day individuals using PCA, featuring only the Mesolithic European samples (see S6 Text for the full plot). Symbols representing newly sequenced individuals have a black contour line. (C) Allele sharing between the SHGs, Latvian Mesolithic hunter-gatherers (Zv) [29], and EHGs versus WHGs measured by the statistic f4(Chimp, SHG; EHG, WHG) calculated for the captured SNPs [20]. Error bars show two block-jackknife standard errors. Data shown in this figure can be found in S1 Data. BP, before present; cal, calibrated; Chimp, Chimpanzee; EHG, eastern hunter-gatherer; PCA, principal component analysis; SHG, Scandinavian hunter-gatherer; WHG, western hunter-gatherer; Zv, Latvian Mesolithic hunter-gatherer from Zvejnieki.
Fig 2
Fig 2. Migration scenarios into postglacial Scandinavia.
Maps showing potential migration routes into Scandinavia. Scenario (a) shows a migration related to the Ahrensburgian tradition from the south (S1 Text). Scenarios (b), (c), and (d) show different possible routes into Scandinavia for the EHG ancestry. The scenarios are discussed in the text and the scenario most consistent with genetic data and stone tools is a combination of routes (a) and (b). All maps were plotted using the R package rworldmap [28]. EHG, eastern hunter-gatherer.
Fig 3
Fig 3. Genetic diversity in prehistoric Europe.
(A) RoH for the six prehistoric humans that have been sequenced to >15× genome coverage, (Kotias is a hunter-gatherer from the Caucasus region [35], NE1 is an early Neolithic individual from modern-day Hungary [27], the other individuals are described in the text), compared to all modern-day non African individuals from the 1000 Genomes Project [32]. (B) LD decay for five prehistoric populations each represented by two individuals (eastern SHGs: SF [SF9 and SF12], western SHGs: Hum [Hum1 and Hum2], CHGs [35]: [Kotias and Satsurblia], WHGs [18,35] [Loschbour and Bichon], and early Neolithic Hungarians [27]: EN_Hungary [NE1 and NE6]). LD was scaled in each distance bin by using the LD for two modern populations [32] as 0 (TSI) and as 1 (PEL). LD was calculated from the covariance of derived allele frequencies of two haploid individuals per population (S7 Text). Error bars show two standard errors estimated during 100 bootstraps across SNP pairs. (C) Effective population size over time as inferred by PSMC’ [44] for four prehistoric humans with high genome coverage. The dashed lines show the effective population sizes for selected modern-day populations. All curves for prehistoric individuals were shifted along the x-axis according to their radiocarbon date. S4 Fig. shows 100 bootstrap replicates per individual. Data shown in this figure can be found in S1 Data. BP, before present; CHG, Caucasus hunter-gatherer; LD, linkage disequilibrium; PEL, modern-day Peruvian individual; PSMC’, pairwise sequentially Markovian coalescent; RoH, runs of homozygosity; SHG, Scandinavian hunter-gatherer; TSI, modern-day Tuscan individual; WHG, western hunter-gatherer.
Fig 4
Fig 4. Adaptation to high-latitude environments.
(A) Plot of similarity between Mesolithic allele frequency and FIN allele frequency in contrast to difference to TSI allele frequency using the statistic Dsel. The figure shows all positive Z scores representing the number of standard deviations each SNP deviates from the mean. The green-highlighted SNPs are all located in the TMEM131 gene. The plot was made with qqman [49]. (B) Derived allele frequencies for three pigmentation-associated SNPs (SLC24A5, SLC45A2, which are associated with skin pigmentation, and OCA2/HERC2, which is associated with eye pigmentation). The dashed line connecting EHG and WHG represents potential allele frequencies if SHG were a linear combination of admixture between EHG and WHG. The solid horizontal line represents the derived allele frequency in SHG. The blue symbols representing SHGs were set on the average genome-wide WHG and EHG mixture proportion (on x-axis) across all SHGs, and the thick black line represents the minimum and maximum admixture proportions across all SHGs. Dashed horizontal lines represent modern European populations (CEU). The p-values were estimated from simulations of SHG allele frequencies based on their genome-wide ancestry proportions (S9 Text). Data shown in this figure can be found in S1 Data. CEU, Utah residents with Central European ancestry; EHG, eastern hunter-gatherer; FIN, modern-day Finnish individual; SHG, Scandinavian hunter-gatherer; TSI, modern-day Tuscan individual; WHG, western hunter-gatherer.

Comment in

  • Northwest passage to Scandinavia.
    Skoglund P. Skoglund P. Nat Ecol Evol. 2018 Apr;2(4):593-594. doi: 10.1038/s41559-018-0505-7. Nat Ecol Evol. 2018. PMID: 29483678 No abstract available.

Similar articles

Cited by

  • Neandertal ancestry through time: Insights from genomes of ancient and present-day humans.
    Iasi LNM, Chintalapati M, Skov L, Mesa AB, Hajdinjak M, Peter BM, Moorjani P. Iasi LNM, et al. bioRxiv [Preprint]. 2024 May 13:2024.05.13.593955. doi: 10.1101/2024.05.13.593955. bioRxiv. 2024. PMID: 38798350 Free PMC article. Preprint.
  • Genomic ancestry and social dynamics of the last hunter-gatherers of Atlantic France.
    Simões LG, Peyroteo-Stjerna R, Marchand G, Bernhardsson C, Vialet A, Chetty D, Ala��amlı E, Edlund H, Bouquin D, Dina C, Garmond N, Günther T, Jakobsson M. Simões LG, et al. Proc Natl Acad Sci U S A. 2024 Mar 5;121(10):e2310545121. doi: 10.1073/pnas.2310545121. Epub 2024 Feb 26. Proc Natl Acad Sci U S A. 2024. PMID: 38408241 Free PMC article.
  • The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes.
    Mallick S, Micco A, Mah M, Ringbauer H, Lazaridis I, Olalde I, Patterson N, Reich D. Mallick S, et al. Sci Data. 2024 Feb 10;11(1):182. doi: 10.1038/s41597-024-03031-7. Sci Data. 2024. PMID: 38341426 Free PMC article.
  • Stable population structure in Europe since the Iron Age, despite high mobility.
    Antonio ML, Weiß CL, Gao Z, Sawyer S, Oberreiter V, Moots HM, Spence JP, Cheronet O, Zagorc B, Praxmarer E, Özdoğan KT, Demetz L, Gelabert P, Fernandes D, Lucci M, Alihodžić T, Amrani S, Avetisyan P, Baillif-Ducros C, Bedić Ž, Bertrand A, Bilić M, Bondioli L, Borówka P, Botte E, Burmaz J, Bužanić D, Candilio F, Cvetko M, De Angelis D, Drnić I, Elschek K, Fantar M, Gaspari A, Gasperetti G, Genchi F, Golubović S, Hukeľová Z, Jankauskas R, Vučković KJ, Jeremić G, Kaić I, Kazek K, Khachatryan H, Khudaverdyan A, Kirchengast S, Korać M, Kozlowski V, Krošláková M, Kušan Špalj D, La Pastina F, Laguardia M, Legrand S, Leleković T, Leskovar T, Lorkiewicz W, Los D, Silva AM, Masaryk R, Matijević V, Cherifi YMS, Meyer N, Mikić I, Miladinović-Radmilović N, Milošević Zakić B, Nacouzi L, Natuniewicz-Sekuła M, Nava A, Neugebauer-Maresch C, Nováček J, Osterholtz A, Paige J, Paraman L, Pieri D, Pieta K, Pop-Lazić S, Ruttkay M, Sanader M, Sołtysiak A, Sperduti A, Stankovic Pesterac T, Teschler-Nicola M, Teul I, Tončinić D, Trapp J, Vulović D, Waliszewski T, Walter D, Živanović M, Filah MEM, Čaušević-Bully M, Šlaus M, Borić D, Novak M, Coppa A, Pinhasi R, Pritchard JK. Antonio ML, et al. Elife. 2024 Jan 30;13:e79714. doi: 10.7554/eLife.79714. Elife. 2024. PMID: 38288729 Free PMC article.
  • Reconstructing DNA methylation maps of ancient populations.
    Barouch A, Mathov Y, Meshorer E, Yakir B, Carmel L. Barouch A, et al. Nucleic Acids Res. 2024 Feb 28;52(4):1602-1612. doi: 10.1093/nar/gkad1232. Nucleic Acids Res. 2024. PMID: 38261973 Free PMC article.

References

    1. Andersen B, Borns HW. The ice age world: an introduction to quaternary history and research with emphasis on North America and Northern Europe during the last 2.5 million years. Oslo [u.a]: Scandinavian Univ. Press; 1994.
    1. Parducci L, Jorgensen T, Tollefsrud MM, Elverland E, Alm T, Fontana SL, et al. Glacial Survival of Boreal Trees in Northern Scandinavia. Science. 2012;335: 1083–1086. doi: 10.1126/science.1216043 - DOI - PubMed
    1. François O, Blum MGB, Jakobsson M, Rosenberg NA. Demographic History of European Populations of Arabidopsis thaliana. PLoS Genet. 2008;4: e1000075 doi: 10.1371/journal.pgen.1000075 - DOI - PMC - PubMed
    1. Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405: 907–913. doi: 10.1038/35016000 - DOI - PubMed
    1. Xenikoudakis G, Ersmark E, Tison J-L, Waits L, Kindberg J, Swenson JE, et al. Consequences of a demographic bottleneck on genetic structure and variation in the Scandinavian brown bear. Mol Ecol. 2015;24: 3441–3454. doi: 10.1111/mec.13239 - DOI - PubMed

Publication types

Grants and funding

Berit Wallenberg foundation (grant number BWS2011.0090) to MF. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. European Research Council starting grant to MJ. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Wenner-Gren Foundations postdoctoral fellowship to TG. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Knut and Alice Wallenberg foundation to MJ, JS, AG. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Riksbankens Jubileumsfond to MJ, JS, AG. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Swedish Research council (grant number 421-2013-730) to JA, JS. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Swedish Research council (grant number 2013-1905) to MJ, JS, AG. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Formas (grant number 2011-1138) to EMS. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

LinkOut - more resources