The role of selection in the evolution of human mitochondrial genomes

A disproportionate role for mtDNA in Dobzhansky-Muller incompatibilities?

Evolution in allopatric populations can lead to incompatibilities that result in reduced hybrid fitness and ultimately reproductive isolation upon secondary contact. The Dobzhansky-Muller (DM) model nicely accounts for the evolution of such incompatibilities. Although DM incompatibilities were originally conceived as resulting of interactions between nuclear genes, recent studies have documented cases where incompatibilities have arisen between nuclear and mitochondrial genomes (mtDNA). Although mtDNA comprises only a tiny component (typically <<0.01%) of an organism's genetic material, several features of mtDNA may lead to a disproportionate contribution to the evolution of hybrid incompatibilities: (i) essentially all functions of mtDNA require interaction with nuclear gene products. All mtDNA-encoded proteins are components of the oxidative phosphorylation (OXPHOS) system and all mtDNA-encoded RNAs are part of the mitochondrial protein synthetic machinery; both processes require interaction with nuclear-encoded proteins for function. (ii) Transcription and replication of mtDNA also involve mitonuclear interactions as nuclear-encoded proteins must bind to regulatory motifs in the mtDNA to initiate these processes. (iii) Although features of mtDNA vary amongst taxa, metazoan mtDNA is typically characterized by high nucleotide substitution rates, lack of recombination and reduced effective population sizes that collectively lead to increased chance fixation of mildly deleterious mutations. Combined, these features create an evolutionary dynamic where rapid mtDNA evolution favours compensatory nuclear gene evolution, ultimately leading to co-adaptation of mitochondrial and nuclear genomes. When previously isolated lineages hybridize in nature or in the lab, intergenomic co-adaptation is disrupted and hybrid breakdown is observed; the role of intergenomic co-adaptation in hybrid breakdown and speciation will generally be most pronounced when rates of mtDNA evolution are high or when restricted gene flow results in significant population differentiation.

Adaptive evolution of the FADS gene cluster within Africa

Long chain polyunsaturated fatty acids (LC-PUFAs) are essential for brain structure, development, and function, and adequate dietary quantities of LC-PUFAs are thought to have been necessary for both brain expansion and the increase in brain complexity observed during modern human evolution. Previous studies conducted in largely European populations suggest that humans have limited capacity to synthesize brain LC-PUFAs such as docosahexaenoic acid (DHA) from plant-based medium chain (MC) PUFAs due to limited desaturase activity. Population-based differences in LC-PUFA levels and their product-to-substrate ratios can, in part, be explained by polymorphisms in the fatty acid desaturase (FADS) gene cluster, which have been associated with increased conversion of MC-PUFAs to LC-PUFAs. Here, we show evidence that these high efficiency converter alleles in the FADS gene cluster were likely driven to near fixation in African populations by positive selection ∼85 kya. We hypothesize that selection at FADS variants, which increase LC-PUFA synthesis from plant-based MC-PUFAs, played an important role in allowing African populations obligatorily tethered to marine sources for LC-PUFAs in isolated geographic regions, to rapidly expand throughout the African continent 60-80 kya.

Results of Mitochondrial DNA Sequence Analysis in Patients with Clinically Diagnosed Leber's Hereditary Optic Neuropathy

OBJECTIVE

To investigate possible mitochondrial DNA (mtDNA) mutations in patients with Leber's hereditary optic neuropathy (LHON) in order to provide a precise diagnosis and genetic counseling.

MATERIAL AND METHODS

Between 1982 and 2007, ten patients were clinically diagnosed with LHON and six of these patients agreed to be involved in this study. Six healthy individuals were also included as a control group. mtDNA was isolated from peripheral blood samples and polymerase chain reaction and mtDNA sequence analysis were performed.

RESULTS

In one of the six patients, a homoplasmic mutant m.11778G>A mutation was detected. All of the clinically diagnosed LHON patients and the control groups had the m.14212C>T and m.14580G>A single nucleotide polymorphisms (SNPs). The m.11719A>G SNP was detected in three of six patients and four of the controls. Two of the six patients had the m.3197T>C SNP and, in addition, the m.14258G>A SNP was found in one of these two patients, while neither of these mutations were present in the control group.

CONCLUSION

The clinical diagnosis of LHON could be supported by molecular genetics only in one patient by the detection of one mutation. The m.3197T>C and m.14258G>A SNPs should be considered as potential mtDNA mutations due to the fact that they were detected in the patient group. These mutations should be investigated further in large case groups for suspected gene loci that could lead to optic neuropathy.

Mitochondrial DNA mutation screening in an ethnically diverse nonsyndromic deafness cohort

Deafness is a heterogeneous trait with many known genetic and environmental causes. Hereditary hearing loss is an extremely common disorder in the general population. Mutations in mitochondrial DNA (mtDNA) are known to be associated with nonsyndromic deafness (NSD) and syndromic deafness. The objective of this article is to investigate the frequency of common mitochondrial mutations (A1555G, G7444A, and A3243G) in an ethnically diverse cohort of probands with NSD from South Florida. These patients were ascertained at the University of Miami. Polymerase chain reaction-restriction fragment length polymorphism analysis and direct sequencing methods were used for mutation screening in a cohort of 217 patients with NSD. The frequency of common mitochondrial mutations is 1.84% (4/217) in this cohort. A1555G and G7444A accounted for four patients with NSD. Our mutation frequencies are comparable with those previously reported in other populations, indicating that mutations in mtDNA are an important cause of NSD in our patient cohort.

Did a discrete event 200,000-100,000 years ago produce modern humans?

Scenarios for modern human origins are often predicated on the assumption that modern humans arose 200,000-100,000 years ago in Africa. This assumption implies that something 'special' happened at this point in time in Africa, such as the speciation that produced Homo sapiens, a severe bottleneck in human population size, or a combination of the two. The common thread is that after the divergence of the modern human and Neandertal evolutionary lineages ∼400,000 years ago, there was another discrete event near in time to the Middle-Late Pleistocene boundary that produced modern humans. Alternatively, modern human origins could have been a lengthy process that lasted from the divergence of the modern human and Neandertal evolutionary lineages to the expansion of modern humans out of Africa, and nothing out of the ordinary happened 200,000-100,000 years ago in Africa. Three pieces of biological (fossil morphology and DNA sequences) evidence are typically cited in support of discrete event models. First, living human mitochondrial DNA haplotypes coalesce ∼200,000 years ago. Second, fossil specimens that are usually classified as 'anatomically modern' seem to appear shortly afterward in the African fossil record. Third, it is argued that these anatomically modern fossils are morphologically quite different from the fossils that preceded them. Here I use theory from population and quantitative genetics to show that lengthy process models are also consistent with current biological evidence. That this class of models is a viable option has implications for how modern human origins is conceptualized.

Characterization of mtDNA variation in a cohort of South African paediatric patients with mitochondrial disease

