Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.     

Coevolution Predicts Direct Interactions between mtDNA-Encoded and nDNA-Encoded Subunits of Oxidative Phosphorylation Complex I

Despite years of research, the structure of the largest mammalian oxidative phosphorylation (OXPHOS) complex, NADH-ubiquinone oxidoreductase (complex I), and the interactions among its 45 subunits are not fully understood. Since complex I harbors subunits encoded by mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) genomes, with the former evolving ∼ 10 times faster than the latter, tight cytonuclear coevolution is expected and observed. Recently, we identified three nDNA-encoded complex I subunits that underwent accelerated amino acid replacement, suggesting their adjustment to the elevated mtDNA rate of change. Hence, they constitute excellent candidates for binding mtDNA-encoded subunits.Here, we further disentangle the network of physical cytonuclear interactions within complex I by analyzing subunits coevolution. Firstly, relying on the bioinformatic analysis of 10 protein complexes possessing solved structures, we show that signals of coevolution identified physically interacting subunits with nearly 90% accuracy, thus lending support to our approach. When applying this approach to cytonuclear interaction within complex I, we predict that the ‘rate-accelerated’ nDNA-encoded subunits of complex I, NDUFC2 and NDUFA1, likely interact with the mtDNA-encoded subunits ND5/ND4 and ND5/ND4/ND1, respectively. Furthermore, we predicted interactions among mtDNA-encoded complex I subunits. Using the yeast two-hybrid system, we experimentally confirmed the predicted interactions of human NDUFC2 with ND4, the interactions of human NDUFA1 with ND1 and ND4, and the lack of interaction of NDUFC2 with ND3 and NDUFA1, thus providing a proof of concept for our approach.Our study shows, for the first time, evidence for direct interactions between nDNA-encoded and mtDNA-encoded subunits of human OXPHOS complex I and paves the path towards deciphering subunit interactions within complexes lacking three-dimensional structures. Our subunit-interactions-predicting method, ComplexCorr, is available at http://webclu.bio.wzw.tum.de/complexcorr.     

Evolution of interacting proteins in the mitochondrial electron transport system in a marine copepod

The extensive interaction between mitochondrial-encoded and nuclear-encoded subunits of electron transport system (ETS) enzymes in mitochondria is expected to lead to intergenomic coadaptation. Whether this coadaptation results from adaptation to the environment or from fixation of deleterious mtDNA mutations followed by compensatory nuclear gene evolution is unknown. The intertidal copepod Tigriopus californicus shows extreme divergence in mtDNA sequence and provides an excellent model system for study of intergenomic coadaptation. Here, we examine genes encoding subunits of complex III of the ETS, including the mtDNA-encoded cytochrome b (CYTB), the nuclear-encoded rieske iron-sulfur protein (RISP), and cytochrome c(1) (CYC1). We compare levels of polymorphism within populations and divergence between populations in these genes to begin to untangle the selective forces that have shaped evolution in these genes. CYTB displays dramatic divergence between populations, but sequence analysis shows no evidence for positive selection driving this divergence. CYC1 and RISP have lower levels of sequence divergence between populations than CYTB, but, again, sequence analysis gives no evidence for positive selection acting on them. However, an examination of variation at cytochrome c (CYC), a nuclear-encoded protein that transfers electrons between complex III and complex IV provides evidence for selective divergence. Hence, it appears that rapid evolution in mitochondrial-encoded subunits is not always associated with rapid divergence in interacting subunits (CYC1 and RISP), but can be in some cases (CYC). Finally, a comparison of nuclear-encoded and mitochondrial-encoded genes from T. californicus suggests that substitution rates in the mitochondrial-encoded genes are dramatically increased relative to nuclear genes.     

Mitochondrial diseases

