Autosomal and uniparental portraits of the native populations of Sakha (Yakutia): implications for the peopling of Northeast Eurasia

High prevalence of m.1555A > G in patients with hearing loss in the Baikal Lake region of Russia as a result of founder effect

Mitochondrial forms account approximately 1-2% of all nonsyndromic cases of hearing loss (HL). One of the most common causative variants of mtDNA is the m.1555A > G variant of the MT-RNR1 gene (OMIM 561000). Currently the detection of the m.1555A > G variant of the MT-RNR1 gene is not included in all research protocols. In this study this variant was screened among 165 patients with HL from the Republic of Buryatia, located in the Baikal Lake region of Russia. In our study, the total contribution of the m.1555A > G variant to the etiology of HL was 12.7% (21/165), while the update global prevalence of this variant is 1.8% (863/47,328). The m.1555A > G variant was notably more prevalent in Buryat (20.2%) than in Russian patients (1.3%). Mitogenome analysis in 14 unrelated Buryat families carrying the m.1555A > G variant revealed a predominant lineage: in 13 families, a cluster affiliated with sub-haplogroup A5b (92.9%) was identified, while one family had the D5a2a1 lineage (7.1%). In a Russian family with the m.1555A > G variant the lineage affiliated with sub-haplogroup F1a1d was found. Considering that more than 90% of Buryat families with the m.1555A > G variant belong to the single maternal lineage cluster we conclude that high prevalence of this variant in patients with HL in the Baikal Lake region can be attributed to a founder effect.

Genetic history of the Koryaks and Evens of the Magadan region based on Y chromosome polymorphism data

In order to clarify the history of gene pool formation of the indigenous populations of the Northern Priokhotye (the northern coast of the Sea of Okhotsk), Y-chromosome polymorphisms were studied in the Koryaks and Evens living in the Magadan region. The results of the study showed that the male gene pool of the Koryaks is represented by haplogroups C-B90-B91, N-B202, and Q-B143, which are also widespread in other peoples of Northeastern Siberia, mainly of Paleo-Asiatic origin. High frequency of haplogroup C-B80, typical of other Tungus-Manchurian peoples, is characteristic of the Evens of the Magadan region. The shared components of the gene pools of the Koryaks and Evens are haplogroups R-M17 and I-P37.2 inherited as a result of admixture with Eastern Europeans (mainly Russians). The high frequency of such Y chromosome haplogroups in the Koryaks (16.7 %) and Evens (37.8 %) is indicative of close interethnic contacts during the last centuries, and most probably especially during the Soviet period. The genetic contribution of the European males' Y chromosome significantly prevails over that of maternally inherited mitochondrial DNA. The study of the Y chromosome haplogroup diversity has shown that only relatively young phylogenetic branches have been preserved in the Koryak gene pool. The age of the oldest component of the Koryak gene pool (haplogroup C-B90-B91) is estimated to be about 3.8 thousand years, the age of the younger haplogroups Q-B143 and N-B202 is about 2.8 and 2.4 thousand years, respectively. Haplogroups C-B90-B91 and N-B202 are Siberian in origin, and haplogroup Q-B143 was apparently inherited by the ancestors of the Koryaks and other Paleo-Asiatic peoples from the Paleo-Eskimos as a result of their migrations to Northeast Asia from the Americas. The analysis of microsatellite loci for haplogroup Q-B143 in the Eskimos of Greenland, Canada and Alaska as well as in the indigenous peoples of Northeastern Siberia showed a decrease in genetic diversity from east to west, pointing to the direction of distribution of the Paleo-Eskimo genetic component in the circumpolar region of America and Asia. At the same time, the Evens appeared in the Northern Priokhotye much later (in the XVII century) as a result of the expansion of the Tungusic tribes, which is confirmed by the results of the analysis of haplogroup C-B80 polymorphisms.

Epigenetics of the far northern Yakutian population

BACKGROUND

Yakuts are one of the indigenous populations of the subarctic and arctic territories of Siberia characterized by a continental subarctic climate with severe winters, with the regular January average temperature in the regional capital city of Yakutsk dipping below - 40 °C. The epigenetic mechanisms of adaptation to such ecologies and environments and, in particular, epigenetic age acceleration in the local population have not been studied before.

RESULTS

This work reports the first epigenetic study of the Yakutian population using whole-blood DNA methylation data, supplemented with the comparison to the residents of Central Russia. Gene set enrichment analysis revealed, among others, geographic region-specific differentially methylated regions associated with adaptation to climatic conditions (water consumption, digestive system regulation), aging processes (actin filament activity, cell fate), and both of them (channel activity, regulation of steroid and corticosteroid hormone secretion). Further, it is demonstrated that the epigenetic age acceleration of the Yakutian representatives is significantly higher than that of Central Russia counterparts. For both geographic regions, we showed that epigenetically males age faster than females, whereas no significant sex differences were found between the regions.

CONCLUSIONS

We performed the first study of the epigenetic data of the Yakutia cohort, paying special attention to region-specific features, aging processes, age acceleration, and sex specificity.

Single nucleotide polymorphism rs527236194 of the cytochrome B gene (MT-CYB) is associated with alterations in sperm parameters

BACKGROUND

The mitochondrial genome is substantially susceptible to mutations and has high polymorphism due to structural features, location, and lack of recombinant variability, as its inheritance is strictly maternal. All of these events can be accompanied by the accumulation of mitochondrial single nucleotide polymorphisms (mtSNPs) in the sperm. The aim of this research was to analyze the influence of mutations in the MT-CYB gene on sperm quality.

METHODS AND RESULTS

