Adaptive selection of mitochondrial complex I subunits during primate radiation

Adverse effect propensity: A new feature of Gulf War illness predicted by environmental exposures

A third of 1990-1 Gulf-deployed personnel developed drug/chemical-induced multisymptom illness, "Gulf War illness" (GWI). Veterans with GWI (VGWI) report increased drug/exposure adverse effects (AEs). Using previously collected data from a case-control study, we evaluated whether the fraction of exposures that engendered AEs ("AE Propensity") is increased in VGWI (it was); whether AE Propensity is related to self-rated "chemical sensitivity" (it did); and whether specific exposures "predicted" AE Propensity (they did). Pesticides and radiation exposure were significant predictors, with copper significantly "protective"-in the total sample (adjusted for GWI-status) and separately in VGWI and controls, on multivariable regression. Mitochondrial impairment and oxidative stress (OS) underlie AEs from many exposures irrespective of nominal specific mechanism. We hypothesize that mitochondrial toxicity and interrelated OS from pesticides and radiation position people on the steep part of the curve of mitochondrial impairment and OS versus symptom/biological disruption, amplifying impact of new exposures. Copper, meanwhile, is involved in critical OS detoxification processes.

Inefficient Batteries in Heart Failure: Metabolic Bottlenecks Disrupting the Mitochondrial Ecosystem

Mitochondrial abnormalities have long been described in the setting of cardiomyopathies and heart failure (HF), yet the mechanisms of mitochondrial dysfunction in cardiac pathophysiology remain poorly understood. Many studies have described HF as an energy-deprived state characterized by a decline in adenosine triphosphate production, largely driven by impaired oxidative phosphorylation. However, impairments in oxidative phosphorylation extend beyond a simple decline in adenosine triphosphate production and, in fact, reflect pervasive metabolic aberrations that cannot be fully appreciated from the isolated, often siloed, interrogation of individual aspects of mitochondrial function. With the application of broader and deeper examinations into mitochondrial and metabolic systems, recent data suggest that HF with preserved ejection fraction is likely metabolically disparate from HF with reduced ejection fraction. In our review, we introduce the concept of the mitochondrial ecosystem, comprising intricate systems of metabolic pathways and dynamic changes in mitochondrial networks and subcellular locations. The mitochondrial ecosystem exists in a delicate balance, and perturbations in one component often have a ripple effect, influencing both upstream and downstream cellular pathways with effects enhanced by mitochondrial genetic variation. Expanding and deepening our vantage of the mitochondrial ecosystem in HF is critical to identifying consistent metabolic perturbations to develop therapeutics aimed at preventing and improving outcomes in HF.

A Naturally Heteroplasmic Clam Provides Clues about the Effects of Genetic Bottleneck on Paternal mtDNA

Mitochondrial DNA (mtDNA) is present in multiple copies within an organism. Since these copies are not identical, a single individual carries a heterogeneous population of mtDNAs, a condition known as heteroplasmy. Several factors play a role in the dynamics of the within-organism mtDNA population: among them, genetic bottlenecks, selection, and strictly maternal inheritance are known to shape the levels of heteroplasmy across mtDNAs.

In Metazoa, the only evolutionarily stable exception to the strictly maternal inheritance of mitochondria is the doubly uniparental inheritance (DUI), reported in 100+ bivalve species. In DUI species, there are two highly divergent mtDNA lineages, one inherited through oocyte mitochondria (F-type) and the other through sperm mitochondria (M-type). Having both parents contributing to the mtDNA pool of the progeny makes DUI a unique system to study the dynamics of mtDNA populations. Since, in bivalves, the spermatozoon has few mitochondria (4–5), M-type mtDNA faces a tight bottleneck during embryo segregation, one of the narrowest mitochondrial bottlenecks investigated so far.

Here, we analyzed the F- and M-type mtDNA variability within individuals of the DUI species Ruditapes philippinarum and investigated for the first time the effects of such a narrow bottleneck affecting mtDNA populations. As a potential consequence of this narrow bottleneck, the M-type mtDNA shows a large variability in different tissues, a condition so pronounced that it leads to genotypes from different tissues of the same individual not to cluster together. We believe that such results may help understanding the effect of low population size on mtDNA bottleneck.

MitImpact 3: modeling the residue interaction network of the Respiratory Chain subunits

Numerous lines of evidence have shown that the interaction between the nuclear and mitochondrial genomes ensures the efficient functioning of the OXPHOS complexes, with substantial implications in bioenergetics, adaptation, and disease. Their interaction is a fascinating and complex trait of the eukaryotic cell that MitImpact explores with its third major release. MitImpact expands its collection of genomic, clinical, and functional annotations of all non-synonymous substitutions of the human mitochondrial genome with new information on putative Compensated Pathogenic Deviations and co-varying amino acid sites of the Respiratory Chain subunits. It further provides evidence of energetic and structural residue compensation by techniques of molecular dynamics simulation. MitImpact is freely accessible at http://mitimpact.css-mendel.it.

