Dissecting the functional specificities of two Hox proteins

The intrinsically disordered transactivation region of HOXA9 regulates its function by auto-inhibition of its DNA-binding activity

HOXA9 transcription factor is expressed in hematopoietic stem cells and is involved in the regulation of their differentiation and maturation to various blood cells. HOXA9 is linked to various leukemia and is a marker for poor prognosis of acute myeloid leukemia (AML). This protein has a conserved DNA-binding homeodomain and a transactivation domain. We show that this N-terminal transactivation domain is intrinsically disordered and inhibits DNA-binding by the homeodomain. Using NMR spectroscopy and molecular dynamics simulation, we show that the hexapeptide 197AANWLH202 in the disordered region transiently occludes the DNA-binding interface. The hexapeptide also forms a rigid segment, as determined by NMR dynamics, in an otherwise flexible disordered region. Interestingly, this hexapeptide is known to mediate the interaction of HOXA9 and its TALE partner proteins, such as PBX1, and help in cooperative DNA binding. Mutation of tryptophan to alanine in the hexapeptide abrogates the DNA-binding auto-inhibition. We propose that the disordered transactivation region plays a dual role in the regulation of HOXA9 function. In the absence of TALE partners, it inhibits DNA binding, and in the presence of TALE partners it interacts with the TALE protein and facilitates the cooperative DNA binding by the HOX-TALE complex.

Stability and dynamics of extradenticle modulates its function

Extradenticle (EXD) is a partner protein of the HOX transcription factors and plays an important role in the development of Drosophila. It confers increased affinity and specificity of DNA-binding to the HOX proteins. However, the DNA-binding homeodomain of EXD has a significantly weaker affinity to DNA compared to the HOX homeodomains. Here, we show that a glycine residue (G290) in the middle of the EXD DNA-binding helix primarily results in this weaker binding. Glycine destabilizes helices. To probe its role in the stability and function of the protein, G290 was mutated to alanine. The intrinsic stability of the DNA-binding helix increased in the G290A mutant as observed by NMR studies and molecular dynamics (MD) simulation. Also, NMR dynamics and MD simulation show that dynamic motions present in the wild-type protein are quenched in the mutant. This in turn resulted in increased stability of the entire homeodomain (ΔΔGG→A of -2.6 kcal/mol). Increased protein stability resulted in three-fold better DNA-binding affinity of the mutant as compared to the wild-type protein. Molecular mechanics with generalized Born and surface area solvation (MMGBSA) analysis of our MD simulation on DNA-bound models of both wild-type and mutant proteins shows that the contribution to binding is enhanced for most of the interface residues in the mutant compared to the wild-type. Interestingly, the flexible N-terminal arm makes more stable contact with the DNA minor groove in the mutant. We found that the two interaction sites i.e. the DNA-binding helix and the unstructured N-terminal arm influence each other via the bound DNA. These results provide an interesting conundrum: alanine at position 290 enhances both the stability and the DNA-binding affinity of the protein, however, evolution prefers glycine at this position. We have provided several plausible explanations for this apparent conundrum. The function of the EXD as a HOX co-factor requires its ability to discriminate similar DNA sequences, which is most likely comprom.

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

Specificity of the Hox member Deformed is determined by transcription factor levels and binding site affinities

Hox proteins have similar binding specificities in vitro, yet they control different morphologies in vivo. This paradox has been partially solved with the identification of Hox low-affinity binding sites. However, anterior Hox proteins are more promiscuous than posterior Hox proteins, raising the question how anterior Hox proteins achieve specificity. We use the AP2x enhancer, which is activated in the maxillary head segment by the Hox TF Deformed (Dfd). This enhancer lacks canonical Dfd-Exd sites but contains several predicted low-affinity sites. Unexpectedly, these sites are strongly bound by Dfd-Exd complexes and their conversion into optimal Dfd-Exd sites results only in a modest increase in binding strength. These small variations in affinity change the sensitivity of the enhancer to different Dfd levels, resulting in perturbed AP-2 expression and maxillary morphogenesis. Thus, Hox-regulated morphogenesis seems to result from the co-evolution of Hox binding affinity and Hox dosage for precise target gene regulation.

Despite the central role of Hox genes in controlling morphogenesis, the DNA binding of different Hox members is relatively similar. Here they show that specificity of Hox member Dfd relies on a precise balance of transcription factors and binding site affinities.

Trimeric complexes of Antp-TBP with TFIIEβ or Exd modulate transcriptional activity

Background

Hox proteins finely coordinate antero-posterior axis during embryonic development and through their action specific target genes are expressed at the right time and space to determine the embryo body plan. As master transcriptional regulators, Hox proteins recognize DNA through the homeodomain (HD) and interact with a multitude of proteins, including general transcription factors and other cofactors. HD binding specificity increases by protein–protein interactions with a diversity of cofactors that outline the Hox interactome and determine the transcriptional landscape of the selected target genes. All these interactions clearly demonstrate Hox-driven transcriptional regulation, but its precise mechanism remains to be elucidated.

Results

Here we report Antennapedia (Antp) Hox protein–protein interaction with the TATA-binding protein (TBP) and the formation of novel trimeric complexes with TFIIEβ and Extradenticle (Exd), as well as its participation in transcriptional regulation. Using Bimolecular Fluorescence Complementation (BiFC), we detected the interaction of Antp-TBP and, in combination with Förster Resonance Energy Transfer (BiFC-FRET), the formation of the trimeric complex with TFIIEβ and Exd in living cells. Mutational analysis showed that Antp interacts with TBP through their N-terminal polyglutamine-stretches. The trimeric complexes of Antp-TBP with TFIIEβ and Exd were validated using different Antp mutations to disrupt the trimeric complexes. Interestingly, the trimeric complex Antp-TBP-TFIIEβ significantly increased the transcriptional activity of Antp, whereas Exd diminished its transactivation.

Conclusions

Our findings provide important insights into the Antp interactome with the direct interaction of Antp with TBP and the two new trimeric complexes with TFIIEβ and Exd. These novel interactions open the possibility to analyze promoter function and gene expression to measure transcription factor binding dynamics at target sites throughout the genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00239-8.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

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.

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.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

Secrets of secretion-How studies of the Drosophila salivary gland have informed our understanding of the cellular networks underlying secretory organ form and function

