Evolution of nuclear- and mitochondrial-encoded subunit interaction in cytochrome c oxidase

Mitochondrial citrate transporters Ctp1-Yhm2 and respiratory chain: A coordinated functional connection in Saccharomyces cerevisiae metabolism

The mitochondrial inner membrane contains some hydrophobic proteins that mediate the exchange of metabolites between the mitochondrial matrix and the cytosol. Ctp1 and Yhm2 are two carrier proteins in the yeast Saccharomyces cerevisiae responsible for the transport of citrate, a tricarboxylate involved in several metabolic pathways. Since these proteins also contribute to respiratory metabolism, in this study we investigated for the first time whether changes in citrate transport can affect the structural organization and functional properties of respiratory complexes. Through experiments in yeast mutant cells in which the gene encoding Ctp1 or Yhm2 was deleted, we found that in the absence of either mitochondrial citrate transporter, mitochondrial respiration was impaired. Structural analysis of the respiratory complexes III and IV revealed different expression levels of the catalytic and supernumerary subunits in the Δctp1 and Δyhm2 strains. In addition, Δyhm2 mitochondria appeared to be more sensitive than Δctp1 to the oxidative damage. Our results provide the first evidence for a coordinated modulation of mitochondrial citrate transport and respiratory chain activity in S. cerevisiae metabolism.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

Mitochondrial effects on fertility and longevity in Tigriopus californicus contradict predictions of the mother's curse hypothesis

Strict maternal inheritance of mitochondria favours the evolutionary accumulation of sex-biased fitness effects, as mitochondrial evolution occurs exclusively in female lineages. The 'mother's curse' hypothesis proposes that male-harming mutations should accumulate in mitochondrial genomes when they have neutral or beneficial effects on female fitness. Rigorous empirical tests have largely focused on Drosophila, where support for the predictions of mother's curse has been mixed. We investigated the impact of mother's curse mutations in Tigriopus californicus, a minute crustacean. Using non-recombinant backcrosses, we introgressed four divergent mitochondrial haplotypes into two nuclear backgrounds and recorded measures of fertility and longevity. We found that the phenotypic effects of mitochondrial mutations were context dependent, being influenced by the nuclear background in which they were expressed, as well as the sex of the individual and rearing temperature. Mitochondrial haplotype effects were greater for fertility than longevity, and temperature effects were greater for longevity. However, in opposition to mother's curse expectations, females had higher mitochondrial genetic variance than males for fertility and longevity, little evidence of sexual antagonism favouring females was found, and the impacts of mitonuclear mismatch harmed females but not males. Together, this indicates that selection on mitochondrial variation has not resulted in the accumulation of male mutation load in Tigriopus californicus.

Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Animals living in extremely high elevations have to adapt to low temperatures and low oxygen availability (hypoxia), but the underlying genetic mechanisms associated with these adaptations are still unclear. The mitochondrial respiratory chain can provide >95% of the ATP in animal cells, and its efficiency is influenced by temperature and oxygen availability. Therefore, the respiratory chain complexes (RCCs) could be important molecular targets for positive selection associated with respiratory adaptation in high-altitude environments. Here, we investigated positive selection in 5 RCCs and their assembly factors by analyzing sequences of 106 genes obtained through RNA-seq of all 15 Chinese Phrynocephalus lizard species, which are distributed from lowlands to the Tibetan plateau (average elevation >4,500 m). Our results indicate that evidence of positive selection on RCC genes is not significantly different from assembly factors, and we found no difference in selective pressures among the 5 complexes. We specifically looked for positive selection in lineages where changes in habitat elevation happened. The group of lineages evolving from low to high altitude show stronger signals of positive selection than lineages evolving from high to low elevations. Lineages evolving from low to high elevation also have more shared codons under positive selection, though the changes are not equivalent at the amino acid level. This study advances our understanding of the genetic basis of animal respiratory metabolism evolution in extreme high environments and provides candidate genes for further confirmation with functional analyses.

An evolutionary, or "Mitocentric" perspective on cellular function and disease

The incidence of common, metabolic diseases (e.g. obesity, cardiovascular disease, diabetes) with complex genetic etiology has been steadily increasing nationally and globally. While identification of a genetic model that explains susceptibility and risk for these diseases has been pursued over several decades, no clear paradigm has yet been found to disentangle the genetic basis of polygenic/complex disease development. Since the evolution of the eukaryotic cell involved a symbiotic interaction between the antecedents of the mitochondrion and nucleus (which itself is a genetic hybrid), we suggest that this history provides a rational basis for investigating whether genetic interaction and co-evolution of these genomes still exists. We propose that both mitochondrial and Mendelian, or "mito-Mendelian" genetics play a significant role in cell function, and thus disease risk. This paradigm contemplates the natural variation and co-evolution of both mitochondrial and nuclear DNA backgrounds on multiple mitochondrial functions that are discussed herein, including energy production, cell signaling and immune response, which collectively can influence disease development. At the nexus of these processes is the economy of mitochondrial metabolism, programmed by both mitochondrial and nuclear genomes.

Recurrent gene duplication in the angiosperm tribe Delphinieae (Ranunculaceae) inferred from intracellular gene transfer events and heteroplasmic mutations in the plastid matK gene

The study of intracellular gene transfer may allow for the detection of interesting evolutionary processes such as ancient polyploidization. We compared 24 plastid genomes (plastomes) from tribe Delphinieae, one from tribe Nigelleae and one from tribe Ranunculeae, including five newly sequenced genomes. The functional transfers of the plastids rpl32 and rps16 to the nucleus in tribe Delphinieae were identified. Unexpectedly, we discovered multiple divergent copies of the nuclear-encoded plastid rpl32 in the genus Aconitum. Phylogenetic and synonymous substitution rate analyses revealed that the nuclear-encoded plastid rpl32 underwent two major duplication events. These ancient gene duplication events probably occurred via multiple polyploidization events in Aconitum between 11.9 and 24.7 Mya. Furthermore, our sequence rate analysis indicated that the eight plastid-encoded rpl subunits in Aconitum had a significantly accelerated evolutionary rate compared to those in other genera, suggesting that highly divergent paralogs targeted to the plastid may contribute to an elevated rate of evolution in plastid rpl genes. In addition, heteroplasmy of the plastid matK from two Aconitum species suggested the existence of potentially functional plastid maturases in its plastome. Our results provide insight into the evolutionary history of the tribe Delphinieae.

Characterization of the Complete Mitochondrial Genome of Harpalus sinicus and Its Implications for Phylogenetic Analyses