Mitochondrial disease can be attributed to both mitochondrial and nuclear gene mutations. It has a heterogeneous clinical and biochemical profile, which is compounded by the diversity of the genetic background. Disease-based epidemiological information has expanded significantly in recent decades, but little information is known that clarifies the aetiology in African patients. The aim of this study was to investigate mitochondrial DNA variation and pathogenic mutations in the muscle of diagnosed paediatric patients from South Africa. A cohort of 71 South African paediatric patients was included and a high-throughput nucleotide sequencing approach was used to sequence full-length muscle mtDNA. The average coverage of the mtDNA genome was 81±26 per position. After assigning haplogroups, it was determined that although the nature of non-haplogroup-defining variants was similar in African and non-African haplogroup patients, the number of substitutions were significantly higher in African patients. We describe previously reported disease-associated and novel variants in this cohort. We observed a general lack of commonly reported syndrome-associated mutations, which supports clinical observations and confirms general observations in African patients when using single mutation screening strategies based on (predominantly non-African) mtDNA disease-based information. It is finally concluded that this first extensive report on muscle mtDNA sequences in African paediatric patients highlights the need for a full-length mtDNA sequencing strategy, which applies to all populations where specific mutations is not present. This, in addition to nuclear DNA gene mutation and pathogenicity evaluations, will be required to better unravel the aetiology of these disorders in African patients.

Impact of fertility transmission and other sociodemographic factors on reproductive success and coalescent trees

Summary Fertility transmission (FT) is a phenomenon with a cultural and/or genetic basis, whereby a positive correlation exists between the number of offspring of an individual and that of his/her parents. Theoretical studies using a haploid individual-based model have shown that FT increases the variance and intergenerational correlation in reproductive success and results in an imbalance in the coalescent tree of sampled genes. This phenomenon has been documented in several demographic studies conducted on the correlation in fertility between generations, or through the reconstruction of the genealogical trees of mitochondrial DNA sequences. However, as mtDNA is a single locus, potentially subject to other forces (e.g. natural selection), it is of interest to extend the theory of FT to nuclear loci. We show that because random mating between individuals leads to a mixing of their fertility profiles, FT in these cases will have less influence on the variance and intergenerational correlation of reproductive success. This, in turn, results in less impact on the shape of the coalescent trees. Nevertheless, in the presence of FT, high heterogeneity in reproductive success and homogamy for family size will increase the imbalance in the coalescent tree. Thus, FT should be easier to detect when occurring in conjunction with these other factors. We also show the utility of analysing different kinds of loci (X-linked, Y-linked, mitochondrial and autosomal) to assess whether FT is matrilineal, patrilineal or biparental. Finally, we demonstrate that the shape of the coalescent tree depends upon population size, in contrast to the classical Kingman's model.

Uniparental genetic markers in South Amerindians

A comprehensive review of uniparental systems in South Amerindians was undertaken. Variability in the Y-chromosome haplogroups were assessed in 68 populations and 1,814 individuals whereas that of Y-STR markers was assessed in 29 populations and 590 subjects. Variability in the mitochondrial DNA (mtDNA) haplogroup was examined in 108 populations and 6,697 persons, and sequencing studies used either the complete mtDNA genome or the highly variable segments 1 and 2. The diversity of the markers made it difficult to establish a general picture of Y-chromosome variability in the populations studied. However, haplogroup Q1a3a* was almost always the most prevalent whereas Q1a3* occurred equally in all regions, which suggested its prevalence among the early colonizers. The STR allele frequencies were used to derive a possible ancient Native American Q-clade chromosome haplotype and five of six STR loci showed significant geographic variation. Geographic and linguistic factors moderately influenced the mtDNA distributions (6% and 7%, respectively) and mtDNA haplogroups A and D correlated positively and negatively, respectively, with latitude. The data analyzed here provide rich material for understanding the biological history of South Amerindians and can serve as a basis for comparative studies involving other types of data, such as cultural data.

Complete mitochondrial genome sequencing reveals novel haplotypes in a Polynesian population

The high risk of metabolic disease traits in Polynesians may be partly explained by elevated prevalence of genetic variants involved in energy metabolism. The genetics of Polynesian populations has been shaped by island hoping migration events which have possibly favoured thrifty genes. The aim of this study was to sequence the mitochondrial genome in a group of Maoris in an effort to characterise genome variation in this Polynesian population for use in future disease association studies. We sequenced the complete mitochondrial genomes of 20 non-admixed Maori subjects using Affymetrix technology. DNA diversity analyses showed the Maori group exhibited reduced mitochondrial genome diversity compared to other worldwide populations, which is consistent with historical bottleneck and founder effects. Global phylogenetic analysis positioned these Maori subjects specifically within mitochondrial haplogroup - B4a1a1. Interestingly, we identified several novel variants that collectively form new and unique Maori motifs – B4a1a1c, B4a1a1a3 and B4a1a1a5. Compared to ancestral populations we observed an increased frequency of non-synonymous coding variants of several mitochondrial genes in the Maori group, which may be a result of positive selection and/or genetic drift effects. In conclusion, this study reports the first complete mitochondrial genome sequence data for a Maori population. Overall, these new data reveal novel mitochondrial genome signatures in this Polynesian population and enhance the phylogenetic picture of maternal ancestry in Oceania. The increased frequency of several mitochondrial coding variants makes them good candidates for future studies aimed at assessment of metabolic disease risk in Polynesian populations.

A "Copernican" reassessment of the human mitochondrial DNA tree from its root

Mutational events along the human mtDNA phylogeny are traditionally identified relative to the revised Cambridge Reference Sequence, a contemporary European sequence published in 1981. This historical choice is a continuous source of inconsistencies, misinterpretations, and errors in medical, forensic, and population genetic studies. Here, after having refined the human mtDNA phylogeny to an unprecedented level by adding information from 8,216 modern mitogenomes, we propose switching the reference to a Reconstructed Sapiens Reference Sequence, which was identified by considering all available mitogenomes from Homo neanderthalensis. This "Copernican" reassessment of the human mtDNA tree from its deepest root should resolve previous problems and will have a substantial practical and educational influence on the scientific and public perception of human evolution by clarifying the core principles of common ancestry for extant descendants.

High mitochondrial mutation rates estimated from deep-rooting Costa Rican pedigrees

Estimates of mutation rates for the noncoding hypervariable Region I (HVR-I) of mitochondrial DNA vary widely, depending on whether they are inferred from phylogenies (assuming that molecular evolution is clock-like) or directly from pedigrees. All pedigree-based studies so far were conducted on populations of European origin. In this article, we analyzed 19 deep-rooting pedigrees in a population of mixed origin in Costa Rica. We calculated two estimates of the HVR-I mutation rate, one considering all apparent mutations, and one disregarding changes at sites known to be mutational hot spots and eliminating genealogy branches which might be suspected to include errors, or unrecognized adoptions along the female lines. At the end of this procedure, we still observed a mutation rate equal to 1.24 × 10(-6) , per site per year, i.e., at least threefold as high as estimates derived from phylogenies. Our results confirm that mutation rates observed in pedigrees are much higher than estimated assuming a neutral model of long-term HVRI evolution. We argue that until the cause of these discrepancies will be fully understood, both lower estimates (i.e., those derived from phylogenetic comparisons) and higher, direct estimates such as those obtained in this study, should be considered when modeling evolutionary and demographic processes.

Revisiting the role of the Himalayas in peopling Nepal: insights from mitochondrial genomes