Secretory organs are critical for organismal survival. Yet, the transcriptional regulatory mechanisms governing their development and maintenance remain unclear for most model secretory organs. The Drosophila embryonic salivary gland (SG) remedies this deficiency as one of the few organs wherein direct connections from the expression of the early patterning genes to cell specification to organ architecture and functional specialization can be made. Few other models of secretion can be accorded this distinction. Studies from the past three decades have made enormous strides in parsing out the roles of distinct transcription factors (TFs) that direct major steps in furnishing this secretory organ. In the first step of specifying the salivary gland, the activity of the Hox factors Sex combs reduced, Extradenticle, and Homothorax activate expression of fork head (fkh), sage, and CrebA, which code for the major suite of TFs that carry forward the task of organ building and maintenance. Then, in the second key step of building the SG, the program for cell fate maintenance and morphogenesis is deployed. Fkh maintains the secretory cell fate by regulating its own expression and that of sage and CrebA. Fkh and Sage maintain secretory cell viability by actively blocking apoptotic cell death. Fkh, along with two other TFs, Hkb and Rib, also coordinates organ morphogenesis, transforming two plates of precursor cells on the embryo surface into elongated internalized epithelial tubes. Acquisition of functional specialization, the third key step, is mediated by CrebA and Fkh working in concert with Sage and yet another TF, Sens. CrebA directly upregulates expression of all of the components of the secretory machinery as well as other genes (e.g., Xbp1) necessary for managing the physiological stress that inexorably accompanies high secretory load. Secretory cargo specificity is controlled by Sage and Sens in collaboration with Fkh. Investigations have also uncovered roles for various signaling pathways, e.g., Dpp signaling, EGF signaling, GPCR signaling, and cytoskeletal signaling, and their interactions within the gene regulatory networks that specify, build, and specialize the SG. Collectively, studies of the SG have expanded our knowledge of secretory dynamics, cell polarity, and cytoskeletal mechanics in the context of organ development and function. Notably, the embryonic SG has made the singular contribution as a model system that revealed the core function of CrebA in scaling up secretory capacity, thus, serving as the pioneer system in which the conserved roles of the mammalian Creb3/3L-family orthologues were first discovered.

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.

The human HOXA9 protein uses paralog-specific residues of the homeodomain to interact with TALE-class cofactors

HOX proteins interact with PBX and MEIS cofactors, which belong to the TALE-class of homeodomain (HD)-containing transcription factors. Although the formation of HOX-PBX complexes depends on a unique conserved HOX motif called hexapeptide (HX), the additional presence of MEIS induces a remodeling of the interaction, leading to a global dispensability of the HX motif for trimeric complex formation in the large majority of HOX proteins. In addition, it was shown that the anterior HOXB3 and central HOXA7 and HOXC8 proteins could use different alternative TALE interaction motifs, with or without the HX motif, depending on the DNA-binding site and cell context. Here we dissected the molecular interaction properties of the human posterior HOXA9 protein with its TALE cofactors, PBX1 and MEIS1. Analysis was performed on different DNA-binding sites in vitro and by doing Bimolecular Fluorescence Complementation (BiFC) in different cell lines. Notably, we observed that the HOXA9-TALE interaction relies consistently on the redundant activity of the HX motif and two paralog-specific residues of the HOXA9 HD. Together with previous work, our results show that HOX proteins interact with their generic TALE cofactors through various modalities, ranging from unique and context-independent to versatile and context-dependent TALE binding interfaces.

Dynamic Studies on Intrinsically Disordered Regions of Two Paralogous Transcription Factors Reveal Rigid Segments with Important Biological Functions

Long stretches of intrinsically disordered regions (IDRs) are abundantly present in eukaryotic transcription factors. Although their biological significance is well appreciated, the underlying structural and dynamic mechanisms of their function are still not clear. Using solution NMR spectroscopy, we have studied the structural and dynamic features of two paralogous HOX transcription factors, SCR and DFD, from Drosophila. Both proteins have a conserved DNA-binding homeodomain and a long stretch of functionally important IDR. Using NMR dynamics, we determined flexibility of each residue in these proteins. The flexibility of the residues in the disordered region is not uniform. In both proteins, the IDRs have short stretches of consecutive residues with relatively less flexibility, that is, higher rigidity. We show that one such rigid segment is specifically recognized by another co-transcription factor, thus highlighting the importance of these rigid segments in IDR-mediated protein-protein interactions. Using molecular dynamics simulation, we further show that the rigid segments sample less conformations compared to the rest of the residues in the disordered region. The restrained conformational sampling of these rigid residues should lower the loss in conformational entropy during their interactions with binding partners resulting in sequence specific binding. This work provides experimental evidence of a "rigid-segment" model of IDRs, where functionally important rigid segments are connected by highly flexible linkers. Furthermore, a comparative study of IDRs in paralogous proteins reveals that in spite of low-sequence conservation, the rigid and flexible segments are sequentially maintained to preserve related functions and regulations of these proteins.

The cis-regulatory logic underlying abdominal Hox-mediated repression versus activation of regulatory elements in Drosophila

During development diverse transcription factor inputs are integrated by cis-regulatory modules (CRMs) to yield cell-specific gene expression. Defining how CRMs recruit the appropriate combinations of factors to either activate or repress gene expression remains a challenge. In this study, we compare and contrast the ability of two CRMs within the Drosophila embryo to recruit functional Hox transcription factor complexes. The DCRE CRM recruits Ultrabithorax (Ubx) and Abdominal-A (Abd-A) Hox complexes that include the Extradenticle (Exd) and Homothorax (Hth) transcription factors to repress the Distal-less leg selector gene, whereas the RhoA CRM selectively recruits Abd-A/Exd/Hth complexes to activate rhomboid and stimulate Epidermal Growth Factor secretion in sensory cell precursors. By swapping binding sites between these elements, we found that the RhoA Exd/Hth/Hox site configuration that mediates Abd-A specific activation can convey transcriptional repression by both Ubx and Abd-A when placed into the DCRE. We further show that the orientation and spacing of Hox sites relative to additional binding sites within the RhoA and DCRE is critical to mediate cell- and segment-specific output. These results indicate that the configuration of Exd, Hth, and Hox site within RhoA is neither Abd-A specific nor activation specific. Instead Hox specific output is largely dependent upon the presence of appropriately spaced and oriented binding sites for additional TF inputs. Taken together, these studies provide insight into the cis-regulatory logic used to generate cell-specific outputs via recruiting Hox transcription factor complexes.

The Generic Facet of Hox Protein Function

Hox transcription factors are essential to promote morphological diversification of the animal body. A substantial number of studies have focused on how Hox proteins reach functional specificity, an issue that arises from the fact that these transcription factors control distinct developmental functions despite sharing similar molecular properties. In this review, we highlight that, besides specific functions, for which these transcription factors are renowned, Hox proteins also often have nonspecific functions. We next discuss some emerging principles of these generic functions and how they relate to specific functions and explore our current grasp of the underlying molecular mechanisms.

Combinatorial action of Grainyhead, Extradenticle and Notch in regulating Hox mediated apoptosis in Drosophila larval CNS

Hox mediated neuroblast apoptosis is a prevalent way to pattern larval central nervous system (CNS) by different Hox genes, but the mechanism of this apoptosis is not understood. Our studies with Abdominal-A (Abd-A) mediated larval neuroblast (pNB) apoptosis suggests that AbdA, its cofactor Extradenticle (Exd), a helix-loop-helix transcription factor Grainyhead (Grh), and Notch signaling transcriptionally contribute to expression of RHG family of apoptotic genes. We find that Grh, AbdA, and Exd function together at multiple motifs on the apoptotic enhancer. In vivo mutagenesis of these motifs suggest that they are important for the maintenance of the activity of the enhancer rather than its initiation. We also find that Exd function is independent of its known partner homothorax in this apoptosis. We extend some of our findings to Deformed expressing region of sub-esophageal ganglia where pNBs undergo a similar Hox dependent apoptosis. We propose a mechanism where common players like Exd-Grh-Notch work with different Hox genes through region specific enhancers to pattern respective segments of larval central nervous system.

