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.

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.

Y chromosome haplotype diversity in Mongolic-speaking populations and gene conversion at the duplicated STR DYS385a,b in haplogroup C3-M407.

Y chromosome microsatellite (Y-STR) diversity has been studied in different Mongolic-speaking populations from South Siberia, Mongolia, North-East China and East Europe. The results obtained indicate that the Mongolic-speaking populations clustered into two groups, with one group including populations from eastern part of South Siberia and Central Asia (the Buryats, Barghuts and Khamnigans) and the other group including populations from western part of Central Asia and East Europe (the Mongols and Kalmyks). High frequency of haplogroup C3-M407 (>50%) is present in the Buryats, Barghuts and Khamnigans, whereas in the Mongols and Kalmyks its frequency is much lower. In addition, two allelic combinations in DYS385a,b loci of C3-M407 haplotypes have been observed: the combination 11,18 (as well as 11,17 and 11,19) is frequent in different Mongolic-speaking populations, but the 11,11 branch is present mainly in the Kalmyks and Mongols. Results of locus-specific sequencing suggest that the action of gene conversion is a more likely explanation for origin of homoallelic 11,11 combination. Moreover, analysis of median networks of Y-STR haplotypes demonstrates that at least two gene conversion events can be revealed-one of them has probably occurred among the Mongols, and the other event occurred in the Barghuts. These two events give an average gene conversion rate range of 0.24-7.1 × 10(-3) per generation.

Mitochondrial DNA Polymerase γ Mutations and Their Implications in mtDNA Alterations in Colorectal Cancer.

Mitochondrial DNA was found to be highly mutated in colorectal cancer cells. One of the key molecules involved in the maintenance of the mitochondrial genome is the nuclear-encoded polymerase gamma. The aim of our study was to determine if there is a link between polymorphisms within the polymerase gamma gene (POLG) and somatic mutations within the mitochondrial genome in cancer cells. We investigated POLG sequence variability in 50 colorectal cancer patients whose complete mitochondrial genome sequences were determined. Relative mtDNA copy number was also determined. We identified 251 sequence variants in the POLG gene. Most of them were germline-specific (∼92%). Twenty-one somatic changes in POLG were found in 10 colorectal cancer patients. We have found no association between the occurrence of mtDNA somatic mutations and the somatically occurring variants in POLG. MtDNA content was reduced in patients carrying somatic variants in POLG or germline nucleotide variants located in the region encoding the POLG polymerase domain, but the difference did not reach statistical significance. Our findings suggest that somatic mtDNA mutations occurring in colorectal cancer are not a consequence of somatic mutations in POLG. Nevertheless, POLG nucleotide variants may lead to a decrease in mtDNA content, and consequently result in mitochondrial dysfunction.

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.

A novel multiplex assay amplifying 13 Y-STRs characterized by rapid and moderate mutation rate.

As microsatellites located on Y chromosome mutate with different rates, they may be exploited in evolutionary studies, genealogical testing of a variety of populations and even, as proven recently, aid individual identification. Currently available commercial Y-STR kits encompass mostly low to moderately mutating loci, making them a perfect choice for the first two applications. Some attempts have been made so far to utilize Y-STRs to provide a discriminatory tool for forensic purposes. Although all 13 rapidly mutating Y-STRs were already multiplexed, no single assay based on single-copy markers allowing at least a portion of close male relatives to be differentiated from one another is available. To fill in the blanks, we constructed and validated an assay comprised of single-copy Y-STR markers only with a mutation rate ranging from 8×10(-3) to 1×10(-2). Performance of the resulting combination of nine RM Y-STRs and four moderately mutating ones was tested on 361 father-son pairs and 1326 males from 9 populations revealing an overall mutation rate of 1.607×10(-1) for the assay as a whole. Application of the proposed 13 Y-STR set to differentiation of haplotypes present among homogenous population of Buryats resulted in a threefold increase of discrimination as compared with 10 Y-STRs from the PowerPlex(®) Y.

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.

Genetic data from Y chromosome STR and SNP loci in Ukrainian population.

