Numerous efforts are underway to determine gene regulatory networks that describe physical relationships between transcription factors (TFs) and their target DNA sequences. Members of paralogous TF families typically recognize similar DNA sequences. Knowledge of the molecular determinants of protein-DNA recognition by paralogous TFs is of central importance for understanding how small differences in DNA specificities can dictate target gene selection. Previously, we determined the in vitro DNA binding specificities of 19 Caenorhabditis elegans basic helix-loop-helix (bHLH) dimers using protein binding microarrays. These TFs bind E-box (CANNTG) and E-box-like sequences. Here, we combine these data with logics, bHLH-DNA co-crystal structures and computational modeling to infer which bHLH monomer can interact with which CAN E-box half-site and we identify a critical residue in the protein that dictates this specificity. Validation experiments using mutant bHLH proteins provide support for our inferences. Our study provides insights into the mechanisms of DNA recognition by bHLH dimers as well as a blueprint for system-level studies of the DNA binding determinants of other TF families in different model organisms and humans.

It has become increasingly clear that transcription factors (TFs) play crucial roles in the development and day-to-day homeostasis that all biological systems experience. TFs target particular genes in a genome, at the appropriate place and time, to regulate their expression so as to elicit the most appropriate biological response from a cell or multicellular organism. TFs can often be grouped into families based on the presence of similar DNA binding domains, and these families are believed to have expanded and diverged throughout evolution by several rounds of gene duplication and mutation. The extent to which TFs within a family have functionally diverged, however, has remained unclear. We propose that systematic analysis of multiple aspects, or parameters, of TF functionality for entire families of TFs could provide clues as to how divergent paralogous TFs really are. We present here a multiparameter integrated network of the activity of the basic helix-loop-helix (bHLH) TFs from the nematode Caenorhabditis elegans. Our data, and the resulting network, indicate that several parameters of bHLH function contribute to their divergence and that many bHLH TFs and their associated parameters exhibit a wide range of connectivity in the network, some being uniquely associated to one another, whereas others are highly connected to multiple parameter associations. We find that 34 bHLH proteins dimerize to form 30 bHLH dimers, which are expressed in a wide range of tissues and cell types, particularly during the development of the nematode. These dimers bind to E-Box DNA sequences and E-Box-like sequences with specificity for nucleotides central to and flanking those E-Boxes and related sequences. Our integrated network is the first such network for a multicellular organism, describing the dimerization specificity, spatiotemporal expression patterns, and DNA binding specificities of an entire family of TFs. The network elucidates the state of bHLH TF divergence in C. elegans with respect to multiple functional parameters and suggests that each bHLH TF, despite many molecular similarities, is distinct from its family members. This functional distinction may indeed explain how TFs from a single family can acquire different biological functions despite descending from common genetic ancestry.

Differences in expression, protein interactions, and DNA binding of paralogous transcription factors ("TF parameters") are thought to be important determinants of regulatory and biological specificity. However, both the extent of TF divergence and the relative contribution of individual TF parameters remain undetermined. We comprehensively identify dimerization partners, spatiotemporal expression patterns, and DNA-binding specificities for the C. elegans bHLH family of TFs, and model these data into an integrated network. This network displays both specificity and promiscuity, as some bHLH proteins, DNA sequences, and tissues are highly connected, whereas others are not. By comparing all bHLH TFs, we find extensive divergence and that all three parameters contribute equally to bHLH divergence. Our approach provides a framework for examining divergence for other protein families in C. elegans and in other complex multicellular organisms, including humans. Cross-species comparisons of integrated networks may provide further insights into molecular features underlying protein family evolution. For a video summary of this article, see the PaperFlick file available with the online Supplemental Data.

Now that numerous high-quality complete genome sequences are available, many efforts are focusing on the "second genomic code", namely the code that determines how the precise temporal and spatial expression of each gene in the genome is achieved. In this regard, the elucidation of transcription regulatory networks that describe combined transcriptional circuits for an organism of interest has become valuable to our understanding of gene expression at a systems level. Such networks describe physical and regulatory interactions between transcription factors (TFs) and the target genes they regulate under different developmental, physiological, or pathological conditions. The mapping of high-quality transcription regulatory networks depends not only on the accuracy of the experimental or computational method chosen, but also relies on the quality of TF predictions. Moreover, the total repertoire of TFs is not only determined by the protein-coding capacity of the genome, but also by different protein properties, including dimerization, co-factor interactions and post-translational modifications. Here, we discuss the factors that influence TF functionality and, hence, the functionality of the networks in which they operate.

