Contrasting patterns of synonymous and nonsynonymous sequence evolution in asexual and sexual freshwater snail lineages

In asexual lineages, both synonymous and nonsynonymous sequence polymorphism may be reduced due to severe founder effects when asexual lineages originate. However, mildly deleterious (nonsynonymous) mutations may accumulate after asexual lineages are formed, because the efficiency of purifying selection is reduced even in the nonrecombining mitochondrial genome. Here we examine patterns of synonymous and nonsynonymous mitochondrial sequence polymorphism in asexual and sexual lineages of the freshwater snail Campeloma. Using clade-specific estimates, we found that synonymous sequence polymorphism was significantly reduced by 75% in asexuals relative to sexuals, whereas nonsynonymous sequence polymorphism did not differ significantly between sexuals and asexuals. Two asexual clades had high negative values for Tajima's D statistic. Coalescent simulations confirmed that various bottleneck scenarios can account for this result. We also used branch-specific estimates of the ratio of amino acid to silent substitutions, K(a)/K(s). Our study revealed that K(a)/K(s) ratios are six times higher in terminal branches of independent asexual lineages compared to sexuals. Coalescent-based reconstruction of gene networks for all sexual and asexual clades indicated that nonsynonymous mutations occurred at a higher frequency in recently derived asexual haplotypes. These findings suggest that patterns of synonymous and nonsynonymous nucleotide polymorphism in asexual snail lineages may be shaped by both severe founder effect and relaxed purifying selection.     

Mitochondrial D-Loop Variation in Persian Multiple Sclerosis Patients: K and A Haplogroups as a Risk Factor!!

Multiple Sclerosis (MS) is a multifocal demyelinating central nervous system disorder in which interplay between genes and the environment are supposed to be involved. Mitochondrial DNA (mtDNA) has the only non-coding regions at the displacement loop (D-loop) region that contains two hypervariable segments (HVS-I and HVS-II) with high polymorphism. mtDNA has already been fully sequenced and many subsequent publications have showed polymorphic sites, haplogroups and haplotypes. Haplogroups could have important implications to understand association between mutability of the mitochondrial genome and disease. To assess relationship between mtDNA haplogroups and MS, we have sequenced the mtDNA HVS-I in 54 MS patients and 100 control subjects. We have found that haplogroups A and K are significantly more abundant in MS patients (P=0.042 for haplogroup A and P=0.0005 for haplogroup K). Thus, these two haplogroups might act synergistically to increase the penetrance of MS disease.     

Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups

The evolution of the human mitochondrial genome is characterized by the emergence of ethnically distinct lineages or haplogroups. Nine European, seven Asian (including Native American), and three African mitochondrial DNA (mtDNA) haplogroups have been identified previously on the basis of the presence or absence of a relatively small number of restriction-enzyme recognition sites or on the basis of nucleotide sequences of the D-loop region. We have used reduced-median-network approaches to analyze 560 complete European, Asian, and African mtDNA coding-region sequences from unrelated individuals to develop a more complete understanding of sequence diversity both within and between haplogroups. A total of 497 haplogroup-associated polymorphisms were identified, 323 (65%) of which were associated with one haplogroup and 174 (35%) of which were associated with two or more haplogroups. Approximately one-half of these polymorphisms are reported for the first time here. Our results confirm and substantially extend the phylogenetic relationships among mitochondrial genomes described elsewhere from the major human ethnic groups. Another important result is that there were numerous instances both of parallel mutations at the same site and of reversion (i.e., homoplasy). It is likely that homoplasy in the coding region will confound evolutionary analysis of small sequence sets. By a linkage-disequilibrium approach, additional evidence for the absence of human mtDNA recombination is presented here.     

The role of mtDNA background in disease expression: a new primary LHON mutation associated with Western Eurasian haplogroup J