Specification of the head to tail axis of the developing central nervous system is carried out by Hox genes. Hox mediated programmed cell death of the neural progenitor cells plays an important role in specification of this axis, but the molecular mechanism of this phenomenon is not well understood. We have studied this phenomenon in abdominal and subesophageal regions of larval central nervous system of Drosophila. We find that different Hox genes use a combination of common players (Extradenticle, Grainyhead and Notch) but employ region specific enhancers to cause progenitor cell death in different segments of developing central nervous system.

Non-specificity of transcription factor function in Drosophila melanogaster

A major problem in developmental genetics is how HOX transcription factors, like Proboscipedia (PB) and Ultrabithorax (UBX), regulate distinct programs of gene expression to result in a proboscis versus a haltere, respectively, when the DNA-binding homeodomain (HD) of HOX transcription factors recognizes similar DNA-binding sequences. Indeed, the lack of DNA-binding specificity is a problem for all transcription factors (TFs), as the DNA-binding domains generally recognize small targets of five to six bases in length. Although not the initial intent of the study, I found extensive non-specificity of TF function. Multiple TFs including HOX and HD-containing and non-HD-containing TFs induced both wingless and eyeless phenotypes. The TFs Labial (LAB), Deformed (DFD), Fushi tarazu (FTZ), and Squeeze (SQZ) induced ectopic larval thoracic (T) 1 beard formation in T2 and T3. The TF Doublesex male (DSX^M) rescued the reduced maxillary palp pb phenotype. These examples of non-specificity of TF function across classes of TFs, combined with previous observations, compromise the implicit, initial assumption often made that an intrinsic mechanism of TF specificity is important for function. Interestingly, the functional complementation of the pb phenotype may suggest a larger role for regulation of expression of TFs in restriction of function as opposed to an intrinsic specificity of TF function. These observations have major ramifications for analysis of functional conservation in evolution and development.

Sensitivity of Allelic Divergence to Genomic Position: Lessons from the Drosophila tan Gene

To identify genetic variants underlying changes in phenotypes within and between species, researchers often utilize transgenic animals to compare the function of alleles in different genetic backgrounds. In Drosophila, targeted integration mediated by the ΦC31 integrase allows activity of alternative alleles to be compared at the same genomic location. By using the same insertion site for each transgene, position effects are generally assumed to be controlled for because both alleles are surrounded by the same genomic context. Here, we test this assumption by comparing the activity of tan alleles from two Drosophila species, D. americana and D. novamexicana, at five different genomic locations in D. melanogaster. We found that the relative effects of these alleles varied among insertion sites, with no difference in activity observed between them at two sites. One of these sites simply silenced both transgenes, but the other allowed expression of both alleles that was sufficient to rescue a mutant phenotype yet failed to reveal the functional differences between the two alleles. These results suggest that more than one insertion site should be used when comparing the activity of transgenes because failing to do so could cause functional differences between alleles to go undetected.

A Hox Transcription Factor Collective Binds a Highly Conserved Distal-less cis-Regulatory Module to Generate Robust Transcriptional Outcomes

cis-regulatory modules (CRMs) generate precise expression patterns by integrating numerous transcription factors (TFs). Surprisingly, CRMs that control essential gene patterns can differ greatly in conservation, suggesting distinct constraints on TF binding sites. Here, we show that a highly conserved Distal-less regulatory element (DCRE) that controls gene expression in leg precursor cells recruits multiple Hox, Extradenticle (Exd) and Homothorax (Hth) complexes to mediate dual outputs: thoracic activation and abdominal repression. Using reporter assays, we found that abdominal repression is particularly robust, as neither individual binding site mutations nor a DNA binding deficient Hth protein abolished cooperative DNA binding and in vivo repression. Moreover, a re-engineered DCRE containing a distinct configuration of Hox, Exd, and Hth sites also mediated abdominal Hox repression. However, the re-engineered DCRE failed to perform additional segment-specific functions such as thoracic activation. These findings are consistent with two emerging concepts in gene regulation: First, the abdominal Hox/Exd/Hth factors utilize protein-protein and protein-DNA interactions to form repression complexes on flexible combinations of sites, consistent with the TF collective model of CRM organization. Second, the conserved DCRE mediates multiple cell-type specific outputs, consistent with recent findings that pleiotropic CRMs are associated with conserved TF binding and added evolutionary constraints.

Enhancers are regulatory elements that interact with transcription factor proteins to control cell-specific gene expression during development. Surprisingly, only a subset of enhancers are highly conserved at the sequence level, even though the expression patterns they control are often conserved and essential for proper development. Why some enhancer sequences are highly conserved whereas others are not is not well understood. In this study, we characterize a highly conserved enhancer that regulates gene expression in leg precursor cells. We find that this enhancer has dual regulatory activities that include gene activation in thoracic segments and gene repression in abdominal segments. Surprisingly, we show that the conserved enhancer can tolerate numerous sequence changes yet mediate robust transcription factor binding and abdominal repression. These findings are consistent with abdominal transcription factors binding numerous different configurations of binding sites. So, why is this enhancer highly conserved? We found that overlapping sequences within the enhancer also contribute to thoracic activation, suggesting the enhancer sequences are under added functional constraints. Altogether, our results provide new insights into why some enhancers are highly conserved at the sequence level while others can tolerate sequence changes.

To Be Specific or Not: The Critical Relationship Between Hox And TALE Proteins

Hox proteins are key regulatory transcription factors that act in different tissues of the embryo to provide specific spatial and temporal coordinates to each cell. These patterning functions often depend on the presence of the TALE-homeodomain class cofactors, which form cooperative DNA-binding complexes with all Hox proteins. How this family of cofactors contributes to the highly diverse and specific functions of Hox proteins in vivo remains an important unsolved question. We review here the most recent advances in understanding the molecular mechanisms underlying Hox-TALE function. In particular, we discuss the role of DNA shape, DNA-binding affinity, and protein-protein interaction flexibility in dictating Hox-TALE specificity. We propose several models to explain how these mechanisms are integrated with each other in the context of the many distinct functions that Hox and TALE factors carry out in vivo.

Cell-Autonomous and Non-cell-autonomous Function of Hox Genes Specify Segmental Neuroblast Identity in the Gnathal Region of the Embryonic CNS in Drosophila

During central nervous system (CNS) development neural stem cells (Neuroblasts, NBs) have to acquire an identity appropriate to their location. In thoracic and abdominal segments of Drosophila, the expression pattern of Bithorax-Complex Hox genes is known to specify the segmental identity of NBs prior to their delamination from the neuroectoderm. Compared to the thoracic, ground state segmental units in the head region are derived to different degrees, and the precise mechanism of segmental specification of NBs in this region is still unclear. We identified and characterized a set of serially homologous NB-lineages in the gnathal segments and used one of them (NB6-4 lineage) as a model to investigate the mechanism conferring segment-specific identities to gnathal NBs. We show that NB6-4 is primarily determined by the cell-autonomous function of the Hox gene Deformed (Dfd). Interestingly, however, it also requires a non-cell-autonomous function of labial and Antennapedia that are expressed in adjacent anterior or posterior compartments. We identify the secreted molecule Amalgam (Ama) as a downstream target of the Antennapedia-Complex Hox genes labial, Dfd, Sex combs reduced and Antennapedia. In conjunction with its receptor Neurotactin (Nrt) and the effector kinase Abelson tyrosine kinase (Abl), Ama is necessary in parallel to the cell-autonomous Dfd pathway for the correct specification of the maxillary identity of NB6-4. Both pathways repress CyclinE (CycE) and loss of function of either of these pathways leads to a partial transformation (40%), whereas simultaneous mutation of both pathways leads to a complete transformation (100%) of NB6-4 segmental identity. Finally, we provide genetic evidences, that the Ama-Nrt-Abl-pathway regulates CycE expression by altering the function of the Hippo effector Yorkie in embryonic NBs. The disclosure of a non-cell-autonomous influence of Hox genes on neural stem cells provides new insight into the process of segmental patterning in the developing CNS.

