Terminal Pleistocene Alaskan genome reveals first founding population of Native Americans

Despite broad agreement that the Americas were initially populated via Beringia, the land bridge that connected far northeast Asia with northwestern North America during the Pleistocene epoch, when and how the peopling of the Americas occurred remains unresolved. Analyses of human remains from Late Pleistocene Alaska are important to resolving the timing and dispersal of these populations. The remains of two infants were recovered at Upward Sun River (USR), and have been dated to around 11.5 thousand years ago (ka). Here, by sequencing the USR1 genome to an average coverage of approximately 17 times, we show that USR1 is most closely related to Native Americans, but falls basal to all previously sequenced contemporary and ancient Native Americans. As such, USR1 represents a distinct Ancient Beringian population. Using demographic modelling, we infer that the Ancient Beringian population and ancestors of other Native Americans descended from a single founding population that initially split from East Asians around 36 ± 1.5 ka, with gene flow persisting until around 25 ± 1.1 ka. Gene flow from ancient north Eurasians into all Native Americans took place 25–20 ka, with Ancient Beringians branching off around 22–18.1 ka. Our findings support a long-term genetic structure in ancestral Native Americans, consistent with the Beringian ‘standstill model’. We show that the basal northern and southern Native American branches, to which all other Native Americans belong, diverged around 17.5–14.6 ka, and that this probably occurred south of the North American ice sheets. We also show that after 11.5 ka, some of the northern Native American populations received gene flow from a Siberian population most closely related to Koryaks, but not Palaeo-Eskimos, Inuits or Kets, and that Native American gene flow into Inuits was through northern and not southern Native American groups. Our findings further suggest that the far-northern North American presence of northern Native Americans is from a back migration that replaced or absorbed the initial founding population of Ancient Beringians.

[1]  A. Dyke,et al.  Deglaciation of North America , 2003 .

[2]  Jonathan Terhorst,et al.  Efficient Computation of the Joint Sample Frequency Spectra for Multiple Populations , 2015, Journal of computational and graphical statistics : a joint publication of American Statistical Association, Institute of Mathematical Statistics, Interface Foundation of North America.

[3]  Cristina E. Valdiosera,et al.  Ancient mitochondrial DNA provides high-resolution time scale of the peopling of the Americas , 2016, Science Advances.

[4]  T. Korneliussen,et al.  Estimating Individual Admixture Proportions from Next Generation Sequencing Data , 2013, Genetics.

[5]  Philip L. F. Johnson,et al.  Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse , 2013, Nature.

[6]  Yun S. Song,et al.  Robust and scalable inference of population history from hundreds of unphased whole genomes , 2016, Nature Genetics.

[7]  Thomas Mailund,et al.  admixturegraph: an R package for admixture graph manipulation and fitting , 2017, Bioinform..

[8]  Stinus Lindgreen,et al.  AdapterRemoval: easy cleaning of next-generation sequencing reads , 2012, BMC Research Notes.

[9]  Cristina E. Valdiosera,et al.  The ancestry and affiliations of Kennewick Man , 2015, Nature.

[10]  B. Potter,et al.  New insights into Eastern Beringian mortuary behavior: A terminal Pleistocene double infant burial at Upward Sun River , 2014, Proceedings of the National Academy of Sciences.

[11]  Ida Moltke,et al.  NgsRelate: a software tool for estimating pairwise relatedness from next-generation sequencing data , 2015, Bioinform..

[12]  Søren Brunak,et al.  Population genomics of Bronze Age Eurasia , 2015, Nature.

[13]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[14]  Noreen von Cramon-Taubadel,et al.  Evolutionary population history of early Paleoamerican cranial morphology , 2017, Science Advances.

[15]  Ludovic Antoine Alexandre,et al.  Improving ancient DNA read mapping against modern reference genomes , 2015 .

[16]  O. Delaneau,et al.  Supplementary Information for ‘ Improved whole chromosome phasing for disease and population genetic studies ’ , 2012 .

[17]  R. Mägi,et al.  Upper Palaeolithic Siberian genome reveals dual ancestry of Native Americans , 2013, Nature.

[18]  Cris E. Hughes,et al.  Patterns of Admixture and Population Structure in Native Populations of Northwest North America , 2014, PLoS genetics.

[19]  Omar E. Cornejo,et al.  The genetic prehistory of the New World Arctic , 2014, Science.

[20]  E. Vajda,et al.  Na-Dene populations descend from the Paleo-Eskimo migration into America , 2016, bioRxiv.

[21]  Mattias Jakobsson,et al.  The genome of a Late Pleistocene human from a Clovis burial site in western Montana , 2014, Nature.

[22]  D. Reich,et al.  Genetic evidence for two founding populations of the Americas , 2015, Nature.

[23]  Eske Willerslev,et al.  Postglacial viability and colonization in North America’s ice-free corridor , 2016, Nature.

[24]  Philip L. F. Johnson,et al.  A Revised Timescale for Human Evolution Based on Ancient Mitochondrial Genomes , 2013, Current Biology.

[25]  Philip L. F. Johnson,et al.  A Draft Sequence of the Neandertal Genome , 2010, Science.

[26]  E. Røe The Oxford Handbook of The Prehistoric Arctic , 2018, Norwegian Archaeological Review.

[27]  Joseph K. Pickrell,et al.  Inference of population splits and mixtures from genome-wide allele frequency data , 2012 .

[28]  Mattias Jakobsson,et al.  Genomic evidence for the Pleistocene and recent population history of Native Americans , 2015, Science.

[29]  Jonas Binladen,et al.  Monte Verde: Seaweed, Food, Medicine, and the Peopling of South America , 2008, Science.

[30]  Philip L. F. Johnson,et al.  A Complete Neandertal Mitochondrial Genome Sequence Determined by High-Throughput Sequencing , 2008, Cell.

[31]  A. Krogh,et al.  Ancient human genome sequence of an extinct Palaeo-Eskimo , 2010, Nature.

[32]  D. O’Rourke,et al.  Beringia and the global dispersal of modern humans , 2016, Evolutionary anthropology.

[33]  David H. Alexander,et al.  Fast model-based estimation of ancestry in unrelated individuals. , 2009, Genome research.

[34]  David B. Witonsky,et al.  Reconstructing Native American Population History , 2012, Nature.

[35]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[36]  Yong Wang,et al.  bammds: a tool for assessing the ancestry of low-depth whole-genome data using multidimensional scaling (MDS) , 2014, Bioinform..

[37]  Anders Albrechtsen,et al.  ANGSD: Analysis of Next Generation Sequencing Data , 2014, BMC Bioinformatics.

[38]  Anders Albrechtsen,et al.  Relatedness mapping and tracts of relatedness for genome‐wide data in the presence of linkage disequilibrium , 2009, Genetic epidemiology.

[39]  Swapan Mallick,et al.  Ancient Admixture in Human History , 2012, Genetics.

[40]  J. Dipierri,et al.  Beringian Standstill and Spread of Native American Founders , 2007, PloS one.

[41]  Yun S. Song,et al.  Inference of complex population histories using whole-genome sequences from multiple populations , 2015, Proceedings of the National Academy of Sciences.

[42]  T. Goebel The “Microblade Adaptation” and Recolonization of Siberia during the Late Upper Pleistocene , 2008 .

[43]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[44]  D. O’Rourke,et al.  Two contemporaneous mitogenomes from terminal Pleistocene burials in eastern Beringia , 2015, Proceedings of the National Academy of Sciences.

[45]  M. Slatkin,et al.  Joint Estimation of Contamination, Error and Demography for Nuclear DNA from Ancient Humans , 2015, bioRxiv.