Signatures of host-pathogen evolutionary conflict reveal MISTR-A conserved MItochondrial STress Response network

Host-pathogen conflicts leave genetic signatures in genes that are critical for host defense functions. Using these "molecular scars" as a guide to discover gene functions, we discovered a vertebrate-specific MItochondrial STress Response (MISTR) circuit. MISTR proteins are associated with electron transport chain (ETC) factors and activated by stress signals such as interferon gamma (IFNγ) and hypoxia. Upon stress, ultraconserved microRNAs (miRNAs) down-regulate MISTR1(NDUFA4) followed by replacement with paralogs MItochondrial STress Response AntiViral (MISTRAV) and/or MItochondrial STress Response Hypoxia (MISTRH). While cells lacking MISTR1(NDUFA4) are more sensitive to chemical and viral apoptotic triggers, cells lacking MISTRAV or expressing the squirrelpox virus-encoded vMISTRAV exhibit resistance to the same insults. Rapid evolution signatures across primate genomes for MISTR1(NDUFA4) and MISTRAV indicate recent and ongoing conflicts with pathogens. MISTR homologs are also found in plants, yeasts, a fish virus, and an algal virus indicating ancient origins and suggesting diverse means of altering mitochondrial function under stress. The discovery of MISTR circuitry highlights the use of evolution-guided studies to reveal fundamental biological processes.

Complex IV - The regulatory center of mitochondrial oxidative phosphorylation

ATP, the universal energy currency in all living cells, is mainly synthesized in mitochondria by oxidative phosphorylation (OXPHOS). The final and rate limiting step of the respiratory chain is cytochrome c oxidase (COX) which represents the regulatory center of OXPHOS. COX is regulated through binding of various effectors to its "supernumerary" subunits, by reversible phosphorylation, and by expression of subunit isoforms. Of particular interest is its feedback inhibition by ATP, the final product of OXPHOS. This "allosteric ATP-inhibition" of phosphorylated and dimeric COX maintains a low and healthy mitochondrial membrane potential (relaxed state), and prevents the formation of ROS (reactive oxygen species) which are known to cause numerous diseases. Excessive work and stress abolish this feedback inhibition of COX by Ca²⁺-activated dephosphorylation which leads to monomerization and movement of NDUFA4 from complex I to COX with higher rates of COX activity and ATP synthesis (active state) but increased ROS formation and decreased efficiency.

Regulation of cytochrome c oxidase contributes to health and optimal life

The generation of cellular energy in the form of ATP occurs mainly in mitochondria by oxidative phosphorylation. Cytochrome c oxidase (CytOx), the oxygen accepting and rate-limiting step of the respiratory chain, regulates the supply of variable ATP demands in cells by “allosteric ATP-inhibition of CytOx.” This mechanism is based on inhibition of oxygen uptake of CytOx at high ATP/ADP ratios and low ferrocytochrome c concentrations in the mitochondrial matrix via cooperative interaction of the two substrate binding sites in dimeric CytOx. The mechanism keeps mitochondrial membrane potential ΔΨ_m and reactive oxygen species (ROS) formation at low healthy values. Stress signals increase cytosolic calcium leading to Ca²⁺-dependent dephosphorylation of CytOx subunit I at the cytosolic side accompanied by switching off the allosteric ATP-inhibition and monomerization of CytOx. This is followed by increase of ΔΨ_m and formation of ROS. A hypothesis is presented suggesting a dynamic change of binding of NDUFA4, originally identified as a subunit of complex I, between monomeric CytOx (active state with high ΔΨ_m, high ROS and low efficiency) and complex I (resting state with low ΔΨ_m, low ROS and high efficiency).

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.

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.

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.

A transgenic mouse expressing miR-210 in proximal tubule cells shows mitochondrial alteration: possible association of miR-210 with a shift in energy metabolism