The central nervous system (CNS) needs to be subdivided into functionally specified regions. In the developing CNS of Drosophila, each neural stem cell, called neuroblasts (NB), acquires a unique identity according to its anterior-posterior and dorso-ventral position to generate a specific cell lineage. Along the anterior-posterior body axis, Hox genes of the Bithorax-Complex convey segmental identities to NBs in the trunk segments. In the derived gnathal and brain segments, the mechanisms specifying segmental NB identities are largely unknown. We investigated the role of Hox genes of the Antennapedia-Complex in the gnathal CNS. In addition to cell-autonomous Hox gene function, we unexpectedly uncovered a parallel non-cell-autonomous pathway in mediating segmental specification of embryonic NBs in gnathal segments. Both pathways restrict the expression of the cell cycle gene CyclinE, ensuring the proper specification of a glial cell lineage. Whereas the Hox gene Deformed mediates this cell-autonomously, labial and Antennapedia influence the identity via transcriptional regulation of the secreted molecule Amalgam (and its downstream pathway) in a non-cell-autonomous manner. These findings shed new light on the role of the highly conserved Hox genes during segmental patterning of neural stem cells in the CNS.

Cellular and molecular insights into Hox protein action

Hox genes encode homeodomain transcription factors that control morphogenesis and have established functions in development and evolution. Hox proteins have remained enigmatic with regard to the molecular mechanisms that endow them with specific and diverse functions, and to the cellular functions that they control. Here, we review recent examples of Hox-controlled cellular functions that highlight their versatile and highly context-dependent activity. This provides the setting to discuss how Hox proteins control morphogenesis and organogenesis. We then summarise the molecular modalities underlying Hox protein function, in particular in light of current models of transcription factor function. Finally, we discuss how functional divergence between Hox proteins might be achieved to give rise to the many facets of their action.

A flexible extension of the Drosophila ultrabithorax homeodomain defines a novel Hox/PBC interaction mode

The patterning function of Hox proteins relies on assembling protein complexes with PBC proteins, which often involves a protein motif found in most Hox proteins, the so-called Hexapeptide (HX). Hox/PBC complexes likely gained functional diversity by acquiring additional modes of interaction. Here, we structurally characterize the first HX alternative interaction mode based on the paralogue-specific UbdA motif and further functionally validate structure-based predictions. The UbdA motif folds as a flexible extension of the homeodomain recognition helix and defines Hox/PBC contacts that occur, compared with those mediated by the HX motif, on the opposing side of the DNA double helix. This provides a new molecular facet to Hox/PBC complex assembly and suggests possible mechanisms for the diversification of Hox protein function.

Probing the kinetic landscape of Hox transcription factor-DNA binding in live cells by massively parallel Fluorescence Correlation Spectroscopy

Hox genes encode transcription factors that control the formation of body structures, segment-specifically along the anterior-posterior axis of metazoans. Hox transcription factors bind nuclear DNA pervasively and regulate a plethora of target genes, deploying various molecular mechanisms that depend on the developmental and cellular context. To analyze quantitatively the dynamics of their DNA-binding behavior we have used confocal laser scanning microscopy (CLSM), single-point fluorescence correlation spectroscopy (FCS), fluorescence cross-correlation spectroscopy (FCCS) and bimolecular fluorescence complementation (BiFC). We show that the Hox transcription factor Sex combs reduced (Scr) forms dimers that strongly associate with its specific fork head binding site (fkh250) in live salivary gland cell nuclei. In contrast, dimers of a constitutively inactive, phospho-mimicking variant of Scr show weak, non-specific DNA-binding. Our studies reveal that nuclear dynamics of Scr is complex, exhibiting a changing landscape of interactions that is difficult to characterize by probing one point at a time. Therefore, we also provide mechanistic evidence using massively parallel FCS (mpFCS). We found that Scr dimers are predominantly formed on the DNA and are equally abundant at the chromosomes and an introduced multimeric fkh250 binding-site, indicating different mobilities, presumably reflecting transient binding with different affinities on the DNA. Our proof-of-principle results emphasize the advantages of mpFCS for quantitative characterization of fast dynamic processes in live cells.

Role of Homothorax in region specific regulation of Deformed in embryonic neuroblasts

The expression and regulation of Hox genes in developing central nervous system (CNS) lack important details like specific cell types where Hox genes are expressed and the transcriptional regulatory players involved in these cells. In this study we have investigated the expression and regulation of Drosophila Hox gene Deformed (Dfd) in specific cell types of embryonic CNS. Using Dfd neural autoregulatory enhancer we find that Dfd autoregulates itself in cells of mandibular neuromere. We have also investigated the role of a Hox cofactor Homothorax (Hth) for its role in regulating Dfd expression in CNS. We find that Hth exhibits a region specific role in controlling the expression of Dfd, but has no direct role in mandibular Dfd neural autoregulatory circuit. Our results also suggest that homeodomain of Hth is not required for regulating Dfd expression in embryonic CNS.

A survey of conservation of sea spider and Drosophila Hox protein activities

Hox proteins have well-established functions in development and evolution, controlling the final morphology of bilaterian animals. The common phylogenetic origin of Hox proteins and the associated evolutionary diversification of protein sequences provide a unique framework to explore the relationship between changes in protein sequence and function. In this study, we aimed at questioning how sequence variation within arthropod Hox proteins influences function. This was achieved by exploring the functional impact of sequence conservation/divergence of the Hox genes, labial, Sex comb reduced, Deformed, Ultrabithorax and abdominalA from two distant arthropods, the sea spider and the well-studied Drosophila. Results highlight a correlation between sequence conservation within the homeodomain and the degree of functional conservation, and identify a novel functional domain in the Labial protein.

Hoxb6 can interfere with somitogenesis in the posterior embryo through a mechanism independent of its rib-promoting activity

Formation of the vertebrate axial skeleton requires coordinated Hox gene activity. Hox group 6 genes are involved in the formation of the thoracic area due to their unique rib-promoting properties. We show here that the linker region (LR) connecting the homeodomain and the hexapeptide is essential for Hoxb6 rib-promoting activity. The LR-defective Hoxb6 protein was still able to bind a target enhancer together with Pax3 producing a dominant negative effect, indicating that the LR brings additional regulatory factors to target DNA elements. We also found an unexpected association between Hoxb6 and segmentation in the paraxial mesoderm. In particular, Hoxb6 can disturb somitogenesis and anterior-posterior somite patterning by deregulating Lfng expression. Interestingly, this interaction occurred differently in thoracic and more caudal embryonic areas, indicating functional differences in somitogenesis before and after the trunk to tail transition. Our results suggest the requirement of precisely regulated Hoxb6 expression for proper segmentation at tailbud stages.

Differential pleiotropy and HOX functional organization

