Comparative genomics and the evolution of human mitochondrial DNA: assessing the effects of selection

This article provides evidence that selection has been a significant force during the evolution of the human mitochondrial genome. Both gene-by-gene and whole-genome approaches were used here to assess selection in the 560 mitochondrial DNA (mtDNA) coding-region sequences that were used previously for reduced-median-network analysis. The results of the present analyses were complex, in that the action of selection was not indicated by all tests, but this is not surprising, in view of the characteristics and limitations of the different analytical methods. Despite these limitations, there is evidence for both gene-specific and lineage-specific variation in selection. Whole-genome sliding-window approaches indicated a lack of selection in large-scale segments of the coding region. In other tests, we analyzed the ratio of nonsynonymous-to-synonymous substitutions in the 13 protein-encoding mtDNA genes. The most straightforward interpretation of those results is that negative selection has acted on the mtDNA during evolution. Single-gene analyses indicated significant departures from neutrality in the CO1, ND4, and ND6 genes, although the data also suggested the possible operation of positive selection on the AT6 gene. Finally, our results and those of other investigators do not support a simple model in which climatic adaptation has been a major force during human mtDNA evolution.     

Lineage-specific selection in human mtDNA: lack of polymorphisms in a segment of MTND5 gene in haplogroup J

Human mitochondrial DNA (mtDNA) is a nonrecombining genome that codes for 13 subunits of the mitochondrial oxidative phosphorylation system, 2 rRNAs, and 22 tRNAs. Mutations have accumulated sequentially in mtDNA lineages that diverged tens of thousands of years ago. The genes in mtDNA are subject to different functional constraints and are therefore expected to evolve at different rates, but the rank order of these rates should be the same in all lineages of a phylogeny. Previous studies have indicated, however, that specific regions of mtDNA may have experienced different histories of selection in different lineages, possibly because of lineage-specific interactions or environmental factors such as climate. We report here on a survey for lineage-specific patterns of nucleotide polymorphism in human mtDNA. We calculated molecular polymorphism indices and neutrality tests for classes of functional sites and genes in 837 human mtDNA sequences, compared the results between continent-specific mtDNA lineages, and used two sliding window methods to identify differences in the patterns of polymorphism between haplogroups. A general correlation between nucleotide position and the level of nucleotide polymorphism was identified in the coding region of the mitochondrial genome. Nucleotide diversity in the protein-coding sequence of mtDNA was generally not much higher than that found for many genes in nuclear DNA. A comparison of nonsynonymous/synonymous rate ratios in the 13 protein-coding genes suggested differences in the relative levels of selection between haplogroups, including the European haplogroup clusters. Interestingly, a segment of the MTND5 gene was found to be almost void of segregating sites and nonsynonymous mutations in haplogroup J, which has been associated with susceptibility to certain complex diseases. Our results suggest that there are haplogroup-specific differences in the intensity of selection against particular regions of the mitochondrial genome, indicating that some mutations may be non-neutral within specific phylogenetic lineages but neutral within others.     

mtDNA variation predicts population size in humans and reveals a major Southern Asian chapter in human prehistory

The relative timing and size of regional human population growth following our expansion from Africa remain unknown. Human mitochondrial DNA (mtDNA) diversity carries a legacy of our population history. Given a set of sequences, we can use coalescent theory to estimate past population size through time and draw inferences about human population history. However, recent work has challenged the validity of using mtDNA diversity to infer species population sizes. Here we use Bayesian coalescent inference methods, together with a global data set of 357 human mtDNA coding-region sequences, to infer human population sizes through time across 8 major geographic regions. Our estimates of relative population sizes show remarkable concordance with the contemporary regional distribution of humans across Africa, Eurasia, and the Americas, indicating that mtDNA diversity is a good predictor of population size in humans. Plots of population size through time show slow growth in sub-Saharan Africa beginning 143-193 kya, followed by a rapid expansion into Eurasia after the emergence of the first non-African mtDNA lineages 50-70 kya. Outside Africa, the earliest and fastest growth is inferred in Southern Asia approximately 52 kya, followed by a succession of growth phases in Northern and Central Asia (approximately 49 kya), Australia (approximately 48 kya), Europe (approximately 42 kya), the Middle East and North Africa (approximately 40 kya), New Guinea (approximately 39 kya), the Americas (approximately 18 kya), and a second expansion in Europe (approximately 10-15 kya). Comparisons of relative regional population sizes through time suggest that between approximately 45 and 20 kya most of humanity lived in Southern Asia. These findings not only support the use of mtDNA data for estimating human population size but also provide a unique picture of human prehistory and demonstrate the importance of Southern Asia to our recent evolutionary past.     