Previously we reported that the microRNA miR-210 is aberrantly upregulated in clear cell renal cell carcinoma (ccRCC) via deregulation of the VHL-HIF pathway. In the present study, to investigate the biological impact of miR-210 in ccRCC tumorigenesis, we developed a transgenic mouse line expressing miR-210 in proximal tubule cells under control of the mouse SGLT2/Slc5a2 promoter. Light microscopy revealed desquamation of the tubule cells and regeneration of the proximal tubule, suggesting that miR-210 expression led to damage of the proximal tubule cells. Electron microscopy revealed alterations to the mitochondria in proximal tubule cells, with marked reduction of the mitochondrial inner membrane, which is the main site of ATP production via oxidative phosphorylation (OxPhos). An additional in vitro study revealed that this loss of the inner membrane was associated with downregulation of Iscu and Ndufa4, the target genes of miR-210, suggesting that the miR-210-ISCU/NDUFA4 axis may affect mitochondrial energy metabolism. Furthermore, metabolome analysis revealed activation of anaerobic glycolysis in miR-210-transfected cells, and consistent with this the secretion of lactate, the final metabolite of anaerobic glycolysis, was significantly increased. Lactate concentration was higher in the kidney cortex of transgenic mice relative to wild-type mice, although the difference was not significant (p = 0.070). On the basis of these findings, we propose that miR-210 may induce a shift of energy metabolism from OxPhos to glycolysis by acting on the mitochondrial inner membrane. In addition to activation of glycolysis, we observed activation of the pentose phosphate pathway (PPP) and an increase in the total amount of amino acids in miR-210-transfected cells. This may help cells synthesize nucleotides and proteins for building new cells. These results suggest that miR-210 may be involved in the metabolic changes in the early stage of ccRCC development, helping the cancer cells to acquire growth and survival advantages. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Evolutionary Rate Correlation between Mitochondrial-Encoded and Mitochondria-Associated Nuclear-Encoded Proteins in Insects

The mitochondrion is a pivotal organelle for energy production, and includes components encoded by both the mitochondrial and nuclear genomes. Functional and evolutionary interactions are expected between the nuclear- and mitochondrial-encoded components. The topic is of broad interest in biology, with implications to genetics, evolution, and medicine. Here, we compare the evolutionary rates of mitochondrial proteins and ribosomal RNAs to rates of mitochondria-associated nuclear-encoded proteins, across the major orders of holometabolous insects. There are significant evolutionary rate correlations (ERCs) between mitochondrial-encoded and mitochondria-associated nuclear-encoded proteins, which are likely driven by different rates of mitochondrial sequence evolution and correlated changes in the interacting nuclear-encoded proteins. The pattern holds after correction for phylogenetic relationships and considering protein conservation levels. Correlations are stronger for both nuclear-encoded OXPHOS proteins that are in contact with mitochondrial OXPHOS proteins and for nuclear-encoded mitochondrial ribosomal amino acids directly contacting the mitochondrial rRNAs. We find that ERC between mitochondrial- and nuclear-encoded proteins is a strong predictor of nuclear-encoded proteins known to interact with mitochondria, and ERC shows promise for identifying new candidate proteins with mitochondrial function. Twenty-three additional candidate nuclear-encoded proteins warrant further study for mitochondrial function based on this approach, including proteins in the minichromosome maintenance helicase complex.

HIF-1α Stabilization Increases miR-210 Eliciting First Trimester Extravillous Trophoblast Mitochondrial Dysfunction

Preeclampsia is associated with first trimester placental dysfunction. miR-210, a small non-coding RNA, is increased in the preeclamptic placenta. The effects of elevated miR-210 on placental function remain unclear. The objectives of this study were to identify targets of miR-210 in first trimester primary extravillous trophoblasts (EVTs) and to investigate functional pathways altered by elevated placental miR-210 during early pregnancy. EVTs isolated from first trimester placentas were exposed to cobalt chloride (CoCl₂), a HIF-1α stabilizer and hypoxia mimetic, and miR-210 expression by qPCR, HIF1α protein levels by western blot and cell invasion were assessed. A custom TruSeq RNA array, including all known/predicted miR-210 targets, was run using miR-210 and miR-negative control transfected EVTs. Mitochondrial function was assessed by high resolution respirometry in transfected EVTs. EVTs exposed to CoCl₂ showed a dose and time-dependent increase in miR-210 and HIF1α and reductions in cell invasion. The TruSeq array identified 49 altered genes in miR-210 transfected EVTs with 27 genes repressed and 22 enhanced. Three of the top six significantly repressed genes, NDUFA4, SDHD, and ISCU, are associated with mitochondrial function. miR-210 transfected EVTs had decreased maximal, complex II and complex I+II mitochondrial respiration. This study suggests that miR-210 alters first trimester trophoblast function. miR-210 overexpression alters EVT mitochondrial function in early pregnancy. Mitochondrial dysfunction may lead to increased reactive oxygen species, trophoblast cell damage and likely contributes to the pathogenesis of preeclampsia.

A novel mutation in the NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 4 (Ndufa4) gene links mitochondrial dysfunction to the development of diabetes in a rodent model