Leber's hereditary optic neuropathy (LHON) is a maternally transmitted form of blindness caused by mitochondrial DNA (mtDNA) mutations. Approximately 90% of LHON cases are caused by 3460A, 11778A, or 14484C mtDNA mutations. These are designated "primary" mutations because they impart a high risk for LHON expression. Although the 11778A and 14484C mutations unequivocally predispose carriers to LHON, they are preferentially associated with mtDNA haplogroup J, one of nine Western Eurasian mtDNA lineages, suggesting a synergistic and deleterious interaction between these LHON mutations and haplogroup J polymorphism(s). We report here the characterization of a new primary LHON mutation in the mtDNA ND4L gene at nucleotide pair 10663. The homoplasmic 10663C mutation has been found in three independent LHON patients who lack a known primary mutation and all of which belong to haplogroup J. This mutation has not been found in a large number of haplotype-matched or non-haplogroup-J control mtDNAs. Phylogenetic analysis with primarily complete mtDNA sequence data demonstrates that the 10663C mutation has arisen at least three independent times in haplogroup J, indicating that it is not a rare lineage-specific polymorphism. Analysis of complex I function in patient lymphoblasts and transmitochondrial cybrids has revealed a partial complex I defect similar in magnitude to the 14484C mutation. Thus, the 10663C mutation appears to be a new primary LHON mutation that is pathogenic when co-occurring with haplogroup J. These results strongly support a role for haplogroup J in the expression of certain LHON mutations.     

Strong purifying selection in transmission of mammalian mitochondrial DNA

There is an intense debate concerning whether selection or demographics has been most important in shaping the sequence variation observed in modern human mitochondrial DNA (mtDNA). Purifying selection is thought to be important in shaping mtDNA sequence evolution, but the strength of this selection has been debated, mainly due to the threshold effect of pathogenic mtDNA mutations and an observed excess of new mtDNA mutations in human population data. We experimentally addressed this issue by studying the maternal transmission of random mtDNA mutations in mtDNA mutator mice expressing a proofreading-deficient mitochondrial DNA polymerase. We report a rapid and strong elimination of nonsynonymous changes in protein-coding genes; the hallmark of purifying selection. There are striking similarities between the mutational patterns in our experimental mouse system and human mtDNA polymorphisms. These data show strong purifying selection against mutations within mtDNA protein-coding genes. To our knowledge, our study presents the first direct experimental observations of the fate of random mtDNA mutations in the mammalian germ line and demonstrates the importance of purifying selection in shaping mitochondrial sequence diversity.     

Dynamics of haplogroup frequencies and survival rates in a contact zone of two mtDNA lineages of the lizard Lacerta vivipara

The analysis of contact zones between lineages that were previously isolated in allopatry can lead to important insights on evolutionary processes such as selection and adaptation. In this paper we conducted a comparative demographic study of two mitochondrial DNA (mtDNA) lineages of the lizard Lacerta vivipara in the western Pyrénées to provide detail on the dynamics of their contact zone. By surveying haplogroup frequency across the contact area, we revealed the existence of a stable and very narrow contact zone between two parapatric lineages, which we infer to demonstrate a role for selection in the maintenance of this contact zone. We suggest these two lineages evolved in allopatry after retreating to different refugia during the Pleistocene glaciations, and subsequently came into secondary contact after the last glacial maximum. Although haplogroup frequencies were stable over time, we found significant age and environment (temperature) dependent survival differences between mtDNA haplogroups in one contact population sampled yearly from 2002 to 2009. Therefore, temperature‐induced demographic differences between the two mtDNA lineages may be responsible for the stability of this narrow contact zone. This is one of the first demographic studies conducted under natural conditions indicating the possibility of selection on mtDNA.     

Nonneutral evolution and differential mutation rate of gender-associated mitochondrial DNA lineages in the marine mussel Mytilus.

Mussels have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally, and the F type is transmitted maternally. To test hypotheses of the molecular evolution of both mtDNA genomes, 50 nucleotide sequences were obtained for 396 bp of the COIII gene of European populations of Mytilus edulis and the Atlantic and Mediterranean forms of M. galloprovincialis. Analysis based on the proportion of synonymous and nonsynonymous substitutions indicate that mtDNA is evolving in a non-neutral and complex fashion. Previous studies on American mussels demonstrated that the F genome experiences a higher purifying selection and that the M genome evolves faster. Here we show that these patterns also hold in European populations. However, in contrast to American populations, where an excess of replacement substitution between F and M lineages has been reported, a significant excess of replacement polymorphism within mtDNA lineages is observed in European populations of M. galloprovincialis. European populations also show an excess of replacement polymorphism within the F but not within the M genome with respect to American M. trossulus, as well as a consistent pattern of excess of rare variants in both F and M genomes. These results are consistent with a nearly neutral model of molecular evolution and a recent relaxation of selective constraints on European mtDNA. Levels of diversity are significantly higher for the M than F genome, and the M genome also accumulates synonymous and nonsynonymous substitutions at a higher rate, in contrast with earlier reports where no difference for the synonymous rate was observed. It is suggested that a subtle balance between relaxed selection and a higher mutation rate explains the faster evolutionary rate of the M lineage.     

