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

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

Construction of Light-Responsive Gene Regulatory Network for Growth, Development and Secondary Metabolite Production in Cordyceps militaris

Cordyceps militaris is an edible fungus that produces many beneficial compounds, including cordycepin and carotenoid. In many fungi, growth, development and secondary metabolite production are controlled by crosstalk between light-signaling pathways and other regulatory cascades. However, little is known about the gene regulation upon light exposure in C. militaris. This study aims to construct a gene regulatory network (GRN) that responds to light in C. militaris. First, a genome-scale GRN was built based on transcription factor (TF)-target gene interactions predicted from the Regulatory Sequence Analysis Tools (RSAT). Then, a light-responsive GRN was extracted by integrating the transcriptomic data onto the genome-scale GRN. The light-responsive network contains 2689 genes and 6837 interactions. From the network, five TFs, Snf21 (CCM_04586), an AT-hook DNA-binding motif TF (CCM_08536), a homeobox TF (CCM_07504), a forkhead box protein L2 (CCM_02646) and a heat shock factor Hsf1 (CCM_05142), were identified as key regulators that co-regulate a large group of growth and developmental genes. The identified regulatory network and expression profiles from our analysis suggested how light may induce the growth and development of C. militaris into a sexual cycle. The light-mediated regulation also couples fungal development with cordycepin and carotenoid production. This study leads to an enhanced understanding of the light-responsive regulation of growth, development and secondary metabolite production in the fungi.

SIX1 transcription factor: A review of cellular functions and regulatory dynamics

Sine Oculis Homeobox 1 (SIX1) is a member of homeobox transcription factor family having pivotal roles in organismal development and differentiation. This protein functionally acts to regulate the expression of different proteins that are involved in organ development during embryogenesis and in disorders like cancer. Aberrant expression of this homeoprotein has therefore been reported in multiple pathological complexities like hearing impairment and renal anomalies during development and tumorigenesis in adult life. Most of the cellular effects mediated by it are mostly due to its role as a transcription factor. This review presents a concise narrative of its structure, interaction partners and cellular functions vis a vis its role in cancer. We thoroughly discuss the reported molecular mechanisms that govern its function in cellular milieu. Its post-translational regulation by phosphorylation and ubiquitination are also discussed with an emphasis on yet to be explored mechanistic insights regulating its molecular dynamics to fully comprehend its role in development and disease.

An Evolutionary Perspective on Hox Binding Site Preferences in Two Different Tissues

Transcription factor (TF) networks define the precise development of multicellular organisms. While many studies focused on TFs expressed in specific cell types to elucidate their contribution to cell specification and differentiation, it is less understood how broadly expressed TFs perform their precise functions in the different cellular contexts. To uncover differences that could explain tissue-specific functions of such TFs, we analyzed here genomic chromatin interactions of the broadly expressed Drosophila Hox TF Ultrabithorax (Ubx) in the mesodermal and neuronal tissues using bioinformatics. Our investigations showed that Ubx preferentially interacts with multiple yet tissue-specific chromatin sites in putative regulatory regions of genes in both tissues. Importantly, we found the classical Hox/Ubx DNA binding motif to be enriched only among the neuronal Ubx chromatin interactions, whereas a novel Ubx-like motif with rather low predicted Hox affinities was identified among the regions bound by Ubx in the mesoderm. Finally, our analysis revealed that tissues-specific Ubx chromatin sites are also different with regards to the distribution of active and repressive histone marks. Based on our data, we propose that the tissue-related differences in Ubx binding behavior could be a result of the emergence of the mesoderm as a new germ layer in triploblastic animals, which might have required the Hox TFs to relax their binding specificity.

Molecular Characterization of HOXA2 and HOXA3 Binding Properties

The highly conserved HOX homeodomain (HD) transcription factors (TFs) establish the identity of different body parts along the antero–posterior axis of bilaterian animals. Segment diversification and the morphogenesis of different structures is achieved by generating precise patterns of HOX expression along the antero–posterior axis and by the ability of different HOX TFs to instruct unique and specific transcriptional programs. However, HOX binding properties in vitro, characterised by the recognition of similar AT-rich binding sequences, do not account for the ability of different HOX to instruct segment-specific transcriptional programs. To address this problem, we previously compared HOXA2 and HOXA3 binding in vivo. Here, we explore if sequence motif enrichments observed in vivo are explained by binding affinities in vitro. Unexpectedly, we found that the highest enriched motif in HOXA2 peaks was not recognised by HOXA2 in vitro, highlighting the importance of investigating HOX binding in its physiological context. We also report the ability of HOXA2 and HOXA3 to heterodimerise, which may have functional consequences for the HOX patterning function in vivo.

Inferring primase-DNA specific recognition using a data driven approach

DNA–protein interactions play essential roles in all living cells. Understanding of how features embedded in the DNA sequence affect specific interactions with proteins is both challenging and important, since it may contribute to finding the means to regulate metabolic pathways involving DNA–protein interactions. Using a massive experimental benchmark dataset of binding scores for DNA sequences and a machine learning workflow, we describe the binding to DNA of T7 primase, as a model system for specific DNA–protein interactions. Effective binding of T7 primase to its specific DNA recognition sequences triggers the formation of RNA primers that serve as Okazaki fragment start sites during DNA replication.

Gamete expression of TALE class HD genes activates the diploid sporophyte program in Marchantia polymorpha

