Cited by 57

The DOI matches its Crossref record; 18 of the 18 references it cites were found in a registry or an index.

Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry

Flegontov, Pavel, Changmai, Piya, Zidkova, Anastassiya, Logacheva, Maria D., Altınışık, N. Ezgi, Flegontova, Olga, Gelfand, Mikhail S., Gerasimov, Evgeny S., Khrameeva, Ekaterina E., Konovalova, Olga P., Neretina, Tatiana, Nikolsky, Yuri V., Starostin, George, Stepanova, Vita V., Travinsky, Igor V., Tříska, Martin, Tříska, Petr, Tatarinova, Tatiana V.·2016·Scientific Reports

Published version of a preprint

Abstract

The Kets, an ethnic group in the Yenisei River basin, Russia, are considered the last nomadic hunter-gatherers of Siberia and Ket language has no transparent affiliation with any language family. We investigated connections between the Kets and Siberian and North American populations, with emphasis on the Mal’ta and Paleo-Eskimo ancient genomes, using original data from 46 unrelated samples of Kets and 42 samples of their neighboring ethnic groups (Uralic-speaking Nganasans, Enets and Selkups). We genotyped over 130,000 autosomal SNPs, identified mitochondrial and Y-chromosomal haplogroups and performed high-coverage genome sequencing of two Ket individuals. We established that Nganasans, Kets, Selkups and Yukaghirs form a cluster of populations most closely related to Paleo-Eskimos in Siberia (not considering indigenous populations of Chukotka and Kamchatka). Kets are closely related to modern Selkups and to some Bronze and Iron Age populations of the Altai region, with all these groups sharing a high degree of Mal’ta ancestry. Implications of these findings for the linguistic hypothesis uniting Ket and Na-Dene languages into a language macrofamily are discussed.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/srep20768

Cites

