* Department of Human Genetics
Department of Anthropology, Howard Hughes Medical Institute, University of Chicago
Correspondence: E-mail: wyckoffg{at}umkc.edu.
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: comparative genomics male-driven evolution mutation rate molecular evolution sex chromosomes synteny bins
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Although most researchers have accepted male-driven evolution in mammals as a biological reality and turned their efforts to estimating its strength in various lineages, two articles have presented arguments that call into question the existence of the phenomenon. McVean and Hurst (1997) argue against male-driven evolution using two primary observations. First, the average silent substitution rate on the rodent X chromosome (KsX) is so dramatically reduced relative to the autosomal rate (KsA) that the ratio of the two (KsX/KsA) is outside the theoretical limits of male-driven evolution (KsX/KsA as reported by McVean and Hurst [1997] produces an estimate of , the strength of male driven evolution, that exceeds infinity). Second, the mutation rate on the Y chromosome is broadly similar to that of the autosomes, yielding an
that is inconsistent with that derived from the X chromosome.
As an alternative to male-driven evolution, McVean and Hurst (1997) invoke a model to which we refer as "modulated mutability." They argue that, because recessive deleterious mutations in X-linked genes are hemizygous in males, their negative effect on fitness is amplified relative to similar mutations in autosomal genes; consequently, selection would favor the systematic reduction of mutation rate, specifically on the X chromosome. Their argument rests on the hypothesis that selection acts at the local level to balance the presumed metabolic costs of suppressing mutations with the fitness advantage of maintaining high genomic fidelity, a model that is biologically plausible (Kimura 1967; Leigh 1970, 1973; Kondrashov 1995; Sniegowski et al. 2000), although there is as yet little empirical evidence to support it, particularly in mammals. In contrast to the conclusions of McVean and Hurst (1997), Lercher et al. (2001) assert that Ks rates are highly heterogeneous among different chromosomes but broadly similar within the same chromosome. Although they did not assess the significance of Ks heterogeneity among different chromosomes, both trends were statistically significant when actual data were compared with simulated genomes. The authors argue that these trends render the mutation rates on the X chromosome indistinguishable from many of the autosomes. Although the empirical findings of Lercher, Williams, and Hurst (2001) directly contradict those of McVean and Hurst (1997), Lercher, Williams, and Hurst (2001) also argue against male-driven evolution, explaining their findings in language very similar to that of McVean and Hurst (1997). They suggest that an entire chromosome represents the main, perhaps least reducible, unit of variation in mutation rates; thus, the differential mutation rates that they observe are the result of a modulation of mutation by selection acting at the individual chromosome level.
Together, these papers have twofold significance: they call into question the applicability of male-driven evolution to the mammalian lineage, and they provide the only existing empirical evidence that mutation rates are locally modulated in mammals. Given the implications of these conclusions for comparative genomic studies, we statistically reexamined the trends in chromosomal mutation rates reported by these researchers. Specifically, we used coding-region data from the same two lineages (mouse-rat and human-mouse) to determine how well the modulated mutability and male-driven evolution models fit the data.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
We obtained orthologs from the two lineages using a two-step process. First, human-mouse orthologs were gathered from synteny maps (Homologene, NCBI) and the sequences for each gene were recovered from RefSeq. Next, mouse-rat orthologs were identified by Blast (Altschul et al. 1990). Rat orthologs to mouse genes were identified as those that were reciprocal best hits when compared with a mouse data file using nucleotide Blast, with an e-value cut off of 1 x 1050. Coding region sequences annotated from RefSeq were used to maximize the speed and utility of the Blast search. In many cases, other supporting information was available to corroborate orthology (e.g., gene name annotation or orthology reported in databases or previous literature).
The coding sequences were aligned in frame using the Pileup and Framealign programs from the GCG package (Genetics Computer Group, Madison, Wis.). Newdiverge from the same suite was used to analyze evolutionary rates (Ka and Ks) according to the Li (1993) method. Gene assignments to human or mouse chromosomes were gathered from Unigene as well as from the human-mouse synteny database.
Statistics
We estimated Ks and the nonsynonymous substitution rate (Ka) based on the Li (1993) method and used these calculations in estimating Ka/Ks. This ratio is commonly employed as a measure of constraint on the evolution of a gene (Li 1993). Values for Ks and Ka calculated by the maximum-likelihood method (Goldman and Yang 1994; Yang 1997; Yang and Nielsen 2000) were nearly indistinguishable from those generated by the Li (1993) method. Thus, we use substitution rates calculated with the Li method for the majority of the following statistical analysis, which affords greater comparability to previous work (McVean and Hurst 1997; Sharp 1997; Fay, Wyckoff, and Wu 2001).
To assess the heterogeneity of Ks rates between chromosomes, we performed pairwise Kolmogorov-Smirnov tests (Hollander and Wolfe 1973). This test evaluates the null hypothesis that two data sets are drawn from the same underlying distribution. In this regard, it is similar to the F-test; however, the parametric F-test loses statistical power when dealing with unequal sample sizes and data that are not normally distributed. The Kolmogorov-Smirnov test, in contrast, is robust under these conditions and was therefore deemed most appropriate to our data, which have these limitations. Kolmogorov-Smirnov tests were performed using the online calculator. (http://www.physics.csbsju.edu/stats/KS-test.html). This program generates a D statistic, which gauges the distance between two data sets, and a P-value that reflects its significance based on a critical table of D values and the degrees of freedom based on the two sample sizes.
At the subchromosomal level, we also investigated the Ks rates and Ka/Ks ratios of genes located on different synteny blocks. Synteny blocks contain widely varying numbers of genes (the smallest synteny block had three genes, whereas the largest had 207, with an average of 20 genes per block). In addition to the issue of unequal sample sizes, Ks rates and Ka/Ks ratios are nearly normally distributed across some blocks but depart radically from a normal distribution within others. Thus, testing the similarity of two blocks is a matter of choosing the test with the most appropriate properties without sacrificing statistical power. We used Kolmogorov-Smirnov tests (as described above) and confirmed the results of the nonparametric test by performing F-tests in Microsoft Excel. To complement the Kolmogorov-Smirnov and F-tests, which examine differences in the distributions of two data sets, we were also interested in examining the differences among synteny block means. To this end, we performed t-tests (technically, a special execution of F-tests) in Microsoft Excel to determine whether or not the mean Ks rates of synteny blocks are significantly different from one another. This confirmatory statistical approach allowed us to ensure that we were not losing statistical power by using the nonparametric Kolmogorov-Smirnov test.
We tested for correlation between synonymous and nonsynonymous substitution rates in two ways. At the chromosome level, we constructed matricies of D statistics from the Kolmogorov-Smirnov tests for synonymous and nonsynonymous substitution rates. We used the Mantel test (Mantel 1967; Legendre and Legendre 1998) (Mantel PPC, available for download at http://www.bio.sdsu.edu/pub/andy/MANTEL.html) with 200 randomizations to determine whether or not these matrices were correlated. Distance matrices describing the relationships among groups can be derived from any of a number of characteristics. The Mantel test determines whether two matrices describing the same groups are correlated with one another by calculating an initial correlation coefficient for the two real matrices and comparing that value to the correlation coefficient between one real matrix and a specified number of random permutations of the other. The significance of the correlation is gauged by how consistently the real correlation coefficient exceeds that of the permutations. In testing for correlations among chromosomal Ka and Ks rates, the Mantel test is appropriate because, unlike a simple linear correlation, which gives a Pearson correlation coefficient, the Mantel test corrects for the problem of nonindependence among group relationships (in this case, chromosomes). We also tested for correlation between Ka and Ks rates estimated by maximum likelihood (Goldman and Yang 1994), using Minitab for Windows XP to generate simple linear correlations.
Lastly, we examined mutation rates in the larger context of relationships among all chromosomes or among synteny blocks using discriminant function analysis (DFA) (Tabachnick and Fidell 2001). In DFA, data points are assigned to true groups and the power of particular parameters to predict the true group to which the data point belongs is assessed. In this case, chromosomes (or synteny blocks) are the true groups for individual genes and Ks is employed as the sole predictor of group affiliation. DFA generates a statistic called the squared Mahalanobis distance (D2). This is a measure of the distance of two groups from one another relative to a regional centroid. This centroid is derived from the characteristics of all groups. The D2 value is then used to calculate an F statistic using the equation given in Droessler (1981). The significance is assessed using a table of critical F values having degrees of freedom equivalent to the number of groups and number of predictors. Thus, in contrast to the Kolmogorov-Smirnov tests, which describe the relationship of two groups in isolation from all other groups, the D2 and F values allow us to view relationships among chromosomal or synteny-block Ks rates in the context of all available groups of data.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Using the formula = (4*[KsA/KsX] 3)/(3 2*[KsA/KsX]), we estimate
(the male-to-female ratio of germline mutation rate) to be 3.52 in the mouse-rat lineage and 2.45 in the human-mouse. Our estimates are in broad agreement with those published for the mammalian lineage by numerous authors (Shimmin, Chang, and Li 1993a; Chang and Li 1995; Makova and Li 2002); however, they are in stark contrast to the value of
calculated by McVean and Hurst (1997) for the mouse-rat lineage, which exceeds infinity. Thus, our results concur with the conclusion of Mural and Venter (2002) and Lercher, Williams, and Hurst (2001) that the reduction of KsX in rodents is much less dramatic than that reported by McVean and Hurst (1997), and we extend this observation to the human-mouse lineage.
McVean and Hurst (1997) also calculated the average mutation rate for the Y chromosome (KsY). They found average KsY to be similar to the average KsA value, rather than exceeding it as predicted by male-driven evolution. They further note a discrepancy between the estimate generated from Y-linked genes and that generated with X-linked genes. Given the dearth of data available for the Y chromosome (McVean's and Hurst's [1997] calculations employ only three genes), the estimates of KsY and the
derived from it are best viewed with caution. Although our analysis does not include Y-derived estimates because of the scarcity of genes, we note that a recent article found that when several pairs of X-Y homologs are compared for the human-mouse lineage, genes on the Y have, on average, Ks values double that of their X-chromosome counterparts (Wyckoff, Li, and Wu 2002).
Our estimates are reported solely for comparison with previous mammalian studies where the strength of
was calculated using coding sequence (see Chang et al. [1994], Agulnik et al. [1997], and Lawson and Hewitt [2002] for examples). Given that selection and variation in Ks across the genome confound calculations of
from coding-region data, we concur with the argument of Makova and Li (2002) that meaningful estimates of the male-mutation bias are best calculated from paired noncoding region data, where selection is minimal.
Mutation Rate of the X Chromosome Is an Outlier Relative to Autosomes
As figure 1 shows, Ks varies considerably from chromosome to chromosome, which is in keeping with the observations of Lercher, Williams, and Hurst (2001). To explore this phenomenon more fully, we statistically examined the differences among chromosomal rates to gauge the degree and significance of this heterogeneity. We used pairwise Kolmogorov-Smirnov tests between chromosomes to determine which chromosomes were significantly different from the others in terms of Ks distributions. In the human-mouse lineage, we find that the distribution of Ks rates on the human X chromosome is statistically significantly different (P < 0.05) from those on 21 of the 22 human autosomes. In the mouse-rat lineage, the mouse X chromosome is significantly different at the same level from 16 of the 19 mouse autosomes (tables 1 and 2). In both lineages, no other chromosome is so uniformly distinct from other autosomes. Thus, in keeping with theories of male-driven evolution, it is clear that the mammalian X chromosome is an outlier to the autosomes (more so than any other autosome is to the rest of the genome) and, therefore, cannot be treated as statistically identical to the autosomes.
|
|
According to Lercher, Williams, and Hurst (2001), the heterogeneity in chromosomal mutation rates is so great that chromosomes are individually distinct and have little relation to one another; however, our data present a more complex web of relationships among chromosomal mutation rates. Squared Mahalanobis distances (D2) were generated from a discriminant function analysis, with human chromosomes as the "true" groups and genic Ks for the human-mouse lineage as the predictor for chromosomal localization (see Materials and Methods). Discriminant function analysis considers the relationships among chromosomes in the context of the variation present throughout the entire data set (including chromosome-to-chromosome variation). Thus, it complements our pairwise Kolmogorov-Smirnov tests, which consider the significance of the difference between chromosomes in isolation from other comparisons. Confirming the Kolmogorov-Smirnov test results, F values derived from pairwise D2 values demonstrate that 47% of the distances between chromosomes are significant (Droessler 1981).
Interestingly, we find that within the human-mouse lineage, several chromosomes (21 and 19) have relatively high Ks rates, forming a small cluster that is significantly different from most other autosomes, even those that also have relatively high Ks rates (e.g., 16 and 4). Chromosomes 15 and 14 form a cluster of relatively low Ks rates, distinct from most other autosomes. This pattern of chromosomal clustering is mirrored in the mouse-rat lineage, although the separation between chromosomal clusters is weaker, and both pairwise distances and D2 values frequently do not achieve significance. Thus, chromosomes seem to cluster into "families" of similar synonymous substitution rates. This pattern motivated us to examine mutation rate heterogeneity further at the subchromosomal level.
Synteny Blocks May Be the Unit of Mutation Rate Variation
Given the fact that the chromosome-to-chromosome heterogeneity in Ks was substantial, if not as pervasive as that reported by Lercher, Williams, and Hurst (2001), we were interested in investigating trends within the autosomes that might help to explain this genome-wide heterogeneity. We gathered information about the locations of genes on chromosomes from NCBI (Homologene). In particular, we identified groups of genes located on evolutionarily conserved pieces of chromosomes (i.e., synteny blocks) between humans and mice. These blocks contain groups of genes that remain physically linked over evolutionary time as chromosomes themselves break and reshuffle. As such, they constitute ideal biologically meaningful units for examining mutation rate heterogeneity at the subchromosomal level (Mural and Venter 2002).
The Kolmogorov-Smirnov tests show that two synteny blocks whose Ks distributions are significantly different from one another are essentially as likely to be on the same chromosome as on different chromosomes. We specifically examined this in the human-mouse lineage, in which synteny blocks are clearly defined and there is high-quality curation for gene locations in both species (table 3). F-tests confirm both the percentage of and trend in positional relationships between significant comparisons. In addition to testing the heterogeneity of Ks distributions among synteny blocks, we used t-tests to determine how many pairwise comparisons of mean Ks rates were significant. Both the percentages and positional relationships of blocks whose mean rates were significantly different from one another mirror the results of the Kolmogorov-Smirnov and F-tests. In summary, pairwise tests of both the means and distributions of Ks rates between synteny blocks demonstrate that considerable variation exists between blocks regardless of whether or not they are positioned on the same chromosome.
|
Although our pairwise analyses in both lineages suggest that positional similarity in Ks rates does not extend across entire chromosomes, an exploratory discriminant function analysis using the human-mouse lineage reveals that the mean Ks rate of a synteny block might not be completely independent of the chromosome on which the block is located. Using human chromosome as the true group, we found that the mean Ks rate of a synteny block correctly predicts the chromosome on which the block is located 10.2% of the time, which is 2.25 times better than by chance alone (data not shown). This effect suggests a slight, yet perhaps real, tendency for blocks with similar substitution rates to be located on the same chromosome. However, this effect is sufficiently weak as to be undetectable by the Kolmogorov-Smirnov test. It is also worth noting that the ability to predict chromosomal locations of synteny blocks based on block-specific Ks is largely attributable to a few outlier chromosomes (i.e., the X chromosome, chromosome 19, and chromosome 21). It is also relevant to note that neighboring synteny blocks have quite different mutation rates, suggesting that Lercher, Williams, and Hurst's (2001) "near-neighbor" approach to local similarity may have obscured sharp differences in substitution rates that exist across synteny boundaries. Our analysis suggests that synteny blocks are more biologically relevant units at which to examine mutation rate variation than are whole chromosomes.
Measures of Constraint Suggest That Selection Acts at the Level of Syntenic Blocks
Another means of testing whether or not selection is acting to adjust mutation rates over the entire length of chromosomes, as hypothesized by Lercher, Williams, and Hurst (2001), is to examine the strength of selective constraint and positive selection at the chromosomal and subchromosomal levels. The Ka/Ks ratio is one such measure; a Ka/Ks near 1 is suggestive of neutrality, whereas most genes have Ka/Ks ratios much less than 1 (suggestive of functional constraint). Rarely, individual genes have a Ka/Ks of greater than 1 (suggestive of positive selection). The average Ka/Ks for a specific synteny block should be indicative of the net effect of constraint, drift, and positive selection for genes on that syntenic segment over a given period of evolutionary time. We compared the means and distributions of Ka/Ks between pairs of synteny blocks either within the same chromosome or between different chromosomes using the Kolmogorov-Smirnov test (see Materials and Methods). It shows that the Ka/Ks ratios of synteny blocks on the same chromosome are essentially just as likely to be significantly different from one another as synteny blocks on different chromosomes. This is inconsistent with the interpretation that selection acts on chromosomes as a single unit.
Correlation Between Ka and Ks May Explain Local Similarity in Ka
Lercher, Williams, and Hurst (2001) also report that Ka rates mirror the chromosomal trends noted for Ks, albeit more weakly. They invoke selection to explain this local similarity, reasoning that it cannot be explained as a function of local Ks similarity, because the two measures are independent when calculated by maximum likelihood. We tested this assumption in two ways.
First, we constructed two distance matrices to describe the relationships among chromosomes in the human-mouse lineage, one employing the D values from the Kolmogorov-Smirnov tests for Ks and the other using those from the Kolomogorov-Smirnov tests for Ka (tables 1 and 2). We then performed a Mantel test with 200 randomizations to assess the significance of the correlation between the two matrices (Mantel 1967; Legendre and Legendre 1998). This confirmed that the two matrices are significantly correlated (P < 0.005). We did not test for Ka and Ks matrix correlation in the mouse-rat lineage for two reasons. First, it seemed likely that the lineage for which we had more data would be more likely to reveal real underlying correlation. Second, we felt that any correlation derived from the mouse-rat lineage, which has a much smaller sample size and represents a considerably shorter evolutionary distance, may be more likely to be spurious.
We also tested the independence of Ka and Ks calculated by maximum likelihood using simple linear correlations and found significant correlation between Ka and Ks in both lineages (P < 0.01, data not shown). Thus, in contrast to the conclusions of Lercher, Williams, and Hurst (2001), as well as other studies contending that correlation between Ka and Ks is an artifact of the Li (1993) method (Smith and Hurst 1998), we find significant correlation between synonymous and nonsynonymous substitution rates as calculated by both methods. Consequently, we caution that there is a strong possibility that any local similarity identified in Ka rates is at least in part a function of the correlation to local Ks rates.
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Our examination of the two lineages in tandem brings our results into much closer agreement with previous literature and provides no impetus for rejecting male-driven evolution or for invoking a modulated mutability model. In contrast to the dramatic reduction of KsX relative to KsA as reported by McVean and Hurst (1997), our larger data set yields a moderate KsX reduction and an estimate of that is in line with those published previously. Although we find the difference between KsX and KsA to be less pronounced than did McVean and Hurst (1997), both pairwise Kolmogorov-Smirnov tests and DFA of chromosomal Ks rates clearly demonstrate that mutation rates on the X chromosome are statistically distinct from autosomal rates, in contradiction of Lercher, Williams, and Hurst's (2001) conclusion that the mammalian X chromosome is not an outlier.
Although our results support the observation of considerable autosomal heterogeneity in Ksrates (Lercher, Williams, and Hurst 2001; Ebersberger et al. 2002), discriminant function analysis suggests that simple pairwise comparisons of chromosomal rates may obscure a possible "familial" relationship among groups of autosomes in terms of mutation rate similarity. This led us to explore the question of within-chromosome heterogeneity in mutation rates. Because genes within synteny blocks share similarities in evolutionary history, physical characteristics (such as chromatin structure), and possibly coordinate expression, we hypothesized that synonymous and nonsynonymous rates might be shared among these smaller units, based in large part on these commonalities. The results from our analysis of mutation rates across synteny blocks is incompatible with the idea that there is wide-spread sharing of evolutionary rates across entire chromosomes. First, we find that synteny blocks on the same chromosome in the human-mouse lineage are essentially as likely to be significantly different in Ks rate as synteny blocks on different chromosomes. The near equality of within- chromosome and between-chromosome heterogeneity is incompatible with the idea that selection modulates the substitution rates of entire chromosomes as a unit. However, discriminant function analysis analysis using chromosomes as true groups and synteny-block-average Ks as the predictor does suggest that the Ks of a block bears some relationship to which chromosome it resides in.
Analysis of Ka/Ks within synteny blocks seems to bear out the implication that selection does not act strongly on the whole-chromosomal level. We find that, in both lineages, Ka/Ks ratios for synteny blocks on the same human or mouse chromosome (using human chromosome in the human-mouse lineage and mouse chromosome in the mouse-rat lineage) are just as likely to be significantly different as Ka/Ks ratio comparisons for synteny blocks between chromosomes when examined by pairwise Kolmogorov-Smirnov tests. This makes it difficult to invoke whole-chromosome selection as a possible explanation for chromosome-to-chromosome rate heterogeneity.
All of the available evidence suggests that male-driven evolution can explain the systematic drop in Ks rates for genes on the X chromosome in mammals. This suggests that global factors (such as the number of cell divisions in the germline) account for a large proportion of mutation rate disparity between the X chromosome and autosomes witnessed in a variety of studies (Shimmin, Chang, and Li 1993b; Chang and Li 1995; Makova and Li 2002), and by extension, perhaps between species as well. We therefore conclude that the high degree of variability in Ks rates across chromosomes must be debated separately from discussions of male-driven evolution. We also suggest that looking at evolutionarily conserved synteny blocks provides a useful biological basis for further analyzing genomic data (Mural and Venter 2002) and offers further insight into rates of substitutions across mammalian chromosomes. Given the current progress on the rat genome project and calls for sequencing of other primates, the data to test hypotheses explaining the large variance in chromosomal Ks rates should soon be available. We submit that, although a possible hypothesis for this variance must include "modulated mutability" as put forth by McVean and Hurst (1997), it is too soon to conclusively argue that this phenomenon exists in mammalian taxa, given the lack of strong empirical evidence and the possibility of alternative explanations.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
![]() |
Literature Cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Agulnik, A., C. Bishop, J. Lerner, S. Agulnik, and V. Solovyev. 1997. Analysis of mutation rates in the SMCY/SMCX genes shows that mammalian evolution is male driven. Mamm. Genome 8:134-138.[CrossRef][ISI][Medline]
Altschul, S. F., W. Gish, W. Miller, E. W. Myers, and D. J. Lipman. 1990. Basic local alignment search tool. J. Mol. Biol. 215:403-410.[CrossRef][ISI][Medline]
Bohossian, H. B., H. Skaletsky, and D. Page. 2000. Unexpectedly similar rates of nucleotide substitution found in male and female hominids. Nature 406:622-625.[CrossRef][ISI][Medline]
Chang, B., and W.-H. Li. 1995. Estimating the intensity of male-driven evolution in rodents by using X-linked and Y-linked UBE 1 genes and pseudogenes. J. Mol. Evol. 40:70-77.[ISI][Medline]
Chang, B., L. Shimmin, S. Shyue, D. Hewett-Emmett, and W. Li. 1994. Weak male-driven molecular evolution in rodents. Proc. Natl. Acad. Sci. USA 91:827-831.
Droessler, J. 1981. Craniometry and biological distance: biocultural continuity and change at the late Woodland-Mississippian interface. Research Series No.1. Center for American Archeology, Northwestern University, Evanston, Ill.
Ebersberger, I., D. Metzler, C. Schwarz, and S. Paabo. 2002. Genomewide comparison of DNA sequences between humans and chimpanzees. Am. J. Hum. Genet. 70:1490-1497.[CrossRef][ISI][Medline]
Fay, J. C., G. J. Wyckoff, and C.-I. Wu. 2001. Positive and negative selection on the human genome. Genetics 158:1227-1234.
Goldman, N., and Z. Yang. 1994. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol. Biol. Evol. 11:725-736.
Hollander, M., and D. A. Wolfe. 1973. Nonparametric statistical methods. John Wiley and Sons, New York.
Kimura, M. 1967. On the evolutionary adjustment of spontaneous mutation rates. Genet. Res. 9:23-34.[ISI]
Kondrashov, A. S. 1995. Modifiers of mutation selection balance. Genet. Res. 66:53-69.[ISI]
Lawson, L., and G. Hewitt. 2002. Comparison of substitution rates in ZFX and ZFY introns of sheep and goat related species supports the hypothesis of male-biased mutation rates. J. Mol. Evol. 54:54-61.[CrossRef][ISI][Medline]
Legendre, P., and L. Legendre. 1998. Numerical ecology. 2nd English edition. Elsevier, New York.
Leigh, E. G., Jr. 1970. Natural selection and mutability. Am. Nat. 104:301-305.[CrossRef][ISI]
Leigh, E. G., Jr. 1973. The evolution of mutation rates. Genetics 73:(suppl): 1-18.
Lercher, M. J., E. J. Williams, and L. D. Hurst. 2001. Local similarity in evolutionary rates extends over whole chromosomes in human-rodent and mouse-rat comparisons: implications for understanding the mechanistic basis of the male mutation bias. Mol. Biol. Evol. 18:2032-2039.
Li, W. H. 1993. Unbiased estimation of the rates of synonymous and nonsynonymous substitution [letter]. J. Mol. Evol. 36:96-99.[ISI][Medline]
Makova, K. D., and W.-H. Li. 2002. Strong male-driven evolution of DNA sequences in humans and apes. Nature 416:624-626.[CrossRef][ISI][Medline]
Mantel, N. 1967. The detection of disease clustering and a generalized regression approach. Cancer Res. 27:209-220.[ISI][Medline]
McVean, G. T., and L. D. Hurst. 1997. Evidence for a selectively favorable reduction in the mutation rate of the X chromosome. Nature 386:388-392.[CrossRef][ISI][Medline]
Miyata, T., H. Hayashida, K. Kuma, K. Mitsuyasu, and T. Yasunaga. 1987. Male-driven molecular evolution: a model and nucleotide sequence analysis. Cold Spring Harbor Symp. Quant. Biol. LII:863-868.
Mural, R. J., and J. C. Venter. 2002. A comparison of whole-genome shotgun-derived mouse chromosome 16 and the human genome. Science 296:1661-1671.
Sharp, P. 1997. In search of molecular Darwinism. Nature 385:111-112.[Medline]
Shimmin, L. C., B. H. Chang, and W. H. Li. 1993a. Male-driven evolution of DNA sequences. Nature 362:745-747.[CrossRef][ISI][Medline]
Shimmin, L. C., B. H. Chang, and W. H. Li. 1993b. Male-driven evolution of DNA sequences. Nature 362:745-747.[CrossRef][ISI][Medline]
Smith, N. G., and L. D. Hurst. 1998. Sensitivity of patterns of molecular evolution to alterations in methodology: a critique of Hughes and Yeager. J. Mol. Evol. 47:493-500.[ISI][Medline]
Sniegowski, P. D., P. J. Garrish, T. Johnson, and A. Shaver. 2000. The evolution of mutation rates: separating causes from conseqeunces. Bioessays 22:1057-1066.[CrossRef][ISI][Medline]
Tabachnick, B. G., and L. S. Fidell. 2001. Using multivariate statistics. 2nd edition. Allyn and Bacon, Boston.
Wyckoff, G. J., J. Li, and C.-I. Wu. 2002. Molecular evolution of functional genes on the mammalian Y chromosome. Mol. Biol. Evol. 19:1633-1636.
Yang, Z. 1997. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13:555-556.[Medline]
Yang, Z., and R. Nielsen. 2000. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evol. 17:32-43.