Eukaryotic life cycles alternate between haploid and diploid phases and in phylogenetically diverse unicellular eukaryotes, expression of paralogous homeodomain genes in gametes primes the haploid-to-diploid transition. In the unicellular chlorophyte alga Chlamydomonas, KNOX and BELL TALE-homeodomain genes mediate this transition. We demonstrate that in the liverwort Marchantia polymorpha, paternal (sperm) expression of three of five phylogenetically diverse BELL genes, MpBELL234, and maternal (egg) expression of both MpKNOX1 and MpBELL34 mediate the haploid-to-diploid transition. Loss-of-function alleles of MpKNOX1 result in zygotic arrest, whereas a loss of either maternal or paternal MpBELL234 results in variable zygotic and early embryonic arrest. Expression of MpKNOX1 and MpBELL34 during diploid sporophyte development is consistent with a later role for these genes in patterning the sporophyte. These results indicate that the ancestral mechanism to activate diploid gene expression was retained in early diverging land plants and subsequently co-opted during evolution of the diploid sporophyte body.

Anterior Hox Genes and the Process of Cephalization

During evolution, bilateral animals have experienced a progressive process of cephalization with the anterior concentration of nervous tissue, sensory organs and the appearance of dedicated feeding structures surrounding the mouth. Cephalization has been achieved by the specialization of the unsegmented anterior end of the body (the acron) and the sequential recruitment to the head of adjacent anterior segments. Here we review the key developmental contribution of Hox1-5 genes to the formation of cephalic structures in vertebrates and arthropods and discuss how this evolved. The appearance of Hox cephalic genes preceded the evolution of a highly specialized head in both groups, indicating that Hox gene involvement in the control of cephalic structures was acquired independently during the evolution of vertebrates and invertebrates to regulate the genes required for head innovation.

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

Structure of the unique tetrameric STENOFOLIA homeodomain bound with target promoter DNA

Homeobox transcription factors are key regulators of morphogenesis and development in both animals and plants. In plants, the WUSCHEL-related homeobox (WOX) family of transcription factors function as central organizers of several developmental programs ranging from embryo patterning to meristematic stem-cell maintenance through transcriptional activation and repression mechanisms. The Medicago truncatula STENOFOLIA (STF) gene is a master regulator of leaf-blade lateral development. Here, the crystal structure of the homeodomain (HD) of STF (STF-HD) in complex with its promoter DNA is reported at 2.1 Å resolution. STF-HD binds DNA as a tetramer, enclosing nearly the entire bound DNA surface. The STF-HD tetramer is partially stabilized by docking of the C-terminal tail of one protomer onto a conserved hydrophobic surface on the head of another protomer in a head-to-tail manner. STF-HD specifically binds TGA motifs, although the promoter sequence also contains TAAT motifs. Helix α3 not only serves a canonical role as a base reader in the major groove, but also provides DNA binding in the minor groove through basic residues located at its C-terminus. The structural and functional data in planta reported here provide new insights into the DNA-binding mechanisms of plant-specific HDs from the WOX family of transcription factors.

Human model of IRX5 mutations reveals key role for this transcription factor in ventricular conduction

AIMS

Several inherited arrhythmic diseases have been linked to single gene mutations in cardiac ion channels and interacting proteins. However, the mechanisms underlying most arrhythmias, are thought to involve altered regulation of the expression of multiple effectors. In this study, we aimed to examine the role of a transcription factor (TF) belonging to the Iroquois homeobox family, IRX5, in cardiac electrical function.

METHODS AND RESULTS

Using human cardiac tissues, transcriptomic correlative analyses between IRX5 and genes involved in cardiac electrical activity showed that in human ventricular compartment, IRX5 expression strongly correlated to the expression of major actors of cardiac conduction, including the sodium channel, Nav1.5, and Connexin 40 (Cx40). We then generated human-induced pluripotent stem cells (hiPSCs) derived from two Hamamy syndrome-affected patients carrying distinct homozygous loss-of-function mutations in IRX5 gene. Cardiomyocytes derived from these hiPSCs showed impaired cardiac gene expression programme, including misregulation in the control of Nav1.5 and Cx40 expression. In accordance with the prolonged QRS interval observed in Hamamy syndrome patients, a slower ventricular action potential depolarization due to sodium current reduction was observed on electrophysiological analyses performed on patient-derived cardiomyocytes, confirming the functional role of IRX5 in electrical conduction. Finally, a cardiac TF complex was newly identified, composed by IRX5 and GATA4, in which IRX5 potentiated GATA4-induction of SCN5A expression.

CONCLUSION

Altogether, this work unveils a key role for IRX5 in the regulation of human ventricular depolarization and cardiac electrical conduction, providing therefore new insights into our understanding of cardiac diseases.

Competition for DNA binding between paralogous transcription factors determines their genomic occupancy and regulatory functions

Most eukaryotic transcription factors (TFs) are part of large protein families, with members of the same family (i.e., paralogous TFs) recognizing similar DNA-binding motifs but performing different regulatory functions. Many TF paralogs are coexpressed in the cell and thus can compete for target sites across the genome. However, this competition is rarely taken into account when studying the in vivo binding patterns of eukaryotic TFs. Here, we show that direct competition for DNA binding between TF paralogs is a major determinant of their genomic binding patterns. Using yeast proteins Cbf1 and Pho4 as our model system, we designed a high-throughput quantitative assay to capture the genomic binding profiles of competing TFs in a cell-free system. Our data show that Cbf1 and Pho4 greatly influence each other's occupancy by competing for their common putative genomic binding sites. The competition is different at different genomic sites, as dictated by the TFs' expression levels and their divergence in DNA-binding specificity and affinity. Analyses of ChIP-seq data show that the biophysical rules that dictate the competitive TF binding patterns in vitro are also followed in vivo, in the complex cellular environment. Furthermore, the Cbf1-Pho4 competition for genomic sites, as characterized in vitro using our new assay, plays a critical role in the specific activation of their target genes in the cell. Overall, our study highlights the importance of direct TF-TF competition for genomic binding and gene regulation by TF paralogs, and proposes an approach for studying this competition in a quantitative and high-throughput manner.

