Relationship between amino acid composition and gene expression in the mouse genome
© Misawa et al; licensee BioMed Central Ltd. 2011
Received: 25 June 2010
Accepted: 27 January 2011
Published: 27 January 2011
Codon bias is a phenomenon that refers to the differences in the frequencies of synonymous codons among different genes. In many organisms, natural selection is considered to be a cause of codon bias because codon usage in highly expressed genes is biased toward optimal codons. Methods have previously been developed to predict the expression level of genes from their nucleotide sequences, which is based on the observation that synonymous codon usage shows an overall bias toward a few codons called major codons. However, the relationship between codon bias and gene expression level, as proposed by the translation-selection model, is less evident in mammals.
We investigated the correlations between the expression levels of 1,182 mouse genes and amino acid composition, as well as between gene expression and codon preference. We found that a weak but significant correlation exists between gene expression levels and amino acid composition in mouse. In total, less than 10% of variation of expression levels is explained by amino acid components. We found the effect of codon preference on gene expression was weaker than the effect of amino acid composition, because no significant correlations were observed with respect to codon preference.
These results suggest that it is difficult to predict expression level from amino acid components or from codon bias in mouse.
Codon bias is a phenomenon that refers to the differences in the frequencies of occurrence of synonymous codons among different genes . In the translation-selection model, natural selection is considered to be a cause of codon bias, because codon usage in highly expressed genes is biased toward "optimal" codons, i.e., codons corresponding to more abundant tRNAs in many organisms [2–6]. Methods have previously been developed to predict the expression level of genes from their nucleotide sequences, which is based on the observation that synonymous codon usage shows an overall bias toward a few codons called major codons [2–12].
Previous studies have provided clear evidence that the translation-selection model applies to some prokaryotes, such as Escherichia coli[6, 13], but not to all bacteria . Additionally, some evidence exists that suggests this model is also applicable to various eukaryotes, including Saccharomyces cerevisiae[14–17], Caenorhabditis elegans[18, 19], and the fruit fly [16, 20], and even to the vertebrate Xenopus laevis. However, the relationship between codon bias and expression level as proposed by the translation-selection model is less evident in mammals [22–30]. Urrutia and Hurst  found a weak correlation between gene expression levels and codon bias in human, but failed to find a relationship between this correlation and tRNA-gene copy numbers.
Amino acid content is also known to be dependent on gene expression level in some bacteria [31, 32], as well as in budding yeast . To determine why the relationship between codon bias and gene expression level, as proposed by the translation-selection model, is less evident in mammals, we investigated the correlations between the expression levels of genes and both the amino acid contents of genes and codon preference, in mouse. Subsequently, we compared the effect of gene expression on codon preference to the effect of gene expression on amino acid composition. We used the expression data of mouse genes contained in the InGap database .
Materials and methods
We obtained cDNA sequences of genes of Mus musculus from the ROUGE (http://www.kazusa.or.jp/rouge/index.html) database . In total, 449,444 codons from 1,182 genes were used. Mouse expression data were retrieved from the InGap database using cDNA microarray .
Amino Acid Contents and Codon Preference
We calculated the proportion of the amino acid contents of all genes. In order to examine the translation-selection model, we classified amino acids into two classes, i.e., C- and T-adapted, on the basis of tRNA-gene copy numbers in the mouse genome, since tRNA-gene copy numbers can be considered as a rough estimate of tRNA abundance . If natural selection is a cause of codon bias, codon usage in C-adapted amino acids of highly expressed genes will be biased toward C-ending codons and vice versa. First, we defined C-ending and T-ending codons; for instance, AGC is a C-ending codon and AGT is a T-ending codon. However, both encode the Ser residue. In the mouse genome, when the number of tRNAs complementary to C-ending codons for an amino acid is larger than the number of tRNAs that are complementary to T-ending codons, the amino acid is defined as a C-adapted amino acid. If the opposite is true, the amino acid is instead classified as a T-adapted amino acid. Furthermore, an amino acid is classified as T-adapted when the number of tRNAs that are complementary to C-ending codons is the same as the number of tRNAs that are complementary to T-ending codons. We obtained the number of tRNAs in the mouse genome from the GtRNAdb database . Ser, Leu, Pro, Arg, Ile, Thr, Val, and Ala are T-adapted amino acids, whereas Phe, Tyr, Cys, His, Asn, Ser, Asp, and Gly are C-adapted amino acids. Of note, Ser is encoded by TCT, TCC, TCA, TCG, AGC, and AGT. The number of tRNAs that are complementary to TCT is larger than the number that are complementary to TCC, whereas the number of tRNAs that are complementary to AGT is smaller than the number that are complementary to AGC. We considered the two types of codons that specifically encoded Ser. We compared the expression levels of genes to the nucleotide composition at the 3rd position of the codons. We conducted this comparison for all amino acids, including the for T-adapted, and C-adapted amino acids
Because of CpG hypermutability, the mutation rates of codons are affected by the 3' adjacent codon [36, 37]. Thus, the frequency of codon occurrence is dependent on the adjacent amino acid [36, 38, 39]. We analyzed the effect of adjacent nucleotides on amino acid composition. Specifically, we calculated the correlation between the proportion of the first and third nucleotides of the 3' adjacent codon in genes and the expression levels of those genes.
The Pearson product-moment correlation coefficients were calculated using R software . Because the probability density functions of the amino acid contents, codon preference, and expression levels are not known, we used a Kendall test, which is a nonparametric correlation test. Some of highly expressed genes might have specific sequences and functions. Thus, we eliminated the outliers from the data; we defined outliers as both the 5% of genes with the highest expression levels and the 5% of genes with the lowest expression levels. We also conducted multiple-regression analysis.
Correlation between amino acid contents and gene expression level
Correlation between amino acid abundance and gene expression level
After Elimination of Outliers
Codon preference and gene expression level
Correlation between the nucleotide composition at the 3rd position of codons and gene expression level
3' adjacent nucleotide
After elimination of outliers
All amino acids
T-adapted amino acids
C-adapted amino acids
Observed number of combinations of nucleotide at the third position and their 3' adjacent nucleotide
The 3' adjacent nucleotide
Variation of gene expression level
There is a large variation among gene expression level . More than 90% of variation of expression levels cannot be explained by amino acid components. Gene expression levels are known to be affected by many factors, such as 3'UTR lengths . Further study must be necessary.
Amino acid contents and gene expression level
To our knowledge, this is the first study that showed amino acid composition depends on the gene expression level in mouse. Previous study has shown that, in the case of budding yeast, some residues showed a positive correlation, and most of these residues were small . Furthermore, Akashi and Gojobori  showed an increase in the abundance of less energetically costly amino acids in highly expressed proteins. This study also suggested that natural selection for energetic efficiency appears to constrain the primary structures of the proteins of Bacillus subtilis and E. coli. Amino acid mutations that do not cause changes in protein functions may result in subtle, but evolutionarily important, fitness consequences through their effects on translation and metabolism.
We compared the estimates of the cost of amino acid synthesis from the above mentioned study  to the correlation coefficients presented in Table 1 (data not shown); however, we determined that the correlation was insignificant. It may be difficult to estimate the accurate metabolic cost of each amino acid, because the mouse obtains amino acids from food. Furthermore, the cost may depend on the environment. In the case of mouse, 10 amino acids, namely, Arg, His, Ile, Leu, Lys, Met, Phe, Thr, Try, and Val, are essential for natural growth . Thus, sparing the incorporation of His and Ile in highly expressed proteins may be advantageous to the mouse. However, Cys is not an essential amino acid, and is negatively correlated with gene expression. Of note, both amino acid composition and gene expression level may be influenced by protein functions. Furthermore, adaptive changes in protein sequences may overcome the increases in the metabolic cost, and the amino acid sequences may not be optimized for metabolic cost. Further study is necessary to elucidate these issues. Our results show that the coefficient of determination is very small so that it would be hard to predict expression level from amino acid contents in mammals.
Codon preference and gene expression level
We determined that the effect of codon preference on gene expression was weaker than the effect of amino acid composition, because no significant correlations were observed with respect to codon preference. This result is consistent with the relationship between codon bias and expression level, as proposed by the translation-selection model, is less evident in mammals [22–30]. In mammals, it would also be hard to predict expression level from codon bias.
Hypermutability of CpG dinucleotides  is one of major causes of codon substitution in mammalian genes [44–48]. CpG dinucleotides are often methylated at sites of cytosine (C); subsequently, the methylated C spontaneously deaminates to thymine (T) with a higher frequency than that of other types of point mutations . It has previously been estimated that approximately 14% of codon substitutions are caused by hypermutations at CpG sites . Furthermore, CpG hypermutation has been shown to affect the rate of amino acid substitution .
Table 2 shows that gene expression levels do not significantly affect codon preference in mouse. Furthermore, Table 3 indicates that the effect of codon preference is weaker than that of CpG hypermutability. Thus, the relationship between codon bias and gene expression level can be explained on the basis of the translation-selection model . This model proposes that codon usage in highly expressed genes is biased toward "optimal" codons, i.e., codons corresponding to more abundant tRNAs. This bias has been demonstrated to affect both elongation rate and accuracy [50, 51]. As shown in Table 2, the calculated negative correlation indicates that the codons used in this study are not optimal. In human and mouse genomes, the most frequently used codons  are not those with the most abundant tRNAs .
Recent studies [36, 39] have shown that CpG mutation rates in the non-coding regions of the human genome negatively correlate with the local GC content [53–56]. Isochores of the human genome  appear to be an influential factor that affects codon composition [53–56], and several studies have shown that this factor is related to gene expression levels [58, 59]. However, additional studies are necessary to confirm the relationship between codon bias and the positional effect of genes.
Plotkin et al. showed that codon usage for tissue-specific genes varies among the tissues in which such genes are expressed, thereby suggesting that this variation may be affected by differential tRNA-gene copy numbers in different tissues. However, this variability in codon usage among tissues is still under debate [22, 60, 61]. Nevertheless, it is noteworthy that codon substitutions are affected by adjacent codons [36, 39], and are therefore indirectly affected by adjacent amino acids . Amino acid frequencies may also be tissue specific, although additional studies are necessary to investigate the effect of CpG hypermutability on tissue-specific codon usage. Furthermore, codon bias in mammalian genomes should also be investigated with regard to the presence of CpG nucleotides [27, 28, 30].
In mouse, the effect of gene expression level on codon bias is weaker than both the effect of gene expression level on amino acid composition and the effect of CpG hypermutability on codon bias. However, to detect the effect of gene expression level on codon bias in mouse, a study of more genes is necessary.
We thank Dr. Osamu Ohara and all members of the Human Genetics Laboratory of Kazusa DNA research institutes for their encouragement and useful comments on the manuscript. The present study was supported by the National Project on "Next-generation Integrated Living Matter Simulation" of the Ministry of Education, Culture, Sports, Science and Technology (MEXT).
- Miyata T, Hayashida H, Yasunaga T, Hasegawa M: The preferential codon usages in variable and constant regions of immunoglobulin genes are quite distinct from each other. Nucleic Acids Res. 1979, 7 (8): 2431-2438. 10.1093/nar/7.8.2431.PubMedPubMed CentralView Article
- Akashi H, Eyre-Walker A: Translational selection and molecular evolution. Curr Opin Genet Dev. 1998, 8 (6): 688-693. 10.1016/S0959-437X(98)80038-5.PubMedView Article
- Willie E, Majewski J: Evidence for codon bias selection at the pre-mRNA level in eukaryotes. Trends Genet. 2004, 20 (11): 534-538. 10.1016/j.tig.2004.08.014.PubMedView Article
- Sharp PM, Bailes E, Grocock RJ, Peden JF, Sockett RE: Variation in the strength of selected codon usage bias among bacteria. Nucleic Acids Res. 2005, 33 (4): 1141-1153. 10.1093/nar/gki242.PubMedPubMed CentralView Article
- Karlin S, Barnett MJ, Campbell AM, Fisher RF, Mrazek J: Predicting gene expression levels from codon biases in alpha-proteobacterial genomes. Proc Natl Acad Sci USA. 2003, 100 (12): 7313-7318. 10.1073/pnas.1232298100.PubMedPubMed CentralView Article
- Roymondal U, Das S, Sahoo S: Predicting gene expression level from relative codon usage bias: an application to Escherichia coli genome. DNA Res. 2009, 16 (1): 13-30. 10.1093/dnares/dsn029.PubMedPubMed CentralView Article
- Henry I, Sharp PM: Predicting gene expression level from codon usage bias. Mol Biol Evol. 2007, 24 (1): 10-12. 10.1093/molbev/msl148.PubMedView Article
- Raghava GP, Han JH: Correlation and prediction of gene expression level from amino acid and dipeptide composition of its protein. BMC Bioinformatics. 2005, 6: 59-10.1186/1471-2105-6-59.PubMedPubMed CentralView Article
- Raghava GP, Han JH, Hwang DJ: ECGpred: Correlation and prediction of gene expression from nucleotide sequence. The Open Bioinformatics Journal. 2008, 2: 64-71. 10.2174/1875036200802010064.View Article
- Jansen R, Bussemaker HJ, Gerstein M: Revisiting the codon adaptation index from a whole-genome perspective: analyzing the relationship between gene expression and codon occurrence in yeast using a variety of models. Nucleic Acids Res. 2003, 31 (8): 2242-2251. 10.1093/nar/gkg306.PubMedPubMed CentralView Article
- Beer MA, Tavazoie S: Predicting gene expression from sequence. Cell. 2004, 117 (2): 185-198. 10.1016/S0092-8674(04)00304-6.PubMedView Article
- Coghlan A, Wolfe KH: Relationship of codon bias to mRNA concentration and protein length in Saccharomyces cerevisiae. Yeast. 2000, 16 (12): 1131-1145. 10.1002/1097-0061(20000915)16:12<1131::AID-YEA609>3.0.CO;2-F.PubMedView Article
- Ikemura T: Codon usage and tRNA content in unicellular and multicellular organisms. Mol Biol Evol. 1985, 2 (1): 13-34.PubMed
- Bennetzen JL, Hall BD: Codon selection in yeast. J Biol Chem. 1982, 257 (6): 3026-3031.PubMed
- Ikemura T: Correlation between the abundance of yeast transfer RNAs and the occurrence of the respective codons in protein genes. Differences in synonymous codon choice patterns of yeast and Escherichia coli with reference to the abundance of isoaccepting transfer RNAs. J Mol Biol. 1982, 158 (4): 573-597. 10.1016/0022-2836(82)90250-9.PubMedView Article
- Akashi H: Inferring weak selection from patterns of polymorphism and divergence at "silent" sites in Drosophila DNA. Genetics. 1995, 139 (2): 1067-1076.PubMedPubMed Central
- Akashi H: Translational selection and yeast proteome evolution. Genetics. 2003, 164 (4): 1291-1303.PubMedPubMed Central
- Duret L: tRNA gene number and codon usage in the C. elegans genome are co-adapted for optimal translation of highly expressed genes. Trends Genet. 2000, 16 (7): 287-289. 10.1016/S0168-9525(00)02041-2.PubMedView Article
- Marais G, Duret L: Synonymous codon usage, accuracy of translation, and gene length in Caenorhabditis elegans. J Mol Evol. 2001, 52 (3): 275-280.PubMed
- Moriyama EN, Powell JR: Codon usage bias and tRNA abundance in Drosophila. J Mol Evol. 1997, 45 (5): 514-523. 10.1007/PL00006256.PubMedView Article
- Musto H, Cruveiller S, D'Onofrio G, Romero H, Bernardi G: Translational selection on codon usage in Xenopus laevis. Mol Biol Evol. 2001, 18 (9): 1703-1707.PubMedView Article
- Urrutia AO, Hurst LD: Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection. Genetics. 2001, 159 (3): 1191-1199.PubMedPubMed Central
- Urrutia AO, Hurst LD: The signature of selection mediated by expression on human genes. Genome Res. 2003, 13 (10): 2260-2264. 10.1101/gr.641103.PubMedPubMed CentralView Article
- Plotkin JB, Robins H, Levine AJ: Tissue-specific codon usage and the expression of human genes. Proc Natl Acad Sci USA. 2004, 101 (34): 12588-12591. 10.1073/pnas.0404957101.PubMedPubMed CentralView Article
- Lavner Y, Kotlar D: Codon bias as a factor in regulating expression via translation rate in the human genome. Gene. 2005, 345 (1): 127-138. 10.1016/j.gene.2004.11.035.PubMedView Article
- Kotlar D, Lavner Y: The action of selection on codon bias in the human genome is related to frequency, complexity, and chronology of amino acids. BMC Genomics. 2006, 7: 67-10.1186/1471-2164-7-67.PubMedPubMed CentralView Article
- Hellmann I, Zollner S, Enard W, Ebersberger I, Nickel B, Paabo S: Selection on human genes as revealed by comparisons to chimpanzee cDNA. Genome Res. 2003, 13 (5): 831-837. 10.1101/gr.944903.PubMedPubMed CentralView Article
- Kondrashov FA, Ogurtsov AY, Kondrashov AS: Selection in favor of nucleotides G and C diversifies evolution rates and levels of polymorphism at mammalian synonymous sites. J Theor Biol. 2006, 240 (4): 616-626. 10.1016/j.jtbi.2005.10.020.PubMedView Article
- Subramanian S: Nearly neutrality and the evolution of codon usage bias in eukaryotic genomes. Genetics. 2008, 178 (4): 2429-2432. 10.1534/genetics.107.086405.PubMedPubMed CentralView Article
- dos Reis M, Wernisch L: Estimating translational selection in eukaryotic genomes. Mol Biol Evol. 2009, 26 (2): 451-461. 10.1093/molbev/msn272.PubMedPubMed CentralView Article
- Akashi H, Gojobori T: Metabolic efficiency and amino acid composition in the proteomes of Escherichia coli and Bacillus subtilis. Proc Natl Acad Sci USA. 2002, 99 (6): 3695-3700. 10.1073/pnas.062526999.PubMedPubMed CentralView Article
- Herbeck JT, Wall DP, Wernegreen JJ: Gene expression level influences amino acid usage, but not codon usage, in the tsetse fly endosymbiont Wigglesworthia. Microbiology. 2003, 149 (Pt 9): 2585-2596. 10.1099/mic.0.26381-0.PubMedView Article
- Koga H, Yuasa S, Nagase T, Shimada K, Nagano M, Imai K, Ohara R, Nakajima D, Murakami M, Kawai M: A comprehensive approach for establishment of the platform to analyze functions of KIAA proteins II: public release of inaugural version of InGaP database containing gene/protein expression profiles for 127 mouse KIAA genes/proteins. DNA Res. 2004, 11 (4): 293-304. 10.1093/dnares/11.4.293.PubMedView Article
- Kikuno R, Nagase T, Nakayama M, Koga H, Okazaki N, Nakajima D, Ohara O: HUGE: a database for human KIAA proteins, a 2004 update integrating HUGEppi and ROUGE. Nucleic Acids Res. 2004, D502-504. 10.1093/nar/gkh035. 32 Database
- Chan PP, Lowe TM: GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res. 2009, D93-97. 10.1093/nar/gkn787. 37 Database
- Misawa K, Kikuno RF: Evaluation of the effect of CpG hypermutability on human codon substitution. Gene. 2009, 431 (1-2): 18-22. 10.1016/j.gene.2008.11.006.PubMedView Article
- Eyre-Walker AC: An analysis of codon usage in mammals: selection or mutation bias?. J Mol Evol. 1991, 33 (5): 442-449. 10.1007/BF02103136.PubMedView Article
- Wang GZ, Chen LL, Zhang HY: Neighboring-site effects of amino acid mutation. Biochem Biophys Res Commun. 2007, 353 (3): 531-534. 10.1016/j.bbrc.2006.12.089.PubMedView Article
- Misawa K, Kamatani N, Kikuno RF: The universal trend of amino acid gain-loss is caused by CpG hypermutability. J Mol Evol. 2008, 67 (4): 334-342. 10.1007/s00239-008-9141-1.PubMedView Article
- R Development Core Team: R: A language and environment for statistical computing. 2008, Vienna, Austria
- Okazaki N, Imai K, Kikuno RF, Misawa K, Kawai M, Inamoto S, Ohara R, Nagase T, Ohara O, Koga H: Influence of the 3'-UTR-length of mKIAA cDNAs and their sequence features to the mRNA expression level in the brain. DNA Res. 2005, 12 (3): 181-189. 10.1093/dnares/dsi001.PubMedView Article
- John AM, Bell JM: Amino acid requirements of the growing mouse. J Nutr. 1976, 106 (9): 1361-1367.PubMed
- Bird AP: DNA methylation and the frequency of CpG in animal DNA. Nucleic Acids Res. 1980, 8 (7): 1499-1504. 10.1093/nar/8.7.1499.PubMedPubMed CentralView Article
- Jukes TH: Codons and nearest-neighbor nucleotide pairs in mammalian messenger RNA. J Mol Evol. 1978, 11 (2): 121-127. 10.1007/BF01733888.PubMedView Article
- Karlin S, Mrazek J: What drives codon choices in human genes?. J Mol Biol. 1996, 262 (4): 459-472. 10.1006/jmbi.1996.0528.PubMedView Article
- Krajewski C, Blacket M, Buckley L, Westerman M: A multigene assessment of phylogenetic relationships within the dasyurid marsupial subfamily Sminthopsinae. Mol Phylogenet Evol. 1997, 8 (2): 236-248. 10.1006/mpev.1997.0421.PubMedView Article
- Huttley GA: Modeling the impact of DNA methylation on the evolution of BRCA1 in mammals. Mol Biol Evol. 2004, 21 (9): 1760-1768. 10.1093/molbev/msh187.PubMedView Article
- Lunter G: Probabilistic whole-genome alignments reveal high indel rates in the human and mouse genomes. Bioinformatics. 2007, 23 (13): i289-296. 10.1093/bioinformatics/btm185.PubMedView Article
- Scarano E, Iaccarino M, Grippo P, Parisi E: The heterogeneity of thymine methyl group origin in DNA pyrimidine isostichs of developing sea urchin embryos. Proc Natl Acad Sci USA. 1967, 57 (5): 1394-1400. 10.1073/pnas.57.5.1394.PubMedPubMed CentralView Article
- Bulmer M: The selection-mutation-drift theory of synonymous codon usage. Genetics. 1991, 129 (3): 897-907.PubMedPubMed Central
- Akashi H: Synonymous codon usage in Drosophila melanogaster: natural selection and translational accuracy. Genetics. 1994, 136 (3): 927-935.PubMedPubMed Central
- Nakamura Y, Gojobori T, Ikemura T: Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res. 2000, 28 (1): 292-10.1093/nar/28.1.292.PubMedPubMed CentralView Article
- Fryxell KJ, Moon WJ: CpG mutation rates in the human genome are highly dependent on local GC content. Mol Biol Evol. 2005, 22 (3): 650-658. 10.1093/molbev/msi043.PubMedView Article
- Taylor J, Tyekucheva S, Zody M, Chiaromonte F, Makova KD: Strong and weak male mutation bias at different sites in the primate genomes: insights from the human-chimpanzee comparison. Mol Biol Evol. 2006, 23 (3): 565-573. 10.1093/molbev/msj060.PubMedView Article
- Tyekucheva S, Makova KD, Karro JE, Hardison RC, Miller W, Chiaromonte F: Human-macaque comparisons illuminate variation in neutral substitution rates. Genome Biol. 2008, 9 (4): R76-10.1186/gb-2008-9-4-r76.PubMedPubMed CentralView Article
- Walser JC, Ponger L, Furano AV: CpG dinucleotides and the mutation rate of non-CpG DNA. Genome Res. 2008, 18 (9): 1403-1414. 10.1101/gr.076455.108.PubMedPubMed CentralView Article
- Bernardi G: The vertebrate genome: isochores and evolution. Mol Biol Evol. 1993, 10 (1): 186-204.PubMed
- Vinogradov AE: Isochores and tissue-specificity. Nucleic Acids Res. 2003, 31 (17): 5212-5220. 10.1093/nar/gkg699.PubMedPubMed CentralView Article
- Vinogradov AE: Noncoding DNA, isochores and gene expression: nucleosome formation potential. Nucleic Acids Res. 2005, 33 (2): 559-563. 10.1093/nar/gki184.PubMedPubMed CentralView Article
- Hsiao LL, Dangond F, Yoshida T, Hong R, Jensen RV, Misra J, Dillon W, Lee KF, Clark KE, Haverty P, et al: A compendium of gene expression in normal human tissues. Physiol Genomics. 2001, 7 (2): 97-104.PubMedView Article
- Semon M, Lobry JR, Duret L: No evidence for tissue-specific adaptation of synonymous codon usage in humans. Mol Biol Evol. 2006, 23 (3): 523-529. 10.1093/molbev/msj053.PubMedView Article
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.