Mitochondrial DNA Analysis of Mongolian Populations and Implications for the Origin of New World Founders

High levels of mitochondrial DNA (mtDNA) diversity were determined for Mongolian populations, represented by the Mongol-speaking Khalkha and Dariganga. Although 103 samples were collected across Mongolia, low levels of genetic substructuring were detected, reflecting the nomadic lifestyle and relatively recent ethnic differentiation of Mongolian populations. mtDNA control region I sequence and seven additional mtDNA polymorphisms were assayed to allow extensive comparison with previous human population studies. Based on a comparative analysis, we propose that indigenous populations in east Central Asia represent the closest genetic link between Old and New World populations. Utilizing restriction/deletion polymorphisms, Mongolian populations were found to carry all four New World founding haplogroups as defined by WALLACE and coworkers. The ubiquitous presence of the four New World haplogroups in the Americas but narrow distribution across Asia weakens support for GREENBERG and coworkers' theory of New World colonization via three independent migrations. The statistical and geographic scarcity of New World haplogroups in Asia makes it improbable that the same four haplotypes would be drawn from one geographic region three independent times. Instead, it is likely that founder effects manifest throughout Asia and the Americas are responsible for differences in mtDNA haplotype frequencies observed in these regions.     

Distinguishing the effects of selection from demographic history in the genetic variation of two sister passerines based on mitochondrial–nuclear comparison

Determining the mechanisms responsible for the distribution of genetic diversity in natural populations has occupied a central role in molecular evolution. Our study was motivated by the unprecedented observation that a widespread Eurasian flycatcher, Ficedula albicilla, exhibited no variation at the mitochondrial DNA (mtDNA) ND2 gene in 75 individuals sampled over a 5000-km distance. In contrast, its sister species, F. parva, had low but considerably higher levels of mtDNA variation. We assessed whether natural selection or demographic factors could explain the absence of mtDNA variation in F. albicilla. Eighteen nuclear genes were sequenced to estimate the two species'' phylogeographic histories, and for comparison to the mtDNA data. Multilocus coalescence analyses suggested that F. albicilla experienced a population expansion perhaps following a population bottleneck. Simulations based on this demographic history, however, did not replicate the extremely low level of mtDNA variation. Historical range changes based on ecological niche models also failed to explain the observed mtDNA patterns. Neutrality tests (DHEW and ML-HKA) suggested a non-neutral pattern in the mtDNA of F. albicilla. We found a transmembrane-skewed distribution of nonsynonymous substitutions between the two species, three of which caused functional change; the results implied that positive selection could have targeted mtDNA. Several lines of evidence support selection rather than demographic history as the main force influencing the patterns of mtDNA variation. Despite the influence of natural selection, many of the phylogeographic inferences derived from mtDNA were robust, including species limits and a high level of gene flow among populations within species.     

中国3个不同地域藏族群体线粒体ATP6、ATP8和Cyt b基因的比较: 探查自然选择在基因组的印记

线粒体基因组(Mitochondrial DNA, mtDNA)具有区别于核基因组的诸多特性(例如缺乏重组等), 在人类进化历史的研究领域具有巨大的应用价值。文章通过对3个不同地域藏族群体线粒体基因组的ATP6、ATP8和Cyt b 基因约2 kb区域的测序研究发现: 随着海拔的增加, 3个基因的全序列、ATP6和ATP8基因逐渐偏离中性模式, 但未呈现显著性差异; 随着海拔的降低, Cyt b 基因受到纯化选择的作用逐渐增大; ATP6基因可能存在适应性选择, 并且随着海拔的增高呈现出适应性选择增强的趋势; 选择的主要因素应该是该群体所处的特殊的地理环境, 即不同的地理环境具有直接的选择作用。     