In this study, we report the complete mitochondrial genome of Harpalus sinicus (occasionally named as the Chinese ground beetle) which is the first mitochondrial genome for Harpalus. The mitogenome is 16,521 bp in length, comprising 37 genes, and a control region. The A + T content of the mitogenome is as high as 80.6%. A mitochondrial origins of light-strand replication (OL)-like region is found firstly in the insect mitogenome, which can form a stem-loop hairpin structure. Thirteen protein-coding genes (PCGs) share high homology, and all of them are under purifying selection. All tRNA genes (tRNAs) can be folded into the classic cloverleaf secondary structures except tRNA-Ser (GCU), which lacks a dihydrouridine (DHU) stem. The secondary structure of two ribosomal RNA genes (rRNAs) is predicted based on previous insect models. Twelve types of tandem repeats and two stem-loop structures are detected in the control region, and two stem-loop structures may be involved in the initiation of replication and transcription. Additionally, phylogenetic analyses based on mitogenomes suggest that Harpalus is an independent lineage in Carabidae, and is closely related to four genera (Abax, Amara, Stomis, and Pterostichus). In general, this study provides meaningful genetic information for Harpalus sinicus and new insights into the phylogenetic relationships within the Carabidae.

Mitochondrial genomes of the stoneflies Mesonemourametafiligera and Mesonemouratritaenia (Plecoptera, Nemouridae), with a phylogenetic analysis of Nemouroidea

In this study, two new mitochondrial genomes (mitogenomes) of Mesonemourametafiligera and Mesonemouratritaenia from the family Nemouridae (Insecta: Plecoptera) were sequenced. The Mesonemourametafiligera mitogenome was a 15,739 bp circular DNA molecule, which was smaller than that of M.tritaenia (15,778 bp) due to differences in the size of the A+T-rich region. Results show that gene content, gene arrangement, base composition, and codon usage were highly conserved in two species. Ka/Ks ratios analyses of protein-coding genes revealed that the highest and lowest rates were found in ND6 and COI and that all these genes were evolving under purifying selection. All tRNA genes in nemourid mitogenomes had a typical cloverleaf secondary structure, except for tRNA^Ser(AGN) which appeared to lack the dihydrouridine arm. The multiple alignments of nemourid lrRNA and srRNA genes showed that sequences of three species were highly conserved. All the A+T-rich region included tandem repeats regions and stem-loop structures. The phylogenetic analyses using Bayesian inference (BI) and maximum likelihood methods (ML) generated identical results. Amphinemurinae and Nemourinae were sister-groups and the family Nemouridae was placed as sister to Capniidae and Taeniopterygidae.

Structure-Related Differences between Cytochrome Oxidase I Proteins in a Stable Heteroplasmic Mitochondrial System

Many bivalve species have two types of mitochondrial DNA passed independently through the female line (F genome) and male line (M genome). Here we study the cytochrome oxidase I protein in such bivalve species and provide evidence for differences between the F and M proteins in amino acid property values, particularly relating to hydrophobicity and helicity. The magnitude of these differences varies between different regions of the protein and the change from the ancestor is most marked in the M protein. The observed changes occur in parallel and in the same direction in the different species studied. Two possible causes are considered, first relaxation of purifying selection with drift and second positive selection. These may operate in different ways in different regions of the protein. Many different amino acid substitutions contribute in a small way to the observed variation, but substitutions involving alanine and serine have a quantitatively large effect. Some of these substitutions are potential targets for phosphorylation and some are close to residues of functional importance in the catalytic mechanism. We propose that the observed changes in the F and M proteins might contribute to functional differences between them relating to ATP production and mitochondrial membrane potential with implications for sperm function.

The Mitonuclear Dimension of Neanderthal and Denisovan Ancestry in Modern Human Genomes

Some human populations interbred with Neanderthals and Denisovans, resulting in substantial contributions to modern-human genomes. Therefore, it is now possible to use genomic data to investigate mechanisms that shaped historical gene flow between humans and our closest hominin relatives. More generally, in eukaryotes, mitonuclear interactions have been argued to play a disproportionate role in generating reproductive isolation. There is no evidence of mtDNA introgression into modern human populations, which means that all introgressed nuclear alleles from archaic hominins must function on a modern-human mitochondrial background. Therefore, mitonuclear interactions are also potentially relevant to hominin evolution. We performed a detailed accounting of mtDNA divergence among hominin lineages and used population-genomic data to test the hypothesis that mitonuclear incompatibilities have preferentially restricted the introgression of nuclear genes with mitochondrial functions. We found a small but significant underrepresentation of introgressed Neanderthal alleles at such nuclear loci. Structural analyses of mitochondrial enzyme complexes revealed that these effects are unlikely to be mediated by physically interacting sites in mitochondrial and nuclear gene products. We did not detect any underrepresentation of introgressed Denisovan alleles at mitochondrial-targeted loci, but this may reflect reduced power because locus-specific estimates of Denisovan introgression are more conservative. Overall, we conclude that genes involved in mitochondrial function may have been subject to distinct selection pressures during the history of introgression from archaic hominins but that mitonuclear incompatibilities have had, at most, a small role in shaping genome-wide introgression patterns, perhaps because of limited functional divergence in mtDNA and interacting nuclear genes.

Compositional and mutational rate heterogeneity in mitochondrial genomes and its effect on the phylogenetic inferences of Cimicomorpha (Hemiptera: Heteroptera)

Background

Mitochondrial genome (mt-genome) data can potentially return artefactual relationships in the higher-level phylogenetic inference of insects due to the biases of accelerated substitution rates and compositional heterogeneity. Previous studies based on mt-genome data alone showed a paraphyly of Cimicomorpha (Insecta, Hemiptera) due to the positions of the families Tingidae and Reduviidae rather than the monophyly that was supported based on morphological characters, morphological and molecular combined data and large scale molecular datasets. Various strategies have been proposed to ameliorate the effects of potential mt-genome biases, including dense taxon sampling, removal of third codon positions or purine-pyrimidine coding and the use of site-heterogeneous models. In this study, we sequenced the mt-genomes of five additional Tingidae species and discussed the compositional and mutational rate heterogeneity in mt-genomes and its effect on the phylogenetic inferences of Cimicomorpha by implementing the bias-reduction strategies mentioned above.

Results

Heterogeneity in nucleotide composition and mutational biases were found in mt protein-coding genes, and the third codon exhibited high levels of saturation. Dense taxon sampling of Tingidae and Reduviidae and the other common strategies mentioned above were insufficient to recover the monophyly of the well-established group Cimicomorpha. When the sites with weak phylogenetic signals in the dataset were removed, the remaining dataset of mt-genomes can support the monophyly of Cimicomorpha; this support demonstrates that mt-genomes possess strong phylogenetic signals for the inference of higher-level phylogeny of this group. Comparison of the ratio of the removal of amino acids for each PCG showed that ATP8 has the highest ratio while CO1 has the lowest. This pattern is largely congruent with the evolutionary rate of 13 PCGs that ATP8 represents the highest evolutionary rate, whereas CO1 appears to be the lowest. Notably, the value of Ka/Ks ratios of all PCGs is less than 1, indicating that these genes are likely evolving under purifying selection.

Conclusions

