Estimating the Influence of Selection on the Variable Amino Acid Sites of the Cytochrome b Protein Functional Domains

David A. McClellan and Kevin G. McCracken

Department of Zoology, Brigham Young University, Provo, Utah
Institute of Statistical Mathematics, Tokyo, Japan, and Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, Yokohama, Japan
Institute of Arctic Biology, and Department of Biology and Wildlife, University of Alaska, Fairbanks, Alaska


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 
We evaluated the effects of selection on the molecular evolution of the functional domains of the mammalian cytochrome b gene as it relates to physicochemical properties shown to correlate with rates of amino acid replacement. Two groups of mammals were considered: pocket gophers of the rodent family Geomyidae, and cetaceans and ungulates of the monophyletic taxon Cetartiodactyla. Several characteristics of cytochrome b evolution were common to both mammal groups. The evolution of the matrix domain reflected the region's relative lack of function. Goodness of fit to neutral expectations indicated that external influences have had very little effect on the evolution of the matrix, although in some cases conservative and moderate changes have been favored. Although rates of synonymous nucleotide substitution have been relatively high, the transmembrane domain exhibited poor goodness of fit to neutral expectations. However, the evolution of the transmembrane domain has been constrained by negative selection, allowing a preponderance of conservative and moderate amino acid replacements. We hypothesize that a high rate of substitution is maintained in spite of negative selection because the codons of the transmembrane coding region are predisposed to conservative changes in all amino acid properties. The evolutionary patterns of the intermembrane domain in pocket gophers and cetartiodactyls, however, were very different. Changes inferred from the pocket gopher phylogenetic tree exhibited a significant fit to neutral expectations for each of the amino acid properties. Changes inferred from the cetartiodactyl tree exhibited significant fit to neutral expectations for polarity and isoelectric point, but not for composition, molecular volume, polar requirement, or hydropathy. In each case, lack of fit was due to selection that promoted conservative or moderate change, with the noteworthy exception of polar requirement. We detected an unexpectedly large change in polar requirement (from aspartic acid to threonine) in two separate lineages (Camelus bactrianus and all cetaceans) at amino acid position 159. This inferred change occurred in a region of the cyt-b protein that directly interacts with external surface proteins of the cytochrome bc1 complex and resulted in a reversion to a more common character state in vertebrates.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 
Cytochrome b (cyt-b) is the central catalytic subunit of the Q-cycle and the only membrane-bound molecule in the cytochrome bc1 complex. The structure and function of cyt-b has been determined in great detail and in many cases down to the importance of individual amino acid residues (Degli Esposti et al. 1993Citation ; Zhang et al. 1998Citation ). Patterns of cyt-b nucleotide and amino acid substitution also have received much attention (e.g., Irwin, Kocher, and Wilson 1991Citation ; DeWalt et al. 1993Citation ; Graybeal 1993Citation ; Xia, Hafner, and Sudman 1996Citation ; Griffiths 1997Citation ). Such studies have evaluated cyt-b in terms of variation at many levels, such as among sites (e.g., codon positions), among base exchange types (i.e., transitions and transversions), and among regions or functional domains of the protein (i.e., intermembrane, transmembrane, and matrix).

Past studies of the molecular evolution of cyt-b indicate that the functional domains evolve at different rates of amino acid replacement (e.g., Irwin, Kocher, and Wilson 1991Citation ; DeWalt et al. 1993Citation ; Griffiths 1997Citation ). The common hypothesis for this phenomenon is that the selection responsible for maintaining the Q-cycle mechanism greatly influences rates of amino acid replacement. However, this hypothesis has yet to be critically evaluated in terms of evolutionary changes in the physicochemical properties that determine the phenotype of the protein. If the functional properties of the domains are being maintained by selection, very few amino acid replacements would be expected in those regions with the most residues implicated in the Q-cycle mechanism.

