Metazoan OXPHOS gene families: evolutionary forces at the level of mitochondrial and nuclear genomes

Mitochondrial and nuclear DNAs contribute to encode the whole mitochondrial protein complement. The two genomes possess highly divergent features and properties, but the forces influencing their evolution, even if different, require strong coordination. The gene content of mitochondrial genome in all Metazoa is in a frozen state with only few exceptions and thus mitochondrial genome plasticity especially concerns some molecular features, i.e. base composition, codon usage, evolutionary rates. In contrast the high plasticity of nuclear genomes is particularly evident at the macroscopic level, since its redundancy represents the main feature able to introduce genetic material for evolutionary innovations. In this context, genes involved in oxidative phosphorylation (OXPHOS) represent a classical example of the different evolutionary behaviour of mitochondrial and nuclear genomes. The simple DNA sequence of Cytochrome c oxidase I (encoded by the mitochondrial genome) seems to be able to distinguish intra- and inter-species relations between organisms (DNA Barcode). Some OXPHOS subunits (cytochrome c, subunit c of ATP synthase and MLRQ) are encoded by several nuclear duplicated genes which still represent the trace of an ancient segmental/genome duplication event at the origin of vertebrates.     

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

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

Regulation of mitochondrial respiration and apoptosis through cell signaling: cytochrome c oxidase and cytochrome c in ischemia/reperfusion injury and inflammation

Cytochrome c (Cytc) and cytochrome c oxidase (COX) catalyze the terminal reaction of the mitochondrial electron transport chain (ETC), the reduction of oxygen to water. This irreversible step is highly regulated, as indicated by the presence of tissue-specific and developmentally expressed isoforms, allosteric regulation, and reversible phosphorylations, which are found in both Cytc and COX. The crucial role of the ETC in health and disease is obvious since it, together with ATP synthase, provides the vast majority of cellular energy, which drives all cellular processes. However, under conditions of stress, the ETC generates reactive oxygen species (ROS), which cause cell damage and trigger death processes. We here discuss current knowledge of the regulation of Cytc and COX with a focus on cell signaling pathways, including cAMP/protein kinase A and tyrosine kinase signaling. Based on the crystal structures we highlight all identified phosphorylation sites on Cytc and COX, and we present a new phosphorylation site, Ser126 on COX subunit II. We conclude with a model that links cell signaling with the phosphorylation state of Cytc and COX. This in turn regulates their enzymatic activities, the mitochondrial membrane potential, and the production of ATP and ROS. Our model is discussed through two distinct human pathologies, acute inflammation as seen in sepsis, where phosphorylation leads to strong COX inhibition followed by energy depletion, and ischemia/reperfusion injury, where hyperactive ETC complexes generate pathologically high mitochondrial membrane potentials, leading to excessive ROS production. Although operating at opposite poles of the ETC activity spectrum, both conditions can lead to cell death through energy deprivation or ROS-triggered apoptosis.     

Genome-wide analysis of cytochrome P450 monooxygenase genes in the silkworm, Bombyx mori

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed. A total of 84 CYP-related sequences were identified and could be classified into 26 families and 47 subfamilies according to standard nomenclature. Seventy eight of the eighty four genes appear to be functional and six are probable pseudogenes. The distribution of Bombyx mori P450s in the genome shows that most of them are tandem arranged on chromosomes, only 34 genes are present as singletons, with 8 clusters including 3 or more than 3 genes. Sequence alignments were used to reconstruct phylogenetic trees and to analyze the intron-exon organizations of the functional genes. The conserved intron positioning agrees perfectly with their common grouping on the tree. The presence of three extremely ancient introns which are conserved across different clans indicates that a few introns are still highly conserved after they have undergone extensive evolutionary changes of B. mori P450 duplication and divergence. Comparison of the P450s from B. mori to the P450s from Drosophila melanogaster shows that the expansion is not uniform across the gene families. Remarkably, two mitochondrial families, the B. mori CYP333 and D. melanogaster Cyp12, formed two orthologous groups in the phylogenetic tree. All CYP333s can be proposed to be related to xenobiotic metabolism in accordance with the D. melanogaster Cyp12s. The characterization and evolutionary analysis of P450s from B. mori in the current study provide useful information for understanding the characteristics and diversity of P450s from B. mori and the baseline for functional analyses of individual P450s in this model Lepidopteran insect.     

The MitoDrome database annotates and compares the OXPHOS nuclear genes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae

The oxidative phosphorylation (OXPHOS) is the primary energy-producing process of all aerobic organisms and the only cellular function under the dual control of both the mitochondrial and the nuclear genomes. Functional characterization and evolutionary study of the OXPHOS system is of great importance for the understanding of many as yet unclear aspects of nucleus-mitochondrion genomic co-evolution and co-regulation gene networks. The MitoDrome database is a web-based database which provides genomic annotations about nuclear genes of Drosophila melanogaster encoding for mitochondrial proteins. Recently, MitoDrome has included a new section annotating genomic information about OXPHOS genes in Drosophila pseudoobscura and Anopheles gambiae and their comparative analysis with their Drosophila melanogaster and human counterparts. The introduction of this new comparative annotation section into MitoDrome is expected to be a useful resource for both functional and structural genomics related to the OXPHOS system.     

