Adaptive evolution signals in mitochondrial genes of Europeans

Evaluating the role of selection in the evolution of mitochondrial genomes of aboriginal peoples of Siberia

Studies of the nature of mitochondrial DNA (mtDNA) variability in human populations have shown that protein-coding genes are under negative (purifying) selection, since their mutation spectra are characterized by a pronounced predominance of synonymous substitutions over non-synonymous ones (Ka/Ks < 1). Meanwhile, a number of studies have shown that the adaptation of populations to various environmental conditions may be accompanied by a relaxation of negative selection in some mtDNA genes. For example, it was previously found that in Arctic populations, negative selection is relaxed in the mitochondrial ATP6 gene, which encodes one of the subunits of ATP synthase. In this work, we performed a Ka/Ks analysis of mitochondrial genes in large samples of three regional population groups in Eurasia: Siberia (N = 803), Western Asia/Transcaucasia (N = 753), and Eastern Europe (N = 707). The main goal of this work is to search for traces of adaptive evolution in the mtDNA genes of aboriginal peoples of Siberia represented by populations of the north (Koryaks, Evens) and the south of Siberia and the adjacent territory of Northeast China (Buryats, Barghuts, Khamnigans). Using standard Ka/Ks analysis, it was found that all mtDNA genes in all studied regional population groups are subject to negative selection. The highest Ka/Ks values in different regional samples were found in almost the same set of genes encoding subunits of ATP synthase (ATP6, ATP8), NADH dehydrogenase complex (ND1, ND2, ND3), and cytochrome bc1 complex (CYB). The highest Ka/Ks value, indicating a relaxation of negative selection, was found in the ATP6 gene in the Siberian group. The results of the analysis performed using the FUBAR method (HyPhy software package) and aimed at searching for mtDNA codons under the influence of selection also showed the predominance of negative selection over positive selection in all population groups. In Siberian populations, nucleotide sites that are under positive selection and associated with mtDNA haplogroups were registered not in the north (which is expected under the assumption of adaptive evolution of mtDNA), but in the south of Siberia.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.