There are many physicochemical properties associated with amino acid residues (e.g., Sneath 1966Citation ; Woese et al. 1966Citation ; Alff-Steinberger 1969Citation ; Grantham 1974Citation ; Kyte and Doolittle 1982Citation ). Several of these have been shown to correlate with the evolution of the genetic code, the amino acid composition of polypeptides, and rates of amino acid replacement (e.g., Grantham 1974Citation ; Hughes, Ota, and Nei 1990Citation ; Haig and Hurst 1991Citation ; Xia and Li 1998Citation ). As a result, similar codons generally specify amino acids with similar physicochemical properties (Sonneborn 1965Citation ; Epstein 1966Citation ; Goldberg and Wittes 1966Citation ; Alff-Steinberger 1969Citation ). Furthermore, these properties contribute to a vast array of structural and functional characteristics of proteins that can be influenced directly by natural selection.

We evaluated the physicochemical changes that result from the molecular evolution of the cyt-b functional domains in pocket gophers and cetartiodactyls. Six amino acid properties shown to correlate with rates of amino acid replacement (Xia and Li 1998Citation ) were considered: composition of the side chain, polarity, and molecular volume (Grantham 1974Citation ), as well as polar requirement (Woese et al. 1966Citation ), hydropathy (Kyte and Doolittle 1982Citation ), and isoelectric point (Alff-Steinberger 1969Citation ). To evaluate the selective influences affecting the molecular evolution of cyt-b, we first calculated the goodness of fit between an observed distribution of physicochemical changes and an expected distribution based on codon composition (briefly outlined in McCracken et al. 1999Citation ) to measure how well the data fit neutral expectations. Next, we compared the inferred numbers of conservative, moderate, radical, and very radical physicochemical changes relative to the number of possible evolutionary pathways using a normal distribution to obtain relative measures (z-scores) of the selective influence on each amino acid property.

The goodness-of-fit score (GF-score) and the set of three z-scores can be interpreted as a measure of the selective influences on each gene region, with the magnitude of the GF-score indicative of the intensity of selection and the z-scores indicative of the direction of selection. For example, an observed distribution with a GF-score > 7.815 (df = 3, {alpha} = 0.05) and a z-score > 1.645 ({alpha} = 0.05) between the conservative and moderate magnitude classes with the conservative class being greater can be said to have experienced negative selection that favors conservative change. Alternatively, an observed distribution with a GF-score > 7.815 and a z-score > 1.645 between radical and very radical magnitude classes with the very radical class being greater can be said to have experienced positive selection that favors very radical change at those sites of amino acid replacement.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 
Pocket Gopher Phylogenetics
Cytochrome b sequences for 15 species of pocket gophers (Geomyidae: Rodentia) were downloaded from the National Center for Biotechnology Information. Phylogenetic analysis was performed using the general time-reversible (GTR) model of maximum-likelihood in PAUP*, version 4.0b1 (Swofford 1998Citation ), with maximum-likelihood estimates of substitution rates, nucleotide frequencies, proportion of invariable sites, and the gamma distribution shape parameter ({alpha}) producing the tree in figure 1A. This tree topology is consistent with trees constructed by Hafner et al. (1994)Citation and Xia, Hafner, and Sudman (1996)Citation using cytochrome oxidase subunit I (COI) and cyt-b.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 1.—Phylogenetic relationships of (A) pocket gophers and (B) cetartiodactyls

 
Cetartiodactyl Phylogenetics
Complete cytochrome b sequences for representatives of eight artiodactyl families and four cetaceans were downloaded from the National Center for Biotechnology Information. Phylogenetic analysis was performed using methods similar to those listed above for pocket gophers, with the exception that the tree search was constrained using the basic relationships published by Nikaido, Rooney, and Okada (1999)Citation based on the inferred occurrence of SINE and LINE insertions, which have been shown to be valuable phylogenetic markers (e.g., Nikaido, Rooney, and Okada 1999Citation ; Shedlock and Okada 2000Citation ; Shedlock, Milinkovitch, and Okada 2000Citation ). This analysis produced the tree in figure 1B.

Substitution Analysis
Ancestral character state reconstructions were optimized using the GTR maximum-likelihood model, with the same parameter estimates used to construct the tree. The sequence at each branch terminus on the tree was compared with its immediate ancestral sequence to estimate the number and types of substitutions (similar to Li 1993Citation ; see also DeWalt et al. 1993Citation ; Xia, Hafner, and Sudman 1996Citation ; Xia 1998Citation ; McClellan 2000Citation ).