The Hox protein conundrum: The "specifics" of DNA binding for Hox proteins and their partners

Homeotic genes (Hox genes) are homeodomain-transcription factors involved in conferring segmental identity along the anterior-posterior body axis. Molecular characterization of HOX protein function raises some interesting questions regarding the source of the binding specificity of the HOX proteins. How do HOX proteins regulate common and unique target specificity across space and time? This review attempts to summarize and interpret findings in this area, largely focused on results from in vitro and in vivo studies in Drosophila and mouse systems. Recent studies related to HOX protein binding specificity compel us to reconsider some of our current models for transcription factor-DNA interactions. It is crucial to study transcription factor binding by incorporating components of more complex, multi-protein interactions in concert with small changes in binding motifs that can significantly impact DNA binding specificity and subsequent alterations in gene expression. To incorporate the multiple elements that can determine HOX protein binding specificity, we propose a more integrative Cooperative Binding model.

Identification of in vivo Hox13-binding sites reveals an essential locus controlling zebrafish brachyury expression

During early embryogenesis, the vertebrate embryo extends from anterior to posterior because of the progressive addition of cells from a posteriorly localized neuromesodermal progenitor (NMp) population. An autoregulatory loop between Wnt and Brachyury/Tbxt is required for NMps to retain mesodermal potential and, hence, normal axis development. We recently showed that Hox13 genes help to support body axis formation and to maintain the autoregulatory loop, although the direct Hox13 target genes were unknown. Here, using a new method for identifying in vivo transcription factor-binding sites, we identified more than 500 potential Hox13 target genes in zebrafish. Importantly, we found two highly conserved Hox13-binding elements far from the tbxta transcription start site that also contain a conserved Tcf7/Lef1 (Wnt response) site. We show that the proximal of the two elements is sufficient to confer somitogenesis-stage expression to a tbxta promoter that, on its own, only drives NMp expression during gastrulation. Importantly, elimination of this proximal element produces shortened embryos due to aberrant formation of the most posterior somites. Our study provides a potential direct connection between Hox13 and regulation of the Wnt/Brachyury loop.

A Computer-Based Methodology to Design Non-Standard Peptides Potentially Able to Prevent HOX-PBX1-Associated Cancer Diseases

In the last decades, HOX proteins have been extensively studied due to their pivotal role in transcriptional events. HOX proteins execute their activity by exploiting a cooperative binding to PBX proteins and DNA. Therefore, an increase or decrease in HOX activity has been associated with both solid and haematological cancer diseases. Thus, inhibiting HOX-PBX interaction represents a potential strategy to prevent these malignancies, as demonstrated by the patented peptide HTL001 that is being studied in clinical trials. In this work, a computational study is described to identify novel potential peptides designed by employing a database of non-natural amino acids. For this purpose, residue scanning of the HOX minimal active sequence was performed to select the mutations to be further processed. According to these results, the peptides were point-mutated and used for Molecular Dynamics (MD) simulations in complex with PBX1 protein and DNA to evaluate complex binding stability. MM-GBSA calculations of the resulting MD trajectories were exploited to guide the selection of the most promising mutations that were exploited to generate twelve combinatorial peptides. Finally, the latter peptides in complex with PBX1 protein and DNA were exploited to run MD simulations and the ΔGbinding average values of the complexes were calculated. Thus, the analysis of the results highlighted eleven combinatorial peptides that will be considered for further assays.

Tlx3 Exerts Direct Control in Specifying Excitatory Over Inhibitory Neurons in the Dorsal Spinal Cord

The spinal cord dorsal horn is a major station for integration and relay of somatosensory information and comprises both excitatory and inhibitory neuronal populations. The homeobox gene Tlx3 acts as a selector gene to control the development of late-born excitatory (dILB) neurons by specifying glutamatergic transmitter fate in dorsal spinal cord. However, since Tlx3 direct transcriptional targets remain largely unknown, it remains to be uncovered how Tlx3 functions to promote excitatory cell fate. Here we combined a genomics approach based on chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq) and expression profiling, with validation experiments in Tlx3 null embryos, to characterize the transcriptional program of Tlx3 in mouse embryonic dorsal spinal cord. We found most dILB neuron specific genes previously identified to be directly activated by Tlx3. Surprisingly, we found Tlx3 also directly represses many genes associated with the alternative inhibitory dILA neuronal fate. In both cases, direct targets include transcription factors and terminal differentiation genes, showing that Tlx3 directly controls cell identity at distinct levels. Our findings provide a molecular frame for the master regulatory role of Tlx3 in developing glutamatergic dILB neurons. In addition, they suggest a novel function for Tlx3 as direct repressor of GABAergic dILA identity, pointing to how generation of the two alternative cell fates being tightly coupled.

Blueprint of human thymopoiesis reveals molecular mechanisms of stage-specific TCR enhancer activation