The mechanisms underlying diabetes remain unresolved. The Cohen diabetic rat represents a model of diet-induced diabetes, in which the disease is induced after exposure to a diabetogenic diet (DD) in the diabetes-sensitive (CDs/y) but not in the -resistant (CDr/y) strain. Diet imposes a metabolic strain that leads to diabetes in the appropriate genetic background. We previously identified, through whole-genome linkage analysis, a diabetes-related quantitative trait locus on rat chromosome 4 (RNO4), which incorporates NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 4 (Ndufa4), a nuclear gene that affects mitochondrial function. Here, we sequenced the gene and found a major deletion in CDs/y that leads to lack of expression of the NDUFA4 protein that has been reported to be involved in the activities of mitochondrial complexes I and IV. In the absence of NDUFA4 in the diabetic CDs/y on DD, complex I activity is reduced in comparison to that in nondiabetic CDs/y on regular diet and CDr/y on either diet; complex IV activity is reduced in both strains provided DD, and thus as a result of diet and unrelated to the gene mutation. ATP fails to increase in diabetic CDs/y in response to DD, in comparison to nondiabetic CDr/y on DD. Plasma malondialdehyde levels are elevated in CDs/y on DD, whereas SOD1 and SOD2 levels fail to increase, indicating increased oxidative stress and inability of the pancreas to generate an appropriate antioxidative stress response. These findings suggest that the Ndufa4 mutation in CDs/y on DD is directly associated with mitochondrial dysfunction, which we attribute to the lack of expression of NDUFA4 and to diet, and which prevents the anticipated increase in ATP production. The resulting enhanced oxidative stress impairs the ability of the pancreas to secrete insulin, leading to the development of diabetes. This is the first demonstration of an inherited mutation in a nuclear gene that adversely affects mitochondrial function and promotes diet-induced diabetes.

Summary: Here, we report, for the first time, a major inherited mutation in a diabetes-prone animal model that adversely affects mitochondrial function and leads, through oxidative stress, to the development of diet-induced diabetes.

The mitochondrial genomes of two walnut pests, Gastrolina depressa depressa and G. depressa thoracica (Coleoptera: Chrysomelidae), and phylogenetic analyses

In this study, the mitochondrial genomes (mitogenomes) of two walnut leaf insect pests, Gastrolina depressa depressa and G. depressa thoracica, were sequenced by Sanger sequencing technology. The mitogenome of G. depressa thoracica was complete at 16,109 bp in length, while the mitogenome of G. depressa depressa (14,277 bp) was partial. The genomic analyses indicated that both mitogenomes have the typical gene content and arrangement. The formerly identified elements, 'TAGTA' between trnSer(UCN) and nad2, and 'ATGATAA' between atp8 and atp6, were more conserved than that between nad4L and nad4, which was 'ATGTTAA' in Coleoptera excluding Polyphaga. Phylogenetic analyses of the 13 protein-coding genes from 36 coleopteran species well supported a close affinity between the subfamily Chrysomelinae including G. depressa thoracica and G. depressa depressa and Galerucinae, as well as a sister relationship of ((Eumolpinae + Cryptocephalinae) + Cassidinae) within Chrysomelidae.

Punctuated invasion of water, ice, snow and terrestrial ecozones by segmented worms (Oligochaeta: Enchytraeidae: Mesenchytraeus)

Segmented worms (Annelida) are among the most successful animal inhabitants of extreme environments worldwide. An unusual group of enchytraeid oligochaetes of genus Mesenchytraeus are abundant in the Pacific northwestern region of North America and occupy geographically proximal ecozones ranging from low elevation rainforests and waterways to high altitude glaciers. Along this altitudinal transect, Mesenchytraeus representatives from disparate habitat types were collected and subjected to deep mitochondrial and nuclear phylogenetic analyses. Our data identify significant topological discordance among gene trees, and near equivalent interspecific divergence levels indicative of a rapid radiation event. Collectively, our results identify a Mesenchytraeus 'explosion' coincident with mountain building in the Pacific northwestern region that gave rise to closely related aquatic, ice, snow and terrestrial worms.

Evolutionary analysis of mitochondrially encoded proteins of toad-headed lizards, Phrynocephalus, along an altitudinal gradient

Background

Animals living at high altitude must adapt to environments with hypoxia and low temperatures, but relatively little is known about underlying genetic changes. Toad-headed lizards of the genus Phrynocephalus cover a broad altitudinal gradient of over 4000 m and are useful models for studies of such adaptive responses. In one of the first studies to have considered selection on mitochondrial protein-coding regions in an ectothermic group distributed over such a wide range of environments, we analysed nineteen complete mitochondrial genomes from all Chinese Phrynocephalus (including eight genomes sequenced for the first time). Initial analyses used site and branch-site model (program: PAML) approaches to examine nonsynonymous: synonymous substitution rates across the mtDNA tree.

