A Machine-Learning-Based Approach to Prediction of Biogeographic Ancestry within Europe.

Data obtained with the use of massive parallel sequencing (MPS) can be valuable in population genetics studies. In particular, such data harbor the potential for distinguishing samples from different populations, especially from those coming from adjacent populations of common origin. Machine learning (ML) techniques seem to be especially well suited for analyzing large datasets obtained using MPS. The Slavic populations constitute about a third of the population of Europe and inhabit a large area of the continent, while being relatively closely related in population genetics terms. In this proof-of-concept study, various ML techniques were used to classify DNA samples from Slavic and non-Slavic individuals. The primary objective of this study was to empirically evaluate the feasibility of discerning the genetic provenance of individuals of Slavic descent who exhibit genetic similarity, with the overarching goal of categorizing DNA specimens derived from diverse Slavic population representatives. Raw sequencing data were pre-processed, to obtain a 1200 character-long binary vector. A total of three classifiers were used-Random Forest, Support Vector Machine (SVM), and XGBoost. The most-promising results were obtained using SVM with a linear kernel, with 99.9% accuracy and F1-scores of 0.9846-1.000 for all classes.

An association between copy number variation of enhancer involved in craniofacial development and biogeographic ancestry.

Human facial morphology is a combination of many complex traits and is determined by a large number of genes and enhancers. Here, we report a Copy Number Variation (CNV) study of enhancer hs1431 in populations of Central European and South Siberian ancestry. Central European samples included 97 Poles, while South Siberian samples included 78 Buryats and 27 Tuvinians. CNVs were detected by real-time PCR, using ViiA™ 7 Real-Time PCR System (Applied Biosystems). We revealed significant differences in CNV of hs1431 enhancer between Polish and Buryat population (p=0.0378), but not between Central European and South Siberian population (p=0.1225). Our results suggest that an increase in copy number variation of hs1431 enhancer is associated with biogeographic ancestry. However, this result needs extending and replicating in larger cohorts. This is the first study revealing the presence of copy number variation of enhancer hs1431 in humans.

Verification of insertion-deletion markers (InDels) and microsatellites (STRs) as subsidiary tools for inferring Slavic population ancestry.

Genetic markers for the prediction of biogeographical ancestry have proved to be effective tools for law enforcement agencies for many years now. In this study, we attempted to assess the potential of insertion-deletion markers (InDel) and microsatellites (STRs) as subsidiary polymorphisms for inference of Slavic population ancestry. For that purpose, we genotyped Slavic-speaking populations samples from Belarus, the Czech Republic, Poland, Serbia, Ukraine and Russia in 46 InDels and 15 STRs by PCR and capillary electrophoresis and analyzed for between-population differentiation with the use of distance-based methods (FST, principal component analysis and multidimensional scaling). Additionally, we studied a sample from a Polish individual of well-documented genealogy whose biogeographic ancestry had previously been inferred by commercial genomic services using autosomal single nucleotide polymorphisms (SNPs), mitochondrial DNA and Y-SNP markers. For comparative purposes, we used genotype data collected in the "forInDel" browser and allele frequencies from previously published papers. The results obtained for InDels and STRs show that the Slavic populations constitute a genetically homogeneous group, with the exception of the Czechs differing clearly from the other tested populations. The analysis of the known Polish sample in the Snipper application proves the usefulness of the InDel markers on the continental level only. Conversely, microsatellites not only improve prediction, but are also informative if considered as an independent set of ancestry markers.

Complete mitochondrial genome and phylogeny of Pleistocene mammoth Mammuthus primigenius.

Phylogenetic relationships between the extinct woolly mammoth (Mammuthus primigenius), and the Asian (Elephas maximus) and African savanna (Loxodonta africana) elephants remain unresolved. Here, we report the sequence of the complete mitochondrial genome (16,842 base pairs) of a woolly mammoth extracted from permafrost-preserved remains from the Pleistocene epoch--the oldest mitochondrial genome sequence determined to date. We demonstrate that well-preserved mitochondrial genome fragments, as long as approximately 1,600-1700 base pairs, can be retrieved from pre-Holocene remains of an extinct species. Phylogenetic reconstruction of the Elephantinae clade suggests that M. primigenius and E. maximus are sister species that diverged soon after their common ancestor split from the L. africana lineage. Low nucleotide diversity found between independently determined mitochondrial genomic sequences of woolly mammoths separated geographically and in time suggests that north-eastern Siberia was occupied by a relatively homogeneous population of M. primigenius throughout the late Pleistocene.

