Coevolution predicts direct interactions between mtDNA-encoded and nDNA-encoded subunits of oxidative phosphorylation complex i

Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions

Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Genomic Signatures of Mitonuclear Coevolution in Mammals

Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

A Primer Genetic Toolkit for Exploring Mitochondrial Biology and Disease Using Zebrafish

Mitochondria are a dynamic eukaryotic innovation that play diverse roles in biology and disease. The mitochondrial genome is remarkably conserved in all vertebrates, encoding the same 37-gene set and overall genomic structure, ranging from 16,596 base pairs (bp) in the teleost zebrafish (Danio rerio) to 16,569 bp in humans. Mitochondrial disorders are amongst the most prevalent inherited diseases, affecting roughly 1 in every 5000 individuals. Currently, few effective treatments exist for those with mitochondrial ailments, representing a major unmet patient need. Mitochondrial dysfunction is also a common component of a wide variety of other human illnesses, ranging from neurodegenerative disorders such as Huntington's disease and Parkinson's disease to autoimmune illnesses such as multiple sclerosis and rheumatoid arthritis. The electron transport chain (ETC) component of mitochondria is critical for mitochondrial biology and defects can lead to many mitochondrial disease symptoms. Here, we present a publicly available collection of genetic mutants created in highly conserved, nuclear-encoded mitochondrial genes in Danio rerio. The zebrafish system represents a potentially powerful new opportunity for the study of mitochondrial biology and disease due to the large number of orthologous genes shared with humans and the many advanced features of this model system, from genetics to imaging. This collection includes 15 mutant lines in 13 different genes created through locus-specific gene editing to induce frameshift or splice acceptor mutations, leading to predicted protein truncation during translation. Additionally, included are 11 lines created by the random insertion of the gene-breaking transposon (GBT) protein trap cassette. All these targeted mutant alleles truncate conserved domains of genes critical to the proper function of the ETC or genes that have been implicated in human mitochondrial disease. This collection is designed to accelerate the use of zebrafish to study many different aspects of mitochondrial function to widen our understanding of their role in biology and human disease.

Identification of a novel m.3955G > A variant in MT-ND1 associated with Leigh syndrome

Leigh syndrome (LS) is one of the most common mitochondrial diseases in children, for which at least 90 causative genes have been identified. However, many LS patients have no genetic diagnosis, indicating that more disease-related genes remain to be identified. In this study, we identified a novel variant, m.3955G > A, in mitochondrially encoded NADH:ubiquinone oxidoreductase core subunit 1 (MT-ND1) in two unrelated LS patients, manifesting as infancy-onset frequent seizures, neurodegeneration, elevated lactate levels, and bilateral symmetrical lesions in the brainstem, basal ganglia, and thalamus. Transfer of the mutant mtDNA with m.3955G > A into cybrids disturbed the MT-ND1 expression and CI assembly, followed by remarkable mitochondrial dysfunction, reactive oxygen species production, and mitochondrial membrane potential reduction. Our findings demonstrated the pathogenicity of the novel m.3955G > A variant, and extend the spectrum of pathogenic mtDNA variants.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Bivalve molluscs as model systems for studying mitochondrial biology

The class Bivalvia is a highly successful and ancient taxon including ∼25,000 living species. During their long evolutionary history bivalves adapted to a wide range of physicochemical conditions, habitats, biological interactions, and feeding habits. Bivalves can have strikingly different size, and despite their apparently simple body plan, they evolved very different shell shapes, and complex anatomic structures. One of the most striking features of this class of animals is their peculiar mitochondrial biology: some bivalves have facultatively anaerobic mitochondria that allow them to survive prolonged periods of anoxia/hypoxia. Moreover, more than 100 species have now been reported showing the only known evolutionarily stable exception to the strictly maternal inheritance of mitochondria in animals, named doubly uniparental inheritance. Mitochondrial activity is fundamental to eukaryotic life, and thanks to their diversity and uncommon features, bivalves represent a great model system to expand our knowledge about mitochondrial biology, so far limited to a few species. We highlight recent works studying mitochondrial biology in bivalves at either genomic or physiological level. A link between these two approaches is still missing, and we believe that an integrated approach and collaborative relationships are the only possible ways to be successful in such endeavour.

Mitochondrial gene expression in single cells shape pancreatic beta cells' sub-populations and explain variation in insulin pathway