Results

Ten positively selected sites were discovered, nine of which corresponded to subunits ND2, ND3, ND4, ND5, and ND6 within the respiratory chain enzyme mitochondrial Complex I (NADH Coenzyme Q oxidoreductase). Four of these sites showed evidence of general long-term selection across the group while the remainder showed evidence of episodic selection across different branches of the tree. Some of these branches corresponded to increases in altitude and/or latitude. Analyses of physicochemical changes in protein structures revealed that residue changes at sites that were under selection corresponded to major functional differences. Analyses of coevolution point to coevolution of selected sites within the ND4 subunit, with key sites associated with proton translocation across the mitochondrial membrane.

Conclusions

Our results identify mitochondrial Complex I as a target for environment-mediated selection in this group of lizards, a complex that frequently appears to be under selection in other organisms. This makes these lizards good candidates for more detailed future studies of molecular evolution.

Electronic supplementary material

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

Regulation of Mammalian 13-Subunit Cytochrome c Oxidase and Binding of other Proteins: Role of NDUFA4

Cytochrome c oxidase (CcO) is the final oxygen accepting enzyme complex (complex IV) of the mitochondrial respiratory chain. In contrast to the other complexes (I, II, and III), CcO is highly regulated via isoforms for six of its ten nuclear-coded subunits, which are differentially expressed in species, tissues, developmental stages, and cellular oxygen concentrations. Recent publications have claimed that NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 4 (NDUFA4), originally identified as subunit of complex I, represents a 14th subunit of CcO. Results on CcO composition in tissues from adult animals and the review of data from recent literature strongly suggest that NDUFA4 is not a 14th subunit of CcO but may represent an assembly factor for CcO or supercomplexes (respirasomes) in mitochondria of growing cells and cancer tissues.

An X-chromosome linked mouse model (Ndufa1S55A) for systemic partial Complex I deficiency for studying predisposition to neurodegeneration and other diseases

The respiratory chain Complex I deficiencies are the most common cause of mitochondrial diseases. Complex I biogenesis is controlled by 58 genes and at least 47 of these cause mitochondrial disease in humans. Two of these are X-chromosome linked nuclear (nDNA) genes (NDUFA1 and NDUFB11), and 7 are mitochondrial (mtDNA, MT-ND1-6, -4L) genes, which may be responsible for sex-dependent variation in the presentation of mitochondrial diseases. In this study, we describe an X-chromosome linked mouse model (Ndufa1^S55A) for systemic partial Complex I deficiency. By homologous recombination, a point mutation T > G within 55th codon of the Ndufa1 gene was introduced. The resulting allele Ndufa1^S55A introduced systemic serine-55-alanine (S55A) mutation within the MWFE protein, which is essential for Complex I assembly and stability. The S55A mutation caused systemic partial Complex I deficiency of ∼50% in both sexes. The mutant males (Ndufa1^S55A/Y) displayed reduced respiratory exchange ratio (RER) and produced less body heat. They were also hypoactive and ate less. They showed age-dependent Purkinje neurons degeneration. Metabolic profiling of brain, liver and serum from males showed reduced heme levels in mutants, which correlated with altered expressions of Fech and Hmox1 mRNAs in tissues. This is the first genuine X-chromosome linked mouse model for systemic partial Complex I deficiency, which shows age-dependent neurodegeneration. The effect of Complex I deficiency on survival patterns of males vs. females was different. We believe this model will be very useful for studying sex-dependent predisposition to both spontaneous and stress-induced neurodegeneration, cancer, diabetes and other diseases.

Mitochondrial dysfunction in ocular disease: Focus on glaucoma

Mitochondrial dysfunction commonly presents with ocular findings as a part of a systemic disorder. These ophthalmic manifestations can be the first sign of a mitochondrial abnormality, which highlights the key role of a comprehensive ophthalmic assessment. On the other hand, a number of visually disabling genetic and acquired eye diseases with no curative treatment show abnormal mitochondrial function. Recent advances in mitochondrial research have improved our understanding of previously unexplained ocular disorders utilising better diagnostic approaches. Further studies on mitochondrial dysfunction and novel modalities of treatment will help to improve outcomes of these conditions. In this review article we discuss the clinical picture of common mitochondrial-related eye diseases, diagnostic approaches and possible treatment options including a very recent interesting report about gene therapy, with a particular focus on glaucoma.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.

Leber Hereditary Optic Neuropathy: Exemplar of an mtDNA Disease

