- Open Access
Between Lake Baikal and the Baltic Sea: genomic history of the gateway to Europe
BMC Genetics volume 18, Article number: 110 (2017)
The history of human populations occupying the plains and mountain ridges separating Europe from Asia has been eventful, as these natural obstacles were crossed westward by multiple waves of Turkic and Uralic-speaking migrants as well as eastward by Europeans. Unfortunately, the material records of history of this region are not dense enough to reconstruct details of population history. These considerations stimulate growing interest to obtain a genetic picture of the demographic history of migrations and admixture in Northern Eurasia.
We genotyped and analyzed 1076 individuals from 30 populations with geographical coverage spanning from Baltic Sea to Baikal Lake. Our dense sampling allowed us to describe in detail the population structure, provide insight into genomic history of numerous European and Asian populations, and significantly increase quantity of genetic data available for modern populations in region of North Eurasia. Our study doubles the amount of genome-wide profiles available for this region.
We detected unusually high amount of shared identical-by-descent (IBD) genomic segments between several Siberian populations, such as Khanty and Ket, providing evidence of genetic relatedness across vast geographic distances and between speakers of different language families. Additionally, we observed excessive IBD sharing between Khanty and Bashkir, a group of Turkic speakers from Southern Urals region. While adding some weight to the “Finno-Ugric” origin of Bashkir, our studies highlighted that the Bashkir genepool lacks the main “core”, being a multi-layered amalgamation of Turkic, Ugric, Finnish and Indo-European contributions, which points at intricacy of genetic interface between Turkic and Uralic populations. Comparison of the genetic structure of Siberian ethnicities and the geography of the region they inhabit point at existence of the “Great Siberian Vortex” directing genetic exchanges in populations across the Siberian part of Asia.
Slavic speakers of Eastern Europe are, in general, very similar in their genetic composition. Ukrainians, Belarusians and Russians have almost identical proportions of Caucasus and Northern European components and have virtually no Asian influence. We capitalized on wide geographic span of our sampling to address intriguing question about the place of origin of Russian Starovers, an enigmatic Eastern Orthodox Old Believers religious group relocated to Siberia in seventeenth century. A comparative reAdmix analysis, complemented by IBD sharing, placed their roots in the region of the Northern European Plain, occupied by North Russians and Finno-Ugric Komi and Karelian people. Russians from Novosibirsk and Russian Starover exhibit ancestral proportions close to that of European Eastern Slavs, however, they also include between five to 10 % of Central Siberian ancestry, not present at this level in their European counterparts.
Our project has patched the hole in the genetic map of Eurasia: we demonstrated complexity of genetic structure of Northern Eurasians, existence of East-West and North-South genetic gradients, and assessed different inputs of ancient populations into modern populations.
The phenotypic diversity of modern humans was shaped under the combined pressure of environment and social relations. Placing the studies of human genetic variation into a geographical context provides powerful insights into how historical events, patterns of migration, and natural selection have led to genetic distinctions between various present-day populations [1, 2]. Moreover, genomic investigations may aid in resolving historic record discrepancies by confirming or rejecting hypotheses of ancient invasions and ethnic intermixing events.
While human genetic diversity has been sampled extensively in many areas of the globe [3,4,5,6], a sizeable gap remains in the region of Northern Eurasia (region including Russia and neighboring countries from the former Soviet Union) which spans from the Arctic Ocean down to Inner Asia, and from Eastern Europe to the Pacific Ocean. Though in total, human populations inhabiting this region were analyzed among others in several genome-wide studies [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22] most of them were focused on other regions and included just a limited number of Northern Eurasian populations. Only five published studies were focused on the areas within North Eurasia: two papers of Yunusbayev [15, 19] investigated genetic composition of the Caucasus and Turkic speaking groups;  focused on Balto-Slavic speakers;  and  studies were even more limited in their geographic coverage. Thus, a panoramic genetic study covering all of Northern Eurasia is still lacking. The Russian Federation represents a unique setting for genetic studies because of its multitude of ethnicities with the evidence for admixture interspersed across several isolated communities. Further, its enormous space and considerable climatic variation created a range of distinct environmental niches which may have contributed to differential shaping of the genomes. However, the limited number of sampled populations in the published datasets translates into significantly less coverage incomparable to that for the western and central regions of Europe.
Here we present high-quality genome-wide analysis of 30 diverse populations from Russia and neighbouring countries (see Table 1). Some of these populations have been previously studied on a smaller scale, and some been sampled here for the first time. Though full genome sequences on population level started to accumulate worldwide extensively, only 246 full genomes were so far published for the Russian populations [10, 11, 16, 17, 23,24,25,26]. In contrast, genome-wide genotypes were published for 963 samples from 51 Russian populations [10, 16]. So, the genotyping arrays remain the most important source of genomic variation within Northern Eurasia. Here we double this aggregate dataset by publishing genome-wide genotype data on 1076 samples (1019 of them unrelated) from 30 populations of Russia and adjacent countries.
Certain unusually diverse areas were given special consideration, such as the Caucasus, where all the major ethnicities, including Abkhaz, Adygei, Chechen, Cherkes, Kabardian, Karachay, Megrel, and Ossetian were profiled. We have also sampled several unique populations, such as the Ket - an isolated, native Siberian people with a distinct language  and the Starover Russians, orthodox Old Believers who left western Russia in the seventeenth century and settled in the dense boreal forests of the banks of Volga and the Russian European North, as well as on the southern outskirts of Siberia . Starovers maintain the liturgical and ritual practices of the Russian Orthodox Church as they existed prior to the reforms of Patriarch Nikon of Moscow between 1652 and 1666. In this work, we studied the descendants of Siberian Starovers, who presumably had limited admixture with other groups.
Several independent groups of researchers [1, 2, 8, 28,29,30,31,32,33] have analyzed the relation between genetic variation and geography, with a variety of biogeographical analysis techniques developed [2, 8, 32,33,34,35,36,37,38,39]. This relationship was extensively studied for European populations [32, 33, 35], for Indian casts [40,41,42], and, more generally, for world-wide populations . Here we present a detailed analysis of Northern Eurasian populations inhabiting the territories of Russian Federation, and neighbouring countries.
Sample collection and quality controls
DNA samples (N = 1076) were collected in course of study expeditions into different parts of Russia, Kazakhstan, Georgia, Uzbekistan and Kyrgyzstan. Samples were genotyped on the Illumina Infinium 370-Duo, 370-Quad, or 610-Quad arrays (https://support.illumina.com/downloads/humancnv370-duo_v10_product_files.html, https://support.illumina.com/downloads/humancnv370-quad_v30_product_files.html, https://support.illumina.com/array/array_kits/human610-quad_beadchip_kit.html). The number of samples per population, source of the samples, and the type of microarray used are given in the Table 1.
All DNA samples were subjected to the following quality control procedures: samples with genotyping success rates <90% were removed, as were male samples with ≥1% heterozygous markers on the X chromosome or female samples with ≤20% heterozygous X chromosome markers. Across retained samples, 95% cut-off for SNP presence was imposed. (Further details are provided in the Additional file 1: Figure S1).
According to the meta-data from questionnaires, all study volunteers were unrelated to each other. Nevertheless, the dataset was analysed for the presence of cryptic relatedness by calculating the kinship coefficients separately in each ethnic group using the King software , assuming the presence of population structure. For 29 related pairs of individuals with a threshold of kinship coefficient set at ≥0.177 , the sample with the lower genotyping-call rate was excluded from further analysis, thus, only 1019 samples remained out of 1076.
To reduce the effect of missing data, the marker panel was limited to autosomal SNPs with genotyping success rates ≥99.5%. For each set of samples, more than 200,000 markers were analysed.
Origin of the samples
To provide a uniform representation of ethnic diversity in Russia we sampled broadly over the country (see Fig. 1 and Table 1) and adjacent territories. In total, genome-wide variation was accessed in 1019 individuals from five countries across the former Soviet Union. All except one of the studied populations (Chinese) were covered by at least 20 samples (see Table 1). To limit influence of recent admixture, we ensured sampling of villagers who reportedly were settled in the sample place for at least three generations (up to grandparents).
The collected data were assembled in three datasets.
The “Extended” dataset includes all individuals genotyped in this study combined with selected previously published modern and ancient samples from Northern Eurasia, which extend geographic span of our study and provide necessary populational context for our analyses (1353 individuals from 55 populations, plus 11 ancient samples, shown in Additional file 3: Table S1a). The “Extended” dataset was used for ADMIXTURE analysis.
The “Core” dataset contains samples genotyped in this study (1019 individuals shown in Additional file 3: Table S1b). The “Core” dataset was used to calculate IBD sharing and f3 statistics.
The “Ancient” dataset includes all individuals genotyped in this study combined with European samples from “1000 Genomes” project as well as previously published ancient DNA samples (1232 individuals shown in Additional file 3: Table S1c). This dataset was used to calculate the f3 outgroup statistics.
ADMIXTURE [44, 45] algorithm was used in unsupervised mode to determine the population structure. The number of components (K) was varied from 2 to 10, and cross-validation errors was recorded for all values of K.
For GPS  and reAdmix  analyses, the reference dataset was obtained from E Elhaik, T Tatarinova, D Chebotarev, IS Piras, C Maria Calò, A De Montis, M Atzori, M Marini, S Tofanelli, P Francalacci, et al. . To enable the comparison with earlier published results, SNPs were converted to the 9-dimensional admixture vectors (“North East Asian”, “Mediterranean”, “South African”, “South West Asian”, “Native American”, Oceanian”, “South East Asian”, “Northern European”, “Sub-Saharan African”) using the ADMIXTURE [44, 45] algorithm in supervised mode.
IBD blocks were identified for every pair of individuals in the “Core” dataset using fastIBD algorithm implemented in BEAGLE 3.3 . For further analysis, the amount of shared IBD was calculated separately in two bins: short 1–3 cM blocks and longer 4–10 cM blocks. Total amount of shared IBD was then averaged per pair of individuals.
The three-population statistics were calculated in threepop software included in Treemix  with the block size parameter (−k) set to 500 SNPs (other parameters set to default values).
GPS and reAdmix provenance identification
To validate compatibility of the genotyping technologies based on different platforms, we applied GPS  and reAdmix  to the combined dataset consisting of populations sampled in our study along with ones previously sampled by the National Genographic project, HapMap, and “1000 genomes”. Depending on the database coverage of a selected region and ethnic group, GPS/reAdmix accuracies differ. For example, in this analysis, self-identifications of all profiled Kyrgyz, 75% of Kabardian, and 60% self-reported Ket were correctly identified by GPS. For populations not covered by reference database, locations were triangulated using neighboring ethnic groups.
In our study, collection points spanned the latitudes from the Northern Europe to Caucasus and longitudes from Eastern Europe to far shores of Siberia. The results of unsupervised analysis of the “extended” dataset with ADMIXTURE [44, 45] varying the number of ancestral clusters (K) from 2 to 10 (Additional file 2: Figure S2) were most informative for K = 6 (Fig.2), which resulted in lowest cross-validation error and meaningful distribution of components between studied populations.
To enable direct comparison with previous world-wide studies by National Genographic Project [7, 8, 10, 31, 39, 48], ADMIXTURE analysis in supervised mode were performed for K = 9. For each chip type, we selected a subset of SNPs (~30 K and ~60 K for the two chip types) that matched National Genographic Project chip Geno 2.0 .
Among the profiled populations, the degrees of ancient admixture varied dramatically. While the populations residing in the Caucasus and East Siberia regions were mostly represented by a single component (dark green and dark blue, correspondingly), the samples from the Volga-Ural region exhibit substantial admixture of European and Asian components (red, dark green, and light blue). Importantly, for each of the populations, the ratio between components was characteristic of that population (Additional files 3: Tables S2 and S3). Moreover, we observed a subdivision of the European genetic component into two clusters, one most prevalent in the Caucasus (dark green) and another with highest frequency in Northern Europe (red). These two European components jointly account for 50% - 90% of admixture vectors in both Turkic and Uralic speakers of Volga-Ural region, while in Finno-Ugric speakers in Northeast Europe and in all Slavic populations these components account for almost 90% of the gene pool.
Slavic speakers – Russians, Ukrainians, and Belarussians - are similar in their genetic composition. Ukrainians and Belarusians have almost identical proportions of the two “European” components and have virtually no “Asian” admixture. Russians from Siberia - both Novosibirsk residents’ and Russian Eastern Orthodox Old Believers’ (Starover) samples - being genetically close to Slavs residing in Europe, also have between five to 10 % of Central Siberian ancestry (light blue).
Asian genetic ancestry of the profiled populations is represented by four components. Two of them, light green (Beringian) and dark blue (East Siberian), are geographically confined to Northeast Asia, while showing only minor impact on populations west of Ural Mountains. On the contrary, East Asian (pink) and Central Siberian (light blue) components are also present in populations to the west of Ural.
In Central Asian Turkic speakers, including Kazakh and Uzbek, East Asian genetic influence is dominant (>35%), while in Bashkir it is detected at somewhat lower levels (~ 20%). Importantly, in Western Turkic speakers, like Chuvash and Volga Tatar, the East Asian component was detected only in low amounts (~ 5%).
The light blue genetic component dominates genetic landscape of populations inhabiting West and Central Siberia: Ugric-speaking Khanty and Mansi, Samoyedic speaking Selkups and linguistically isolated Ket. However, this ancestry component is present not only in Siberia, but also on the western side of Ural Mountains, though at somewhat lower frequencies - 20-30% in Komi (16% on average) and Udmurt (27% on average) who belong to the Permic branch of Uralic languages. Interestingly, similar levels of this ancestry component (16–23%) are also exhibited by Turkic speaking Chuvash (20% on average) and Bashkir (17% on average), while Tatar, who also reside in the Volga region and have related linguistic and cultural profiles, only show at most 15% (10% on average) of this genetic component. Even lower levels of this ancestry component (<5%) were observed in Turkic speakers of Central Asia.
The Beringian component (light green) is confined exclusive to indigenous populations of Eskimo, Chukchi and Koryak. The East Siberian (dark blue) component is represented by Turkic and Samoyedic speakers of Central Siberian plateau: Yakut, Dolgan and Nganasan. This component is also found at moderate frequencies in Mongolic and Turkic speakers in Baikal region and Central Asia (5–15%), and, at low but discernible frequencies (1–5%), in Turkic speakers residing in Volga-Ural region.
f3 population test
ADMIXTURE-guided ancestry clustering suggests that, in most of the populations studied, the genetic background is complex, as it includes at least three hypothetical ancestral components: European, East Asian and North Asian (Siberian). Because the ADMIXTURE-guided ancestry clustering cannot be used in lieu of the formal test of admixture, three-population test  was conducted to confirm proposed admixture events. Surrogate populations were selected for each of three ancestral components, followed by f3 test of admixture, to find out whether a target population is admixed between two source populations. The combinations of source and target populations and the Z score in the f3 test are reported in Additional file 3: Table S4.
Admixture scenarios were tested in sampled populations geographically grouped as Northern Europe, Volga-Ural region and Central Asia and, while assuming two-way mixture between European and East Asian, or between European and Central Siberian surrogate populations. Significant Z scores (Z < −3) for admixture between East European (Belarus, Ukrainian) and either East Asian (Tuva, Yakut) or Central Siberian (Khanty, Ket) populations were obtained for all populations of Northern Europe and Volga-Ural region.
Note, that the value of f3(X; A, B) (and the corresponding Z score) is negative if X is a mixture of A and B. Among all pairs compared, the most negative f3 values were obtained for populations of Volga-Ural region using East Slavs and Siberians as source groups (Additional file 3: Table S4). This suggests that both, Western (including ancestors of East Slavs) and Eastern (Siberian) sources of the formation of Volga-Ural populations, which could be also seen from the ADMIXTURE plot (Fig. 2). Although Ukrainians and Belarus received similar scores in the f3 test, Belarus turned out to be a slightly better proxy for East European component in populations of Volga-Ural region.
Populations from Central Asia also received significant negative Z scores in the f3 test using East Asia and Europe populations as a source, although their best supported surrogate populations were different. For Uzbek and Kazakh, the best surrogate for a European source was Abkhaz, while for Volga-Ural populations it was either Belarus or Ukrainian. For Uzbek, the best East Asian surrogate was Yakut, while for Kazakh and Kyrgyz it was Tuva (followed by Yakut).
Analysis of IBD blocks
To identify shared identical-by-descent (IBD) blocks, we used the fastIBD algorithm implemented in the Beagle package . We first calculated total amount of IBD in centimorgans (cM) shared between the populations, which was then averaged per pair of individuals (Additional file 3: Table S6). Since the length of IBD blocks is anti-correlated with their respective age, analysing the distribution of length of blocks allows us to examine patterns in ancestry sharing with temporal resolution . We analysed the amount of shared IBD in two length bins: 1–3 cM (ancient blocks) and 4–10 cM (recent blocks) (Fig. 3 and Additional file 3: Table S5). We focused on three regions, providing the densest sample coverage: Caucasus, Volga basin and Siberia. Populations from the Caucasus share most of IBD blocks between themselves and the amount of shared IBD ranges between 3.26 cM to 12.39 cM per pair, which is roughly comparable to the amount of IBD blocks shared between Eastern European populations in our dataset (Additional file 3: Table S5). A conspicuous exception are the Chechens, who share almost all detected IBD blocks among themselves and only a scant amount of IBD with neighbouring Caucasus populations. This may have happened because the Chechen sample is the only representative of the North-East Caucasus in our dataset. Low amount of IBD shared outside the cluster of Caucasian populations suggests a lack of recent ancestry links with Uralic, Slavic or Turkic people (except for Turkic from the Caucasus) present in our dataset.
The Volga-Ural region is populated by three major language and cultural groups: Uralic, Turkic and Slavic speakers. Bashkir and Tatar are major Turkic groups in the region. Although both ethnic groups live in the same region and their languages are mutually intelligible, we surprisingly detected only a limited amount of ancient IBD blocks shared between them, and their overall IBD sharing pattern is different: Tatar share moderate amount of IBD (3.55–7.35 cM per pair) with all neighbouring populations, while Bashkir share most of their ancient blocks (on average 8.62 cM per pair) with Khanty, a group of Uralic speakers from Western Siberia. We speculate that this disparity between cultural and genetic affinities of Tatar and Bashkir can be attributed to a phenomenon of cultural dominance: the population ancestral to Bashkir adopted the Turkic language during Turkic expansion from the east (language replacement event).
European and Siberian Uralic speakers are separated by the Ural Mountain range. This separation affects sharing of recent IBD blocks: Komi and Udmurt share in the interval of 4–10 cM almost double the amount of IBD they share with Khanty. Interestingly, the situation is different when we look at sharing of ancient IBD blocks: the amount of ancient IBD blocks that Udmurt and Komi share with Khanty (10.63 and 8.62 cM per pair, respectively) is comparable to the amount of ancient IBD shared between them (11.38 cM per pair). These data agree well with the ethnic history of European and Asian Uralic speakers . Udmurt and Komi belong to the Permic branch of the Uralic language family and share a recent origin, as demonstrated by common short IBD blocks, their split from their ancestral Uralic population is dated back to the first half of the 1st millennium BC. Split between ancestors of Asian Uralic people, represented in this paper by Khanty, and ancestors of modern European Uralic ethnic groups (Udmurt and Komi) is much older and dated back approximately to 3rd millennium BC. All modern Uralic populations share common genetic substrate inherited from some ancient Uralic people, reflected in long and similar size ancient IBD blocks shared between Udmurt, Komi and Khanty. All analysed native Siberian populations exhibit high levels of intrapopulation sharing of IBD (Fig 3), which is in line with observed long runs of homozygosity in these populations. High rates of shared IBD blocks were detected also between pairs of Siberian populations, particularly between Ket and Khanty.
Following B Yunusbayev, M Metspalu, E Metspalu, A Valeev, S Litvinov, R Valiev, V Akhmetova, E Balanovska, O Balanovsky, S Turdikulova, et al. , who observed that the number of shared IBD blocks decline exponentially with the distance between populations, we have calculated linear regression between geographic distance and logarithm of IBD for all pairs of populations from the “Core” dataset (see Additional file 3: Tables S7 and S8). This resulted in the following equation: log 10 IBD = 1.772 − 0.0005975 × (Distance in kilometers), adjusted R 2 = 0.4923, p-value <10−16.
Although Russian Starovers appeared in Siberia only about 300 years ago, details of their geographic origin are not clear. Starovers share roughly twice as much of ancient blocks with Komi and Karelian than between themselves, or with other Slavic speaking groups (Additional file 3: Table S5). This finding, corroborated by results of GPS and reAdmix (described below), strongly points to Northern European ancestral ties of Russian Starovers. The most pronounced difference between predicted and observed IBD rate was found for Siberian Russians, the residents of Novosibirsk and the Starovers, and for the Chinese outgroup (Additional file 3: Tables S7 and S8, Fig. 4). This was expected due to relatively recent relocation of Eastern European Slavic individuals to Siberia. Interestingly, one of the Caucasian ethnicities, the Chechens, also share fewer IBD blocks with Chinese than would be expected from the geographical distance separating these populations. This outlier could be possibly explained by underestimation of the degree of apparent geographical isolation of Chechens, occupying hard to reach highlands of Caucasus range, since elevation is not considered by the regression. From patterns of IBD-sharing and from ADMIXTURE-based analysis we see that Caucasus populations differ from their neighbors in the European part of the Russian Federation. One of the factors being that some Caucasus tribes reside at high altitudes of 2000 m above sea level or more, where they have been genetically isolated for centuries [51,52,53,54].
At the next stage of analysis, Novosibirsk Russian and Chinese populations were excluded, leaving us with 325 population pairs. This analysis produced following equation: log 10 IBD = 1.394 − 0.0004121 × (Distance in kilometers), adjusted R 2 = 0.2774, p–value <10−16. Using this equation, we have identified pairs of populations (using the “core” dataset) that were more than two standard deviations away from the predicted log 10 IBD (Fig. 4). In this analysis, the departure from the regression line suggests unusual gene flow events, unaccounted for by the calculation of geographic proximity. Most significant departures were observed for Yakut with Kyrgyz, Kazakh, Buryat and Tuva, as well as for Ket with Khanty and Tuva combinations. Starover, who recently migrated to Siberia, again show genetic proximity to Karelian (Fig. 5).
GPS and reAdmix analyses of self-identified Russians
In our dataset, ethnic Russians constituted the largest population (N = 80). These samples came from two groups: residents of Novosibirsk district (39 samples, designated “NSK”) and the Siberian Starover (41 samples, designated “STV”). The Novosibirsk residents were thoroughly surveyed about their ancestors and selected only if reporting at least three preceding Russian generations, while members of the Starover cohort were all assumed to be “authentic” Russians. Since their resettlement from the European part of Russia in the seventeenth century, Russian Starovers deliberately adhered to a strict religious routine and avoided contact with neighbouring Native Siberian populations. Both the Starovers and most Novosibirsk residents are informally considered as “canonical Russians”. Nevertheless, only the Novosibirsk group represents a uniform sample from the modern Russian gene pool.
For all populations sampled in this study, SNP array data were compared to the worldwide collection of populations using Geno 2.0130 K ancestry-informative markers (AIMs) . Both SNP platforms used in our analysis contain a subset of these markers. Chip 370 includes 60,730 AIMs, while chip 610 includes 90,231 AIMs. As demonstrated earlier , even in case if admixture vectors are determined with as little as 40,000 AIMs, the difference between the admixture vectors obtained from the complete set of AIMs and the reduced set does not exceed 3%. Therefore, reduction of AIMs to 60,000 or more resides within the range of natural variation and does not affect the accuracy of population assignment.
In the following analysis, we used admixture vectors obtained by ADMIXTURE software run with reference dataset from E Elhaik, T Tatarinova, D Chebotarev, IS Piras, C Maria Calò, A De Montis, M Atzori, M Marini, S Tofanelli, P Francalacci, et al.  in supervised mode (К = 9). When admixture vectors for Novosibirsk and Starover Russians were compared, relative weights of the “Northern European” component were found to differ by 2% (t-test p-value <0.009). The provenances of the samples were inferred by two algorithms, GPS  and reAdmix . For each tested individual, GPS algorithm determines a location on a world map, where people with similar genotypes are most likely to reside. Notably, this algorithm is not suitable for analysis of recently mixed individuals, such as children of parents from two different ethnic groups. When subjected to GPS, a recently admixed sample would result in a report of high uncertainty of prediction.
To analyze modern Russians from Novosibirsk and Starover Russians, we used Russian diversity panel data genotyped on Geno 2.0 chip (Balanovsky et al., unpublished data). Nearly 37% of self-identified Russians from Novosibirsk were mapped to various Russian populations from European part of Russia: 13% to Tver region, 13% to Arkhangelsk region, 5% to Ryazan region and 3% to Don Cossacks and Vologda region each. Not surprisingly, as many as 27% of Novosibirsk residents were identified as Mordva: 24% as Erzya and 3% as Moksha (Fig. 6). These two subethnic groups were followed by Chuvash (16%), Karelian and Evenki (5% each). Many singular representations of other ethnic groups of the Russian Federation were also reported.
Starover Russians appear to be more closely related to European Russians than Russians in Novosibirsk, with 58% identified as descendants of the migrants from various cities and villages in European part of Russia: for 45% of them, the provenance was traced to Arkhangelsk region, for 7% to Vologda region, for 2% to Yaroslavl region, and for 2% each to Tver region and Don Cossacks. Other notable ethnic component groups include 23% of Erzya and Moksha Mordva collectively, 12% of Karelian, and 5% of Veps.
Altogether, GPS analysis of Starovers suggests that most of them came from northern areas of European Russia. This agrees with the slightly higher value of the Northern European component in Starovers as compared to Novosibirsk Russians.
In addition to the proposed population and geographic location, the GPS algorithm also reports prediction uncertainty calculated from the distance to the nearest reference population. One of the Starover individuals was identified by GPS as a Khakas, a Turkic ethnicity living in the Republic of Khakassia located in southern Siberia, Russia. The same individual had the largest prediction uncertainty (7%) as compared to the average 3% prediction uncertainty for other Starovers samples. Typically, the prediction uncertainties which exceed 4% indicate mixed origin of an individual. For these cases, GPS algorithm should not be used.
Therefore, for further analysis of Starovers and Novosibirsk individuals, we used reAdmix , which represents each individual as weighted sums of modern reference populations (see Fig. 7). In agreement with the GPS results, self-identified Russians from Novosibirsk appear to be more admixed than the Starovers. In Novosibirsk, 37% of genetic input came from ethnic Russians (15% from Northern Russia and 23% from Southern Russia), 25% Finno-Ugric (Veps, Karelian, Mordva), and 38% to other (Buriat, Chukchi, Chuvash, Dolgan, Evenki, Ket, Nenets, Nganasan, Selkup, Tatar, Tuvinian, Yakut, Yukaghir). Among the Starovers, 50% of ancestry was attributed to Russians (with 41% from Northern Russian and 9% of Southern Russia), 33% to Finno-Ugric (Veps, Karelian, Mordva), and 17% to other, including native Siberian populations (such as Tuva, Buryat, Yakut, Ket, Khanty). This observation supports the notion that Siberian Starovers represent relatively large heterogeneous group, which did not stay entirely isolated.
Since the strict religious rules prevented Russian Starovers from marrying members of other ethnic groups, they are commonly believed to be less admixed with native Siberians than other Russian communities in the region. However, our ADMIXTURE analysis showed that the admixture profiles of Russian Starovers and Russians from Novosibirsk are similar (Fig. 2) and that both groups experienced comparable gene flow.
This genetic input can be attributed to multiple known and unknown events in the history of Starovers. We can summarize our observations as follows. According to both GPS and reAdmix analyses, Starover Russians have more significant input from Northern Russians and Finno-Ugric populations than from the South of Russia. Novosibirsk Russians represent a typical mixed Russian population of the early twenty-first century; lesser degree of admixture in the genomes of Starovers point at an increase in the rates of admixture of Russian populations with neighboring ethnicities that occurred in the last 300–400 years.
f3 outgroup analysis of relatedness to ancient genomes
Earlier comparative studies of ancient and modern human DNA have helped to delineate human migration routes around the world [56,57,58,59,60]. We used f3 outgroup statistics  to test for shared genetic drift between our studied populations and selected ancient populations, namely East European hunter gatherers, Caucasus hunter gatherers, Anatolian farmers and Mal’ta (See Additional file 3: Table S1c). It was demonstrated that f3 is positive if and only if the branch supporting the population tree is longer than the two branches discordant with the population tree . Therefore, large positive values of f3 show that the two tested populations had a large amount of shared population drift.
All populations from the “Extended” dataset were used as test populations, with Mal’ta , Eastern European hunter-gatherers [56, 58], Caucasus hunter-gatherers (Jones et al. 2015) and Neolithic samples (Mathieson et al. 2015) as the reference and Yoruba as an outgroup (Additional file 3: Table S9). The summary of the findings is shown in the Fig. 8.
This analysis confirms local inheritance of genetic structure between ancient and modern populations, as evident from consideration of aDNA samples from the Caucasus and Europe. We did not find the “source” population for our Eastern Siberian samples (Yakut and Buryat). We also confirmed that modern European population is an amalgamation of ancient European Hunter-Gatherers with Neolithic Farmers. [57, 63,64,65]. Neolithic Farmers’ genetic influence is present in a wide range of modern Eurasian populations (from the Iberian Peninsula in the West to the Altay mountains in the East). East of Altay the signal fades. Genetic signal from European Hunter Gatherers is present across several Northern Eurasian populations. The modern populations of the Caucasus show a strong signal from Caucasus hunter gatherers, that is almost absent elsewhere. Ancient North Europeans (represented by Mal’ta boy) left their genetic mark on several genomes of modern Northern Eurasians, without affecting Western or Southern Europeans or Eastern Siberians or Central Asians.
Since the pioneering effort by the HapMap Consortium made in 2003 , multiple studies were conducted to investigate human genetic diversity, population structure, migration routes, and genotype-phenotype association [2, 8, 16, 32, 33, 55, 66,67,68,69,70,71,72,73,74,75]. These studies produced a variety of computational tools and reference datasets, leaving just a few blind spots.
One of these blind spots is in Russia, where only a handful of genome-wide human variation studies were conducted to date [10, 16, 17, 19, 23]. In this work, using the whole-genome SNP analysis, we surveyed 1019 individuals from Northern Eurasia for their genetic diversity. Newly acquired genome-wide high-density coverage for almost 30 ethnic groups in Russia enabled us to perform both inter-population and population-specific analyses. Combined with genome sequencing data available for a limited number of individuals, such as described in , our study provides one of the most comprehensive datasets covering genetic variation in Russia.
The relationship between genetics and geography was analysed by a combination of ADMIXTURE-based and IBD sharing approaches. We showed that Russian populations of diverse demographic histories and geographic localization share many genetic features, as reflected in their relatively tight ADMIXTURE groupings and outputs of GPS and reAdmix. The apparent positioning of some Russian samples in the genetic space of Caucasus and Siberian populations may reflect either traces of historical assimilation of these groups during the expansion of the Russians, or a recent contribution from neighbouring ethnic groups to the genomes of specific individuals. When we compared Starover and Novosibirsk Russians, representing snapshots of historical (as old as the seventeenth century) and modern (twenty-first century) Russian population, respectively, an apparent recent increase in the rates of admixture with various neighbouring population was evident. Admixture profiles of Modern Novosibirsk Russians have a lower percentage of Northern European components compared to Starovers Russians. In addition, various analyses including GPS, IBD and reAdmix suggest that Starover Russians were genetically influenced by Finno-Ugric people; this hypothesis agrees with the historical record concerning the patterns of Starover migrations within Russian Empire.
One of most curious findings involved the Bashkir, an ethnicity with an extremely complex historical background. There are three main theories describing Bashkir origins: “Turkic”, “Finno-Ugric”, and “Iranian” [76, 77]. According to the “Turkic” theory, most Bashkir genetic ancestry was formed by Turkic tribes migrating from Central Asia in the first millennium AD. The “Finno-Ugric” theory stipulates that the nucleus of Bashkir ancestry was formed by the Magyar (Hungarians), who were later assimilated by Turkic tribes and adopted a Turkic language, while the “Iranian” theory considers Bashkir to be descendants of Sarmatians from the southern Ural.
Speaking generally, our findings add weight to “Finno-Ugric” theory of the origin of Bashkir. A majority of Bashkir IBD fragments were shared with Khanty, an ethnicity related to Magyar. Interestingly, some works point out that before the thirteenth century the Hungarians were commonly called Bashkir (, pp. 289–294). It is surmised that the Magyar ethnicity was formed in the region between Volga and the Ural Mountains, then, at the end of the sixth century AD, moved to the Don-Kuban steppes abandoned by the Proto-Bulgarians followed by the move to their present location between Dnieper and Danube somewhat later.
Further analyses (ADMIXTURE and recent IBD) pointed to proximity of Bashkir to Turkic-speaking Tatar and Chuvash as well as to Finno-Ugric Udmurt and Khanty. In addition, results of f3 outgroup analysis indicate that Bashkir, in contrary to other Turkic speakers, were strongly influenced by Ancient Northern Eurasians, highlighting a mismatch of their cultural background and genetic ancestry and an intricacy of the historic interface between Turkic and Uralic populations. As a general pattern, the Eastern European speakers of Uralic languages share large amounts of IBD with Khanty and Ket, with Turkic speaking Bashkir being added to this rule.
It is noteworthy that the genomes of closest linguistic relatives of Bashkir, Volga Tatar, bears very little traces of East Asian or Central Siberian ancestry. Volga Tatar are a mix between Bulgar who carried a large Finno-Ugric component, Pecheneg, Kuman, Khazar, local Finno-Ugric tribes, and even Alan. Therefore, Volga Tatars are predominantly European ethnicity with a tiny contribution of East-Asian component. As most Tatar’ IBD is shared with various Turkic and Uralic populations from Volga-Ural region, an amalgamation of various cultures is evident. When the original Finno-Ugric speaking people were conquered by Turkic tribes, both Tatar and Chuvash are likely to have experience language replacement, while retaining their genetic core. Most likely, these events took place sometime around VIII century AD, after the relocation of Bulgar tribes to Volga and Kama river basins, and expansion of Turkic people.
We speculate that Bashkir, Tatar, Chuvash and Finno-Ugric speakers from Volga basin has a common Turkic component, which could have been acquired as a result of Turkic expansion to Volga-Urals region. However, the original Finno-Ugric substrate was not homogeneous: Tatar and Chuvash genomes carry mainly “Finno-Permic” component, while Bashkir carry the “Magyar” one. The fraction of the Turkic component in Bashkir is, undoubtedly, quite significant, and larger than that in Tatar and Chuvash. This component reflects the South Siberian influence on Bashkir, which makes them related to Altai, Kyrgyz, Tuvinian, and Kazakh people.
As a standalone approach, an analysis of shared IBD is not sufficient to support the Finno-Ugric hypothesis of Bashkir origin as a sole source, while pointing at temporal separation of genetic components in Bashkir. Hence, we demonstrated that Bashkir genepool is a multifaceted, multicomponent system, lacking the main “core”; it is an amalgamation of Turkic, Ugric, Finnish and Indo-European contributions. In this mosaic, it is impossible to identify the leading element. Therefore, Bashkir are the most genetically diverse ethnic group of the Volga-Urals region.
Many Siberian populations share an unusually high amount of IBD, which may be explained by a combination of the following factors: 1) shared origin, 2) relative isolation from outside world, 3) rapid recent population growth and strong founder effect in Yakut, Buryat, and Tuva, or 4) gene flow facilitated by some migrating population. The structure of these population also reflects the role of multiple South-North travel routes along the great waterways of Ob, Yenisei, and Lena, while the Siberian taiga, which is notoriously hard to traverse, to some degree prevented lateral access. On the other hand, Southern Siberia, where the steppes border the forests, is easier to travel. The same is true for the Northern Siberia, where the cold, flat tundra is suitable for travel by deer herders. These geographical limitations corralled the East-West migration to either “northern” or “southern” corridors and North-South migrations to the banks of great Siberian rivers. The footprints of these geographical restrictions could be seen in the patterns of IBD sharing between the Siberian populations studied. We christened it as the “Siberian genetic vortex”.
High IBD between West Siberian Ket and Khanty populations may reflect their relatively recent admixture with Selkup. Close genetic relationships between Ket and Tuvan can be explained by the existence of an ancient pre-Turkic and pre-Samoyedic Yenisei substrate which constitutes the main genetic component in Ket and still present in Tuva due to assimilation of extinct Yeniseian peoples (such as Kott, Arin, and Pumpokol)  inhabited Yenisei source area in the Southern Siberia .
High levels of shared IBD blocks in Altaic-speaking populations from Southern Siberia (Tuva, Buryat), North Asia (Yakut) and Central Asia (Kyrgyz) supports their recently formed common genetic core, which is geographically related to the Altay-Sayan Mountains region in Southern Siberia. Yakut and Kyrgyz populations which are now distant from this region were resettled from Southern Siberia relatively recently. It is accepted that ancestors of Yakuts (Kurykan) migrated from the Southern Yenisei to Lake Baikal area in seventh century AD, and then travelled the Lena river North in 12th -14th centuries AD , while Kyrgyz, who until recently were known as Yenisei Kyrgyz, migrated from Southern Siberia to Central Asia in 13th - 15th centuries AD after the collapse of the Mongol Empire .
The discovery of long runs of homozygosity in native Siberian populations (such as Tuva, Buryat, Yakut, Ket, Khanty) supports the earlier finding of pronounced founder effects and low genetic diversity in Siberians due to genetic drift, isolation by distance and recent population expansion events, that were made using the Y-chromosome analysis [83,84,85,86,87,88,89].
Comparative analysis of modern and ancient genomes suggests that Western Siberians have more Ancient North European ancestry (represented by Mal’ta) than other populations of the Russian Federation. Other studied populations show genetic affinity to various ancient genomes, either co-located with modern inhabitants, pointing to direct gene flow and relatively sessile population, or geographically removed, pointing to their migration to currently occupied locations.
We see that the shared genetic drift associated with hunter gatherers (Fig. 8) is correlated with Northern European ancestry of studied individuals. At the same time, the shared genetic drift of farmers has a pronounced gradient: it is large in the areas suitable for agriculture and drops to zero in Ket and Khanty-inhabited boreal forest areas of Siberia, where the climate is harsh and summers are too short for a sustainable harvest. In Siberian forests, the signal of Neolithic ancestry is no longer detected, but the ancient northern Eurasian (ANE) signal predominates instead. Possibly, the ancient Northern Eurasians met with more western groups of ancient hunters or with ancient farmers in the steppe, formed a certain population resembling the steppe samples of Yamnaya and Afanasyevo cultures, which then spread this North Eurasian component across and beyond the boreal forests of Siberia. This suggests an extensive westward migration from the steppe, discussed in detail elsewhere . It is also possible that there was wave of northern or western Europeans migrating to the steppes from an opposite direction.
Our project has filled an important lacuna in the genetic map of Eurasia. We revealed the complexity of genetic structure of Northern Eurasians, the existence of East-West and North-South genetic gradients, and varying inputs of ancient populations into modern populations. In particular, we have collected evidence in support of Finno-Ugric influence on the formation of Bashkir, shed light on
to the genetic make-up of Russian Starovers (Old Believers), and postulated the existence of a Great Siberian Vortex directing genetic exchanges in populations across the Siberian part of Asia.
Cavalli-Sforza LL, Menozzi P, Piazza A, Princeton NJ. The history and geography of human genes, vol. 541: Princeton University Press; SRC - GoogleScholar. 1994.
Novembre J, Johnson T, Bryc K, Kutalik Z, Boyko AR, Auton A, Indap A, King KS, Bergmann S, Nelson MR, et al. Genes mirror geography within Europe. Nature. 2008;456(7218):98–101.
Cavalli-Sforza LL. The human genome diversity project: past, present and future. Nat Rev Genet. 2005;6(4):333–40.
Sabeti P, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne E, McCarroll S, Gaudet R, et al. The international HapMap project. Nature. 2003;426(6968):789–96.
Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422 SRC - GoogleScholar):56–65.
Abdulla MA, Ahmed I, Assawamakin A, Bhak J, Brahmachari SK, Calacal GC. Mapping human genetic diversity in Asia. Science. 2009;326(5959 SRC - GoogleScholar):1541–5.
Elhaik E, Greenspan E, Staats S, Krahn T, Tyler-Smith C, Xue Y, Tofanelli S, Francalacci P, Cucca F, Pagani L, et al. The GenoChip: a new tool for genetic anthropology. Gen Biol Evol. 2013;5(5):1021–31.
Elhaik E, Tatarinova T, Chebotarev D, Piras IS, Maria Calò C, De Montis A, Atzori M, Marini M, Tofanelli S, Francalacci P, et al. Geographic population structure analysis of worldwide human populations infers their biogeographical origins. Nat Commun. 2014;5:3513.
Raghavan M, Skoglund P, Graf KE, Metspalu M, Albrechtsen A, Moltke I, Rasmussen S, Stafford TW Jr, Orlando L, Metspalu E, et al. Upper Palaeolithic Siberian genome reveals dual ancestry of native Americans. Nature. 2014;505(7481):87–91.
Flegontov P, Changmai P, Zidkova A, Logacheva MD, Flegontova O, Gelfand MS, Gerasimov ES, Khrameeva E, Konovalova OP, Neretina T, et al. Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient north Eurasian ancestry. Sci Rep. 2016;
Raghavan M, Steinrucken M, Harris K, Schiffels S, Rasmussen S, DeGiorgio M, Albrechtsen A, Valdiosera C, Avila-Arcos MC, Malaspinas AS, et al. POPULATION GENETICS. Genomic evidence for the Pleistocene and recent population history of native Americans. Science. 2015;349(6250):aab3884.
Skoglund P, Malmstrom H, Raghavan M, Stora J, Hall P, Willerslev E, Gilbert MT, Gotherstrom A, Jakobsson M. Origins and genetic legacy of Neolithic farmers and hunter-gatherers in Europe. Science. 2012;336(6080):466–9.
Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, et al. Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008;319(5866):1100–4.
Nelis M, Esko T, Magi R, Zimprich F, Zimprich A, Toncheva D, Karachanak S, Piskackova T, Balascak I, Peltonen L, et al. Genetic structure of Europeans: a view from the north-east. PLoS One. 2009;4(5):e5472.
Yunusbayev B, Metspalu M, Jarve M, Kutuev I, Rootsi S, Metspalu E, Behar DM, Varendi K, Sahakyan H, Khusainova R, et al. The Caucasus as an asymmetric semipermeable barrier to ancient human migrations. Mol Biol Evol. 2012;29(1):359–65.
Mallick S, Li H, Lipson M, Mathieson I, Gymrek M, Racimo F, Zhao M, Chennagiri N, Nordenfelt S, Tandon A, et al. The Simons genome diversity project: 300 genomes from 142 diverse populations. Nature. 2016;538(7624):201–6.
Pagani L, Lawson DJ, Jagoda E, Morseburg A, Eriksson A, Mitt M, Clemente F, Hudjashov G, DeGiorgio M, Saag L, et al. Genomic analyses inform on migration events during the peopling of Eurasia. Nature. 2016;538(7624):238–42.
Kushniarevich A, Utevska O, Chuhryaeva M, Agdzhoyan A, Dibirova K, Uktveryte I, Mols M, Mulahasanovic L, Pshenichnov A, Frolova S, et al. Genetic heritage of the Balto-Slavic speaking populations: a synthesis of Autosomal, mitochondrial and Y-chromosomal data. PLoS One. 2015;10(9):e0135820.
Yunusbayev B, Metspalu M, Metspalu E, Valeev A, Litvinov S, Valiev R, Akhmetova V, Balanovska E, Balanovsky O, Turdikulova S, et al. The genetic legacy of the expansion of Turkic-speaking nomads across Eurasia. PLoS Genet. 2015;11(4):e1005068.
Fedorova SA, Reidla M, Metspalu E, Metspalu M, Rootsi S, Tambets K, Trofimova N, Zhadanov SI, Hooshiar Kashani B, Olivieri A, et al. Autosomal and uniparental portraits of the native populations of Sakha (Yakutia): implications for the peopling of northeast Eurasia. BMC Evol Biol. 2013;13:127.
Behar DM, Metspalu M, Baran Y, Kopelman NM, Yunusbayev B, Gladstein A, Tzur S, Sahakyan H, Bahmanimehr A, Yepiskoposyan L, et al. No evidence from genome-wide data of a Khazar origin for the Ashkenazi Jews. Hum Biol. 2013;85(6):859–900.
Behar DM, Yunusbayev B, Metspalu M, Metspalu E, Rosset S, Parik J, Rootsi S, Chaubey G, Kutuev I, Yudkovsky G, et al. The genome-wide structure of the Jewish people. Nature. 2010;466(7303):238–42.
Wong EH, Khrunin A, Nichols L, Pushkarev D, Khokhrin D, Verbenko D, Evgrafov O, Knowles J, Novembre J, Limborska S, et al. Reconstructing genetic history of Siberian and northeastern European populations. Genome Res. 2017;27(1):1–14.
Clemente FJ, Cardona A, Inchley CE, Peter BM, Jacobs G, Pagani L, Lawson DJ, Antao T, Vicente M, Mitt M, et al. A selective sweep on a deleterious mutation in CPT1A in Arctic populations. Am J Hum Genet. 2014;95(5):584–9.
Raghavan M, DeGiorgio M, Albrechtsen A, Moltke I, Skoglund P, Korneliussen TS, Gronnow B, Appelt M, Gullov HC, Friesen TM, et al. The genetic prehistory of the new world Arctic. Science. 2014;345(6200):1255832.
Skryabin KG, Prokhortchouk EB, Mazur AM, Boulygina ES, Tsygankova SV, Nedoluzhko AV, Rastorguev SM, Matveev VB, Chekanov NN, AG D, et al. Combining two technologies for full genome sequencing of human. Acta Nat. 2009;1(3):102–7.
Zenkovsky SA: Russian Starovers in XVII-XIX [Зеньковский С.А. Русское старообрядчество: в 2 т. - 3-е изд., испр. и доп. - М.: Институт ДИ-ДИК, 2016. - 712 с., 70х100/16] ISBN 978-5-93311-013-2. 2016.
Eller E. Population substructure and isolation by distance in three continental regions. Am J Phys Anthropol. 1999;108(2):147–59.
Relethford JH. Global analysis of regional differences in craniometric diversity and population substructure. Hum Biol. 2001;73(5):629–36.
Ramachandran S, Deshpande O, Roseman CC, Rosenberg NA, Feldman MW, Cavalli-Sforza LL. Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc Natl Acad Sci U S A. 2005;102(44):15942–7.
ArunKumar G, Tatarinova TV, Duty J, Rollo D, Syama A, Arun VS, Kavitha VJ, Triska P, Greenspan B, Wells RS, et al. Genome-wide signatures of male-mediated migration shaping the Indian gene pool. J Hum Genet. 2015;60(9):493–9.
Yang WY, Platt A, Chiang CW, Eskin E, Novembre J, Pasaniuc B. Spatial localization of recent ancestors for admixed individuals. G3 (Bethesda). 2014;4(12):2505–18.
Yang WY, Novembre J, Eskin E, Halperin E. A model-based approach for analysis of spatial structure in genetic data. Nat Genet. 2012;44(6):725–31.
François O, Currat M, Ray N, Han E, Excoffier L, Novembre J. Principal component analysis under population genetic models of range expansion and admixture. Mol Biol Evol. 2010;27(6):1257–68.
Novembre J, Ramachandran S. Perspectives on human population structure at the cusp of the sequencing era. Annu Rev Genomics Hum Genet. 2011;12:245–74.
Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, Srinivasan BS, Barsh GS, Myers RM, Feldman MW, et al. Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009;19(5):826–37.
Royal CD, Novembre J, Fullerton SM, Goldstein DB, Long JC, Bamshad MJ, Clark AG. Inferring genetic ancestry: opportunities, challenges, and implications. Am J Hum Genet. 2010;86(5):661–73.
Kozlov KN, Samsonov AM, Samsonova MG. Method of entirely parallel differential evolution for model adaptation in systems biology. Biofizika. 2015;60(6):1219–20.
Kozlov K, Chebotarev D, Hassan M, Triska M, Triska P, Flegontov P, Tatarinova TV. Differential evolution approach to detect recent admixture. BMC Genomics. 2015;16(Suppl 8):S9.
Reich D, Thangaraj K, Patterson N, Price AL, Singh L. Reconstructing Indian population history. Nature. 2009;461(7263):489–94.
ArunKumar G, David F, Soria-Hernanz, Kavitha VJ, Arun VS, Syama A, Ashokan KS, Gandhirajan KT, Vijayakumar K, Narayanan M, Jayalakshmi M, et al. Population differentiation of southern Indian male lineages correlates with agricultural expansions predating the caste system. PLoS One. 2012;7(11):e50269.
ArunKumar G, Tatarinova TV, Duty J, Rollo D, Syama A, Arun VS, Kavitha VJ, Triska P, Greenspan B, Wells RS, Pitchappan R. Genographic Consortium. Genome-wide signatures of male-mediated migration shaping the Indian gene pool. J Hum Genet. 2015;60(9):493–9.
Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen W-M. Robust relationship inference in genome-wide association studies. Bioinformatics (Oxford, England). 2010;26(22):2867–73.
Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19(9):1655–64.
Alexander DH, Lange K, Admixture C. BM: enhancements to the for individual ancestry estimation. BMC Bioinformatics. 2011;12(1 SRC - GoogleScholar):246.
Browning BL, Browning SR. A fast, powerful method for detecting identity by descent. Am J Hum Genet. 2011;88(2):173–82.
Pickrell JK, Pritchard JK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012;8(11):e1002967.
Elhaik E, Tatarinova TV, Klyosov AA, Graur D. The ‘extremely ancient’ chromosome that isn't: a forensic bioinformatic investigation of Albert Perry’s X-degenerate portion of the Y chromosome. Eur J Hum Genet. 2014;22(9):1111–6.
Ringbauer H, Coop G, Barton NH. Inferring recent demography from isolation by distance of long shared sequence blocks. Genetics. 2017;205(3):1335–51.
Napolskikh VV. [Vvedenije v istoricheskuju uralistiku] Introduction to Historical Uralistics. Izhevsk. 1997. http://elibrary.unatlib.ru/handle/123456789/22075.
Karafet TM, Bulayeva KB, Nichols J, Bulayev OA, Gurgenova F, Omarova J, Yepiskoposyan L, Savina OV, Rodrigue BH, Hammer MF. Coevolution of genes and languages and high levels of population structure among the highland populations of Daghestan. J Hum Genet. 2016;61(3):181–91. doi:10.1038/jhg.2015.132. Epub 2015 Nov 26. PMID: 26607180.
Karafet TM, Bulayeva KB, Bulayev OA, Gurgenova F, Omarova J, Yepiskoposyan L, Savina OV, Veeramah KR, Hammer MF. Extensive genome-wide autozygosity in the population isolates of Daghestan. Eur J Hum Genet. 2015;23(10):1405–12.
Haber M, Mezzavilla M, Xue Y, Comas D, Gasparini P, Zalloua P, Tyler-Smith C. Genetic evidence for an origin of the Armenians from bronze age mixing of multiple populations. Eur J Hum Genet. 2016;24(6):931–6.
Balanovsky O, Dibirova K, Dybo A, Mudrak O, Frolova S, Pocheshkhova E, Haber M, Platt D, Schurr T, Haak W, et al. Parallel evolution of genes and languages in the Caucasus region. Mol Biol Evol. 2011;28(10):2905–20.
Elhaik E, Greenspan E, Staats S, Krahn T, Tyler-Smith C, Xue Y, Tofanelli S, Francalacci P, Cucca F, Pagani L, et al. The GenoChip: a new tool for genetic anthropology. Genome Biol Evol. 2013;
Haak W, Lazaridis I, Patterson N, Rohland N, Mallick S, Llamas B, Brandt G, Nordenfelt S, Harney E, Stewardson K, et al. Massive migration from the steppe was a source for indo-European languages in Europe. Nature. 2015;
Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, Harney E, Stewardson K, Fernandes D, Novak M, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015;528(7583):499–503.
Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, Sudmant PH, Schraiber JG, Castellano S, Lipson M, et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014;513(7518):409–13.
Morozova I, Flegontov P, Mikheyev AS, Bruskin S, Asgharian H, Ponomarenko P, Klyuchnikov V, ArunKumar G, Prokhortchouk E, Gankin Y, et al. Toward high-resolution population genomics using archaeological samples. DNA Res. 2016;23(4):295–310.
Der Sarkissian C, Balanovsky O, Brandt G, Khartanovich V, Buzhilova A, Koshel S, Zaporozhchenko V, Gronenborn D, Moiseyev V, Kolpakov E, et al. Ancient DNA reveals prehistoric gene-flow from siberia in the complex human population history of north East Europe. PLoS Genet. 2013;9(2):e1003296.
Patterson N, Moorjani P, Luo Y, Mallick S, Rohland N, Zhan Y, Genschoreck T, Webster T, Reich D. Ancient admixture in human history. Genetics. 2012;192(3):1065–93.
Peter BM. Admixture, population structure, and F-statistics. Genetics. 2016;202(4):1485–501.
Bramanti B, Thomas MG, Haak W, Unterlaender M, Jores P, Tambets K, Antanaitis-Jacobs I, Haidle MN, Jankauskas R, Kind CJ, et al. Genetic discontinuity between local hunter-gatherers and central Europe's first farmers. Science. 2009;326(5949):137–40.
Pereira JB, Costa MD, Vieira D, Pala M, Bamford L, Harich N, Cherni L, Alshamali F, Hatina J, Rychkov S, Stefanescu G, King T, Torroni A, Soares P, Pereira L, Richards MB. Reconciling evidence from ancient and contemporary genomes: a major source for the European Neolithic within Mediterranean Europe. Proc Biol Sci. 2017;284(1851). doi:10.1098/rspb.2016.1976.
Hofmanova Z, Kreutzer S, Hellenthal G, Sell C, Diekmann Y, Diez-Del-Molino D, van Dorp L, Lopez S, Kousathanas A, Link V, et al. Early farmers from across Europe directly descended from Neolithic Aegeans. Proc Natl Acad Sci U S A. 2016;113(25):6886–91.
Heath SC, Gut IG, Brennan P, McKay JD, Bencko V, Fabianova E, Foretova L, Georges M, Janout V, Kabesch M, et al. Investigation of the fine structure of European populations with applications to disease association studies. Eur J Hum Genet : EJHG. 2008;16(12):1413–29.
Seldin MF, Shigeta R, Villoslada P, Selmi C, Tuomilehto J, Silva G, Belmont JW, Klareskog L, Gregersen PK. European population substructure: clustering of northern and southern populations. PLoS Genet. 2006;2(9):e143.
Tian C, Plenge RM, Ransom M, Lee A, Villoslada P, Selmi C, Klareskog L, Pulver AE, Qi L, Gregersen PK, et al. Analysis and application of European genetic substructure using 300 K SNP information. PLoS Genet. 2008;4(1):e4.
Price AL, Butler J, Patterson N, Capelli C, Pascali VL, Scarnicci F, Ruiz-Linares A, Groop L, Saetta AA, Korkolopoulou P, et al. Discerning the ancestry of European Americans in genetic association studies. PLoS Genet. 2008;4(1):e236.
Bauchet M, McEvoy B, Pearson LN, Quillen EE, Sarkisian T, Hovhannesyan K, Deka R, Bradley DG, Shriver MD. Measuring European population stratification with microarray genotype data. Am J Hum Genet. 2007;80(5):948–56.
International HapMap C, Altshuler DM, Gibbs RA, Peltonen L, Altshuler DM, Gibbs RA, Peltonen L, Dermitzakis E, Schaffner SF, Yu F, et al. Integrating common and rare genetic variation in diverse human populations. Nature. 2010;467(7311):52–8.
McVean G. A genealogical interpretation of principal components analysis. PLoS Genet. 2009;5(10):e1000686.
Rodriguez-Flores JL, Fakhro K, Agosto-Perez F, Ramstetter MD, Arbiza L, Vincent TL, Robay A, Malek JA, Suhre K, Chouchane L, et al. Indigenous Arabs are descendants of the earliest split from ancient Eurasian populations. Genome Res. 2016;26(2):151–62.
Wall JD, Yang MA, Jay F, Kim SK, Durand EY, Stevison LS, Gignoux C, Woerner A, Hammer MF, Slatkin M. Higher levels of neanderthal ancestry in east Asians than in Europeans. Genetics. 2013;194(1):199–209.
Petkova D, Novembre J, Stephens M. Visualizing spatial population structure with estimated effective migration surfaces. Nat Genet. 2016;48(1):94–100.
Yamaeva LA. On ehnogenesis of Bashkir: culturological approach--K voprosu ob etnogenese bashkir: kulturologiheskij podxod. In: International scientific conference dedicted to the memory of R M Yusupov. Ufa: Russian Academy of Science; 2011. p. 354–8.
Bashkir ethnogenesis, Yanguzin R., Vatandash Online Journal, 2004, available at http://vatandash.ru/index.php?article=1079.
Róna-Tas A. Hungarians and Europe in the early middle ages: an introduction to early Hungarian history. Budapest, Hungary: Central European University Press; 1999.
Vajda EJ: Languages and prehistory of Central Siberia. Current issues in linguistic theory (presentation of the Yeniseian family and its speakers, together with neighboring languages and their speakers, in linguistic, historical and archeological view). John Benjamin Publishing Company; 2004.
Funk DA, Alexeev NA. Turkic people of eastern Siberia. Moscow: Nauka; 2008.
Alekseev NA, Romanova EN, Sokolova ZP. Yakuts (Saha). In: Yakuts (Saha). Moscow: Nauka; 2012. p. 599.
Butanayev VY, Khudyakov YS. History of Yeniseyan Kyrgyz. Abakan: Katanov Khakass State University; 2000.
Karafet TM, Osipova LP, Gubina MA, Posukh OL, Zegura SL, Hammer MF. High levels of Y-chromosome differentiation among native Siberian populations and the genetic signature of a boreal hunter-gatherer way of life. Hum Biol. 2002;74(6):761–89.
Har'kov VN, Hamina KV, Medvedeva OF, Simonova KV, Eremina ER, Stepanov VA. Gene pool of Buryats: clinal variability and territorial subdivision based on data of Y-chromosome markers. Genetika. 2014;50(2):203–13.
Kharkov V, Khamina K, Medvedeva O, Simonova K, Khitrinskaya I, Stepanov V. Gene-pool structure of Tuvinians inferred from Y-chromosome marker data. Genetika. 2013;49(12):1416–25.
Derenko MV, Malyarchuk BA, Denisova GA, Dorzhu ChM, Karamchakova ON, Luzina FA, Lotosh EA, Dambueva IK, Ondar UN, Zakharov IA. Polymorphism of the Y-Chromosome Diallelic Loci in Ethnic Groups of the Altai–Sayan Region, Russian Journal of Genetics. 2002;38(3):309–14.
Khar'kov VN, Khamina KV, Medvedeva OF, Stepanov VA, Shtygasheva OV. Genetic diversity of the Khakass gene pool: subethnic differentiation and the structure of Y-chromosome haplogroups. Mol Biol. 2011;45(3):404–16.
Khar'kov VN, Stepanov VA, Medvedev OF, Spiridonova MG, Maksimova NR, Nogovitsyna AN, Puzyrev VP. The origin of Yakuts: analysis of Y-chromosome haplotypes. Mol Biol. 2008;42(2):226–37.
Kharkov VN, Stepanov VA, Medvedeva OF, Spiridonova MG, Puzyrev VP, Maksimova NR, Nogovitsina AN. The origin of Yakuts: analysis of the Y-chromosome haplotypes. Mol Biol. 2008;42(2):198–208.
We thank Prof. Roger Jelliffe for careful reading of the manuscript and Prof. Mikhail Kovalchuk for comprehensive assistance with the work.
TVT, PT, MT were supported by grants from National Institute of General Medical Sciences (GM068968), Eunice Kennedy Shriver National Institute of Child Health and Human Development (HD070996), National Science Foundation Division of Evolutionary Biology (1456634), and National Science Foundation Small Business Technology Transfer (1622840). TVT and OB were supported by Russian Scientific Foundation (17–14-01345). AB was supported by National Science Foundation Small Business Technology Transfer (1622840). VS was supported by Russian Scientific Foundation (16–15-00020), Russian Foundation for Basic Research (15–04-02442). VK was supported by Russian Foundation for Basic Research (16–34-60,222) and Grant of President of Russian Federation (MD-8886.2016.4). Publication costs were funded by the authors themselves.
Availability of data and materials
The genotype data (core dataset) are available (after the paper is published) from www.RussianGenome.ru.
About this supplement
This article has been published as part of BMC Genetics Volume 18 Supplement 1, 2017: Selected articles from Belyaev Conference 2017: genetics. The full contents of the supplement are available online at https://bmcgenet.biomedcentral.com/articles/supplements/volume-18-supplement-1.
Ethics approval and consent to participate
The study protocol was approved by the Ethics Committee of the Research Institute of Medical Genetics, Tomsk National Medical Research Center of the Russian Academy of Science, the Ethics Committee of the Institute of Biochemistry and Genetics of Ufa Scientific Centre Russian Academy of Science, and the Ethics Committee of the Center “Bioengineering” (BNG) Russian Academy of Science. Written informed consent was obtained from all participants in the study. For each participant, information about geographic and ethnic origins of their parents and grandparents was recorded.
Consent for publication
On behalf of all authors, the corresponding author states that there is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Quality control process. (PDF 25 kb)
Results of ADMIXTURE for K = 2–10. (PNG 157 kb)
List of samples included in “Extended” dataset. Table S1b List of samples included in “Core” dataset. Table S1c List of samples included in “Ancient” dataset. Table S2 Results of ADMIXTURE for K = 9. Table S3 Results of ADMIXTURE for K = 6, 7, 8. Table S4 Results of f3 test. Table S5 Results of IBD sharing analysis in 1–3 cM and 4–10 cM bins. Table S6 Total amount of shared IBD between populations. Table S7 Standard residue of linear regression analysis of distance-IBD sharing. Table S8 Distance and shared IBD between pairs of populations. Table S9: Results of f3 outgroup test with ancient samples. (XLSX 482 kb)