Information pertaining to each inferred nonsynonymous substitution was recorded, including codon position, type of base exchange (transition or transversion), and exact location in the protein. Substitutions were categorized further by the functional domain in which they occurred (Degli Esposti et al. 1993Citation ; Zhang et al. 1998Citation ). The magnitudes of physicochemical changes in composition, polarity, and molecular volume (Grantham 1974Citation ), as well as polar requirement (Woese et al. 1966Citation ), hydropathy (Kyte and Doolittle 1982Citation ), and isoelectric point (Alff-Steinberger 1969Citation ), resulting from each inferred nonsynonymous substitution were also recorded.

Expected Distribution of Physicochemical Changes
Distributions of changes in physicochemical properties of inferred amino acid replacements were compared with distributions calculated based on the assumption of completely random amino acid replacement (McCracken et al. 1999Citation ; Xia 2000Citation ) expected under the condition of selective neutrality. Expected frequencies per physicochemical property magnitude class were calculated by


where Fi is the expected relative frequency of amino acid replacements of magnitude class i (i = 1, 2, ... , m, where m is the greatest magnitude class in which the condition Fi > 0 prevails), rij is the number of possible evolutionary pathways of physicochemical property magnitude class i across all possible nonsynonymous substitutions of codon j (j = 1, 2, ... , 60), and nj is the average number of times that codon j occurs in the extant DNA sequences. Equation (1) thus describes an average distribution of possible physicochemical differences. Thus, as selection decreases, the observed distribution of physicochemical differences of inferred amino acid replacements is predicted to take on a shape similar to the expected distribution.

The calculation of this expected distribution differs from that presented in McCracken et al. (1999)Citation in that the magnitude class frequencies are not transformed by the nonsynonymous substitution frequencies expected under the codon degeneracy model (McClellan 2000Citation ). This is a null model of amino acid replacement and, as it applies to the present study, assumes completely random replacement. Thus, modification of the expected distribution of physicochemical changes with regard to patterns of codon degeneracy and an inferred transition bias is unnecessary.

The distributions of physicochemical changes inferred from the phylogenetic analysis of the pocket gopher cyt-b DNA sequences were compared with the expected distribution calculated with equation (1) using a {chi}2 GF test of four magnitude categories: conservative, moderate, radical, and very radical. These categories are similar to those used by Wyckoff, Wang, and Wu (2000Citation ). Good fit (P > 0.05) between expected and observed distributions indicates that observed evolution of the DNA sequences is consistent with the assumption of selective neutrality. This is not to say, however, that the evolution of the DNA sequences is not influenced by selection. Good fit means only that the preponderance of sites experiencing amino acid replacements are accumulating changes as predicted by the model.

Detecting the Causes of Deviation from Expectations
The method developed by Hughes, Ota, and Nei (1990)Citation extends Nei and Gojobori's (1986)Citation method of estimating rates of synonymous and nonsynonymous nucleotide substitution. Nei and Gojobori (1986)Citation categorized each nucleotide site as either synonymous, nonsynonymous, or a fraction thereof. Hughes, Ota, and Nei (1990)Citation took this one step further by dividing each nonsynonymous site (or fractional nonsynonymous site) into either conservative, radical, or a fraction thereof based on the extent to which the possible evolutionary pathways result in either conservative or radical changes relative to a certain amino acid property. The number of inferred conservative substitutions is then divided by the mean number of conservative sites, and the number of inferred radical substitutions is divided by the mean number of radical sites. If these two ratios are equal, the amino acid replacements are said to take place at random relative to that particular property. If the ratio for more conservative replacements is greater than the ratio for more radical replacements, however, the property is being conserved, as would be expected under conditions of purifying selection. If the conservative ratio is less than the radical ratio, the replacements are said to promote radical changes, as would be expected under conditions of positive selection (Hughes, Ota, and Nei 1990Citation ).