Our results demonstrate that mt-genomes have sites with strong phylogenetic signals for the inference of higher-level phylogeny of Cimicomorpha. Consequently, bioinformatic approaches to removing sites with weak phylogenetic signals in mt-genome without relying on an a priori tree topology would greatly improve this field.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4650-9) contains supplementary material, which is available to authorized users.

Role of hypoxia in Diffuse Large B-cell Lymphoma: Metabolic repression and selective translation of HK2 facilitates development of DLBCL

Published molecular profiling studies in patients with lymphoma suggested the influence of hypoxia inducible factor-1 alpha (HIF1α) targets in prognosis of DLBCL. Yet, the role of hypoxia in hematological malignancies remains unclear. We observed that activation of HIF1α resulted in global translation repression during hypoxic stress in DLBCL. Protein translation efficiency as measured using 35S-labeled methionine incorporation revealed a ≥50% reduction in translation upon activation of HIF1α. Importantly, translation was not completely inhibited and expression of clinically correlated hypoxia targets such as GLUT1, HK2, and CYT-C was found to be refractory to translational repression under hypoxia in DLBCL cells. Notably, hypoxic induction of these genes was not observed in normal primary B-cells. Translational repression was coupled with a decrease in mitochondrial function. Screening of primary DLBCL patient samples revealed that expression of HK2, which encodes for the enzyme hexokinase 2, was significantly correlated with DLBCL phenotype. Genetic knockdown studies demonstrated that HK2 is required for promoting growth of DLBCL under hypoxic stress. Altogether, our findings provide strong support for the direct contribution of HK2 in B-cell lymphoma development and suggest that HK2 is a key metabolic driver of the DLBCL phenotype.

The molecular evolutionary dynamics of oxidative phosphorylation (OXPHOS) genes in Hymenoptera

BACKGROUND

The primary energy-producing pathway in eukaryotic cells, the oxidative phosphorylation (OXPHOS) system, comprises proteins encoded by both mitochondrial and nuclear genes. To maintain the function of the OXPHOS system, the pattern of substitutions in mitochondrial and nuclear genes may not be completely independent. It has been suggested that slightly deleterious substitutions in mitochondrial genes are compensated by substitutions in the interacting nuclear genes due to positive selection. Among the four largest insect orders, Coleoptera (beetles), Hymenoptera (sawflies, wasps, ants, and bees), Diptera (midges, mosquitoes, and flies) and Lepidoptera (moths and butterflies), the mitochondrial genes of Hymenoptera exhibit an exceptionally high amino acid substitution rate while the evolution of nuclear OXPHOS genes is largely unknown. Therefore, Hymenoptera is an excellent model group for testing the hypothesis of positive selection driving the substitution rate of nuclear OXPHOS genes. In this study, we report the evolutionary rate of OXPHOS genes in Hymenoptera and test for evidence of positive selection in nuclear OXPHOS genes of Hymenoptera.

RESULTS

Our analyses revealed that the amino acid substitution rate of mitochondrial and nuclear OXPHOS genes in Hymenoptera is higher than that in other studied insect orders. In contrast, the amino acid substitution rate of non-OXPHOS genes in Hymenoptera is lower than the rate in other insect orders. Overall, we found the dN/dS ratio of the nuclear OXPHOS genes to be higher in Hymenoptera than in other insect orders. However, nuclear OXPHOS genes with high dN/dS ratio did not always exhibit a high amino acid substitution rate. Using branch-site and site model tests, we identified various codon sites that evolved under positive selection in nuclear OXPHOS genes.

CONCLUSIONS

Our results showed that nuclear OXPHOS genes in Hymenoptera are evolving faster than the genes in other three insect orders. The branch test suggested that while some nuclear OXPHOS genes in Hymenoptera show a signature of positive selection, the pattern is not consistent across all nuclear OXPHOS genes. As only few codon sites were under positive selection, we suggested that positive selection might not be the only factor contributing to the rapid evolution of nuclear OXPHOS genes in Hymenoptera.

Subfunctionalization of COX4 paralogs in fish

Cytochrome c oxidase (COX) subunit 4 has two paralogs in most vertebrates. The mammalian COX4-2 gene is hypoxia responsive, and the protein has a disrupted ATP-binding site that confers kinetic properties on COX that distinguish it from COX4-1. The structure-function of COX4-2 orthologs in other vertebrates remains uncertain. Phylogenetic analyses suggest the two paralogs arose in basal vertebrates, but COX4-2 orthologs diverged faster than COX4-1 orthologs. COX4-1/4-2 protein levels in tilapia tracked mRNA levels across tissues, and did not change in hypoxia, arguing against a role for differential post-translational regulation of paralogs. The heart, and to a lesser extent the brain, showed a size-dependent shift from COX4-1 to COX4-2 (transcript and protein). ATP allosterically inhibited both velocity and affinity for oxygen in COX assayed from both muscle (predominantly COX4-2) and gill (predominantly COX4-1). We saw some evidence of cellular and subcellular discrimination of COX4 paralogs in heart. In cardiac ventricle, some non-cardiomyocyte cells were COX positive but lacked detectible COX4-2. Within heart, the two proteins partitioned to different mitochondrial subpopulations. Cardiac subsarcolemmal mitochondria had mostly COX4-1 and intermyofibrillar mitochondria had mostly COX4-2. Collectively, these data argue that, despite common evolutionary origins, COX4-2 orthologs of fish show unique patterns of subfunctionalization with respect to transcriptional and posttranslation regulation relative to the rodents and primates that have been studied to date.

Mitochondrial-nuclear co-evolution leads to hybrid incompatibility through pentatricopeptide repeat proteins

Mitochondrial-nuclear incompatibility has a major role in reproductive isolation between species. However, the underlying mechanism and driving force of mitochondrial-nuclear incompatibility remain elusive. Here, we report a pentatricopeptide repeat-containing (PPR) protein, Ccm1, and its interacting partner, 15S rRNA, to be involved in hybrid incompatibility between two yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus S. bayanus-Ccm1 has reduced binding affinity for S. cerevisiae-15S rRNA, leading to respiratory defects in hybrid cells. This incompatibility can be rescued by single mutations on several individual PPR motifs, demonstrating the highly evolvable nature of PPR proteins. When we examined other PPR proteins in the closely related Saccharomyces sensu stricto yeasts, about two-thirds of them showed detectable incompatibility. Our results suggest that fast co-evolution between flexible PPR proteins and their mitochondrial RNA substrates may be a common driving force in the development of mitochondrial-nuclear hybrid incompatibility.

A Mitochondrial Genome of Rhyparochromidae (Hemiptera: Heteroptera) and a Comparative Analysis of Related Mitochondrial Genomes