Optimization of the Y831C mutation detection in human DNA polymerase gamma by allelic discrimination assay.

Many well-defined mutations in the gene for the catalytic subunit of polymerase gamma (POLG1) have been found to be associated with disease, whereas the status of several mutations remains unresolved due to the conflicting reports on their frequencies in populations of healthy individuals. Here, we have developed a highly sensitive, real-time allelic discrimination assay enabling detection of the Y831C mutation in the POLG1 gene. The Y831C mutation is present in the Polish population at a frequency of 2.25%. The new assay is well suited to both extensive population studies and molecular diagnostics of POLG1.

Genetic variation of 15 STR loci (D3S1358, vWA, FGA, TH01, TPOX, CSF1PO, D5S818, D13S317, D7S820, D16S539, D2S1338, D8S1179, D21S11, D18S51, and D19S433) in populations of north and central Poland.

Allele frequencies for 15 short tandem repeats (STRs) included in AmpFlSTR Profiler and AmpFlSTR SGM Plus kits were obtained from populations of Pomorze Gdanskie, Wielkopolska, Kujawy, Pomorze Zachodnie, Mazury and Mazowsze regions of Poland.

Simple and cost-effective 14-loci SNP assay designed for differentiation of European, East Asian and African samples.

During every criminal investigation, it is vital to extract as much information as possible from every piece of evidence. When it comes to DNA testing, simple short tandem repeat (STR) typing may soon become a relic because it is now possible to genotype more characteristics. Ancestry informative markers are receiving attention from the forensic community because individuals can be assigned to their population or territory of origin based on their analysis. Many panels of this kind have been proposed so far, yet most of them require typing of a large number of loci. In many cases it is crucial to pick a minimal set of the most informative markers due to the limited amount of material available for analysis. In this study, we demonstrate that 14 carefully picked SNPs combined in two multiplex assays are capable of fast, robust and cost effective three-way differentiation of East Asians, Europeans and Africans.

Mitochondrial DNA variation in the aboriginal populations of the Altai-Baikal region: implications for the genetic history of North Asia and America.

The discovery of mtDNA types common to Asians and Amerinds (types A, B, C, and D) forced investigators to search for those nations of Asia which, though not considered the ancestors of the Amerinds, have retained a close genetic resemblance with them. We collected samples and studied the gene pools of the Turkic-speaking nations of South Siberia: Altaians, Khakassians, Shorians, Tuvinians, Todjins, Tofalars, Sojots, as well as Mongolian-speaking Buryats. The data indicate that nearly all Turkic-speaking nations of Siberia and Central Asia, as well as the Buryats, have types A, B, C, and D in their gene pool. The highest total frequency of these types is observed in the Tuvinians and Sojots. They, as well as the Buryats, also have the lowest frequency of the europeoid types. The most mixed Asian-Europeoid gene pool examined turned out to be that of the Shorians. An important finding was the presence of type X in the Altaians, which had not yet been detected in Asia. As shown by computer analysis, this DNA sequence is not a late European admixture. Rather, the Altai variant X is ancient and can be close to the ancestral form of the variants of contemporary Europeans and Amerinds. The presented results prove that of all nations in Asia, the Turkic-speaking nations living between Altai and Baikal along the Sayan mountains are genetically closest to the Amerinds.

Genome-wide analysis of cold adaptation in indigenous Siberian populations.

