Analysis of homeodomain specificities allows the family-wide prediction of preferred recognition sites

Exploring the DNA-recognition potential of homeodomains

The recognition potential of most families of DNA-binding domains (DBDs) remains relatively unexplored. Homeodomains (HDs), like many other families of DBDs, display limited diversity in their preferred recognition sequences. To explore the recognition potential of HDs, we utilized a bacterial selection system to isolate HD variants, from a randomized library, that are compatible with each of the 64 possible 3' triplet sites (i.e., TAANNN). The majority of these selections yielded sets of HDs with overrepresented residues at specific recognition positions, implying the selection of specific binders. The DNA-binding specificity of 151 representative HD variants was subsequently characterized, identifying HDs that preferentially recognize 44 of these target sites. Many of these variants contain novel combinations of specificity determinants that are uncommon or absent in extant HDs. These novel determinants, when grafted into different HD backbones, produce a corresponding alteration in specificity. This information was used to create more explicit HD recognition models, which can inform the prediction of transcriptional regulatory networks for extant HDs or the engineering of HDs with novel DNA-recognition potential. The diversity of recovered HD recognition sequences raises important questions about the fitness barrier that restricts the evolution of alternate recognition modalities in natural systems.

Emx2 is a dose-dependent negative regulator of Sox2 telencephalic enhancers

The transcription factor Sox2 is essential for neural stem cells (NSC) maintenance in the hippocampus and in vitro. The transcription factor Emx2 is also critical for hippocampal development and NSC self-renewal. Searching for ‘modifier’ genes affecting the Sox2 deficiency phenotype in mouse, we observed that loss of one Emx2 allele substantially increased the telencephalic β-geo (LacZ) expression of a transgene driven by the 5′ or 3′ Sox2 enhancer. Reciprocally, Emx2 overexpression in NSC cultures inhibited the activity of the same transgene. In vivo, loss of one Emx2 allele increased Sox2 levels in the medial telencephalic wall, including the hippocampal primordium. In hypomorphic Sox2 mutants, retaining a single ‘weak’ Sox2 allele, Emx2 deficiency substantially rescued hippocampal radial glia stem cells and neurogenesis, indicating that Emx2 functionally interacts with Sox2 at the stem cell level. Electrophoresis mobility shift assays and transfection indicated that Emx2 represses the activities of both Sox2 enhancers. Emx2 bound to overlapping Emx2/POU-binding sites, preventing binding of the POU transcriptional activator Brn2. Additionally, Emx2 directly interacted with Brn2 without binding to DNA. These data imply that Emx2 may perform part of its functions by negatively modulating Sox2 in specific brain areas, thus controlling important aspects of NSC function in development.

Integration of an abdominal Hox complex with Pax2 yields cell-specific EGF secretion from Drosophila sensory precursor cells

Cis-regulatory modules (CRMs) ensure specific developmental outcomes by mediating both proper spatiotemporal gene expression patterns and appropriate transcriptional levels. In Drosophila, the precise transcriptional control of the serine protease rhomboid regulates EGF signaling to specify distinct cell types. Recently, we identified a CRM that activates rhomboid expression and thereby EGF secretion from a subset of abdominal sensory organ precursor cells (SOPs) to induce an appropriate number of lipid-processing cells called oenocytes. Here, we use scanning mutagenesis coupled with reporter assays, biochemistry and genetics to dissect the transcriptional mechanisms regulating SOP-specific rhomboid activation. Our results show that proper spatial activity of the rhomboid CRM is dependent upon direct integration of the abdomen-specific Hox factor Abdominal-A and the SOP-restricted Pax2 factor. In addition, we show that the Extradenticle and Homothorax Hox co-factors are differentially integrated on the rhomboid CRM by abdominal versus thoracic Hox proteins in the presence of Pax2. Last, we show that Abdominal-A uses both Pax2-dependent and Pax2-independent mechanisms to stimulate rhomboid CRM activity to induce proper oenocyte numbers. Thus, these data demonstrate how a CRM integrates Hox and neural transcriptional inputs to regulate the appropriate spatial pattern and levels of EGF secretion to specify an essential cell fate.

Antagonistic regulation of apoptosis and differentiation by the Cut transcription factor represents a tumor-suppressing mechanism in Drosophila

Apoptosis is essential to prevent oncogenic transformation by triggering self-destruction of harmful cells, including those unable to differentiate. However, the mechanisms linking impaired cell differentiation and apoptosis during development and disease are not well understood. Here we report that the Drosophila transcription factor Cut coordinately controls differentiation and repression of apoptosis via direct regulation of the pro-apoptotic gene reaper. We also demonstrate that this regulatory circuit acts in diverse cell lineages to remove uncommitted precursor cells in status nascendi and thereby interferes with their potential to develop into cancer cells. Consistent with the role of Cut homologues in controlling cell death in vertebrates, we find repression of apoptosis regulators by Cux1 in human cancer cells. Finally, we present evidence that suggests that other lineage-restricted specification factors employ a similar mechanism to put the brakes on the oncogenic process.

Apoptosis is a highly conserved cellular function to remove excessive or unstable cells in diverse developmental processes and disease-responses. An important example is the elimination of cells unable to differentiate, which have the potential to generate tumors. Despite the significance of this process, the mechanisms coupling loss of differentiation and apoptosis have remained elusive. Using cell-type specification in Drosophila as a model, we now identify a conserved regulatory logic that underlies cell-type specific removal of uncommitted cells by apoptosis. We find that the transcription factor Cut activates differentiation, while it simultaneously represses cell death via the direct regulation of a pro-apoptotic gene. We show that this regulatory interaction occurs in many diverse cell types and is essential for normal development. Using in vivo Drosophila cancer models, we demonstrate that apoptosis activation in differentiation-compromised cells is an immediate-early cancer prevention mechanism. Importantly, we show that this type of regulatory wiring is also found in vertebrates and that other cell-type specification factors might employ a similar mechanism for tumor suppression. Thus, our findings suggest that the coupling of differentiation and apoptosis by individual transcription factors is a widely used and evolutionarily conserved cancer prevention module, which is hard-wired into the developmental program.

A machine learning approach for identifying novel cell type-specific transcriptional regulators of myogenesis