The Rhyparochromidae, the largest family of Lygaeoidea, encompasses more than 1,850 described species, but no mitochondrial genome has been sequenced to date. Here we describe the first mitochondrial genome for Rhyparochromidae: a complete mitochondrial genome of Panaorus albomaculatus (Scott, 1874). This mitochondrial genome is comprised of 16,345 bp, and contains the expected 37 genes and control region. The majority of the control region is made up of a large tandem-repeat region, which has a novel pattern not previously observed in other insects. The tandem-repeats region of P. albomaculatus consists of 53 tandem duplications (including one partial repeat), which is the largest number of tandem repeats among all the known insect mitochondrial genomes. Slipped-strand mispairing during replication is likely to have generated this novel pattern of tandem repeats. Comparative analysis of tRNA gene families in sequenced Pentatomomorpha and Lygaeoidea species shows that the pattern of nucleotide conservation is markedly higher on the J-strand. Phylogenetic reconstruction based on mitochondrial genomes suggests that Rhyparochromidae is not the sister group to all the remaining Lygaeoidea, and supports the monophyly of Lygaeoidea.

COX7AR is a Stress-inducible Mitochondrial COX Subunit that Promotes Breast Cancer Malignancy

Cytochrome c oxidase (COX), the terminal enzyme of the mitochondrial respiratory chain, plays a key role in regulating mitochondrial energy production and cell survival. COX subunit VIIa polypeptide 2-like protein (COX7AR) is a novel COX subunit that was recently found to be involved in mitochondrial supercomplex assembly and mitochondrial respiration activity. Here, we report that COX7AR is expressed in high energy-demanding tissues, such as brain, heart, liver, and aggressive forms of human breast cancer cells. Under cellular stress that stimulates energy metabolism, COX7AR is induced and incorporated into the mitochondrial COX complex. Functionally, COX7AR promotes cellular energy production in human mammary epithelial cells. Gain- and loss-of-function analysis demonstrates that COX7AR is required for human breast cancer cells to maintain higher rates of proliferation, clone formation, and invasion. In summary, our study revealed that COX7AR is a stress-inducible mitochondrial COX subunit that facilitates human breast cancer malignancy. These findings have important implications in the understanding and treatment of human breast cancer and the diseases associated with mitochondrial energy metabolism.

Plastid-Nuclear Interaction and Accelerated Coevolution in Plastid Ribosomal Genes in Geraniaceae

Plastids and mitochondria have many protein complexes that include subunits encoded by organelle and nuclear genomes. In animal cells, compensatory evolution between mitochondrial and nuclear-encoded subunits was identified and the high mitochondrial mutation rates were hypothesized to drive compensatory evolution in nuclear genomes. In plant cells, compensatory evolution between plastid and nucleus has rarely been investigated in a phylogenetic framework. To investigate plastid–nuclear coevolution, we focused on plastid ribosomal protein genes that are encoded by plastid and nuclear genomes from 27 Geraniales species. Substitution rates were compared for five sets of genes representing plastid- and nuclear-encoded ribosomal subunit proteins targeted to the cytosol or the plastid as well as nonribosomal protein controls. We found that nonsynonymous substitution rates (dN) and the ratios of nonsynonymous to synonymous substitution rates (ω) were accelerated in both plastid- (CpRP) and nuclear-encoded subunits (NuCpRP) of the plastid ribosome relative to control sequences. Our analyses revealed strong signals of cytonuclear coevolution between plastid- and nuclear-encoded subunits, in which nonsynonymous substitutions in CpRP and NuCpRP tend to occur along the same branches in the Geraniaceae phylogeny. This coevolution pattern cannot be explained by physical interaction between amino acid residues. The forces driving accelerated coevolution varied with cellular compartment of the sequence. Increased ω in CpRP was mainly due to intensified positive selection whereas increased ω in NuCpRP was caused by relaxed purifying selection. In addition, the many indels identified in plastid rRNA genes in Geraniaceae may have contributed to changes in plastid subunits.

Conservative and compensatory evolution in oxidative phosphorylation complexes of angiosperms with highly divergent rates of mitochondrial genome evolution

Interactions between nuclear and mitochondrial gene products are critical for eukaryotic cell function. Nuclear genes encoding mitochondrial-targeted proteins (N-mt genes) experience elevated rates of evolution, which has often been interpreted as evidence of nuclear compensation in response to elevated mitochondrial mutation rates. However, N-mt genes may be under relaxed functional constraints, which could also explain observed increases in their evolutionary rate. To disentangle these hypotheses, we examined patterns of sequence and structural evolution in nuclear- and mitochondrial-encoded oxidative phosphorylation proteins from species in the angiosperm genus Silene with vastly different mitochondrial mutation rates. We found correlated increases in N-mt gene evolution in species with fast-evolving mitochondrial DNA. Structural modeling revealed an overrepresentation of N-mt substitutions at positions that directly contact mutated residues in mitochondrial-encoded proteins, despite overall patterns of conservative structural evolution. These findings support the hypothesis that selection for compensatory changes in response to mitochondrial mutations contributes to the elevated rate of evolution in N-mt genes. We discuss these results in light of theories implicating mitochondrial mutation rates and mitonuclear coevolution as drivers of speciation and suggest comparative and experimental approaches that could take advantage of heterogeneity in rates of mtDNA evolution across eukaryotes to evaluate such theories.

Heterogeneous natural selection on oxidative phosphorylation genes among fishes with extreme high and low aerobic performance

Background

Oxidative phosphorylation (OXPHOS) is the primary source of ATP in eukaryotes and serves as a mechanistic link between variation in genotypes and energetic phenotypes. While several physiological and anatomical factors may lead to increased aerobic capacity, variation in OXPHOS proteins may influence OXPHOS efficiency and facilitate adaptation in organisms with varied energy demands. Although there is evidence that natural selection acts on OXPHOS genes, the focus has been on detection of directional (positive) selection on specific phylogenetic branches where traits that increase energetic demands appear to have evolved. We examined patterns of selection in a broader evolutionary context, i.e., on multiple lineages of fishes with extreme high and low aerobic performance.

Results

We found that patterns of natural selection on mitochondrial OXPHOS genes are complex among fishes with different swimming performance. Positive selection is not consistently associated with high performance taxa and appears to be strongest on lineages containing low performance taxa. In contrast, within high performance lineages, purifying (negative) selection appears to predominate.

Conclusions

We provide evidence that selection on OXPHOS varies in both form and intensity within and among lineages through evolutionary time. These results provide evidence for fluctuating selection on OXPHOS associated with divergence in aerobic performance. However, in contrast to previous studies, positive selection was strongest on low performance taxa suggesting that adaptation of OXPHOS involves many factors beyond enhancing ATP production in high performance taxa. The broader pattern indicates a complex interplay between organismal adaptations, ATP demand, and OXPHOS function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0453-7) contains supplementary material, which is available to authorized users.

Sex-specific effects of sympatric mitonuclear variation on fitness in Drosophila subobscura

Background