18 of 18 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 (A) Admixture coefficients plotted for dataset ‘GenoChip + Illumina arrays’. Abbreviated names of admixture components are shown on the left as follows: SAM, South American; NAM, North American; ESK, Eskimo (Beringian); SEA, South-East Asian; SIB, Siberian; NEU, North European; ME, Middle Eastern; CAU, Caucasian; SAS, South Asian; OCE, Oceanian; AFR, African. The Ket-Uralic (‘Ket’) admixture component appears at K ≥ 11, and admixture coefficients are plotted for K = 4, 10, 11, and 19. Although K = 20 demonstrates the lowest average cross-validation error, the Ket-Uralic component splits in two at this K value, therefore K = 19 was chosen for the final analysis. Only populations containing at least one individual with >5% of the Ket-Uralic component at K = 19 are plotted, and individuals are sorted according to values of the Ket-Uralic component. Admixture coefficients for the Saqqaq ancient genome are shown separately on the right, and for two reference Kets and two Ket individuals from this study - on the left. (B) Average cross-validation (CV) error graph with standard deviations plotted. Ten-fold cross-validation was performed. The graph has a minimum at K = 20. (C) Color-coded values of the Ket-Uralic admixture component at K = 19 plotted on the world map using QGIS v.2.8. Maximum values in each population are taken, and only values >5% are plotted. Top five values of the component are shown in the bottom left corner, and the value for Saqqaq is shown on the map.
Figure 1 (A) Admixture coefficients plotted for dataset ‘GenoChip + Illumina arrays’. Abbreviated names of admixture components are shown on the left as follows: SAM, South American; NAM, North American; ESK, Eskimo (Beringian); SEA, South-East Asian; SIB, Siberian; NEU, North European; ME, Middle Eastern; CAU, Caucasian; SAS, South Asian; OCE, Oceanian; AFR, African. The Ket-Uralic (‘Ket’) admixture component appears at K ≥ 11, and admixture coefficients are plotted for K = 4, 10, 11, and 19. Although K = 20 demonstrates the lowest average cross-validation error, the Ket-Uralic component splits in two at this K value, therefore K = 19 was chosen for the final analysis. Only populations containing at least one individual with >5% of the Ket-Uralic component at K = 19 are plotted, and individuals are sorted according to values of the Ket-Uralic component. Admixture coefficients for the Saqqaq ancient genome are shown separately on the right, and for two reference Kets and two Ket individuals from this study - on the left. (B) Average cross-validation (CV) error graph with standard deviations plotted. Ten-fold cross-validation was performed. The graph has a minimum at K = 20. (C) Color-coded values of the Ket-Uralic admixture component at K = 19 plotted on the world map using QGIS v.2.8. Maximum values in each population are taken, and only values >5% are plotted. Top five values of the component are shown in the bottom left corner, and the value for Saqqaq is shown on the map.CC BY · PubMed Central
Figure 2 (A) A maximum likelihood tree with 6 migration edges computed on the dataset ‘Ket genomes + HumanOrigins’ with selected populations (194,750 SNPs, 39 populations, 527 individuals). Drift parameter is shown on the x-axis. (B) Residuals from the fit of the model to the data visualized. 98% of variance is explained by the tree.
Figure 2 (A) A maximum likelihood tree with 6 migration edges computed on the dataset ‘Ket genomes + HumanOrigins’ with selected populations (194,750 SNPs, 39 populations, 527 individuals). Drift parameter is shown on the x-axis. (B) Residuals from the fit of the model to the data visualized. 98% of variance is explained by the tree.CC BY · PubMed Central
Figure 3 (A) A maximum likelihood tree with 7 migration edges computed on the genome-based dataset without transitions. Edge weight and bootstrap support values are shown in the table, the drift parameter is shown on the x-axis, and bootstrap support values for tree nodes are indicated. Migration edges are numbered according to the order of their appearance in the sequence of trees from m = 0 to m = 8. Note to the figure: as migration edges and tree topology are inter-dependent in bootstrapped trees, bootstrap support for the edges in the original tree was calculated by summing up support for closely similar edges in bootstrapped trees. Below these edge groups are listed for edges #1–7: 1/ Australian and/or Papuan the (Nivkh, Han, Dai, Kinh) clade or any of its members; 2/ Greenlander Inuit or the (Greenlander, Aleutian) clade Saqqaq and/or Late Dorset (optionally a wider clade with Nivkh); 3/ any clade containing African populations any clade composed of Nivkh/Han/Dai/Kinh (optionally a wider clade with Late Dorset and/or Saqqaq and/or Iron Age Altai); 4/ any clade composed of Mal’ta/Afanasievo/Andronovo (optionally a wider clade with Aleutian and/or Mari) Karasuk; 5/ Mal’ta (optionally a wider clade with Motala12/Afanasievo/Andronovo/Aleutian) any clade composed exclusively of Native Americans and/or Greenlander; 6/ any clade composed exclusively of populations with European ancestry Aleutian; 7/ Ket (optionally a wider clade with Karasuk and/or Iron Age Altai and/or Iron Age Russia) Saqqaq and/or Late Dorset. (B) Residuals from the fit of the model to the data visualized. 96.72% variance is explained by the tree.
Figure 3 (A) A maximum likelihood tree with 7 migration edges computed on the genome-based dataset without transitions. Edge weight and bootstrap support values are shown in the table, the drift parameter is shown on the x-axis, and bootstrap support values for tree nodes are indicated. Migration edges are numbered according to the order of their appearance in the sequence of trees from m = 0 to m = 8. Note to the figure: as migration edges and tree topology are inter-dependent in bootstrapped trees, bootstrap support for the edges in the original tree was calculated by summing up support for closely similar edges in bootstrapped trees. Below these edge groups are listed for edges #1–7: 1/ Australian and/or Papuan the (Nivkh, Han, Dai, Kinh) clade or any of its members; 2/ Greenlander Inuit or the (Greenlander, Aleutian) clade Saqqaq and/or Late Dorset (optionally a wider clade with Nivkh); 3/ any clade containing African populations any clade composed of Nivkh/Han/Dai/Kinh (optionally a wider clade with Late Dorset and/or Saqqaq and/or Iron Age Altai); 4/ any clade composed of Mal’ta/Afanasievo/Andronovo (optionally a wider clade with Aleutian and/or Mari) Karasuk; 5/ Mal’ta (optionally a wider clade with Motala12/Afanasievo/Andronovo/Aleutian) any clade composed exclusively of Native Americans and/or Greenlander; 6/ any clade composed exclusively of populations with European ancestry Aleutian; 7/ Ket (optionally a wider clade with Karasuk and/or Iron Age Altai and/or Iron Age Russia) Saqqaq and/or Late Dorset. (B) Residuals from the fit of the model to the data visualized. 96.72% variance is explained by the tree.CC BY · PubMed Central
Figure 4 (A) PC3 vs. PC4 plot for the dataset ‘Ket genomes + HumanOrigins array’. African populations are not shown. Populations are color-coded by geographic region or language affiliation (in the case of Siberian and Central Asian populations), and most relevant populations are differentiated by marker shapes. Ancient genomes are shown in black. For the corresponding PC1 vs. PC2 plot see Suppl. Fig. 6.7. (B) PC3 vs. PC4 plot, zoom on the Ket individuals. Here is a list of populations closest to Saqqaq based on the average Euclidean distances in the multi-dimensional space of ten principal components (distances in parentheses): Ket (0.022), Nganasan (0.025), Selkup (0.026), Yukaghir (0.028), Eskimo (0.032), Koryak (0.032), Mansi (0.032), Itelmen (0.033), Chukchi (0.033), Dolgan (0.035).
Figure 4 (A) PC3 vs. PC4 plot for the dataset ‘Ket genomes + HumanOrigins array’. African populations are not shown. Populations are color-coded by geographic region or language affiliation (in the case of Siberian and Central Asian populations), and most relevant populations are differentiated by marker shapes. Ancient genomes are shown in black. For the corresponding PC1 vs. PC2 plot see Suppl. Fig. 6.7. (B) PC3 vs. PC4 plot, zoom on the Ket individuals. Here is a list of populations closest to Saqqaq based on the average Euclidean distances in the multi-dimensional space of ten principal components (distances in parentheses): Ket (0.022), Nganasan (0.025), Selkup (0.026), Yukaghir (0.028), Eskimo (0.032), Koryak (0.032), Mansi (0.032), Itelmen (0.033), Chukchi (0.033), Dolgan (0.035).CC BY · PubMed Central
Figure 5 Statistics f4 (Mal'ta, Yoruba; Y, X) (A), f4 (Ket, Yoruba; Y, X) (B) computed on the genome-based dataset with African, Australian and Papuan populations excluded. See the corresponding results for the dataset without transitions in Suppl. Figs 8.15 and 8.16, respectively. A matrix of color-coded Z-scores is shown, and ancient genomes are marked with asterisks. Z-score equals the number of standard errors by which the statistic differs from zero, and |Z| > 2.9 demonstrates that the statistic is significantly different from zero using Bonferroni correction for 27 independent tests (threshold p-value of 0.001852). Rows show Z-scores for f4 (Mal'ta, Yoruba; row, column) or f4 (Ket, Yoruba; row, column), vice versa for columns.
Figure 5 Statistics f4 (Mal'ta, Yoruba; Y, X) (A), f4 (Ket, Yoruba; Y, X) (B) computed on the genome-based dataset with African, Australian and Papuan populations excluded. See the corresponding results for the dataset without transitions in Suppl. Figs 8.15 and 8.16, respectively. A matrix of color-coded Z-scores is shown, and ancient genomes are marked with asterisks. Z-score equals the number of standard errors by which the statistic differs from zero, and |Z| > 2.9 demonstrates that the statistic is significantly different from zero using Bonferroni correction for 27 independent tests (threshold p-value of 0.001852). Rows show Z-scores for f4 (Mal'ta, Yoruba; row, column) or f4 (Ket, Yoruba; row, column), vice versa for columns.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 18 references