Mitochondrial gene expression is pivotal to cell metabolism. Nevertheless, it is unknown whether it diverges within a given cell type. Here, we analysed single-cell RNA-seq experiments from human pancreatic alpha (N = 3471) and beta cells (N = 1989), as well as mouse beta cells (N = 1094). Cluster analysis revealed two distinct human beta cells populations, which diverged by mitochondrial (mtDNA) and nuclear DNA (nDNA)-encoded oxidative phosphorylation (OXPHOS) gene expression in healthy and diabetic individuals, and in newborn but not in adult mice. Insulin gene expression was elevated in beta cells with higher mtDNA gene expression in humans and in young mice. Such human beta cell populations also diverged in mitochondrial RNA mutational repertoire, and in their selective signature, thus implying the existence of two previously overlooked distinct and conserved beta cell populations. While applying our approach to human alpha cells, two sub-populations of cells were identified which diverged in mtDNA gene expression, yet these cellular populations did not consistently diverge in nDNA OXPHOS genes expression, nor did they correlate with the expression of glucagon, the hallmark of alpha cells. Thus, pancreatic beta cells within an individual are divided into distinct groups with unique metabolic-mitochondrial signature.

Mitonuclear Interactions in the Maintenance of Mitochondrial Integrity

In eukaryotic cells, mitochondria originated in an α-proteobacterial endosymbiont. Although these organelles harbor their own genome, the large majority of genes, originally encoded in the endosymbiont, were either lost or transferred to the nucleus. As a consequence, mitochondria have become semi-autonomous and most of their processes require the import of nuclear-encoded components to be functional. Therefore, the mitochondrial-specific translation has evolved to be coordinated by mitonuclear interactions to respond to the energetic demands of the cell, acquiring unique and mosaic features. However, mitochondrial-DNA-encoded genes are essential for the assembly of the respiratory chain complexes. Impaired mitochondrial function due to oxidative damage and mutations has been associated with numerous human pathologies, the aging process, and cancer. In this review, we highlight the unique features of mitochondrial protein synthesis and provide a comprehensive insight into the mitonuclear crosstalk and its co-evolution, as well as the vulnerabilities of the animal mitochondrial genome.

A novel m.11406 T > A mutation in mitochondrial ND4 gene causes MELAS syndrome

Pathogenic point mutations of mitochondrial DNA (mtDNA) are associated with a large number of heterogeneous diseases involving multiple systems with which patients may present with a wide range of clinical phenotypes. In this study, we describe a novel heteroplasmic missense mutation, m.11406 T > A, of the ND4 gene encoding the subunit 4 of mitochondrial complex I in a 32-year-old woman with recurrent epileptic seizure, headache and bilateral hearing loss. Skeletal muscle histochemistry demonstrated that approximately 20% of fibers were cytochrome C oxidase (COX) deficient with increased activity of succinate dehydrogenase (SDH). Further investigations in muscle specimens showed significantly reduced level of ND4 protein. It is interesting that the subunits of complex I (ND1 and NDFUB8) and complex IV(CO1) were also remarkably decreased. These findings indicate that ND1, NDFUB8 and CO1 are more susceptible than other subunits to mutations in the mitochondrial ND4 gene.

Mitonuclear Compensatory Coevolution

In bilaterian animals, the mitochondrial genome is small, haploid, does not typically recombine, and is subject to accumulation of deleterious alleles via Muller's ratchet. These basic features of the genomic architecture present a paradox: mutational erosion of these genomes should lead to decline in mitochondrial function over time, yet no such decline is observed. Compensatory coevolution, whereby the nuclear genome evolves to compensate for the deleterious alleles in the mitochondrial genome, presents a potential solution to the paradox of Muller's ratchet without loss of function. Here, I review different proposed forms of mitonuclear compensatory coevolution. Empirical evidence from diverse eukaryotic taxa supports the mitonuclear compensatory coevolution hypothesis, but the ubiquity and importance of such compensatory coevolution remains a topic of debate.

Admixture Effects on Coevolved Metabolic Systems

Oxidative phosphorylation (OXPHOS) is the primary energy generating system in eukaryotic organisms. The complexes within the OXPHOS pathway are of mixed genomic origin. Although most subunit-coding genes are located within the nuclear genome, several genes are coded for in the mitochondrial genome. There is strong evidence to support coadaptation between the two genomes in these OXPHOS gene regions in order to create tight protein interactions necessary for a functional energetics system. In this study, we begin to assess the physiological impact of separating coevolved protein motifs that make up the highly conserved energy production pathway, as we hypothesize that divergent matings will significantly diminish the protein interactions and therefore hinder efficient OXPHOS activity We measured mitochondrial activity in high energy-demanding tissues from six strains of Mus musculus with varying degrees of mixed ancestral background. Mice with divergent mitochondrial and nuclear backgrounds consistently yielded lower mitochondrial activity. Bioinformatic analysis of common single nucleotide variants across the nuclear and mitochondrial genomes failed to identify any non-synonymous variants that could account for the energetic differences, suggesting that interpopulational mating between ancestrally distinct groups influences energy production efficiency.

Codon based co-occurrence network motifs in human mitochondria