A number of recent studies have shown that the pattern of mitochondrial DNA variation and evolution is at odds with a neutral equilibrium model. Theory has suggested that selection on mitonuclear genotypes can act to maintain stable mitonuclear polymorphism within populations. However, this effect largely relies upon selection being either sex-specific or frequency dependent. Here, we use mitonuclear introgression lines to assess differences in a series of key life-history traits (egg-to-adult developmental time, viability, offspring sex-ratio, adult longevity and resistance to desiccation) in Drosophila subobscura fruit flies carrying one of three different sympatric mtDNA haplotypes.

Results

We found functional differences between these sympatric mtDNA haplotypes, but these effects were contingent upon the nuclear genome with which they were co-expressed. Further, we demonstrate a significant mitonuclear genetic effect on adult sex ratio, as well as a sex × mtDNA × nuDNA interaction for adult longevity.

Conclusions

The observed effects suggest that sex specific mitonuclear selection contributes to the maintenance of mtDNA polymorphism and to mitonuclear linkage disequilibrium in this model system.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0421-2) contains supplementary material, which is available to authorized users.

Fast fish face fewer mitochondrial mutations: Patterns of dN/dS across fish mitogenomes

Mitochondrial DNA is routinely used to answer a variety of biological questions; and there is growing evidence suggesting that its accumulation of mutations is influenced by life history, effective population size and cellular energy requirements. This study examines the influence of phylogenetic patterns of metabolic activity on the evolution of mitochondrial DNA in fishes, given energy requirements associated with high performance versus sedentary life histories. It was determined that all 13 protein coding genes of the mitogenome experience a relaxation of purifying selection in sedentary fishes. This phenomenon was not detected in nuclear housekeeping genes, suggesting that it can be explained by the energy requirements of these groups, and possibly their effective population sizes. This study also examined the subunit binding sites of two subunits of cytochrome c oxidase (COXI and COXIII), and did not detect any differences in selection between these groups of fishes. These cytochrome c oxidase subunits interact with subunits that are encoded by the nuclear genome and it has been suggested that a unique form of coevolution occurs between these genomes in order to maintain function, and may have implications for speciation. Although this was not a main focus of this study, our preliminary results suggest that substitutions in subunit binding site regions are rare. The results from this study add to the growing literature on the complex relationship between mitochondrial DNA and the evolution of life histories across the tree of life.

Hemipteran mitochondrial genomes: features, structures and implications for phylogeny

The study of Hemipteran mitochondrial genomes (mitogenomes) began with the Chagas disease vector, Triatoma dimidiata, in 2001. At present, 90 complete Hemipteran mitogenomes have been sequenced and annotated. This review examines the history of Hemipteran mitogenomes research and summarizes the main features of them including genome organization, nucleotide composition, protein-coding genes, tRNAs and rRNAs, and non-coding regions. Special attention is given to the comparative analysis of repeat regions. Gene rearrangements are an additional data type for a few families, and most mitogenomes are arranged in the same order to the proposed ancestral insect. We also discuss and provide insights on the phylogenetic analyses of a variety of taxonomic levels. This review is expected to further expand our understanding of research in this field and serve as a valuable reference resource.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

Mito-nuclear co-evolution: the positive and negative sides of functional ancient mutations

Most cell functions are carried out by interacting factors, thus underlying the functional importance of genetic interactions between genes, termed epistasis. Epistasis could be under strong selective pressures especially in conditions where the mutation rate of one of the interacting partners notably differs from the other. Accordingly, the order of magnitude higher mitochondrial DNA (mtDNA) mutation rate as compared to the nuclear DNA (nDNA) of all tested animals, should influence systems involving mitochondrial-nuclear (mito-nuclear) interactions. Such is the case of the energy producing oxidative phosphorylation (OXPHOS) and mitochondrial translational machineries which are comprised of factors encoded by both the mtDNA and the nDNA. Additionally, the mitochondrial RNA transcription and mtDNA replication systems are operated by nDNA-encoded proteins that bind mtDNA regulatory elements. As these systems are central to cell life there is strong selection toward mito-nuclear co-evolution to maintain their function. However, it is unclear whether (A) mito-nuclear co-evolution befalls only to retain mitochondrial functions during evolution or, also, (B) serves as an adaptive tool to adjust for the evolving energetic demands as species' complexity increases. As the first step to answer these questions we discuss evidence of both negative and adaptive (positive) selection acting on the mtDNA and nDNA-encoded genes and the effect of both types of selection on mito-nuclear interacting factors. Emphasis is given to the crucial role of recurrent ancient (nodal) mutations in such selective events. We apply this point-of-view to the three available types of mito-nuclear co-evolution: protein-protein (within the OXPHOS system), protein-RNA (mainly within the mitochondrial ribosome), and protein-DNA (at the mitochondrial replication and transcription machineries).

Protein-protein interfaces from cytochrome c oxidase I evolve faster than nonbinding surfaces, yet negative selection is the driving force

Respiratory complexes are encoded by two genomes (mitochondrial DNA [mtDNA] and nuclear DNA [nDNA]). Although the importance of intergenomic coadaptation is acknowledged, the forces and constraints shaping such coevolution are largely unknown. Previous works using cytochrome c oxidase (COX) as a model enzyme have led to the so-called “optimizing interaction” hypothesis. According to this view, mtDNA-encoded residues close to nDNA-encoded residues evolve faster than the rest of positions, favoring the optimization of protein–protein interfaces. Herein, using evolutionary data in combination with structural information of COX, we show that failing to discern the effects of interaction from other structural and functional effects can lead to deceptive conclusions such as the “optimizing hypothesis.” Once spurious factors have been accounted for, data analysis shows that mtDNA-encoded residues engaged in contacts are, in general, more constrained than their noncontact counterparts. Nevertheless, noncontact residues from the surface of COX I subunit are a remarkable exception, being subjected to an exceptionally high purifying selection that may be related to the maintenance of a suitable heme environment. We also report that mtDNA-encoded residues involved in contacts with other mtDNA-encoded subunits are more constrained than mtDNA-encoded residues interacting with nDNA-encoded polypeptides. This differential behavior cannot be explained on the basis of predicted thermodynamic stability, as interactions between mtDNA-encoded subunits contribute more weakly to the complex stability than those interactions between subunits encoded by different genomes. Therefore, the higher conservation observed among mtDNA-encoded residues involved in intragenome interactions is likely due to factors other than structural stability.

Evolutionary perspectives on the links between mitochondrial genotype and disease phenotype

BACKGROUND

Disorders of the mitochondrial respiratory chain are heterogeneous in their symptoms and underlying genetics. Simple links between candidate mutations and expression of disease phenotype typically do not exist. It thus remains unclear how the genetic variation in the mitochondrial genome contributes to the phenotypic expression of complex traits and disease phenotypes.

SCOPE OF REVIEW

I summarize the basic genetic processes known to underpin mitochondrial disease. I highlight other plausible processes, drawn from the evolutionary biological literature, whose contribution to mitochondrial disease expression remains largely empirically unexplored. I highlight recent advances to the field, and discuss common-ground and -goals shared by researchers across medical and evolutionary domains.

MAJOR CONCLUSIONS

