Prediction of similarly acting cis-regulatory modules by subsequence profiling and comparative genomics in Drosophila melanogaster and D.pseudoobscura.
BACKGROUND
To date, computational searches for cis-regulatory modules (CRMs) have relied on two methods. The first, phylogenetic footprinting, has been used to find CRMs in non-coding sequence, but does not directly link DNA sequence with spatio-temporal patterns of expression. The second, based on searches for combinations of transcription factor (TF) binding motifs, has been employed in genome-wide discovery of similarly acting enhancers, but requires prior knowledge of the set of TFs acting at the CRM and the TFs' binding motifs.


RESULTS
We propose a method for CRM discovery that combines aspects of both approaches in an effort to overcome their individual limitations. By treating phylogenetically footprinted non-coding regions (PFRs) as proxies for CRMs, we endeavor to find PFRs near co-regulated genes that are comprised of similar short, conserved sequences. Using Markov chains as a convenient formulation to assess similarity, we develop a sampling algorithm to search a large group of PFRs for the most similar subset. When starting with a set of genes involved in Drosophila early blastoderm development and using phylogenetic comparisons of Drosophila melanogaster and D.pseudoobscura genomes, we show here that our algorithm successfully detects known CRMs. Further, we use our similarity metric, based on Markov chain discrimination, in a genome-wide search, and uncover additional known and many candidate early blastoderm CRMs.


BACKGROUND
Software is available via http://arep.med.harvard.edu/enhancer
Grad
Yonatan H
YH
The Lipper Center for Computational Genetics, Department of Genetics, Harvard Medical School, 77 Avenue Louis Pasteur, Boston, Massachusetts, 02115, USA.
Roth
Frederick P
FP
Halfon
Marc S
MS
Church
George M
GM
D000465:Algorithms; D000818:Animals; D002874:Chromosome Mapping / Q000379:methods*; D016384:Consensus Sequence; D018983:DNA Footprinting / Q000379:methods*; D004330:Drosophila / Q000145:classification / Q000235:genetics*; D020869:Gene Expression Profiling / Q000379:methods*; D018507:Gene Expression Regulation, Developmental / Q000235:genetics; D005809:Genes, Regulator / Q000235:genetics; D010802:Phylogeny; D012045:Regulatory Sequences, Nucleic Acid / Q000235:genetics*; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods; D012689:Sequence Homology, Nucleic Acid; D012984:Software; D013045:Species Specificity