The local origin of the Tibetan pig and additional insights into the origin of Asian pigs

Background

The domestic pig currently indigenous to the Tibetan highlands is supposed to have been introduced during a continuous period of colonization by the ancestors of modern Tibetans. However, there is no direct genetic evidence of either the local origin or exotic migration of the Tibetan pig.

Methods and Findings

We analyzed mtDNA hypervariable segment I (HVI) variation of 218 individuals from seven Tibetan pig populations and 1,737 reported mtDNA sequences from domestic pigs and wild boars across Asia. The Bayesian consensus tree revealed a main haplogroup M and twelve minor haplogroups, which suggested a large number of small scale in situ domestication episodes. In particular, haplogroups D1 and D6 represented two highly divergent lineages in the Tibetan highlands and Island Southeastern Asia, respectively. Network analysis of haplogroup M further revealed one main subhaplogroup M1 and two minor subhaplogroups M2 and M3. Intriguingly, M2 was mainly distributed in Southeastern Asia, suggesting for a local origin. Similar with haplogroup D6, M3 was mainly restricted in Island Southeastern Asia. This pattern suggested that Island Southeastern Asia, but not Southeastern Asia, might be the center of domestication of the so-called Pacific clade (M3 and D6 here) described in previous studies. Diversity gradient analysis of major subhaplogroup M1 suggested three local origins in Southeastern Asia, the middle and downstream regions of the Yangtze River, and the Tibetan highlands, respectively.

Conclusions

We identified two new origin centers for domestic pigs in the Tibetan highlands and in the Island Southeastern Asian region.     

To the problem of differentiation of north Mongoloids,Caucasoids, and Amerinds on the territory of Eurasia (genetic data)

The matrix of genetic distances of 11 human populations of Europe, Asia, and America was calculated using data on frequencies of 28 alleles of 12 loci of proteins, enzymes, and blood groups. This matrix and constructed on their basis dendrogram permit one to suggest that the place of division of paleolitic population on ancestors of North Mongoloids, Caucasoids, and Amerinds was in South Siberia and possibly in neighbouring regions of Central Asia. Published data on polymorphism of mtDNA in the human populations of Eurasia and America were analysed. These data support the author's hypothesis.     

NONRANDOM FACTORS IN MODERN HUMAN MORPHOLOGICAL DIVERSIFICATION: A STUDY OF CRANIOFACIAL VARIATION IN SOUTHERN SOUTH AMERICAN POPULATIONS

The causes of craniofacial variation among human populations have been the subject of controversy. In this work, we studied aboriginal populations from southern South America, the last continental region peopled by humans and with a wide range of ecological conditions. Because of these characteristics, southern South America provides a unique opportunity to study the relative importance of random and nonrandom factors in human diversification. Previous craniometric studies recognized remarkable differences among populations from this region, usually resorting to random factors as the main explanation. In contrast, here we suggest, using tests based on quantitative genetic models, that: (1) the rate of craniofacial divergence among these populations is too high and (2) the patterns of variation within and between populations are too different to be explained by genetic drift alone. In addition, the among-sample craniofacial variation is correlated with climate and diet but not with mtDNA variation. We suggest that the influence of nonrandom factors (e.g., plasticity, selection) on human craniofacial diversification in regions with large ecological variation is more important than generally acknowledged and capable to generate large craniofacial divergence in a short period of time. These results bring nonrandom factors into focus for the interpretation of human craniofacial variation.     

Origins of domestic dog in southern East Asia is supported by analysis of Y-chromosome DNA