Himalayas was believed to be a formidably geographical barrier between South and East Asia. The observed high frequency of the East Eurasian paternal lineages in Nepal led some researchers to suggest that these lineages were introduced into Nepal from Tibet directly; however, it is also possible that the East Eurasian genetic components might trace their origins to northeast India where abundant East Eurasian maternal lineages have been detected. To trace the origin of the Nepalese maternal genetic components, especially those of East Eurasian ancestry, and then to better understand the role of the Himalayas in peopling Nepal, we have studied the matenal genetic composition extensively, especially the East Eurasian lineages, in Nepalese and its surrounding populations. Our results revealed the closer affinity between the Nepalese and the Tibetans, specifically, the Nepalese lineages of the East Eurasian ancestry generally are phylogenetically closer with the ones from Tibet, albeit a few mitochondrial DNA haplotypes, likely resulted from recent gene flow, were shared between the Nepalese and northeast Indians. It seems that Tibet was most likely to be the homeland for most of the East Eurasian in the Nepalese. Taking into account the previous observation on Y chromosome, now it is convincing that bearer of the East Eurasian genetic components had entered Nepal across the Himalayas around 6 kilo years ago (kya), a scenario in good agreement with the previous results from linguistics and archeology.

How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.

Association of mitochondrial DNA variations with lung cancer risk in a Han Chinese population from southwestern China

Mitochondrial DNA (mtDNA) is particularly susceptible to oxidative damage and mutation due to the high rate of reactive oxygen species (ROS) production and limited DNA-repair capacity in mitochondrial. Previous studies demonstrated that the increased mtDNA copy number for compensation for damage, which was associated with cigarette smoking, has been found to be associated with lung cancer risk among heavy smokers. Given that the common and “non-pathological” mtDNA variations determine differences in oxidative phosphorylation performance and ROS production, an important determinant of lung cancer risk, we hypothesize that the mtDNA variations may play roles in lung cancer risk. To test this hypothesis, we conducted a case-control study to compare the frequencies of mtDNA haplogroups and an 822 bp mtDNA deletion between 422 lung cancer patients and 504 controls. Multivariate logistic regression analysis revealed that haplogroups D and F were related to individual lung cancer resistance (OR = 0.465, 95%CI = 0.329–0.656, p<0.001; and OR = 0.622, 95%CI = 0.425–0.909, p = 0.014, respectively), while haplogroups G and M7 might be risk factors for lung cancer (OR = 3.924, 95%CI = 1.757–6.689, p<0.001; and OR = 2.037, 95%CI = 1.253–3.312, p = 0.004, respectively). Additionally, multivariate logistic regression analysis revealed that cigarette smoking was a risk factor for the 822 bp mtDNA deletion. Furthermore, the increased frequencies of the mtDNA deletion in male cigarette smoking subjects of combined cases and controls with haplogroup D indicated that the haplogroup D might be susceptible to DNA damage from external ROS caused by heavy cigarette smoking.

Somatic mitochondrial DNA mutations in cancer escape purifying selection and high pathogenicity mutations lead to the oncocytic phenotype: pathogenicity analysis of reported somatic mtDNA mutations in tumors

BACKGROUND

The presence of somatic mitochondrial DNA (mtDNA) mutations in cancer cells has been interpreted in controversial ways, ranging from random neutral accumulation of mutations, to positive selection for high pathogenicity, or conversely to purifying selection against high pathogenicity variants as occurs at the population level.

METHODS

Here we evaluated the predicted pathogenicity of somatic mtDNA mutations described in cancer and compare these to the distribution of variations observed in the global human population and all possible protein variations that could occur in human mtDNA. We focus on oncocytic tumors, which are clearly associated with mitochondrial dysfunction. The protein variant pathogenicity was predicted using two computational methods, MutPred and SNPs&GO.

RESULTS

The pathogenicity score of the somatic mtDNA variants were significantly higher in oncocytic tumors compared to non-oncocytic tumors. Variations in subunits of Complex I of the electron transfer chain were significantly more common in tumors with the oncocytic phenotype, while variations in Complex V subunits were significantly more common in non-oncocytic tumors.

CONCLUSIONS

Our results show that the somatic mtDNA mutations reported over all tumors are indistinguishable from a random selection from the set of all possible amino acid variations, and have therefore escaped the effects of purifying selection that act strongly at the population level. We show that the pathogenicity of somatic mtDNA mutations is a determining factor for the oncocytic phenotype. The opposite associations of the Complex I and Complex V variants with the oncocytic and non-oncocytic tumors implies that low mitochondrial membrane potential may play an important role in determining the oncocytic phenotype.

Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication

Archaeological and genetic evidence concerning the time and mode of wild horse (Equus ferus) domestication is still debated. High levels of genetic diversity in horse mtDNA have been detected when analyzing the control region; recurrent mutations, however, tend to blur the structure of the phylogenetic tree. Here, we brought the horse mtDNA phylogeny to the highest level of molecular resolution by analyzing 83 mitochondrial genomes from modern horses across Asia, Europe, the Middle East, and the Americas. Our data reveal 18 major haplogroups (A-R) with radiation times that are mostly confined to the Neolithic and later periods and place the root of the phylogeny corresponding to the Ancestral Mare Mitogenome at ~130-160 thousand years ago. All haplogroups were detected in modern horses from Asia, but F was only found in E. przewalskii--the only remaining wild horse. Therefore, a wide range of matrilineal lineages from the extinct E. ferus underwent domestication in the Eurasian steppes during the Eneolithic period and were transmitted to modern E. caballus breeds. Importantly, now that the major horse haplogroups have been defined, each with diagnostic mutational motifs (in both the coding and control regions), these haplotypes could be easily used to (i) classify well-preserved ancient remains, (ii) (re)assess the haplogroup variation of modern breeds, including Thoroughbreds, and (iii) evaluate the possible role of mtDNA backgrounds in racehorse performance.

The abundance of deleterious polymorphisms in humans

Here I show a gradual decline in the proportion of deleterious nonsynonymous SNPs (nSNPs) from tip to root of the human population tree. This study reveals that up to 48% of nSNPs specific to a single genome are deleterious in nature, which underscores the abundance of deleterious polymorphisms in humans.

Genetic uniqueness of the Waorani tribe from the Ecuadorian Amazon

South America and especially the Amazon basin is known to be home to some of the most isolated human groups in the world. Here, we report on a study of mitochondrial DNA (mtDNA) in the Waorani from Ecuador, probably the most warlike human population known to date. Seeking to look in more depth at the characterization of the genetic diversity of this Native American tribe, molecular markers from the X and Y chromosomes were also analyzed. Only three different mtDNA haplotypes were detected among the Waorani sample. One of them, assigned to Native American haplogroup A2, accounted for more than 94% of the total diversity of the maternal gene pool. Our results for sex chromosome molecular markers failed to find close genetic kinship between individuals, further emphasizing the low genetic diversity of the mtDNA. Bearing in mind the results obtained for both the analysis of the mtDNA control region and complete mitochondrial genomes, we suggest the existence of a 'Waorani-specific' mtDNA lineage. According to current knowledge on the phylogeny of haplogroup A2, we propose that this lineage could be designated as subhaplogroup A2s. Its wide predominance among the Waorani people might have been conditioned by severe genetic drift episodes resulting from founding events, long-term isolation and a traditionally small population size most likely associated with the striking ethnography of this Amazonian community. In all, the Waorani constitute a fine example of how genetic imprint may mirror ethnopsychology and sociocultural features in human populations.