One questionable aspect of this line of thought is the fractional assignment of a single nucleotide site as part conservative and part radical. We prefer to consider potential evolution by the number of possible nonsynonymous evolutionary pathways available to a given codon (Xia 1998Citation ; Xia and Li 1998Citation ). Another questionable aspect of the Hughes, Ota, and Nei (1990)Citation model is the consideration of just two categories of amino acid replacement (conservative and radical) based on the qualitative properties described by Taylor (1986)Citation . We prefer to use four categories of replacement with respect to quantifiable amino acid properties (as in Wyckoff, Wang, and Wu 2000Citation ). Considering four categories (conservative, p1; moderate, p2; radical, p3; and very radical, p4) allows a more detailed analysis of more dynamic evolutionary trends. Each ratio is calculated by


where rij and nj are as defined previously, and Di is the number of inferred nonsynonymous substitutions resulting in a physicochemical change of magnitude class i (i = 1, 2, 3, 4, depending on its relative conservative or radical nature). The standard error (SE) for each proportion is calculated by


and evaluated using a standard z test for comparing two binomial proportions.


    Results and Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 
Cytochrome b Functional Domain Evolution
The cyt-b protein can be divided into three functional domains: intermembrane, matrix, and transmembrane. The intermembrane domain, which extends between the inner and outer mitochondrial membranes, evolves significantly more slowly than the other two functional domains (Irwin, Kocher, and Wilson 1991Citation ; Griffiths 1997Citation ). The standard hypothesis for this phenomenon is that evolutionary constraints are imposed by the function of the Qo redox center (Irwin, Kocher, and Wilson 1991Citation ; DeWalt et al. 1993Citation ; Griffiths 1997Citation ). The intermembrane contains the greatest proportion of residues implicated in this function (Degli Esposti et al. 1993Citation ) and consists of 105 amino acid residues in mammals (Zhang et al. 1998Citation ), with about 29% of these either completely or mostly conserved in metazoans (Degli Esposti et al. 1993Citation ).

The matrix domain, which is located entirely within the inner surface of the inner mitochondrial membrane, is composed of 65 amino acid residues and has a relatively high proportion of polar and basic amino acid residues (Griffiths 1997Citation ). There are relatively few residues in the matrix domain that are conserved or have been implicated in the proton-input function of a Qi redox center (Degli Esposti et al. 1993Citation ). Most regions of the matrix domain have no known function (Irwin, Kocher, and Wilson 1991Citation ) or conserved properties in metazoans (Degli Esposti et al. 1993Citation ), and thus variable sites within this domain are expected to evolve in a nearly neutral manner relative to the other two functional domains.

The transmembrane domain, which consists of that portion of the cyt-b protein that transverses the inner mitochondrial membrane, is composed of 209 amino acid residues in mammals (Zhang et al. 1998Citation ) and is characterized by its hydrophobic properties (Irwin, Kocher, and Wilson 1991Citation ; Griffiths 1997Citation ). Furthermore, most of the observed amino acid replacements in the transmembrane region are among the hydrophobic amino acids leucine, isoleucine, and valine (Irwin, Kocher, and Wilson 1991Citation ; Kornegay et al. 1993Citation ). In addition, about 19% of the residues that compose the transmembrane domain are completely or nearly conserved in the vast majority of metazoans (Degli Esposti et al. 1993Citation ). Many of these residues have been implicated in the function of heme ligation, redox activity, or structural stability.

Summary of Cyt-b Data
Rates of synonymous substitution, rsyn, did not differ significantly between domains (table 1 ), with the exception being the matrix domain in cetartiodactyls. Rates of nonsynonymous substitution, rns, however, did differ significantly. Both the matrix and the transmembrane domains showed amino acid replacements significantly more frequently than the intermembrane domain in both geomyids and cetartiodactyls. The proportion of invariable amino acid sites, pinvar, corresponded with the pattern of rates of nonsynonymous substitution among the functional domains. The intermembrane domain exhibited a greater proportion of invariable sites than either the matrix or the transmembrane domain in both geomyids and cetartiodactyls, indicating that purifying selection is much greater in the intermembrane than in the other two domains. The number of nonsynonymous substitutions per variable site, however, did not appear to be greatly affected by selection in either group. These results suggest that variable amino acid sites in the intermembrane domain are influenced by much less selective constraint than the invariable sites, which experience an extreme form of purifying selection.


View this table:
[in this window]
[in a new window]
 
