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.

Searching for improvements in predicting human eye colour from DNA.

Increasing understanding of human genome variability allows for better use of the predictive potential of DNA. An obvious direct application is the prediction of the physical phenotypes. Significant success has been achieved, especially in predicting pigmentation characteristics, but the inference of some phenotypes is still challenging. In search of further improvements in predicting human eye colour, we conducted whole-exome (enriched in regulome) sequencing of 150 Polish samples to discover new markers. For this, we adopted quantitative characterization of eye colour phenotypes using high-resolution photographic images of the iris in combination with DIAT software analysis. An independent set of 849 samples was used for subsequent predictive modelling. Newly identified candidates and 114 additional literature-based selected SNPs, previously associated with pigmentation, and advanced machine learning algorithms were used. Whole-exome sequencing analysis found 27 previously unreported candidate SNP markers for eye colour. The highest overall prediction accuracies were achieved with LASSO-regularized and BIC-based selected regression models. A new candidate variant, rs2253104, located in the ARFIP2 gene and identified with the HyperLasso method, revealed predictive potential and was included in the best-performing regression models. Advanced machine learning approaches showed a significant increase in sensitivity of intermediate eye colour prediction (up to 39%) compared to 0% obtained for the original IrisPlex model. We identified a new potential predictor of eye colour and evaluated several widely used advanced machine learning algorithms in predictive analysis of this trait. Our results provide useful hints for developing future predictive models for eye colour in forensic and anthropological studies.

Variability of the rs333 in Polish patients with lupus erythematosus.

INTRODUCTION: Lupus erythematosus (LE) is an autoimmune disease with a strong influence of genetic and environmental factors. C-C motif chemokine receptor 5 (CCR5) gene expression may affect the development and intensity of LE. AIM: To evaluate the possible association between the 32bp deletion in rs333 locus located within the CCR5 gene and the development of LE or the occurrence of various clinical symptoms in the course of the disease. MATERIAL AND METHODS: One hundred and twenty patients with LE (77 with systemic lupus erythematosus (SLE) and 43 with discoid lupus erythematosus (DLE)) and 100 healthy controls from the Polish population were genotyped for deletion in rs333. RESULTS: 32 bp deletion in the rs333 was significantly more frequent among healthy individuals than DLE patients. Moreover, heterozygotes and homozygotes with deletion in rs333 were significantly more frequent within the control group than the group of patients with discoid lupus erythematosus. In contrast, any statistically significant differences in allele or genotype frequencies between healthy persons and SLE patients were observed. Furthermore, nucleotide sequence variability of rs333 was not associated with certain clinical symptoms of LE patients. CONCLUSIONS: Deletion in the rs333 might be a protective factor for DLE, but not SLE in the Polish population. Nevertheless further studies performed on larger populations are needed to confirm these observations.

Exploring the possibility of predicting human head hair greying from DNA using whole-exome and targeted NGS data.

BACKGROUND: Greying of the hair is an obvious sign of human aging. In addition to age, sex- and ancestry-specific patterns of hair greying are also observed and the progression of greying may be affected by environmental factors. However, little is known about the genetic control of this process. This study aimed to assess the potential of genetic data to predict hair greying in a population of nearly 1000 individuals from Poland. RESULTS: The study involved whole-exome sequencing followed by targeted analysis of 378 exome-wide and literature-based selected SNPs. For the selection of predictors, the minimum redundancy maximum relevance (mRMRe) method was used, and then two prediction models were developed. The models included age, sex and 13 unique SNPs. Two SNPs of the highest mRMRe score included whole-exome identified KIF1A rs59733750 and previously linked with hair loss FGF5 rs7680591. The model for greying vs. no greying prediction achieved accuracy of cross-validated AUC = 0.873. In the 3-grade classification cross-validated AUC equalled 0.864 for no greying, 0.791 for mild greying and 0.875 for severe greying. Although these values present fairly accurate prediction, most of the prediction information was brought by age alone. Genetic variants explained < 10% of hair greying variation and the impact of particular SNPs on prediction accuracy was found to be small. CONCLUSIONS: The rate of changes in human progressive traits shows inter-individual variation, therefore they are perceived as biomarkers of the biological age of the organism. The knowledge on the mechanisms underlying phenotypic aging can be of special interest to the medicine, cosmetics industry and forensics. Our study improves the knowledge on the genetics underlying hair greying processes, presents prototype models for prediction and proves hair greying being genetically a very complex trait. Finally, we propose a four-step approach based on genetic and epigenetic data analysis allowing for i) sex determination; ii) genetic ancestry inference; iii) greying-associated SNPs assignment and iv) epigenetic age estimation, all needed for a final prediction of greying.

