Response to Wyckelsma et al.: Loss of α-actinin-3 during human evolution provides superior cold resilience and muscle heat generation.

The common loss-of-function mutation R577X in the structural muscle protein ACTN3 emerged as a potential target of positive selection from early studies and has been the focus of insightful physiological work suggesting a significant impact on muscle metabolism. Adaptation to cold climates has been proposed as a key adaptive mechanism explaining its global allele frequency patterns. Here, we re-examine this hypothesis analyzing modern (n = 3,626) and ancient (n = 1,651) genomic data by using allele-frequency as well as haplotype homozygosity-based methods. The presented results are more consistent with genetic drift rather than selection in cold climates as the main driver of the ACTN3 R577X frequency distribution in human populations across the world. This Matters Arising paper is in response to Wyckelsma et al. (2021),1 published in The American Journal of Human Genetics. See also the response by Wyckelsma et al. (2022),2 published in this issue.

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.

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.

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.

Mitogenomics of modern Mongolic-speaking populations.

Here, we present a comprehensive data set of 489 complete mitogenomes (211 of which are new) from four Mongolic-speaking populations (Mongols, Barghuts, Khamnigans, and Buryats) to investigate their matrilineal genetic structure, ancestry and relationship with other ethnic groups. We show that along with very high levels of genetic diversity and lack of genetic differentiation, Mongolic-speaking populations exhibit strong genetic resemblance to East Asian populations of Chinese, Japanese, and Uyghurs. Phylogeographic analysis of complete mitogenomes reveals the presence of different components in the gene pools of modern Mongolic-speaking populations-the main East Eurasian component is represented by mtDNA lineages of East Asian, Siberian and autochthonous (the Baikal region/Mongolian) ancestry, whereas the less pronounced West Eurasian component can be ascribed to Europe and West Asia/Caucasus. We also observed that up to one third of the mtDNA subhaplogroups identified in Mongolic-speaking populations can be considered as Mongolic-specific with the coalescence age of most of them not exceeding 1.7 kya. This coincides well with the population size growth which started around 1.1 kya and is detectable only in the Bayesian Skyline Plot constructed based on Mongolic-specific mitogenomes. Our data suggest that the genetic structure established during the Mongol empire is still retained in present-day Mongolic-speaking populations.

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.

Insights into matrilineal genetic structure, differentiation and ancestry of Armenians based on complete mitogenome data.

Distinctive peculiarities of Armenians such as their millennia-long genetic isolation and strong national identity attract a keen interest while studying the demographic history of the West Asia. Here, to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 536 complete mitogenomes (141 of which are novel) from 8 geographically different Armenian populations, covering the whole stretch of historical Armenia. The observed patterns highlight a remarkable degree of matrilineal genetic heterogeneity and weak population structuring of Armenians. Moreover, our phylogeographic analysis reveals common ancestries for some mtDNA lineages shared by West Asians, Transcaucasians, Europeans, Central Asians and Armenians. About third of the mtDNA subhaplogroups found in Armenian gene pool might be considered as Armenian-specific, as these are virtually absent elsewhere in Europe, West Asia and Transcaucasia. Coalescence ages of most of these lineages do not exceed 3.1 kya and coincide well with the population size growth started around 1.8-2.8 kya detectable only in the Bayesian Skyline Plots based on the Armenian-specific mtDNA haplotypes.

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).

Mitogenomic diversity and differentiation of the Buryats.

In this paper we present a results of first comprehensive study of the complete mitogenomes in the Buryats with regard to their belonging to the main regional (eastern and western Buryats); tribal (Khori, Ekhirid, Bulagad, and Khongodor), and ethno-territorial (Aginsk, Alar, Balagansk, Barguzin, Ida, Khorinsk, Kuda, Selenga, Verkholensk, Olkhon, Tunka, and Shenehen Buryats) groups. The analysis of molecular variation performed using regional, tribal, and ethno-territorial divisions of the Buryats showed lack of genetic differentiation at all levels. Nonetheless, the complete mitogenome analysis revealed a very high level of genetic diversity in the Buryats which is the highest among Siberian populations and comparable to that in populations of eastern and western Asia. The AMOVA and MDS analyses results imply to a strong genetic similarity between the Buryats and eastern Asian populations of Chinese and Japanese, suggesting their origin on the basis of common maternal ancestry components. Several new Buryat-specific branches of haplogroup G (G2a2a, G2a1i, G2a5a) display signals of dispersals dating to 2.6-6.6 kya with a possible origin in eastern Asia, thus testifying Bronze Age and Neolithic arrival of ancestral eastern Asian component to the South Siberia region.

The genetic legacy of the expansion of Turkic-speaking nomads across Eurasia.