Mitochondrial genetic variance is linked to phenotypic variance across a variety of traits (e.g. reproductive function, life expectancy) fundamental to the upkeep of good health. Evolutionary theory predicts that mitochondrial genomes are destined to accumulate male-harming (but female-friendly) mutations, and this prediction has received proof-of-principle support. Furthermore, mitochondrial effects on the phenotype are typically manifested via interactions between mitochondrial and nuclear genes. Thus, whether a mitochondrial mutation is pathogenic in effect can depend on the nuclear genotype in which is it expressed.

GENERAL SIGNIFICANCE

Many disease phenotypes associated with OXPHOS malfunction might be determined by the outcomes of mitochondrial-nuclear interactions, and by the evolutionary forces that historically shaped mitochondrial DNA (mtDNA) sequences. Concepts and results drawn from the evolutionary sciences can have broad, but currently under-utilized, applicability to the medical sciences and provide new insights into understanding the complex genetics of mitochondrial disease. This article is part of a Special Issue entitled Frontiers of Mitochondrial Research.

Bioenergetics in human evolution and disease: implications for the origins of biological complexity and the missing genetic variation of common diseases

Two major inconsistencies exist in the current neo-Darwinian evolutionary theory that random chromosomal mutations acted on by natural selection generate new species. First, natural selection does not require the evolution of ever increasing complexity, yet this is the hallmark of biology. Second, human chromosomal DNA sequence variation is predominantly either neutral or deleterious and is insufficient to provide the variation required for speciation or for predilection to common diseases. Complexity is explained by the continuous flow of energy through the biosphere that drives the accumulation of nucleic acids and information. Information then encodes complex forms. In animals, energy flow is primarily mediated by mitochondria whose maternally inherited mitochondrial DNA (mtDNA) codes for key genes for energy metabolism. In mammals, the mtDNA has a very high mutation rate, but the deleterious mutations are removed by an ovarian selection system. Hence, new mutations that subtly alter energy metabolism are continuously introduced into the species, permitting adaptation to regional differences in energy environments. Therefore, the most phenotypically significant gene variants arise in the mtDNA, are regional, and permit animals to occupy peripheral energy environments where rarer nuclear DNA (nDNA) variants can accumulate, leading to speciation. The neutralist-selectionist debate is then a consequence of mammals having two different evolutionary strategies: a fast mtDNA strategy for intra-specific radiation and a slow nDNA strategy for speciation. Furthermore, the missing genetic variation for common human diseases is primarily mtDNA variation plus regional nDNA variants, both of which have been missed by large, inter-population association studies.

Structural analysis of mitochondrial mutations reveals a role for bigenomic protein interactions in human disease

Mitochondria are the energy producing organelles of the cell, and mutations within their genome can cause numerous and often severe human diseases. At the heart of every mitochondrion is a set of five large multi-protein machines collectively known as the mitochondrial respiratory chain (MRC). This cellular machinery is central to several processes important for maintaining homeostasis within cells, including the production of ATP. The MRC is unique due to the bigenomic origin of its interacting proteins, which are encoded in the nucleus and mitochondria. It is this, in combination with the sheer number of protein-protein interactions that occur both within and between the MRC complexes, which makes the prediction of function and pathological outcome from primary sequence mutation data extremely challenging. Here we demonstrate how 3D structural analysis can be employed to predict the functional importance of mutations in mtDNA protein-coding genes. We mined the MITOMAP database and, utilizing the latest structural data, classified mutation sites based on their location within the MRC complexes III and IV. Using this approach, four structural classes of mutation were identified, including one underexplored class that interferes with nuclear-mitochondrial protein interactions. We demonstrate that this class currently eludes existing predictive approaches that do not take into account the quaternary structural organization inherent within and between the MRC complexes. The systematic and detailed structural analysis of disease-associated mutations in the mitochondrial Complex III and IV genes significantly enhances the predictive power of existing approaches and our understanding of how such mutations contribute to various pathologies. Given the general lack of any successful therapeutic approaches for disorders of the MRC, these findings may inform the development of new diagnostic and prognostic biomarkers, as well as new drugs and targets for gene therapy.

Animal mitochondria, positive selection and cyto-nuclear coevolution: insights from pulmonates

Pulmonate snails have remarkably high levels of mtDNA polymorphism within species and divergence between species, making them an interesting group for the study of mutation and selection on mitochondrial genomes. The availability of sequence data from most major lineages – collected largely for studies of phylogeography - provides an opportunity to perform several tests of selection that may provide general insights into the evolutionary forces that have produced this unusual pattern. Several protein coding mtDNA datasets of pulmonates were analyzed towards this direction. Two different methods for the detection of positive selection were used, one based on phylogeny, and the other on the McDonald-Kreitman test. The cyto-nuclear coevolution hypothesis, often implicated to account for the high levels of mtDNA divergence of some organisms, was also addressed by assessing the divergence pattern exhibited by a nuclear gene. The McDonald-Kreitman test indicated multiple signs of positive selection in the mtDNA genes, but was significantly biased when sequence divergence was high. The phylogenetic method identified five mtDNA datasets as affected by positive selection. In the nuclear gene, the McDonald-Kreitman test provided no significant results, whereas the phylogenetic method identified positive selection as likely present. Overall, our findings indicate that: 1) slim support for the cyto-nuclear coevolution hypothesis is present, 2) the elevated rates of mtDNA polymorphims and divergence in pulmonates do not appear to be due to pervasive positive selection, 3) more stringent tests show that spurious positive selection is uncovered when distant taxa are compared and 4) there are significant examples of positive selection acting in some cases, so it appears that mtDNA evolution in pulmonates can escape from strict deleterious evolution suggested by the Muller’s ratchet effect.

Comparative mitogenomic analysis of damsel bugs representing three tribes in the family Nabidae (Insecta: Hemiptera)

BACKGROUND

Nabidae, a family of predatory heteropterans, includes two subfamilies and five tribes. We previously reported the complete mitogenome of Alloeorhynchus bakeri, a representative of the tribe Prostemmatini in the subfamily Prostemmatinae. To gain a better understanding of architecture and evolution of mitogenome in Nabidae, mitogenomes of five species representing two tribes (Gorpini and Nabini) in the subfamily Nabinae were sequenced, and a comparative mitogenomic analysis of three nabid tribes in two subfamilies was carried out.

METHODOLOGY/PRINCIPAL FINDINGS

Nabid mitogenomes share a similar nucleotide composition and base bias, except for the control region, where differences are observed at the subfamily level. In addition, the pattern of codon usage is influenced by the GC content and consistent with the standard invertebrate mitochondrial genetic code and the preference for A+T-rich codons. The comparison among orthologous protein-coding genes shows that different genes have been subject to different rates of molecular evolution correlated with the GC content. The stems and anticodon loops of tRNAs are extremely conserved, and the nucleotide substitutions are largely restricted to TψC and DHU loops and extra arms, with insertion-deletion polymorphisms. Comparative analysis shows similar rates of substitution between the two rRNAs. Long non-coding regions are observed in most Gorpini and Nabini mtDNAs in-between trnI-trnQ and/or trnS2-nad1. The lone exception, Nabis apicalis, however, has lost three tRNAs. Overall, phylogenetic analysis using mitogenomic data is consistent with phylogenies constructed mainly form morphological traits.