The nucleotide polymorphism in the human mitochondrial genome (mtDNA) tolled by codon position bias plays an indispensable role in human population dispersion and expansion. Herein, genome-wide nucleotide co-occurrence networks were constructed using data comprised of five different geographical regions and around 3000 samples for each region. We developed a powerful network model to describe complex mitochondrial evolutionary patterns among codon and non-codon positions. We found evidence that the evolution of human mitochondria DNA is dominated by adaptive forces, particularly mutation and selection, which was supported by many previous studies. The diversity observed in the mtDNA was compared with mutations, co-occurring mutations, network motifs considering codon positions as causing agent. This comparison showed that long-range nucleotide co-occurrences have a large effect on genomic diversity. Most notably, codon motifs apparently underpinned the preferences among codon positions for co-evolution which is probably highly biased during the origin of the genetic code. Our analysis also showed that variable nucleotide positions of different human sub-populations implemented the independent mtDNA evolution to its geographical dispensation. Ergo, this study has provided both a network framework and a codon glance to investigate co-occurring genomic variations that are critical in underlying complex mitochondrial evolution.

Extracellular haemoglobin upregulates and binds to tissue factor on macrophages: Implications for coagulation and oxidative stress

The mechanisms of crosstalk between haemolysis, coagulation and innate immunity are evolutionarily conserved from the invertebrate haemocyanin to the vertebrate haemoglobin (Hb). In vertebrates, extracellular Hb resulting from haemolytic infections binds bacterial lipopolysaccharide (LPS) to unleash the antimicrobial redox activity of Hb. Because bacterial invasion also upregulates tissue factor (TF), the vertebrate coagulation initiator, we asked whether there may be functional interplay between the redox activity of Hb and the procoagulant activity of TF. Using real-time PCR, TF-specific ELISA, flow cytometry and TF activity assay, we found that Hb upregulated the expression of functional TF in macrophages. ELISA, flow cytometry and immunofluorescence microscopy showed binding between Hb and TF, in isolation and in situ. Bioinformatic analysis of Hb and TF protein sequences showed co-evolution across species, suggesting that Hbβ binds TF. Empirically, TF suppressed the LPS-induced activation of Hb redox activity. Furthermore, Hb desensitised TF to the effects of antioxidants like glutathione or serum. This bi-directional regulation between Hb and TF constitutes a novel link between coagulation and innate immunity. In addition, induction of TF by Hb is a potentially central mechanism for haemolysis to trigger coagulation.

Unraveling the complexity of mitochondrial complex I assembly: A dynamic process

Mammalian complex I is composed of 44 different subunits and its assembly requires at least 13 specific assembly factors. Proper function of the mitochondrial respiratory chain enzyme is of crucial importance for cell survival due to its major participation in energy production and cell signaling. Complex I assembly depends on the coordination of several crucial processes that need to be tightly interconnected and orchestrated by a number of assembly factors. The understanding of complex I assembly evolved from simple sequential concept to the more sophisticated modular assembly model describing a convoluted process. According to this model, the different modules assemble independently and associate afterwards with each other to form the final enzyme. In this review, we aim to unravel the complexity of complex I assembly and provide the latest insights in this fundamental and fascinating process. This article is part of a Special Issue entitled Respiratory complex I, edited by Volker Zickermann and Ulrich Brandt.

Loss of mitochondrial DNA-encoded protein ND1 results in disruption of complex I biogenesis during early stages of assembly

Mitochondrial complex I (NADH:ubiquinone oxidoreductase) must be assembled precisely from 45 protein subunits for it to function correctly. One of its mitochondrial DNA (mtDNA) encoded subunits, ND1, is incorporated during the early stages of complex I assembly. However, little is known about how mutations in ND1 affect this assembly process. We found that in human 143B cybrid cells carrying a homoplasmic MT-ND1 mutation, ND1 protein could not be translated. As a result, the early stages of complex I assembly were disrupted, with mature complex I undetectable and complex I-linked respiration severely reduced to 2.0% of control levels. Interestingly, complex IV (ferrocytochrome c:oxygen oxidoreductase) steady-state levels were also reduced to 40.3%, possibly due to its diminished stability in the absence of respiratory supercomplex formation. This was in comparison with 143B cybrid controls (that contained wild-type mtDNA on the same nuclear background), which exhibited normal complex I, complex IV, and supercomplex assembly. We conclude that the loss of ND1 stalls complex I assembly during the early stages of its biogenesis, which not only results in the loss of mature complex I but also disrupts the stability of complex IV and the respiratory supercomplex to cause mitochondrial dysfunction.-Lim, S. C., Hroudová, J., Van Bergen, N. J., Lopez Sanchez, M. I. G., Trounce, I. A., McKenzie, M. Loss of mitochondrial DNA-encoded protein ND1 results in disruption of complex I biogenesis during early stages of assembly.

Mitochondrial DNA variation in human radiation and disease

Environmental adaptation, predisposition to common diseases, and, potentially, speciation may all be linked through the adaptive potential of mitochondrial DNA (mtDNA) alterations of bioenergetics. This Perspective synthesizes evidence that human mtDNA variants may be adaptive or deleterious depending on environmental context and proposes that the accrual of mtDNA variation could contribute to animal speciation via adaptation to marginal environments.