We have tested a sample of 154 unrelated males from Lviv region (Ukraine) for 11 Y-chromosomal single nucleotide polymorphisms (SNPs) and 17 Y-chromosomal STR loci (DYS19, DYS385a, DYS385b, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, YGATA_H4.1). Haplotype and haplogroup diversity values were calculated for the population under study. Genetic distances (R(ST)) to 9 other Slavic populations were calculated based on 12 Y-STR loci. Haplotype frequencies and MDS plots were constructed based on genetic distances. Haplogroup frequency patterns revealed in Ukraine are similar to those characteristic of other European populations. However, it also allowed for identification a specific genetic component in Ukrainian sample which seems to originate from areas dwelled by Western Slavs, i.e. subhaplogroup R1a1a7, at frequency of 13.65%. Analysis of R(ST) distances and AMOVA revealed high level of heterogeneity between Slavic populations inhabiting the south and north part of Europe, determined geographically rather than by linguistic factors. It has also been found a closer similarity (in the values of R(ST)) between Ukrainian and Slovak populations than between Ukrainians and other Slavic population samples.

On the Y-chromosome haplogroup C3c classification.

As there are ambiguities in classification of the Y-chromosome haplogroup C3c, relatively frequent in populations of Northern Asia, we analyzed all three haplogroup-defining markers M48, M77 and M86 in C3-M217-individuals from Siberia, Eastern Asia and Eastern Europe. We have found that haplogroup C3c is characterized by the derived state at M48, whereas mutations at both M77 and M86 define subhaplogroup C3c1. The branch defined by M48 alone would belong to subhaplogroup C3c*, characteristic for some populations of Central and Eastern Siberia, such as Koryaks, Evens, Evenks and Yukaghirs. Subhaplogroup C3c* individuals could be considered as remnants of the Neolithic population of Siberia, based on the age of C3c*-short tandem repeat variation amounting to 4.5 ± 2.4 thousand years.

The landscape of mitochondrial DNA variation in human colorectal cancer on the background of phylogenetic knowledge.

Recently, an increasing number of studies indicate that mutations in mitochondrial genome may contribute to cancer development or metastasis. Hence, it is important to determine whether the mitochondrial DNA might be a good, clinically applicable marker of cancer. This review describes hereditary as well as somatic mutations reported in mitochondrial DNA of colorectal cancer cells. We showed here that the entire mitochondrial genome mutational spectra are different in colorectal cancer and non-tumor cells. We also placed the described mutations on the phylogenetic context, which highlighted the recurrent problem of data quality. Therefore, the most important rules for adequately assessing the quality of mitochondrial DNA sequence analysis in cancer have been summarized. As follows from this review, neither the reliable spectrum of mtDNA somatic mutations nor the association between hereditary mutations and colorectal cancer risk have been resolved. This indicates that only high resolution studies on mtDNA variability, followed by a proper data interpretation employing phylogenetic knowledge may finally verify the utility of mtDNA sequence (if any) in clinical practice.

Diversity of 15 human X chromosome microsatellite loci in Polish population.

