Search

Cytochrome c (CYC) oxidase (COX), a multisubunit enzyme that functions in mitochondrial aerobic energy production, catalyzes the transfer of electrons from CYC to oxygen and participates in creating the electrochemical gradient used for ATP synthesis. Modeling three-dimensional structural data on COX and CYC reveals that 57 of the >1,500 COX residues can be implicated in binding CYC. Because of the functional importance of the transfer of electrons to oxygen, it might be expected that natural selection would drastically constrain amino acid replacement rates of CYC and COX. Instead, in anthropoid primates, although not in other mammals, CYC and COX show markedly accelerated amino acid replacement rates, with the COX acceleration being much greater at the positions that bind CYC than at those that do not. Specifically, in the anthropoid lineage descending from the last common ancestor of haplorhines (tarsiers and anthropoids) to that of anthropoids (New World monkeys and catarrhines) and that of catarrhines (Old World monkeys and apes, including humans), a minimum of 27 of the 57 COX amino acid residues that bind CYC were replaced, most frequently from electrostatically charged to noncharged residues. Of the COX charge-bearing residues involved in binding CYC, half (11 of 22) have been replaced with uncharged residues. CYC residues that interact with COX residues also frequently changed, but only two of the CYC changes altered charge. We suggest that reducing the electrostatic interaction between COX and CYC was part of the adaptive evolution underlying the emergence of anthropoid primates.

Cytochrome c (CYC) and 9 of the 13 subunits of cytochrome c oxidase (complex IV; COX) were previously shown to have accelerated rates of nonsynonymous substitution in anthropoid primates. Cytochrome b, the mtDNA encoded subunit of ubiquinol-cytochrome c reductase (complex III), also showed an accelerated nonsynonymous substitution rate in anthropoid primates but rate information about the nuclear encoded subunits of complex III has been lacking. We now report that phylogenetic and relative rates analysis of a nuclear encoded catalytically active subunit of complex III, the iron-sulfur protein (ISP), shows an accelerated rate of amino acid replacement similar to cytochrome b. Because both ISP and subunit 9, whose function is not directly related to electron transport, are produced by cleavage into two subunits of the initial translation product of a single gene, it is probable that these two subunits of complex III have essentially identical underlying rates of mutation. Nevertheless, we find that the catalytically active ISP has an accelerated rate of amino acid replacement in anthropoid primates whereas the catalytically inactive subunit 9 does not.

Cytochrome c oxidase (COX) is a 13-subunit protein complex that catalyzes the last step in mitochondrial electron transfer in mammals. Of the 10 subunits encoded by nuclear DNA (three are mtDNA products), some are expressed as tissue- and/or development-specific isoforms. For COX subunit VIII, previous work showed that expression of the contractile muscle-specific isoform gene, COX8H, is absent in humans and Old World monkeys, and the other isoform gene, COX8L, is expressed ubiquitously. Here, we show that COX8H is transcribed in most primate clades, but its expression is absent in catarrhines, that is, in Old World monkeys and hominids (apes, including humans), having become a pseudogene in the stem of the catarrhines. The ubiquitously expressed isoform, COX8L, underwent nonsynonymous rate acceleration and elevation in the ratio of nonsynonymous/synonymous changes in the stem of anthropoid primates (New World monkeys and catarrhines), possibly setting the stage for loss of the heart-type (H) isoform. The most rapidly evolving region of VIII-L is one that interacts with COX I, suggesting that the changes are functionally coadaptive. Because accelerated rates of nonsynonymous substitutions in anthropoids such as observed for COX8L are also shown by genes for at least 13 other electron transport chain components, these encoded amino acid replacements may be viewed as part of a series of coadaptive changes that optimized the anthropoid biochemical machinery for aerobic energy metabolism. We argue that these changes were linked to the evolution of an expanded neocortex in anthropoid primates.

Mitochondrial DNA (mtDNA)-encoded proteins function in eukaryotes as subunits of respiratory complexes that also contain nuclear DNA (nDNA)-encoded subunits. The importance of functional interactions between mtDNA- and nDNA-encoded proteins was previously demonstrated by testing the survivability of cybrid cells or individuals containing nDNA and mtDNA from different populations or species. This report focuses on the multisubunit respiratory complex cytochrome c oxidase (COX), made up of both mtDNA-encoded and nDNA-encoded subunits. A combination of evolutionary and crystallographic data is employed to determine whether rates of nonsynonymous substitutions have been higher, the same, or lower for residues in close proximity that are encoded by a different genome (nDNA or mtDNA). This determination is performed by simply taking the ratio, called the interaction ratio i, of the nonsynonymous substitution rate of the close-contact residues to the nonsynonymous substitution rate of the noncontact residues. We assume that the close-contact residues (which are more likely to interact) are functionally important and that, therefore, amino acid replacements among these residues cannot escape the scrutiny of natural selection. i = 1 indicates that the close-contact residues have been under neither greater purifying selection nor greater positive selection than the noncontact residues as a specific consequence of their being encoded by separate genomes. i < 1 indicates that the close-contact residues have been under greater purifying selection but less positive selection than have the noncontact residues. Conversely, i > 1 indicates that the close-contact residues have been under less purifying but greater positive selection than have the noncontact residues. i < 1 may be referred to as a constraining interaction; i.e., the close-contact residues compared with the noncontact residues appear to be under greater structural-functional constraints. On the other hand, i > 1 may be referred to as an optimizing interaction; i.e., apparently many different amino acid replacements are required to optimize this subunit's interaction with the other subunit. A major finding is that the nDNA-encoded residues in close physical proximity to mtDNA-encoded residues evolve more slowly than the other nuclear-encoded residues (and thus display a constraining interaction), whereas the mtDNA-encoded residues in close physical proximity to nDNA-encoded residues evolve more rapidly than the other mitochondrial-encoded residues (and thus display an optimizing interaction). A possible reason for this striking difference between the nuclear- and mitochondrial-encoded COX subunits in how their functional interaction evolves is discussed.