Mitochondrial Cytochrome c Oxidase subunit 1 Sequence Variation in Prostate Cancer

PURPOSE

Mitochondrial DNA (mtDNA) gene mutations have been described in nearly every adult solid neoplasm including prostate cancer. There are marked racial differences in specific inherited mutations within the cytochrome c oxidase subunit 1 (COI) gene in individuals with prostate cancer (PCa). The purpose of this study was to identify the variation in COI gene sequence in prostate cancer patients and to compare the mutations in African and Caucasian Americans.

MATERIALS AND METHODS

We sequenced the COI gene in DNA derived from peripheral blood in 482 prostate cancer patients and 189 controls. All bases that differed from the revised Cambridge Reference Sequence (rCRS) were classified as either silent (non-amino acid altering) or missense (amino acid altering) and the compiled alterations were then compared between races and published reports of mutations in this gene in both Caucasian and African-Americans.

RESULTS AND CONCLUSIONS

We found inherited mtDNA COI missense variants in 8.8% of Caucasian prostate cancer patients (vs. 0.0% controls) and 72.8 % of African-American prostate cancer patients (vs. 64.3% controls) A total of 144 COI variants were identified, of which 30 were missense mutations. Of 482 PCa patients, 116 (24.1%) had one or more missense mutations. Further evaluation of this gene and these mutations may allow for the identification of genetically at-risk populations. The high rate of COI mutations in African-Americans may account for some of the racial disparity observed in prostate cancer.

Evolution of the couple cytochrome c and cytochrome c oxidase in primates

Mitochondrial energy metabolism has been affected by a broad set of ancient and recent evolutionary events. The oldest example is the endosymbiosis theory that led to mitochondria and a recently proposed example is adaptation to cold climate by anatomically modern human lineages. Mitochondrial energy metabolism has also been associated with an important area in anthropology and evolutionary biology, brain enlargement in human evolution. Indeed, several studies have pointed to the need for a major metabolic rearrangement to supply a sufficient amount of energy for brain development in primates.The genes encoding for the coupled cytochrome c (Cyt c) and cytochrome c oxidase (COX, complex IV, EC 1.9.3.1) seem to have an exceptional pattern of evolution in the anthropoid lineage. It has been proposed that this evolution was linked to the rearrangement of energy metabolism needed for brain enlargement. This hypothesis is reinforced by the fact that the COX enzyme was proposed to have a large role in control of the respiratory chain and thereby global energy production.After summarizing major events that occurred during the evolution of COX and cytochrome c on the primate lineage, we review the different evolutionary forces that could have influenced primate COX evolution and discuss the probable causes and consequences of this evolution. Finally, we discuss and review the co-occurring primate phenotypic evolution.

Mitochondria in anthropology and forensic medicine

Mitochondria's role in crucial metabolic pathways is probably the first answer which comes to our minds for the question: what do these tiny organelles serve for? However, specific features of their DNA made them extremely useful also in the field of anthropology and forensics. MtDNA analyses became a milestone in the complex task of unraveling earliest human migrations. Evidence provided by these experiments left no doubts on modern humans origins pointing to Africa being our cradle. It also contributed to interpretation of putative ways of our dispersal around Asia and Americas thousands years ago. On the other hand, analysis of mtDNA is well established and valuable tool in forensic genetics. When other definitely more popular markers give no answer on identity, it is the time to employ information carried by mitochondria. This chapter summarizes not only current reports on the role of mitochondria in forensics and reconstruction of modern humans phylogeny, but also calls one's attention to a broad range of difficulties and constraints associated with mtDNA analyses.

Differences in mtDNA whole sequence between Tibetan and Han populations suggesting adaptive selection to high altitude

We performed a mitochondrial whole-genome comparison study in 40 Tibetan and 50 Han Chinese. All subjects could be classified into 13 haplogroups pertained to the Macrohaplogroup M and N that pitched different quadrants by principal component analysis. We observed a difference in the M9 haplogroup and identified 18 significant variants by comparing whole sequences between Tibetan and Han populations. Variants in ND2, COX2, tRNA alanine and 12S rRNA were predicted to confer increased protein stability in Tibetans. We compared the base substitutions of nonsynonymous (NS) versus synonymous (S) of 13 protein-encoding genes and found the NS/S values of the ATP6, ATP8, and Cyt b genes were larger (>1) in Tibetans than that in Han population. Our findings provide clues for the existence of adaptive selection for the ATP6, ATP8, Cyt b, ND2, COX2, tRNA alanine and 12S rRNA genes in Tibetans which likely contributed to adaptation to their specific geographic environment, such as high altitude.

Native Americans experienced a strong population bottleneck coincident with European contact

The genetic and demographic impact of European contact with Native Americans has remained unclear despite recent interest. Whereas archeological and historical records indicate that European contact resulted in widespread mortality from various sources, genetic studies have found little evidence of a recent contraction in Native American population size. In this study we use a large dataset including both ancient and contemporary mitochondrial DNA to construct a high-resolution portrait of the Holocene and late Pleistocene population size of indigenous Americans. Our reconstruction suggests that Native Americans suffered a significant, although transient, contraction in population size some 500 y before the present, during which female effective size was reduced by ∼50%. These results support analyses of historical records indicating that European colonization induced widespread mortality among indigenous Americans.

Mitochondrial DNA variation in human metabolic rate and energy expenditure

The role of climate in driving selection of mtDNA as Homo sapiens migrated out of Africa into Eurasia remains controversial. We evaluated the role of mtDNA variation in resting metabolic rate (RMR) and total energy expenditure (TEE) among 294 older, community-dwelling African and European American adults from the Health, Aging and Body Composition Study. Common African haplogroups L0, L2 and L3 had significantly lower RMRs than European haplogroups H, JT and UK with haplogroup L1 RMR being intermediate to these groups. This study links mitochondrial haplogroups with ancestry-associated differences in metabolic rate and energy expenditure.

Large scale mitochondrial sequencing in Mexican Americans suggests a reappraisal of Native American origins

Background

The Asian origin of Native Americans is largely accepted. However uncertainties persist regarding the source population(s) within Asia, the divergence and arrival time(s) of the founder groups, the number of expansion events, and migration routes into the New World. mtDNA data, presented over the past two decades, have been used to suggest a single-migration model for which the Beringian land mass plays an important role.

Results

In our analysis of 568 mitochondrial genomes, the coalescent age estimates of shared roots between Native American and Siberian-Asian lineages, calculated using two different mutation rates, are A4 (27.5 ± 6.8 kya/22.7 ± 7.4 kya), C1 (21.4 ± 2.7 kya/16.4 ± 1.5 kya), C4 (21.0 ± 4.6 kya/20.0 ± 6.4 kya), and D4e1 (24.1 ± 9.0 kya/17.9 ± 10.0 kya). The coalescent age estimates of pan-American haplogroups calculated using the same two mutation rates (A2:19.5 ± 1.3 kya/16.1 ± 1.5 kya, B2:20.8 ± 2.0 kya/18.1 ± 2.4 kya, C1:21.4 ± 2.7 kya/16.4 ± 1.5 kya and D1:17.2 ± 2.0 kya/14.9 ± 2.2 kya) and estimates of population expansions within America (~21-16 kya), support the pre-Clovis occupation of the New World. The phylogeography of sublineages within American haplogroups A2, B2, D1 and the C1b, C1c andC1d subhaplogroups of C1 are complex and largely specific to geographical North, Central and South America. However some sub-branches (B2b, C1b, C1c, C1d and D1f) already existed in American founder haplogroups before expansion into the America.