Mitochondrial Involvement in Vertebrate Speciation? The Case of Mito-nuclear Genetic Divergence in Chameleons

Compatibility between the nuclear (nDNA) and mitochondrial (mtDNA) genomes is important for organismal health. However, its significance for major evolutionary processes such as speciation is unclear, especially in vertebrates. We previously identified a sharp mtDNA-specific sequence divergence between morphologically indistinguishable chameleon populations (Chamaeleo chamaeleon recticrista) across an ancient Israeli marine barrier (Jezreel Valley). Because mtDNA introgression and gender-based dispersal were ruled out, we hypothesized that mtDNA spatial division was maintained by mito-nuclear functional compensation. Here, we studied RNA-seq generated from each of ten chameleons representing the north and south populations and identified candidate nonsynonymous substitutions (NSSs) matching the mtDNA spatial distribution. The most prominent NSS occurred in 14 nDNA-encoded mitochondrial proteins. Increased chameleon sample size (N = 70) confirmed the geographic differentiation in POLRMT, NDUFA5, ACO1, LYRM4, MARS2, and ACAD9. Structural and functionality evaluation of these NSSs revealed high functionality. Mathematical modeling suggested that this mito-nuclear spatial divergence is consistent with hybrid breakdown. We conclude that our presented evidence and mathematical model underline mito-nuclear interactions as a likely role player in incipient speciation in vertebrates.

Fifty years of co-evolution and beyond: integrating co-evolution from molecules to species

Fifty years after Ehrlich and Raven's seminal paper, the idea of co-evolution continues to grow as a key concept in our understanding of organic evolution. This concept has not only provided a compelling synthesis between evolutionary biology and community ecology, but has also inspired research that extends beyond its original scope. In this article, we identify unresolved questions about the co-evolutionary process and advocate for the integration of co-evolutionary research from molecular to interspecific interactions. We address two basic questions: (i) What is co-evolution and how common is it? (ii) What is the unit of co-evolution? Both questions aim to explore the heart of the co-evolutionary process. Despite the claim that co-evolution is ubiquitous, we argue that there is in fact little evidence to support the view that reciprocal natural selection and coadaptation are common in nature. We also challenge the traditional view that co-evolution only occurs between traits of interacting species. Co-evolution has the potential to explain evolutionary processes and patterns that result from intra- and intermolecular biochemical interactions within cells, intergenomic interactions (e.g. nuclear-cytoplasmic) within species, as well as intergenomic interactions mediated by phenotypic traits between species. Research that bridges across these levels of organization will help to advance our understanding of the importance of the co-evolutionary processes in shaping the diversity of life on Earth.

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.

Coevolutionary analyses require phylogenetically deep alignments and better null models to accurately detect inter-protein contacts within and between species

Background

When biomolecules physically interact, natural selection operates on them jointly. Contacting positions in protein and RNA structures exhibit correlated patterns of sequence evolution due to constraints imposed by the interaction, and molecular arms races can develop between interacting proteins in pathogens and their hosts. To evaluate how well methods developed to detect coevolving residues within proteins can be adapted for cross-species, inter-protein analysis, we used statistical criteria to quantify the performance of these methods in detecting inter-protein residues within 8 angstroms of each other in the co-crystal structures of 33 bacterial protein interactions. We also evaluated their performance for detecting known residues at the interface of a host-virus protein complex with a partially solved structure.

Results

Our quantitative benchmarking showed that all coevolutionary methods clearly benefit from alignments with many sequences. Methods that aim to detect direct correlations generally outperform other approaches. However, faster mutual information based methods are occasionally competitive in small alignments and with relaxed false positive rates. Two commonly used null distributions are anti-conservative and have high false positive rates in some scenarios, although the empirical distribution of scores performs reasonably well with deep alignments.

Conclusions

We conclude that coevolutionary analysis of cross-species protein interactions holds great promise but requires sequencing many more species pairs.

Electronic supplementary material

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

Patterns of natural selection acting on the mitochondrial genome of a locally adapted fish species

Background

Mitochondrial DNA (mtDNA) is frequently used in population genetic studies and is usually considered as a neutral marker. However, given the functional importance of the proteins encoded by the mitochondrial genome, and the prominent role of mitochondria in cellular energy production, the assumption of neutrality is increasingly being questioned.

Results

We tested for evidence of selection on the mitochondrial genome of the Atlantic salmon, which is a locally adapted and widely farmed species and is distributed across a large latitudinal cline. We analysed 20 independent regions of the salmon mtDNA that represented nine genes (ND1, ND2, ND3, COX1, COX2, ATP6, ND4, ND5, and CYTB). These 20 mtDNA regions were sequenced using a 454 approach from samples collected across the entire European range of this species. We found evidence of positive selection at the ND1, ND3 and ND4 genes, which is supported by at least two different codon-based methods and also by differences in the chemical properties of the amino acids involved. The geographical distribution of some of the mutations indicated to be under selection was not random, and some mutations were private to artic populations. We discuss the possibility that selection acting on the Atlantic salmon mtDNA genome might be related to the need for increased metabolic efficiency at low temperatures.