The report in 1988 that Leber Hereditary Optic Neuropathy (LHON) was the product of mitochondrial DNA (mtDNA) mutations provided the first demonstration of the clinical relevance of inherited mtDNA variation. From LHON studies, the medical importance was demonstrated for the mtDNA showing its coding for the most important energy genes, its maternal inheritance, its high mutation rate, its presence in hundreds to thousands of copies per cell, its quantitatively segregation of biallelic genotypes during both mitosis and meiosis, its preferential effect on the most energetic tissues including the eye and brain, its wide range of functional polymorphisms that predispose to common diseases, and its accumulation of mutations within somatic tissues providing the aging clock. These features of mtDNA genetics, in combination with the genetics of the 1-2000 nuclear DNA (nDNA) coded mitochondrial genes, is not only explaining the genetics of LHON but also providing a model for understanding the complexity of many common diseases. With the maturation of LHON biology and genetics, novel animal models for complex disease have been developed and new therapeutic targets and strategies envisioned, both pharmacological and genetic. Multiple somatic gene therapy approaches are being developed for LHON which are applicable to other mtDNA diseases. Moreover, the unique cytoplasmic genetics of the mtDNA has permitted the first successful human germline gene therapy via spindle nDNA transfer from mtDNA mutant oocytes to enucleated normal mtDNA oocytes. Such LHON lessons are actively being applied to common ophthalmological diseases like glaucoma and neurological diseases like Parkinsonism.

Mitonuclear linkage disequilibrium in human populations

There is extensive evidence from model systems that disrupting associations between co-adapted mitochondrial and nuclear genotypes can lead to deleterious and even lethal consequences. While it is tempting to extrapolate from these observations and make inferences about the human-health effects of altering mitonuclear associations, the importance of such associations may vary greatly among species, depending on population genetics, demographic history and other factors. Remarkably, despite the extensive study of human population genetics, the statistical associations between nuclear and mitochondrial alleles remain largely uninvestigated. We analysed published population genomic data to test for signatures of historical selection to maintain mitonuclear associations, particularly those involving nuclear genes that encode mitochondrial-localized proteins (N-mt genes). We found that significant mitonuclear linkage disequilibrium (LD) exists throughout the human genome, but these associations were generally weak, which is consistent with the paucity of population genetic structure in humans. Although mitonuclear LD varied among genomic regions (with especially high levels on the X chromosome), N-mt genes were statistically indistinguishable from background levels, suggesting that selection on mitonuclear epistasis has not preferentially maintained associations involving this set of loci at a species-wide level. We discuss these findings in the context of the ongoing debate over mitochondrial replacement therapy.

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population-specific mitonuclear coadaptations create between-population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.

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.

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.

The complete mitochondrial genome of the yellow-browed bunting, Emberiza chrysophrys (Passeriformes: Emberizidae), and phylogenetic relationships within the genus Emberiza

Mitochondrial genomes have proved to be powerful tools in resolving phylogenetic relationships. Emberiza chrysophrys (least concern species: IUCN 2013) is a passerine bird in the bunting family, Emberizidae. The complete mitochondrial genome of E. chrysophrys was sequenced. This circular mitochondrial genome was 16,803 bp in length, with an A+T content of 52.26%, containing 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs and a putative control region (CR). The CR of E. chrysophrys was divided into three conserved domains. Six conserved sequence boxes in the central conserved domain II were identified as F, E, D, C, b and B. An obvious positive AT-skew and negative GC-skew bias were found for all 28 genes encoded by the H strand, whereas it was the reverse in the remaining nine genes encoded by the L strand. Remarkable rate heterogeneity was present in the mitochondrial genome of E. chrysophrys. Notably, unusual slow rate of evolution in the mitochondrial CR of E. chrysophrys was detected, which is rarely seen in other birds. Phylogenetic analyses were carried out based on 13 PCGs that showed E. pusilla was the sister group of E. rustica, and the monophyly of Emberiza was established.

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

Comparison of whole mitochondrial genome sequences from two clades of the invasive ascidian, Didemnum vexillum