Yeast one-hybrid (Y1H) assays provide a gene-centered method for the identification of interactions between gene promoters and regulatory transcription factors (TFs). To date, Y1H assays have involved library screens that are relatively expensive and laborious. We present two Y1H strategies that allow immediate prey identification: matrix assays that use an array of 755 individual Caenorhabditis elegans TFs, and smart-pool assays that use TF multiplexing. Both strategies simplify the Y1H pipeline and reduce the cost of protein-DNA interaction identification. We used a Steiner triple system (STS) to create smart pools of 4-25 TFs. Notably, we uniplexed a small number of highly connected TFs to allow efficient assay deconvolution. Both strategies outperform library screens in terms of coverage, confidence and throughput. These versatile strategies can be adapted both to TFs in other systems and, likely, to other biomolecules and assays as well.

BACKGROUND: The C. elegans Promoterome is a powerful resource for revealing the regulatory mechanisms by which transcription is controlled pan-genomically. Transcription factors will form the core of any systems biology model of genome control and therefore the promoter activity of Promoterome inserts for C. elegans transcription factor genes was examined, in vivo, with a reporter gene approach. RESULTS: Transgenic C. elegans strains were generated for 366 transcription factor promoter/gfp reporter gene fusions. GFP distributions were determined, and then summarized with reference to developmental stage and cell type. Reliability of these data was demonstrated by comparison to previously described gene product distributions. A detailed consideration of the results for one C. elegans transcription factor gene family, the Six family, comprising ceh-32, ceh-33, ceh-34 and unc-39 illustrates the value of these analyses. The high proportion of Promoterome reporter fusions that drove GFP expression, compared to previous studies, led to the hypothesis that transcription factor genes might be involved in local gene duplication events less frequently than other genes. Comparison of transcription factor genes of C. elegans and Caenorhabditis briggsae was therefore carried out and revealed very few examples of functional gene duplication since the divergence of these species for most, but not all, transcription factor gene families. CONCLUSION: Examining reporter expression patterns for hundreds of promoters informs, and thereby improves, interpretation of this data type. Genes encoding transcription factors involved in intrinsic developmental control processes appear acutely sensitive to changes in gene dosage through local gene duplication, on an evolutionary time scale.

Transcription regulatory networks consist of physical and functional interactions between transcription factors (TFs) and their target genes. The systematic mapping of TF-target gene interactions has been pioneered in unicellular systems, using "TF-centered" methods (e.g., chromatin immunoprecipitation). However, metazoan systems are less amenable to such methods. Here, we used "gene-centered" high-throughput yeast one-hybrid (Y1H) assays to identify 283 interactions between 72 C. elegans digestive tract gene promoters and 117 proteins. The resulting protein-DNA interaction (PDI) network is highly connected and enriched for TFs that are expressed in the digestive tract. We provide functional annotations for approximately 10% of all worm TFs, many of which were previously uncharacterized, and find ten novel putative TFs, illustrating the power of a gene-centered approach. We provide additional in vivo evidence for multiple PDIs and illustrate how the PDI network provides insights into metazoan differential gene expression at a systems level. Y1H dataset can be found as a supplemental file to this paper. See Additional Files below. Legend at bottom of Excel spreadsheet.

BACKGROUND: Transcription regulatory networks are composed of interactions between transcription factors and their target genes. Whereas unicellular networks have been studied extensively, metazoan transcription regulatory networks remain largely unexplored. Caenorhabditis elegans provides a powerful model to study such metazoan networks because its genome is completely sequenced and many functional genomic tools are available. While C. elegans gene predictions have undergone continuous refinement, this is not true for the annotation of functional transcription factors. The comprehensive identification of transcription factors is essential for the systematic mapping of transcription regulatory networks because it enables the creation of physical transcription factor resources that can be used in assays to map interactions between transcription factors and their target genes. RESULTS: By computational searches and extensive manual curation, we have identified a compendium of 934 transcription factor genes (referred to as wTF2.0). We find that manual curation drastically reduces the number of both false positive and false negative transcription factor predictions. We discuss how transcription factor splice variants and dimer formation may affect the total number of functional transcription factors. In contrast to mouse transcription factor genes, we find that C. elegans transcription factor genes do not undergo significantly more splicing than other genes. This difference may contribute to differences in organism complexity. We identify candidate redundant worm transcription factor genes and orthologous worm and human transcription factor pairs. Finally, we discuss how wTF2.0 can be used together with physical transcription factor clone resources to facilitate the systematic mapping of C. elegans transcription regulatory networks. CONCLUSION: wTF2.0 provides a starting point to decipher the transcription regulatory networks that control metazoan development and function.