Conclusions

The analysis of sequences representing nine mitochondrial genes that are involved in the OXPHOS pathway revealed signatures of positive selection in the mitochondrial genome of the Atlantic salmon. The properties of the amino acids involved suggest that some of the mutations that were identified to be under positive selection might have functional implications, possibly in relation to metabolic efficiency. Experimental evidence, and better understanding of regional phylogeographic structuring, are needed to clarify the potential role of selection acting on the mitochondrial genome of Atlantic salmon and other locally adapted fishes.

Electronic supplementary material

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

Differential Expression of Genes that Control Respiration Contribute to Thermal Adaptation in Redband Trout (Oncorhynchus mykiss gairdneri)

Organisms can adapt to local environmental conditions as a plastic response or become adapted through natural selection on genetic variation. The ability to adapt to increased water temperatures will be of paramount importance for many fish species as the climate continues to warm and water resources become limited. Because increased water temperatures will reduce the dissolved oxygen available for fish, we hypothesized that adaptation to low oxygen environments would involve improved respiration through oxidative phosphorylation (OXPHOS). To test this hypothesis, we subjected individuals from two ecologically divergent populations of inland (redband) rainbow trout (Oncorhynchus mykiss gairdneri) with historically different temperature regimes (desert and montane) and their F1 progeny to diel cycles of temperature stress and then examined gene expression data for 80 nuclear- and mitochondrial-encoded OXPHOS subunits that participate in respiration. Of the 80 transcripts, 7 showed ≥ 2-fold difference in expression levels in gill tissue from desert fish under heat stress whereas the montane fish had none and the F1 only had one differentially expressed gene. A structural analysis of the proteins encoded by those genes suggests that the response could coordinate the formation of supercomplexes and oligomers. Supercomplexes may increase the efficiency of respiration because complexes I, III, and IV are brought into close proximity and oligomerization of complex V alters the macrostructure of mitochondria to improve respiration. Significant differences in gene expression patterns in response to heat stress in a common environment indicate that the response was not due to plasticity but had a genetic basis.

Redefining the roles of mitochondrial DNA-encoded subunits in respiratory Complex I assembly

Respiratory Complex I deficiency is implicated in numerous degenerative and metabolic diseases. In particular, mutations in several mitochondrial DNA (mtDNA)-encoded Complex I subunits including ND4, ND5 and ND6 have been identified in several neurological diseases. We previously demonstrated that these subunits played essential roles in Complex I assembly which in turn affected mitochondrial function. Here, we carried out a comprehensive study of the Complex I assembly pathway. We identified a new Complex I intermediate containing both membrane and matrix arms at an early assembly stage. We find that lack of the ND6 subunit does not hinder membrane arm formation; instead it recruits ND1 and ND5 enters the intermediate. While ND4 is important for the formation of the newly identified intermediate, the addition of ND5 stabilizes the complex and is required for the critical transition from Complex I to supercomplex assembly. As a result, the Complex I assembly pathway has been redefined in this study.

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).

Disrupting mitochondrial-nuclear coevolution affects OXPHOS complex I integrity and impacts human health

The mutation rate of the mitochondrial DNA (mtDNA), which is higher by an order of magnitude as compared with the nuclear genome, enforces tight mitonuclear coevolution to maintain mitochondrial activities. Interruption of such coevolution plays a role in interpopulation hybrid breakdown, speciation events, and disease susceptibility. Previously, we found an elevated amino acid replacement rate and positive selection in the nuclear DNA-encoded oxidative phosphorylation (OXPHOS) complex I subunit NDUFC2, a phenomenon important for the direct interaction of NDUFC2 with the mtDNA-encoded complex I subunit ND4. This finding underlines the importance of mitonuclear coevolution to physical interactions between mtDNA and nuclear DNA-encoded factors. Nevertheless, it remains unclear whether this interaction is important for the stability and activity of complex I. Here, we show that siRNA silencing of NDUFC2 reduced growth of human D-407 retinal pigment epithelial cells, significantly diminished mitochondrial membrane potential, and interfered with complex I integrity. Moreover, site-directed mutagenesis of a positively selected amino acid in NDUFC2 significantly interfered with the interaction of NDUFC2 with its mtDNA-encoded partner ND4. Finally, we show that a genotype combination involving this amino acid (NDUFC2 residue 46) and the mtDNA haplogroup HV likely altered susceptibility to type 2 diabetes mellitus in Ashkenazi Jews. Therefore, mitonuclear coevolution is important for maintaining mitonuclear factor interactions, OXPHOS, and for human health.

An update on complex I assembly: the assembly of players

