Recombination rate and protein evolution in yeast
© Connallon and Knowles; licensee BioMed Central Ltd. 2007
Received: 22 February 2007
Accepted: 27 November 2007
Published: 27 November 2007
Theory and artificial selection experiments show that recombination can promote adaptation by enhancing the efficacy of natural selection, but the extent to which recombination affects levels of adaptation across the genome is still an open question. Because patterns of molecular evolution reflect long-term processes of mutation and selection in nature, interactions between recombination rate and genetic differentiation between species can be used to test the benefits of recombination. However, this approach faces a major difficulty: different evolutionary processes (i.e. negative versus positive selection) produce opposing relationships between recombination rate and genetic divergence, and obscure patterns predicted by individual benefits of recombination.
We use a combination of polymorphism and genomic data from the yeast Saccharomyces cerevisiae to infer the relative importance of nearly-neutral (i.e. slightly deleterious) evolution in different gene categories. For genes with high opportunities for slightly deleterious substitution, recombination substantially reduces the rate of molecular evolution, whereas divergence in genes with little opportunity for slightly deleterious substitution is not strongly affected by recombination.
These patterns indicate that adaptation throughout the genome can be strongly influenced by each gene's recombinational environment, and suggest substantial long-term fitness benefits of enhanced purifying selection associated with sexual recombination.
Genetic drift is expected to overpower natural selection when selection is weak and effective population size (N e ) is small [1–3]. Recombination increases the effective population size in which genes evolve by reducing interference between linked loci under selection [4, 5]. As a result, recombination is expected to facilitate the spread of beneficial mutations and the elimination of deleterious mutations [6, 7]. Because recombination rates vary between different regions of a genome [e.g. yeast: ; Drosophila: ; Mammals: ; plants: ], adaptation at the molecular level might be strongly affected by each gene's recombinational environment – genes evolving in low recombination regions are expected to be poorly adapted relative to those in high recombination regions [12, 13].
Different processes of molecular evolution produce different correlations between recombination rate and nucleotide divergence between species1.
Model of evolution
Fitness effect of substitutions2
Correlation between recombination rate and divergence
s = 0
-1/(2N e ) <s < 0
s > 0
The relationship between recombination and divergence will be shaped by the predominant process of molecular evolution (i.e. neutral; slightly deleterious; adaptive). To test whether recombination facilitates adaptation throughout the genome (by enhancing purifying and positive selection), genes evolving under purifying selection and those evolving via positive selection should be analyzed separately, as each predicts a different relationship between divergence and recombination rate. Unfortunately, inferring the processes causing molecular divergence has traditionally been problematic without detailed within- and between-species genetic data , which limits the extent of the genome that can be analyzed.
Results & Discussion
Inferring the fitness effects of deleterious mutations
Highly expressed genes appear to evolve under stronger purifying selection than low-expressed genes (although the mechanistic basis of this pattern is still debated; [24–27]). Consequently, gene expression level is a good predictor of the average fitness effect of deleterious mutations. Experimental gene knockouts have also identified suites of genes that are essential for survival, while many others are nonessential. To the extent that whole gene knockout phenotypes reflect the fitness effects of individual mutations, mutations in essential genes are predicted to have larger fitness effects than mutations in nonessential genes [25, 28–30]. Lastly, proteins have variable numbers of interaction partners (protein-protein interactions per gene – PPI – range up to nearly 300 PPI in yeast; Connallon & Knowles unpub.), which indicate the level of constraint due to pleiotropy . Because individual mutations are likely to disrupt more cellular processes in genes with many PPI compared to those with few PPI, purifying selection is expected to be stronger in genes with many PPI .
Previous inferences of the strength of purifying selection acting on different gene categories are based on an observed elevated rate of nonsynonymous substitution in low-expressed nonessential genes with few PPI , which assumes that elevated rates of substitution are caused by genetic drift. An alternative possibility is that rapidly evolving genes undergo frequent bouts of positive selection. To test this assumption, we analyzed available polymorphism data from S. cerevisiae genes (see Methods and Additional file 1). Under a neutral/nearly-neutral model, patterns of within species polymorphism are expected to mirror patterns of interspecific substitution (i.e. genes with high substitutions rates also exhibit high levels of polymorphism ). Positive selection decouples patterns of polymorphism and divergence and is expected to increase the number of substitutions relative to polymorphisms .
Recombination and protein divergence
These results clearly show that recombination can influence the rate of protein evolution at a genome wide scale and that the impact of recombination rate variation is strongest for low-expressed, nonessential genes with few PPI. Associations between recombination and divergence rate cannot be explained by covariation between recombination rate and several variables that independently affect protein evolution (the effects were controlled for; see Methods). Estimates of the relative recombination rate between genes are coarse and limited by the quality of the S. cerevisiae recombination map, and there are potential evolutionary changes in recombination between S. cerevisiae and S. paradoxus. However, both of these factors will decrease the strength of associations between divergence and recombination, and will cause our test to be conservative.
Mutation bias is also unlikely to account for the effect of recombination on protein evolution. We present associations between recombination and divergence at nonsynonymous sites (dN) rather than between recombination and dN/dS ratios because synonymous sites are under selection in yeast [e.g. ]. Indeed, codon usage bias (F op ) is positively correlated with recombination (as previously reported ), most strongly for highly expressed genes, which presumably have stronger selection for optimal codons (highest 50% expression partial rrec-Fop= 0.214 vs. lowest 50% rrec-Fop= 0.100; highest 25% expression partial rrec-Fop= 0.249 vs. lowest 25% rrec-Fop= 0.089; P < 0.0001 for both comparisons). As a consequence, dS is negatively correlated with recombination (highest 50% expression partial rrec-Fop= -0.131 vs. lowest 50% rrec-Fop= -0.077, P = 0.031; highest 25% expression partial rrec-Fop= -0.122 vs. lowest 25% rrec-Fop= -0.075, P = 0.123). Furthermore, direct estimates indicate higher mutation rates in regions of high recombination [36, 37], which should make our tests conservative. Despite these caveats with respect to using dS to estimate underlying mutational dynamics across the genome, dN/dS produces nearly-identical patterns of covariation with recombination (see Additional file 1, Fig. S3).
The results are consistent with evolutionary theory suggesting that recombination enhances the efficacy of selection [e.g. [6, 12]]. Mutations with weak fitness effects respond to selection when the effective population size (N e ) is large, but evolve via genetic drift when N e is small. By increasing N e , recombination enhances the power of selection and minimizes genetic drift. Furthermore, the adaptive consequences of recombination may be extreme in yeast since most genes in the yeast genome can be defined by weak purifying selection (~75% of genes are nonessential; ~50% have one PPI; TC & LLK unpub.). Such genes also tend to reside in genomic regions with relatively low recombination frequencies (see Additional file 1, Fig. S4). The correlations revealed by the data are particularly striking when one considers the method by which the genes are partitioned. Functional genomic data permits classification of genes according to their relative opportunities for slightly deleterious evolution. However, multiple types of substitutions (i.e. slightly deleterious, neutral, and beneficial) are likely to contribute to each gene's total genetic divergence between species. This plurality should dampen patterns predicted by any single processes, and will cause the conclusions presented here to be conservative.
Why does highly-expressed gene divergence show no correlation with recombination rate? There are two major possibilities. If genes under stronger selection tend to experience fewer beneficial mutations, their overall divergence rate might be relatively unaffected by local recombination. This might occur because strongly-selected genes are closer to perfection than weakly-selected genes and therefore have less opportunity for improvement, or because tradeoffs via pleiotropy (as indicated by high PPI ) limit the opportunity for beneficial mutations [38–40]. Secondly, selection for beneficial mutations might be very strong. The adaptive impact of varying recombination rate (and thus N e ) is expected to decrease with the strength of selection (i.e. s; [3, 41]). If beneficial mutations tend to be strongly advantageous, they will tend to become fixed in low or in high recombinational environments.
This study shows that recombination reduces evolutionary divergence in genes under relatively weak purifying selection (e.g. low-expressed, nonessential, few PPI), and at best, marginally increases divergence in genes under strong purifying selection (e.g. highly expressed, essential, many PPI). This pattern suggests that enhanced purifying selection is a primary long-term benefit of recombination in nature. The efficient removal of deleterious mutations might increase the competitive ability of sexual species and contribute to the observed ubiquity of sexual reproduction in eukaryotes [6, 42].
While this interpretation is appealing and supported by both theory and data from other taxa [e.g. [12, 13]], it should be noted that inferences about the processes driving nucleotide divergence between species are tentative and reflect a major limitation of molecular divergence data. Future studies using entire-genome polymorphism and divergence data can add resolution by estimating the proportion of adaptive substitutions per gene [see ]. Such estimates, combined with inferences about the slightly deleterious substitution rate (derived from expression, essentiality and PPI data), will permit a much improved analysis of the benefits of recombination.
Publically available polymorphism data was obtained via  and . Genes with at least four samples from S. cerevisiae and at least one polymorphic site were included in the analysis, resulting in a dataset of 35 genes (lower 25% expression, n = 11; lower 50% expression, n = 12; upper 50% expression, n = 23; upper 25% expression, n = 17; essential genes, n = 7; nonessential genes, n = 28), comprising 34443 nonsynonymous, and 9975 synonymous nucleotide sites. The mean number of samples per gene was = 22. Orthologous sequences from S. paradoxus were obtained via BLAST search at .
Per gene recombination rates for S. cerevisiae are those reported by Gerton et al.  [data available at ]. These estimates refer to recombination rates per sexual generation. The total, per generation recombination rate during the evolutionary history of each gene is the product of the rate under sexual reproduction (R sex ) and the frequency of outcrossing (O c ); R TOT = R sex O c (modified from ). Because the exact value of O c will be the same for all genes in the genome, R sex per gene i, relative to R sex for other genes, will be the same as R TOT per gene i, relative to R TOT for other genes , and R sex should accurately capture relative rates of total recombination for each gene. Furthermore, because 0 <O c < 1 , the critical transition between R TOT ≈ 0 to R TOT > 0, predicted to most strongly impact the efficacy of selection , will be represented with yeast genes.
Protein divergence data for S. cerevisiae and S. paradoxus orthologous genes were kindly provided by D. Allan Drummond (see  for details). Gene expression values for S. cerevisiae were calculated from seven time period estimates during the diauxic shift by DeRisi et al.  [data available at ]. Average and maximum expression levels for the seven time periods were calculated using the method of Kliman et al. . Results for average gene expression (E avg ) across the time periods are presented here; results do not differ when maximum expression values (E max ) are used. Essential genes were identified through the GeneMerge database . Gene length, dispensability, protein-protein interactions and genome map positions were obtained from the Saccharomyces Genome Database. Genes with no known interaction partner were excluded from analyses involving PPI as a variable. Space between genes (SBG) was calculated by the method of Hey & Kliman . Recombination, expression, length, SBG and divergence estimates were available for 4786 genes in total, including essential genes with two or more PPI (n = 329), essential genes with 1 PPI (n = 195), nonessential genes with two or more PPI (n = 762), and nonessential genes with 1 PPI (n = 835). Our dataset is available upon request.
Population samples for each gene were aligned with ClustalW , available online, and manually adjusted. P n , P s , D n , and D s values were calculated with DnaSP, Version 4.10 . Watterson's estimate of silent nucleotide diversity (theta) was calculated by hand (as described in ). The complete polymorphism dataset is provided in Additional file 1.
Genes were classified into low and high expression level categories, based on quantile partitions. These correspond with ranges of lower 25%: 2.25 to 3.15 (n = 1197); lower 50%: 2.25 to 3.31 (n = 2393); upper 50%: 3.31 to 3.53 (n = 2393); and upper 25%: 3.53 to 4.54 (n = 1197) log mRNA abundance.
All divergence estimates were log10 transformed to facilitate linear comparisons, which are presented (values of dN = 0 were converted to dN = 0.0001 prior to log transformation); the results are robust and also obtained with nonparametric comparisons (TC & LLK unpub.). Partial correlation analysis was used to compare recombination rate with dN (the rate of nonsynonymous substitutions); the same results were obtained for the comparison between recombination rate and dN/dS. The partial r statistic reported here reflects the association between recombination rate and protein divergence after associations between gene expression, gene length, and SBG were removed. These factors are known to influence patterns of protein evolution [; TC & LLK unpub.], are all correlated with one another (i.e. recombination is positively correlated with expression and gene density, but negatively correlated with length), and can therefore give rise to spurious correlations between the variables of interest. All statistical analyses were carried out with JMP (SAS Institute). Statistical comparisons between r for different gene categories were carried out with software available online http://department.obg.cuhk.edu.hk/researchsupport/Correlation.asp. Bonferroni corrections for multiple comparisons (α/5; because of the 5 categories explored in Figure 2) were used to adjust P values of statistical significance.
We are grateful to D. Allan Drummond for providing genetic divergence estimates, and to Alex Kondrashov, Daven Presgraves, Akane Uesugi, Jianzhi Zhang, and two anonymous reviewers for valuable comments on an earlier version of the manuscript. This research was funded by a grant from the University of Michigan Rackham Graduate School and a training grant from the National Human Genome Research Institute (T32 HG00040) that supported T. C., and by a National Science Foundation grant (DEB-04-47224) to L. L. K.
- Fisher RA: The distribution of gene ratios for rare mutations. Proc Roy Soc Edinburgh. 1930, 50: 204-219.Google Scholar
- Wright S: Evolution in Mendelian populations. Genetics. 1931, 16: 97-159.PubMed CentralPubMedGoogle Scholar
- Kimura M: On the probability of fixation of mutant genes in a population. Genetics. 1962, 47: 713-719.PubMed CentralPubMedGoogle Scholar
- Hill WG, Robertson A: The effect of linkage on limits to artificial selection. Genet Res. 1966, 8: 269-294.View ArticlePubMedGoogle Scholar
- Felsenstein J: The evolutionary advantage of recombination. Genetics. 1974, 78: 737-756.PubMed CentralPubMedGoogle Scholar
- Rice WR: Experimental tests of the adaptive significance of sexual recombination. Nat Rev Genet. 2002, 3: 241-251. 10.1038/nrg760.View ArticlePubMedGoogle Scholar
- Otto SP, Lenormand T: Resolving the paradox of sex and recombination. Nat Rev Genet. 2002, 3: 252-261. 10.1038/nrg761.View ArticlePubMedGoogle Scholar
- Gerton JL, DeRisi J, Schroff R, Lichten M, Brown PO, Petes TD: Global mapping of meiotic recombination hotspots and coldspots in the yeast Saccharomyces cerevisiae. Proc Natl Acad Sci USA. 2000, 97: 11383-11390. 10.1073/pnas.97.21.11383.PubMed CentralView ArticlePubMedGoogle Scholar
- Hey J, Kliman RM: Interactions between natural selection, recombination and gene density in the genes of Drosophila. Genetics. 2002, 160: 595-608.PubMed CentralPubMedGoogle Scholar
- Jensen-Seaman MI, Furey TS, Payseur BA, Lu YT, Roskin KM, Chen CF, Thomas MA, Haussler D, Jacob HJ: Comparative recombination rates in the rat, mouse, and human genomes. Genome Res. 2004, 14: 528-538. 10.1101/gr.1970304.PubMed CentralView ArticlePubMedGoogle Scholar
- Mezard C: Meiotic recombination hotspots in plants. Biochem Soc Trans. 2006, 34: 531-534. 10.1042/BST0340531.View ArticlePubMedGoogle Scholar
- Barton NH: Linkage and the limits to natural selection. Genetics. 1995, 140: 821-841.PubMed CentralPubMedGoogle Scholar
- Presgraves DC: Recombination enhances protein adaptation in Drosophila melanogaster. Curr Biol. 2005, 15: 1651-1656. 10.1016/j.cub.2005.07.065.View ArticlePubMedGoogle Scholar
- Pal C, Papp B, Hurst LD: Does the recombination rate affect the efficiency of purifying selection? The yeast genome provides a partial answer. Mol Biol Evol. 2001, 18: 2323-2326.View ArticlePubMedGoogle Scholar
- Betancourt AJ, Presgraves DC: Linkage limits the power of natural selection in Drosophila. Proc Nat Acad Sci USA. 2002, 99: 13616-13620. 10.1073/pnas.212277199.PubMed CentralView ArticlePubMedGoogle Scholar
- Bachtrog D: Protein evolution and codon usage bias on the neo-sex chromosomes of Drosophila miranda. Genetics. 2003, 165 (3): 1221-1232.PubMed CentralPubMedGoogle Scholar
- Marais G, Domazet-Loso T, Tautz D, Charlesworth B: Correlated evolution of synonymous and nonsynonymous sites in Drosophila. J Mol Evol. 2004, 59: 771-779. 10.1007/s00239-004-2671-2.View ArticlePubMedGoogle Scholar
- Zhang Z, Parsch J: Positive correlation between evolutionary rate and recombination rate in Drosophila genes with male-biased expression. Mol Biol Evol. 2005, 22: 1945-1947. 10.1093/molbev/msi189.View ArticlePubMedGoogle Scholar
- Paland S, Lynch M: Transitions to asexuality result in excess amino acid substitutions. Science. 2006, 311: 990-992. 10.1126/science.1118152.View ArticlePubMedGoogle Scholar
- Birky CW, Walsh JB: Effects of linkage on rates of molecular evolution. Proc Natl Acad Sci USA. 1988, 85: 6414-6418. 10.1073/pnas.85.17.6414.PubMed CentralView ArticlePubMedGoogle Scholar
- Charlesworth B: The effect of background selection against deleterious mutations on weakly selected, linked variants. Genet Res. 1994, 63: 213-227.View ArticlePubMedGoogle Scholar
- Peck JR: A ruby in the rubbish: beneficial mutations, deleterious mutations and the evolution of sex. Genetics. 1994, 137: 597-606.PubMed CentralPubMedGoogle Scholar
- Eyre-Walker A: The genomic rate of adaptive evolution. Trends Ecol Evol. 2006, 21: 569-575. 10.1016/j.tree.2006.06.015.View ArticlePubMedGoogle Scholar
- Drummond DA, Bloom JD, Adami C, Wilke CO, Arnold FH: Why highly expressed proteins evolve slowly. Proc Natl Acad Sci USA. 2005, 102: 14338-14343. 10.1073/pnas.0504070102.PubMed CentralView ArticlePubMedGoogle Scholar
- Pal C, Papp B, Lercher MJ: An integrated view of protein evolution. Nat Rev Genet. 2006, 7: 337-348. 10.1038/nrg1838.View ArticlePubMedGoogle Scholar
- McInerney JO: The causes of protein evolutionary rate variation. Trends Ecol Evol. 2006, 21: 230-232. 10.1016/j.tree.2006.03.008.View ArticlePubMedGoogle Scholar
- Rocha EPC: The quest for the universals of protein evolution. Trends Genet. 2006, 22: 412-416. 10.1016/j.tig.2006.06.004.View ArticlePubMedGoogle Scholar
- Hirsh AE, Fraser HB: Protein dispensability and rate of evolution. Nature. 2001, 411: 1046-1049. 10.1038/35082561.View ArticlePubMedGoogle Scholar
- Wall DP, Hirsh AE, Fraser HB, Kumm J, Giaever G, Eisen MB, Feldman MW: Functional genomic analysis of the rates of protein evolution. Proc Natl Acad Sci USA. 2005, 102: 5483-5488. 10.1073/pnas.0501761102.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang J, He X: Significant impact of protein dispensability on the instantaneous rate of protein evolution. Mol Biol Evol. 2005, 22: 1147-1155. 10.1093/molbev/msi101.View ArticlePubMedGoogle Scholar
- He X, Zhang J: Toward a molecular understanding of pleiotropy. Genetics. 2006, 173: 1885-1891. 10.1534/genetics.106.060269.PubMed CentralView ArticlePubMedGoogle Scholar
- Kimura M: The Neutral Theory of Molecular Evolution. 1983, Cambridge, Cambridge University PressView ArticleGoogle Scholar
- McDonald JH, Kreitman M: Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991, 351: 652-654. 10.1038/351652a0.View ArticlePubMedGoogle Scholar
- Bierne N, Eyre-Walker A: The genomic rate of adaptive amino acid substitution in Drosophila. Mol Biol Evol. 2004, 21: 1350-1360. 10.1093/molbev/msh134.View ArticlePubMedGoogle Scholar
- Kliman RM, Irving N, Santiago M: Selection conflicts, gene expression, and codon usage trends in yeast. J Mol Evol. 2003, 57: 98-109. 10.1007/s00239-003-2459-9.View ArticlePubMedGoogle Scholar
- Datta A, Jinksrobertson S: Association of increased spontaneous mutation rates with high levels of transcription in yeast. Science. 1995, 268: 1616-1619. 10.1126/science.7777859.View ArticlePubMedGoogle Scholar
- Holbeck SL, Strathern JN: A role for REV3 in mutagenesisduring double stranded break repair in Saccharomyces cerevisiae. Genetics. 1997, 147: 1017-1024.PubMed CentralPubMedGoogle Scholar
- Fisher RA: The genetical theory of natural selection. 1930, Oxford: Clarendon PressView ArticleGoogle Scholar
- Orr HA: The population genetics of adaptation: the distribution of factors fixed during adaptive evolution. Evolution. 1998, 52: 935-949. 10.2307/2411226.View ArticleGoogle Scholar
- Otto SP: Two steps forward, one step back: the pleiotropic effects of favoured alleles. Proc Roy Soc Lond B. 2004, 271: 705-714. 10.1098/rspb.2003.2635.View ArticleGoogle Scholar
- Hedrick PW: Genetics of populations. 1983, Boston: Science Books InternationalGoogle Scholar
- Bell G: The masterpiece of nature. 1982, Berkeley: University of California PressGoogle Scholar
- Genbank. [http://www.ncbi.nlm.nih.gov/sites/entrez?db=Popset]
- Polymorphix Database. [http://pbil.univ-lyon1.fr/polymorphix/]
- The Saccharomyces Genome Database. [http://www.yeastgenome.org/]
- Data from Gerton et al. 2000. [http://derisilab14.ucsf.edu/hotspots/]
- Ruderfer DM, Pratt SC, Seidel HS, Kruglyak L: Population genomic analysis of outcrossing and recombination in yeast. Nat Genet. 2006, 38: 1077-1081. 10.1038/ng1859.View ArticlePubMedGoogle Scholar
- McVean GAT, Charlesworth B: The effects of Hill-Robertson interference between weakly selected mutations on patterns of molecular evolution and variation. Genetics. 2000, 155: 929-944.PubMed CentralPubMedGoogle Scholar
- DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686. 10.1126/science.278.5338.680.View ArticlePubMedGoogle Scholar
- Data from DeRisi et al. 1997. [http://cmgm.stanford.edu/pbrown/explore/array.txt]
- Castillo-Davis CI, Hartl DL: GeneMerge – post-genomic analysis, data mining, and hypothesis testing. Bioinformatics. 2003, 19: 891-892. 10.1093/bioinformatics/btg114. [http://www.oeb.harvard.edu/hartl/lab/publications/GeneMerge/GeneMerge.html]View ArticlePubMedGoogle Scholar
- ClustalW online. [http://www.ch.embnet.org/software/ClustalW-XXL.html]
- Rozas J, Sánchez-DelBarrio JC, Messegyer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Hein J, Schierup MH, Wiuf C: Gene Genealogies, Variation and Evolution: A Primer in Coalescent Theory. 2005, New York: Oxford University PressGoogle Scholar
- Wloch DM, Szafraniec K, Borts RH, Korona R: Direct estimate of the mutation rate and the distribution of fitness effects in the yeast Saccaromyces cerevisiae. Genetics. 2001, 159: 441-452.PubMed CentralPubMedGoogle Scholar
- Keightley PD, Lynch M: Toward a realistic model of mutations affecting fitness. Evolution. 2003, 57: 683-685.View ArticlePubMedGoogle Scholar
- Joseph SB, Hall DW: Spontaneous mutations in diploid Saccharomyces cerevisiae: More beneficial than expected. Genetics. 2004, 168: 1817-1825. 10.1534/genetics.104.033761.PubMed CentralView ArticlePubMedGoogle Scholar
- Shaw FH, Geyer CJ, Shaw RG: A comprehensive model of mutations affecting fitness and inferences for Arabidobsis thaliana. Evolution. 2002, 56: 453-463.View ArticlePubMedGoogle Scholar
- Shaw RG, Shaw FH, Geyer C: What fraction of mutations reduces fitness? A reply to Keightley and Lynch. Evolution. 2003, 57: 686-689.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.