Key studies led to the idea that transcription factors are composed of defined modular protein motifs or domains, each with separable, unique function. During evolution, the recombination of these modular domains could give rise to transcription factors with new properties, as has been shown using recombinant molecules. This archetypic, modular view of transcription factor organization is based on the analyses of a few transcription factors such as GAL4, which may represent extreme exemplars rather than an archetype or the norm. Recent work with a set of Homeotic selector (HOX) proteins has revealed differential pleiotropy: the observation that highly-conserved HOX protein motifs and domains make small, additive, tissue specific contributions to HOX activity. Many of these differentially pleiotropic HOX motifs may represent plastic sequence elements called short linear motifs (SLiMs). The coupling of differential pleiotropy with SLiMs, suggests that protein sequence changes in HOX transcription factors may have had a greater impact on morphological diversity during evolution than previously believed. Furthermore, differential pleiotropy may be the genetic consequence of an ensemble nature of HOX transcription factor allostery, where HOX proteins exist as an ensemble of states with the capacity to integrate an extensive array of developmental information. Given a new structural model for HOX functional domain organization, the properties of the archetypic TF may require reassessment.

Intrinsically disordered proteins and multicellular organisms

Intrinsically disordered proteins (IDPs) and IDP regions lack stable tertiary structure yet carry out numerous biological functions, especially those associated with signaling, transcription regulation, DNA condensation, cell division, and cellular differentiation. Both post-translational modifications (PTMs) and alternative splicing (AS) expand the functional repertoire of IDPs. Here we propose that an "IDP-based developmental toolkit," which is comprised of IDP regions, PTMs, especially multiple PTMs, within these IDP regions, and AS events within segments of pre-mRNA that code for these same IDP regions, allows functional diversification and environmental responsiveness for molecules that direct the development of complex metazoans.

Drosophila melanogaster Hox transcription factors access the RNA polymerase II machinery through direct homeodomain binding to a conserved motif of mediator subunit Med19

Hox genes in species across the metazoa encode transcription factors (TFs) containing highly-conserved homeodomains that bind target DNA sequences to regulate batteries of developmental target genes. DNA-bound Hox proteins, together with other TF partners, induce an appropriate transcriptional response by RNA Polymerase II (PolII) and its associated general transcription factors. How the evolutionarily conserved Hox TFs interface with this general machinery to generate finely regulated transcriptional responses remains obscure. One major component of the PolII machinery, the Mediator (MED) transcription complex, is composed of roughly 30 protein subunits organized in modules that bridge the PolII enzyme to DNA-bound TFs. Here, we investigate the physical and functional interplay between Drosophila melanogaster Hox developmental TFs and MED complex proteins. We find that the Med19 subunit directly binds Hox homeodomains, in vitro and in vivo. Loss-of-function Med19 mutations act as dose-sensitive genetic modifiers that synergistically modulate Hox-directed developmental outcomes. Using clonal analysis, we identify a role for Med19 in Hox-dependent target gene activation. We identify a conserved, animal-specific motif that is required for Med19 homeodomain binding, and for activation of a specific Ultrabithorax target. These results provide the first direct molecular link between Hox homeodomain proteins and the general PolII machinery. They support a role for Med19 as a PolII holoenzyme-embedded “co-factor” that acts together with Hox proteins through their homeodomains in regulated developmental transcription.

Mutations of Hox developmental genes in the fruit fly Drosophila melanogaster may provoke spectacular changes in form: transformations of one body part into another, or loss of organs. This attribute identifies them as important developmental genes. Insect and vertebrate Hox proteins contain highly related homeodomain motifs used to bind to regulatory DNA and influence expression of developmental target genes. This occurs at the level of transcription of target gene DNA to messenger RNA by RNA polymerase II and its associated protein machinery (>50 proteins). How Hox homeodomain proteins induce fine-tuned transcription remains an open question. We provide an initial response, finding that Hox proteins also use their homeodomains to bind one machinery protein, Mediator complex subunit 19 (Med19) through a Med19 sequence that is highly conserved in animal phyla. Med19 mutants isolated in this work (the first animal mutants) show that Med19 assists Hox protein functions. Further, they indicate that homeodomain binding to the Med19 motif is required for normal expression of a Hox target gene. Our work provides new clues for understanding how the specific transcriptional inputs of the highly conserved Hox class of transcription factors are integrated at the level of the whole transcription machinery.

A combination of activation and repression by a colinear Hox code controls forelimb-restricted expression of Tbx5 and reveals Hox protein specificity

Tight control over gene expression is essential for precision in embryonic development and acquisition of the regulatory elements responsible is the predominant driver for evolution of new structures. Tbx5 and Tbx4, two genes expressed in forelimb and hindlimb-forming regions respectively, play crucial roles in the initiation of limb outgrowth. Evolution of regulatory elements that activate Tbx5 in rostral LPM was essential for the acquisition of forelimbs in vertebrates. We identified such a regulatory element for Tbx5 and demonstrated Hox genes are essential, direct regulators. While the importance of Hox genes in regulating embryonic development is clear, Hox targets and the ways in which each protein executes its specific function are not known. We reveal how nested Hox expression along the rostro-caudal axis restricts Tbx5 expression to forelimb. We demonstrate that Hoxc9, which is expressed in caudal LPM where Tbx5 is not expressed, can form a repressive complex on the Tbx5 forelimb regulatory element. This repressive capacity is limited to Hox proteins expressed in caudal LPM and carried out by two separate protein domains in Hoxc9. Forelimb-restricted expression of Tbx5 and ultimately forelimb formation is therefore achieved through co-option of two characteristics of Hox genes; their colinear expression along the body axis and the functional specificity of different paralogs. Active complexes can be formed by Hox PG proteins present throughout the rostral-caudal LPM while restriction of Tbx5 expression is achieved by superimposing a dominant repressive (Hoxc9) complex that determines the caudal boundary of Tbx5 expression. Our results reveal the regulatory mechanism that ensures emergence of the forelimbs at the correct position along the body. Acquisition of this regulatory element would have been critical for the evolution of limbs in vertebrates and modulation of the factors we have identified can be molecular drivers of the diversity in limb morphology.

The acquisition of limbs during vertebrate evolution was a very successful innovation that enabled this group of species to diversify and colonise land. It has become clear recently that the primary driver behind the evolution of new structures, such as limbs, is the acquisition of novel regulatory elements that control when and where genes are activated rather than the proteins encoded by the genes themselves acquiring novel functions. We have identified the regulatory element from a gene, Tbx5. Activation of Tbx5 in the forelimb-forming region of the developing embryos is essential for forelimbs to form and disruption of human TBX5 causes limb abnormalities. We show that activation of Tbx5 in a restricted territory is achieved through a combination of activation inputs that are present broadly throughout the embryo flank and dominant, repressive inputs present only in more caudal regions of the flank. The sum of these inputs yields restricted activation in the rostral, forelimb-forming flank. Our results explain how the regulatory switches that were harnessed for the acquisition of limbs during evolution operate and how they can be turned off during the evolution of limblessness in species such as the snake.

Molecular insights into the origin of the Hox-TALE patterning system