Cell differentiation is accompanied by epigenetic changes leading to precise lineage definition and cell identity. Here we present a comprehensive resource of epigenomic data of human T cell precursors along with an integrative analysis of other hematopoietic populations. Although T cell commitment is accompanied by large scale epigenetic changes, we observed that the majority of distal regulatory elements are constitutively unmethylated throughout T cell differentiation, irrespective of their activation status. Among these, the TCRA gene enhancer (Eα) is in an open and unmethylated chromatin structure well before activation. Integrative analyses revealed that the HOXA5-9 transcription factors repress the Eα enhancer at early stages of T cell differentiation, while their decommission is required for TCRA locus activation and enforced αβ T lineage differentiation. Remarkably, the HOXA-mediated repression of Eα is paralleled by the ectopic expression of homeodomain-related oncogenes in T cell acute lymphoblastic leukemia. These results highlight an analogous enhancer repression mechanism at play in normal and cancer conditions, but imposing distinct developmental constraints.

Genome-Wide Binding Analyses of HOXB1 Revealed a Novel DNA Binding Motif Associated with Gene Repression

Knowledge of the diverse DNA binding specificities of transcription factors is important for understanding their specific regulatory functions in animal development and evolution. We have examined the genome-wide binding properties of the mouse HOXB1 protein in embryonic stem cells differentiated into neural fates. Unexpectedly, only a small number of HOXB1 bound regions (7%) correlate with binding of the known HOX cofactors PBX and MEIS. In contrast, 22% of the HOXB1 binding peaks display co-occupancy with the transcriptional repressor REST. Analyses revealed that co-binding of HOXB1 with PBX correlates with active histone marks and high levels of expression, while co-occupancy with REST correlates with repressive histone marks and repression of the target genes. Analysis of HOXB1 bound regions uncovered enrichment of a novel 15 base pair HOXB1 binding motif HB1RE (HOXB1 response element). In vitro template binding assays showed that HOXB1, PBX1, and MEIS can bind to this motif. In vivo, this motif is sufficient for direct expression of a reporter gene and over-expression of HOXB1 selectively represses this activity. Our analyses suggest that HOXB1 has evolved an association with REST in gene regulation and the novel HB1RE motif contributes to HOXB1 function in part through a repressive role in gene expression.

Contrasting DNA-binding behaviour by ISL1 and LHX3 underpins differential gene targeting in neuronal cell specification

•

The mechanisms by which ISL1 and LHX3 specify neuronal cell identity are unknown.

•

EMSA/SPR data show ISL1 and LHX3 have markedly different DNA-binding behaviours.

•

SAXS shows ISL1/LHX3:DNA complexes are flexible in nature.

•

ISL1 binds DNA poorly but appears to modulate the DNA-binding specificity of LHX3.

The roles of ISL1 and LHX3 in the development of spinal motor neurons have been well established. Whereas LHX3 triggers differentiation into interneurons, the additional expression of ISL1 in developing neuronal cells is sufficient to redirect their developmental trajectory towards spinal motor neurons. However, the underlying mechanism of this action by these transcription factors is less well understood. Here, we used electrophoretic mobility shift assays (EMSAs) and surface plasmon resonance (SPR) to probe the different DNA-binding behaviours of these two proteins, both alone and in complexes mimicking those found in developing neurons, and found that ISL1 shows markedly different binding properties to LHX3. We used small angle X-ray scattering (SAXS) to structurally characterise DNA-bound species containing ISL1 and LHX3. Taken together, these results have allowed us to develop a model of how these two DNA-binding modules coordinate to regulate gene expression and direct development of spinal motor neurons.

Conserved Gsx2/Ind homeodomain monomer versus homodimer DNA binding defines regulatory outcomes in flies and mice

How homeodomain proteins gain sufficient specificity to control different cell fates has been a long-standing problem in developmental biology. The conserved Gsx homeodomain proteins regulate specific aspects of neural development in animals from flies to mammals, and yet they belong to a large transcription factor family that bind nearly identical DNA sequences in vitro. Here, we show that the mouse and fly Gsx factors unexpectedly gain DNA binding specificity by forming cooperative homodimers on precisely spaced and oriented DNA sites. High-resolution genomic binding assays revealed that Gsx2 binds both monomer and homodimer sites in the developing mouse ventral telencephalon. Importantly, reporter assays showed that Gsx2 mediates opposing outcomes in a DNA binding site-dependent manner: Monomer Gsx2 binding represses transcription, whereas homodimer binding stimulates gene expression. In Drosophila, the Gsx homolog, Ind, similarly represses or stimulates transcription in a site-dependent manner via an autoregulatory enhancer containing a combination of monomer and homodimer sites. Integrating these findings, we test a model showing how the homodimer to monomer site ratio and the Gsx protein levels defines gene up-regulation versus down-regulation. Altogether, these data serve as a new paradigm for how cooperative homeodomain transcription factor binding can increase target specificity and alter regulatory outcomes.

HOX paralogs selectively convert binding of ubiquitous transcription factors into tissue-specific patterns of enhancer activation

Gene expression programs determine cell fate in embryonic development and their dysregulation results in disease. Transcription factors (TFs) control gene expression by binding to enhancers, but how TFs select and activate their target enhancers is still unclear. HOX TFs share conserved homeodomains with highly similar sequence recognition properties, yet they impart the identity of different animal body parts. To understand how HOX TFs control their specific transcriptional programs in vivo, we compared HOXA2 and HOXA3 binding profiles in the mouse embryo. HOXA2 and HOXA3 directly cooperate with TALE TFs and selectively target different subsets of a broad TALE chromatin platform. Binding of HOX and tissue-specific TFs convert low affinity TALE binding into high confidence, tissue-specific binding events, which bear the mark of active enhancers. We propose that HOX paralogs, alone and in combination with tissue-specific TFs, generate tissue-specific transcriptional outputs by modulating the activity of TALE TFs at selected enhancers.

