{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,6,19]],"date-time":"2026-06-19T07:34:21Z","timestamp":1781854461527,"version":"3.54.5"},"reference-count":53,"publisher":"Oxford University Press (OUP)","issue":"24","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2012,12,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>Motivation: Statistical methods for comparing relative rates of synonymous and non-synonymous substitutions maintain a central role in detecting positive selection. To identify selection, researchers often estimate the ratio of these relative rates () at individual alignment sites. Fitting a codon substitution model that captures heterogeneity in across sites provides a reliable way to perform such estimation, but it remains computationally prohibitive for massive datasets. By using crude estimates of the numbers of synonymous and non-synonymous substitutions at each site, counting approaches scale well to large datasets, but they fail to account for ancestral state reconstruction uncertainty and to provide site-specific estimates.<\/jats:p>\n               <jats:p>Results: We propose a hybrid solution that borrows the computational strength of counting methods, but augments these methods with empirical Bayes modeling to produce a relatively fast and reliable method capable of estimating site-specific values in large datasets. Importantly, our hybrid approach, set in a Bayesian framework, integrates over the posterior distribution of phylogenies and ancestral reconstructions to quantify uncertainty about site-specific estimates. Simulations demonstrate that this method competes well with more-principled statistical procedures and, in some cases, even outperforms them. We illustrate the utility of our method using human immunodeficiency virus, feline panleukopenia and canine parvovirus evolution examples.<\/jats:p>\n               <jats:p>Availability: Renaissance counting is implemented in the development branch of BEAST, freely available at http:\/\/code.google.com\/p\/beast-mcmc\/. The method will be made available in the next public release of the package, including support to set up analyses in BEAUti.<\/jats:p>\n               <jats:p>Contact: \u00a0philippe.lemey@rega.kuleuven.be or msuchard@ucla.edu<\/jats:p>\n               <jats:p>Supplementary information: \u00a0Supplementary data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/bts580","type":"journal-article","created":{"date-parts":[[2012,10,13]],"date-time":"2012-10-13T04:11:30Z","timestamp":1350101490000},"page":"3248-3256","source":"Crossref","is-referenced-by-count":89,"title":["A counting renaissance: combining stochastic mapping and empirical Bayes to quickly detect amino acid sites under positive selection"],"prefix":"10.1093","volume":"28","author":[{"given":"Philippe","family":"Lemey","sequence":"first","affiliation":[{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Vladimir N.","family":"Minin","sequence":"additional","affiliation":[{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Filip","family":"Bielejec","sequence":"additional","affiliation":[{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Sergei L.","family":"Kosakovsky Pond","sequence":"additional","affiliation":[{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Marc A.","family":"Suchard","sequence":"additional","affiliation":[{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"},{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"},{"name":"1 Department of Microbiology and Immunology, Rega Institute, KU Leuven, B-3000 Leuven, Belgium, 2Department of Statistics, University of Washington, Seattle, WA 98195, 3Department of Medicine, University of California, San Diego, CA 92103, 4Department of Biomathematics, 5Department of Human Genetics, David Geffen School of Medicine and 6Department of Biostatistics, School of Public Health, University of California, Los Angeles, CA 90095, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2012,10,12]]},"reference":[{"key":"2023012513254008700_bts580-B1","doi-asserted-by":"crossref","first-page":"170","DOI":"10.1093\/sysbio\/syr100","article-title":"BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics","volume":"61","author":"Ayres","year":"2012","journal-title":"Syst. Biol."},{"key":"2023012513254008700_bts580-B2","doi-asserted-by":"crossref","first-page":"447","DOI":"10.1097\/QAD.0b013e3280129964","article-title":"Identification of accessory mutations associated with high-level resistance in HIV-1 reverse transcriptase","volume":"21","author":"Cane","year":"2007","journal-title":"AIDS"},{"key":"2023012513254008700_bts580-B3","first-page":"157","article-title":"Polymorphisms in the viral reverse transcriptase predict the evolution towards distinct thymidine analogue mutational patterns: a longitudinal analysis","volume":"11","author":"De Luca","year":"2006","journal-title":"Antivir. Ther."},{"key":"2023012513254008700_bts580-B4","doi-asserted-by":"crossref","first-page":"1969","DOI":"10.1093\/molbev\/mss075","article-title":"Bayesian Phylogenetics with BEAUti and the BEAST 1.7","volume":"29","author":"Drummond","year":"2012","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B5","doi-asserted-by":"crossref","first-page":"265","DOI":"10.1089\/aid.1991.7.265","article-title":"In vivo prevalence of azidothymidine (AZT) resistance mutations in an aids patient before and after AZT therapy","volume":"7","author":"Fitzgibbon","year":"1991","journal-title":"AIDS Res. Hum. Retroviruses"},{"key":"2023012513254008700_bts580-B6","first-page":"733","article-title":"Posterior predictive assessment of model fitness via realized discrepancies","volume":"6","author":"Gelman","year":"1996","journal-title":"Stat. Sin."},{"key":"2023012513254008700_bts580-B7","first-page":"725","article-title":"A codon-based model of nucleotide substitution for protein-coding DNA sequences","volume":"11","author":"Goldman","year":"1994","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B8","doi-asserted-by":"crossref","first-page":"995","DOI":"10.7326\/0003-4819-130-12-199906150-00008","article-title":"Toward evidence-based medical statistics. 1: The P value fallacy","volume":"130","author":"Goodman","year":"1999","journal-title":"Ann. Intern. Med."},{"key":"2023012513254008700_bts580-B9","doi-asserted-by":"crossref","first-page":"327","DOI":"10.1126\/science.1090727","article-title":"Unifying the epidemiological and evolutionary dynamics of pathogens","volume":"303","author":"Grenfell","year":"2004","journal-title":"Science"},{"key":"2023012513254008700_bts580-B10","doi-asserted-by":"crossref","first-page":"160","DOI":"10.1007\/BF02101694","article-title":"Dating the human-ape splitting by a molecular clock of mitochondrial DNA","volume":"22","author":"Hasegawa","year":"1985","journal-title":"J. Mol. Evol."},{"key":"2023012513254008700_bts580-B11","doi-asserted-by":"crossref","first-page":"1157","DOI":"10.1198\/016214505000001393","article-title":"Post-processing posterior predictive p values","volume":"101","author":"Hjort","year":"2006","journal-title":"J Am. Stat. Assoc."},{"key":"2023012513254008700_bts580-B12","doi-asserted-by":"crossref","first-page":"1204","DOI":"10.1214\/09-AOAS247","article-title":"Simulation from endpoint-conditioned, continuous-time Markov chains on a finite state space, with applications to molecular evolution","volume":"3","author":"Hobolth","year":"2009","journal-title":"Ann. Appl. Stat."},{"key":"2023012513254008700_bts580-B13","doi-asserted-by":"crossref","first-page":"661","DOI":"10.1007\/s00239-004-2588-9","article-title":"Bayesian estimation of positively selected sites","volume":"58","author":"Huelsenbeck","year":"2004","journal-title":"J. Mol. Evol."},{"key":"2023012513254008700_bts580-B14","doi-asserted-by":"crossref","first-page":"167","DOI":"10.1038\/335167a0","article-title":"Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection","volume":"335","author":"Hughes","year":"1988","journal-title":"Nature"},{"key":"2023012513254008700_bts580-B15","first-page":"113","article-title":"Compensatory fixation explains long term persistence of the m41l in HIV-1 reverse transcriptase in a large transmission cluster","volume":"11","author":"Huigen","year":"2006","journal-title":"Antivir. Ther."},{"key":"2023012513254008700_bts580-B16","doi-asserted-by":"crossref","first-page":"e225","DOI":"10.1371\/journal.pone.0000225","article-title":"Sequential broadening of CTL responses in early HIV-1 infection is associated with viral escape","volume":"2","author":"Karlsson","year":"2007","journal-title":"PLoS One"},{"key":"2023012513254008700_bts580-B17","doi-asserted-by":"crossref","first-page":"5093","DOI":"10.1128\/JVI.72.6.5093-5098.1998","article-title":"A novel polymorphism at codon 333 of human immunodeficiency virus type 1 reverse transcriptase can facilitate dual resistance to zidovudine and L-2\u2019,3\u2019-dideoxy-3\u2019-thiacytidine","volume":"72","author":"Kemp","year":"1998","journal-title":"J. Virol."},{"key":"2023012513254008700_bts580-B18","doi-asserted-by":"crossref","first-page":"1208","DOI":"10.1093\/molbev\/msi105","article-title":"Not so different after all: a comparison of methods for detecting amino acid sites under selection","volume":"22","author":"Kosakovsky Pond","year":"2005","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B19","doi-asserted-by":"crossref","first-page":"676","DOI":"10.1093\/bioinformatics\/bti079","article-title":"HyPhy: hypothesis testing using phylogenies","volume":"21","author":"Kosakovsky Pond","year":"2005","journal-title":"Bioinformatics"},{"key":"2023012513254008700_bts580-B20","doi-asserted-by":"crossref","first-page":"e1000304","DOI":"10.1371\/journal.pgen.1000304","article-title":"The population genetics of dN\/dS","volume":"4","author":"Kryazhimskiy","year":"2008","journal-title":"PLoS Genet."},{"key":"2023012513254008700_bts580-B21","doi-asserted-by":"crossref","first-page":"1155","DOI":"10.1126\/science.2479983","article-title":"Multiple mutations in HIV-1 reverse transcriptase confer high-level resistance to zidovudine (AZT)","volume":"246","author":"Larder","year":"1989","journal-title":"Science"},{"key":"2023012513254008700_bts580-B22","doi-asserted-by":"crossref","first-page":"1701","DOI":"10.1089\/cmb.2006.13.1701","article-title":"Conjugate Gibbs sampling for Bayesian phylogenetic models","volume":"13","author":"Lartillot","year":"2006","journal-title":"J. Comput. Biol."},{"key":"2023012513254008700_bts580-B23","doi-asserted-by":"crossref","first-page":"e29","DOI":"10.1371\/journal.pcbi.0030029","article-title":"Synonymous substitution rates predict HIV disease progression as a result of underlying replication dynamics","volume":"3","author":"Lemey","year":"2007","journal-title":"PLoS Comput. Biol."},{"key":"2023012513254008700_bts580-B24","doi-asserted-by":"crossref","first-page":"e1000520","DOI":"10.1371\/journal.pcbi.1000520","article-title":"Bayesian phylogeography finds its root","volume":"5","author":"Lemey","year":"2009","journal-title":"PLoS Comput. Biol."},{"key":"2023012513254008700_bts580-B25","doi-asserted-by":"crossref","first-page":"1877","DOI":"10.1093\/molbev\/msq067","article-title":"Phylogeography takes a relaxed random walk in continuous space and time","volume":"27","author":"Lemey","year":"2010","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B26","doi-asserted-by":"crossref","first-page":"20","DOI":"10.1097\/01.qai.0000172369.82456.36","article-title":"Effect of the Q207D mutation in HIV type 1 reverse transcriptase on zidovudine susceptibility and replicative fitness","volume":"40","author":"Lu","year":"2005","journal-title":"J. Acquir. Immune Defic. Syndr."},{"key":"2023012513254008700_bts580-B27","doi-asserted-by":"crossref","first-page":"349","DOI":"10.1093\/biomet\/56.2.349","article-title":"Empirical Bayes estimation for the Poisson distribution","volume":"56","author":"Maritz","year":"1969","journal-title":"Biometrika"},{"key":"2023012513254008700_bts580-B28","doi-asserted-by":"crossref","first-page":"151","DOI":"10.1038\/385151a0","article-title":"Episodic adaptive evolution of primate lysozymes","volume":"385","author":"Messier","year":"1997","journal-title":"Nature"},{"key":"2023012513254008700_bts580-B29","doi-asserted-by":"crossref","first-page":"391","DOI":"10.1007\/s00285-007-0120-8","article-title":"Counting labeled transitions in continuous-time Markov models of evolution","volume":"56","author":"Minin","year":"2008","journal-title":"J. Math. Biol."},{"key":"2023012513254008700_bts580-B30","doi-asserted-by":"crossref","first-page":"17","DOI":"10.2202\/1544-6115.1650","article-title":"Imputation estimators partially correct for model misspecification","volume":"10","author":"Minin","year":"2011","journal-title":"Stat. Appl. Genet. Mol. Biol."},{"key":"2023012513254008700_bts580-B31","doi-asserted-by":"crossref","first-page":"3985","DOI":"10.1098\/rstb.2008.0176","article-title":"Fast, accurate and simulation-free stochastic mapping of discrete traits","volume":"363","author":"Minin","year":"2008","journal-title":"Philos. Trans. R. Soc. Lond. B: Biol. Sci."},{"key":"2023012513254008700_bts580-B32","first-page":"715","article-title":"A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome","volume":"11","author":"Muse","year":"1994","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B33","doi-asserted-by":"crossref","first-page":"729","DOI":"10.1080\/10635150290102393","article-title":"Mapping mutations on phylogenies","volume":"51","author":"Nielsen","year":"2002","journal-title":"Syst. Biol."},{"key":"2023012513254008700_bts580-B34","doi-asserted-by":"crossref","first-page":"929","DOI":"10.1093\/genetics\/148.3.929","article-title":"Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene","volume":"148","author":"Nielsen","year":"1998","journal-title":"Genetics"},{"key":"2023012513254008700_bts580-B35","doi-asserted-by":"crossref","first-page":"801","DOI":"10.1093\/molbev\/msp003","article-title":"Learning to count: robust estimates for labeled distances between molecular sequences","volume":"26","author":"O\u2019Brien","year":"2009","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B36","doi-asserted-by":"crossref","first-page":"2375","DOI":"10.1093\/molbev\/msi232","article-title":"Site-to-site variation of synonymous substitution rates","volume":"22","author":"Pond","year":"2005","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B37","doi-asserted-by":"crossref","first-page":"e62","DOI":"10.1371\/journal.pcbi.0020062","article-title":"Adaptation to different human populations by HIV-1 revealed by codon-based analyses","volume":"2","author":"Pond","year":"2006","journal-title":"PLoS Comput. Biol."},{"key":"2023012513254008700_bts580-B38","doi-asserted-by":"crossref","first-page":"540","DOI":"10.1038\/nrg2583","article-title":"Evolutionary analysis of the dynamics of viral infectious disease","volume":"10","author":"Pybus","year":"2009","journal-title":"Nat. Rev. Genet."},{"key":"2023012513254008700_bts580-B39","first-page":"157","article-title":"An empirical Bayes approach to statistics","volume-title":"Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability, Volume 1: Contributions to the Theory of Statistics","author":"Robbins","year":"1956"},{"key":"2023012513254008700_bts580-B40","doi-asserted-by":"crossref","first-page":"56","DOI":"10.1093\/bioinformatics\/btm532","article-title":"Uniformization for sampling realizations of Markov processes: applications to Bayesian implementations of codon substitution models","volume":"24","author":"Rodrigue","year":"2008","journal-title":"Bioinformatics"},{"key":"2023012513254008700_bts580-B41","doi-asserted-by":"crossref","first-page":"9","DOI":"10.1002\/jmv.20500","article-title":"Impact of unreported HIV-1 reverse transcriptase mutations on phenotypic resistance to nucleoside and non-nucleoside inhibitors","volume":"78","author":"Saracino","year":"2006","journal-title":"J. Med. Virol."},{"key":"2023012513254008700_bts580-B42","doi-asserted-by":"crossref","first-page":"379","DOI":"10.1073\/pnas.0406765102","article-title":"High rate of viral evolution associated with the emergence of carnivore parvovirus","volume":"102","author":"Shackelton","year":"2005","journal-title":"Proc. Natl Acad. Sci. USA"},{"key":"2023012513254008700_bts580-B43","first-page":"190","article-title":"New methods for detecting lineage-specific selection","volume-title":"Proceedings of the 10th International Conference on Research in Computational Molecular Biology","author":"Siepel","year":"2006"},{"key":"2023012513254008700_bts580-B44","doi-asserted-by":"crossref","first-page":"4000","DOI":"10.1128\/AAC.46.12.4000-4003.2002","article-title":"Phenotypic and genotypic analysis of biologically cloned human immunodeficiency virus type 1 isolates from patients treated with zidovudine and lamivudine","volume":"46","author":"Stoeckli","year":"2002","journal-title":"Antimicrob. Agents Chemother."},{"key":"2023012513254008700_bts580-B45","doi-asserted-by":"crossref","first-page":"1370","DOI":"10.1093\/bioinformatics\/btp244","article-title":"Many-core algorithms for statistical phylogenetics","volume":"25","author":"Suchard","year":"2009","journal-title":"Bioinformatics"},{"key":"2023012513254008700_bts580-B46","doi-asserted-by":"crossref","first-page":"1315","DOI":"10.1093\/oxfordjournals.molbev.a026042","article-title":"A method for detecting positive selection at single amino acid sites","volume":"16","author":"Suzuki","year":"1999","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B47","doi-asserted-by":"crossref","first-page":"7186","DOI":"10.1128\/JVI.02084-05","article-title":"Involvement of novel human immunodeficiency virus type 1 reverse transcriptase mutations in the regulation of resistance to nucleoside inhibitors","volume":"80","author":"Svicher","year":"2006","journal-title":"J. Virol."},{"key":"2023012513254008700_bts580-B48","doi-asserted-by":"crossref","first-page":"2276","DOI":"10.1128\/AAC.45.8.2276-2279.2001","article-title":"Variants other than aspartic acid at codon 69 of the human immunodeficiency virus type 1 reverse transcriptase gene affect susceptibility to nuleoside analogs","volume":"45","author":"Winters","year":"2001","journal-title":"Antimicrob. Agents Chemother."},{"key":"2023012513254008700_bts580-B49","doi-asserted-by":"crossref","first-page":"496","DOI":"10.1016\/S0169-5347(00)01994-7","article-title":"Statistical methods for detecting molecular adaptation","volume":"15","author":"Yang","year":"2000","journal-title":"Trends Ecol. Evol."},{"key":"2023012513254008700_bts580-B50","doi-asserted-by":"crossref","first-page":"587","DOI":"10.1007\/BF02352289","article-title":"Maximum-likelihood models for combined analyses of multiple sequence data","volume":"42","author":"Yang","year":"1996","journal-title":"J. Mol. Evol."},{"key":"2023012513254008700_bts580-B51","first-page":"555","article-title":"PAML: a program package for phylogenetic analysis by maximum likelihood","volume":"13","author":"Yang","year":"1997","journal-title":"Comput. Appl. Biosc."},{"key":"2023012513254008700_bts580-B52","doi-asserted-by":"crossref","first-page":"1107","DOI":"10.1093\/molbev\/msi097","article-title":"Bayes empirical Bayes inference of amino acid sites under positive selection","volume":"22","author":"Yang","year":"2005","journal-title":"Mol. Biol. Evol."},{"key":"2023012513254008700_bts580-B53","doi-asserted-by":"crossref","first-page":"340","DOI":"10.1007\/s00239-007-9019-7","article-title":"Exploring variation in the","volume":"65","author":"Zhai","year":"2007","journal-title":"J. Mol. Evol."}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/28\/24\/3248\/48883263\/bioinformatics_28_24_3248.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/28\/24\/3248\/48883263\/bioinformatics_28_24_3248.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,1,25]],"date-time":"2023-01-25T19:23:03Z","timestamp":1674674583000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/28\/24\/3248\/243893"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2012,10,12]]},"references-count":53,"journal-issue":{"issue":"24","published-print":{"date-parts":[[2012,12,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bts580","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2012,12]]},"published":{"date-parts":[[2012,10,12]]}}}