CONCLUSIONS/SIGNIFICANCE

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the family Nabidae. Nucleotide diversity and mitogenomic traits are phylogenetically informative at subfamily level. Furthermore, inclusion of a broader range of samples representing various taxonomic levels is critical for the understanding of mitogenomic evolution in damsel bugs.

Mutational bias plays an important role in shaping longevity-related amino acid content in mammalian mtDNA-encoded proteins

During the course of evolution, amino acid shifts might have resulted in mitochondrial proteomes better endowed to resist oxidative stress. However, owing to the problem of distinguishing between functional constraints/adaptations in protein sequences and mutation-driven biases in the composition of these sequences, the adaptive value of such amino acid shifts remains under discussion. We have analyzed the coding sequences of mtDNA from 173 mammalian species, dissecting the effect of nucleotide composition on amino acid usages. We found remarkable cysteine avoidance in mtDNA-encoded proteins. However, no effect of longevity on cysteine content could be detected. On the other hand, nucleotide compositional shifts fully accounted for threonine usages. In spite of a strong effect of mutational bias on methionine abundances, our results suggest a role of selection in determining the composition of methionine. Whether this selective effect is linked or not to protection against oxidative stress is still a subject of debate.

Thermodynamic stability explains the differential evolutionary dynamics of cytochrome b and COX I in mammals

By using a combination of evolutionary and structural data from 231 species, we have addressed the relationship between evolution and structural features of cytochrome b and COX I, two mtDNA-encoded proteins. The interior of cytochrome b, in contrast to that of COX I, exhibits a remarkable tolerance to changes. The higher evolvability of cytochrome b contrasts with the lower rate of synonymous substitutions of its gene when compared to that of COX I, suggesting that the latter is subjected to a stronger purifying selection. We present evidences that the stability effect of mutations (ΔΔG) may be behind these differential behaviour.

Mitochondrial DNA sequence variation associated with dementia and cognitive function in the elderly

Mitochondrial dysfunction is a prominent hallmark of Alzheimer's disease (AD). Mitochondrial DNA (mtDNA) damage may be a major cause of abnormal reactive oxidative species production in AD or increased neuronal susceptibility to oxidative injury during aging. The purpose of this study was to assess the influence of mtDNA sequence variation on clinically significant cognitive impairment and dementia risk in the population-based Health, Aging, and Body Composition (Health ABC) Study. We first investigated the role of common mtDNA haplogroups and individual variants on dementia risk and 8-year change on the Modified Mini-Mental State Examination (3MS) and Digit Symbol Substitution Test (DSST) among 1,631 participants of European genetic ancestry. Participants were free of dementia at baseline and incidence was determined in 273 cases from hospital and medication records over 10-12 follow-up years. Participants from haplogroup T had a statistically significant increased risk of developing dementia (OR = 1.86, 95% CI = 1.23, 2.82, p = 0.0008) and haplogroup J participants experienced a statistically significant 8-year decline in 3MS (β = -0.14, 95% CI = -0.27, -0.03, p = 0.0006), both compared with common haplogroup H. The m.15244A>G, p.G166G, CytB variant was associated with a significant decline in DSST score (β = -0.58, 95% CI -0.89, -0.28, p = 0.00019) and the m.14178T>C, p.I166V, ND6 variant was associated with a significant decline in 3MS score (β = -0.87, 95% CI -1.31, -3.86, p = 0.00012). Finally, we sequenced the complete ~16.5 kb mtDNA from 135 Health ABC participants and identified several highly conserved and potentially functional nonsynonymous variants unique to 22 dementia cases and aggregate sequence variation across the hypervariable 2-3 regions that influences 3MS and DSST scores.

Is the COI barcoding gene involved in speciation through intergenomic conflict?

We here test the proposition that changes in the barcoding region of COI are commonly involved in speciation through intergenomic conflict. We demonstrate that this is unlikely given that even with incomplete taxon sampling, 78-90% of closely-related animal species have identical COI amino acid sequences. In addition, in those cases where amino acid substitutions between closely related species are observed, the inter- and intra-specific substitution patterns are very similar and/or lack consistent differences in the number, position and type of amino acid change. Overall, we conclude that there is little evidence for a widespread involvement of the barcoding gene in speciation.

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.

Silencing, positive selection and parallel evolution: busy history of primate cytochromes C

Cytochrome c (cyt c) participates in two crucial cellular processes, energy production and apoptosis, and unsurprisingly is a highly conserved protein. However, previous studies have reported for the primate lineage (i) loss of the paralogous testis isoform, (ii) an acceleration and then a deceleration of the amino acid replacement rate of the cyt c somatic isoform, and (iii) atypical biochemical behavior of human cyt c. To gain insight into the cause of these major evolutionary events, we have retraced the history of cyt c loci among primates. For testis cyt c, all primate sequences examined carry the same nonsense mutation, which suggests that silencing occurred before the primates diversified. For somatic cyt c, maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses yielded the same tree topology. The evolutionary analyses show that a fast accumulation of non-synonymous mutations (suggesting positive selection) occurred specifically on the anthropoid lineage root and then continued in parallel on the early catarrhini and platyrrhini stems. Analysis of evolutionary changes using the 3D structure suggests they are focused on the respiratory chain rather than on apoptosis or other cyt c functions. In agreement with previous biochemical studies, our results suggest that silencing of the cyt c testis isoform could be linked with the decrease of primate reproduction rate. Finally, the evolution of cyt c in the two sister anthropoid groups leads us to propose that somatic cyt c evolution may be related both to COX evolution and to the convergent brain and body mass enlargement in these two anthropoid clades.

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.

Cytochrome c oxidase: evolution of control via nuclear subunit addition

According to theory, present eukaryotic cells originated from a beneficial association between two free-living cells. Due to this endosymbiotic event the pre-eukaryotic cell gained access to oxidative phosphorylation (OXPHOS), which produces more than 15 times as much ATP as glycolysis. Because cellular ATP needs fluctuate and OXPHOS both requires and produces entities that can be toxic for eukaryotic cells such as ROS or NADH, we propose that the success of endosymbiosis has largely depended on the regulation of endosymbiont OXPHOS. Several studies have presented cytochrome c oxidase as a key regulator of OXPHOS; for example, COX is the only complex of mammalian OXPHOS with known tissue-specific isoforms of nuclear encoded subunits. We here discuss current knowledge about the origin of nuclear encoded subunits and the appearance of different isozymes promoted by tissue and cellular environments such as hypoxia. We also review evidence for recent selective pressure acting on COX among vertebrates, particularly in primate lineages, and discuss the unique pattern of co-evolution between the nuclear and mitochondrial genomes. Finally, even though the addition of nuclear encoded subunits was a major event in eukaryotic COX evolution, this does not lead to emergence of a more efficient COX, as might be expected from an anthropocentric point of view, for the "higher" organism possessing large brains and muscles. The main function of these subunits appears to be "only" to control the activity of the mitochondrial subunits. We propose that this control function is an as yet under appreciated key point of evolution. Moreover, the importance of regulating energy supply may have caused the addition of subunits encoded by the nucleus in a process comparable to a "domestication scenario" such that the host tends to control more and more tightly the ancestral activity of COX performed by the mtDNA encoded subunits.