Complete mitogenome data for the Serbian population: the contribution to high-quality forensic databases.

Mitochondrial genome (mtDNA) is a valuable resource in resolving various human forensic casework. The usage of variability of complete mtDNA genomes increases their discriminatory power to the maximum and enables ultimate resolution of distinct maternal lineages. However, their wider employment in forensic casework is nowadays limited by the lack of appropriate reference database. In order to fill in the gap in the reference data, which, considering Slavic-speaking populations, currently comprises only mitogenomes of East and West Slavs, we present mitogenome data for 226 Serbians, representatives of South Slavs from the Balkan Peninsula. We found 143 (sub)haplogroups among which West Eurasian ones were dominant. The percentage of unique haplotypes was 85%, and the random match probability was as low as 0.53%. We support previous findings on both high levels of genetic diversity in the Serbian population and patterns of genetic differentiation among this and ten studied European populations. However, our high-resolution data supported more pronounced genetic differentiation among Serbians and two Slavic populations (Russians and Poles) as well as expansion of the Serbian population after the Last Glacial Maximum and during the Migration period (fourth to ninth century A.D.), as inferred from the Bayesian skyline analysis. Phylogenetic analysis of haplotypes found in Serbians contributed towards the improvement of the worldwide mtDNA phylogeny, which is essential for the interpretation of the mtDNA casework.

The microbiome and its impact on food allergy and atopic dermatitis in children.

Food allergy (FA) affects 4-10% of children, especially children with atopic dermatitis (AD). During infancy the gut microbiome may determine both the course of FA and tolerance to food allergens. Analogically, the skin microbiome changes in the course of AD. Most studies have associated FA with a lower abundance and diversity of Lactobacillales and Clostridiales, but greater numbers of Enterobacterales, while AD in children has been associated with lower numbers of Staphylococcus epidermidis and S. hominis but an abundance of S. aureus and Streptococcus species. An understanding of the impact of the microbiome on the clinical course of FA and AD may allow for the development of new models of allergy treatment and prevention.

TP53 somatic mutations are associated with somatic mitogenome substitutions but not indels in colorectal cancer cells.

BACKGROUND: p53 is a tumour suppressor protein that is involved in many cancer-related processes. Growing evidence suggests that p53 also plays an important role in mitochondrial (mtDNA) maintenance. Somatic mitogenome mutations are frequently observed in colorectal cancer (CC) cells. Thus, it was important to determine whether somatic mtDNA changes are associated with TP53 mutational status. METHODS: In the present study, we analysed the TP53 gene in 67 CC patients, for whom mitogenome haplotypes were previously described. In total, 134 TP53 sequences (of cancer and matched normal specimens) were determined using the dideoxy method. RESULTS: Nine hereditary polymorphisms in the TP53 gene were detected in normal colon cells. None of them (neither alleles, nor genotypes) was associated with somatic mitogenome mutations in CC cells. Moreover, 42 somatic TP53 mutations were found in approximately 36% of CC tissues. These somatic changes were significantly more frequent in CC cells with somatic mtDNA mutations (p = 0.0069). Furthermore, we show that only mitochondrial somatic substitutions (p = 0.0017), but not indels (p > 0.05), were associated with somatic TP53 mutations. CONCLUSIONS: The results of the present study suggest that changes in TP53 may modify p53 properties, which may result in the accumulation of somatic substitutions in CC mitogenomes.

Mitogenomic differences between the normal and tumor cells of colorectal cancer patients.