Global mitochondrial DNA (mtDNA) data indicates that the dog originates from domestication of wolf in Asia South of Yangtze River (ASY), with minor genetic contributions from dog-wolf hybridisation elsewhere. Archaeological data and autosomal single nucleotide polymorphism data have instead suggested that dogs originate from Europe and/or South West Asia but, because these datasets lack data from ASY, evidence pointing to ASY may have been overlooked. Analyses of additional markers for global datasets, including ASY, are therefore necessary to test if mtDNA phylogeography reflects the actual dog history and not merely stochastic events or selection. Here, we analyse 14,437 bp of Y-chromosome DNA sequence in 151 dogs sampled worldwide. We found 28 haplotypes distributed in five haplogroups. Two haplogroups were universally shared and included three haplotypes carried by 46% of all dogs, but two other haplogroups were primarily restricted to East Asia. Highest genetic diversity and virtually complete phylogenetic coverage was found within ASY. The 151 dogs were estimated to originate from 13-24 wolf founders, but there was no indication of post-domestication dog-wolf hybridisations. Thus, Y-chromosome and mtDNA data give strikingly similar pictures of dog phylogeography, most importantly that roughly 50% of the gene pools are shared universally but only ASY has nearly the full range of genetic diversity, such that the gene pools in all other regions may derive from ASY. This corroborates that ASY was the principal, and possibly sole region of wolf domestication, that a large number of wolves were domesticated, and that subsequent dog-wolf hybridisation contributed modestly to the dog gene pool.     

Genetic data concerning the problem of differentiation of Northern Mongoloids, American Indians and Caucasoids in the northern territory of Eurasia

Basing on the frequencies of 28 alleles of 12 polymorphic loci of blood groups, serum proteins and red cell enzymes the matrix of genetic distances between 11 populations of Europe, Asia and America was calculated. This matrix and the dendrogram based on it permitted to suggest that the region of South Siberia and the neighbouring regions of Central Asia was the place, where the paleolithic populations were divided into the ancestors of the Northern Mongoloids, Caucasoids and American Indians. The published data concerning the human mtDNA polymorphisms support the hypothesis of the author.     

Mitochondrial DNA Content Contributes to Climate Adaptation Using Chinese Populations as a Model

Maintaining a balance between ATP synthesis and heat generation is crucial for adapting to changes in climate. Variation in the mitochondrial DNA (mtDNA), which encodes 13 subunits of the respiratory chain complexes, may contribute to climate adaptation by regulating thermogenesis and the use of bioenergy. However, studies looking for a relationship between mtDNA haplogroups and climate have obtained mixed results, leaving unresolved the role of mtDNA in climate adaptation. Since mtDNA content can regulate human bioenergy processes and is known to influence many physiological traits and diseases, it is possible that mtDNA content contributes to climate adaptation in human populations. Here, we analyze the distribution of mtDNA content among 27 Chinese ethnic populations residing across China and find a significant association between mtDNA content and climate, with northern populations having significantly higher mtDNA content than southern populations. Functional studies have shown that high mtDNA content correlates with an increase in the expression of energy metabolism enzymes, which may accelerate thermogenesis. This suggests that the significantly higher mtDNA content observed in northern populations may confer a selective advantage in adapting to colder northern climates     

Genetic studies of human diversity in East Asia

East Asia is one of the most important regions for studying evolution and genetic diversity of human populations. Recognizing the relevance of characterizing the genetic diversity and structure of East Asian populations for understanding their genetic history and designing and interpreting genetic studies of human diseases, in recent years researchers in China have made substantial efforts to collect samples and generate data especially for markers on Y chromosomes and mtDNA. The hallmark of these efforts is the discovery and confirmation of consistent distinction between northern and southern East Asian populations at genetic markers across the genome. With the confirmation of an African origin for East Asian populations and the observation of a dominating impact of the gene flow entering East Asia from the south in early human settlement, interpretation of the north-south division in this context poses the challenge to the field. Other areas of interest that have been studied include the gene flow between East Asia and its neighbouring regions (i.e. Central Asia, the Sub-continent, America and the Pacific Islands), the origin of Sino-Tibetan populations and expansion of the Chinese.     