Table 1 Rates (mean ± SE) of Cytochrome b Gene Evolution

 
Comparisons of the relative proportions of the several amino acid residues that have undergone replacement in the cyt-b protein were indicative of the relative patterns of nonsynonymous substitution. For example, in pocket gophers, residues with a positive hydropathy (Kyte and Doolittle 1982Citation ) constituted about 32% of the matrix domain and 34% of the intermembrane domain, but 59% of the transmembrane domain. Furthermore, amino acid replacements at these sites constituted about 35%, 25%, and 88%, respectively, of all inferred replacements in these domains. Similar relationships between replacement sites and the magnitude of the physicochemical property were observed for composition (amino acid replacements resulting in no change in composition of the side chain: 24% in the matrix, 25% in the intermembrane, and 88% in the transmembrane), polarity (changes in polarity <= 6.2: 35%, 25%, and 86%, respectively), and polar requirement (changes in polar requirement <= 5.6: 35%, 25%, and 86%, respectively) and, to a lesser extent, for molecular volume (changes in molecular volume > 105: 35%, 38%, and 69%, respectively). The relationship between the magnitude of changes in these amino acid properties and the general rate of amino acid replacement underscores the importance of partitioning gene sequence data by functional domain when analyzing the molecular evolution of cyt-b.

Evolution Relative to Individual Amino Acid Properties
With regard to changes in physicochemical properties, the evolution of cyt-b in the two mammalian groups was similar (figs. 2 and 3 ). The matrix domain exhibited a good fit to expected magnitude distributions (P > 0.05) in relation to every property considered. The transmembrane domain, however, exhibited a poor fit (P < 0.05) to every expected distribution. Intermembrane domain evolution was not nearly as conservative across taxonomic groups. Pocket gophers exhibited a good fit to expected distributions for every amino acid property, whereas cetartiodactyls exhibited a good fit for polarity and isoelectric point, but not composition of the side chain, molecular volume, polar requirement, or hydropathy. The good fit of the pocket gopher data may be somewhat artificial, being the result of the extremely low rate of nonsynonymous substitution in the intermembrane domain (table 1 ).



View larger version (54K):
[in this window]
[in a new window]
 