So far, a reliable spectrum of mitochondrial DNA mutations in colorectal cancer cells is still unknown, and neither is their significance in carcinogenesis. Indeed, it remains debatable whether mtDNA mutations are "drivers" or "passengers" of colorectal carcinogenesis. Thus, we analyzed 200 mitogenomes from normal and cancer tissues of 100 colorectal cancer patients. Minority variant mutations were detected at the 1% level. We showed that somatic mutations frequently occur in colorectal cancer cells (75%) and are randomly distributed across the mitochondrial genome. Mutational signatures of somatic mitogenome mutations suggest that they might arise through nucleotide deamination due to oxidative stress. The majority of somatic mutations localized within the coding region (in positions not known from the human phylogeny) and was potentially pathogenic to cell metabolism. Further analysis suggested that the relaxation of negative selection in the mitogenomes of colorectal cancer cells may allow accumulation of somatic mutations. Thus, a shift in glucose metabolism from oxidative phosphorylation to glycolysis may create advantageous conditions for accumulation of mtDNA mutations. Considering the fact that the presence of somatic mtDNA mutations was not associated with any clinicopathological features, we suggested that mtDNA somatic mutations are "passengers" rather than the cause of colorectal carcinogenesis.

Whole mitochondrial genome diversity in two Hungarian populations.

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

An Analysis of the Association between Epilepsy-Related Genes and Vertigo in the Polish Population.

Considering the possibility of a common genetic background of vertigo and epilepsy, we genotyped an affected group of individuals with vertigo and an unaffected group, by studying 26 single-nucleotide polymorphisms (SNPs) in 14 genes which were previously reported to be of particular importance for epilepsy. Significant differences were found between the patients and the control group (χ2 = 38.3, df = 3, p = 1.6 × 10-7) for the frequencies of haplotypes consist ing of 2 SNPs located in chromosome 11 (rs1939012 and rs1783901 within genes MMP8 and SCN3B, respectively). The haplotype rs1939012:C-rs1783901:A, consisting of the minor-frequency alleles was found to be associated with a higher risk of vertigo (OR = 5.0143, 95% CI = 1.6991-14.7980, p = 0.0035). In contrast, the haplotype rs1939012:T-rs1783901:A showed a significant association with a decreased risk of the disease (OR = 0.0597, 95% CI = 0.0136-0.2620, p = 0.0002). Our results suggest that the SNPs rs1939012 and rs1783901 may play a potential role of gene regulation and/or epistasis in a complex etiology of vertigo.

The Association of Apolipoprotein E Gene Polymorphism With Cognitive Performance in Nondemented Polish Adults Aged 55 to 75.

The ε4 allele of the apolipoprotein E (APOE) gene is known as a risk factor for dementia. How APOE ε polymorphism affects cognitive performance in nondemented aging subjects remains less clear. In this study, the relationship between APOE status and cognitive performance across various cognitive domains in adults aged 55 to 75 years ( n = 74) without dementia was investigated. E4 carriers ( n = 11) performed worse versus noncarriers on forward Digit Span and delayed recall of the Rey-Osterrieth complex figure. General linear model analysis revealed a small but significant main effect of ε4 on Rey-Osterrieth complex figure delayed recall. Comparing ε2 carriers, ε3 homozygotes, and ε4 carriers, ε3/ε3 performed significantly better on Trail Making Test part B and derived score Trail Making Test B-A. The findings support the relation between the APOE ε polymorphism and visual memory, short-term auditory memory, visuospatial attention, and executive functions in an aging sample without dementia.

X-linked TLR7 gene polymorphisms are associated with diverse immunological conditions but not with discoid lupus erythematosus in Polish patients.

INTRODUCTION: Toll-like receptor 7 (TLR7) is an important molecule involved in the development of autoimmunity and the response to different pathogens. Several polymorphisms within the TLR7 gene were previously found to be associated with systemic lupus erythematosus (SLE). However, none of those studies investigated the TLR7 promoter flanking variants rs1634318 and rs1616583. TLR7 gene diversity has not been analyzed with respect to discoid lupus erythematosus (DLE) development, while its role in the human immunological response to fungal infection is not fully known. AIM: To clarify the potential involvement of two novel single-nucleotide polymorphisms (SNPs) located in the TLR7 gene (rs1634318 and rs1616583) in a variety of immune-related conditions, we studied the variability of these loci in patients from a Polish population with SLE and DLE, as well as in immunocompromised patients who were affected by invasive aspergillosis (IA) and those who were not affected. MATERIAL AND METHODS: Real-time polymerase chain reaction was used to genotype SNPs. Statistically significant differences between case and control groups for both allele and genotype frequencies were assessed using the χ2 test with Yates' correction or two-tailed Fisher's exact test. The results were Bonferroni-corrected for multiple comparisons and odds ratios were calculated. RESULTS: Two polymorphisms located in TLR7 might be associated with the development of SLE but not DLE within the Polish population. Moreover, variation of the two investigated SNPs was found to be associated with IA in immunocompromised Polish patients. CONCLUSIONS: In Polish patients, TLR7 promoter flanking gene polymorphisms might be associated with IA and SLE but not DLE.