Differential abilities to engage inaccessible chromatin diversify vertebrate Hox binding patterns

Although Hox genes encode for conserved transcription factors (TFs), they are further divided into anterior, central and posterior groups based on their DNA-binding domain similarity. The posterior Hox group expanded in the deuterostome clade and patterns caudal and distal structures. We aimed to address how similar Hox TFs diverge to induce different positional identities. We studied Hox TF DNA-binding and regulatory activity during an in vitro motor neuron differentiation system that recapitulates embryonic development. We found diversity in the genomic binding profiles of different Hox TFs, even among the posterior group paralogs that share similar DNA-binding domains. These differences in genomic binding were explained by differing abilities to bind to previously inaccessible sites. For example, the posterior group HOXC9 had a greater ability to bind occluded sites than the posterior HOXC10, producing different binding patterns and driving differential gene expression programs. From these results, we propose that the differential abilities of posterior Hox TFs to bind to previously inaccessible chromatin drive patterning diversification.This article has an associated 'The people behind the papers' interview.

A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

Structural basis for designing an array of engrailed homeodomains

Small DNA-binding proteins that target desired sequences have the potential to act as a scaffold for molecular tools such as genome editing. In this study, an engrailed homeodomain (EHD) was chosen and it was evaluated whether it could be used as a molecular module that can connect to itself to recognize a longer target sequence. It was previously shown that two EHDs connected by a linker (EHD₂) recognize a target sequence twice as long as that recognized by a single EHD in cells only when Arg53 in each EHD in the tandem protein is mutated to alanine {(EHD[R53A])₂}. To investigate the recognition mechanism of (EHD[R53A])₂, the crystal structure of the (EHD[R53A])₂-DNA complex was determined at 1.6 Å resolution. The individual EHDs were found to adopt the typical homeodomain fold. Most importantly, the base-specific interactions in the major groove necessary for the affinity/specificity of wild-type EHD were preserved in (EHD[R53A])₂. Bacterial assays confirmed that the base-specific interactions are retained under cellular conditions. These observations indicate that the R53A mutation only causes a loss of the arginine-phosphate interaction at the protein-DNA interface, which reduces the DNA-binding affinity compared with the wild type. It is therefore concluded that (EHD[R53A])₂ precisely recognizes tandem target sites within cells, enabling the individual EHDs to concurrently bind to the target sites with modest binding affinity. This suggests that modulation of the binding activity of each EHD is vital to construct a protein array that can precisely recognize a sequence with multiple target sites.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

The geometric influence on the Cys2His2 zinc finger domain and functional plasticity

The Cys₂His₂ zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein–DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.

Constraint and conservation of paired-type homeodomains predicts the clinical outcome of missense variants of uncertain significance

The need to interpret the pathogenicity of novel missense variants of unknown significance identified in the homeodomain of X-chromosome aristaless-related homeobox (ARX) gene prompted us to assess the utility of conservation and constraint across these domains in multiple genes compared to conventional in vitro functional analysis. Pathogenic missense variants clustered in the homeodomain of ARX contribute to intellectual disability (ID) and epilepsy, with and without brain malformation in affected males. Here we report novel c.1112G>A, p.Arg371Gln and c.1150C>T, p.Arg384Cys variants in male patients with ID and severe seizures. The third case of a male patient with a c.1109C>T, p.Ala370Val variant is perhaps the first example of ID and autism spectrum disorder (ASD), without seizures or brain malformation. We compiled data sets of pathogenic variants from ClinVar and presumed benign variation from gnomAD and demonstrated that the high levels of sequence conservation and constraint of benign variation within the homeodomain impacts upon the ability of publicly available in silico prediction tools to accurately discern likely benign from likely pathogenic variants in these data sets. Despite this, considering the inheritance patterns of the genes and disease variants with the conservation and constraint of disease variants affecting the homeodomain in conjunction with current clinical assessments may assist in predicting the pathogenicity of missense variants, particularly for genes with autosomal recessive and X-linked patterns of disease inheritance, such as ARX. In vitro functional analysis demonstrates that the transcriptional activity of all three variants was diminished compared to ARX-Wt. We review the associated phenotypes of the published cases of patients with ARX homeodomain variants and propose expansion of the ARX-related phenotype to include severe ID and ASD without brain malformations or seizures. We propose that the use of the constraint and conservation data in conjunction with consideration of the patient phenotype and inheritance pattern may negate the need for the experimental functional validation currently required to achieve a diagnosis.

I-Block: a simple Escherichia coli-based assay for studying sequence-specific DNA binding of proteins

We have developed a simple method called I-Block assay, which can detect sequence-specific binding of proteins to DNA in Escherichia coli. The method works by detecting competition between the protein of interest and RNA polymerase for binding to overlapping target sites in a plasmid-borne lacI promoter variant. The assay utilizes two plasmids and an E. coli host strain, from which the gene of the Lac repressor (lacI) has been deleted. One of the plasmids carries the lacI gene with a unique NheI restriction site created in the lacI promoter. The potential recognition sequences of the tested protein are inserted into the NheI site. Introduction of the plasmids into the E. coliΔlacI host represses the constitutive β-galactosidase synthesis of the host bacterium. If the studied protein expressed from a compatible plasmid binds to its target site in the lacI promoter, it will interfere with lacI transcription and lead to increased β-galactosidase activity. The method was tested with two zinc finger proteins, with the lambda phage cI857 repressor, and with CRISPR-dCas9 targeted to the lacI promoter. The I-Block assay was shown to work with standard liquid cultures, with cultures grown in microplate and with colonies on X-gal indicator plates.