Transcriptional enhancers integrate the contributions of multiple classes of transcription factors (TFs) to orchestrate the myriad spatio-temporal gene expression programs that occur during development. A molecular understanding of enhancers with similar activities requires the identification of both their unique and their shared sequence features. To address this problem, we combined phylogenetic profiling with a DNA-based enhancer sequence classifier that analyzes the TF binding sites (TFBSs) governing the transcription of a co-expressed gene set. We first assembled a small number of enhancers that are active in Drosophila melanogaster muscle founder cells (FCs) and other mesodermal cell types. Using phylogenetic profiling, we increased the number of enhancers by incorporating orthologous but divergent sequences from other Drosophila species. Functional assays revealed that the diverged enhancer orthologs were active in largely similar patterns as their D. melanogaster counterparts, although there was extensive evolutionary shuffling of known TFBSs. We then built and trained a classifier using this enhancer set and identified additional related enhancers based on the presence or absence of known and putative TFBSs. Predicted FC enhancers were over-represented in proximity to known FC genes; and many of the TFBSs learned by the classifier were found to be critical for enhancer activity, including POU homeodomain, Myb, Ets, Forkhead, and T-box motifs. Empirical testing also revealed that the T-box TF encoded by org-1 is a previously uncharacterized regulator of muscle cell identity. Finally, we found extensive diversity in the composition of TFBSs within known FC enhancers, suggesting that motif combinatorics plays an essential role in the cellular specificity exhibited by such enhancers. In summary, machine learning combined with evolutionary sequence analysis is useful for recognizing novel TFBSs and for facilitating the identification of cognate TFs that coordinate cell type-specific developmental gene expression patterns.

Cofactor binding evokes latent differences in DNA binding specificity between Hox proteins

Members of transcription factor families typically have similar DNA binding specificities yet execute unique functions in vivo. Transcription factors often bind DNA as multiprotein complexes, raising the possibility that complex formation might modify their DNA binding specificities. To test this hypothesis, we developed an experimental and computational platform, SELEX-seq, that can be used to determine the relative affinities to any DNA sequence for any transcription factor complex. Applying this method to all eight Drosophila Hox proteins, we show that they obtain novel recognition properties when they bind DNA with the dimeric cofactor Extradenticle-Homothorax (Exd). Exd-Hox specificities group into three main classes that obey Hox gene collinearity rules and DNA structure predictions suggest that anterior and posterior Hox proteins prefer DNA sequences with distinct minor groove topographies. Together, these data suggest that emergent DNA recognition properties revealed by interactions with cofactors contribute to transcription factor specificities in vivo.

Molecular mechanism underlying the regulatory specificity of a Drosophila homeodomain protein that specifies myoblast identity

A subfamily of Drosophila homeodomain (HD) transcription factors (TFs) controls the identities of individual muscle founder cells (FCs). However, the molecular mechanisms by which these TFs generate unique FC genetic programs remain unknown. To investigate this problem, we first applied genome-wide mRNA expression profiling to identify genes that are activated or repressed by the muscle HD TFs Slouch (Slou) and Muscle segment homeobox (Msh). Next, we used protein-binding microarrays to define the sequences that are bound by Slou, Msh and other HD TFs that have mesodermal expression. These studies revealed that a large class of HDs, including Slou and Msh, predominantly recognize TAAT core sequences but that each HD also binds to unique sites that deviate from this canonical motif. To understand better the regulatory specificity of an individual FC identity HD, we evaluated the functions of atypical binding sites that are preferentially bound by Slou relative to other HDs within muscle enhancers that are either activated or repressed by this TF. These studies showed that Slou regulates the activities of particular myoblast enhancers through Slou-preferred sequences, whereas swapping these sequences for sites that are capable of binding to multiple HD family members does not support the normal regulatory functions of Slou. Moreover, atypical Slou-binding sites are overrepresented in putative enhancers associated with additional Slou-responsive FC genes. Collectively, these studies provide new insights into the roles of individual HD TFs in determining cellular identity, and suggest that the diversity of HD binding preferences can confer regulatory specificity.

Inferring divergence of context-dependent substitution rates in Drosophila genomes with applications to comparative genomics

Nucleotide substitution is a major evolutionary driving force that can incrementally and stochastically give rise to broad divergence patterns among species. The substitution process at each genomic position is frequently modeled independently of the other positions, although complex interactions between nearby bases are known to significantly affect mutation rates. Here, we study the evolution of 12 fly genomes using new algorithms for accurate inference of parameter-rich substitution models. By comparing models between lineages, we reveal the evolutionary histories of substitution rates at different flanking nucleotide contexts. We demonstrate these driving forces of molecular evolution to be constantly changing, suggesting that neutral drift of mutation rates is an important factor in the evolution of genomes and their sequence composition. This observation is used to develop a scalable approach for parameter-rich comparative genomics. By screening short DNA sequences, we demonstrate how homeoboxes and other transcription factor binding motifs are highly conserved based on our parameter-rich models but not according to standard conservation assays. With the increasing availability of genome sequences, rich substitution models become an attractive and practical approach for evolutionary analysis in general and comparative genomics in particular.

Distinct DNA binding and transcriptional repression characteristics related to different ARX mutations

Mutations in the Aristaless-related homeobox gene (ARX) are associated with a wide variety of neurologic disorders including lissencephaly, hydrocephaly, West syndrome, Partington syndrome, and X-linked intellectual disability with or without epilepsy. A genotype-phenotype correlation exists for ARX mutations; however, the molecular basis for this association has not been investigated. To begin understanding the molecular basis for ARX mutations, we tested the DNA binding sequence preference and transcriptional repression activity for Arx, deletion mutants and mutants associated with various neurologic disorders. We found DNA binding preferences of Arx are influenced by the amino acid sequences adjacent to the homeodomain. Mutations in the homeodomain show a loss of DNA binding activity, while the T333N and P353R homeodomain mutants still possess DNA binding activities, although less than the wild type. Transcription repression activity, the primary function of ARX, is reduced in all mutants except the L343Q, which has no DNA binding activity and does not functionally repress Arx targets. These data indicate that mutations in the homeodomain result in not only a loss of DNA binding activity but also loss of transcriptional repression activity. Our results provide novel insights into the pathogenesis of ARX-related disorders and possible directions to pursue potential therapeutic interventions.

Genome-wide occupancy links Hoxa2 to Wnt-β-catenin signaling in mouse embryonic development

The regulation of gene expression is central to developmental programs and largely depends on the binding of sequence-specific transcription factors with cis-regulatory elements in the genome. Hox transcription factors specify the spatial coordinates of the body axis in all animals with bilateral symmetry, but a detailed knowledge of their molecular function in instructing cell fates is lacking. Here, we used chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) to identify Hoxa2 genomic locations in a time and space when it is actively instructing embryonic development in mouse. Our data reveals that Hoxa2 has large genome coverage and potentially regulates thousands of genes. Sequence analysis of Hoxa2-bound regions identifies high occurrence of two main classes of motifs, corresponding to Hox and Pbx-Hox recognition sequences. Examination of the binding targets of Hoxa2 faithfully captures the processes regulated by Hoxa2 during embryonic development; in addition, it uncovers a large cluster of potential targets involved in the Wnt-signaling pathway. In vivo examination of canonical Wnt-β-catenin signaling reveals activity specifically in Hoxa2 domain of expression, and this is undetectable in Hoxa2 mutant embryos. The comprehensive mapping of Hoxa2-binding sites provides a framework to study Hox regulatory networks in vertebrate developmental processes.