The mitochondria are the main source of cellular energy production and have an important role in development, fertility, and thermal limitations. Adaptive mitochondrial DNA mutations have the potential to be of great importance in determining aspects of the life history of an organism. Phylogenetic analyses of the globally invasive marine ascidian Didemnum vexillum using the mitochondrial cytochrome c oxidase 1 (COX1) coding region, revealed two distinct clades. Representatives of one clade (denoted by 'B') are geographically restricted to D. vexillum's native region (north-west Pacific Ocean, including Japan), whereas members of the other clade (denoted by 'A') have been introduced and become invasive in temperate coastal areas around the world. Persistence of clade B's restricted distribution may reflect it being inherently less invasive than clade A. To investigate this we sought to determine if the two clades differ significantly in other mitochondrial genes of functional significance, specifically, alterations in amino acids encoded in mitochondrial enzyme subunits. Differences in functional mitochondrial genes could indicate an increased ability for clade A colonies to tolerate a wider range of environmental temperature. Full mitochondrial genomic sequences from D. vexillum clades A and B were obtained and they predict significant sequence differences in genes encoding for enzymes involved in oxidative phosphorylation. Diversity levels were relatively high and showed divergence across almost all genes, with p-distance values between the two clades indicating recent divergence. Both clades showed an excess of rare variants, which is consistent with balancing selection or a recent population expansion. Results presented here will inform future research focusing on examining the functional properties of the corresponding mitochondrial respiration enzymes, of A and B clade enzymes. By comparing closely related taxa that have differing distributions it is possible to identify genes and phenotypes suited to particular environments. The examination of mitochondrial genotypes, and associated enzyme functioning, across populations may aid in our understanding of thermal tolerance and environmental adaptation.

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.

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.

Gene expression profiling of mitochondrial oxidative phosphorylation (OXPHOS) complex I in Friedreich ataxia (FRDA) patients

Friedreich ataxia (FRDA) is the most frequent progressive autosomal recessive disorder associated with unstable expansion of GAA trinucleotide repeats in the first intron of the FXN gene, which encodes for the mitochondrial frataxin protein. The number of repeats correlates with disease severity, where impaired transcription of the FXN gene results in reduced expression of the frataxin protein. Gene expression studies provide insights into disease pathogenicity and identify potential biomarkers, an important goal of translational research in neurodegenerative diseases. Here, using real-time PCR (RT-PCR), the expression profiles of mitochondrial (mtDNA) and nuclear DNA (nDNA) genes that encode for the mitochondrial subunits of respiratory oxidative phosphorylation (OXPHOS) complex I in the blood panels of 21 FRDA patients and 24 healthy controls were investigated. Here, the expression pattern of mtDNA-encoded complex I subunits was distinctly different from the expression pattern of nDNA-encoded complex I subunits, where significant (p<0.05) down-regulation of the mitochondrial ND2, ND4L, and ND6 complex I genes, compared to controls, were observed. In addition, the expression pattern of one nDNA-encoded gene, NDUFA1, was significantly (p<0.05) down-regulated compared to control. These findings suggest, for the first time, that the regulation of complex I subunit expression in FRDA is complex, rather than merely being a reflection of global co-regulation, and may provide important clues toward novel therapeutic strategies for FRDA and mitochondrial complex I deficiency.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease

The unorthodox genetics of the mtDNA is providing new perspectives on the etiology of the common "complex" diseases. The maternally inherited mtDNA codes for essential energy genes, is present in thousands of copies per cell, and has a very high mutation rate. New mtDNA mutations arise among thousands of other mtDNAs. The mechanisms by which these "heteroplasmic" mtDNA mutations come to predominate in the female germline and somatic tissues is poorly understood, but essential for understanding the clinical variability of a range of diseases. Maternal inheritance and heteroplasmy also pose major challengers for the diagnosis and prevention of mtDNA disease.

Dynamics of mitochondrial inheritance in the evolution of binary mating types and two sexes

The uniparental inheritance (UPI) of mitochondria is thought to explain the evolution of two mating types or even true sexes with anisogametes. However, the exact role of UPI is not clearly understood. Here, we develop a new model, which considers the spread of UPI mutants within a biparental inheritance (BPI) population. Our model explicitly considers mitochondrial mutation and selection in parallel with the spread of UPI mutants and self-incompatible mating types. In line with earlier work, we find that UPI improves fitness under mitochondrial mutation accumulation, selfish conflict and mitonuclear coadaptation. However, we find that as UPI increases in the population its relative fitness advantage diminishes in a frequency-dependent manner. The fitness benefits of UPI ‘leak’ into the biparentally reproducing part of the population through successive matings, limiting the spread of UPI. Critically, while this process favours some degree of UPI, it neither leads to the establishment of linked mating types nor the collapse of multiple mating types to two. Only when two mating types exist beforehand can associated UPI mutants spread to fixation under the pressure of high mitochondrial mutation rate, large mitochondrial population size and selfish mutants. Variation in these parameters could account for the range of UPI actually observed in nature, from strict UPI in some Chlamydomonas species to BPI in yeast. We conclude that UPI of mitochondria alone is unlikely to have driven the evolution of two mating types in unicellular eukaryotes.

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.

Novel genomic resources for a climate change sensitive mammal: characterization of the American pika transcriptome

BACKGROUND