The Turkic peoples represent a diverse collection of ethnic groups defined by the Turkic languages. These groups have dispersed across a vast area, including Siberia, Northwest China, Central Asia, East Europe, the Caucasus, Anatolia, the Middle East, and Afghanistan. The origin and early dispersal history of the Turkic peoples is disputed, with candidates for their ancient homeland ranging from the Transcaspian steppe to Manchuria in Northeast Asia. Previous genetic studies have not identified a clear-cut unifying genetic signal for the Turkic peoples, which lends support for language replacement rather than demic diffusion as the model for the Turkic language's expansion. We addressed the genetic origin of 373 individuals from 22 Turkic-speaking populations, representing their current geographic range, by analyzing genome-wide high-density genotype data. In agreement with the elite dominance model of language expansion most of the Turkic peoples studied genetically resemble their geographic neighbors. However, western Turkic peoples sampled across West Eurasia shared an excess of long chromosomal tracts that are identical by descent (IBD) with populations from present-day South Siberia and Mongolia (SSM), an area where historians center a series of early Turkic and non-Turkic steppe polities. While SSM matching IBD tracts (> 1cM) are also observed in non-Turkic populations, Turkic peoples demonstrate a higher percentage of such tracts (p-values ≤ 0.01) compared to their non-Turkic neighbors. Finally, we used the ALDER method and inferred admixture dates (~9th-17th centuries) that overlap with the Turkic migrations of the 5th-16th centuries. Thus, our results indicate historical admixture among Turkic peoples, and the recent shared ancestry with modern populations in SSM supports one of the hypothesized homelands for their nomadic Turkic and related Mongolic ancestors.

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

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.

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.

Mitochondrial DNA perspective of Serbian genetic diversity.

Although south-Slavic populations have been studied to date from various aspects, the population of Serbia, occupying the central part of the Balkan Peninsula, is still genetically understudied at least at the level of mitochondrial DNA (mtDNA) variation. We analyzed polymorphisms of the first and the second mtDNA hypervariable segments (HVS-I and HVS-II) and informative coding-region markers in 139 Serbians to shed more light on their mtDNA variability, and used available data on other Slavic and neighboring non-Slavic populations to assess their interrelations in a broader European context. The contemporary Serbian mtDNA profile is consistent with the general European maternal landscape having a substantial proportion of shared haplotypes with eastern, central, and southern European populations. Serbian population was characterized as an important link between easternmost and westernmost south-Slavic populations due to the observed lack of genetic differentiation with all other south-Slavic populations and its geographical positioning within the Balkan Peninsula. An increased heterogeneity of south Slavs, most likely mirroring turbulent demographic events within the Balkan Peninsula over time (i.e., frequent admixture and differential introgression of various gene pools), and a marked geographical stratification of Slavs to south-, east-, and west-Slavic groups, were also found. A phylogeographic analyses of 20 completely sequenced Serbian mitochondrial genomes revealed not only the presence of mtDNA lineages predominantly found within the Slavic gene pool (U4a2a*, U4a2a1, U4a2c, U4a2g, HV10), supporting a common Slavic origin, but also lineages that may have originated within the southern Europe (H5*, H5e1, H5a1v) and the Balkan Peninsula in particular (H6a2b and L2a1k).

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.

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.

Y-chromosome diversity in the Kalmyks at the ethnical and tribal levels.

The Mongolic-speaking Kalmyks currently inhabiting the steppes of the Volga region have Central Asian ancestry and are organized into the tribal groups. The genetic relationships among these tribes and their origin have remained obscure. We analyzed 17 short tandem repeat and 44 binary polymorphisms of Y-chromosome in 426 individuals mainly from three major tribes of the Kalmyks (the Torguuds, Dörwöds and Khoshuuds). Among these tribes, the Dörwöds and Torguuds, as well as the Kalmyks collectively as an ethnic group, showed relatively close genetic affinities to each other and to the Mongols and Altaian Kazakhs, whereas the Khoshuuds were clearly separated from all of them, gathering with the Manchu, Tibetans or Evenks (depending on the algorithm used to calculate genetic distances). The genetic results also indicate that paternal gene flow from East Europeans to the Kalmyks is very little, despite their cohabitation in the North Caspian Steppe during the last 380 years. The occurrence of unique cluster of N1c-Tat haplotypes in the Khoshuuds, which dates to about 340 years and is likely to have East European ancestry, is considered as a result of interethnic contacts occurred soon after the appearance of the Kalmyk tribes in the Volga-Ural region.

Phylogeny and genetic history of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870) inferred from complete mitochondrial genomes.

We assessed phylogeny of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870), the most northern ectothermic, terrestrial vertebrate in Eurasia, by sequence analysis of complete mitochondrial genomes in 26 specimens from different localities (China, Khabarovsk region, Sakhalin, Yakutia, Magadan region, Chukotka, Kamchatka, Ural, European part of Russia). In addition, a complete mitochondrial genome of the Schrenck salamander, Salamandrella schrenckii, was determined for the first time. Bayesian phylogenetic analysis of the entire mtDNA genomes of S. keyserlingii demonstrates that two haplotype clades, AB and C, radiated about 1.4 million years ago (Mya). Bayesian skyline plots of population size change through time show an expansion around 250 thousand years ago (kya) and then a decline around the Last Glacial Maximum (25 kya) with subsequent restoration of population size. Climatic changes during the Quaternary period have dramatically affected the population genetic structure of the Siberian salamanders. In addition, complete mtDNA sequence analysis allowed us to recognize that the vast area of Northern Eurasia was colonized only by the Siberian salamander clade C1b during the last 150 kya. Meanwhile, we were unable to find evidence of molecular adaptation in this clade by analyzing the whole mitochondrial genomes of the Siberian salamanders.