Analysis of specific protein-DNA interactions by bacterial one-hybrid assay

The DNA-binding specificity of transcription factors allows the prediction of regulatory targets in a genome. However, very few factor specificities have been characterized and still too little is known about how these proteins interact with their targets to make predictions a priori. To provide a greater understanding of how proteins and DNA interact, we have developed a bacterial one-hybrid system that allows the sensitive, high-throughput, and cost-effective assay of the interaction at the protein-DNA interface. This system makes survival of the bacteria dependent on activation of the reporter gene and therefore dependent on the protein-DNA interaction that recruits the polymerase. We have used this system to characterize DNA-binding specificities for representative members of the most common DNA-binding domain (DBD) families. We have also been able to engineer DBDs with novel specificity to be used as artificial transcription factors and zinc finger nucleases. The B1H assay provides a simple and inexpensive method to investigate protein-DNA interactions that is accessible to essentially any laboratory.

Comparative genomics-based identification and analysis of cis-regulatory elements

Identification of cis-regulatory elements, such as enhancers and promoters, is very important not only for analysis of gene regulatory networks but also as a tool for targeted gene expression experiments. In this chapter, we introduce an easy but reliable approach to predict enhancers of a gene of interest by comparing mammalian and Xenopus genome sequences, and to examine their activity using a co-transgenesis technique in Xenopus embryos. Since the bioinformatics analysis utilizes publically available web tools, bench biologists can easily perform it without any need for special computing capability. The co-transgenesis assay, which directly uses polymerase chain reaction products, quickly screens for the activity of the candidate elements in a cloning-free manner.

Roles for intrinsic disorder and fuzziness in generating context-specific function in Ultrabithorax, a Hox transcription factor

Surprisingly few transcription factors drive animal development relative to the number and diversity of final tissues and body structures. Therefore, most transcription factors must function in more than one tissue. In a famous example, members of the Hox transcription factor family are expressed in contiguous stripes along the anterior/posterior axis during animal development. Individual Hox transcription factors specify all tissues within their expression domain and thus must respond to cellular cues to instigate the correct tissue-specific gene regulatory cascade. We describe how, in the Drosophila Hox protein Ultrabithorax, intrinsically disordered regions implement, regulate and co-ordinate multiple functions, potentially enabling context-specific gene regulation. The large N-terminal disordered domain encodes most of the transcription activation domain and directly impacts DNA binding affinity by the Ubx homeodomain. Similarly, the C-terminal disordered domain alters DNA binding affinity and specificity, interaction with a Hox binding protein and strongly influences both transcription activation and repression. Phosphorylation of the N-terminal disordered domain and alternative splicing of the C-terminal disordered domain could allow the cell to both regulate and co-ordinate DNA binding, protein interactions and transcription regulation. For regulatory mechanisms relying on disorder to continue to be available when Ubx is bound to other proteins or DNA, fuzziness would need to be preserved in these macromolecular complexes. The intrinsically disordered domains in Hox proteins are predicted to be on the very dynamic end of the disorder spectrum, potentially allowing disorder to persist when Ubx is bound to proteins or DNA to regulate the function of these "fuzzy" complexes. Because both intrinsically disordered regions within Ubx have multiple roles, each region may implement several different regulatory mechanisms identified in fuzzy complexes. These intrinsic disorder-based regulatory mechanisms are likely to be critical for allowing Ubx to sense tissue identity and respond by implementing a context-specific gene regulatory cascade.

Genome-level identification of targets of Hox protein Ultrabithorax in Drosophila: novel mechanisms for target selection

Hox proteins are transcription factors and key regulators of segmental identity along the anterior posterior axis across all bilaterian animals. Despite decades of research, the mechanisms by which Hox proteins select and regulate their targets remain elusive. We have carried out whole-genome ChIP-chip experiments to identify direct targets of Hox protein Ultrabithorax (Ubx) during haltere development in Drosophila. Direct targets identified include upstream regulators or cofactors of Ubx. Homothorax, a cofactor of Ubx during embryonic development, is one such target and is required for normal specification of haltere. Although Ubx bound sequences are conserved amongst various insect genomes, no consensus Ubx-specific motif was detected. Surprisingly, binding motifs for certain transcription factors that function either upstream or downstream to Ubx are enriched in these sequences suggesting complex regulatory loops governing Ubx function. Our data supports the hypothesis that specificity during Hox target selection is achieved by associating with other transcription factors.

Variable motif utilization in homeotic selector (Hox)-cofactor complex formation controls specificity

Homeotic selector (Hox) proteins often bind DNA cooperatively with cofactors such as Extradenticle (Exd) and Homothorax (Hth) to achieve functional specificity in vivo. Previous studies identified the Hox YPWM motif as an important Exd interaction motif. Using a comparative approach, we characterize the contribution of this and additional conserved sequence motifs to the regulation of specific target genes for three Drosophila Hox proteins. We find that Sex combs reduced (Scr) uses a simple interaction mechanism, where a single tryptophan-containing motif is necessary for Exd-dependent DNA-binding and in vivo functions. Abdominal-A (AbdA) is more complex, using multiple conserved motifs in a context-dependent manner. Lastly, Ultrabithorax (Ubx) is the most flexible, in that it uses multiple conserved motifs that function in parallel to regulate target genes in vivo. We propose that using different binding mechanisms with the same cofactor may be one strategy to achieve functional specificity in vivo.

Functional dissection of the Hox protein Abdominal-B in Drosophila cell culture

Hox transcription factors regulate the morphogenesis along the anterior-posterior (A/P) body axis through the interaction with small cis-regulatory modules (CRMs) of their target gene, however so far very few Hox CRMs are known and have been analyzed in detail. In this study we have identified a new Hox CRM, ct340, which guides the expression of the cell type specification gene cut (ct) in the posterior spiracle under the direct control of the Hox protein Abdominal-B (Abd-B). Using the ct340 enhancer activity as readout, an efficient cloning system to generate VP16 activation domain fusion protein was developed to unambiguously test protein-DNA interaction in Drosophila cell culture. By functionally dissecting the Abd-B protein, new features of Abd-B dependent target gene regulation were detected. Due to its easy adaptability, this system can be generally used to map functional domains within sequence-specific transcriptional factors in Drosophila cell culture, and thus provide preliminary knowledge of the protein functional domain structure for further in vivo analysis.