Do Haplogroups H and U Act to Increase the Penetrance of Alzheimer’s Disease?

1. Alzheimer’s disease (AD) is the most common form of dementia in the elderly in which interplay between genes and the environment is supposed to be involved. Mitochondrial DNA (mtDNA) has the only noncoding regions at the displacement loop (D-loop) region that contains two hypervariable segments (HVS-I and HVS-II) with high polymorphism. mtDNA has already been fully sequenced and many subsequent publications have shown polymorphic sites, haplogroups, and haplotypes. Haplogroups could have important implications to understand the association between mutability of the mitochondrial genome and the disease. 2. To assess the relationship between mtDNA haplogroup and AD, we sequenced the mtDNA HVS-I in 30 AD patients and 100 control subjects. We could find that haplogroups H and U are significantly more abundant in AD patients (P = 0.016 for haplogroup H and P = 0.0003 for haplogroup U), Thus, these two haplogroups might act synergistically to increase the penetrance of AD disease.     

Gene Expression Pattern in Transmitochondrial Cytoplasmic Hybrid Cells Harboring Type 2 Diabetes-Associated Mitochondrial DNA Haplogroups

Decreased mitochondrial function plays a pivotal role in the pathogenesis of type 2 diabetes mellitus (T2DM). Recently, it was reported that mitochondrial DNA (mtDNA) haplogroups confer genetic susceptibility to T2DM in Koreans and Japanese. Particularly, mtDNA haplogroup N9a is associated with a decreased risk of T2DM, whereas haplogroups D5 and F are associated with an increased risk. To examine functional consequences of these haplogroups without being confounded by the heterogeneous nuclear genomic backgrounds of different subjects, we constructed transmitochondrial cytoplasmic hybrid (cybrid) cells harboring each of the three haplogroups (N9a, D5, and F) in a background of a shared nuclear genome. We compared the functional consequences of the three haplogroups using cell-based assays and gene expression microarrays. Cell-based assays did not detect differences in mitochondrial functions among the haplogroups in terms of ATP generation, reactive oxygen species production, mitochondrial membrane potential, and cellular dehydrogenase activity. However, differential expression and clustering analyses of microarray data revealed that the three haplogroups exhibit a distinctive nuclear gene expression pattern that correlates with their susceptibility to T2DM. Pathway analysis of microarray data identified several differentially regulated metabolic pathways. Notably, compared to the T2DM-resistant haplogroup N9a, the T2DM-susceptible haplogroup F showed down-regulation of oxidative phosphorylation and up-regulation of glycolysis. These results suggest that variations in mtDNA can affect the expression of nuclear genes regulating mitochondrial functions or cellular energetics. Given that impaired mitochondrial function caused by T2DM-associated mtDNA haplogroups is compensated by the nuclear genome, we speculate that defective nuclear compensation, under certain circumstances, might lead to the development of T2DM.     

High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria

When a population size is reduced, genetic drift may fix slightly deleterious mutations, and an increase in nonsynonymous substitution is expected. It has been suggested that past aridity has seriously affected and decreased the populations of cichlid fishes in Lake Victoria, while geographical studies have shown that the water levels in Lake Tanganyika and Lake Malawi have remained fairly constant. The comparably stable environments in the latter two lakes might have kept the populations of cichlid fishes large enough to remove slightly deleterious mutations. The difference in the stability of cichlid fish population sizes between Lake Victoria and the Lakes Tanganyika and Malawi is expected to have caused differences in the nonsynonymous/synonymous ratio, ω (= d_N/d_S), of the evolutionary rate. Here, we estimated ω and compared it between the cichlids of the three lakes for 13 mitochondrial protein-coding genes using maximum likelihood methods. We found that the lineages of the cichlids in Lake Victoria had a significantly higher ω for several mitochondrial loci. Moreover, positive selection was indicated for several codons in the mtDNA of the Lake Victoria cichlid lineage. Our results indicate that both adaptive and slightly deleterious molecular evolution has taken place in the Lake Victoria cichlids' mtDNA genes, whose nonsynonymous sites are generally conserved.     

High divergence among Drosophila simulans mitochondrial haplogroups arose in midst of long term purifying selection