Searching for evidence of selection in avian DNA barcodes

The barcode of life project has assembled a tremendous number of mitochondrial cytochrome c oxidase I (COI) sequences. Although these sequences were gathered to develop a DNA-based system for species identification, it has been suggested that further biological inferences may also be derived from this wealth of data. Recurrent selective sweeps have been invoked as an evolutionary mechanism to explain limited intraspecific COI diversity, particularly in birds, but this hypothesis has not been formally tested. In this study, I collated COI sequences from previous barcoding studies on birds and tested them for evidence of selection. Using this expanded data set, I re-examined the relationships between intraspecific diversity and interspecific divergence and sampling effort, respectively. I employed the McDonald-Kreitman test to test for neutrality in sequence evolution between closely related pairs of species. Because amino acid sequences were generally constrained between closely related pairs, I also included broader intra-order comparisons to quantify patterns of protein variation in avian COI sequences. Lastly, using 22 published whole mitochondrial genomes, I compared the evolutionary rate of COI against the other 12 protein-coding mitochondrial genes to assess intragenomic variability. I found no conclusive evidence of selective sweeps. Most evidence pointed to an overall trend of strong purifying selection and functional constraint. The COI protein did vary across the class Aves, but to a very limited extent. COI was the least variable gene in the mitochondrial genome, suggesting that other genes might be more informative for probing factors constraining mitochondrial variation within species.

Evolutionary patterns of the mitochondrial genome in Metazoa: exploring the role of mutation and selection in mitochondrial protein coding genes

The mitochondrial genome is a fundamental component of the eukaryotic domain of life, encoding for several important subunits of the respiratory chain, the main energy production system in cells. The processes by means of which mitochondrial DNA (mtDNA) replicates, expresses itself and evolves have been explored over the years, although various aspects are still debated. In this review, we present several key points in modern research on the role of evolutionary forces in affecting mitochondrial genomes in Metazoa. In particular, we assemble the main data on their evolution, describing the contributions of mutational pressure, purifying, and adaptive selection, and how they are related. We also provide data on the evolutionary fate of the mitochondrial synonymous variation, related to the nonsynonymous variation, in comparison with the pattern detected in the nucleus.

Elevated mutational pressure characterizes the evolution of the mitochondrial synonymous variation, whereas purging selection, physiologically due to phenomena such as cell atresia and intracellular mtDNA selection, guarantees coding sequence functionality. This enables mitochondrial adaptive mutations to emerge and fix in the population, promoting mitonuclear coevolution.

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.

Genetical and molecular analysis reveals a cooperating relationship between cytoplasmic male sterility- and fertility restoration-related genes in Oryza species

Although the characterization of genes associated with cytoplasmic male sterility (CMS) and fertility restoration (Rf) has been well documented, the evolutionary relationship between nuclear Rf and CMS factors in mitochondria in Oryza species is still less understood. Here, 41 accessions from 7 Oryza species with AA genome were employed for analyzing the evolutionary relationships between the CMS factors and Rf candidates on chromosome 10. The phylogenetic tree based on restriction fragment length polymorphism patterns of CMS-associated mitochondrial genes showed that these 41 Oryza accessions fell into 3 distinct groups. Another phylogenetic tree based on PCR profiles of the nuclear Rf candidates on chromosome 10 was also established, and three groups were distinctively grouped. The accessions in each subgroup/group of the two phylogenetic trees are well parallel to each other. Furthermore, the 41 investigated accessions were test-crossed with Honglian (gametophytic type) and Wild-abortive (sporophytic type) CMS, and 5 groups were classified according to their restoring ability. The accessions in the same subgroup of the two phylogenetic trees shared similar fertility restoring pattern. Therefore, we conclude that the CMS-associated mitotypes are compatible to the Rf candidate-related nucleotypes, CMS and Rf have a parallel evolutionary relation in the Oryza species.

Distinct patterns of mitochondrial genome diversity in bonobos (Pan paniscus) and humans

BACKGROUND

We have analyzed the complete mitochondrial genomes of 22 Pan paniscus (bonobo, pygmy chimpanzee) individuals to assess the detailed mitochondrial DNA (mtDNA) phylogeny of this close relative of Homo sapiens.

RESULTS

We identified three major clades among bonobos that separated approximately 540,000 years ago, as suggested by Bayesian analysis. Incidentally, we discovered that the current reference sequence for bonobo likely is a hybrid of the mitochondrial genomes of two distant individuals. When comparing spectra of polymorphic mtDNA sites in bonobos and humans, we observed two major differences: (i) Of all 31 bonobo mtDNA homoplasies, i.e. nucleotide changes that occurred independently on separate branches of the phylogenetic tree, 13 were not homoplasic in humans. This indicates that at least a part of the unstable sites of the mitochondrial genome is species-specific and difficult to be explained on the basis of a mutational hotspot concept. (ii) A comparison of the ratios of non-synonymous to synonymous changes (dN/dS) among polymorphic positions in bonobos and in 4902 Homo sapiens mitochondrial genomes revealed a remarkable difference in the strength of purifying selection in the mitochondrial genes of the F0F1-ATPase complex. While in bonobos this complex showed a similar low value as complexes I and IV, human haplogroups displayed 2.2 to 7.6 times increased dN/dS ratios when compared to bonobos.

CONCLUSIONS

Some variants of mitochondrially encoded subunits of the ATPase complex in humans very likely decrease the efficiency of energy conversion leading to production of extra heat. Thus, we hypothesize that the species-specific release of evolutionary constraints for the mitochondrial genes of the proton-translocating ATPase is a consequence of altered heat homeostasis in modern humans.

Recent studies of iron deficiency during brain development in nonhuman primates

Recent studies of the effects of developmental iron deficiency (ID) and iron deficiency anemia in nonhuman primates have provided new insights into this widespread and well-recognized human nutritional deficiency. The rhesus monkey was the animal model in these experiments, which used extensive hematological and behavioral evaluations in addition to noninvasive brain measures. Two important findings were as follows: 1) different behavioral consequences depending on the timing of ID relative to brain developmental stages and 2) the potential for long-lasting changes in brain iron regulatory systems. Further work in this model, including integration with studies in humans and in laboratory rodents, is ongoing.