By convention, the term "mitochondrial diseases" refers to disorders of the mitochondrial respiratory chain, which is the only metabolic pathway in the cell that is under the dual control of the mitochondrial genome (mtDNA) and the nuclear genome (nDNA). Therefore, a genetic classification of the mitochondrial diseases distinguishes disorders due to mutations in mtDNA, which are governed by the relatively lax rules of mitochondrial genetics, and disorders due to mutations in nDNA, which are governed by the stricter rules of mendelian genetics. Mutations in mtDNA can be divided into those that impair mitochondrial protein synthesis in toto and those that affect any one of the 13 respiratory chain subunits encoded by mtDNA. Essential clinical features for each group of diseases are reviewed. Disorders due to mutations in nDNA are more abundant not only because most respiratory chain subunits are nucleus-encoded but also because correct assembly and functioning of the respiratory chain require numerous steps, all of which are under the control of nDNA. These steps (and related diseases) include: (i) synthesis of assembly proteins; (ii) intergenomic signaling; (iii) mitochondrial importation of nDNA-encoded proteins; (iv) synthesis of inner mitochondrial membrane phospholipids; (v) mitochondrial motility and fission.     

Pharmacological targeting of mitochondrial complex I deficiency: the cellular level and beyond

Complex I (CI) represents a major entry point of electrons in the mitochondrial electron transport chain (ETC). It consists of 45 different subunits, encoded by the mitochondrial (mtDNA) and nuclear DNA (nDNA). In humans, mutations in nDNA-encoded subunits cause severe neurodegenerative disorders like Leigh Syndrome with onset in early childhood. The pathophysiological mechanism of these disorders is still poorly understood. Here we summarize the current knowledge concerning the consequences of nDNA-encoded CI mutations in patient-derived cells, present mouse models for human CI deficiency, and discuss potential treatment strategies for CI deficiency.     

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.     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.     

Neoplastic transformation is associated with coordinate induction of nuclear and cytoplasmic oxidative phosphorylation genes

Neoplastic transformation was found to have a marked effect on the expression of nuclear DNA (nDNA)- and mitochondrial DNA (mtDNA)-encoded oxidative phosphorylation (OXPHOS) genes. Examining three pairs of human diploid fibroblasts and their SV 40-transformed counterparts revealed that mRNAs for the nuclear-encoded ATP synthase beta and the adenine nucleotide translocator (ANT) isoform 1 and 2 genes were markedly induced, whereas the mRNA for the ANT isoform 3 gene remained unchanged. The mRNA levels for the mtDNA-encoded 12 S rRNA, ND2, ATPase6+8, COIII, ND5+6, and Cytb genes were also increased, whereas the mtDNA number declined. Similar analysis of a cervical carcinoma (HeLa), fibrosarcoma (HT1080), and an Epstein-Barr virus (EBV)-transformed lymphoblastoid line (EBV-L) revealed that all three ANT isoforms were also expressed in these cells. Hence, changes in the expression of OXPHOS genes may be a common feature of transformed cells.     

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.     

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.     

Evaluating Purifying Selection in the Mitochondrial DNA of Various Mammalian Species

Mitochondrial DNA (mtDNA), the circular DNA molecule inside the mitochondria of all eukaryotic cells, has been shown to be under the effect of purifying selection in several species. Traditional testing of purifying selection has been based simply on ratios of nonsynonymous to synonymous mutations, without considering the relative age of each mutation, which can be determined by phylogenetic analysis of this non-recombining molecule. The incorporation of a mutation time-ordering from phylogeny and of predicted pathogenicity scores for nonsynonymous mutations allow a quantitative evaluation of the effects of purifying selection in human mtDNA. Here, by using this additional information, we show that purifying selection undoubtedly acts upon the mtDNA of other mammalian species/genera, namely Bos sp., Canis lupus, Mus musculus, Orcinus orca, Pan sp. and Sus scrofa. The effects of purifying selection were comparable in all species, leading to a significant major proportion of nonsynonymous variants with higher pathogenicity scores in the younger branches of the tree. We also derive recalibrated mutation rates for age estimates of ancestors of these various species and proposed a correction curve in order to take into account the effects of selection. Understanding this selection is fundamental to evolutionary studies and to the identification of deleterious mutations.     

Molecular Evolution of Cytochrome c Oxidase in High-Performance Fish (Teleostei: Scombroidei)

