  Genes are often regulated in living cells by proteins called transcription
factors (TFs) that bind directly to short segments of DNA in close proximity to
specific genes. These binding sites have a conserved nucleotide appearance,
which is called a motif. Several recent studies of transcriptional regulation
require the reduction of a large collection of motifs into clusters based on
the similarity of their nucleotide composition. We present a principled
approach to this clustering problem based upon a Bayesian hierarchical model
that accounts for both within- and between-motif variability. We use a
Dirichlet process prior distribution that allows the number of clusters to vary
and we also present a novel generalization that allows the core width of each
motif to vary. This clustering model is implemented, using a Gibbs sampling
strategy, on several collections of transcription factor motif matrices. Our
clusters provide a means by which to organize transcription factors based on
binding motif similarities, which can be used to reduce motif redundancy within
large databases such as JASPAR and TRANSFAC. Finally, our clustering procedure
has been used in combination with discovery of evolutionarily-conserved motifs
to predict co-regulated genes. An alternative to our Dirichlet process prior
distribution is explored but shows no substantive difference in the clustering
results for our datasets. Our Bayesian clustering model based on the Dirichlet
process has several advantages over traditional clustering methods that could
make our procedure appropriate and useful for many clustering applications.