Defects in Complex I assembly is one of the emerging underlying causes of severe mitochondrial disorders. The assembly of Complex I has been difficult to understand due to its large size, dual genetic control and the number of proteins involved. Mutations in Complex I subunits as well as assembly factors have been reported to hinder its assembly and give rise to a range of mitochondria disorders. In this review, we summarize the recent progress made in understanding the Complex I assembly pathway. In particularly, we focus on the known as well as novel assembly factors and their role in assembly of Complex I and human disease.

Bioenergetic constraints on the evolution of complex life

All morphologically complex life on Earth, beyond the level of cyanobacteria, is eukaryotic. All eukaryotes share a common ancestor that was already a complex cell. Despite their biochemical virtuosity, prokaryotes show little tendency to evolve eukaryotic traits or large genomes. Here I argue that prokaryotes are constrained by their membrane bioenergetics, for fundamental reasons relating to the origin of life. Eukaryotes arose in a rare endosymbiosis between two prokaryotes, which broke the energetic constraints on prokaryotes and gave rise to mitochondria. Loss of almost all mitochondrial genes produced an extreme genomic asymmetry, in which tiny mitochondrial genomes support, energetically, a massive nuclear genome, giving eukaryotes three to five orders of magnitude more energy per gene than prokaryotes. The requirement for endosymbiosis radically altered selection on eukaryotes, potentially explaining the evolution of unique traits, including the nucleus, sex, two sexes, speciation, and aging.

Practical aspects of protein co-evolution

Co-evolution is a fundamental aspect of Evolutionary Theory. At the molecular level, co-evolutionary linkages between protein families have been used as indicators of protein interactions and functional relationships from long ago. Due to the complexity of the problem and the amount of genomic data required for these approaches to achieve good performances, it took a relatively long time from the appearance of the first ideas and concepts to the quotidian application of these approaches and their incorporation to the standard toolboxes of bioinformaticians and molecular biologists. Today, these methodologies are mature (both in terms of performance and usability/implementation), and the genomic information that feeds them large enough to allow their general application. This review tries to summarize the current landscape of co-evolution-based methodologies, with a strong emphasis on describing interesting cases where their application to important biological systems, alone or in combination with other computational and experimental approaches, allowed getting new insight into these.

Functional recurrent mutations in the human mitochondrial phylogeny: dual roles in evolution and disease

Mutations frequently reoccur in the human mitochondrial DNA (mtDNA). However, it is unclear whether recurrent mtDNA nodal mutations (RNMs), that is, recurrent mutations in stems of unrelated phylogenetic nodes, are functional and hence selectively constrained. To answer this question, we performed comprehensive parsimony and maximum likelihood analyses of 9,868 publicly available whole human mtDNAs revealing 1,606 single nodal mutations (SNMs) and 679 RNMs. We then evaluated the potential functionality of synonymous, nonsynonymous and RNA SNMs and RNMs. For synonymous mutations, we have implemented the Codon Adaptation Index. For nonsynonymous mutations, we assessed evolutionary conservation, and employed previously described pathogenicity score assessment tools. For RNA genes’ mutations, we designed a bioinformatic tool which compiled evolutionary conservation and potential effect on RNA structure. While comparing the functionality scores of nonsynonymous and RNA SNMs and RNMs with those of disease-causing mtDNA mutations, we found significant difference (P < 0.001). However, 24 RNMs and 67 SNMs had comparable values with disease-causing mutations reflecting their potential function thus being the best candidates to participate in adaptive events of unrelated lineages. Strikingly, some functional RNMs occurred in unrelated mtDNA lineages that independently altered susceptibility to the same diseases, thus suggesting common functionality. To our knowledge, this is the most comprehensive analysis of selective signatures in the mtDNA not only within proteins but also within RNA genes. For the first time, we discover virtually all positively selected RNMs in our phylogeny while emphasizing their dual role in past evolutionary events and in disease today.

Emerging methods in protein co-evolution

Co-evolution is a fundamental component of the theory of evolution and is essential for understanding the relationships between species in complex ecological networks. A wide range of co-evolution-inspired computational methods has been designed to predict molecular interactions, but it is only recently that important advances have been made. Breakthroughs in the handling of phylogenetic information and in disentangling indirect relationships have resulted in an improved capacity to predict interactions between proteins and contacts between different protein residues. Here, we review the main co-evolution-based computational approaches, their theoretical basis, potential applications and foreseeable developments.

Protein co-evolution: how do we combine bioinformatics and experimental approaches?

Molecular co-evolution is manifested by compensatory changes in proteins designed to enable adaptation to their natural environment. In recent years, bioinformatics approaches allowed for the detection of co-evolution at the level of the whole protein or of specific residues. Such efforts enabled prediction of protein-protein interactions, functional assignments of proteins and the identification of interacting residues, thereby providing information on protein structure. Still, despite such advances, relatively little is known regarding the functional implications of sequence divergence resulting from protein co-evolution. While bioinformatics approaches usually analyze thousands of proteins to obtain a broad view of protein co-evolution, experimental evaluation of protein co-evolution serves to study only individual proteins. In this review, we describe recent advances in bioinformatics and experimental efforts aimed at examining protein co-evolution. Accordingly, we discuss possible modes of crosstalk between the bioinformatics and experimental approaches to facilitate the identification of co-evolutionary signals in proteins and to understand their implications for the structure and function of proteins.