Fig. 2.—Proportions of pocket gopher amino acid replacements per evolutionary pathway for the three functional domains of the cytochrome b protein and six physicochemical properties shown to correlate with rates of amino acid replacement. The x-axis gives the number of amino acid replacements per magnitude class (CONS. = conservative changes; MOD. = moderate changes; RAD. = radical changes; V. RAD. = very radical changes), indicated by the bars. Associated standard errors are also indicated. The y-axis gives the value of the z-score, indicated by the lines and circles. Levels of significance are also indicated: white circles, P > 0.05; gray circles, P < 0.05; black circles, P < 0.01. Arrows indicate direction of significant difference

 
Despite the good fit of the matrix domain to expected magnitude distributions, statistical comparisons between magnitude classes (conservative, moderate, radical, and very radical) revealed selection that has constrained evolution such that conservative and moderate changes predominate with regard to molecular volume in pocket gophers and with regard to all properties but hydropathy in cetartiodactyls. This trend may reflect a general evolutionary maintenance of some amino acid properties in the matrix domain. Furthermore, although the low GF scores denote an evolutionary profile that is not significantly different than neutral expectations (a conclusion that may correspond to the domain's general lack of function), overall conservation of some amino acid properties may nonetheless be important to the overall evolution of the cyt-b protein in mammals.

The relatively high rates of nonsynonymous substitution in the transmembrane domain (table 1 ) are not necessarily indicative of low levels of selection. The poor fit of the resultant amino acid replacements to neutral expectations supports this assertion, especially when compared with the GF scores calculated for the matrix domain, which exhibits similar rates of substitution. Poor GF scores with high rates of substitution in the transmembrane domain probably are due to the effect of codon composition.

Each codon possesses an inherent suite of possible amino acid changes precluding all others. As a result, some codons are inherently conservative by nature, whereas others are inherently radical. Codons that encode phenylalanine, leucine, isoleucine, methionine, and valine (codons with thymine at the second position) are examples of inherently conservative amino acids. Of the 104 possible evolutionary pathways for these 16 codons, only 12 result in nonconservative (i.e., moderate, radical, or very radical) changes with regard to composition of the side chain, 19 with regard to polarity, 16 with regard to molecular volume, 12 with regard to polar requirement, 26 with regard to hydropathy, and 5 with regard to isoelectric point. Collectively, these codons make up more than 52% of the transmembrane domain in both pocket gophers and cetartiodactyls. Phenylalanine, leucine, isoleucine, methionine, and valine are each hydrophobic and thus act to anchor the transmembrane domain into position through the membrane. The codons for these amino acids were the targets of a disproportionate number of nonsynonymous substitutions (72.5%) in the transmembrane domain, not because there are necessarily more mutations at these sites, but most probably because mutations at these sites are more likely to result in conservative changes, which in turn are more likely to become fixed. The correspondence between the magnitude of changes in amino acid properties and the general rate of amino acid replacement noted in the previous section may be the result of the genetic code playing a role that has been called "a guardrail set near the sharp ridge of the protein landscape" (Aita, Urata, and Husimi 2000). The parameters of the genetic code may constrain evolution in some regions, yet promote high, nearly neutral rates of evolution in others, and in both cases the effects of negative selection tend to amplify the results, in this case maintaining a rate of nonsynonymous substitution in the transmembrane domain that is comparable to the nonsynonymous rate of the matrix domain, which is affected by selection relatively very little (figs. 2 and 3 ).

As stated above, the pocket gopher intermembrane domain exhibited good fit to neutral expectations; the only possible exception was isoelectric point (fig. 2 ). Just one of the 17 amino acid replacements in this domain was "very radical," whereas the rest were conservative. Although these numbers were insufficient to constitute a lack of global fit, the z-scores comparing magnitude classes indicated a significant difference between conservative and moderate categories and between radical and very radical categories. This may be interpreted as a state of general negative selection, with enough positive selection for the replacement glutamine (Q) -> arginine (R) at amino acid position 137 to become fixed in a single taxon (Geomys bursarius) contrary to neutral expectations and conserved cyt-b primary structure (Degli Esposti et al. 1993Citation ).

The evolution of the cetartiodactyl intermembrane domain exhibited good fit for polarity and isoelectric point, although the latter appeared somewhat constrained. Composition, molecular volume, polar requirement, and hydropathy exhibited poor fits to expected distributions, indicating that there have been other influences on the evolution of these properties. Comparisons of observed magnitude classes suggested that conservative changes in molecular volume, as well as moderate changes in both composition and hydropathy, have been promoted by selection in this functional domain (fig. 3 ). Analysis of inferred changes in polar requirement suggested that not only has there been negative selection promoting conservative changes, but there has also been enough positive selection to promote the same "very radical" amino acid replacement, aspartic acid (D) -> threonine (T) at amino acid site 159, in two separate lineages (Camelus bactrianus and, most notably, all cetaceans). Although this amino acid residue has not been implicated in the Q-cycle mechanism (Degli Esposti et al. 1993Citation ), this change represents a large local decrease in polar requirement and may be of substantial evolutionary importance. The location of this residue is near the amino-terminal end in the second helix of the large cd-loop, which interacts directly with the other external surface domains of the bc1 complex. The effect of such a replacement has yet to be determined. However, a T at this position is a common characteristic shared by many vertebrates, whereas D is shared by relatively few (Degli Esposti et al. 1993Citation ).



View larger version (56K):
[in this window]
[in a new window]
 
Fig. 3.—Proportions of cetartiodactyl amino acid replacements per evolutionary pathway for the three functional domains of the cytochrome b protein and six physicochemical properties shown to correlate with rates of amino acid replacement. The x-axis gives the number of amino acid replacements per magnitude class (CONS. = conservative changes; MOD. = moderate changes; RAD. = radical changes; V. RAD. = very radical changes), indicated by the bars. Associated standard errors are also indicated. The y-axis gives the value of the z-score, indicated by the lines and circles. Levels of significance are also indicated: white circles, P > 0.05; gray circles, P < 0.05; black circles, P < 0.01. Arrows indicate direction of significant difference

 

    Conclusions
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 
Our results support previous hypotheses for different evolutionary rates of substitution in the cyt-b functional domains and propose mechanisms for their maintenance. The good fit of amino acid replacements to neutral expectations in the cyt-b matrix domain reflects this region's lack of function. Although the transmembrane domain shows a nonsynonymous substitution rate comparable to the matrix domain, transmembrane amino acid replacements do not fit neutral expectations for every amino acid property considered. Our analyses suggest that this phenomenon results from interaction between codon composition and purifying selection. As a result, purifying selection serves to "amplify" transmembrane amino acid replacement rates because the most plentiful amino acids are predisposed to conservative change. Not surprisingly, the intermembrane domain, which has the greatest number of amino acids implicated in the Q-cycle mechanism of the electron transport chain, evolves most slowly. Amino acid replacements in the pocket gopher intermembrane domain fit neutral expectations very well, whereas the cetartiodactyl intermembrane domain does not fit the model for four of the six amino acid properties. Positive selection resulting in "very radical" amino acid replacements was detected in both groups.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 
Nora J. Besansky, David W. Foltz, Mark S. Hafner, James V. Remsen, Frederick H. Sheldon, and two anonymous reviewers made helpful comments on the manuscript. Special thanks go to Norihiro Okada and Masami Hasegawa for office space and resources. This research was supported by grants from the Louisiana State University Museum of Natural Science, the American Museum of Natural History Theodore Roosevelt Fund, the American Society of Mammalogists, and the Japan Society for the Promotion of Science Monbusho Fund.


    Footnotes
 
David Irwin, Reviewing Editor

1 Keywords: codon composition physicochemical properties amino acid replacement cytochrome b functional domains Geomyidae Cetartiodactyla Back

2 Address for correspondence and reprints: David A. McClellan, Department of Zoology, 574 WIDB, Brigham Young University, Provo, Utah 84602. david_mcclellan{at}byu.edu Back


    literature cited
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Conclusions
 Acknowledgements
 literature cited
 

    Aita, T., S. Urata, and Y. Husimi. 2000. From amino acid landscape to protein landscape: analysis of genetic codes in terms of fitness landscape. J. Mol. Evol. 50:313–323[ISI][Medline]

    Alff-Steinberger, C.. 1969. The genetic code and error transmission. Proc. Natl. Acad. Sci. USA. 64:584–591[Abstract]

    Degli Esposti, M., S. De Vries, M. Crimi, A. Ghelli, T. Paternello, and A. Meyer. 1993. Mitochondrial cytochrome b: evolution and structure of the protein. Biochem. Biophys. Acta. 1143:243–271[ISI][Medline]

    DeWalt, T. S., P. D. Sudman, M. S. Hafner, and S. K. Davis. 1993. Phylogenetic relationships of pocket gophers (Cratogeomys and Pappogeomys) based on mitochondrial DNA cytochrome b sequences. Mol. Phylogenet. Evol. 2:193–204[Medline]

    Epstein, C. J.. 1966. Role of the amino-acid ‘code’ and of selection for conformation in the evolution of proteins. Nature. 210:25–28[ISI][Medline]

    Goldberg, A. L., and R. E. Wittes. 1966. Genetic code: aspects of organization. Science. 153:420–424[ISI][Medline]

    Grantham, R.. 1974. Amino acid difference formula to help explain protein evolution. Science. 185:862–864[ISI][Medline]

    Graybeal, A.. 1993. The phylogenetic utility of cytochrome b: lessons from bufonid frogs. Mol. Phylogenet. Evol. 2:256–269[Medline]

    Griffiths, C. S.. 1997. Correlation of functional domains and rates of nucleotide substitution in cytochrome b.. Mol. Phylogenet. Evol. 7:352–365[ISI][Medline]

    Hafner, M. S., P. D. Sudman, F. X. Villablanca, T. A. Spradling, J. W. Demastes, and S. A. Nadler. 1994. Disparate rates of molecular evolution in cospeciating hosts and parasites. Science. 265:1087–1090[ISI][Medline]

    Haig, D., and L. D. Hurst. 1991. A quantitative measure of error minimization in the genetic code. J. Mol. Evol. 33:412–417[ISI][Medline]

    Hughes, A. L., T. Ota, and M. Nei. 1990. Positive Darwinian selection promotes charge profile diversity in the antigen-binding cleft of class I major-histocompatibility-complex molecules. Mol. Biol. Evol. 7:515–524[Abstract]

    Irwin, D. M., T. D. Kocher, and A. C. Wilson. 1991. Evolution of the cytochrome b gene of mammals. J. Mol. Evol. 32:128–144[ISI][Medline]

    Kornegay, J. R., T. D. Kocher, L. A. Williams, and A. C. Wilson. 1993. Pathways of lysozyme evolution inferred from the sequences of cytochrome b in birds. J. Mol. Evol. 37:367–379[ISI][Medline]

    Kyte, J., and R. F. Doolittle. 1982. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 157:105–132[ISI][Medline]

    Li, W.-H.. 1993. Unbiased estimation of the rates of synonymous and nonsynonymous substitution. J. Mol. Evol. 36:96–99[ISI][Medline]

    McClellan, D. A.. 2000. The codon-degeneracy model of molecular evolution. J. Mol. Evol. 50:131–140[ISI][Medline]

    McCracken, K. G., J. Harshman, D. A. McClellan, and A. D. Afton. 1999. Data set incongruence and correlated character evolution: an example of functional convergence in the hind-limbs of stifftail diving ducks. Syst. Biol. 48:683–714[ISI][Medline]

    Nei, M., T. Gojobori. 1986. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol. Biol. Evol. 3:418–426[Abstract]

    Nikaido, M., A. P. Rooney, and N. Okada. 1999. Phylogenetic relationships among cetartiodactyls based on insertions of short and long interspersed elements: hippopotamuses are the closest extant relatives of whales. Proc. Natl. Acad. Sci. USA. 96:10261–10266[Abstract/Free Full Text]

    Shedlock, A. M., M. C. Milinkovitch, and N. Okada. 2000. SINE evolution, missing data, and the origin of whales. Syst. Biol. 49:1–10[ISI]

    Shedlock, A. M., and N. Okada. 2000. SINE insertions: powerful tools for molecular systematics. BioEssays. 22:148–160[ISI][Medline]

    Sneath, P. H. A.. 1966. Relations between chemical structure and biological activity. J. Theor. Biol. 12:157–195[ISI][Medline]

    Sonneborn, T. M.. 1965. Degeneracy of the genetic code: extent, nature, and genetic implicationsPp. 377–397 in V. Bryson and H. J. Vogel, eds. Evolving genes and proteins. Academic Press, New York

    Swofford, D. L.. 1998. PAUP*Phylogenetic analysis using parsimony (*and other methods). Version 4. Sinauer, Sunderland, Mass

    Taylor, W. R.. 1986. The classification of amino acid conservation. J. Theor. Biol. 119:205–218[ISI][Medline]

    Woese, C. R., D. H. Dugre, S. A. Dugre, M. Kondo, and W. C. Saxinger. 1966. On the fundamental nature and evolution of the genetic code. Cold Spring Harb. Symp. Quant. Biol. 31:723–736[ISI][Medline]

    Wyckoff, G. J., W. Wang, and C.-I. Wu. 2000. Rapid evolution of male reproductive genes in the descent of man. Nature. 403:304–308[ISI][Medline]

    Xia, X.. 1998. The rate heterogeneity of nonsynonymous substitutions in mammalian mitochondrial genes. Mol. Biol. Evol. 15:336–344[Abstract]

    ———.. 2000. Phylogenetic relationships among horseshoe crab species: effect of substitution models on phylogenetic analysis. Syst. Biol. 49:87–100[ISI][Medline]

    Xia, X., M. S. Hafner, and P. D. Sudman. 1996. On transition bias in mitochondrial genes of pocket gophers. J. Mol. Evol. 43:32–40[ISI][Medline]

    Xia, X., and W.-H. Li. 1998. What amino acid properties affect protein evolution?. J. Mol. Evol. 47:557–564[ISI][Medline]

    Zhang, Z., L. Huang, V. M. Shulmeister, Y.-I. Chi, K. K. Kim, L.-W. Hung, A. R. Crofts, E. A. Berry, and S.-H. Kim. 1998. Electron transfer by domain movement in cytochrome bc1. Nature. 392:677–684[ISI][Medline]

Accepted for publication February 5, 2001.