Reengineering Cro protein functional specificity with an evolutionary code

Cro proteins from different lambdoid bacteriophages are extremely variable in their target consensus DNA sequences and constitute an excellent model for evolution of transcription factor specificity. We experimentally tested a bioinformatically derived evolutionary code relating switches between pairs of amino acids at three recognition helix sites in Cro proteins to switches between pairs of nucleotide bases in the cognate consensus DNA half-sites. We generated all eight possible code variants of bacteriophage λ Cro and used electrophoretic mobility shift assays to compare binding of each variant to its own putative cognate site and to the wild-type cognate site; we also tested the wild-type protein against all eight DNA sites. Each code variant showed stronger binding to its putative cognate site than to the wild-type site, except some variants containing proline at position 27; each also bound its cognate site better than wild-type Cro bound the same site. Most code variants, however, displayed poorer affinity and specificity than wild-type λ Cro. Fluorescence anisotropy assays on λ Cro and the triple code variant (PSQ) against the two cognate sites confirmed the switch in specificity and showed larger apparent effects on binding affinity and specificity. Bacterial one-hybrid assays of λ Cro and PSQ against libraries of sequences with a single randomized half-site showed the expected switches in specificity at two of three coded positions and no clear switches in specificity at noncoded positions. With a few caveats, these results confirm that the proposed Cro evolutionary code can be used to reengineer Cro specificity.

How to make stripes: deciphering the transition from non-periodic to periodic patterns in Drosophila segmentation

The generation of metameric body plans is a key process in development. In Drosophila segmentation, periodicity is established rapidly through the complex transcriptional regulation of the pair-rule genes. The 'primary' pair-rule genes generate their 7-stripe expression through stripe-specific cis-regulatory elements controlled by the preceding non-periodic maternal and gap gene patterns, whereas 'secondary' pair-rule genes are thought to rely on 7-stripe elements that read off the already periodic primary pair-rule patterns. Using a combination of computational and experimental approaches, we have conducted a comprehensive systems-level examination of the regulatory architecture underlying pair-rule stripe formation. We find that runt (run), fushi tarazu (ftz) and odd skipped (odd) establish most of their pattern through stripe-specific elements, arguing for a reclassification of ftz and odd as primary pair-rule genes. In the case of run, we observe long-range cis-regulation across multiple intervening genes. The 7-stripe elements of run, ftz and odd are active concurrently with the stripe-specific elements, indicating that maternal/gap-mediated control and pair-rule gene cross-regulation are closely integrated. Stripe-specific elements fall into three distinct classes based on their principal repressive gap factor input; stripe positions along the gap gradients correlate with the strength of predicted input. The prevalence of cis-elements that generate two stripes and their genomic organization suggest that single-stripe elements arose by splitting and subfunctionalization of ancestral dual-stripe elements. Overall, our study provides a greatly improved understanding of how periodic patterns are established in the Drosophila embryo.

Correlated evolution of transcription factors and their binding sites

MOTIVATION

The interaction between transcription factor (TF) and transcription factor binding site (TFBS) is essential for gene regulation. Mutation in either the TF or the TFBS may weaken their interaction and thus result in abnormalities. To maintain such vital interaction, a mutation in one of the interacting partners might be compensated by a corresponding mutation in its binding partner during the course of evolution. Confirming this co-evolutionary relationship will guide us in designing protein sequences to target a specific DNA sequence or in predicting TFBS for poorly studied proteins, or even correcting and rescuing disease mutations in clinical applications.

RESULTS

Based on six, publicly available, experimentally validated TF-TFBS binding datasets for the basic Helix-Loop-Helix (bHLH) family, Homeo family, High-Mobility Group (HMG) family and Transient Receptor Potential channels (TRP) family, we showed that the evolutions of the TFs and their TFBSs are significantly correlated across eukaryotes. We further developed a mutual information-based method to identify co-evolved protein residues and DNA bases. This research sheds light on the dynamic relationship between TF and TFBS during their evolution. The same principle and strategy can be applied to co-evolutionary studies on protein-DNA interactions in other protein families.

AVAILABILITY

All the datasets, scripts and other related files have been made freely available at: http://jjwanglab.org/co-evo.

CONTACT

junwen@uw.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Iroquois homeobox gene 3 establishes fast conduction in the cardiac His-Purkinje network

Rapid electrical conduction in the His-Purkinje system tightly controls spatiotemporal activation of the ventricles. Although recent work has shed much light on the regulation of early specification and morphogenesis of the His-Purkinje system, less is known about how transcriptional regulation establishes impulse conduction properties of the constituent cells. Here we show that Iroquois homeobox gene 3 (Irx3) is critical for efficient conduction in this specialized tissue by antithetically regulating two gap junction-forming connexins (Cxs). Loss of Irx3 resulted in disruption of the rapid coordinated spread of ventricular excitation, reduced levels of Cx40, and ectopic Cx43 expression in the proximal bundle branches. Irx3 directly represses Cx43 transcription and indirectly activates Cx40 transcription. Our results reveal a critical role for Irx3 in the precise regulation of intercellular gap junction coupling and impulse propagation in the heart.

Spatial regulation of lag-2 transcription during vulval precursor cell fate patterning in Caenorhabditis elegans

lag-2 encodes a ligand for LIN-12/Notch and is a component of the lateral signal that activates LIN-12/Notch during Caenorhabditis elegans vulval precursor cell (VPC) fate patterning. lag-2 is specifically transcribed in one VPC, named P6.p, in response to activation of EGFR/Ras/MAPK by the inductive signal that initiates vulval development. Here, we show that a critical molecular event linking inductive and lateral signaling is the relief of VPC-wide lag-2 repression in P6.p. We find that the lag-2 promoter contains an element, VPCrep, which mediates repression in all VPCs when the inductive signal is absent, and another promoter element, VPCact, which is required for activation when repression is relieved by the inductive signal. We show that repression through VPCrep is mediated by the Elk1 ortholog LIN-1, and that the level and subcellular accumulation of a functional LIN-1::GFP protein is similar in all six VPCs before and after vulval induction, suggesting that relief of LIN-1-mediated repression in P6.p is likely due to the known MAPK-dependent phosphorylation of LIN-1. We also provide evidence that the factor(s) acting through VPCact is present in all VPCs but is not modulated by the inductive signal, and that transcription of lag-2 requires the Hth/Meis ortholog UNC-62 and the Mediator complex component SUR-2. Relief of repression of lag-2 in P6.p offers a plausible mechanistic basis for spatial restriction of lag-2 in generating the precise spatial pattern of VPC fates.

Mathematical modeling of gene expression: a guide for the perplexed biologist

