Although the proteins that read the gene regulatory code, transcription factors (TFs), have been largely identified, it is not well known which sequences TFs can recognize. We have analyzed the sequence-specific binding of human TFs using high-throughput SELEX and ChIP sequencing. A total of 830 binding profiles were obtained, describing 239 distinctly different binding specificities. The models represent the majority of human TFs, approximately doubling the coverage compared to existing systematic studies. Our results reveal additional specificity determinants for a large number of factors for which a partial specificity was known, including a commonly observed A- or T-rich stretch that flanks the core motifs. Global analysis of the data revealed that homodimer orientation and spacing preferences, and base-stacking interactions, have a larger role in TF-DNA binding than previously appreciated. We further describe a binding model incorporating these features that is required to understand binding of TFs to DNA.

The majority of CpG dinucleotides in the human genome are methylated at cytosine bases. However, active gene regulatory elements are generally hypomethylated relative to their flanking regions, and the binding of some transcription factors (TFs) is diminished by methylation of their target sequences. By analysis of 542 human TFs with methylation-sensitive SELEX (systematic evolution of ligands by exponential enrichment), we found that there are also many TFs that prefer CpG-methylated sequences. Most of these are in the extended homeodomain family. Structural analysis showed that homeodomain specificity for methylcytosine depends on direct hydrophobic interactions with the methylcytosine 5-methyl group. This study provides a systematic examination of the effect of an epigenetic DNA modification on human TF binding specificity and reveals that many developmentally important proteins display preference for mCpG-containing sequences.

Uterine leiomyomas, or fibroids, are benign tumors that affect millions of women worldwide and that can cause considerable morbidity. To study the genetic basis of this tumor type, we examined 18 uterine leiomyomas derived from 17 different patients by exome sequencing and identified tumor-specific mutations in the mediator complex subunit 12 (MED12) gene in 10. Through analysis of 207 additional tumors, we determined that MED12 is altered in 70% (159 of 225) of tumors from a total of 80 patients. The Mediator complex is a 26-subunit transcriptional regulator that bridges DNA regulatory sequences to the RNA polymerase II initiation complex. All mutations resided in exon 2, suggesting that aberrant function of this region of MED12 contributes to tumorigenesis.