We characterize the type of selection acting within and among mitochondrial lineages in five closely related Drosophila species. We focus on D. simulans, where three genetically distinct mitochondrial haplogroups show high interhaplogroup divergence and low intrahaplogroup polymorphism. Using maximum likelihood models we find that the branches leading to these three distinct mitochondrial groups show a significantly reduced rate of nonsynonymous relative to synonymous substitution. This interhaplogroup rate is significantly reduced compared to the intrahaplogroup rate, and closely resembles the rate observed between distinct species. The data suggest that slightly deleterious mutations segregating within D. simulans haplogroups are removed by selection prior to their fixation among haplogroups. We explore several hypotheses to explain how lineages within a single species can be compatible with this model of slightly deleterious mutation. The most likely hypothesis is that D. simulans haplogroups have persisted in isolation, perhaps due to association with the bacterial symbiont Wolbachia and/or demographic history, introducing a bias against the fixation of slightly deleterious mutations.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

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.     

Mitochondrial haplogroup B is negatively associated with elite Korean endurance athlete status

Mitochondrial DNA (mtDNA) mutations could contribute to aerobic performance, since they provide information to generate aerobic ATP energy by oxidative phosphorylation (OXPHOS) in their respiratory chain. Owing to haploid and absence of recombination, specific mutations in the mtDNA genome associated with human exercise tolerance or intolerance arise and remain in particular genetic lineages referred to as haplogroups. Recent studies have suggested that certain mtDNA lineages were associated with individual differences in trainability and physical capacity, but their origins and replicate test remain to be elucidated. We have therefore analyzed mtDNA haplogroup B variants to assess the possible role of the lineages to differences in elite athletic performance in a population-based case-control study in Cheonan, Korea. We demonstrated a significant negative association between mtDNA haplogroup B and the status of elite endurance athlete [n=378, odds ratio = 0.37 (95% CI: 0.14?C0.97, p = 0.016)]. Thus, our results imply that specific mtDNA lineages may provide a significant effect on elite Korean endurance status, although larger sample sizes functional studies are necessary to further substantiate these findings.     

Mitochondrial Haplogroups Modify the Risk of Developing Hypertrophic Cardiomyopathy in a Danish Population

Hypertrophic cardiomyopathy (HCM) is a genetic disorder caused by mutations in genes coding for proteins involved in sarcomere function. The disease is associated with mitochondrial dysfunction. Evolutionarily developed variation in mitochondrial DNA (mtDNA), defining mtDNA haplogroups and haplogroup clusters, is associated with functional differences in mitochondrial function and susceptibility to various diseases, including ischemic cardiomyopathy. We hypothesized that mtDNA haplogroups, in particular H, J and K, might modify disease susceptibility to HCM. Mitochondrial DNA, isolated from blood, was sequenced and haplogroups identified in 91 probands with HCM. The association with HCM was ascertained using two Danish control populations. Haplogroup H was more prevalent in HCM patients, 60% versus 46% (p = 0.006) and 41% (p = 0.003), in the two control populations. Haplogroup J was less prevalent, 3% vs. 12.4% (p = 0.017) and 9.1%, (p = 0.06). Likewise, the UK haplogroup cluster was less prevalent in HCM, 11% vs. 22.1% (p = 0.02) and 22.8% (p = 0.04). These results indicate that haplogroup H constitutes a susceptibility factor and that haplogroup J and haplogroup cluster UK are protective factors in the development of HCM. Thus, constitutive differences in mitochondrial function may influence the occurrence and clinical presentation of HCM. This could explain some of the phenotypic variability in HCM. The fact that haplogroup H and J are also modifying factors in ischemic cardiomyopathy suggests that mtDNA haplotypes may be of significance in determining whether a physiological hypertrophy develops into myopathy. mtDNA haplotypes may have the potential of becoming significant biomarkers in cardiomyopathy.     

Mitochondrial genome of the colorless green alga Polytomella parva: two linear DNA molecules with homologous inverted repeat Termini

Most of the well-characterized mitochondrial genomes from diverse green algal lineages are circular mapping DNA molecules; however, Chlamydomonas reinhardtii has a linear 15.8 kb unit mitochondrial genome with 580 or 581 bp inverted repeat ends. In mitochondrial-enriched fractions prepared from Polytomella parva (=P. agilis), a colorless, naturally wall-less relative of C. reinhardtii, we have detected two linear mitochondrial DNA (mtDNA) components with sizes of 13.5 and 3.5 kb. Sequences spanning 97% and 86% of the 13.5- and 3.5-kb mtDNAs, respectively, reveal that these molecules contain long, at least 1.3 kb, homologous inverted repeat sequences at their termini. The 3.5-kb mtDNA has only one coding region (nad6), the functionality of which is supported by both the relative rate at which it has accumulated nonsynonymous nucleotide substitutions and its absence from the 13.5-kb mtDNA which encodes nine genes (i.e., large and small subunit rRNA [LSU and SSU rRNA] genes, one tRNA gene, and six protein-coding genes). On the basis of DNA sequence data, we propose that a variant start codon, GTG, is utilized by the P. parva 13.5-kb mtDNA-encoded gene, nad5. Using the relative rate test with Chlamydomonas moewusii (=C. eugametos) as the outgroup, we conclude that the nonsynonymous nucleotide substitution rate in the mitochondrial protein-coding genes of P. parva is on an average about 3.3 times that of the C. reinhardtii counterparts.     