The detailed analysis of transcriptional networks holds a key for understanding central biological processes, and interest in this field has exploded due to new large-scale data acquisition techniques. Mathematical modeling can provide essential insights, but the diversity of modeling approaches can be a daunting prospect to investigators new to this area. For those interested in beginning a transcriptional mathematical modeling project, we provide here an overview of major types of models and their applications to transcriptional networks. In this discussion of recent literature on thermodynamic, Boolean, and differential equation models, we focus on considerations critical for choosing and validating a modeling approach that will be useful for quantitative understanding of biological systems.

Direct measurement of DNA affinity landscapes on a high-throughput sequencing instrument

Several methods for characterizing DNA-protein interactions are available, but none have demonstrated both high throughput and quantitative measurement of affinity. Here we describe 'high-throughput sequencing'-'fluorescent ligand interaction profiling' (HiTS-FLIP), a technique for measuring quantitative protein-DNA binding affinity at unprecedented depth. In this approach, the optics built into a high-throughput sequencer are used to visualize in vitro binding of a protein to sequenced DNA in a flow cell. Application of HiTS-FLIP to the protein Gcn4 (Gcn4p), the master regulator of the yeast amino acid starvation response, yielded ~440 million binding measurements, enabling determination of dissociation constants for all 12-mer sequences having submicromolar affinity. These data revealed a complex interdependency between motif positions, allowed improved discrimination of in vivo Gcn4p binding sites and regulatory targets relative to previous methods and showed that sets of genes with different promoter affinities to Gcn4p have distinct functions and expression kinetics. Broad application of this approach should increase understanding of the interactions that drive transcription.

Establishment of medial fates along the proximodistal axis of the Drosophila leg through direct activation of dachshund by Distalless

The proximodistal (PD) axis of the Drosophila leg is thought to be established by the combined gradients of two secreted morphogens, Wingless (Wg) and Decapentaplegic (Dpp). According to this model, high [Wg+Dpp] activates Distalless (Dll) and represses dachshund (dac) in the distal cells of the leg disc, while intermediate [Wg+Dpp] activates dac in medial tissue. To test this model we identified and characterized a dac cis-regulatory element (dac RE) that recapitulates dac's medial expression domain during leg development. Counter to the gradient model, we find that Wg and Dpp do not act in a graded manner to activate RE. Instead, dac RE is activated directly by Dll and repressed distally by a combination of factors, including the homeodomain protein Bar. Thus, medial leg fates are established via a regulatory cascade in which Wg+Dpp activate Dll and then Dll directly activates dac, with Wg+Dpp as less critical, permissive inputs.

Systematic characterization of protein-DNA interactions

Sequence-specific protein-DNA interactions (PDIs) are critical for regulating many cellular processes, including transcription, DNA replication, repair, and rearrangement. We review recent experimental advances in high-throughput technologies designed to characterize PDIs and discuss recent studies that use these tools, including ChIP-chip/seq, SELEX-based approaches, yeast one-hybrid, bacterial one-hybrid, protein binding microarray, and protein microarray. The results of these studies have challenged some long-standing concepts of PDI and provide valuable insights into the complex transcriptional regulatory networks.

Characterization of protein-DNA interactions using protein microarrays

INTRODUCTION

Protein-DNA interactions (PDIs) are critical for many cellular processes. We present here a protocol for the identification of PDIs in vitro using protein microarray technology. The procedure involves double-stranding synthesized DNA oligonucleotides with a fluorescent-labeled primer, binding the labeled double-stranded DNA directly to the protein microarray, and analyzing binding of the resulting PDIs. This approach provides simultaneous identification of PDIs for thousands of proteins, and multiple carefully designed DNA probes can be tested in parallel, which enables a rapid mapping of PDIs on a proteome-wide scale.

A modified bacterial one-hybrid system yields improved quantitative models of transcription factor specificity

We examine the use of high-throughput sequencing on binding sites recovered using a bacterial one-hybrid (B1H) system and find that improved models of transcription factor (TF) binding specificity can be obtained compared to standard methods of sequencing a small subset of the selected clones. We can obtain even more accurate binding models using a modified version of B1H selection method with constrained variation (CV-B1H). However, achieving these improved models using CV-B1H data required the development of a new method of analysis--GRaMS (Growth Rate Modeling of Specificity)--that estimates bacterial growth rates as a function of the quality of the recognition sequence. We benchmark these different methods of motif discovery using Zif268, a well-characterized C(2)H(2) zinc-finger TF on both a 28 bp randomized library for the standard B1H method and on 6 bp randomized library for the CV-B1H method for which 45 different experimental conditions were tested: five time points and three different IPTG and 3-AT concentrations. We find that GRaMS analysis is robust to the different experimental parameters whereas other analysis methods give widely varying results depending on the conditions of the experiment. Finally, we demonstrate that the CV-B1H assay can be performed in liquid media, which produces recognition models that are similar in quality to sequences recovered from selection on solid media.

Differential selection within the Drosophila retinal determination network and evidence for functional divergence between paralog pairs

The retinal determination (RD) network in Drosophila comprises 14 known nuclear proteins that include DNA-binding proteins, transcriptional coactivators, kinases, and phosphatases. The composition of the network varies considerably throughout the animal kingdom, with the network in several basal insects having fewer members and with vertebrates having potentially significantly higher numbers of RD genes. One important contributing factor for the variation in gene number within the network is gene duplication. For example, 10 members of the RD network in Drosophila are derived from duplication events. Here we present an analysis of the coding regions of the five pairs of duplicate genes from within the RD network of several different Drosophila species. We demonstrate that there is differential selection across the coding regions of all RD genes. Additionally, some of the most significant differences in ratios of non-silent-to-silent site substitutions (d(N)/d(S)) between paralog pairs are found within regions that have no ascribed function. Previous structure/function analyses of several duplicate genes have identified areas within one gene that contain novel activities when compared with its paralog. The evolutionary analysis presented here identifies these same areas in the paralogs as being under high levels of relaxed selection. We suggest that sequence divergence between paralogs and selection signatures can be used as a reasonable predictor of functional changes in rapidly evolving motifs.

Genome-wide analysis of the binding of the Hox protein Ultrabithorax and the Hox cofactor Homothorax in Drosophila