Structural basis for the complex DNA binding behavior of the plant stem cell regulator WUSCHEL

Stem cells are one of the foundational evolutionary novelties that allowed the independent emergence of multicellularity in the plant and animal lineages. In plants, the homeodomain (HD) transcription factor WUSCHEL (WUS) is essential for the maintenance of stem cells in the shoot apical meristem. WUS has been reported to bind to diverse DNA motifs and to act as transcriptional activator and repressor. However, the mechanisms underlying this remarkable behavior have remained unclear. Here, we quantitatively delineate WUS binding to three divergent DNA motifs and resolve the relevant structural underpinnings. We show that WUS exhibits a strong binding preference for TGAA repeat sequences, while retaining the ability to weakly bind to TAAT elements. This behavior is attributable to the formation of dimers through interactions of specific residues in the HD that stabilize WUS DNA interaction. Our results provide a mechanistic basis for dissecting WUS dependent regulatory networks in plant stem cell control.

Developmental regulation of cell type-specific transcription by novel promoter-proximal sequence elements

Cell type-specific transcriptional programs that drive differentiation of specialized cell types are key players in development and tissue regeneration. One of the most dramatic changes in the transcription program in Drosophila occurs with the transition from proliferating spermatogonia to differentiating spermatocytes, with >3000 genes either newly expressed or expressed from new alternative promoters in spermatocytes. Here we show that opening of these promoters from their closed state in precursor cells requires function of the spermatocyte-specific tMAC complex, localized at the promoters. The spermatocyte-specific promoters lack the previously identified canonical core promoter elements except for the Inr. Instead, these promoters are enriched for the binding site for the TALE-class homeodomain transcription factors Achi/Vis and for a motif originally identified under tMAC ChIP-seq peaks. The tMAC motif resembles part of the previously identified 14-bp β2UE1 element critical for spermatocyte-specific expression. Analysis of downstream sequences relative to transcription start site usage suggested that ACA and CNAAATT motifs at specific positions can help promote efficient transcription initiation. Our results reveal how promoter-proximal sequence elements that recruit and are acted upon by cell type-specific chromatin binding complexes help establish a robust, cell type-specific transcription program for terminal differentiation.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Discovering human transcription factor physical interactions with genetic variants, novel DNA motifs, and repetitive elements using enhanced yeast one-hybrid assays

Identifying transcription factor (TF) binding to noncoding variants, uncharacterized DNA motifs, and repetitive genomic elements has been technically and computationally challenging. Current experimental methods, such as chromatin immunoprecipitation, generally test one TF at a time, and computational motif algorithms often lead to false-positive and -negative predictions. To address these limitations, we developed an experimental approach based on enhanced yeast one-hybrid assays. The first variation of this approach interrogates the binding of >1000 human TFs to repetitive DNA elements, while the second evaluates TF binding to single nucleotide variants, short insertions and deletions (indels), and novel DNA motifs. Using this approach, we detected the binding of 75 TFs, including several nuclear hormone receptors and ETS factors, to the highly repetitive Alu elements. Further, we identified cancer-associated changes in TF binding, including gain of interactions involving ETS TFs and loss of interactions involving KLF TFs to different mutations in the TERT promoter, and gain of a MYB interaction with an 18-bp indel in the TAL1 superenhancer. Additionally, we identified TFs that bind to three uncharacterized DNA motifs identified in DNase footprinting assays. We anticipate that these enhanced yeast one-hybrid approaches will expand our capabilities to study genetic variation and undercharacterized genomic regions.

Genome-wide analysis of HOXC4 and HOXC6 regulated genes and binding sites in prostate cancer cells

Aberrant expression of HOXC6 and HOXC4 is commonly detected in prostate cancer. The high expression of these transcription factors is associated with aggressive prostate cancer and can predict cancer recurrence after treatment. Thus, HOXC4 and HOXC6 are clinically relevant biomarkers of aggressive prostate cancer. However, the molecular mechanisms by which these HOXC genes contribute to prostate cancer is not yet understood. To begin to address the role of HOXC4 and HOXC6 in prostate cancer, we performed RNA-seq analyses before and after siRNA-mediated knockdown of HOXC4 and/or HOXC6 and also performed ChIP-seq to identify genomic binding sites for both of these transcription factors. Our studies demonstrate that HOXC4 and HOXC6 co-localize with HOXB13, FOXA1 and AR, three transcription factors previously shown to contribute to the development of prostate cancer. We suggest that the aberrantly upregulated HOXC4 and HOXC6 proteins may compete with HOXB13 for binding sites, thus altering the prostate transcriptome. This competition model may be applicable to many different human cancers that display increased expression of a HOX transcription factor.

Sharing DNA-binding information across structurally similar proteins enables accurate specificity determination

We are now in an era where protein-DNA interactions have been experimentally assayed for thousands of DNA-binding proteins. In order to infer DNA-binding specificities from these data, numerous sophisticated computational methods have been developed. These approaches typically infer DNA-binding specificities by considering interactions for each protein independently, ignoring related and potentially valuable interaction information across other proteins that bind DNA via the same structural domain. Here we introduce a framework for inferring DNA-binding specificities by considering protein-DNA interactions for entire groups of structurally similar proteins simultaneously. We devise both constrained optimization and label propagation algorithms for this task, each balancing observations at the individual protein level against dataset-wide consistency of interaction preferences. We test our approaches on two large, independent Cys2His2 zinc finger protein-DNA interaction datasets. We demonstrate that jointly inferring specificities within each dataset individually dramatically improves accuracy, leading to increased agreement both between these two datasets and with a fixed external standard. Overall, our results suggest that sharing protein-DNA interaction information across structurally similar proteins is a powerful means to enable accurate inference of DNA-binding specificities.