Mitochondrial DNA polymorphisms as risk factors for Parkinson’s disease and Parkinson’s disease dementia

The activity of complex I of the mitochondrial respiratory chain has been found to be decreased in patients with Parkinsons disease (PD), but no mutations have been identified in genes encoding complex I subunits. Recent studies have suggested that polymorphisms in mitochondrial DNA (mtDNA)-encoded complex I genes (MTND) modify susceptibility to PD. We hypothesize that the risk of PD is conveyed by the total number of nonsynonymous substitutions in the MTND genes in various mtDNA lineages rather than by single mutations. To test this possibility, we determined the number of nonsynonymous substitutions of the seven MTND genes from 183 Finns. The differences in the total number of nonsynonymous substitutions and the nonsynonymous to synonymous substitution rate ratio (K_a/K_s) of MTND genes between the European mtDNA haplogroup clusters (HV, JT, KU, IWX) were analysed by using a statistical approach. Patients with PD (n=238) underwent clinical examination together with mtDNA haplogroup analysis and the clinical features between patient groups defined by the number of nonsynonymous substitutions were compared. Our analysis revealed that the haplogroup clusters HV and KU had a lower average number of amino acid replacements and a lower K_a/K_s ratio in the MTND genes than clusters JT and IWX. Supercluster JTIWX with the highest number of amino acid replacements was more frequent among PD patients and even more frequent among patients with PD who developed dementia. Our results suggest that a relative excess of nonsynonymous mutations in MTND genes in supercluster JTWIX is associated with an increased risk of PD and the disease progression to dementia.     

Analysis of distribution of "mongoloid" haplogroups of mitochondrial DNA among indigenous population of the Tuva Republic

Variation of Mongoloid-specific restriction sites of mitochondrial genome was analyzed in three territorial groups of Tuvinians. Distribution of mitochondrial DNA haplogroups A, B, C, and D on the territory of the Tuva Republic was estimated. The populations studied did not display distinct differentiation in respect to the mtDNA polymorphism. The specific feature of Tuvinian mitochondrial gene pool was the prevalence of only one haplogroup C (over 40%), mainly represented by two mitotypes. The high frequency of this haplogroup makes Tuvinians similar to more northern Siberian populations. On the other hand, the presence of haplogroup B indicates that Tuvinians have affinity to ethnic groups of Central Asia.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Whole-mtDNA genome sequence analysis of ancient African lineages

Studies of human mitochondrial (mt) DNA genomes demonstrate that the root of the human phylogenetic tree occurs in Africa. Although 2 mtDNA lineages with an African origin (haplogroups M and N) were the progenitors of all non-African haplogroups, macrohaplogroup L (including haplogroups L0-L6) is limited to sub-Saharan Africa. Several L haplogroup lineages occur most frequently in eastern Africa (e.g., L0a, L0f, L5, and L3g), but some are specific to certain ethnic groups, such as haplogroup lineages L0d and L0k that previously have been found nearly exclusively among southern African "click" speakers. Few studies have included multiple mtDNA genome samples belonging to haplogroups that occur in eastern and southern Africa but are rare or absent elsewhere. This lack of sampling in eastern Africa makes it difficult to infer relationships among mtDNA haplogroups or to examine events that occurred early in human history. We sequenced 62 complete mtDNA genomes of ethnically diverse Tanzanians, southern African Khoisan speakers, and Bakola Pygmies and compared them with a global pool of 226 mtDNA genomes. From these, we infer phylogenetic relationships amongst mtDNA haplogroups and estimate the time to most recent common ancestor (TMRCA) for haplogroup lineages. These data suggest that Tanzanians have high genetic diversity and possess ancient mtDNA haplogroups, some of which are either rare (L0d and L5) or absent (L0f) in other regions of Africa. We propose that a large and diverse human population has persisted in eastern Africa and that eastern Africa may have been an ancient source of dispersion of modern humans both within and outside of Africa.