MULTILOCUS COALESCENCE ANALYSES SUPPORT A mtDNA-BASED PHYLOGEOGRAPHIC HISTORY FOR A WIDESPREAD PALEARCTIC PASSERINE BIRD, SITTA EUROPAEA

Our understanding of species phylogeography in much of the Palearctic is incomplete. In addition, many existing studies based solely on mitochondrial DNA (mtDNA) can provide a biased view of phylogeographic history because of the effects of lineage sorting, natural selection, or hybridization. We analyzed 13 introns to assess a mtDNA study of the Eurasian nuthatch (Sitta europaea) that suggested a seemingly contemporaneous origin of distinct taxa in the Caucasus, Europe, and Asia. Neutrality tests showed no evidence of selection on either the mtDNA or nuclear sequences. Most nuclear gene trees, except for Z-linked ones, did not recover the three lineages, which we attribute to recent splitting. Analyses of the 13 introns combined revealed the same three groups as did the mtDNA and suggested that nuthatches experienced a trichotomous (or two indistinguishable) split(s) 1-2 million years ago (Mya) and have remained isolated with trifling if not zero gene flow since then, and the Asian group increased in population size. This result demonstrates the usefulness of mtDNA in discovering phylogeographic patterns. The use of multiple nuclear loci facilitated detection of an introgressed individual and improved estimates of process parameters such as divergence time and population expansion. We recommend that phylogeographic studies should be based on both mtDNA and nuclear genes.     

Climate shaped the worldwide distribution of human mitochondrial DNA sequence variation

There is an ongoing discussion in the literature on whether human mitochondrial DNA (mtDNA) evolves neutrally. There have been previous claims for natural selection on human mtDNA based on an excess of non-synonymous mutations and higher evolutionary persistence of specific mitochondrial mutations in Arctic populations. However, these findings were not supported by the reanalysis of larger datasets. Using a geographical framework, we perform the first direct test of the relative extent to which climate and past demography have shaped the current spatial distribution of mtDNA sequences worldwide. We show that populations living in colder environments have lower mitochondrial diversity and that the genetic differentiation between pairs of populations correlates with difference in temperature. These associations were unique to mtDNA; we could not find a similar pattern in any other genetic marker. We were able to identify two correlated non-synonymous point mutations in the ND3 and ATP6 genes characterized by a clear association with temperature, which appear to be plausible targets of natural selection producing the association with climate. The same mutations have been previously shown to be associated with variation in mitochondrial pH and calcium dynamics. Our results indicate that natural selection mediated by climate has contributed to shape the current distribution of mtDNA sequences in humans.     

Mitochondrial diversity and phylogeographic structure of Chinese domestic goats

China has numerous native domestic goat breeds, but so far there has been no extensive study on genetic diversity, population demographic history, and origin of Chinese goats. Here, we examined the genetic diversity and phylogeographic structure of Chinese domestic goats by determining a 481-bp fragment of the first hypervariable region of mitochondrial DNA (mtDNA) control region from 368 individuals representing 18 indigenous breeds. Phylogenetic analyses revealed that there were four mtDNA lineages (A-D) identified in Chinese goats, in which lineage A was predominant, lineage B was moderate, and lineages C and D were at low frequency. These results further support the multiple maternal origins of domestic goats. The pattern of genetic variation in goat mtDNA sequences indicated that the two larger lineages A and B had undergone population expansion events. In a combined analysis of previously reported sequences and our sequences belonging to lineage B, we detected two subclades, in which one was unique to eastern Asia and another was shared between eastern and southern Asia. A larger genetic variation in eastern Asia than southern Asia and the pattern of phylogeographic variation in lineage B suggest that at least one subclade of lineage B originated from eastern Asia. There was no significant geographical structuring in Chinese goat populations, which suggested that there existed strong gene flow among goat populations caused by extensive transportation of goats in history.     

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans,but not in chimpanzees.

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA''s D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.     

Multiple maternal origins of chickens: out of the Asian jungles