Following the dispersal out of Africa, where hominins evolved in warm environments for millions of years, our species has colonised different climate zones of the world, including high latitudes and cold environments. The extent to which human habitation in (sub-)Arctic regions has been enabled by cultural buffering, short-term acclimatization and genetic adaptations is not clearly understood. Present day indigenous populations of Siberia show a number of phenotypic features, such as increased basal metabolic rate, low serum lipid levels and increased blood pressure that have been attributed to adaptation to the extreme cold climate. In this study we introduce a dataset of 200 individuals from ten indigenous Siberian populations that were genotyped for 730,525 SNPs across the genome to identify genes and non-coding regions that have undergone unusually rapid allele frequency and long-range haplotype homozygosity change in the recent past. At least three distinct population clusters could be identified among the Siberians, each of which showed a number of unique signals of selection. A region on chromosome 11 (chr11:66-69 Mb) contained the largest amount of clustering of significant signals and also the strongest signals in all the different selection tests performed. We present a list of candidate cold adaption genes that showed significant signals of positive selection with our strongest signals associated with genes involved in energy regulation and metabolism (CPT1A, LRP5, THADA) and vascular smooth muscle contraction (PRKG1). By employing a new method that paints phased chromosome chunks by their ancestry we distinguish local Siberian-specific long-range haplotype signals from those introduced by admixture.

Developing STR databases on structured populations: the native South Siberian population versus the Russian population.

Developing a forensic DNA database on a population that consists of local ethnic groups separated by physical and cultural barriers is questionable as it can be genetically subdivided. On the other side, small sizes of ethnic groups, especially in alpine regions where they are sub-structured further into small villages, prevent collecting a large sample from each ethnic group. For such situations, we suggest to obtain both a total population database on allele frequencies across ethnic groups and a list of theta-values between the groups and the total data. We have genotyped 558 individuals from the native population of South Siberia, consisting of nine ethnic groups, at 17 autosomal STR loci of the kit packages AmpFlSTR SGM Plus i, Cyrillic AmpFlSTR Profiler Plus. The groups differentiate from each other with average theta-values of around 1.1%, and some reach up to three to four percent at certain loci. There exists between-village differentiation as well. Therefore, a database for the population of South Siberia is composed of data on allele frequencies in the pool of ethnic groups and data on theta-values that indicate variation in allele frequencies across the groups. Comparison to additional data on northeastern Asia (the Chukchi and Koryak) shows that differentiation in allele frequencies among small groups that are separated by large geographic distance can be even greater. In contrast, populations of Russians that live in large cities of the European part of Russia are homogeneous in allele frequencies, despite large geographic distance between them, and thus can be described by a database on allele frequencies alone, without any specific information on theta-values.

Complex interactions of the Eastern and Western Slavic populations with other European groups as revealed by mitochondrial DNA analysis.

Mitochondrial DNA sequence variation was examined by the control region sequencing (HVS I and HVS II) and RFLP analysis of haplogroup-diagnostic coding region sites in 570 individuals from four regional populations of Poles and two Russian groups from northwestern part of the country. Additionally, sequences of complete mitochondrial genomes representing K1a1b1a subclade in Polish and Polish Roma populations have been determined. Haplogroup frequency patterns revealed in Poles and Russians are similar to those characteristic of other Europeans. However, there are several features of Slavic mtDNA pools seen on the level of regional populations which are helpful in the understanding of complex interactions of the Eastern and Western Slavic populations with other European groups. One of the most important is the presence of subhaplogroups U5b1b1, D5, Z1 and U8a with simultaneous scarcity of haplogroup K in populations of northwestern Russia suggesting the participation of Finno-Ugrian tribes in the formation of mtDNA pools of Russians from this region. The results of genetic structure analyses suggest that Russians from Velikii Novgorod area (northwestern Russia) and Poles from Suwalszczyzna (northeastern Poland) differ from all remaining Polish and Russian samples. Simultaneously, northwestern Russians and northeastern Poles bear some similarities to Baltic (Latvians) and Finno-Ugrian groups (Estonians) of northeastern Europe, especially on the level of U5 haplogroup frequencies. The occurrence of K1a1b1a subcluster in Poles and Polish Roma is one of the first direct proofs of the presence of Ashkenazi-specific mtDNA lineages in non-Jewish European populations.

Mitochondrial DNA variability in the Czech population, with application to the ethnic history of Slavs.