X-STR analysis is a powerful tool in both phylogeny reconstruction and forensic investigation. Hereby, we provide new population data concerning 15 X-STR loci (included in commercially available typing kit Mentype Argus X-8 (Biotype AG, Dresden, Germany) (DXS10135, DXS8378, DXS7132, DXS10074, HPRTB, DXS10101, DXS10134 and DXS7423) and another seven (DXS6807, DXS9898, DXS101, DXS7424, DXS7133, DXS8377 and DXS10011) that were previously described by Poetsch et al. [1] obtained from a sample of 311 individuals from Poland and compared to the results previously obtained from other populations of European, Asian and African origin [2-4]. Numerous experiments seem to prove that X-STRs are valuable markers for human identification, kinship testing and even phylogenetic research - thus serving as a complement for autosomal microsatellites, Y-STRs and mtDNA [5-7].

High frequency of somatic mutations in rat liver mitochondrial DNA.

Somatic mutations in mitochondrial DNA (mtDNA) are thought to play an important role in both aging and neurodegenerative diseases although their specific contributions remain a subject of intense debate. We analyzed somatic mutations in the mtDNA control regions in the liver of Wistar rats. The mutation rate was found to be high and increased with age from 5.3x10(-4) mutations per position to 4.48x10(-3) mutations per position at 3 and 12 months of age, respectively. The vast majority of nucleotide substitutions are transitions ( approximately 95%) with A:T>G:C transitions being the most frequent type of substitution (>50%). In 3-month-old Wistar rats, approximately 40% of somatic mutations in the control region of mtDNA are significantly consistent with the model of dislocation mutagenesis which is a signature of error-prone DNA synthesis by mtDNA polymerase gamma. The results are consistent with the previous hypothesis that normal intramitochondrial dNTP pool asymmetries, which have been shown to reduce the fidelity of mtDNA polymerase gamma, substantially contribute to somatic mutagenesis of the rat mtDNA.

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.

Analysis of forensically used autosomal short tandem repeat markers in Polish and neighboring populations.

The purpose of this study was to evaluate the homogeneity of Polish populations with respect to STRs chosen as core markers of the Polish Forensic National DNA Intelligence Database, and to provide reference allele frequencies and to explore the genetic interrelationship between Poland and neighboring countries. The allele frequency distribution of 10 STRs included in the SGMplus kit was analyzed among 2176 unrelated individuals from 6 regional Polish populations and among 4321 individuals from Germany (three samples), Austria, The Netherlands, Sweden, Czech Republic, Slovakia, Belarus, Ukraine and the Russian Federation (six samples). The statistical approach consisted of AMOVA, calculation of pairwise Rst values and analysis by multidimensional scaling. We found homogeneity of present day Poland and consistent differences between Polish and German populations which contrasted with relative similarities between Russian and German populations. These discrepancies between genetic and geographic distances were confirmed by analysis of an independent data set on Y chromosome STRs. Migrations of Goths, Viking influences, German settlements in the region of Volga river and/or forced population resettlements and other events related to World War II are the historic events which might have caused these finding.

Reconstructing the phylogeny of African mitochondrial DNA lineages in Slavs.

To elucidate the origin of African-specific mtDNA lineages, revealed previously in Slavonic populations (at frequency of about 0.4%), we completely sequenced eight African genomes belonging to haplogroups L1b, L2a, L3b, L3d and M1 gathered from Russians, Czechs, Slovaks and Poles. Results of phylogeographic analysis suggest that at least part of the African mtDNA lineages found in Slavs (such as L1b, L3b1, L3d) appears to be of West African origin, testifying to an opportunity of their occurrence as a result of migrations to Eastern Europe through Iberia. However, a prehistoric introgression of African mtDNA lineages into Eastern Europe (approximately 10 000 years ago) seems to be probable only for European-specific subclade L2a1a, defined by coding region mutations at positions 6722 and 12903 and detected in Czechs and Slovaks. Further studies of the nature of African admixture in gene pools of Europeans require the essential enlargement of databases of African complete mitochondrial genomes.

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.

[Phylogeographic approach in the interpretation of mitochondrial DNA sequencing results in forensics]. / Zastosowanie analizy filogeograficznej w interpretacji wyników sekwencjonowania mitochondrialnego DNA dla celów sadowych.

In recent years, forensic mitochondrial DNA analysis has been undertaken from an evolutionary perspective. In particular, the phylogeographic approach based on a phylogenetic analysis of the spatial distribution of mitochondrial haplotypes and haplogroups appears to be a useful tool in the interpretation of identification cases. In this study, the phylogeographic approach has been employed in the analysis of three difficult forensic cases, where single nucleotide, homoplasmic differences were found between the reference and evidentiary haplotypes. mtDNA sequence variation has been examined by the control region (HVS I and HVS II) direct sequencing. Additionally, in order to clarify the subhaplogroup status of the selected haplotypes, DNA sequences of entire mitochondrial genomes obtained from two samples representing J1b subclade have been analyzed.