Association of genetic variation in the mitochondrial genome with blood pressure and metabolic traits

Elevated blood pressure (BP) is a major risk factor for cardiovascular disease. Several studies have noted a consistent maternal effect on BP; consequently, mitochondrial DNA variation has become an additional target of investigation of the missing BP heritability. Analyses of common mitochondrial DNA polymorphisms, however, have not found evidence of association with hypertension. To explore associations of uncommon (frequency>5%) mitochon drial DNA variants with BP, we identified uncommon/rare variants through sequencing the entire mitochondrial genome in 32 unrelated individuals with extreme-high BP in the Framingham Heart Study and genotyped 40 mitochondrial single nucleotide polymorphisms in 7219 Framingham Heart Study participants. The nonsynonymous mitochondrial single nucleotide polymorphism 5913G>A (Asp4Asn) in the cytochrome c oxidase subunit 1 of respiratory complex IV demonstrated significant associations with BP and fasting blood glucose (FBG) levels. Individuals with the rare 5913A allele had, on average, 7-mm Hg higher systolic BP at baseline (Pempirical=0.05) and 17-mg/dL higher mean FBG over 25 years of follow-up (Pempirical=0.009). Significant associations with FBG levels were also detected for nonsynonymous mitochondrial single nucleotide polymorphism 3316G>A (Ala4Thr) in the NADH dehydrogenase subunit 1 of complex I. On average, individuals with rare allele 3316A had 17- and 25-mg/dL higher FBG at baseline (Pempirical=0.01) and over 25 years of follow-up (Pempirical=0.007). Our findings provide the first evidence of putative association of variants in the mitochondrial genome with systolic BP and FBG in the general population. Replication in independent samples, however, is needed to confirm these putative associations.

A three-dimensional topology of complex I inferred from evolutionary correlations

Background

The quaternary structure of eukaryotic NADH:ubiquinone oxidoreductase (complex I), the largest complex of the oxidative phosphorylation, is still mostly unresolved. Furthermore, it is unknown where transiently bound assembly factors interact with complex I. We therefore asked whether the evolution of complex I contains information about its 3D topology and the binding positions of its assembly factors. We approached these questions by correlating the evolutionary rates of eukaryotic complex I subunits using the mirror-tree method and mapping the results into a 3D representation by multidimensional scaling.

Results

More than 60% of the evolutionary correlation among the conserved seven subunits of the complex I matrix arm can be explained by the physical distance between the subunits. The three-dimensional evolutionary model of the eukaryotic conserved matrix arm has a striking similarity to the matrix arm quaternary structure in the bacterium Thermus thermophilus (rmsd=19 Å) and supports the previous finding that in eukaryotes the N-module is turned relative to the Q-module when compared to bacteria. By contrast, the evolutionary rates contained little information about the structure of the membrane arm. A large evolutionary model of 45 subunits and assembly factors allows to predict subunit positions and interactions (rmsd = 52.6 Å). The model supports an interaction of NDUFAF3, C8orf38 and C2orf56 during the assembly of the proximal matrix arm and the membrane arm. The model further suggests a tight relationship between the assembly factor NUBPL and NDUFA2, which both have been linked to iron-sulfur cluster assembly, as well as between NDUFA12 and its paralog, the assembly factor NDUFAF2.

Conclusions

The physical distance between subunits of complex I is a major correlate of the rate of protein evolution in the complex I matrix arm and is sufficient to infer parts of the complex’s structure with high accuracy. The resulting evolutionary model predicts the positions of a number of subunits and assembly factors.

Mitochondrial DNA variation, but not nuclear DNA, sharply divides morphologically identical chameleons along an ancient geographic barrier

The Levant is an important migration bridge, harboring border-zones between Afrotropical and palearctic species. Accordingly, Chameleo chameleon, a common species throughout the Mediterranean basin, is morphologically divided in the southern Levant (Israel) into two subspecies, Chamaeleo chamaeleon recticrista (CCR) and C. c. musae (CCM). CCR mostly inhabits the Mediterranean climate (northern Israel), while CCM inhabits the sands of the north-western Negev Desert (southern Israel). AFLP analysis of 94 geographically well dispersed specimens indicated moderate genetic differentiation (PhiPT = 0.097), consistent with the classical division into the two subspecies, CCR and CCM. In contrast, sequence analysis of a 637 bp coding mitochondrial DNA (mtDNA) fragment revealed two distinct phylogenetic clusters which were not consistent with the morphological division: one mtDNA cluster consisted of CCR specimens collected in regions northern of the Jezreel Valley and another mtDNA cluster harboring specimens pertaining to both the CCR and CCM subspecies but collected southern of the Jezreel Valley. AMOVA indicated clear mtDNA differentiation between specimens collected northern and southern to the Jezreel Valley (PhiPT = 0.79), which was further supported by a very low coalescent-based estimate of effective migration rates. Whole chameleon mtDNA sequencing (∼17,400 bp) generated from 11 well dispersed geographic locations revealed 325 mutations sharply differentiating the two mtDNA clusters, suggesting a long allopatric history further supported by BEAST. This separation correlated temporally with the existence of an at least 1 million year old marine barrier at the Jezreel Valley exactly where the mtDNA clusters meet. We discuss possible involvement of gender-dependent life history differences in maintaining such mtDNA genetic differentiation and suggest that it reflects (ancient) local adaptation to mitochondrial-related traits.