Mitochondrial DNA (mtDNA) variability was studied in a sample of 179 individuals representing the Czech population of Western Bohemia. Sequencing of two hypervariable segments, HVS I and HVS II, in combination with screening of coding-region haplogroup-specific RFLP markers revealed that most Czech mtDNAs belong to the common West Eurasian mitochondrial haplogroups (H, pre-V HV*, J, T, U, N1, W, and X). However, about 3% of Czech mtDNAs encompass East Eurasian lineages (A, N9a, D4, M*). A comparative analysis with published data showed that different Slavonic populations in Central and Eastern Europe contain small but marked amounts of East Eurasian mtDNAs. We suggest that the presence of East Eurasian mtDNA haplotypes is not an original feature of the gene pool of the proto-Slavs but rather may be mostly a consequence of admixture with Central Asian nomadic tribes, who migrated into Central and Eastern Europe in the early Middle Ages.

Low variability of the POLG (CAG)n repeat in north Eurasian populations.

We investigated the frequency of different repeat-length alleles of the trinucleotide CAG microsatellite repeat in the coding sequence of the nuclear gene for the catalytic subunit of mitochondrial DNA polymerase gamma (POLG) in 12 ethnic groups from northern Eurasia. The population sample consisted of 1,330 individuals from 3 large geographic areas: Europe, Southwest Asia, and Siberia/East Asia. We found that the 10-repeat allele of the POLG gene is the most frequent in all analyzed populations, with a frequency of 88-96%. The heterozygosity level ranges from 22% in Europe to 13.6% in Southwest Asia with the lowest value of 7.4% in Siberia/East Asia. The present study provides evidence of clinal distribution of POLG gene heterozygosity in North Eurasian populations. In general, we found an extremely low variability of the trinucleotide CAG microsatellite repeat, suggesting that purifying selection acts against deleterious alleles, although low mutability of the repeated region cannot be ruled out.

Analysis of phylogenetically reconstructed mutational spectra in human mitochondrial DNA control region.

Analysis of mutations in mitochondrial DNA is an important issue in population and evolutionary genetics. To study spontaneous base substitutions in human mitochondrial DNA we reconstructed the mutational spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. An excess of pyrimidine transitions was found both in HVS I and II regions. By means of classification analysis numerous mutational hotspots were revealed in these spectra. Context analysis of hotspots revealed a complex influence of neighboring bases on mutagenesis in the HVS I region. Further statistical analysis suggested that a transient misalignment dislocation mutagenesis operating in monotonous runs of nucleotides play an important role for generating base substitutions in mitochondrial DNA and define context properties of mtDNA. Our results suggest that dislocation mutagenesis in HVS I and II is a fingerprint of errors produced by DNA polymerase gamma in the course of human mitochondrial DNA replication

Disuniting uniformity: a pied cladistic canvas of mtDNA haplogroup H in Eurasia.

It has been often stated that the overall pattern of human maternal lineages in Europe is largely uniform. Yet this uniformity may also result from an insufficient depth and width of the phylogenetic analysis, in particular of the predominant western Eurasian haplogroup (Hg) H that comprises nearly a half of the European mitochondrial DNA (mtDNA) pool. Making use of the coding sequence information from 267 mtDNA Hg H sequences, we have analyzed 830 mtDNA genomes, from 11 European, Near and Middle Eastern, Central Asian, and Altaian populations. In addition to the seven previously specified subhaplogroups, we define fifteen novel subclades of Hg H present in the extant human populations of western Eurasia. The refinement of the phylogenetic resolution has allowed us to resolve a large number of homoplasies in phylogenetic trees of Hg H based on the first hypervariable segment (HVS-I) of mtDNA. As many as 50 out of 125 polymorphic positions in HVS-I were found to be mutated in more than one subcluster of Hg H. The phylogeographic analysis revealed that sub-Hgs H1*, H1b, H1f, H2a, H3, H6a, H6b, and H8 demonstrate distinct phylogeographic patterns. The monophyletic subhaplogroups of Hg H provide means for further progress in the understanding of the (pre)historic movements of women in Eurasia and for the understanding of the present-day genetic diversity of western Eurasians in general.