When faced with climate change, species must either shift their home range or adapt in situ in order to maintain optimal physiological balance with their environment. The American pika (Ochotona princeps) is a small alpine mammal with limited dispersal capacity and low tolerance for thermal stress. As a result, pikas have become an important system for examining biotic responses to changing climatic conditions. Previous research using amplified fragment length polymorphisms (AFLPs) has revealed evidence for environmental-mediated selection in O. princeps populations distributed along elevation gradients, yet the anonymity of AFLP loci and lack of available genomic resources precluded the identification of associated gene regions. Here, we harnessed next-generation sequencing technology in order to characterize the American pika transcriptome and identify a large suite of single nucleotide polymorphisms (SNPs), which can be used to elucidate elevation- and site-specific patterns of sequence variation.

RESULTS

We constructed pooled cDNA libraries of O. princeps from high (1400 m) and low (300 m) elevation sites along a previously established transect in British Columbia. Transcriptome sequencing using the Roche 454 GS FLX titanium platform generated 780 million base pairs of data, which were assembled into 7,325 high coverage contigs. These contigs were used to identify 24,261 novel SNP loci. Using high resolution melt analysis, we developed 17 of these SNPs into genotyping assays, which were validated with independent DNA samples from British Columbia Canada and Oregon State USA. In addition, we detected haplotypes in the NADH dehydrogenase subunit 5 of the mitochondrial genome that were fixed and different among elevations, suggesting that this may be an informative target gene for studying the role of cellular respiration in local adaptation. We also identified contigs that were unique to each elevation, including a high elevation-specific contig that was a positive match with the hemoglobin alpha chain from the plateau pika, a species restricted to high elevation steppes in Asia. Elevation-specific contigs may represent candidate regions subject to differential levels of gene expression along this elevation gradient.

CONCLUSIONS

To our knowledge, this is the first broad-scale, transcriptome-level study conducted within the Ochotonidae, providing novel genomic resources for studying pika ecology, behaviour and population history.

Complete mitochondrial genomes from four subspecies of common chaffinch (Fringilla coelebs): new inferences about mitochondrial rate heterogeneity, neutral theory, and phylogenetic relationships within the order Passeriformes

We describe whole mitochondrial genome sequences from four subspecies of the common chaffinch (Fringilla coelebs), and compare them to 31 publicly available mitochondrial genome sequences from other Passeriformes. Rates and patterns of mitochondrial gene evolution are analyzed at different taxonomic levels within this avian order, and evidence is adduced for and against the nearly neutral theory of molecular evolution and the role of positive selection in shaping genetic variation of this small but critical genome. We find evidence of mitochondrial rate heterogeneity in birds as in other vertebrates, likely due to differences in mutational pressure across the genome. Unlike in gadine fish and some of the human mitochondrial work we do not observe strong support for the nearly neutral theory of molecular evolution; instead evidence from molecular clocks, distribution of dN/dS ratios at different levels of the taxonomic hierarchy and in different lineages, McDonald-Kreitman tests within Fringillidae, and site-specific tests of selection within Passeriformes, all point to a role for positive selection, especially for the complex I NADH dehydrogenase genes. The protein-coding mitogenome phylogeny of the order Passeriformes is broadly consistent with previously-reported molecular findings, but provides support for a sister relationship between the superfamilies Muscicapoidea and Passeroidea on a short basal internode of the Passerida where relationships have been difficult to resolve. An unexpected placement of the Paridae (represented by Hume's groundpecker) within the Muscicapoidea was observed. Consistent with other molecular studies the mtDNA phylogeny reveals paraphyly within the Muscicapoidea and a sister relationship of Fringilla with Carduelis rather than Emberiza.

Accumulation of slightly deleterious mutations in the mitochondrial genome: a hallmark of animal domestication

The hypothesis that domestication leads to a relaxation of purifying selection on mitochondrial (mt) genomes was tested by comparative analysis of mt genes from dog, pig, chicken, and silkworm. The three vertebrate species showed mt genome phylogenies in which domestic and wild isolates were intermingled, whereas the domestic silkworm (Bombyx mori) formed a distinct cluster nested within its closest wild relative (Bombyx mandarina). In spite of these differences in phylogenetic pattern, significantly greater proportions of nonsynonymous SNPs than of synonymous SNPs were unique to the domestic populations of all four species. Likewise, in all four species, significantly greater proportions of RNA-encoding SNPs than of synonymous SNPs were unique to the domestic populations. Thus, domestic populations were characterized by an excess of unique polymorphisms in two categories generally subject to purifying selection: nonsynonymous sites and RNA-encoding sites. Many of these unique polymorphisms thus seem likely to be slightly deleterious; the latter hypothesis was supported by the generally lower gene diversities of polymorphisms unique to domestic populations in comparison to those of polymorphisms shared by domestic and wild populations.

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

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.