The 13 peptides encoded by vertebrate mitochondrial DNA (mtDNA) are essential subunits of oxidative phosphorylation (OXPHOS) enzymes. These genes normally experience purifying selection and also coevolve with nuclear-encoded subunits of OXPHOS complexes. However, the role of positive selection on mtDNA evolution is still unclear, as most examples of intergenomic coevolution appear to be the result of compensation by nuclear-encoded genes for mildly deleterious mtDNA mutations, and not simultaneous positive selection in both genomes. Organisms that have experienced strong selective pressures to increase aerobic capacity or adapt to changes in thermal environment may be better candidates in which to examine the impact of positively selected changes on mtDNA evolution. The tuna (suborder Scombroidei, family Scombridae) and billfish (suborder Scombroidei, families Xiphiidae and Istiophoridae) are highly aerobic fish with multiple specializations in muscle energetics, including a high mitochondrial content and regional endothermy. We examined the role of positively selected mtDNA substitutions in the production of these unique phenotypes. Focusing on a catalytic subunit of cytochrome c oxidase (COX II), we found that the rate ratio of nonsynonymous (d(N); amino acid changing)-to-synonymous (d(S); silent) substitutions was not increased in lineages leading to the tuna but was significantly increased in the lineage preceding the billfish. Furthermore, there are a number of individual positively selected sites that, when mapped onto the COX crystal structure, appear to interact with other COX subunits and may affect OXPHOS function and regulation in billfish.     

The age of nonsynonymous and synonymous mutations in animal mtDNA and implications for the mildly deleterious theory.

McDonald/Kreitman tests performed on animal mtDNA consistently reveal significant deviations from strict neutrality in the direction of an excess number of polymorphic nonsynonymous sites, which is consistent with purifying selection acting on nonsynonymous sites. We show that under models of recurrent neutral and deleterious mutations, the mean age of segregating neutral mutations is greater than the mean age of segregating selected mutations, even in the absence of recombination. We develop a test of the hypothesis that the mean age of segregating synonymous mutations equals the mean age of segregating nonsynonymous mutations in a sample of DNA sequences. The power of this age-of-mutation test and the power of the McDonald/Kreitman test are explored by computer simulations. We apply the new test to 25 previously published mitochondrial data sets and find weak evidence for selection against nonsynonymous 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.     

The role of climate in human mitochondrial DNA evolution: a reappraisal

Previous studies have proposed that selection has been involved in the differentiation of human mitochondrial DNA (mtDNA) and climate was the main driving force. This viewpoint, however, gets no support from the subsequent studies and remains controversial thus far. To clarify this issue, a total of 237 complete mtDNA sequences belonging to autochthonous lineages from South Asia, Oceania, and East Asia were collected to seek for the imprint of selection. Based on nonsynonymous (N) and synonymous (S) substitutions analysis, our results confirmed that purifying selection was the predominant force during the evolution of human mtDNA. However, no significant and extensive difference was detected among these three regions, which did not support the climate adaptation hypothesis but preferred random genetic drift to be the main factor in shaping the current landscape of human mtDNA, at least those from Asian and Oceanian regions.     

Mitochondrial Genetic Control of Assembly and Function of Complex I in Mammalian Cells

Sixteen years ago, we demonstrated, by immunological and biochemical approaches, that seven subunits of complex I are encoded in mitochondrial DNA (mtDNA) and synthesized on mitochondrial ribosomes in mammalian cells. More recently, we carried out a biochemical, molecular, and cellular analysis of a mutation in the gene for one of these subunits, ND4, that causes Leber's hereditary optic neuropathy (LHON). We demonstrated that, in cells carrying this mutation, the mtDNA-encoded subunits of complex I are assembled into a complex, but the rate of complex I-dependent respiration is decreased. Subsequently, we isolated several mutants affected in one or another of the mtDNA-encoded subunits of complex I by exposing established cell lines to high concentrations of rotenone. Our analyses of these mtDNA mutations affecting subunits of complex I have shown that at least two of these subunits, ND4 and ND6, are essential for the assembly of the enzyme. ND5 appears to be located at the periphery of the enzyme and, while it is not essential for assembly of the other mtDNA-encoded subunits into a complex, it is essential for complex I activity. In fact, the synthesis of the ND5 polypeptide is rate limiting for the activity of the enzyme.