Conclusions

Our results suggest that Native American founders diverged from their Siberian-Asian progenitors sometime during the last glacial maximum (LGM) and expanded into America soon after the LGM peak (~20-16 kya). The phylogeography of haplogroup C1 suggest that this American founder haplogroup differentiated in Siberia-Asia. The situation is less clear for haplogroup B2, however haplogroups A2 and D1 may have differentiated soon after the Native American founders divergence. A moderate population bottle neck in American founder populations just before the expansion most plausibly resulted in few founder types in America. The similar estimates of the diversity indices and Bayesian skyline analysis in North America, Central America and South America suggest almost simultaneous (~ 2.0 ky from South to North America) colonization of these geographical regions with rapid population expansion differentiating into more or less regional branches across the pan-American haplogroups.

Major population expansion of East Asians began before neolithic time: evidence of mtDNA genomes

It is a major question in archaeology and anthropology whether human populations started to grow primarily after the advent of agriculture, i.e., the Neolithic time, especially in East Asia, which was one of the centers of ancient agricultural civilization. To answer this question requires an accurate estimation of the time of lineage expansion as well as that of population expansion in a population sample without ascertainment bias. In this study, we analyzed all available mtDNA genomes of East Asians ascertained by random sampling, a total of 367 complete mtDNA sequences generated by the 1000 Genome Project, including 249 Chinese (CHB, CHD, and CHS) and 118 Japanese (JPT). We found that major mtDNA lineages underwent expansions, all of which, except for two JPT-specific lineages, including D4, D4b2b, D4a, D4j, D5a2a, A, N9a, F1a1'4, F2, B4, B4a, G2a1 and M7b1'2'4, occurred before 10 kya, i.e., before the Neolithic time (symbolized by Dadiwan Culture at 7.9 kya) in East Asia. Consistent to this observation, the further analysis showed that the population expansion in East Asia started at 13 kya and lasted until 4 kya. The results suggest that the population growth in East Asia constituted a need for the introduction of agriculture and might be one of the driving forces that led to the further development of agriculture.

Mitogroup: continent-specific clusters of mitochondrial OXPHOS complexes based on nuclear non-synonymous polymorphisms

OXPHOS polymorphisms are known to be population specific and to influence disease. Previous studies have focused on mtDNA polymorphisms. Based on a world sampling of 629 unrelated individuals, we have now studied the polymorphisms of the 80 genes encoding OXPHOS nuclear subunits. We have shown that (i) amino-acid replacement frequencies are significantly correlated with their pathogenicity probability, and (ii) populations can be distinguished based only on amino-acid replacements in nuclear encoded OXPHOS subunits. These results are congruent with the major mtDNA haplogroups, which suggests that OXPHOS complexes are different across the populations in both nuclear and in mitochondrial encoded subunits.

Positive Darwinian selection in the piston that powers proton pumps in complex I of the mitochondria of Pacific salmon

The mechanism of oxidative phosphorylation is well understood, but evolution of the proteins involved is not. We combined phylogenetic, genomic, and structural biology analyses to examine the evolution of twelve mitochondrial encoded proteins of closely related, yet phenotypically diverse, Pacific salmon. Two separate analyses identified the same seven positively selected sites in ND5. A strong signal was also detected at three sites of ND2. An energetic coupling analysis revealed several structures in the ND5 protein that may have co-evolved with the selected sites. These data implicate Complex I, specifically the piston arm of ND5 where it connects the proton pumps, as important in the evolution of Pacific salmon. Lastly, the lineage to Chinook experienced rapid evolution at the piston arm.

MtDNA SNP multiplexes for efficient inference of matrilineal genetic ancestry within Oceania

Human mitochondrial DNA (mtDNA) is a convenient marker for tracing matrilineal bio-geographic ancestry and is widely applied in forensic, genealogical and anthropological studies. In forensic applications, DNA-based ancestry inference can be useful for finding unknown suspects by concentrating police investigations in cases where autosomal STR profiling was unable to provide a match, or can help provide clues in missing person identification. Although multiplexed mtDNA single nucleotide polymorphism (SNP) assays to infer matrilineal ancestry at a (near) continental level are already available, such tools are lacking for the Oceania region. Here, we have developed a hierarchical system of three SNaPshot multiplexes for genotyping 26 SNPs defining all major mtDNA haplogroups for Oceania (including Australia, Near Oceania and Remote Oceania). With this system, it was possible to conclusively assign 74% of Oceanian individuals to their Oceanian matrilineal ancestry in an established literature database (after correcting for obvious external admixture). Furthermore, in a set of 161 genotyped individuals collected in Australia, Papua New Guinea and Fiji, 87.6% were conclusively assigned an Oceanian matrilineal origin. For the remaining 12.4% of the genotyped samples either a Eurasian origin was detected indicating likely European admixture (1.9%), the identified haplogroups are shared between Oceania and S/SE-Asia (5%), or the SNPs applied did not allow a geographic inference to be assigned (5.6%). Sub-regional assignment within Oceania was possible for 32.9% of the individuals genotyped: 49.5% of Australians were assigned an Australian origin and 13.7% of the Papua New Guineans were assigned a Near Oceanian origin, although none of the Fijians could be assigned a specific Remote Oceanian origin. The low assignment rates of Near and Remote Oceania are explained by recent migrations from Asia via Near Oceania into Remote Oceania. Combining the mtDNA multiplexes for Oceania introduced here with those we developed earlier for all other continental regions, global matrilineal bio-geographic ancestry assignment from DNA is now achievable in a highly efficient way that is also suitable for applications with limited material such as forensic case work.

The impact of modern migrations on present-day multi-ethnic Argentina as recorded on the mitochondrial DNA genome

BACKGROUND

The genetic background of Argentineans is a mosaic of different continental ancestries. From colonial to present times, the genetic contribution of Europeans and sub-Saharan Africans has superposed to or replaced the indigenous genetic 'stratum'. A sample of 384 individuals representing different Argentinean provinces was collected and genotyped for the first and the second mitochondrial DNA (mtDNA) hypervariable regions, and selectively genotyped for mtDNA SNPs. This data was analyzed together with additional 440 profiles from rural and urban populations plus 304 from Native American Argentineans, all available from the literature. A worldwide database was used for phylogeographic inferences, inter-population comparisons, and admixture analysis. Samples identified as belonging to hg (hg) H2a5 were sequenced for the entire mtDNA genome.

RESULTS