Molecular evolution of cytochrome c oxidase: rate variation among subunit VIa isoforms

Cytochrome c oxidase (COX) consists of 13 subunits, 3 encoded in the mitochondrial genome and 10 in the nucleus. Little is known of the role of the nuclear-encoded subunits, some of which exhibit tissue-specific isoforms. Subunit VIa is unique in having tissue-specific isoforms in all mammalian species examined. We examined relative evolutionary rates for the COX6A heart (H) and liver (L) isoform genes along the length of the molecule, specifically in relation to the tissue-specific function(s) of the two isoforms. Nonsynonymous (amino acid replacement) substitutions in the COX6AH gene occurred more frequently than in the ubiquitously expressed COX6AL gene. Maximum-parsimony analysis and sequence divergences from reconstructed ancestral sequences revealed that after the ancestral COX6A gene duplicated to yield the genes for the H and L isoforms, the sequences encoding the mitochondrial matrix region of the COX VIa protein experienced an elevated rate of nonsynonymous substitutions relative to synonymous substitutions. This is expected for relaxed selective constraints after gene duplication followed by purifying selection to preserve the replacements with tissue-specific functions.      

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the ‘core’ proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome‐wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.     

The multiple functions of cytochrome c and their regulation in life and death decisions of the mammalian cell: From respiration to apoptosis

Cytochrome c (Cytc) is essential in mitochondrial electron transport and intrinsic type II apoptosis. Mammalian Cytc also scavenges reactive oxygen species (ROS) under healthy conditions, produces ROS with the co-factor p66(Shc), and oxidizes cardiolipin during apoptosis. The recent finding that Cytc is phosphorylated in vivo underpins a model for the pivotal role of Cytc regulation in making life and death decisions. An apoptotic sequence of events is proposed involving changes in Cytc phosphorylation, increased ROS via increased mitochondrial membrane potentials or the p66(Shc) pathway, and oxidation of cardiolipin by Cytc followed by its release from the mitochondria. Cytc regulation in respiration and cell death is discussed in a human disease context including neurodegenerative and cardiovascular diseases, cancer, and sepsis.     

Comparison of the oxidative phosphorylation (OXPHOS) nuclear genes in the genomes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae

Background

In eukaryotic cells, oxidative phosphorylation (OXPHOS) uses the products of both nuclear and mitochondrial genes to generate cellular ATP. Interspecies comparative analysis of these genes, which appear to be under strong functional constraints, may shed light on the evolutionary mechanisms that act on a set of genes correlated by function and subcellular localization of their products.

Results

We have identified and annotated the Drosophila melanogaster, D. pseudoobscura and Anopheles gambiae orthologs of 78 nuclear genes encoding mitochondrial proteins involved in oxidative phosphorylation by a comparative analysis of their genomic sequences and organization. We have also identified 47 genes in these three dipteran species each of which shares significant sequence homology with one of the above-mentioned OXPHOS orthologs, and which are likely to have originated by duplication during evolution. Gene structure and intron length are essentially conserved in the three species, although gain or loss of introns is common in A. gambiae. In most tissues of D. melanogaster and A. gambiae the expression level of the duplicate gene is much lower than that of the original gene, and in D. melanogaster at least, its expression is almost always strongly testis-biased, in contrast to the soma-biased expression of the parent gene.

Conclusions

Quickly achieving an expression pattern different from the parent genes may be required for new OXPHOS gene duplicates to be maintained in the genome. This may be a general evolutionary mechanism for originating phenotypic changes that could lead to species differentiation.     

Molecular genetic and clinical aspects of mitochondrial disorders in childhood

Mitochondrial OXPHOS disorders are caused by mutations in mitochondrial or nuclear genes, which directly or indirectly affect mitochondrial oxidative phosphorylation (OXPHOS). Primary mtDNA abnormalities in children are due to rearrangements (deletions or duplications) and point mutations or insertions. Mutations in the nuclear-encoded polypeptide subunits of OXPHOS result in complex I and II deficiency, whereas mutations in the nuclear proteins involved in the assembly of OXPHOS subunits cause defects in complexes I, III, IV, and V. Here, we review recent progress in the identification of mitochondrial and nuclear gene defects and the associated clinical manifestations of these disorders in childhood.     

The mitonuclear compatibility hypothesis of sexual selection

Why females assess ornaments when choosing mates remains a central question in evolutionary biology. We hypothesize that the imperative for a choosing female to find a mate with nuclear oxidative phosphorylation (OXPHOS) genes that are compatible with her mitochondrial OXPHOS genes drives the evolution of ornaments. Indicator traits are proposed to signal the efficiency of OXPHOS function thus enabling females to select mates with nuclear genes that are compatible with maternal mitochondrial genes in the formation of OXPHOS complexes. Species-typical pattern of ornamentation is proposed to serve as a marker of mitochondrial type ensuring that females assess prospective mates with a shared mitochondrial background. The mitonuclear compatibility hypothesis predicts that the production of ornaments will be closely linked to OXPHOS pathways, and that sexual selection for compatible mates will be strongest when genes for nuclear components of OXPHOS complexes are Z-linked. The implications of this hypothesis are that sexual selection may serve as a driver for the evolution of more efficient cellular respiration.     

