Genome Biology and Evolution Advance Access published online on May 25, 2009
Genome Biology and Evolution, doi:10.1093/gbe/evp010
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Measuring Transcription Factor Binding Site Turnover: A Maximum Likelihood Approach using Phylogenies
1 Lehrstuhl für Bioinformatik, Institut für Informatik, Universität Leipzig, Härtelstraße 16-18, D-04107 Leipzig, Germany
2 Department of Ecology and Evolutionary Biology, Yale University, 165 Prospect Street, New Haven, CT 06520-8106, USA
* Author for Correspondence: Günter P. Wagner, Department of Ecology and Evolutionary Biology, Yale University , 165 Prospect Street, New Haven, CT 06520-8106, USA, Tel: +203-432-9998, Fax: +203-432-3870, Email: gunter.wagner{at}yale.edu
| Abstract |
|---|
A major mode of gene expression evolution is based on changes in cis-regulatory elements (CREs) whose function critically depends on the presence of transcription factor binding sites (TFBS). Since CREs experience extensive TFBS turnover even with conserved function, alignment-based studies of CRE sequence evolution are limited to very closely related species. Here, we propose an alternative approach based on a stochastic model of TFBS turnover. We implemented a maximum likelihood model that permits variable turnover rates in different parts of the species tree. This model can be used to detect changes in turnover rate as a proxy for differences in the selective pressures acting on TFBS in different clades. We applied this method to five TFBS in the fungi methionine biosynthesis pathway and three TFBS in the HoxA clusters of vertebrates. We find that the estimated turnover rate is generally high, with half-life ranging between
5 million and 150 million years, with a mode around 10s of millions of years. This rate is consistent with the finding that even functionally conserved enhancers can show very low sequence similarity. We also detect statistically significant differences in the equilibrium densities of estrogen- and progesterone response elements in the HoxA clusters between mammal and non-mammal vertebrates. Even more extreme clade-specific differences were found in the fungal data. We conclude that stochastic models of transcription factor binding site turnover enable the detection of shifts in the selective pressures acting on CREs in different organisms.
The analysis tool, called CRETO (Cis-Regulatory Element Turn-Over) can be downloaded from http://www.bioinf.uni-leipzig.de/Software/creto/.
Keywords: cis-regulatory evolution, non-coding sequences, evolution of gene regulation, enhancer evolution, promoter evolution, evolution of development
Received April 3, 2009; Accepted May 15, 2009