Cis-regulatory basis of sister cell type divergence in the vertebrate retina

Multicellular organisms evolved via repeated functional divergence of transcriptionally related sister cell types, but the mechanisms underlying sister cell type divergence are not well understood. Here, we study a canonical pair of sister cell types, retinal photoreceptors and bipolar cells, to identify the key cis-regulatory features that distinguish them. By comparing open chromatin maps and transcriptomic profiles, we found that while photoreceptor and bipolar cells have divergent transcriptomes, they share remarkably similar cis-regulatory grammars, marked by enrichment of K50 homeodomain binding sites. However, cell class-specific enhancers are distinguished by enrichment of E-box motifs in bipolar cells, and Q50 homeodomain motifs in photoreceptors. We show that converting K50 motifs to Q50 motifs represses reporter expression in bipolar cells, while photoreceptor expression is maintained. These findings suggest that partitioning of Q50 motifs within cell type-specific cis-regulatory elements was a critical step in the evolutionary divergence of the bipolar transcriptome from that of photoreceptors.

Systematic domain-based aggregation of protein structures highlights DNA-, RNA- and other ligand-binding positions

Domains are fundamental subunits of proteins, and while they play major roles in facilitating protein–DNA, protein–RNA and other protein–ligand interactions, a systematic assessment of their various interaction modes is still lacking. A comprehensive resource identifying positions within domains that tend to interact with nucleic acids, small molecules and other ligands would expand our knowledge of domain functionality as well as aid in detecting ligand-binding sites within structurally uncharacterized proteins. Here, we introduce an approach to identify per-domain-position interaction ‘frequencies’ by aggregating protein co-complex structures by domain and ascertaining how often residues mapping to each domain position interact with ligands. We perform this domain-based analysis on ∼91000 co-complex structures, and infer positions involved in binding DNA, RNA, peptides, ions or small molecules across 4128 domains, which we refer to collectively as the InteracDome. Cross-validation testing reveals that ligand-binding positions for 2152 domains are highly consistent and can be used to identify residues facilitating interactions in ∼63–69% of human genes. Our resource of domain-inferred ligand-binding sites should be a great aid in understanding disease etiology: whereas these sites are enriched in Mendelian-associated and cancer somatic mutations, they are depleted in polymorphisms observed across healthy populations. The InteracDome is available at http://interacdome.princeton.edu.

In vivo Hox binding specificity revealed by systematic changes to a single cis regulatory module

Hox proteins belong to a family of transcription factors with similar DNA binding specificities that control animal differentiation along the antero-posterior body axis. Hox proteins are expressed in partially overlapping regions where each one is responsible for the formation of particular organs and structures through the regulation of specific direct downstream targets. Thus, explaining how each Hox protein can selectively control its direct targets from those of another Hox protein is fundamental to understand animal development. Here we analyse a cis regulatory module directly regulated by seven different Drosophila Hox proteins and uncover how different Hox class proteins differentially control its expression. We find that regulation by one or another Hox protein depends on the combination of three modes: Hox-cofactor dependent DNA-binding specificity; Hox-monomer binding sites; and interaction with positive and negative Hox-collaborator proteins. Additionally, we find that similar regulation can be achieved by Amphioxus orthologs, suggesting these three mechanisms are conserved from insects to chordates.

Hox proteins are expressed in partially overlapping regions to inform development along the embryo’s head-tail axis. Here the authors analyse a cis regulatory module directly regulated by seven different Drosophila Hox proteins to uncover how different Hox class proteins differentially control its expression.

Balance between DNA-binding affinity and specificity enables selective recognition of longer target sequences in vivo

Although genome-editing enzymes such as TALEN and CRISPR/Cas9 are being widely used, they have an essential limitation in that their relatively high-molecular weight makes them difficult to be delivered to cells. To develop a novel genome-editing enzyme with a smaller molecular weight, we focused on the engrailed homeodomain (EHD). We designed and constructed proteins composed of two EHDs connected by a linker to increase sequence specificity. In bacterial one-hybrid assays and electrophoresis mobility shift assay analyses, the created proteins exhibited good affinity for DNA sequences consisting of two tandemly aligned EHD target sequences. However, they also bound to individual EHD targets. To avoid binding to single target sites, we introduced amino acid mutations to reduce the protein-DNA affinity of each EHD monomer and successfully created a small protein with high specificity for tandem EHD target sequences.

Multi-enhancer transcriptional hubs confer phenotypic robustness

We previously showed in Drosophila melanogaster embryos that low-affinity Ultrabithorax (Ubx)-responsive shavenbaby (svb) enhancers drive expression using localized transcriptional environments and that active svb enhancers on different chromosomes tended to colocalize (Tsai et al., 2017). Here, we test the hypothesis that these multi-enhancer 'hubs' improve phenotypic resilience to stress by buffering against decreases in transcription factor concentrations and transcriptional output. Deleting a redundant enhancer from the svb locus led to reduced trichome numbers in embryos raised at elevated temperatures. Using high-resolution fluorescence microscopy, we observed lower Ubx concentration and transcriptional output in this deletion allele. Transcription sites of the full svb cis-regulatory region inserted into a different chromosome colocalized with the svb locus, increasing Ubx concentration, the transcriptional output of svb, and partially rescuing the phenotype. Thus, multiple enhancers could reinforce a local transcriptional hub to buffer against environmental stresses and genetic perturbations, providing a mechanism for phenotypical robustness.