Mitochondrial-nuclear co-evolution and its effects on OXPHOS activity and regulation

Factors required for mitochondrial function are encoded both by the nuclear and mitochondrial genomes. The order of magnitude higher mutation rate of animal mitochondrial DNA (mtDNA) enforces tight co-evolution of mtDNA and nuclear DNA encoded factors. In this essay we argue that such co evolution exists at the population and inter-specific levels and affect disease susceptibility. We also argue for the existence of three modes of co-evolution in the mitochondrial genetic system, which include the interaction of mtDNA and nuclear DNA encoded proteins, nuclear protein - mtDNA-encoded RNA interaction within the mitochondrial translation machinery and nuclear DNA encoded proteins-mtDNA binging sites interaction in the frame of the mtDNA replication and transcription machineries. These modes of co evolution require co-regulation of the interacting factors encoded by the two genomes. Thus co evolution plays an important role in modulating mitochondrial activity. This article is part of a Special Issue entitled: Mitochondrial Gene Expression.

Mitonuclear match: optimizing fitness and fertility over generations drives ageing within generations

Many conserved eukaryotic traits, including apoptosis, two sexes, speciation and ageing, can be causally linked to a bioenergetic requirement for mitochondrial genes. Mitochondrial genes encode proteins involved in cell respiration, which interact closely with proteins encoded by nuclear genes. Functional respiration requires the coadaptation of mitochondrial and nuclear genes, despite divergent tempi and modes of evolution. Free-radical signals emerge directly from the biophysics of mosaic respiratory chains encoded by two genomes prone to mismatch, with apoptosis being the default penalty for compromised respiration. Selection for genomic matching is facilitated by two sexes, and optimizes fitness, adaptability and fertility in youth. Mismatches cause infertility, low fitness, hybrid breakdown, and potentially speciation. The dynamics of selection for mitonuclear function optimize fitness over generations, but the same selective processes also operate within generations, driving ageing and age-related diseases. This coherent view of eukaryotic energetics offers striking insights into infertility and age-related diseases.

Understanding mitochondrial complex I assembly in health and disease

Complex I (NADH:ubiquinone oxidoreductase) is the largest multimeric enzyme complex of the mitochondrial respiratory chain, which is responsible for electron transport and the generation of a proton gradient across the mitochondrial inner membrane to drive ATP production. Eukaryotic complex I consists of 14 conserved subunits, which are homologous to the bacterial subunits, and more than 26 accessory subunits. In mammals, complex I consists of 45 subunits, which must be assembled correctly to form the properly functioning mature complex. Complex I dysfunction is the most common oxidative phosphorylation (OXPHOS) disorder in humans and defects in the complex I assembly process are often observed. This assembly process has been difficult to characterize because of its large size, the lack of a high resolution structure for complex I, and its dual control by nuclear and mitochondrial DNA. However, in recent years, some of the atomic structure of the complex has been resolved and new insights into complex I assembly have been generated. Furthermore, a number of proteins have been identified as assembly factors for complex I biogenesis and many patients carrying mutations in genes associated with complex I deficiency and mitochondrial diseases have been discovered. Here, we review the current knowledge of the eukaryotic complex I assembly process and new insights from the identification of novel assembly factors. This article is part of a Special Issue entitled: Biogenesis/Assembly of Respiratory Enzyme Complexes.

Molecular evolution of protein conformational changes revealed by a network of evolutionarily coupled residues

An improved understanding of protein conformational changes has broad implications for elucidating the mechanisms of various biological processes and for the design of protein engineering experiments. Understanding rearrangements of residue interactions is a key component in the challenge of describing structural transitions. Evolutionary properties of protein sequences and structures are extensively studied; however, evolution of protein motions, especially with respect to interaction rearrangements, has yet to be explored. Here, we investigated the relationship between sequence evolution and protein conformational changes and discovered that structural transitions are encoded in amino acid sequences as coevolving residue pairs. Furthermore, we found that highly coevolving residues are clustered in the flexible regions of proteins and facilitate structural transitions by forming and disrupting their interactions cooperatively. Our results provide insight into the evolution of protein conformational changes and help to identify residues important for structural transitions.