title: The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC.
abstract: Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio ω. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the ω ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to ∼80%) and low GC (∼30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.
authors.affiliation: Department of Ecology and Evolution, Biophore, Lausanne University, Lausanne, Switzerland.
authors.forename: Walid H
authors.initials: WH
authors.lastname: Gharib
authors.forename: Marc
authors.initials: M
authors.lastname: Robinson-Rechavi
mesh_terms: D019943:Amino Acid Substitution / Q000235:genetics*; D001482:Base Composition / Q000235:genetics; D003062:Codon / Q000235:genetics*; D003198:Computer Simulation; D019143:Evolution, Molecular; D008957:Models, Genetic; D016366:Open Reading Frames / Q000235:genetics*; D010802:Phylogeny; D012641:Selection, Genetic
chemical_list: D003062:Codon
keywords: adaptive evolution; base composition; codon model
references.citation: Genome Res. 2008 Sep;18(9):1393-402
references.pmid: 18562677
references.citation: Genetics. 2001 Oct;159(2):907-11
references.pmid: 11693127
references.citation: Annu Rev Genomics Hum Genet. 2009;10:285-311
references.pmid: 19630562
references.citation: Nucleic Acids Res. 2002 Jan 1;30(1):38-41
references.pmid: 11752248
references.citation: Mol Biol Evol. 2007 Jul;24(7):1464-79
references.pmid: 17400572
references.citation: Mol Biol Evol. 2011 Nov;28(11):3033-43
references.pmid: 21670087
references.citation: Gene. 2000 Jan 4;241(1):3-17
references.pmid: 10607893
references.citation: Mol Biol Evol. 2012 Jul;29(7):1861-74
references.pmid: 22319139
references.citation: Philos Trans R Soc Lond B Biol Sci. 2010 Aug 27;365(1552):2571-80
references.pmid: 20643747
references.citation: Genome Biol Evol. 2011;3:974-84
references.pmid: 21795750
references.citation: PLoS Genet. 2010 Sep 09;6(9):e1001107
references.pmid: 20838593
references.citation: Immunogenetics. 2008 Sep;60(9):495-506
references.pmid: 18581108
references.citation: Heredity (Edinb). 2012 Apr;108(4):347-53
references.pmid: 22045380
references.citation: Mol Biol Evol. 2012 Oct;29(10):2889-93
references.pmid: 22490825
references.citation: Mol Biol Evol. 2013 Jan;30(1):177-90
references.pmid: 22936721
references.citation: Mol Biol Evol. 2001 Aug;18(8):1585-92
references.pmid: 11470850
references.citation: Genome Biol Evol. 2009 Jun 05;1:114-8
references.pmid: 20333182
references.citation: Mol Biol Evol. 2000 Jan;17(1):32-43
references.pmid: 10666704
references.citation: Mol Biol Evol. 2010 Oct;27(10):2257-67
references.pmid: 20447933
references.citation: Mol Biol Evol. 2009 Feb;26(2):255-71
references.pmid: 18922761
references.citation: Mol Biol Evol. 2011 Dec;28(12):3297-308
references.pmid: 21690564
references.citation: Genetics. 2003 Jul;164(3):1229-36
references.pmid: 12871927
references.citation: Mol Biol Evol. 2002 Jun;19(6):950-8
references.pmid: 12032251
references.citation: Mol Biol Evol. 2012 Jan;29(1):1-5
references.pmid: 21772063
references.citation: Proc Natl Acad Sci U S A. 2009 Apr 21;106(16):6700-5
references.pmid: 19339501
references.citation: Mol Biol Evol. 2007 Aug;24(8):1586-91
references.pmid: 17483113
references.citation: Science. 1985 May 24;228(4702):953-8
references.pmid: 4001930
references.citation: Mol Biol Evol. 2005 Apr;22(4):1107-18
references.pmid: 15689528
references.citation: Mol Biol Evol. 2007 May;24(5):1219-28
references.pmid: 17339634
references.citation: Genetics. 2000 May;155(1):431-49
references.pmid: 10790415
references.citation: Science. 2008 Jan 25;319(5862):473-6
references.pmid: 18218900
references.citation: PLoS Genet. 2008 Aug 01;4(8):e1000144
references.pmid: 18670650
references.citation: Syst Biol. 1998 Mar;47(1):125-33
references.pmid: 12064232
references.citation: Genome Res. 2009 May;19(5):922-33
references.pmid: 19411606
references.citation: Mol Biol Evol. 2005 Dec;22(12):2472-9
references.pmid: 16107592
references.citation: Mol Biol Evol. 2012 Apr;29(4):1125-39
references.pmid: 22049066
references.citation: Proc Natl Acad Sci U S A. 2003 Aug 5;100(16):9440-5
references.pmid: 12883005
references.citation: Mol Phylogenet Evol. 2007 Feb;42(2):388-93
references.pmid: 16959502
references.citation: Nat Rev Genet. 2001 Jul;2(7):549-55
references.pmid: 11433361
references.citation: Heredity (Edinb). 2007 Oct;99(4):364-73
references.pmid: 17622265
references.citation: Mol Biol Evol. 2011 Mar;28(3):1217-28
references.pmid: 21087944
references.citation: BMC Bioinformatics. 2009 Jun 16;10 Suppl 6:S3
references.pmid: 19534752
references.citation: Mol Biol Evol. 2012 Mar;29(3):1047-57
references.pmid: 22075116