As part of our goal to reconstruct human evolution at the DNA level, we have been examining changes in the biochemical machinery for aerobic energy metabolism. We find that protein subunits of two of the electron transfer complexes, complex III and complex IV, and cytochrome c, the protein carrier that connects them, have all undergone a period of rapid protein evolution in the anthropoid lineage that ultimately led to humans. Indeed, subunit IV of cytochrome c oxidase (COX; complex IV) provides one of the best examples of positively selected changes of any protein studied. The rate of subunit IV evolution accelerated in our catarrhine ancestors in the period between 40 to 18 million years ago and then decelerated in the descendant hominid lineages, a pattern of rate changes indicative of positive selection of adaptive changes followed by purifying selection acting against further changes. Besides clear evidence that adaptive evolution occurred for cytochrome c and subunits of complexes III (e.g., cytochrome c(1)) and IV (e.g., COX2 and COX4), modest rate accelerations in the lineage that led to humans are seen for other subunits of both complexes. In addition the contractile muscle-specific isoform of COX subunit VIII became a pseudogene in an anthropoid ancestor of humans but appears to be a functional gene in the nonanthropoid primates. These changes in the aerobic energy complexes coincide with the expansion of the energy-dependent neocortex during the emergence of the higher primates. Discovering the biochemical adaptations suggested by molecular evolutionary analysis will be an exciting challenge., (Copyright 2001 Academic Press.)

Phylogenetic analyses carried out on cytochrome c oxidase (COX) subunit I mitochondrial genes from 14 primates representing the major branches of the order and four outgroup nonprimate eutherians revealed that transversions and amino acid replacements (i.e., the more slowly occurring sequence changes) contained lower levels of homoplasy and thus provided more accurate information on cladistic relationships than transitions (i.e., the more rapidly occurring sequence changes). Several amino acids, each with a high likelihood of functionality involving the binding of cytochrome c or interaction with COX VIII, have changed in Anthropoidea, the primate suborder grouping New World monkey, Old World monkey, ape, and human lineages. They are conserved in other mammalian lineages and in nonanthropoid primates. Maximum-likelihood ancestral COX I nucleotide sequences were determined utilizing a near most parsimonious branching arrangement for the primate sequences that was consistent with previously hypothesized primate cladistic relationships based on larger and more diverse data sets. Relative rate tests of COX I mitochondrial sequences showed an elevated nonsynonymous (N) substitution rate for anthropoid-nonanthropoid comparisons. This finding for the largest mitochondrial (mt) DNA-encoded subunit is consistent with previous observations of elevated nonsynonymous substitution/synonymous substitution (S) rates in primates for mt-encoded COX II and for the nuclear-encoded COX IV and COX VIIa-H. Other COX-related proteins, including cytochrome c and cytochrome b, also show elevated amino acid replacement rates or N/S during similar time frames, suggesting that this group of interacting genes is likely to have coevolved during primate evolution., (Copyright 2000 Academic Press.)

COX VIIa is one of 10 nuclear-encoded subunits of the COX holoenzyme, and one of three that have isoforms with tissue-specific differences in expression. Analysis of nucleotide substitution rates revealed an accelerated rate of nonsynonymous substitutions relative to that of synonymous substitutions for the heart isoform gene (COX7AH) in six primate lineages. Rate accelerations have been noted for four other COX-related genes in this time period, suggesting that the COX holoenzyme has experienced an episode of adaptive evolution. A third member of the gene family, COX7AR, has recently been described. Although its function is currently unknown, low nonsynonymous substitution/synonymous substitution (N/S) ratios in mammalian evolution suggest that COX7AR is of functional importance. When the COX7A isoforms were divided into domains, examination of nucleotide substitution rates suggested that mitochondrial targeting residues experienced an accelerated nonsynonymous substitution rate in the period following gene duplication. In contrast, paralogous comparisons of the targeting residues of each isoform show they have been relatively conserved in mammalian evolution. This pattern is consistent with the evolution of tissue-specific function.

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.