Hox genes encode a family of transcription factors that are key developmental regulators with a highly conserved role in specifying segmental diversity along the metazoan body axis. Although they have been shown to regulate a wide variety of downstream processes, direct transcriptional targets have been difficult to identify and this has been a major obstacle to our understanding of Hox gene function. We report the identification of genome-wide binding sites for the Hox protein Ultrabithorax (Ubx) using a YFP-tagged Drosophila protein-trap line together with chromatin immunoprecipitation and microarray analysis. We identify 1,147 genes bound by Ubx at high confidence in chromatin from the haltere imaginal disc, a prominent site of Ubx function where it specifies haltere versus wing development. The functional relevance of these genes is supported by their overlap with genes differentially expressed between wing and haltere imaginal discs. The Ubx-bound gene set is highly enriched in genes involved in developmental processes and contains both high-level regulators as well as genes involved in more basic cellular functions. Several signalling pathways are highly enriched in the Ubx target gene set and our analysis supports the view that Hox genes regulate many levels of developmental pathways and have targets distributed throughout the gene network. We also performed genome-wide analysis of the binding sites for the Hox cofactor Homothorax (Hth), revealing a striking similarity with the Ubx binding profile. We suggest that these binding profiles may be strongly influenced by chromatin accessibility and provide evidence of a link between Ubx/Hth binding and chromatin state at genes regulated by Polycomb silencing. Overall, we define a set of direct Ubx targets in the haltere imaginal disc and suggest that chromatin accessibility has important implications for Hox target selection and for transcription factor binding in general.

Genome-wide tissue-specific occupancy of the Hox protein Ultrabithorax and Hox cofactor Homothorax in Drosophila

The Hox genes are responsible for generating morphological diversity along the anterior-posterior axis during animal development. The Drosophila Hox gene Ultrabithorax (Ubx), for example, is required for specifying the identity of the third thoracic (T3) segment of the adult, which includes the dorsal haltere, an appendage required for flight, and the ventral T3 leg. Ubx mutants show homeotic transformations of the T3 leg towards the identity of the T2 leg and the haltere towards the wing. All Hox genes, including Ubx, encode homeodomain containing transcription factors, raising the question of what target genes Ubx regulates to generate these adult structures. To address this question, we carried out whole genome ChIP-chip studies to identify all of the Ubx bound regions in the haltere and T3 leg imaginal discs, which are the precursors to these adult structures. In addition, we used ChIP-chip to identify the sites bound by the Hox cofactor, Homothorax (Hth). In contrast to previous ChIP-chip studies carried out in Drosophila embryos, these binding studies reveal that there is a remarkable amount of tissue- and transcription factor-specific binding. Analyses of the putative target genes bound and regulated by these factors suggest that Ubx regulates many downstream transcription factors and developmental pathways in the haltere and T3 leg. Finally, we discovered additional DNA sequence motifs that in some cases are specific for individual data sets, arguing that Ubx and/or Hth work together with many regionally expressed transcription factors to execute their functions. Together, these data provide the first whole-genome analysis of the binding sites and target genes regulated by Ubx to specify the morphologies of the adult T3 segment of the fly.

Using the recognition code to swap homeodomain target specificity in cell culture

The homeodomain (HD) is a 60 amino acid-long DNA-binding domain. A large fraction of HDs binds with high affinity sequences containing the 5'-TAAT-3' core motif. However, NK-2 class HDs recognizes sequences containing the 5'-CAAG-3' core motif. By using a cell transfection approach, here we show that modification of residues located in the N-terminal arm (at positions 6, 7 and 8) and in the recognition helix (at position 54) is enough to swap the "in vivo" binding specificity of TTF-1 HD (which is a member of the NK-2 class HD) from 5'-CAAG-3' to 5'-TAAT-3'-containing targets. The role of residue at position 54 is also supported by data obtained with the HD of the Drosophila engrailed protein. These data support the notion that DNA-binding specificity "in vivo" is dictated by few critical residues.

Novel computational analysis of protein binding array data identifies direct targets of Nkx2.2 in the pancreas

Background

The creation of a complete genome-wide map of transcription factor binding sites is essential for understanding gene regulatory networks in vivo. However, current prediction methods generally rely on statistical models that imperfectly model transcription factor binding. Generation of new prediction methods that are based on protein binding data, but do not rely on these models may improve prediction sensitivity and specificity.

Results

We propose a method for predicting transcription factor binding sites in the genome by directly mapping data generated from protein binding microarrays (PBM) to the genome and calculating a moving average of several overlapping octamers. Using this unique algorithm, we predicted binding sites for the essential pancreatic islet transcription factor Nkx2.2 in the mouse genome and confirmed >90% of the tested sites by EMSA and ChIP. Scores generated from this method more accurately predicted relative binding affinity than PWM based methods. We have also identified an alternative core sequence recognized by the Nkx2.2 homeodomain. Furthermore, we have shown that this method correctly identified binding sites in the promoters of two critical pancreatic islet β-cell genes, NeuroD1 and insulin2, that were not predicted by traditional methods. Finally, we show evidence that the algorithm can also be applied to predict binding sites for the nuclear receptor Hnf4α.

Conclusions

PBM-mapping is an accurate method for predicting Nkx2.2 binding sites and may be widely applicable for the creation of genome-wide maps of transcription factor binding sites.

Using a structural and logics systems approach to infer bHLH-DNA binding specificity determinants

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.

Retinal determination the beginning of eye development

The road to producing an eye begins with the decision to commit a population of cells to adopting an eye tissue fate, the process of retinal determination. Over the past decade and a half, a network of transcription factors has been found to mediate this process in all seeing animals. This retinal determination network is known to regulate not only tissue fate but also cell proliferation, pattern formation, compartment boundary establishment, and even retinal cell specification. The compound eye of the fruit fly, Drosophila melanogaster, has proven to be an excellent experimental system to study the mechanisms by which this network regulates organogenesis and tissue patterning. In fact the founding members of most of the gene families that make up this network were first isolated in Drosophila based on loss-of-function phenotypes that affect the eye. This chapter will highlight the history of discovery of the retinal determination network and will draw attention to the molecular and biochemical mechanisms that underlie our understanding of how the fate of the retina is determined.

Transient activation of meox1 is an early component of the gene regulatory network downstream of hoxa2

Hox genes encode transcription factors that regulate morphogenesis in all animals with bilateral symmetry. Although Hox genes have been extensively studied, their molecular function is not clear in vertebrates, and only a limited number of genes regulated by Hox transcription factors have been identified. Hoxa2 is required for correct development of the second branchial arch, its major domain of expression. We now show that Meox1 is genetically downstream from Hoxa2 and is a direct target. Meox1 expression is downregulated in the second arch of Hoxa2 mouse mutant embryos. In chromatin immunoprecipitation (ChIP), Hoxa2 binds to the Meox1 proximal promoter. Two highly conserved binding sites contained in this sequence are required for Hoxa2-dependent activation of the Meox1 promoter. Remarkably, in the absence of Meox1 and its close homolog Meox2, the second branchial arch develops abnormally and two of the three skeletal elements patterned by Hoxa2 are malformed. Finally, we show that Meox1 can specifically bind the DNA sequences recognized by Hoxa2 on its functional target genes. These results provide new insight into the Hoxa2 regulatory network that controls branchial arch identity.