Research works printed but not linked: 37

  1. Vajda, E. J. Loanwords in Ket. The Typology of Loanwords, ed. Haspelmath, M., Tadmoor, U. Oxford: Oxford University Press, 125–139 (2009). · DOI
  2. Vajda, E. J. Yeniseian Peoples and Languages: a History of Their Study with an Annotated Bibliography and a Source Guide. Surrey, England: Curzon Press, 389 p. (2001).
  3. Dul’zon, A. P. Ketskie toponimy Zapadnoy Sibiri [Ket toponyms of Western Siberia]. Uchenye Zapisky Tomskogo Gosudarstvennogo Pedagogicheskogo Instituta [Scholarly Proceedings of Tomsk State Pedagogical Institute] 18, 91–111 (1959).
  4. Vajda, E. J. Siberian landscapes in Ket traditional culture. Landscape and Culture in Northern Eurasia, ed. Jordan, P. Walnut Creek, CA: Left Coast Press, 297–304 (2011).
  5. Vajda, E. J. Siberian link with Na-Dene languages. The Dene-Yeniseian Connection, ed. Kari, J., Potter, B. A. Anthropological Papers of the University of Alaska: New Series 5, 33–99 (2010).
  6. Comrie, B. The Dene-Yeniseian hypothesis: an introduction. The Dene-Yeniseian Connection, ed. Kari, J., Potter, B. A. Anthropological Papers of the University of Alaska: New Series 5, 25–32 (2010).
  7. Ruhlen, M. The origin of the Na-Dene. Proc. Natl. Acad. Sci. USA 95, 13994–13996 (1998). · DOI
  8. Potter, B. A. Archaeological patterning in Northeast Asia and Northwest North America: an examination of the Dene-Yeniseian hypothesis. The Dene-Yeniseian Connection, ed. Kari, J., Potter, B. A. Anthropological Papers of the University of Alaska: New Series 5, 138–167 (2010).
  9. Scott, R. G. & O’Rourke, D. Genes across Beringia: a physical anthropological perspective on the Dene-Yeniseian hypothesis. The Dene-Yeniseian Connection, ed. Kari, J., Potter, B. A. Anthropological Papers of the University of Alaska: New Series 5, 119–137 (2010).
  10. Ives, J. W. Dene-Yeniseian, migration and prehistory. The Dene-Yeniseian Connection, ed. Kari, J., Potter, B. A. Anthropological Papers of the University of Alaska: New Series 5, 324–334 (2010).
  11. Fedorova, S. a. et al. Autosomal and uniparental portraits of the native populations of Sakha (Yakutia): implications for the peopling of Northeast Eurasia. BMC Evol. Biol. 13, 127 (2013). · DOI
  12. Yunusbayev, B. et al. The genetic legacy of the expansion of Turkic-speaking nomads across Eurasia. PLoS Genet. 11, e1005068 (2015). · DOI · Paper page
  13. Elhaik, E. et al. The GenoChip: a new tool for genetic anthropology. Genome Biol. Evol. 5, 1021–1031 (2013). · DOI
  14. Raghavan, M. et al. The genetic prehistory of the New World Arctic. Science 345, 1255832 (2014).
  15. McGhee, R. Ancient People of the Arctic. Vancouver: UBC Press (1996).
  16. Reich, D. et al. Reconstructing Native American population history. Nature 488, 370–374 (2012). · DOI · Paper page
  17. Tamm, E. et al. Beringian standstill and spread of Native American founders. PLoS One 2, e829 (2007). · DOI · Paper page
  18. Gilbert, M. T. P. et al. Paleo-Eskimo mtDNA genome reveals matrilineal discontinuity in Greenland. Science 320, 1787–1789 (2008). · DOI
  19. Malhi, R. S. et al. Native American mtDNA prehistory in the American Southwest. Am. J. Phys. Anthropol. 120, 108–124 (2003). · DOI
  20. Sicoli, M. A. & Holton, G. Linguistic phylogenies support back-migration from Beringia to Asia. PLoS One 9, e91722 (2014). · DOI
  21. Pringle, H. Welcome to Beringia. Science 343, 961–963 (2014). · DOI
  22. Elhaik, E. et al. Geographic population structure analysis of worldwide human populations infers their biogeographical origins. Nat. Commun. 5, 3513 (2014). · DOI
  23. Kozlov, K. et al. Differential Evolution approach to detect recent admixture. BMC Genomics 16 Suppl 8, S9 (2015). · DOI
  24. Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009). · DOI · Paper page
  25. Quinque, D., Kittler, R., Kayser, M., Stoneking, M. & Nasidze, I. Evaluation of saliva as a source of human DNA for population and association studies. Anal. Biochem. 353, 272–277 (2006). · DOI
  26. Abecasis, G. R. et al. An integrated map of genetic variation from 1,092 human genomes. Nature 491, 56–65 (2012). · DOI · Paper page
  27. Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007). · DOI · Paper page
  28. Behar, D. M. et al. The genome-wide structure of the Jewish people. Nature 466, 238–242 (2010). · DOI
  29. Cardona, A. et al. Genome-wide analysis of cold adaptation in indigenous Siberian populations. PLoS One 9, e98076 (2014). · DOI
  30. Li, J. Z. et al. Worldwide human relationships inferred from genome-wide patterns of variation. Science 319, 1100–1104 (2008). · DOI
  31. Kidd, J. R. et al. Single nucleotide polymorphisms and haplotypes in Native American populations. Am. J. Phys. Anthropol. 146, 495–502 (2011). · DOI
  32. Silva-Zolezzi, I. et al. Analysis of genomic diversity in Mexican Mestizo populations to develop genomic medicine in Mexico. Proc. Natl. Acad. Sci. 106, 8611–8616 (2009). · DOI · Paper page
  33. Surakka, I. et al. Founder population-specific HapMap panel increases power in GWA studies through improved imputation accuracy and CNV tagging. Genome Res. 20, 1344–1351 (2010). · DOI
  34. Yunusbayev, B. et al. The Caucasus as an Asymmetric Semipermeable Barrier to Ancient Human Migrations. Mol. Biol. Evol. 29, 359–365 (2011). · DOI
  35. Behar, D. M. et al. No Evidence from Genome-Wide Data of a Khazar Origin for the Ashkenazi Jews. Hum. Biol. 85, 859–900 (2013). · DOI
  36. Verdu, P. et al. Patterns of Admixture and Population Structure in Native Populations of Northwest North America. PLoS Genet. 10, e1004530 (2014). · DOI
  37. Axelsson, E., Willerslev, E., Gilbert, M. T. P. & Nielsen, R. The Effect of Ancient DNA Damage on Inferences of Demographic Histories. Mol. Biol. Evol. 25, 2181–2187 (2008). · DOI · Paper page