Despite tremendous body form diversity in nature, bilaterian animals share common sets of developmental genes that display conserved expression patterns in the embryo. Among them are the Hox genes, which define different identities along the anterior–posterior axis. Hox proteins exert their function by interaction with TALE transcription factors. Hox and TALE members are also present in some but not all non-bilaterian phyla, raising the question of how Hox–TALE interactions evolved to provide positional information. By using proteins from unicellular and multicellular lineages, we showed that these networks emerged from an ancestral generic motif present in Hox and other related protein families. Interestingly, Hox-TALE networks experienced additional and extensive molecular innovations that were likely crucial for differentiating Hox functions along body plans. Together our results highlight how homeobox gene families evolved during eukaryote evolution to eventually constitute a major patterning system in Eumetazoans.

DOI:http://dx.doi.org/10.7554/eLife.01939.001

Any animal with a body that is symmetric about an imaginary line that runs from its head to its tail is known as a bilaterian. Humans and most animals are bilateral, whereas jellyfish and starfish are not. Bilateral symmetry can take many forms—as demonstrated by the differences between flies, frogs and humans—but all bilaterians express many of the same genes during development.

One of these groups of genes is known as the Hox family. The expression of specific Hox genes at specific times instructs cells in the developing embryo to adopt different fates according to their position along the anterior–posterior (head to tail) axis. The patterning function of Hox genes relies on the presence of two additional cofactors that belong to the so-called TALE family. Although both Hox and TALE proteins were present early on during animal evolution, it is unclear how and when the interactions between them first began to generate symmetrical body plans.

Now, Hudry et al. have provided insights into the origin of the Hox-TALE network by analysing the expression and molecular properties of Hox and TALE proteins from various multicellular and unicellular organisms. These experiments revealed that Hox and TALE proteins of the sea anemone Nematostella, which belongs to a group of animals called cnidarians that have radial rather than bilateral symmetry, interact with one another in a similar manner to the interactions seen in bilaterians.

Hudry et al. then showed that two Nematostella Hox genes were able to substitute for their bilaterian equivalents in fruit flies, and that a Nematostella TALE gene was able to take over neuronal functions of its equivalent in Xenopus frogs. This striking conservation of function between species suggests that Hox and TALE genes were already working together in the common ancestor of all bilaterian and cnidarian animals.

By contrast, TALE members from a unicellular amoeba were unable to interact with Hox proteins, suggesting that Hox–TALE interactions first emerged in multicellular animals. In addition to increasing our knowledge of highly conserved Hox signalling, these data provide insight into the molecular mechanisms that gave rise to the symmetrical body plan that has been adopted, and adapted, by the majority of animals since.

DOI:http://dx.doi.org/10.7554/eLife.01939.002

Onychophoran Hox genes and the evolution of arthropod Hox gene expression

Introduction

Onychophora is a relatively small phylum within Ecdysozoa, and is considered to be the sister group to Arthropoda. Compared to the arthropods, that have radiated into countless divergent forms, the onychophoran body plan is overall comparably simple and does not display much in-phylum variation. An important component of arthropod morphological diversity consists of variation of tagmosis, i.e. the grouping of segments into functional units (tagmata), and this in turn is correlated with differences in expression patterns of the Hox genes. How these genes are expressed in the simpler onychophorans, the subject of this paper, would therefore be of interest in understanding their subsequent evolution in the arthropods, especially if an argument can be made for the onychophoran system broadly reflecting the ancestral state in the arthropods.

Results

The sequences and embryonic expression patterns of the complete set of ten Hox genes of an onychophoran (Euperipatoides kanangrensis) are described for the first time. We find that they are all expressed in characteristic patterns that suggest a function as classical Hox genes. The onychophoran Hox genes obey spatial colinearity, and with the exception of Ultrabithorax (Ubx), they all have different and distinct anterior expression borders. Notably, Ubx transcripts form a posterior to anterior gradient in the onychophoran trunk. Expression of all onychophoran Hox genes extends continuously from their anterior border to the rear end of the embryo.

Conclusions

The spatial expression pattern of the onychophoran Hox genes may contribute to a combinatorial Hox code that is involved in giving each segment its identity. This patterning of segments in the uniform trunk, however, apparently predates the evolution of distinct segmental differences in external morphology seen in arthropods. The gradient-like expression of Ubx may give posterior segments their specific identity, even though they otherwise express the same set of Hox genes. We suggest that the confined domains of Hox gene expression seen in arthropods evolved from an ancestral onychophoran-like Hox gene pattern. Reconstruction of the ancestral arthropod Hox pattern and comparison with the patterns in the different arthropod classes reveals phylogenetic support for Mandibulata and Tetraconata, but not Myriochelata and Atelocerata.

Measuring Hox-DNA binding by electrophoretic mobility shift analysis

Understanding gene regulation by Hox transcription factors requires understanding the forces that underlie DNA binding by these proteins. Electrophoretic mobility shift analysis (EMSA) not only allows measurement of protein affinity and cooperativity but also permits visualization of differently migrating protein-DNA complexes, including complexes with different compositions or complexes with identical compositions yet assembled in different geometries. Furthermore, protein activity can be measured, allowing correction of binding constants for the percentage of protein that is properly folded and capable of binding DNA. Protocols for measuring protein activity and the equilibrium DNA-binding dissociation constant (K d) are provided. This versatile assay system can be adjusted based on specific needs to measure other parameters, including the kinetic association and dissociation constants (k a and k d) and the formation of heterologous protein-protein interactions.

Binding polymorphism in the DNA bound state of the Pdx1 homeodomain

The subtle effects of DNA-protein recognition are illustrated in the homeodomain fold. This is one of several small DNA binding motifs that, in spite of limited DNA binding specificity, adopts crucial, specific roles when incorporated in a transcription factor. The homeodomain is composed of a 3-helix domain and a mobile N-terminal arm. Helix 3 (the recognition helix) interacts with the DNA bases through the major groove, while the N-terminal arm becomes ordered upon binding a specific sequence through the minor groove. Although many structural studies have characterized the DNA binding properties of homeodomains, the factors behind the binding specificity are still difficult to elucidate. A crystal structure of the Pdx1 homeodomain bound to DNA (PDB 2H1K) obtained previously in our lab shows two complexes with differences in the conformation of the N-terminal arm, major groove contacts, and backbone contacts, raising new questions about the DNA recognition process by homeodomains. Here, we carry out fully atomistic Molecular Dynamics simulations both in crystal and aqueous environments in order to elucidate the nature of the difference in binding contacts. The crystal simulations reproduce the X-ray experimental structures well. In the absence of crystal packing constraints, the differences between the two complexes increase during the solution simulations. Thus, the conformational differences are not an artifact of crystal packing. In solution, the homeodomain with a disordered N-terminal arm repositions to a partially specific orientation. Both the crystal and aqueous simulations support the existence of different stable binding conformers identified in the original crystallographic data with different degrees of specificity. We propose that protein-protein and protein-DNA interactions favor a subset of the possible conformations. This flexibility in DNA binding may facilitate multiple functions for the same transcription factor.

Hox transcription factors influence motoneuron identity through the integrated actions of both homeodomain and non-homeodomain regions

BACKGROUND

