Divergence in DNA specificity among paralogous transcription factors contributes to their differential in vivo binding [RG_Runx1Runx2_v1]

Members of transcription factor (TF) families, i.e. paralogous TFs, are oftentimes reported to have identical DNA-binding motifs, despite the fact that they perform distinct regulatory functions in the cell. Differential genomic targeting by paralogous TFs is generally assumed to be due to interactions with protein cofactors or the chromatin environment. Contrary to previous assumptions, we find that paralogous TFs have different intrinsic preferences for DNA, not captured by current motif models, and these differences partly explain differential genomic binding and functional specificity. Our finding was possible due to a unique combination of carefully designed high-throughput assays and rigorous computation modeling, integrated into a unified framework called iMADS. We used iMADS to quantity, model, and analyze specificity differences between 11 paralogous TFs from 4 distinct human TF families. Our finding of differential specificity between closely related TFs has important implications for the interpretation of the regulatory effects of non-coding genetic variants. Four genomic-context protein binding microarray (gcPBM) experiments were performed for recombinant, full-length, human transcription factors Runx1 and Runx2. Briefly, the PBMs involved binding of GST-tagged transcription factors Runx1 and Runx2 to double-stranded 180K Agilent microarrays in order to determine their binding specificity for putative DNA binding sites in native genomic context. Each genomic DNA sequence represented on the array is present in 6 replicate spots. We report the gcPBM signal intensity for each spot.