We conducted a case‒control study to identify mutations in the mitochondrial cytochrome B (MT-CYB) gene in men with asthenoteratozoospermia (89 cases) and oligoasthenoteratozoospermia (65 cases). The comparison group consisted of 164 fertile men. Somatic cell lysis followed by mtDNA extraction was conducted to analyze three mtDNA polymorphisms, rs28357373 (T15629C (Leu295=), rs527236194 (T15784C (p.Pro346=), rs2853506 (A15218G, p.Thr158Ala). Detection and genotyping of polymorphic loci in the MT-CYB gene was performed using the TaqMan allelic discrimination assay. To verify mutations in the MT-CYB gene, automated Sanger DNA sequencing was used. We found that rs527236194 was associated with asthenoteratozoospermia. rs28357373 in the MT-CYB gene did not show any polymorphism in the analyzed groups, which indicates a rare frequency of the TT genotype in our region. Rs28357373 and rs2853506 are not associated with male sperm abnormalities in the Volga-Ural region.

CONCLUSION

The association of the rs527236194 polymorphic variant with sperm parameter alterations suggests its role in the pathophysiology of male infertility and requires further investigation in larger samples.

Mitogenomics of the Koryaks and Evens of the northern coast of the Sea of Okhotsk

Due to the geographical proximity of the northern coast of the Sea of Okhotsk and Kamchatka Peninsula to the Beringia, the indigenous populations of these territories are of great interest for elucidating the human settlement history of northern Asia and America. Meanwhile, there is a clear shortage of genetic studies of the indigenous populations of the northern coast of the Sea of Okhotsk. Here, in order to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 203 complete mitogenomes (174 of which are new) from population samples of the Koryaks and Evens of the northern coast of the Sea of Okhotsk and the Chukchi of the extreme northeast Asia. The patterns observed underscore the reduced level of genetic diversity found in the Koryak, Even, and Chukchi populations, which, along with the high degree of interpopulation differentiation, may be the result of genetic drift. Our phylogeographic analysis reveals common Paleo-Asiatic ancestry for 51.1% of the Koryaks and 17.8% of the Evens. About third of the mitogenomes found in the Koryaks and Evens might be considered as ethno-specific, as these are virtually absent elsewhere in North, Central and East Asia. Coalescence ages of most of these lineages coincide well with the emergence and development of the Tokarev and Old Koryak archaeological cultures associated with the formation of the Koryaks, as well as with the period of separation and split of the North Tungusic groups migrated northwards from the Lake Baikal or the Amur River area.

Lactase deficiency in Russia: multiethnic genetic study

BACKGROUND

Lactase persistence-the ability to digest lactose through adulthood-is closely related to evolutionary adaptations and has affected many populations since the beginning of cattle breeding. Nevertheless, the contrast initial phenotype, lactase non-persistence or adult lactase deficiency, is still observed in large numbers of people worldwide.

METHODS

We performed a multiethnic genetic study of lactase deficiency on 24,439 people, the largest in Russia to date. The percent of each population group was estimated according to the local ancestry inference results. Additionally, we calculated frequencies of rs4988235 GG genotype in Russian regions using the information of current location and birthplace data from the client's questionnaire.

RESULTS

The attained results show that among all studied population groups, the frequency of GG genotype in rs4988235 is higher than the average in the European populations. In particular, the prevalence of lactase deficiency genotype in the East Slavs group was 42.8% (95% CI: 42.1-43.4%). We also investigated the regional prevalence of lactase deficiency based on the current place of residence.

CONCLUSIONS

Our study emphasizes the significance of genetic testing for diagnostics, i.e., specifically for lactose intolerance parameter, as well as the scale of the problem of lactase deficiency in Russia which needs to be addressed by the healthcare and food sectors.

The matrilineal ancestry of Nepali populations

The Tibetan plateau and high mountain ranges of Nepal are one of the challenging geographical regions inhabited by modern humans. While much of the ethnographic and population-based genetic studies were carried out to investigate the Tibetan and Sherpa highlanders, little is known about the demographic processes that enabled the colonization of the hilly areas of Nepal. Thus, the present study aimed to investigate the past demographic events that shaped the extant Nepalese genetic diversity using mitochondrial DNA (mtDNA) variations from ethnic Nepalese groups. We have analyzed mtDNA sequences of 999 Nepalese and compared data with 38,622 published mtDNA sequences from rest of the world. Our analysis revealed that the genomic landscapes of prehistoric Himalayan settlers of Nepal were similar to that of the low-altitude extant Nepalese (LAN), especially Newar and Magar population groups, but differ from contemporary high-altitude Sherpas. LAN might have derived their East Eurasian ancestry mainly from low-altitude Tibeto-Burmans, who likely have migrated from East Asia and assimilated across the Eastern Himalayas extended from the Eastern Nepal to the North-East of India, Bhutan, Tibet and Northern Myanmar. We also identified a clear genetic sub-structure across different ethnic groups of Nepal based on mtDNA haplogroups and ectodysplasin-A receptor (EDAR) gene polymorphism. Our comprehensive high-resolution mtDNA-based genetic study of Tibeto-Burman communities reconstructs the maternal origins of prehistoric Himalayan populations and sheds light on migration events that have brought most of the East Eurasian ancestry to the present-day Nepalese population.

Transferability of the PRS estimates for height and BMI obtained from the European ethnic groups to the Western Russian populations

Genetic data plays an increasingly important role in modern medicine. Decrease in the cost of sequencing with subsequent increase in imputation accuracy, and the accumulation of large amounts of high-quality genetic data enable the creation of polygenic risk scores (PRSs) to perform genotype-phenotype associations. The accuracy of phenotype prediction primarily depends on the overall trait heritability, Genome-wide association studies cohort size, and the similarity of genetic background between the base and the target cohort. Here we utilized 8,664 high coverage genomic samples collected across Russia by "Evogen", a Russian biomedical company, to evaluate the predictive power of PRSs based on summary statistics established on cohorts of European ancestry for basic phenotypic traits, namely height and BMI. We have demonstrated that the PRSs calculated for selected traits in three distinct Russian populations, recapitulate the predictive power from the original studies. This is evidence that GWAS summary statistics calculated on cohorts of European ancestry are transferable onto at least some ethnic groups in Russia.

Genetic admixture and language shift in the medieval Volga-Oka interfluve

The Volga-Oka interfluve in northwestern Russia has an intriguing history of population influx and language shift during the Common Era. Today, most inhabitants of the region speak Russian, but until medieval times, northwestern Russia was inhabited by Uralic-speaking peoples.^1,2,3 A gradual shift to Slavic languages started in the second half of the first millennium with the expansion of Slavic tribes, which led to the foundation of the Kievan Rus' state in the late 9^th century CE. The medieval Rus' was multicultural and multilingual-historical records suggest that its northern regions comprised Slavic and Uralic peoples ruled by Scandinavian settlers.^4,5,6 In the 10^th-11^th centuries, the introduction of Christianity and Cyrillic literature raised the prestige status of Slavic, driving a language shift from Uralic to Slavic.³ This eventually led to the disappearance of the Uralic languages from northwestern Russia. Here, we study a 1,500-year time transect of 30 ancient genomes and stable isotope values from the Suzdal region in the Volga-Oka interfluve. We describe a previously unsampled local Iron Age population and a gradual genetic turnover in the following centuries. Our time transect captures the population shift associated with the spread of Slavic languages and illustrates the ethnically mixed state of medieval Suzdal principality, eventually leading to the formation of the admixed but fully Slavic-speaking population that inhabits the area today. We also observe genetic outliers that highlight the importance of the Suzdal region in medieval times as a hub of long-reaching contacts via trade and warfare.

North Asian population relationships in a global context

Population genetic studies of North Asian ethnic groups have focused on genetic variation of sex chromosomes and mitochondria. Studies of the extensive variation available from autosomal variation have appeared infrequently. We focus on relationships among population samples using new North Asia microhaplotype data. We combined genotypes from our laboratory on 58 microhaplotypes, distributed across 18 autosomes, on 3945 individuals from 75 populations with corresponding data extracted for 26 populations from the Thousand Genomes consortium and for 22 populations from the GenomeAsia 100 K project. A total of 7107 individuals in 122 total populations are analyzed using STRUCTURE, Principal Component Analysis, and phylogenetic tree analyses. North Asia populations sampled in Mongolia include: Buryats, Mongolians, Altai Kazakhs, and Tsaatans. Available Siberians include samples of Yakut, Khanty, and Komi Zyriane. Analyses of all 122 populations confirm many known relationships and show that most populations from North Asia form a cluster distinct from all other groups. Refinement of analyses on smaller subsets of populations reinforces the distinctiveness of North Asia and shows that the North Asia cluster identifies a region that is ancestral to Native Americans.

Agent-Based Modeling of Autosomal Recessive Deafness 1A (DFNB1A) Prevalence with Regard to Intensity of Selection Pressure in Isolated Human Population

An increase in the prevalence of autosomal recessive deafness 1A (DFNB1A) in populations of European descent was shown to be promoted by assortative marriages among deaf people. Assortative marriages became possible with the widespread introduction of sign language, resulting in increased genetic fitness of deaf individuals and, thereby, relaxing selection against deafness. However, the effect of this phenomenon was not previously studied in populations with different genetic structures. We developed an agent-based computer model for the analysis of the spread of DFNB1A. Using this model, we tested the impact of different intensities of selection pressure against deafness in an isolated human population over 400 years. Modeling of the "purifying" selection pressure on deafness ("No deaf mating" scenario) resulted in a decrease in the proportion of deaf individuals and the pathogenic allele frequency. Modeling of the "relaxed" selection ("Assortative mating" scenario) resulted in an increase in the proportion of deaf individuals in the first four generations, which then quickly plateaued with a subsequent decline and a decrease in the pathogenic allele frequency. The results of neutral selection pressure modeling ("Random mating" scenario) showed no significant changes in the proportion of deaf individuals or the pathogenic allele frequency after 400 years.

A common founder effect of the splice site variant c.-23 + 1G > A in GJB2 gene causing autosomal recessive deafness 1A (DFNB1A) in Eurasia

Mutations in the GJB2 gene are known to be a major cause of autosomal recessive deafness 1A (OMIM 220290). The most common pathogenic variants of the GJB2 gene have a high ethno-geographic specificity in their distribution, being attributed to a founder effect related to the Neolithic migration routes of Homo sapiens. The c.-23 + 1G > A splice site variant is frequently found among deaf patients of both Caucasian and Asian origins. It is currently unknown whether the spread of this mutation across Eurasia is a result of the founder effect or if it could have multiple local centers of origin. To determine the origin of c.-23 + 1G > A, we reconstructed haplotypes by genotyping SNPs on an Illumina OmniExpress 730 K platform of 23 deaf individuals homozygous for this variant from different populations of Eurasia. The analyses revealed the presence of common regions of homozygosity in different individual genomes in the sample. These data support the hypothesis of the common founder effect in the distribution of the c.-23 + 1G > A variant of the GJB2 gene. Based on the published data on the c.-23 + 1G > A prevalence among 16,177 deaf people and the calculation of the TMRCA of the modified f2-haplotypes carrying this variant, we reconstructed the potential migration routes of the carriers of this mutation around the world. This analysis indicates that the c.-23 + 1G > A variant in the GJB2 gene may have originated approximately 6000 years ago in the territory of the Caucasus or the Middle East then spread throughout Europe, South and Central Asia and other regions of the world.

Variation of Genomic Sites Associated with Severe Covid-19 Across Populations: Global and National Patterns

Background

Information about the distribution of clinically significant genetic markers in different populations may be helpful in elaborating personalized approaches to the clinical management of COVID-19 in the absence of consensus guidelines.

Aim

Analyze frequencies and distribution patterns of two markers associated with severe COVID-19 (rs11385942 and rs657152) and look for potential correlations between these markers and deaths from COVID-19 among populations in Russia and across the world.

Methods

We genotyped 1883 samples from 91 ethnic groups pooled into 28 populations representing Russia and its neighbor states. We also compiled a dataset on 32 populations from other regions using genotypes extracted or imputed from the available databases. Geographic maps showing the frequency distribution of the analyzed markers were constructed using the obtained data.

Results

The cartographic analysis revealed that rs11385942 distribution follows the West Eurasian pattern: the marker is frequent among the populations of Europe, West Asia and South Asia but rare or absent in all other parts of the globe. Notably, the transition from high to low rs11385942 frequencies across Eurasia is not abrupt but follows the clinal variation pattern instead. The distribution of rs657152 is more homogeneous. The analysis of correlations between the frequencies of the studied markers and the epidemiological characteristics of COVID-19 in a population revealed that higher frequencies of both risk alleles correlated positively with mortality from this disease. For rs657152, the correlation was especially strong (r = 0.59, p = 0.02). These reasonable correlations were observed for the "Russian" dataset only: no such correlations were established for the "world" dataset. This could be attributed to the differences in methodology used to collect COVID-19 statistics in different countries.

Conclusion

Our findings suggest that genetic differences between populations make a small yet tangible contribution to the heterogeneity of the pandemic worldwide.

Whole-Genome Sequencing of a 900-Year-Old Human Skeleton Supports Two Past Migration Events from the Russian Far East to Northern Japan

Recent studies on paleogenomics have reported some Paleolithic and Neolithic genomes that have provided new insights into the human population history in East and Northeast Asia. However, there remain some cases where more recent migration events need to be examined to elucidate the detailed formation process of local populations. Although the area around northern Japan is one of the regions archaeologically suggested to have been affected by migration waves after the Neolithic period, the genetic source of these migrations are still unclear. Thus, genomic data from such past migrant populations would be highly informative to clarify the detailed formation process of local populations in this region. Here, we report the genome sequence of a 900-year-old adult female (NAT002) belonging to the prehistoric Okhotsk people, who have been considered to be the past migrants to northern Japan after the Neolithic period. We found a close relationship between NAT002 and modern Lower Amur populations and past admixture events between the Amur, Jomon, and Kamchatka ancestries. The admixture dating suggested migration of Amur-related ancestry at approximately 1,600 BP, which is compatible with the archaeological evidence regarding the settlement of the Okhotsk people. Our results also imply migration of Kamchatka-related ancestry at approximately 2,000 BP. In addition, human leukocyte antigen (HLA) typing detected the HLA-B*40 allele, which is reported to increase the risk of arthritis, suggesting the genetic vulnerability of NAT002 to hyperostosis, which was observed around her chest clavicle.

Exploring correlations in genetic and cultural variation across language families in northeast Asia

Culture evolves in ways that are analogous to, but distinct from, genomes. Previous studies examined similarities between cultural variation and genetic variation (population history) at small scales within language families, but few studies have empirically investigated these parallels across language families using diverse cultural data. We report an analysis comparing culture and genomes from in and around northeast Asia spanning 11 language families. We extract and summarize the variation in language (grammar, phonology, lexicon), music (song structure, performance style), and genomes (genome-wide SNPs) and test for correlations. We find that grammatical structure correlates with population history (genetic history). Recent contact and shared descent fail to explain the signal, suggesting relationships that arose before the formation of current families. Our results suggest that grammar might be a cultural indicator of population history while also demonstrating differences among cultural and genetic relationships that highlight the complex nature of human history.

[N-aсetyltransferase (NAT2) gene polymorphism and gene network analysis]

To search for new targets of therapy, it is necessary to reconstruct the gene network of the disease, and identify the interaction of genes, proteins, and drug compounds. Using the online bioinformatics tools we have analyzed the current data set related to the metabolism of xenobiotics, mediated by the N-acetyltransferase 2 (NAT2) gene. The study of allelic polymorphism of the NAT2 gene has a prognostic value, allowing to determine the risk of a number of oncological diseases, the degree of increased risk due to smoking and exposure to chemical carcinogens, including drugs. The aim of this study was to determine the frequencies of two important "slow" variants of the NAT2 gene (NAT2*5, rs1801280 and NAT2*7, rs1799931), which significantly affected the rate of xenobiotic acetylation among the indigenous Nenets population of Northern Siberia. The obtained frequencies of polymorphic variants among the Nenets occupy an intermediate value between those for Europeans and Asians, which might indicate specific features of adaptation. We present a model of the distribution of two polymorphic variants of the NAT2 gene involved in the biotransformation of xenobiotics to study the characteristics of their metabolism in the indigenous inhabitants of Yamal.

Whole-Exome Sequencing Analysis of Human Semen Quality in Russian Multiethnic Population

The global trend toward the reduction of human spermatogenic function observed in many countries, including Russia, raised the problem of extensive screening and monitoring of male fertility and elucidation of its genetic and ethnic mechanisms. Recently, whole-exome sequencing (WES) was developed as a powerful tool for genetic analysis of complex traits. We present here the first Russian WES study for identification of new genes associated with semen quality. The experimental 3 × 2 design of the WES study was based on the analysis of 157 samples including three ethnic groups—Slavs (59), Buryats (n = 49), and Yakuts (n = 49), and two different semen quality groups—pathozoospermia (n = 95) and normospermia (n = 62). Additionally, our WES study group was negative for complete AZF microdeletions of the Y-chromosome. The normospermia group included men with normal sperm parameters in accordance with the WHO-recommended reference limit. The pathozoospermia group included men with impaired semen quality, namely, with any combined parameters of sperm concentration <15 × 10⁶/ml, and/or progressive motility <32%, and/or normal morphology <4%. The WES was performed for all 157 samples. Subsequent calling and filtering of variants were carried out according to the GATK Best Practices recommendations. On the genotyping stage, the samples were combined into four cohorts: three sets corresponded to three ethnic groups, and the fourth set contained all the 157 whole-exome samples. Association of the obtained polymorphisms with semen quality parameters was investigated using the χ2 test. To prioritize the obtained variants associated with pathozoospermia, their effects were determined using Ensembl Variant Effect Predictor. Moreover, polymorphisms located in genes expressed in the testis were revealed based on the genomic annotation. As a result, the nine potential SNP markers rs6971091, rs557806, rs610308, rs556052, rs1289658, rs278981, rs1129172, rs12268007, and rs17228441 were selected for subsequent verification on our previously collected population sample (about 1,500 males). The selected variants located in seven genes FAM71F1, PPP1R15A, TRIM45, PRAME, RBM47, WDFY4, and FSIP2 that are expressed in the testis and play an important role in cell proliferation, meiosis, and apoptosis.

Autosomal recessive cataract (CTRCT18) in the Yakut population isolate of Eastern Siberia: a novel founder variant in the FYCO1 gene

Congenital autosomal recessive cataract with unknown genetic etiology is one of the most common Mendelian diseases among the Turkic-speaking Yakut population (Eastern Siberia, Russia). To identify the genetic cause of congenital cataract spread in this population, we performed whole-exome sequencing (Illumina NextSeq 500) in one Yakut family with three affected siblings whose parents had preserved vision. We have revealed the novel homozygous c.1621C>T transition leading to premature stop codon p.(Gln541*) in exon 8 of the FYCO1 gene (NM_024513.4). Subsequent screening of c.1621C>T p.(Gln541*) revealed this variant in a homozygous state in 25 out of 29 Yakut families with congenital cataract (86%). Among 424 healthy individuals from seven populations of Eastern Siberia (Russians, Yakuts, Evenks, Evens, Dolgans, Chukchi, and Yukaghirs), the highest carrier frequency of c.1621C>T p.(Gln541*) was found in the Yakut population (7.9%). DNA samples of 25 homozygous for c.1621C>T p.(Gln541*) patients with congenital cataract and 114 unaffected unrelated individuals without this variant were used for a haplotype analysis based on the genotyping of six STR markers (D3S3512, D3S3685, D3S3582, D3S3561, D3S1289, and D3S3698). The structure of the identified haplotypes indicates a common origin for all of the studied mutant chromosomes bearing c.1621C>T p.(Gln541*). The age of the с.1621C>T p.(Gln541*) founder haplotype was estimated to be approximately 260 ± 65 years (10 generations). These findings characterize Eastern Siberia as the region of the world with the most extensive accumulation of the unique variant c.1621C>T p.(Gln541*) in the FYCO1 gene as a result of the founder effect.

Integrating Linguistics, Social Structure, and Geography to Model Genetic Diversity within India

India represents an intricate tapestry of population substructure shaped by geography, language, culture, and social stratification. Although geography closely correlates with genetic structure in other parts of the world, the strict endogamy imposed by the Indian caste system and the large number of spoken languages add further levels of complexity to understand Indian population structure. To date, no study has attempted to model and evaluate how these factors have interacted to shape the patterns of genetic diversity within India. We merged all publicly available data from the Indian subcontinent into a data set of 891 individuals from 90 well-defined groups. Bringing together geography, genetics, and demographic factors, we developed Correlation Optimization of Genetics and Geodemographics to build a model that explains the observed population genetic substructure. We show that shared language along with social structure have been the most powerful forces in creating paths of gene flow in the subcontinent. Furthermore, we discover the ethnic groups that best capture the diverse genetic substructure using a ridge leverage score statistic. Integrating data from India with a data set of additional 1,323 individuals from 50 Eurasian populations, we find that Indo-European and Dravidian speakers of India show shared genetic drift with Europeans, whereas the Tibeto-Burman speaking tribal groups have maximum shared genetic drift with East Asians.

Dissecting dynamics and differences of selective pressures in the evolution of human pigmentation

Human pigmentation is a highly diverse and complex trait among populations and has drawn particular attention from both academic and non-academic investigators for thousands of years. Previous studies detected selection signals in several human pigmentation genes, but few studies have integrated contribution from multiple genes to the evolution of human pigmentation. Moreover, none has quantified selective pressures on human pigmentation over epochs and between populations. Here, we dissect dynamics and differences of selective pressures during different periods and between distinct populations with new approaches. We use genotype data of 19 genes associated with human pigmentation from 17 publicly available datasets and obtain data for 2346 individuals of six representative population groups from across the world. Our results quantify the strength of natural selection on light pigmentation not only in modern Europeans (0.0259/generation) but also in proto-Eurasians (0.00650/generation). Our results also suggest that several derived alleles associated with human dark pigmentation may be under positive directional selection in some African populations. Our study provides the first attempt to quantitatively investigate the dynamics of selective pressures during different time periods in the evolution of human pigmentation.This article has an associated First Person interview with the first author of the article.

Mitochondrial genome diversity on the Central Siberian Plateau with particular reference to the prehistory of northernmost Eurasia

The Central Siberian Plateau was the last geographic area in Eurasia to become habitable by modern humans after the Last Glacial Maximum (LGM). Through a comprehensive dataset of mitochondrial DNA (mtDNA) genomes retained in the remnats of earlier ("Old") Siberians, primarily the Ket, Tofalar, and Todzhi, we explored genetic links between the Yenisei-Sayan region and Northeast Eurasia (best represented by the Yukaghir) over the last 10,000 years. We generated 218 new complete mtDNA sequences and placed them into compound phylogenies with 7 newly obtained and 70 published ancient mitochondrial genomes. We have considerably extended the mtDNA sequence diversity (at the entire mtDNA genome level) of autochthonous Siberians, which remain poorly sampled, and these new data may have a broad impact on the study of human migration. We compared present-day mtDNA diversity in these groups with complete mitochondrial genomes from ancient samples from the region and placed the samples into combined genealogical trees. The resulting components were used to clarify the origins and expansion history of mtDNA lineages that evolved in the refugia of south-central Siberia and beyond, as well as multiple phases of connection between this region and distant parts of Eurasia.

Cartographic atlas of frequency variation for 45 pharmacogenetic markers in populations of Russia and its neighbor states

The lack of information about the frequency of pharmacogenetic markers in Russia impedes the adoption of personalized treatment algorithms originally developed for West European populations. The aim of this paper was to study the distribution of some clinically significant pharmacogenetic markers across Russia. A total of 45 pharmacogenetic markers were selected from a few population genetic datasets, including ADME, drug target and hemostasis-controlling genes. The total number of donors genotyped for these markers was 2,197. The frequencies of these markers were determined for 50 different populations, comprised of 137 ethnic and subethnic groups. A comprehensive pharmacogenetic atlas was created, i.e. a systematic collection of gene geographic maps of frequency variation for 45 pharmacogenetic DNA markers in Russia and its neighbor states. The maps revealed 3 patterns of geographic variation. Clinal variation (a gradient change in frequency along the East-West axis) is observed in the pharmacogenetic markers that follow the main pattern of variation for North Eurasia (13% of the maps). Uniform distribution singles out a group of markers that occur at average frequency in most Russian regions (27% of the maps). Focal variation is observed in the markers that are specific to a certain group of populations and are absent in other regions (60% of the maps). The atlas reveals that the average frequency of the marker and its frequency in individual populations do not indicate the type of its distribution in Russia: a gene geographic map is needed to uncover the pattern of its variation.

Gut microbiome signature of Viliuisk encephalomyelitis in Yakuts includes an increase in microbes linked to lean body mass and eating behaviour

BACKGROUND

Viliuisk encephalomyelitis (VE) is a rare endemic neurodegenerative disease occurring in the Yakut population of Northeastern Siberia. The main clinical features of VE are spasticity, dysarthria, dementia, central paresis and paralysis, and cortical atrophy observed via MRI. Many hypotheses have been proposed regarding its etiology, including infectious agents, genetics, environmental factors, and immunopathology. Each of these hypotheses has been supported to some extent by epidemiological and experimental data. Nevertheless, none of them has been decisively proven. Gut microbiome is one of the factors that might be involved in VE pathogenesis.

RESULTS

Here we performed a pilot survey of the stool microbiomes of Yakut subjects with VE (n = 6) and without VE (n = 11). 16S rRNA sequencing showed that in comparison with the control group, the Yakuts with VE had increased proportions of Methanobrevibacter and Christensenella, which are reported to be linked to body mass index, metabolism, dietary habits and potentially to neurodegenerative disorders. The identified associations suggest that the microbiome may be involved in VE. Overall, the Yakut microbiome was quite specific in comparison with other populations, such as metropolitan Russians and native inhabitants of the Canadian Arctic.

CONCLUSIONS

Describing the gut microbiome of indigenous human populations will help to elucidate the impact of dietary and environmental factors on microbial community structure and identify risks linked to the lifestyles of such groups as well as endemic diseases.

The genetic legacy of legendary and historical Siberian chieftains

Seventeen years of archaeological and anthropological expeditions in North-Eastern Siberia (in the Sakha Republic, Yakutia) have permitted the genetic analysis of 150 ancient (15th-19th century) and 510 modern individuals. Almost all males were successfully analysed (Y-STR) and this allowed us to identify paternal lineages and their geographical expansion through time. This genetic data was confronted with mythological, historical and material evidence to establish the sequence of events that built the modern Yakut genetic diversity. We show that the ancient Yakuts recovered from this large collection of graves are not representative of an ancient population. Uncommonly, we were also able to demonstrate that the funerary preference observed here involved three specific male lineages, especially in the 18^th century. Moreover, this dominance was likely caused by the Russian conquest of Siberia which allowed some male clans to rise to new levels of power. Finally, we give indications that some mythical and historical figures might have been the actors of those genetic changes. These results help us reconsider the genetic dynamics of colonization in some regions, question the distinction between fact and myth in national histories and provide a rare insight into a funerary ensemble by revealing the biased process of its composition.

High Rates of Three Common GJB2 Mutations c.516G>C, c.-23+1G>A, c.235delC in Deaf Patients from Southern Siberia Are Due to the Founder Effect

The mutations in the GJB2 gene (13q12.11, MIM 121011) encoding transmembrane protein connexin 26 (Cx26) account for a significant portion of hereditary hearing loss worldwide. Earlier we found a high prevalence of recessive GJB2 mutations c.516G>C, c.-23+1G>A, c.235delC in indigenous Turkic-speaking Siberian peoples (Tuvinians and Altaians) from the Tyva Republic and Altai Republic (Southern Siberia, Russia) and proposed the founder effect as a cause for their high rates in these populations. To reconstruct the haplotypes associated with each of these mutations, the genotyping of polymorphic genetic markers both within and flanking the GJB2 gene was performed in 28 unrelated individuals homozygous for c.516G>C (n = 18), c.-23+1G>A (n = 6), or c.235delC (n = 4) as well as in the ethnically matched controls (62 Tuvinians and 55 Altaians) without these mutations. The common haplotypes specific for mutations c.516G>C, c.-23+1G>A, or c.235delC were revealed implying a single origin of each of these mutations. The age of mutations estimated by the DMLE+ v2.3 software and the single marker method is discussed in relation to ethnic history of Tuvinians and Altaians. The data obtained in this study support a crucial role of the founder effect in the high prevalence of GJB2 mutations c.516G>C, c.-23+1G>A, c.235delC in indigenous populations of Southern Siberia.

Genetic Landscapes Reveal How Human Genetic Diversity Aligns with Geography

Geographic patterns in human genetic diversity carry footprints of population history and provide insights for genetic medicine and its application across human populations. Summarizing and visually representing these patterns of diversity has been a persistent goal for human geneticists, and has revealed that genetic differentiation is frequently correlated with geographic distance. However, most analytical methods to represent population structure do not incorporate geography directly, and it must be considered post hoc alongside a visual summary of the genetic structure. Here, we estimate "effective migration" surfaces to visualize how human genetic diversity is geographically structured. The results reveal local patterns of differentiation in detail and emphasize that while genetic similarity generally decays with geographic distance, the relationship is often subtly distorted. Overall, the visualizations provide a new perspective on genetics and geography in humans and insight to the geographic distribution of human genetic variation.

PGG.SNV: understanding the evolutionary and medical implications of human single nucleotide variations in diverse populations

Despite the tremendous growth of the DNA sequencing data in the last decade, our understanding of the human genome is still in its infancy. To understand the implications of genetic variants in the light of population genetics and molecular evolution, we developed a database, PGG.SNV ( https://www.pggsnv.org ), which gives much higher weight to previously under-investigated indigenous populations in Asia. PGG.SNV archives 265 million SNVs across 220,147 present-day genomes and 1018 ancient genomes, including 1009 newly sequenced genomes, representing 977 global populations. Moreover, estimation of population genetic diversity and evolutionary parameters is available in PGG.SNV, a unique feature compared with other databases.

PGG.Population: a database for understanding the genomic diversity and genetic ancestry of human populations

There are a growing number of studies focusing on delineating genetic variations that are associated with complex human traits and diseases due to recent advances in next-generation sequencing technologies. However, identifying and prioritizing disease-associated causal variants relies on understanding the distribution of genetic variations within and among populations. The PGG.Population database documents 7122 genomes representing 356 global populations from 107 countries and provides essential information for researchers to understand human genomic diversity and genetic ancestry. These data and information can facilitate the design of research studies and the interpretation of results of both evolutionary and medical studies involving human populations. The database is carefully maintained and constantly updated when new data are available. We included miscellaneous functions and a user-friendly graphical interface for visualization of genomic diversity, population relationships (genetic affinity), ancestral makeup, footprints of natural selection, and population history etc. Moreover, PGG.Population provides a useful feature for users to analyze data and visualize results in a dynamic style via online illustration. The long-term ambition of the PGG.Population, together with the joint efforts from other researchers who contribute their data to our database, is to create a comprehensive depository of geographic and ethnic variation of human genome, as well as a platform bringing influence on future practitioners of medicine and clinical investigators. PGG.Population is available at https://www.pggpopulation.org.

Shifts in the Genetic Landscape of the Western Eurasian Steppe Associated with the Beginning and End of the Scythian Dominance

The Early Iron Age nomadic Scythians have been described as a confederation of tribes of different origins, based on ancient DNA evidence [1-3]. It is still unclear how much of the Scythian dominance in the Eurasian Steppe was due to movements of people and how much reflected cultural diffusion and elite dominance. We present new whole-genome sequences of 31 ancient Western and Eastern Steppe individuals, including Scythians as well as samples pre- and postdating them, allowing us to set the Scythians in a temporal context (in the Western, i.e., Ponto-Caspian Steppe). We detect an increase of eastern (Altaian) affinity along with a decrease in eastern hunter-gatherer (EHG) ancestry in the Early Iron Age Ponto-Caspian gene pool at the start of the Scythian dominance. On the other hand, samples of the Chernyakhiv culture postdating the Scythians in Ukraine have a significantly higher proportion of Near Eastern ancestry than other samples of this study. Our results agree with the Gothic source of the Chernyakhiv culture and support the hypothesis that the Scythian dominance did involve a demic component.

The population history of northeastern Siberia since the Pleistocene

Northeastern Siberia has been inhabited by humans for more than 40,000 years but its deep population history remains poorly understood. Here we investigate the late Pleistocene population history of northeastern Siberia through analyses of 34 newly recovered ancient genomes that date to between 31,000 and 600 years ago. We document complex population dynamics during this period, including at least three major migration events: an initial peopling by a previously unknown Palaeolithic population of 'Ancient North Siberians' who are distantly related to early West Eurasian hunter-gatherers; the arrival of East Asian-related peoples, which gave rise to 'Ancient Palaeo-Siberians' who are closely related to contemporary communities from far-northeastern Siberia (such as the Koryaks), as well as Native Americans; and a Holocene migration of other East Asian-related peoples, who we name 'Neo-Siberians', and from whom many contemporary Siberians are descended. Each of these population expansions largely replaced the earlier inhabitants, and ultimately generated the mosaic genetic make-up of contemporary peoples who inhabit a vast area across northern Eurasia and the Americas.

The Arrival of Siberian Ancestry Connecting the Eastern Baltic to Uralic Speakers further East

In this study, we compare the genetic ancestry of individuals from two as yet genetically unstudied cultural traditions in Estonia in the context of available modern and ancient datasets: 15 from the Late Bronze Age stone-cist graves (1200-400 BC) (EstBA) and 6 from the Pre-Roman Iron Age tarand cemeteries (800/500 BC-50 AD) (EstIA). We also included 5 Pre-Roman to Roman Iron Age Ingrian (500 BC-450 AD) (IngIA) and 7 Middle Age Estonian (1200-1600 AD) (EstMA) individuals to build a dataset for studying the demographic history of the northern parts of the Eastern Baltic from the earliest layer of Mesolithic to modern times. Our findings are consistent with EstBA receiving gene flow from regions with strong Western hunter-gatherer (WHG) affinities and EstIA from populations related to modern Siberians. The latter inference is in accordance with Y chromosome (chrY) distributions in present day populations of the Eastern Baltic, as well as patterns of autosomal variation in the majority of the westernmost Uralic speakers [1-5]. This ancestry reached the coasts of the Baltic Sea no later than the mid-first millennium BC; i.e., in the same time window as the diversification of west Uralic (Finnic) languages [6]. Furthermore, phenotypic traits often associated with modern Northern Europeans, like light eyes, hair, and skin, as well as lactose tolerance, can be traced back to the Bronze Age in the Eastern Baltic. VIDEO ABSTRACT.

Genome-wide sequence analyses of ethnic populations across Russia

The Russian Federation is the largest and one of the most ethnically diverse countries in the world, however no centralized reference database of genetic variation exists to date. Such data are crucial for medical genetics and essential for studying population history. The Genome Russia Project aims at filling this gap by performing whole genome sequencing and analysis of peoples of the Russian Federation. Here we report the characterization of genome-wide variation of 264 healthy adults, including 60 newly sequenced samples. People of Russia carry known and novel genetic variants of adaptive, clinical and functional consequence that in many cases show allele frequency divergence from neighboring populations. Population genetics analyses revealed six phylogeographic partitions among indigenous ethnicities corresponding to their geographic locales. This study presents a characterization of population-specific genomic variation in Russia with results important for medical genetics and for understanding the dynamic population history of the world's largest country.

Medieval mummies of Zeleny Yar burial ground in the Arctic Zone of Western Siberia

Notwithstanding the pioneering achievements of studies on arctic mummies in Siberia, there are insufficient data for any comprehensive understanding of the bio-cultural details of medieval people living in the region. In the Western Siberian arctic, permafrost mummies have been found in 12th to 13th century graves located in the Zeleny Yar (Z-Y) burial ground (66°19'4.54"С; 67°21'13.54"В). In 2013-2016, we were fortunate to be able to excavate that cemetery, locating a total of 47 burials, including cases of mummification. Some of these mummies had been wrapped in a multi-layered birch-bark cocoon. After removal of the cocoon, we conducted interdisciplinary studies using various scientific techniques. Gross anatomical examination and CT radiography showed that the internal organs were still well preserved inside the body cavities. Under light and electron microscopy, the histological findings were very similar to those for naturally mummified specimens discovered in other countries. Ancient DNA analysis showed that the Z-Y mummies' mtDNA haplotypes belong to five different haplogroups, namely U5a (#34), H3ao (#53), D (#67-1), U4b1b1 (#67-2), and D4j8 (#68), which distinguish them for their unique combination of Western- and Eastern Siberia-specific mtDNA haplogroups. Our interdisciplinary study obtained fundamental information that will form the foundation of successful future investigations on medieval mummies found in the Western Siberian arctic.

The Genetic Ancestry of Modern Indus Valley Populations from Northwest India

The Indus Valley has been the backdrop for several historic and prehistoric population movements between South Asia and West Eurasia. However, the genetic structure of present-day populations from Northwest India is poorly characterized. Here we report new genome-wide genotype data for 45 modern individuals from four Northwest Indian populations, including the Ror, whose long-term occupation of the region can be traced back to the early Vedic scriptures. Our results suggest that although the genetic architecture of most Northwest Indian populations fits well on the broader North-South Indian genetic cline, culturally distinct groups such as the Ror stand out by being genetically more akin to populations living west of India; such populations include prehistorical and early historical ancient individuals from the Swat Valley near the Indus Valley. We argue that this affinity is more likely a result of genetic continuity since the Bronze Age migrations from the Steppe Belt than a result of recent admixture. The observed patterns of genetic relationships both with modern and ancient West Eurasians suggest that the Ror can be used as a proxy for a population descended from the Ancestral North Indian (ANI) population. Collectively, our results show that the Indus Valley populations are characterized by considerable genetic heterogeneity that has persisted over thousands of years.

Siberian genetic diversity reveals complex origins of the Samoyedic-speaking populations

OBJECTIVES

We examined autosomal genome-wide SNPs and Y-chromosome data from 15 Siberian and 12 reference populations to study the affinities of Siberian populations, and to address hypotheses about the origin of the Samoyed peoples.

METHODS

Samples were genotyped for 567 096 autosomal SNPs and 147 Y-chromosome polymorphic sites. For several analyses, we used 281 093 SNPs from the intersection of our data with publicly available ancient Siberian samples. To examine genetic relatedness among populations, we applied PCA, F_ST , TreeMix, and ADMIXTURE analyses. To explore the potential effect of demography and evolutionary processes, the distribution of ROH and IBD sharing within population were studied.

RESULTS

Analyses of autosomal and Y-chromosome data reveal high differentiation of the Siberian groups. The Siberian populations have a large proportion of their genome in ROH and IBD segments. Several populations (ie, Nganasans, Evenks, Yukagirs, and Koryaks) do not appear to have experienced admixture with other Siberian populations (ie, producing only positive f3), while for the other tested populations the composition of mixing sources always included Nganasans or Evenks. The Nganasans from the Taymyr Peninsula demonstrate the greatest level of shared shorter ROH and IBD with nearly all other Siberian populations.

CONCLUSIONS

Autosomal SNP and Y-chromosome data demonstrate that Samoyedic populations differ significantly in their genetic composition. Genetic relationship is observed only between Forest and Tundra Nentsi. Selkups are affiliated with the Kets from the Yenisey River, while the Nganasans are separated from their linguistic neighbors, showing closer affinities with the Evenks and Yukagirs.

Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya origin in the conquering Hungarians

It has been widely accepted that the Finno-Ugric Hungarian language, originated from proto Uralic people, was brought into the Carpathian Basin by the conquering Hungarians. From the middle of the 19th century this view prevailed against the deep-rooted Hungarian Hun tradition, maintained in folk memory as well as in Hungarian and foreign written medieval sources, which claimed that Hungarians were kinsfolk of the Huns. In order to shed light on the genetic origin of the Conquerors we sequenced 102 mitogenomes from early Conqueror cemeteries and compared them to sequences of all available databases. We applied novel population genetic algorithms, named Shared Haplogroup Distance and MITOMIX, to reveal past admixture of maternal lineages. Our results show that the Conquerors assembled from various nomadic groups of the Eurasian steppe. Population genetic results indicate that they had closest connection to the Onogur-Bulgar ancestors of Volga Tatars. Phylogenetic results reveal that more than one third of the Conqueror maternal lineages were derived from Central-Inner Asia and their most probable ultimate sources were the Asian Scythians and Asian Huns, giving support to the Hungarian Hun tradition. The rest of the lineages most likely originated from the Bronze Age Potapovka-Poltavka-Srubnaya cultures of the Pontic-Caspian steppe. Available data imply that the Conquerors did not have a major contribution to the gene pool of the Carpathian Basin.

Genes reveal traces of common recent demographic history for most of the Uralic-speaking populations

BACKGROUND

The genetic origins of Uralic speakers from across a vast territory in the temperate zone of North Eurasia have remained elusive. Previous studies have shown contrasting proportions of Eastern and Western Eurasian ancestry in their mitochondrial and Y chromosomal gene pools. While the maternal lineages reflect by and large the geographic background of a given Uralic-speaking population, the frequency of Y chromosomes of Eastern Eurasian origin is distinctively high among European Uralic speakers. The autosomal variation of Uralic speakers, however, has not yet been studied comprehensively.

RESULTS

Here, we present a genome-wide analysis of 15 Uralic-speaking populations which cover all main groups of the linguistic family. We show that contemporary Uralic speakers are genetically very similar to their local geographical neighbours. However, when studying relationships among geographically distant populations, we find that most of the Uralic speakers and some of their neighbours share a genetic component of possibly Siberian origin. Additionally, we show that most Uralic speakers share significantly more genomic segments identity-by-descent with each other than with geographically equidistant speakers of other languages. We find that correlated genome-wide genetic and lexical distances among Uralic speakers suggest co-dispersion of genes and languages. Yet, we do not find long-range genetic ties between Estonians and Hungarians with their linguistic sisters that would distinguish them from their non-Uralic-speaking neighbours.

CONCLUSIONS

We show that most Uralic speakers share a distinct ancestry component of likely Siberian origin, which suggests that the spread of Uralic languages involved at least some demic component.

Ancient human parallel lineages within North America contributed to a coastal expansion

Little is known regarding the first people to enter the Americas and their genetic legacy. Genomic analysis of the oldest human remains from the Americas showed a direct relationship between a Clovis-related ancestral population and all modern Central and South Americans as well as a deep split separating them from North Americans in Canada. We present 91 ancient human genomes from California and Southwestern Ontario and demonstrate the existence of two distinct ancestries in North America, which possibly split south of the ice sheets. A contribution from both of these ancestral populations is found in all modern Central and South Americans. The proportions of these two ancestries in ancient and modern populations are consistent with a coastal dispersal and multiple admixture events.

The first horse herders and the impact of early Bronze Age steppe expansions into Asia

The Yamnaya expansions from the western steppe into Europe and Asia during the Early Bronze Age (~3000 BCE) are believed to have brought with them Indo-European languages and possibly horse husbandry. We analyzed 74 ancient whole-genome sequences from across Inner Asia and Anatolia and show that the Botai people associated with the earliest horse husbandry derived from a hunter-gatherer population deeply diverged from the Yamnaya. Our results also suggest distinct migrations bringing West Eurasian ancestry into South Asia before and after, but not at the time of, Yamnaya culture. We find no evidence of steppe ancestry in Bronze Age Anatolia from when Indo-European languages are attested there. Thus, in contrast to Europe, Early Bronze Age Yamnaya-related migrations had limited direct genetic impact in Asia.

Investigating Holocene human population history in North Asia using ancient mitogenomes

Archaeogenomic studies have largely elucidated human population history in West Eurasia during the Stone Age. However, despite being a broad geographical region of significant cultural and linguistic diversity, little is known about the population history in North Asia. We present complete mitochondrial genome sequences together with stable isotope data for 41 serially sampled ancient individuals from North Asia, dated between c.13,790 BP and c.1,380 BP extending from the Palaeolithic to the Iron Age. Analyses of mitochondrial DNA sequences and haplogroup data of these individuals revealed the highest genetic affinity to present-day North Asian populations of the same geographical region suggesting a possible long-term maternal genetic continuity in the region. We observed a decrease in genetic diversity over time and a reduction of maternal effective population size (N_e) approximately seven thousand years before present. Coalescent simulations were consistent with genetic continuity between present day individuals and individuals dating to 7,000 BP, 4,800 BP or 3,000 BP. Meanwhile, genetic differences observed between 7,000 BP and 3,000 BP as well as between 4,800 BP and 3,000 BP were inconsistent with genetic drift alone, suggesting gene flow into the region from distant gene pools or structure within the population. These results indicate that despite some level of continuity between ancient groups and present-day populations, the region exhibits a complex demographic history during the Holocene.

Exome Sequencing Provides Evidence of Polygenic Adaptation to a Fat-Rich Animal Diet in Indigenous Siberian Populations

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.

Reconstructing the demographic history of the Himalayan and adjoining populations

The rugged topography of the Himalayan region has hindered large-scale human migrations, population admixture and assimilation. Such complexity in geographical structure might have facilitated the existence of several small isolated communities in this region. We have genotyped about 850,000 autosomal markers among 35 individuals belonging to the four major populations inhabiting the Himalaya and adjoining regions. In addition, we have genotyped 794 individuals belonging to 16 ethnic groups from the same region, for uniparental (mitochondrial and Y chromosomal DNA) markers. Our results in the light of various statistical analyses suggest a closer link of the Himalayan and adjoining populations to East Asia than their immediate geographical neighbours in South Asia. Allele frequency-based analyses likely support the existence of a specific ancestry component in the Himalayan and adjoining populations. The admixture time estimate suggests a recent westward migration of populations living to the East of the Himalaya. Furthermore, the uniparental marker analysis among the Himalayan and adjoining populations reveal the presence of East, Southeast and South Asian genetic signatures. Interestingly, we observed an antagonistic association of Y chromosomal haplogroups O3 and D clines with the longitudinal distance. Thus, we summarise that studying the Himalayan and adjoining populations is essential for a comprehensive reconstruction of the human evolutionary and ethnolinguistic history of eastern Eurasia.

Between Lake Baikal and the Baltic Sea: genomic history of the gateway to Europe

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Population specific analysis of Yakut exomes

We studied the genetic diversity of the Yakut population using exome sequencing. We performed comparative analysis of the Yakut population and the populations that are included in the "1000 Genomes" project and we identified the alleles specific to the Yakut population. We showed, that the Yakuts population is a separate cluster between Europeans and East Asians.

An earlier revolution: genetic and genomic analyses reveal pre-existing cultural differences leading to Neolithization

Archaeological evidence shows that, in the long run, Neolitization (the transition from foraging to food production) was associated with demographic growth. We used two methods (patterns of linkage disequilibrium from whole-genome SNPs and MSMC estimates on genomes) to reconstruct the demographic profiles for respectively 64 and 24 modern-day populations with contrasting lifestyles across the Old World (sub-Saharan Africa, south-eastern Asia, Siberia). Surprisingly, in all regions, food producers had larger effective population sizes (N_e) than foragers already 20 k years ago, well before the Neolithic revolution. As expected, this difference further increased ~12–10 k years ago, around or just before the onset of food production. Using paleoclimate reconstructions, we show that the early difference in N_e cannot be explained by food producers inhabiting more favorable regions. A number of mechanisms, including ancestral differences in census size, sedentism, exploitation of the natural resources, social stratification or connectivity between groups, might have led to the early differences in Ne detected in our analyses. Irrespective of the specific mechanisms involved, our results provide further evidence that long term cultural differences among populations of Palaeolithic hunter-gatherers are likely to have played an important role in the later Neolithization process.

Human ancestry correlates with language and reveals that race is not an objective genomic classifier

Genetic and archaeological studies have established a sub-Saharan African origin for anatomically modern humans with subsequent migrations out of Africa. Using the largest multi-locus data set known to date, we investigated genetic differentiation of early modern humans, human admixture and migration events, and relationships among ancestries and language groups. We compiled publicly available genome-wide genotype data on 5,966 individuals from 282 global samples, representing 30 primary language families. The best evidence supports 21 ancestries that delineate genetic structure of present-day human populations. Independent of self-identified ethno-linguistic labels, the vast majority (97.3%) of individuals have mixed ancestry, with evidence of multiple ancestries in 96.8% of samples and on all continents. The data indicate that continents, ethno-linguistic groups, races, ethnicities, and individuals all show substantial ancestral heterogeneity. We estimated correlation coefficients ranging from 0.522 to 0.962 between ancestries and language families or branches. Ancestry data support the grouping of Kwadi-Khoe, Kx’a, and Tuu languages, support the exclusion of Omotic languages from the Afroasiatic language family, and do not support the proposed Dené-Yeniseian language family as a genetically valid grouping. Ancestry data yield insight into a deeper past than linguistic data can, while linguistic data provide clarity to ancestry data.

Reconstructing genetic history of Siberian and Northeastern European populations

Siberia and Northwestern Russia are home to over 40 culturally and linguistically diverse indigenous ethnic groups, yet genetic variation and histories of peoples from this region are largely uncharacterized. We present deep whole-genome sequencing data (∼38×) from 28 individuals belonging to 14 distinct indigenous populations from that region. We combined these data sets with additional 32 modern-day and 46 ancient human genomes to reconstruct genetic histories of several indigenous Northern Eurasian populations. We found that Siberian and East Asian populations shared 38% of their ancestry with a 45,000-yr-old Ust'-Ishim individual who was previously believed to have no modern-day descendants. Western Siberians trace 57% of their ancestry to ancient North Eurasians, represented by the 24,000-yr-old Siberian Mal'ta boy MA-1. Eastern Siberian populations formed a distinct sublineage that separated from other East Asian populations ∼10,000 yr ago. In addition, we uncovered admixtures between Siberians and Eastern European hunter-gatherers from Samara, Karelia, Hungary, and Sweden (from 8000-6600 yr ago); Yamnaya people (5300-4700 yr ago); and modern-day Northeastern Europeans. Our results provide new insights into genetic histories of Siberian and Northeastern European populations and evidence of ancient gene flow from Siberia into Europe.

Human Y Chromosome Haplogroup N: A Non-trivial Time-Resolved Phylogeography that Cuts across Language Families

The paternal haplogroup (hg) N is distributed from southeast Asia to eastern Europe. The demographic processes that have shaped the vast extent of this major Y chromosome lineage across numerous linguistically and autosomally divergent populations have previously been unresolved. On the basis of 94 high-coverage re-sequenced Y chromosomes, we establish and date a detailed hg N phylogeny. We evaluate geographic structure by using 16 distinguishing binary markers in 1,631 hg N Y chromosomes from a collection of 6,521 samples from 56 populations. The more southerly distributed sub-clade N4 emerged before N2a1 and N3, found mostly in the north, but the latter two display more elaborate branching patterns, indicative of regional contrasts in recent expansions. In particular, a number of prominent and well-defined clades with common N3a3'6 ancestry occur in regionally dissimilar northern Eurasian populations, indicating almost simultaneous regional diversification and expansion within the last 5,000 years. This patrilineal genetic affinity is decoupled from the associated higher degree of language diversity.

The frequency of SLCO1B1*5 polymorphism genotypes among Russian and Sakha (Yakutia) patients with hypercholesterolemia

Introduction

Statins are the most commonly prescribed medicines for treatment of hypercholesterolemia. At the same time, up to 25% of patients cannot tolerate or have to discontinue the statin therapy due to statin-induced myopathy. In a majority of cases, statin-induced myopathy is attributed to SLCO1B1 gene polymorphism. The strongest association between statin-induced myopathy and SLCO1B1 gene polymorphism was described for simvastatin. Our research was focused on the frequency of SLCO1B1*5 genetic variant in the Russian population and in the native population of Sakha (Yakutia).

Materials and methods

A total of 1,071 hyperlipidemic Russian and 76 hyperlipidemic Sakha (Yakutian) patients were included in the study. Genotypes of SLCO1B1*5 (c.521T>C, rs4149056) were determined with polymerase chain reaction amplification. The results of our study were compared with data about hyperlipidemic patients in available publications.

Results

In the Russian population 665 (62%) patients had TT genotype of SLCO1B1*5, 346 (32%) patients had TC genotype, and in 60 patients (6%) CC variant was found (Hardy–Weinberg’s chi-square test was 3.1 P=0.21). In comparison with Brazil, France, the People’s Republic of China, Japan, and the native population of Sakha (Yakutia), C-allele, which causes an increased risk of statin-induced myopathy, was found significantly more often in the Russian population. In the native population of Sakha (Yakutia) SLCO1B1 polymorphism was TT – 62 (82%), TC – 11 (14%), CC – 3 (4%) (Hardy–Weinberg’s chi-square test was 5.13 P=0.077). In comparison with data from Brazil, France, the People’s Republic of China, and Japan, C-allele frequency in the Sakha (Yakutian) population was not significantly different.

Conclusion

Thus, we have studied the incidence of pathologic SLCO1B1 c.521C-allele in Russian and Sakha hyperlipidemic patients. The presence of SLCO1B1 C-allele in patients with hyperlipidemia forces us to be more careful in hypolipidemic drug prescription, especially statins, according to a higher risk of statin-induced myopathy development. The fact that SLCO1B1 C-allele is rarer among Sakha patients, could be interesting from the point of studying adverse drug effects frequency and statins’ effectiveness.