Hox transcription factors play a critical role in the specification of motoneuron subtypes within the spinal cord. Our previous work showed that two orthologous members of this family, Hoxd10 and Hoxd11, exert opposing effects on motoneuron development in the lumbosacral (LS) spinal cord of the embryonic chick: Hoxd10 promotes the development of lateral motoneuron subtypes that project to dorsal limb muscles, while Hoxd11 represses the development of lateral subtypes in favor of medial subtypes that innervate ventral limb muscles and axial muscles. The striking degree of homology between the DNA-binding homeodomains of Hoxd10 and Hoxd11 suggested that non-homeodomain regions mediate their divergent effects. In the present study, we investigate the relative contributions of homeodomain and non-homeodomain regions of Hoxd10 and Hoxd11 to motoneuron specification.

RESULTS

Using in ovo electroporation to express chimeric and mutant constructs in LS motoneurons, we find that both the homeodomain and non-homeodomain regions of Hoxd10 are necessary to specify lateral motoneurons. In contrast, non-homeodomain regions of Hoxd11 are sufficient to repress lateral motoneuron fates in favor of medial fates.

CONCLUSIONS

Together, our data demonstrate that even closely related Hox orthologues rely on distinct combinations of homeodomain-dependent and -independent mechanisms to specify motoneuron identity.

Dual role for Hox genes and Hox co-factors in conferring leg motoneuron survival and identity in Drosophila

Adult Drosophila walk using six multi-jointed legs, each controlled by ∼50 leg motoneurons (MNs). Although MNs have stereotyped morphologies, little is known about how they are specified. Here, we describe the function of Hox genes and homothorax (hth), which encodes a Hox co-factor, in Drosophila leg MN development. Removing either Hox or Hth function from a single neuroblast (NB) lineage results in MN apoptosis. A single Hox gene, Antennapedia (Antp), is primarily responsible for MN survival in all three thoracic segments. When cell death is blocked, partially penetrant axon branching errors are observed in Hox mutant MNs. When single MNs are mutant, errors in both dendritic and axon arborizations are observed. Our data also suggest that Antp levels in post-mitotic MNs are important for specifying their identities. Thus, in addition to being essential for survival, Hox and hth are required to specify accurate MN morphologies in a level-dependent manner.

Distinct molecular strategies for Hox-mediated limb suppression in Drosophila: from cooperativity to dispensability/antagonism in TALE partnership

The emergence following gene duplication of a large repertoire of Hox paralogue proteins underlies the importance taken by Hox proteins in controlling animal body plans in development and evolution. Sequence divergence of paralogous proteins accounts for functional specialization, promoting axial morphological diversification in bilaterian animals. Yet functionally specialized paralogous Hox proteins also continue performing ancient common functions. In this study, we investigate how highly divergent Hox proteins perform an identical function. This was achieved by comparing in Drosophila the mode of limb suppression by the central (Ultrabithorax and AbdominalA) and posterior class (AbdominalB) Hox proteins. Results highlight that Hox-mediated limb suppression relies on distinct modes of DNA binding and a distinct use of TALE cofactors. Control of common functions by divergent Hox proteins, at least in the case studied, relies on evolving novel molecular properties. Thus, changes in protein sequences not only provide the driving force for functional specialization of Hox paralogue proteins, but also provide means to perform common ancient functions in distinct ways.

Genetic and functional modularity of Hox activities in the specification of limb-innervating motor neurons

A critical step in the assembly of the neural circuits that control tetrapod locomotion is the specification of the lateral motor column (LMC), a diverse motor neuron population targeting limb musculature. Hox6 paralog group genes have been implicated as key determinants of LMC fate at forelimb levels of the spinal cord, through their ability to promote expression of the LMC-restricted genes Foxp1 and Raldh2 and to suppress thoracic fates through exclusion of Hoxc9. The specific roles and mechanisms of Hox6 gene function in LMC neurons, however, are not known. We show that Hox6 genes are critical for diverse facets of LMC identity and define motifs required for their in vivo specificities. Although Hox6 genes are necessary for generating the appropriate number of LMC neurons, they are not absolutely required for the induction of forelimb LMC molecular determinants. In the absence of Hox6 activity, LMC identity appears to be preserved through a diverse array of Hox5–Hox8 paralogs, which are sufficient to reprogram thoracic motor neurons to an LMC fate. In contrast to the apparently permissive Hox inputs to early LMC gene programs, individual Hox genes, such as Hoxc6, have specific roles in promoting motor neuron pool diversity within the LMC. Dissection of motifs required for Hox in vivo specificities reveals that either cross-repressive interactions or cooperativity with Pbx cofactors are sufficient to induce LMC identity, with the N-terminus capable of promoting columnar, but not pool, identity when transferred to a heterologous homeodomain. These results indicate that Hox proteins orchestrate diverse aspects of cell fate specification through both the convergent regulation of gene programs regulated by many paralogs and also more restricted actions encoded through specificity determinants in the N-terminus.

Coordinated motor behaviors—as complex as playing a musical instrument or as simple as walking—rely on the ability of motor neurons within the spinal cord to navigate towards and establish specific connections with muscles in the limbs. The establishment of connections between motor neurons and limb muscles is mediated through the actions of genes encoding Hox proteins, a large family of transcription factors conserved amongst all metazoans. However, the specific requirements for Hox genes in motor neuron specification and patterns of muscle connectivity are poorly understood. We have found that members of the Hox6 gene paralog group (Hoxa6, Hoxc6, and Hoxb6) contribute to diverse aspects of motor neuron subtype differentiation. Hox6 gene activity is required during two critical phases of motor neuron development: first as motor axons select a trajectory toward the forelimb and second as they choose specific muscles to innervate. At the molecular level, these two functions are encoded by distinct peptide domains within Hox proteins. This work indicates that Hox proteins execute their critical functions in motor neurons through intrinsic modules that confer distinct specificities and that these activities are central in the genetic network required for motor neuron differentiation.

The transcriptional specificity of NF-κB dimers is coded within the κB DNA response elements

Nuclear factor κB (NF-κB) regulates gene expression by binding to specific DNA elements, known collectively as κB sites, that are contained within the promoters/enhancers of target genes. We found that the identity of the central base pair (bp) of κB sites profoundly affects the transcriptional activity of NF-κB dimers. RelA dimers prefer an A/T bp at this position for optimal transcriptional activation (A/T-centric) and discriminate against G/C-centric κB sites. The p52 homodimer, in contrast, activates transcription from G/C-centric κB sites in complex with Bcl3 but represses transcription from the A/T-centric sites. The p52:Bcl3 complex binds to these two classes of κB sites in distinct modes, permitting the recruitment of coactivator, corepressor, or both coactivator and corepressor complexes in promoters that contain G/C-, A/T-, or both G/C- and A/T-centric sites. Therefore, through sensing of bp differences within κB sites, NF-κB dimers modulate biological programs by activating, repressing, and altering the expression of effector genes.

Disentangling the many layers of eukaryotic transcriptional regulation

Regulation of gene expression in eukaryotes is an extremely complex process. In this review, we break down several critical steps, emphasizing new data and techniques that have expanded current gene regulatory models. We begin at the level of DNA sequence where cis-regulatory modules (CRMs) provide important regulatory information in the form of transcription factor (TF) binding sites. In this respect, CRMs function as instructional platforms for the assembly of gene regulatory complexes. We discuss multiple mechanisms controlling complex assembly, including cooperative DNA binding, combinatorial codes, and CRM architecture. The second section of this review places CRM assembly in the context of nucleosomes and condensed chromatin. We discuss how DNA accessibility and histone modifications contribute to TF function. Lastly, new advances in chromosomal mapping techniques have provided increased understanding of intra- and interchromosomal interactions. We discuss how these topological maps influence gene regulatory models.