Phylogenetic and admixture analyses indicate that only half of the Native American component in urban Argentineans might be attributed to the legacy of extinct ancestral Argentineans and that the Spanish genetic contribution is slightly higher than the Italian one. Entire H2a5 genomes linked these Argentinean mtDNAs to the Basque Country and improved the phylogeny of this Basque autochthonous clade. The fingerprint of African slaves in urban Argentinean mtDNAs was low and it can be phylogeographically attributed predominantly to western African. The European component is significantly more prevalent in the Buenos Aires province, the main gate of entrance for Atlantic immigration to Argentina, while the Native American component is larger in North and South Argentina. AMOVA, Principal Component Analysis and hgs/haplotype patterns in Argentina revealed an important level of genetic sub-structure in the country.

CONCLUSIONS

Studies aimed to compare mtDNA frequency profiles from different Argentinean geographical regions (e.g., forensic and case-control studies) should take into account the important genetic heterogeneity of the country in order to prevent false positive claims of association in disease studies or inadequate evaluation of forensic evidence.

Inferring the population expansions in peopling of Japan

BACKGROUND

Extensive studies in different fields have been performed to reconstruct the prehistory of populations in the Japanese archipelago. Estimates the ancestral population dynamics based on Japanese molecular sequences can extend our understanding about the colonization of Japan and the ethnogenesis of modern Japanese.

METHODOLOGY/PRINCIPAL FINDINGS

We applied Bayesian skyline plot (BSP) with a dataset based on 952 Japanese mitochondrial DNA (mtDNA) genomes to depict the female effective population size (N(ef)) through time for the total Japanese and each of the major mtDNA haplogroups in Japanese. Our results revealed a rapid N(ef) growth since ∼5 thousand years ago had left ∼72% Japanese mtDNA lineages with a salient signature. The BSP for the major mtDNA haplogroups indicated some different demographic history.

CONCLUSIONS/SIGNIFICANCE

The results suggested that the rapid population expansion acted as a major force in shaping current maternal pool of Japanese. It supported a model for population dynamics in Japan in which the prehistoric population growth initiated in the Middle Jomon Period experienced a smooth and swift transition from Jomon to Yayoi, and then continued through the Yayoi Period. The confounding demographic backgrounds of different mtDNA haplogroups could also have some implications for some related studies in future.

Mutation rate switch inside Eurasian mitochondrial haplogroups: impact of selection and consequences for dating settlement in Europe

R-lineage mitochondrial DNA represents over 90% of the European population and is significantly present all around the planet (North Africa, Asia, Oceania, and America). This lineage played a major role in migration “out of Africa” and colonization in Europe. In order to determine an accurate dating of the R lineage and its sublineages, we analyzed 1173 individuals and complete mtDNA sequences from Mitomap. This analysis revealed a new coalescence age for R at 54.500 years, as well as several limitations of standard dating methods, likely to lead to false interpretations. These findings highlight the association of a striking under-accumulation of synonymous mutations, an over-accumulation of non-synonymous mutations, and the phenotypic effect on haplogroup J. Consequently, haplogroup J is apparently not a Neolithic group but an older haplogroup (Paleolithic) that was subjected to an underestimated selective force. These findings also indicated an under-accumulation of synonymous and non-synonymous mutations localized on coding and non-coding (HVS1) sequences for haplogroup R0, which contains the major haplogroups H and V. These new dates are likely to impact the present colonization model for Europe and confirm the late glacial resettlement scenario.

Mitochondrial complex 1 gene analysis in keratoconus

PURPOSE

Keratoconus is characterized by the thinning of corneal stroma, resulting in reduced vision. The exact etiology of keratoconus (KC) is still unknown. The involvement of oxidative stress (OS) in this disease has been reported. However, the exact mechanism of OS in keratoconus is still unknown. Thus we planned this study to screen mitochondrial complex I genes for sequence changes in keratoconus patients and controls, as mitochondrial complex I is the chief source of reactive oxygen species (ROS) production.

METHODS

A total of 20 keratoconus cases and 20 healthy controls without any ocular disorder were enrolled in this study. Mitochondrial complex I genes (ND1, 2, 3, 4, 4L, 5, and 6) were amplified in all patients and controls using 12 pairs of primers by PCR. After sequencing, DNA sequences were analyzed against the mitochondrial reference sequence NC_012920. Haplogroup frequency based Principle Component Analysis (PCA) was constructed to determine whether the gene pool of keratoconus patients is closer to major populations in India.

RESULTS

DNA sequencing revealed a total 84 nucleotide variations in patients and 29 in controls. Of 84 nucleotide changes, 18 variations were non-synonymous and two novel frame-shift mutations were detected in cases. Non-synonymous mtDNA sequence variations may account for increased ROS and decreased ATP production. This ultimately leads to OS; which is a known cause for variety of corneal abnormalities. Haplotype analysis showed that most of the patients were clustered under the haplogroups: T, C4a2a, R2'TJ, M21'Q1a, M12'G2a2a, M8'CZ and M7a2a, which are present as negligible frequency in normal Indian population, whereas only few patients were found to be a part of the other haplogroups like U7 (Indo-European), R2 and R31, whose origin is contentious.

CONCLUSIONS

Mt complex I sequence variations are the main cause of elevated ROS production which leads oxidative stress. This oxidative stress then starts a cascade of events which ultimately can lead to keratoconus. Prompt antioxidant therapy should be initiated in keratoconus patients to minimize ROS related damage.

Adaptive evolution signals in mitochondrial genes of Europeans

Since modern Europeans appear to be descendants of the Late Pleistocene European peoples who survived the last glacial period, it is quite reasonable to expect the presence of adaptive genetic variants that originated in the Ice Age in the modern gene pool of Europeans. To find such adaptive variants, mitochondrial genomes have been analyzed of the modern population from Eastern and Central Europe belonging to haplogroups U4, U5, and V, that diversified during the Late Pleistocene and Holocene periods. Analysis of distribution of nonsynonymous and synonymous substitutions, as well as results of search for radical amino acid changes that arose under the influence of adaptation (positive destabilizing selection) allowed us to detect signals of molecular adaptation in different mitochondrial genes and haplogroups of mtDNA. However, there were very few strong adaptive signals (z > 3.09, P < 0.001) that could be due to the loss of adaptive mtDNA haplotypes during the Holocene warming.

New insights into the Lake Chad Basin population structure revealed by high-throughput genotyping of mitochondrial DNA coding SNPs

BACKGROUND

Located in the Sudan belt, the Chad Basin forms a remarkable ecosystem, where several unique agricultural and pastoral techniques have been developed. Both from an archaeological and a genetic point of view, this region has been interpreted to be the center of a bidirectional corridor connecting West and East Africa, as well as a meeting point for populations coming from North Africa through the Saharan desert.

METHODOLOGY/PRINCIPAL FINDINGS

Samples from twelve ethnic groups from the Chad Basin (n = 542) have been high-throughput genotyped for 230 coding region mitochondrial DNA (mtDNA) Single Nucleotide Polymorphisms (mtSNPs) using Matrix-Assisted Laser Desorption/Ionization Time-Of-Flight (MALDI-TOF) mass spectrometry. This set of mtSNPs allowed for much better phylogenetic resolution than previous studies of this geographic region, enabling new insights into its population history. Notable haplogroup (hg) heterogeneity has been observed in the Chad Basin mirroring the different demographic histories of these ethnic groups. As estimated using a Bayesian framework, nomadic populations showed negative growth which was not always correlated to their estimated effective population sizes. Nomads also showed lower diversity values than sedentary groups.