Domestic chickens have long been important to human societies for food, religion, entertainment, and decorative uses, yet the origins and phylogeography of chickens through Eurasia remain uncertain. Here, we assessed their origins and phylogeographic history by analyzing the mitochondrial DNA hypervariable segment I (HVS-I) for 834 domestic chickens (Gallus gallus domesticus) across Eurasia as well as 66 wild red jungle fowls (Gallus gallus) from Southeast Asia and China. Phylogenetic analyses revealed nine highly divergent mtDNA clades (A-I) in which seven clades contained both the red jungle fowls and domestic chickens. There was no breed-specific clade in the chickens. The clades A, B, and E are distributed ubiquitously in Eurasia, while the other clades were restricted to South and Southeast Asia. Clade C was mainly distributed in Japan and Southeast China, while clades F and G were exclusive to Yunnan, China. The geographic distribution of clade D was closely related to the distribution of the pastime of cock fighting. Statistical tests detect population expansion within each subclade. These distinct distribution patterns and expansion signatures suggest that different clades may originate from different regions, such as Yunnan, South and Southwest China and/or surrounding areas (i.e., Vietnam, Burma, and Thailand), and the Indian subcontinent, respectively, which support the theory of multiple origins in South and Southeast Asia.     

Comprehensive study of mtDNA among Southwest Asian dogs contradicts independent domestication of wolf, but implies dog-wolf hybridization

Studies of mitochondrial DNA (mtDNA) diversity indicate explicitly that dogs were domesticated, probably exclusively, in southern East Asia. However, Southwest Asia (SwAsia) has had poor representation and geographical coverage in these studies. Other studies based on archaeological and genome-wide SNP data have suggested an origin of dogs in SwAsia. Hence, it has been suspected that mtDNA evidence for this scenario may have remained undetected. In the first comprehensive investigation of genetic diversity among SwAsian dogs, we analyzed 582 bp of mtDNA for 345 indigenous dogs from across SwAsia, and compared with 1556 dogs across the Old World. We show that 97.4% of SwAsian dogs carry haplotypes belonging to a universal mtDNA gene pool, but that only a subset of this pool, five of the 10 principal haplogroups, is represented in SwAsia. A high frequency of haplogroup B, potentially signifying a local origin, was not paralleled with the high genetic diversity expected for a center of origin. Meanwhile, 2.6% of the SwAsian dogs carried the rare non-universal haplogroup d2. Thus, mtDNA data give no indication that dogs originated in SwAsia through independent domestication of wolf, but dog-wolf hybridization may have formed the local haplogroup d2 within this region. Southern East Asia remains the only region with virtually full extent of genetic variation, strongly indicating it to be the primary and probably sole center of wolf domestication. An origin of dogs in southern East Asia may have been overlooked by other studies due to a substantial lack of samples from this region.     

Unexpected island effects at an extreme: reduced Y chromosome and mitochondrial DNA diversity in Nias

The amount of genetic diversity in a population is determined by demographic and selection events in its history. Human populations which exhibit greatly reduced overall genetic diversity, presumably resulting from severe bottlenecks or founder events, are particularly interesting, not least because of their potential to serve as valuable resources for health studies. Here, we present an unexpected case, the human population of Nias Island in Indonesia, that exhibits severely reduced Y chromosome (non-recombining portion of the Y chromosome [NRY]) and to a lesser extent also reduced mitochondrial DNA (mtDNA) diversity as compared with most other populations from the Asia/Oceania region. Our genetic data, collected from more than 400 individuals from across the island, suggest a strong previously undetected bottleneck or founder event in the human population history of Nias, more pronounced for males than for females, followed by subsequent genetic isolation. Our findings are unexpected given the island's geographic proximity to the genetically highly diverse Southeast Asian world, as well as our previous knowledge about the human history of Nias. Furthermore, all NRY and virtually all mtDNA haplogroups observed in Nias can be attributed to the Austronesian expansion, in line with linguistic data, and in contrast with archaeological evidence for a pre-Austronesian occupation of Nias that, as we show here, left no significant genetic footprints in the contemporary population. Our work underlines the importance of human genetic diversity studies not only for a better understanding of human population history but also because of the potential relevance for genetic disease-mapping studies.