Regulatory role for a conserved motif adjacent to the homeodomain of Hox10 proteins

Development of the vertebrate axial skeleton requires the concerted activity of several Hox genes. Among them, Hox genes belonging to the paralog group 10 are essential for the formation of the lumbar region of the vertebral column, owing to their capacity to block rib formation. In this work, we explored the basis for the rib-repressing activity of Hox10 proteins. Because genetic experiments in mice demonstrated that Hox10 proteins are strongly redundant in this function, we first searched for common motifs among the group members. We identified the presence of two small sequences flanking the homeodomain that are phylogenetically conserved among Hox10 proteins and that seem to be specific for this group. We show here that one of these motifs is required but not sufficient for the rib-repressing activity of Hox10 proteins. This motif includes two potential phosphorylation sites, which are essential for protein activity as their mutation to alanines resulted in a total loss of rib-repressing properties. Our data indicates that this motif has a significant regulatory function, modulating interactions with more N-terminal parts of the Hox protein, eventually triggering the rib-repressing program. In addition, this motif might also regulate protein activity by alteration of the protein's DNA-binding affinity through changes in the phosphorylation state of two conserved tyrosine residues within the homeodomain.

The cis-regulatory code of Hox function in Drosophila

Precise gene expression is a fundamental aspect of organismal function and depends on the combinatorial interplay of transcription factors (TFs) with cis-regulatory DNA elements. While much is known about TF function in general, our understanding of their cell type-specific activities is still poor. To address how widely expressed transcriptional regulators modulate downstream gene activity with high cellular specificity, we have identified binding regions for the Hox TF Deformed (Dfd) in the Drosophila genome. Our analysis of architectural features within Hox cis-regulatory response elements (HREs) shows that HRE structure is essential for cell type-specific gene expression. We also find that Dfd and Ultrabithorax (Ubx), another Hox TF specifying different morphological traits, interact with non-overlapping regions in vivo, despite their similar DNA binding preferences. While Dfd and Ubx HREs exhibit comparable design principles, their motif compositions and motif-pair associations are distinct, explaining the highly selective interaction of these Hox proteins with the regulatory environment. Thus, our results uncover the regulatory code imprinted in Hox enhancers and elucidate the mechanisms underlying functional specificity of TFs in vivo.

Hox proteins display a common and ancestral ability to diversify their interaction mode with the PBC class cofactors

Hox protein function during development and evolution relies on conserved multiple interaction modes with cofactors of the PBC and Meis families.

Hox transcription factors control a number of developmental processes with the help of the PBC class proteins. In vitro analyses have established that the formation of Hox/PBC complexes relies on a short conserved Hox protein motif called the hexapeptide (HX). This paradigm is at the basis of the vast majority of experimental approaches dedicated to the study of Hox protein function. Here we questioned the unique and general use of the HX for PBC recruitment by using the Bimolecular Fluorescence Complementation (BiFC) assay. This method allows analyzing Hox-PBC interactions in vivo and at a genome-wide scale. We found that the HX is dispensable for PBC recruitment in the majority of investigated Drosophila and mouse Hox proteins. We showed that HX-independent interaction modes are uncovered by the presence of Meis class cofactors, a property which was also observed with Hox proteins of the cnidarian sea anemone Nematostella vectensis. Finally, we revealed that paralog-specific motifs convey major PBC-recruiting functions in Drosophila Hox proteins. Altogether, our results highlight that flexibility in Hox-PBC interactions is an ancestral and evolutionary conserved character, which has strong implications for the understanding of Hox protein functions during normal development and pathologic processes.

Hox proteins are key transcriptional regulators of animal development, famously helping to determine identity along the anterior-posterior body axis. Although their evolution and developmental roles are well established, the molecular mechanisms underlying their specific functions remain poorly characterized. The current dominant view is that interaction with different members of the PBC family of transcription factors confers specific DNA-binding properties on different Hox proteins. However, this idea conflicts with in vitro evidence that a short “hexapeptide” (HX) motif shared by most Hox proteins is solely responsible for generic PBC recruitment. Here we have used the BiFC (bimolecular fluorescence complementation) method to address the global importance of the HX motif for Hox-PBC interactions in living cells and living animals including fruit flies and chick embryos. We observe that most interactions between Hox and PBC proteins do not depend on HX, and that alternative protein motifs are widely used for PBC recruitment in vivo. We also show that DNA binding by a second family of cofactors, the Meis proteins, unmasks these alternative interaction modes and that this property is conserved not only across Bilateria, but also in the basal animal phylum Cnidaria. Taken together, our results demonstrate that Hox-PBC partnership relies on multiple interaction modes, which can be influenced by additional transcriptional partners. We propose that this ancestral feature has been essential for ensuring Hox functional plasticity during development and evolution.

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.

Competition for cofactor-dependent DNA binding underlies Hox phenotypic suppression

Hox transcription factors exhibit an evolutionarily conserved functional hierarchy, termed phenotypic suppression, in which the activity of posterior Hox proteins dominates over more anterior Hox proteins. Using directly regulated Hox targeted reporter genes in Drosophila, we show that posterior Hox proteins suppress the activities of anterior ones by competing for cofactor-dependent DNA binding. Furthermore, we map a motif in the posterior Hox protein Abdominal-A (AbdA) that is required for phenotypic suppression and facilitates cooperative DNA binding with the Hox cofactor Extradenticle (Exd). Together, these results suggest that Hox-specific motifs endow posterior Hox proteins with the ability to dominate over more anterior ones via a cofactor-dependent DNA-binding mechanism.

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.

Evolution of a derived protein-protein interaction between HoxA11 and Foxo1a in mammals caused by changes in intramolecular regulation

Current models of developmental evolution suggest changes in gene regulation underlie the evolution of morphology. Despite the fact that protein complexes regulate gene expression, the evolution of regulatory protein complexes is rarely studied. Here, we investigate the evolution of a protein-protein interaction (PPI) between Homeobox A11 (HoxA11) and Forkhead box 01A (Foxo1a). Using extant and "resurrected" ancestral proteins, we show that the physical interaction between HoxA11 and Foxo1a originated in the mammalian stem lineage. Functional divergence tests and coimmunoprecipitation with heterologous protein pairs indicate that the evolution of interaction was attributable to changes in HoxA11, and deletion studies demonstrate that the interaction interface is located in the homeodomain region of HoxA11. However, there are no changes in amino acid sequence in the homeodomain region during this time period, indicating that the origin of the derived PPI was attributable to changes outside the binding interface. We infer that the amino acid substitutions in HoxA11 altered Foxo1a's access to the conserved binding interface at the HoxA11 homeodomain. We also found an expansion in the number of paired Hox/Fox binding sites in the genomes of mammalian lineage species suggesting the complex has a biological function. Our data indicate that the physical interaction between HoxA11 and Foxo1a evolved through noninterface changes that facilitate the PPI, which prevents inappropriate interactions, rather than through the evolution of a novel binding interface. We speculate that evolutionary changes of intramolecular regulation have limited pleiotropic effects compared with changes to interaction domains themselves.