Novel structural features drive DNA binding properties of Cmr, a CRP family protein in TB complex mycobacteria

Mycobacterium tuberculosis (Mtb) encodes two CRP/FNR family transcription factors (TF) that contribute to virulence, Cmr (Rv1675c) and CRP_Mt (Rv3676). Prior studies identified distinct chromosomal binding profiles for each TF despite their recognizing overlapping DNA motifs. The present study shows that Cmr binding specificity is determined by discriminator nucleotides at motif positions 4 and 13. X-ray crystallography and targeted mutational analyses identified an arginine-rich loop that expands Cmr’s DNA interactions beyond the classical helix-turn-helix contacts common to all CRP/FNR family members and facilitates binding to imperfect DNA sequences. Cmr binding to DNA results in a pronounced asymmetric bending of the DNA and its high level of cooperativity is consistent with DNA-facilitated dimerization. A unique N-terminal extension inserts between the DNA binding and dimerization domains, partially occluding the site where the canonical cAMP binding pocket is found. However, an unstructured region of this N-terminus may help modulate Cmr activity in response to cellular signals. Cmr’s multiple levels of DNA interaction likely enhance its ability to integrate diverse gene regulatory signals, while its novel structural features establish Cmr as an atypical CRP/FNR family member.

Low-Affinity Binding Sites and the Transcription Factor Specificity Paradox in Eukaryotes

Eukaryotic transcription factors (TFs) from the same structural family tend to bind similar DNA sequences, despite the ability of these TFs to execute distinct functions in vivo. The cell partly resolves this specificity paradox through combinatorial strategies and the use of low-affinity binding sites, which are better able to distinguish between similar TFs. However, because these sites have low affinity, it is challenging to understand how TFs recognize them in vivo. Here, we summarize recent findings and technological advancements that allow for the quantification and mechanistic interpretation of TF recognition across a wide range of affinities. We propose a model that integrates insights from the fields of genetics and cell biology to provide further conceptual understanding of TF binding specificity. We argue that in eukaryotes, target specificity is driven by an inhomogeneous 3D nuclear distribution of TFs and by variation in DNA binding affinity such that locally elevated TF concentration allows low-affinity binding sites to be functional.

Dissecting the sharp response of a canonical developmental enhancer reveals multiple sources of cooperativity

Developmental enhancers integrate graded concentrations of transcription factors (TFs) to create sharp gene expression boundaries. Here we examine the hunchback P2 (HbP2) enhancer which drives a sharp expression pattern in the Drosophila blastoderm embryo in response to the transcriptional activator Bicoid (Bcd). We systematically interrogate cis and trans factors that influence the shape and position of expression driven by HbP2, and find that the prevailing model, based on pairwise cooperative binding of Bcd to HbP2 is not adequate. We demonstrate that other proteins, such as pioneer factors, Mediator and histone modifiers influence the shape and position of the HbP2 expression pattern. Comparing our results to theory reveals how higher-order cooperativity and energy expenditure impact boundary location and sharpness. Our results emphasize that the bacterial view of transcription regulation, where pairwise interactions between regulatory proteins dominate, must be reexamined in animals, where multiple molecular mechanisms collaborate to shape the gene regulatory function.

Chromatin accessibility plays a key role in selective targeting of Hox proteins

BACKGROUND

Hox transcription factors specify segmental diversity along the anterior-posterior body axis in metazoans. While the different Hox family members show clear functional specificity in vivo, they all show similar binding specificity in vitro and a satisfactory understanding of in vivo Hox target selectivity is still lacking.

RESULTS

Using transient transfection in Kc167 cells, we systematically analyze the binding of all eight Drosophila Hox proteins. We find that Hox proteins show considerable binding selectivity in vivo even in the absence of canonical Hox cofactors Extradenticle and Homothorax. Hox binding selectivity is strongly associated with chromatin accessibility, being highest in less accessible chromatin. Individual Hox proteins exhibit different propensities to bind less accessible chromatin, and high binding selectivity is associated with high-affinity binding regions, leading to a model where Hox proteins derive binding selectivity through affinity-based competition with nucleosomes. Extradenticle/Homothorax cofactors generally facilitate Hox binding, promoting binding to regions in less accessible chromatin but with little effect on the overall selectivity of Hox targeting. These cofactors collaborate with Hox proteins in opening chromatin, in contrast to the pioneer factor, Glial cells missing, which facilitates Hox binding by independently generating accessible chromatin regions.

CONCLUSIONS

These studies indicate that chromatin accessibility plays a key role in Hox selectivity. We propose that relative chromatin accessibility provides a basis for subtle differences in binding specificity and affinity to generate significantly different sets of in vivo genomic targets for different Hox proteins.

Similarity regression predicts evolution of transcription factor sequence specificity

Transcription factor (TF) binding specificities (motifs) are essential for the analysis of gene regulation. Accurate prediction of TF motifs is critical, because it is infeasible to assay all TFs in all sequenced eukaryotic genomes. There is ongoing controversy regarding the degree of motif diversification among related species that is, in part, because of uncertainty in motif prediction methods. Here we describe similarity regression, a significantly improved method for predicting motifs, which we use to update and expand the Cis-BP database. Similarity regression inherently quantifies TF motif evolution, and shows that previous claims of near-complete conservation of motifs between human and Drosophila are inflated, with nearly half of the motifs in each species absent from the other, largely due to extensive divergence in C2H2 zinc finger proteins. We conclude that diversification in DNA-binding motifs is pervasive, and present a new tool and updated resource to study TF diversity and gene regulation across eukaryotes.