Mitochondrial super-haplogroup U diversity in Serbians.

BACKGROUND: Available mitochondrial (mtDNA) data demonstrate genetic differentiation among South Slavs inhabiting the Balkan Peninsula. However, their resolution is insufficient to elucidate the female-specific aspects of the genetic history of South Slavs, including the genetic impact of various migrations which were rather common within the Balkans, a region having a turbulent demographic history. AIM: The aim was to thoroughly study complete mitogenomes of Serbians, a population linking westward and eastward South Slavs. SUBJECTS AND METHODS: Forty-six predominantly Serbian super-haplogroup U complete mitogenomes were analysed phylogenetically against ∼4000 available complete mtDNAs of modern and ancient Western Eurasians. RESULTS: Serbians share a number of U mtDNA lineages with Southern, Eastern-Central and North-Western Europeans. Putative Balkan-specific lineages (e.g. U1a1c2, U4c1b1, U5b3j, K1a4l and K1a13a1) and lineages shared among Serbians (South Slavs) and West and East Slavs were detected (e.g. U2e1b1, U2e2a1d, U4a2a, U4a2c, U4a2g1, U4d2b and U5b1a1). CONCLUSION: The exceptional diversity of maternal lineages found in Serbians may be associated with the genetic impact of both autochthonous pre-Slavic Balkan populations whose mtDNA gene pool was affected by migrations of various populations over time (e.g. Bronze Age pastoralists) and Slavic and Germanic newcomers in the early Middle Ages.

Massively parallel targeted resequencing reveals novel genetic variants associated with aspergillosis in paediatric patients with haematological malignancies.

This study aimed to find novel genetic variants of susceptibility to aspaergillosis in paediatric patients with haematological malignancies. Complete sequences of fifteen genes of human innate immunity (CCL2, CCR2, CD209, CLEC6A, CLEC7A and ten TLR genes) were studied in 40 patients diagnosed with haematological disorders (20 unaffected and 20 affected by aspergillosis). All samples were sequenced with MiSeq (Illumina) and 454 (Roche Diagnostics) technologies. Statistical significance of the differences between studied groups was determined using the two-tailed Fisher's exact test. Sixty variants of potential importance were identified, the vast majority of which are located in non-coding parts of the targeted genes. At the threshold of p < 0.000005, one intergenic (TLR2 rs4585282) and one intronic variant (CLEC6A rs12099687) were found significant between the case and control groups for genotype and allele frequencies, respectively. Rs12099687 in CLEC6A was predicted to constitute an alternative isoform or cryptic splice site, which potentially changes activity of the Dectin-2 protein. Overall, we assume that the two strongest associations reported in this study are expected to be reproducible even in the absence of other evidence, while another twelve associations may be strong enough to justify additional research in larger cohorts.

Genetic similarities and differences between discoid and systemic lupus erythematosus patients within the Polish population.

INTRODUCTION: Many studies have shown that some SNPs might be a risk factor for systemic lupus erythematosus (SLE), but little is known about potential susceptibility loci of the skin types of the disease. Discoid lupus erythematosus (DLE) is the most common form of the cutaneous lupus erythematosus. Nevertheless, a genetic contribution to DLE is not fully recognized. AIM: We aimed to analyze three SNPs located in the STAT4 (rs7574865), ITGAM (rs1143679) and TNXB (rs1150754) genes in both DLE and SLE patients from Poland. MATERIAL AND METHODS: SNPs were genotyped using real-time polymerase chain reaction (PCR). Statistical significance of the differences between patient and control groups in both allele and genotype frequencies were calculated using two tailed Fisher's exact test. The correction for multiple testing by the Bonferroni adjustment and odds ratio were also calculated. RESULTS: For the first time, we have shown that the polymorphisms located in the STAT4 (rs7574865), but not in the ITGAM (rs1143679) nor the TNXB (rs1150754) genes, might be associated with the development of DLE within the Polish population. The variation of the three investigated SNPs was found to be associated with SLE in our dataset. CONCLUSIONS: The results of our study suggest differences in the molecular background between DLE and SLE within the Polish population.

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.