A gene regulatory network controlling hhex transcription in the anterior endoderm of the organizer

The homeobox gene hhex is one of the earliest markers of the anterior endoderm, which gives rise to foregut organs such as the liver, ventral pancreas, thyroid, and lungs. The regulatory networks controlling hhex transcription are poorly understood. In an extensive cis-regulatory analysis of the Xenopus hhex promoter, we determined how the Nodal, Wnt, and BMP pathways and their downstream transcription factors regulate hhex expression in the gastrula organizer. We show that Nodal signaling, present throughout the endoderm, directly activates hhex transcription via FoxH1/Smad2 binding sites in the proximal -0.44 Kb promoter. This positive action of Nodal is suppressed in the ventral-posterior endoderm by Vent 1 and Vent2, homeodomain repressors that are induced by BMP signaling. Maternal Wnt/β-catenin on the dorsal side of the embryo cooperates with Nodal and indirectly activates hhex expression via the homeodomain activators Siamois and Twin. Siamois/Twin stimulate hhex transcription through two mechanisms: (1) they induce the expression of Otx2 and Lim1 and together Siamois, Twin, Otx2, and Lim1 appear to promote hhex transcription through homeobox sites in a Wnt-responsive element located between -0.65 to -0.55 Kb of the hhex promoter. (2) Siamois/Twin also induce the expression of the BMP-antagonists Chordin and Noggin, which are required to exclude Vents from the organizer allowing hhex transcription. This study reveals a complex network regulating anterior endoderm transcription in the early embryo.

Quantitative approaches to monitor protein-nucleic acid interactions using fluorescent probes

Sequence-specific recognition of nucleic acids by proteins is required for nearly every aspect of gene expression. Quantitative binding experiments are a useful tool to measure the ability of a protein to distinguish between multiple sequences. Here, we describe the use of fluorophore-labeled oligonucleotide probes to quantitatively monitor protein/nucleic acid interactions. We review two complementary experimental methods, fluorescence polarization and fluorescence electrophoretic mobility shift assays, that enable the quantitative measurement of binding affinity. We also present two strategies for post-synthetic end-labeling of DNA or RNA oligonucleotides with fluorescent dyes. The approaches discussed here are efficient and sensitive, providing a safe and accessible alternative to the more commonly used radio-isotopic methods.

The gap gene network

Gap genes are involved in segment determination during the early development of the fruit fly Drosophila melanogaster as well as in other insects. This review attempts to synthesize the current knowledge of the gap gene network through a comprehensive survey of the experimental literature. I focus on genetic and molecular evidence, which provides us with an almost-complete picture of the regulatory interactions responsible for trunk gap gene expression. I discuss the regulatory mechanisms involved, and highlight the remaining ambiguities and gaps in the evidence. This is followed by a brief discussion of molecular regulatory mechanisms for transcriptional regulation, as well as precision and size-regulation provided by the system. Finally, I discuss evidence on the evolution of gap gene expression from species other than Drosophila. My survey concludes that studies of the gap gene system continue to reveal interesting and important new insights into the role of gene regulatory networks in development and evolution.

Genomic approaches to understanding Hox gene function

For many years, biologists have sought to understand how the homeodomain-containing transcriptional regulators encoded by Hox genes are able to control the development of animal morphology. Almost a century of genetics and several decades of molecular biology have defined the conserved organization of homeotic gene clusters in animals and the basic molecular properties of Hox transcription factors. In contrast to these successes, we remain relatively ignorant of how Hox proteins find their target genes in the genome or what sets of genes a Hox protein regulates to direct morphogenesis. The recent deployment of genomic methods, such as whole transcriptome mRNA expression profiling and genome-wide analysis of protein-DNA interactions, begins to shed light on these issues. Results from such studies, principally in the fruit fly, indicate that Hox proteins control the expression of hundreds, if not thousands, of genes throughout the gene regulatory network and that, in many cases, the effects on the expression of individual genes may be quite subtle. Hox proteins regulate both high-level effectors, including other transcription factors and signaling molecules, as well as the cytodifferentiation genes or Realizators at the bottom of regulatory hierarchies. Insights emerging from mapping Hox binding sites in the genome begin to suggest that Hox binding may be strongly influenced by chromatin accessibility rather than binding site affinity. If this is the case, it indicates we need to refocus our efforts at understanding Hox function toward the dynamics of gene regulatory networks and chromatin epigenetics.

Methods for Analysis of Transcription Factor DNA-Binding Specificity In Vitro

Transcription of genes during development and in response to environmental stimuli is determined by genomic DNA sequence. The DNA sequences regulating transcription are read by sequence-specific transcription factors (TFs) that recognize relatively short sequences, generally between four and twenty base pairs in length. Transcriptional regulation generally requires binding of multiple TFs in close proximity to each other. Mechanistic understanding of transcription in an organism thus requires detailed knowledge of binding affinities of all its TFs to all possible DNA sequences, and the co-operative interactions between the TFs. However, very little is known about such co-operative binding interactions, and even the simple TF-DNA binding information exists only for a very small proportion of all TFs - for example, mammals have approximately 1,300-2,000 TFs [1, 2], yet the largest public databases for TF binding specificity, Jaspar and Uniprobe [3, 4] currently list only approximately 500 moderate to high resolution profiles for human or mouse. This lack of knowledge is in part due to the fact that analysis of TF DNA binding has been laborious and expensive. In this chapter, we review methods that can be used to determine binding specificity of TFs to DNA, mainly focusing on recently developed assays that allow high-resolution analysis of TF binding specificity in relatively high throughput.

Refolding the engrailed homeodomain: structural basis for the accumulation of a folding intermediate

The ultrafast folding pathway of the engrailed homeodomain has been exceptionally well characterized by experiment and simulation. Helices II and III of the three-helix bundle protein form the native helix-turn-helix motif as an on-pathway intermediate within a few microseconds. The slow step is then the proper docking of the helices in approximately 15 mus. However, there is still the unexplained puzzle of why helix docking is relatively slow, which is part of the more general question as to why rearrangements of intermediates occur slowly. To address this problem, we performed 46 all-atom molecular dynamics refolding simulations in explicit water, for a total of 15 micros of simulation time. The simulations started from an intermediate state structure that was generated in an unfolding simulation at 498 K and was then quenched to folding-permissive temperatures. The protein refolded successfully in only one of the 46 simulations, and in that case the refolding pathway mirrored the unfolding pathway at high temperature. In the 45 simulations in which the protein did not fully fold, nonnative salt bridges trapped the protein, which explains why the protein folds relatively slowly from the intermediate state.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