CONCLUSIONS/SIGNIFICANCE

Compared to sedentary population, nomads showed signals of stronger genetic drift occurring in their ancestral populations. These populations, however, retained more haplotype diversity in their hypervariable segments I (HVS-I), but not their mtSNPs, suggesting a more ancestral ethnogenesis. Whereas the nomadic population showed a higher Mediterranean influence signaled mainly by sub-lineages of M1, R0, U6, and U5, the other populations showed a more consistent sub-Saharan pattern. Although lifestyle may have an influence on diversity patterns and hg composition, analysis of molecular variance has not identified these differences. The present study indicates that analysis of mtSNPs at high resolution could be a fast and extensive approach for screening variation in population studies where labor-intensive techniques such as entire genome sequencing remain unfeasible.

Purifying selection of mtDNA and its implications for understanding evolution and mitochondrial disease

Mutations of mitochondrial DNA (mtDNA) are frequent in humans and are implicated in many different types of pathology. The high substitution rate and the maternal, asexual mode of transmission of mtDNA make it more likely to accumulate deleterious mutations. Here, we discuss recent evidence that mtDNA transmission is subject to strong purifying selection in the mammalian female germ line, limiting the accumulation of such mutations. This process shapes mitochondrial sequence diversity and is therefore probably of fundamental importance for animal evolution and in human mitochondrial disease.

Mitochondrial genome sequences illuminate maternal lineages of conservation concern in a rare carnivore

Background

Science-based wildlife management relies on genetic information to infer population connectivity and identify conservation units. The most commonly used genetic marker for characterizing animal biodiversity and identifying maternal lineages is the mitochondrial genome. Mitochondrial genotyping figures prominently in conservation and management plans, with much of the attention focused on the non-coding displacement ("D") loop. We used massively parallel multiplexed sequencing to sequence complete mitochondrial genomes from 40 fishers, a threatened carnivore that possesses low mitogenomic diversity. This allowed us to test a key assumption of conservation genetics, specifically, that the D-loop accurately reflects genealogical relationships and variation of the larger mitochondrial genome.

Results

Overall mitogenomic divergence in fishers is exceedingly low, with 66 segregating sites and an average pairwise distance between genomes of 0.00088 across their aligned length (16,290 bp). Estimates of variation and genealogical relationships from the displacement (D) loop region (299 bp) are contradicted by the complete mitochondrial genome, as well as the protein coding fraction of the mitochondrial genome. The sources of this contradiction trace primarily to the near-absence of mutations marking the D-loop region of one of the most divergent lineages, and secondarily to independent (recurrent) mutations at two nucleotide position in the D-loop amplicon.

Conclusions

Our study has two important implications. First, inferred genealogical reconstructions based on the fisher D-loop region contradict inferences based on the entire mitogenome to the point that the populations of greatest conservation concern cannot be accurately resolved. Whole-genome analysis identifies Californian haplotypes from the northern-most populations as highly distinctive, with a significant excess of amino acid changes that may be indicative of molecular adaptation; D-loop sequences fail to identify this unique mitochondrial lineage. Second, the impact of recurrent mutation appears most acute in closely related haplotypes, due to the low level of evolutionary signal (unique mutations that mark lineages) relative to evolutionary noise (recurrent, shared mutation in unrelated haplotypes). For wildlife managers, this means that the populations of greatest conservation concern may be at the highest risk of being misidentified by D-loop haplotyping. This message is timely because it highlights the new opportunities for basing conservation decisions on more accurate genetic information.

Evolutionary analyses of entire genomes do not support the association of mtDNA mutations with Ras/MAPK pathway syndromes

Background

There are several known autosomal genes responsible for Ras/MAPK pathway syndromes, including Noonan syndrome (NS) and related disorders (such as LEOPARD, neurofibromatosis type 1), although mutations of these genes do not explain all cases. Due to the important role played by the mitochondrion in the energetic metabolism of cardiac muscle, it was recently proposed that variation in the mitochondrial DNA (mtDNA) genome could be a risk factor in the Noonan phenotype and in hypertrophic cardiomyopathy (HCM), which is a common clinical feature in Ras/MAPK pathway syndromes. In order to test these hypotheses, we sequenced entire mtDNA genomes in the largest series of patients suffering from Ras/MAPK pathway syndromes analyzed to date (n = 45), most of them classified as NS patients (n = 42).

Methods/Principal Findings

The results indicate that the observed mtDNA lineages were mostly of European ancestry, reproducing in a nutshell the expected haplogroup (hg) patterns of a typical Iberian dataset (including hgs H, T, J, and U). Three new branches of the mtDNA phylogeny (H1j1, U5b1e, and L2a5) are described for the first time, but none of these are likely to be related to NS or Ras/MAPK pathway syndromes when observed under an evolutionary perspective. Patterns of variation in tRNA and protein genes, as well as redundant, private and heteroplasmic variants, in the mtDNA genomes of patients were as expected when compared with the patterns inferred from a worldwide mtDNA phylogeny based on more than 8700 entire genomes. Moreover, most of the mtDNA variants found in patients had already been reported in healthy individuals and constitute common polymorphisms in human population groups.

Conclusions/Significance

As a whole, the observed mtDNA genome variation in the NS patients was difficult to reconcile with previous findings that indicated a pathogenic role of mtDNA variants in NS.

Comparing phylogeny and the predicted pathogenicity of protein variations reveals equal purifying selection across the global human mtDNA diversity

We used detailed phylogenetic trees for human mtDNA, combined with pathogenicity predictions for each amino acid change, to evaluate selection on mtDNA-encoded protein variants. Protein variants with high pathogenicity scores were significantly rarer in the older branches of the tree. Variants that have formed and survived multiple times in the human phylogenetics tree had significantly lower pathogenicity scores than those that only appear once in the tree. We compared the distribution of pathogenicity scores observed on the human phylogenetic tree to the distribution of all possible protein variations to define a measure of the effect of selection on these protein variations. The measured effect of selection increased exponentially with increasing pathogenicity score. We found no measurable difference in this measure of purifying selection in mtDNA across the global population, represented by the macrohaplogroups L, M, and N. We provide a list of all possible single amino acid variations for the human mtDNA-encoded proteins with their predicted pathogenicity scores and our measured selection effect as a tool for assessing novel protein variations that are often reported in patients with mitochondrial disease of unknown origin or for assessing somatic mutations acquired through aging or detected in tumors.

Mitochondrial-nuclear epistasis: implications for human aging and longevity

There is substantial evidence that mitochondria are involved in the aging process. Mitochondrial function requires the coordinated expression of hundreds of nuclear genes and a few dozen mitochondrial genes, many of which have been associated with either extended or shortened life span. Impaired mitochondrial function resulting from mtDNA and nuclear DNA variation is likely to contribute to an imbalance in cellular energy homeostasis, increased vulnerability to oxidative stress, and an increased rate of cellular senescence and aging. The complex genetic architecture of mitochondria suggests that there may be an equally complex set of gene interactions (epistases) involving genetic variation in the nuclear and mitochondrial genomes. Results from Drosophila suggest that the effects of mtDNA haplotypes on longevity vary among different nuclear allelic backgrounds, which could account for the inconsistent associations that have been observed between mitochondrial DNA (mtDNA) haplogroups and survival in humans. A diversity of pathways may influence the way mitochondria and nuclear-mitochondrial interactions modulate longevity, including: oxidative phosphorylation; mitochondrial uncoupling; antioxidant defenses; mitochondrial fission and fusion; and sirtuin regulation of mitochondrial genes. We hypothesize that aging and longevity, as complex traits having a significant genetic component, are likely to be controlled by nuclear gene variants interacting with both inherited and somatic mtDNA variability.