RNA isolation from bloodstains collected on FTA cards - application in clinical and forensic genetics.

Aim of the study: In recent years, RNA analysis has been increasingly used in clinical and forensic genetics. Nevertheless, a major limitation of RNA-based applications is very low RNA stability in biological material, due to the RNAse activity. This highlights the need for improving the methods of RNA collection and storage. Technological approaches such as FTA Classic Cards (Whatman) could provide a solution for the problem of RNA degradation. However, different methods of RNA isolation from FTA cards could have diverse effects on RNA quantity and quality. The purpose of this research was to analyze the utility of three different methods of RNA isolation from peripheral blood collected on FTA Classic Cards (Whatman). The study also aimed at assessing RNA stability in bloodstains deposited on FTA cards. Material and methods: The study was performed on peripheral bloodstains collected from 59 individuals on FTA Classic Cards (Whatman). RNA was isolated with High Pure RNA Isolation Kit (Roche Diagnostics), Universal RNA/miRNA Purification (EURx) and TRIzol Reagent (Life Technologies). RNA was subjected to quantitative analysis followed by reverse transcription and Real - Time PCR reaction. Results: The study has shown that FTA Classic Cards (Whatman) are useful tools for storing bloodstains at room temperature for RNA analysis. Moreover, the method of RNA extraction employing TRIzol Reagent (Life Technologies) provides the highest efficiency and reproducibility for samples stored for no more than 2 years. Conclusions: The FTA cards are suitable for collecting and storing bloodstains for RNA analysis in clinical and forensic genetics.

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.

On the origin of mongrels: evolutionary history of free-breeding dogs in Eurasia.

Although a large part of the global domestic dog population is free-ranging and free-breeding, knowledge of genetic diversity in these free-breeding dogs (FBDs) and their ancestry relations to pure-breed dogs is limited, and the indigenous status of FBDs in Asia is still uncertain. We analyse genome-wide SNP variability of FBDs across Eurasia, and show that they display weak genetic structure and are genetically distinct from pure-breed dogs rather than constituting an admixture of breeds. Our results suggest that modern European breeds originated locally from European FBDs. East Asian and Arctic breeds show closest affinity to East Asian FBDs, and they both represent the earliest branching lineages in the phylogeny of extant Eurasian dogs. Our biogeographic reconstruction of ancestral distributions indicates a gradual westward expansion of East Asian indigenous dogs to the Middle East and Europe through Central and West Asia, providing evidence for a major expansion that shaped the patterns of genetic differentiation in modern dogs. This expansion was probably secondary and could have led to the replacement of earlier resident populations in Western Eurasia. This could explain why earlier studies based on modern DNA suggest East Asia as the region of dog origin, while ancient DNA and archaeological data point to Western Eurasia.

Western Eurasian ancestry in modern Siberians based on mitogenomic data.

BACKGROUND: Although the genetic heritage of aboriginal Siberians is mostly of eastern Asian ancestry, a substantial western Eurasian component is observed in the majority of northern Asian populations. Traces of at least two migrations into southern Siberia, one from eastern Europe and the other from western Asia/the Caucasus have been detected previously in mitochondrial gene pools of modern Siberians. RESULTS: We report here 166 new complete mitochondrial DNA (mtDNA) sequences that allow us to expand and re-analyze the available data sets of western Eurasian lineages found in northern Asian populations, define the phylogenetic status of Siberian-specific subclades and search for links between mtDNA haplotypes/subclades and events of human migrations. From a survey of 158 western Eurasian mtDNA genomes found in Siberia we estimate that nearly 40% of them most likely have western Asian and another 29% European ancestry. It is striking that 65 of northern Asian mitogenomes, i.e. ~41%, fall into 19 branches and subclades which can be considered as Siberian-specific being found so far only in Siberian populations. From the coalescence analysis it is evident that the sequence divergence of Siberian-specific subclades was relatively small, corresponding to only 0.6-9.5 kya (using the complete mtDNA rate) and 1-6 kya (coding region rate). CONCLUSIONS: The phylogeographic analysis implies that the western Eurasian founders, giving rise to Siberian specific subclades, may trace their ancestry only to the early and mid-Holocene, though some of genetic lineages may trace their ancestry back to the end of Last Glacial Maximum (LGM). We have not found the modern northern Asians to have western Eurasian genetic components of sufficient antiquity to indicate traces of pre-LGM expansions.