Cited by

See all 57

Metadata

Authors
Flegontov, P., Changmai, P., Zidkova, A., Logacheva, M. D., Altınışık, N. E., Flegontova, O., Gelfand, M. S., Gerasimov, E. S., Khrameeva, E. E., Konovalova, O. P., Neretina, T., Nikolsky, Y. V., Starostin, G., Stepanova, V. V., Travinsky, I. V., Tříska, M., Tříska, P., & Tatarinova, T. V.
Venue
Scientific Reports
Published
2016-02-11
Volume
6
Issue
1
ISSN
2045-2322
Cited by
57
References
18
Language
en

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

How to cite

Cite CoChat for the checks on this page, and the paper for its own findings.

This page
CoChat. (2026). Citation check: Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry. https://app.cochat.ai/papers/10.1038/srep20768
The paper
Flegontov, P., Changmai, P., Zidkova, A., Logacheva, M. D., Altınışık, N. E., Flegontova, O., Gelfand, M. S., Gerasimov, E. S., Khrameeva, E. E., Konovalova, O. P., Neretina, T., Nikolsky, Y. V., Starostin, G., Stepanova, V. V., Travinsky, I. V., Tříska, M., Tříska, P., & Tatarinova, T. V. (2016). Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry. Scientific Reports. https://doi.org/10.1038/srep20768