Detection of nonneutral substitution rates on mammalian phylogenies.
Methods for detecting nucleotide substitution rates that are faster or slower than expected under neutral drift are widely used to identify candidate functional elements in genomic sequences. However, most existing methods consider either reductions (conservation) or increases (acceleration) in rate but not both, or assume that selection acts uniformly across the branches of a phylogeny. Here we examine the more general problem of detecting departures from the neutral rate of substitution in either direction, possibly in a clade-specific manner. We consider four statistical, phylogenetic tests for addressing this problem: a likelihood ratio test, a score test, a test based on exact distributions of numbers of substitutions, and the genomic evolutionary rate profiling (GERP) test. All four tests have been implemented in a freely available program called phyloP. Based on extensive simulation experiments, these tests are remarkably similar in statistical power. With 36 mammalian species, they all appear to be capable of fairly good sensitivity with low false-positive rates in detecting strong selection at individual nucleotides, moderate selection in 3-bp elements, and weaker or clade-specific selection in longer elements. By applying phyloP to mammalian multiple alignments from the ENCODE project, we shed light on patterns of conservation/acceleration in known and predicted functional elements, approximate fractions of sites subject to constraint, and differences in clade-specific selection in the primate and glires clades. We also describe new "Conservation" tracks in the UCSC Genome Browser that display both phyloP and phastCons scores for genome-wide alignments of 44 vertebrate species.
Pollard
Katherine S
KS
Gladstone Institutes, University of California, San Francisco, San Francisco, California 94158, USA. kpollard@gladstone.ucsf.edu
Hubisz
Melissa J
MJ
Rosenbloom
Kate R
KR
Siepel
Adam
A
D000818:Animals; D001483:Base Sequence; D003198:Computer Simulation; D017124:Conserved Sequence; D019143:Evolution, Molecular; D006801:Humans; D016013:Likelihood Functions; D008322:Mammals / Q000145:classification / Q000235:genetics*; D008957:Models, Genetic; D015233:Models, Statistical; D010802:Phylogeny; D011323:Primates / Q000235:genetics; D012641:Selection, Genetic; D016415:Sequence Alignment; D012984:Software; D013045:Species Specificity