Mitochondrial DNA evidence supports northeast Indian origin of the aboriginal Andamanese in the Late Paleolithic

In view of the geographically closest location to Andaman archipelago, Myanmar was suggested to be the origin place of aboriginal Andamanese. However, for lacking any genetic information from this region, which has prevented to resolve the dispute on whether the aboriginal Andamanese were originated from mainland India or Myanmar. To solve this question and better understand the origin of the aboriginal Andamanese, we screened for haplogroups M31 (from which Andaman-specific lineage M31a1 branched off) and M32 among 846 mitochondrial DNAs (mtDNAs) sampled across Myanmar. As a result, two Myanmar individuals belonging to haplogroup M31 were identified, and completely sequencing the entire mtDNA genomes of both samples testified that the two M31 individuals observed in Myanmar were probably attributed to the recent gene flow from northeast India populations. Since no root lineages of haplogroup M31 or M32 were observed in Myanmar, it is unlikely that Myanmar may serve as the source place of the aboriginal Andamanese. To get further insight into the origin of this unique population, the detailed phylogenetic and phylogeographic analyses were performed by including additional 7 new entire mtDNA genomes and 113 M31 mtDNAs pinpointed from South Asian populations, and the results suggested that Andaman-specific M31a1 could in fact trace its origin to northeast India. Time estimation results further indicated that the Andaman archipelago was likely settled by modern humans from northeast India via the land-bridge which connected the Andaman archipelago and Myanmar around the Last Glacial Maximum (LGM), a scenario in well agreement with the evidence from linguistic and palaeoclimate studies.

Rapid, global demographic expansions after the origins of agriculture

The invention of agriculture is widely assumed to have driven recent human population growth. However, direct genetic evidence for population growth after independent agricultural origins has been elusive. We estimated population sizes through time from a set of globally distributed whole mitochondrial genomes, after separating lineages associated with agricultural populations from those associated with hunter-gatherers. The coalescent-based analysis revealed strong evidence for distinct demographic expansions in Europe, southeastern Asia, and sub-Saharan Africa within the past 10,000 y. Estimates of the timing of population growth based on genetic data correspond neatly to dates for the initial origins of agriculture derived from archaeological evidence. Comparisons of rates of population growth through time reveal that the invention of agriculture facilitated a fivefold increase in population growth relative to more ancient expansions of hunter-gatherers.

Rapid identification of mtDNA somatic mutations in gastric cancer tissues based on the mtDNA phylogeny

Mitochondrial DNA (mtDNA) somatic mutations have been identified in nearly all kinds of cancer during the past decade. Normally one need to determine the complete mtDNA sequences from both cancerous and normal tissues of the same patient to score the somatic mutation in cancer. In this study, we intended to explore a strategy to quickly identify somatic mutations in the entire mtDNA genome based on its phylogeny. The principal assumption for this strategy is that somatic mutations, as recently accumulated in cancerous tissue, have younger age and will be located in the terminal branches of mtDNA phylogenetic tree. In contrast, the haplogroup-specific variants, which appear as germline variants and have ancient age, will be located in the basal or intermediate-node branches of the tree, depending on their relative age. When the complete mtDNA sequence of the cancerous tissue is determined and is classified relative to the available mtDNA phylogeny, we only need to screen the variants that are located in the terminal branch in the paracancerous tissue or other normal tissue from the same patient to identify somatic mutations in cancer. We validated this strategy by using paired gastric cancer tissue and paracancerous tissue or blood from 10 Chinese patients (including one with gastric stromal tumor). A total of seven somatic mutations were identified in the cancerous tissues from four patients. Our result suggests that employing mtDNA phylogenetic knowledge facilitates rapid identification of mitochondrial genome somatic mutations in cancer.

Evaluating self-declared ancestry of U.S. Americans with autosomal, Y-chromosomal and mitochondrial DNA

The current U.S. population represents an amalgam of individuals originating mainly from four continental regions (Africa, Europe, Asia and America). To study the genetic ancestry and compare with self-declared ancestry we have analyzed paternally, maternally and bi-parentally inherited DNA markers sensitive for indicating continental genetic ancestry in all four major U.S. American groups. We found that self-declared U.S. Hispanics and U.S. African Americans tend to show variable degrees of continental genetic admixture among the three genetic systems, with evidence for a marked sex-biased admixture history. Moreover, for these two groups we observed significant regional variation across the country in genetic admixture. In contrast, self-declared U.S. European and U.S. Asian Americans were genetically more homogeneous at the continental ancestry level. Two autosomal ancestry-sensitive markers located in skin pigmentation candidate genes showed significant differences in self-declared U.S. African Americans or U.S. European Americans, relative to their assumed parental populations from Africa or Europe. This provides genetic support for the importance of skin color in the complex process of ancestry identification.

Genetic affinities between the Yami tribe people of Orchid Island and the Philippine Islanders of the Batanes archipelago

Background

Yami and Ivatan islanders are Austronesian speakers from Orchid Island and the Batanes archipelago that are located between Taiwan and the Philippines. The paternal genealogies of the Yami tribe from 1962 monograph of Wei and Liu were compared with our dataset of non-recombining Y (NRY) chromosomes from the corresponding families. Then mitochondrial DNA polymorphism was also analyzed to determine the matrilineal relationships between Yami, Ivatan, and other East Asian populations.

Results

The family relationships inferred from the NRY Phylogeny suggested a low number of paternal founders and agreed with the genealogy of Wei and Liu (P < 0.01). Except for one Y short tandem repeat lineage (Y-STR), seen in two unrelated Yami families, no other Y-STR lineages were shared between villages, whereas mtDNA haplotypes were indiscriminately distributed throughout Orchid Island.

The genetic affinity seen between Yami and Taiwanese aborigines or between Ivatan and the Philippine people was closer than that between Yami and Ivatan, suggesting that the Orchid islanders were colonized separately by their nearest neighbors and bred in isolation. However a northward gene flow to Orchid Island from the Philippines was suspected as Yami and Ivatan peoples both speak Western Malayo-Polynesian languages which are not spoken in Taiwan. Actually, only very little gene flow was observed between Yami and Ivatan or between Yami and the Philippines as indicated by the sharing of mtDNA haplogroup B4a1a4 and one O1a1* Y-STR lineage.

Conclusions

The NRY and mtDNA genetic information among Yami tribe peoples fitted well the patrilocal society model proposed by Wei and Liu. In this proposal, there were likely few genetic exchanges among Yami and the Philippine people. Trading activities may have contributed to the diffusion of Malayo-Polynesian languages among them.

Finally, artifacts dating 4,000 YBP, found on Orchid Island and indicating association with the Out of Taiwan hypothesis might be related to a pioneering stage of settlement, as most dating estimates inferred from DNA variation in our data set ranged between 100-3,000 YBP.