Control of DNA replication: a new facet of Hox proteins?

Hox proteins are well-known as developmental transcription factors controlling cell and tissue identity, but recent findings suggest that they are also part of the cell replication machinery. Hox-mediated control of transcription and replication may ensure coordinated control of cell growth and differentiation, two processes that need to be tightly and precisely coordinated to allow proper organ formation and patterning. In this review we summarize the available data linking Hox proteins to the replication machinery and discuss the developmental and pathological implications of this new facet of Hox protein function.

FlyFactorSurvey: a database of Drosophila transcription factor binding specificities determined using the bacterial one-hybrid system

FlyFactorSurvey (http://pgfe.umassmed.edu/TFDBS/) is a database of DNA binding specificities for Drosophila transcription factors (TFs) primarily determined using the bacterial one-hybrid system. The database provides community access to over 400 recognition motifs and position weight matrices for over 200 TFs, including many unpublished motifs. Search tools and flat file downloads are provided to retrieve binding site information (as sequences, matrices and sequence logos) for individual TFs, groups of TFs or for all TFs with characterized binding specificities. Linked analysis tools allow users to identify motifs within our database that share similarity to a query matrix or to view the distribution of occurrences of an individual motif throughout the Drosophila genome. Together, this database and its associated tools provide computational and experimental biologists with resources to predict interactions between Drosophila TFs and target cis-regulatory sequences.

Homeoprotein Six1 increases TGF-beta type I receptor and converts TGF-beta signaling from suppressive to supportive for tumor growth

The Six1 homeodomain protein is a developmental transcription factor that has been implicated in tumor onset and progression. Our recent work shows that Six1 overexpression in human breast cancer cell lines is sufficient to induce epithelial-to-mesenchymal transition (EMT) and metastasis. Importantly, Six1-induced EMT and metastasis are dependent on TGF-β signaling. The TGF-β pathway plays a dual role in cancer, acting as a tumor suppressor in early lesions but enhancing metastatic spread in more advanced tumors. Our previous work indicated that Six1 may be a critical mediator of the switch in TGF-β signaling from tumor suppressive to tumor promotional. However, the mechanism by which Six1 impinges on the TGF-β pathway was, until now, unclear. In this work, we identify the TGF-β type I receptor (TβRI) as a target of Six1 and a critical effector of Six1-induced TGF-β signaling and EMT. We show that Six1-induced upregulation of TβRI is both necessary and sufficient to activate TGF-β signaling and induce properties of EMT. Interestingly, increased TβRI expression is not sufficient to induce experimental metastasis, providing in vivo evidence that Six1 overexpression is required to switch TGF-β signaling to the prometastatic phenotype and showing that induction of EMT is not sufficient to induce experimental metastasis. Together, these results show a novel mechanism for the activation of TGF-β signaling, identify TβRI as a new target of Six1, and implicate Six1 as a determinant of TGF-β function in breast cancer.

Identification and genome-wide prediction of DNA binding specificities for the ApiAP2 family of regulators from the malaria parasite

The molecular mechanisms underlying transcriptional regulation in apicomplexan parasites remain poorly understood. Recently, the Apicomplexan AP2 (ApiAP2) family of DNA binding proteins was identified as a major class of transcriptional regulators that are found across all Apicomplexa. To gain insight into the regulatory role of these proteins in the malaria parasite, we have comprehensively surveyed the DNA-binding specificities of all 27 members of the ApiAP2 protein family from Plasmodium falciparum revealing unique binding preferences for the majority of these DNA binding proteins. In addition to high affinity primary motif interactions, we also observe interactions with secondary motifs. The ability of a number of ApiAP2 proteins to bind multiple, distinct motifs significantly increases the potential complexity of the transcriptional regulatory networks governed by the ApiAP2 family. Using these newly identified sequence motifs, we infer the trans-factors associated with previously reported plasmodial cis-elements and provide evidence that ApiAP2 proteins modulate key regulatory decisions at all stages of parasite development. Our results offer a detailed view of ApiAP2 DNA binding specificity and take the first step toward inferring comprehensive gene regulatory networks for P. falciparum.

Determining the specificity of protein-DNA interactions

Proteins, such as many transcription factors, that bind to specific DNA sequences are essential for the proper regulation of gene expression. Identifying the specific sequences that each factor binds can help to elucidate regulatory networks within cells and how genetic variation can cause disruption of normal gene expression, which is often associated with disease. Traditional methods for determining the specificity of DNA-binding proteins are slow and laborious, but several new high-throughput methods can provide comprehensive binding information much more rapidly. Combined with in vivo determinations of transcription factor binding locations, this information provides more detailed views of the regulatory circuitry of cells and the effects of variation on gene expression.

Experimental strategies for studying transcription factor-DNA binding specificities

Specific binding of transcription factors (TFs) determines in a large part the connectivity of gene regulatory networks as well as the quantitative level of gene expression. A multiplicity of both experimental and computational methods is currently used to discover and characterize the underlying TF-DNA interactions. Experimental methods can be further subdivided into in vitro- and in vivo-based approaches, each accenting different aspects of TF-binding events. In this review we summarize the flexibility and performance of a selection of both types of experimental methods. In conclusion, we argue that a serial combination of methods with different throughput and data type constitutes an optimal experimental strategy.

Foxa2 regulates the expression of Nato3 in the floor plate by a novel evolutionarily conserved promoter

The development of the neural tube into a complex central nervous system involves morphological, cellular and molecular changes, all of which are tightly regulated. The floor plate (FP) is a critical organizing center located at the ventral-most midline of the neural tube. FP cells regulate dorsoventral patterning, differentiation and axon guidance by secreting morphogens. Here we show that the bHLH transcription factor Nato3 (Ferd3l) is specifically expressed in the spinal FP of chick and mouse embryos. Using in ovo electroporation to understand the regulation of the FP-specific expression of Nato3, we have identified an evolutionarily conserved 204 bp genomic region, which is necessary and sufficient to drive expression to the chick FP. This promoter contains two Foxa2-binding sites, which are highly conserved among distant phyla. The two sites can bind Foxa2 in vitro, and are necessary for the expression in the FP in vivo. Gain and loss of Foxa2 function in vivo further emphasize its role in Nato3 promoter activity. Thus, our data suggest that Nato3 is a direct target of Foxa2, a transcription activator and effector of Sonic hedgehog, the hallmark regulator of FP induction and spinal cord development. The identification of the FP-specific promoter is an important step towards a better understanding of the molecular mechanisms through which Nato3 transcription is regulated and for uncovering its function during nervous system development. Moreover, the promoter provides us with a powerful tool for conditional genetic manipulations in the FP.