Mouse histocompatibility-related genes are not conserved in other mammals

In addition to the genes for classical H-2 antigens, the H-2 complex of the mouse contains numerous homologous genes belonging to several distinct families. It is not known whether they have any functions. To address this question, I have investigated whether these genes are separately conserved in evolution. Subcloned 5' gene segments, encoding the variable domains, were used as hybridisation probes on genomic DNA blots of various mammals. Only the largest gene family, which includes the classical H-2/HLA genes, is detectable in humans and other mammals. The other gene families, including Qa-2 and T1a, are not conserved even in rodents. Most or all of their coding sequences are therefore redundant.     

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

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

CRIF1 Is Essential for the Synthesis and Insertion of Oxidative Phosphorylation Polypeptides in the Mammalian Mitochondrial Membrane

Although substantial progress has been made in understanding the mechanisms underlying the expression of mtDNA-encoded polypeptides, the regulatory factors involved in mitoribosome-mediated synthesis and simultaneous insertion of mitochondrial oxidative phosphorylation (OXPHOS) polypeptides into the inner membrane of mitochondria are still unclear. In the present study, disruption of the mouse Crif1 gene, which encodes a mitochondrial protein, resulted in a profound deficiency in OXPHOS caused by the disappearance of OXPHOS subunits and complexes in?vivo. CRIF1 was associated with large mitoribosomal subunits that were located close to the polypeptide exit tunnel, and the elimination of CRIF1 led to both aberrant synthesis and defective insertion of mtDNA-encoded nascent OXPHOS polypeptides into the inner membrane. CRIF1 interacted with nascent OXPHOS polypeptides and molecular chaperones, e.g., Tid1. Taken together, these results suggest that CRIF1 plays a critical role in the integration of OXPHOS polypeptides into the mitochondrial membrane in mammals.     

Phylogenetic analysis of vertebrate kininogen genes

Kininogens, the precursors of bradykinins, vary extremely in both structure and function among different taxa of animals, in particular between mammals and amphibians. This includes even the most conserved bradykinin domain in terms of biosynthesis mode and structure. To elucidate the evolutionary dynamics of kininogen genes, we have identified 19 novel amino acid sequences from EST and genomic databases (for mammals, birds, and fishes) and explored their phylogenetic relationships using combined amino acid sequence and gene structure as markers. Our results show that there were initially two paralogous kininogen genes in vertebrates. During their evolution, the original gene was saved with frequent multiplication in amphibians, but lost in fishes, birds, and mammals, while the novel gene was saved with multiple functions in fishes, birds, and mammals, but became a pseudogene in amphibians. We also propose that the defense mechanism against specific predators in amphibian skin secretions has been bradykinin receptor dependent. Our findings may provide a foundation for identification and structural, functional, and evolutionary analyses of more kininogen genes and other gene families.     

Association of Genes,Pathways, and Haplogroups of the Mitochondrial Genome with the Risk of Colorectal Cancer: The Multiethnic Cohort

The mitochondrial genome encodes for the synthesis of 13 proteins that are essential for the oxidative phosphorylation (OXPHOS) system. Inherited variation in mitochondrial genes may influence cancer development through changes in mitochondrial proteins, altering the OXPHOS process, and promoting the production of reactive oxidative species. To investigate the role of the OXPHOS pathway and mitochondrial genes in colorectal cancer (CRC) risk, we tested 185 mitochondrial SNPs (mtSNPs), located in 13 genes that comprise four complexes of the OXPHOS pathway and mtSNP groupings for rRNA and tRNA, in 2,453 colorectal cancer cases and 11,930 controls from the Multiethnic Cohort Study. Using the sequence kernel association test, we examined the collective set of 185 mtSNPs, as well as subsets of mtSNPs grouped by mitochondrial pathways, complexes, and genes, adjusting for age, sex, principal components of global ancestry, and self-reported maternal race/ethnicity. We also tested for haplogroup associations using unconditional logistic regression, adjusting for the same covariates. Stratified analyses were conducted by self-reported maternal race/ethnicity. In European Americans, a global test of all genetic variants of the mitochondrial genome identified an association with CRC risk (P = 0.04). In mtSNP-subset analysis, the NADH dehydrogenase 2 (MT-ND2) gene in Complex I was associated with CRC risk at a P-value of 0.001 (q = 0.015). In addition, haplogroup T was associated with CRC risk (OR = 1.66, 95% CI: 1.19–2.33, P = 0.003). No significant mitochondrial pathway and gene associations were observed in the remaining four racial/ethnic groups—African Americans, Asian Americans, Latinos, and Native Hawaiians. In summary, our findings suggest that variations in the mitochondrial genome and particularly in the MT-ND2 gene may play a role in CRC risk among European Americans, but not in other maternal racial/ethnic groups. Further replication is warranted and future studies should evaluate the contribution of mitochondrial proteins encoded by both the nuclear and mitochondrial genomes to CRC risk.