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

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.

Expanding the Iroquois genes repertoire: a non-transcriptional function in cell cycle progression

Drosophila Iroquois (Iro) proteins are components of the TALE homeodomain family of transcriptional regulators. They play key roles in territorial specification and pattern formation. A recent study has disclosed a novel developmental function of the Iro proteins. In the eye and wing imaginal discs, they can regulate the size of the territories that they specify. They do so by cell-autonomously controlling cell cycle progression. Indeed, Iro proteins down-regulate the activity of the CyclinE/Cdk2 complex by a transcription-independent mechanism. This novel function is executed mainly through 2 evolutionarily conserved domains of the Iro proteins: the Cyclin Binding Domain and the IRO-box, which mediate their binding to CyclinE-containing protein complexes. Here we discuss the functional implications of the control of the cell cycle by Iro proteins for development and oncogenesis.

Roles of cofactors and chromatin accessibility in Hox protein target specificity

Background

The regulation of specific target genes by transcription factors is central to our understanding of gene network control in developmental and physiological processes yet how target specificity is achieved is still poorly understood. This is well illustrated by the Hox family of transcription factors as their limited in vitro DNA-binding specificity contrasts with their clear in vivo functional specificity.

Results

We generated genome-wide binding profiles for three Hox proteins, Ubx, Abd-A and Abd-B, following transient expression in Drosophila Kc167 cells, revealing clear target specificity and a striking influence of chromatin accessibility. In the absence of the TALE class homeodomain cofactors Exd and Hth, Ubx and Abd-A bind at a very similar set of target sites in accessible chromatin, whereas Abd-B binds at an additional specific set of targets. Provision of Hox cofactors Exd and Hth considerably modifies the Ubx genome-wide binding profile enabling Ubx to bind at an additional novel set of targets. Both the Abd-B specific targets and the cofactor-dependent Ubx targets are in chromatin that is relatively DNase1 inaccessible prior to the expression of Hox proteins/Hox cofactors.

Conclusions

Our experiments demonstrate a strong role for chromatin accessibility in Hox protein binding and suggest that Hox protein competition with nucleosomes has a major role in Hox protein target specificity in vivo.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0049-x) contains supplementary material, which is available to authorized users.

Multi-reporter selection for the design of active and more specific zinc-finger nucleases for genome editing

Engineered nucleases have transformed biological research and offer great therapeutic potential by enabling the straightforward modification of desired genomic sequences. While many nuclease platforms have proven functional, all can produce unanticipated off-target lesions and have difficulty discriminating between homologous sequences, limiting their therapeutic application. Here we describe a multi-reporter selection system that allows the screening of large protein libraries to uncover variants able to discriminate between sequences with substantial homology. We have used this system to identify zinc-finger nucleases that exhibit high cleavage activity (up to 60% indels) at their targets within the CCR5 and HBB genes and strong discrimination against homologous sequences within CCR2 and HBD. An unbiased screen for off-target lesions using a novel set of CCR5-targeting nucleases confirms negligible CCR2 activity and demonstrates minimal off-target activity genome wide. This system offers a straightforward approach to generate nucleases that discriminate between similar targets and provide exceptional genome-wide specificity.

Zinc finger nucleases have an established role in genome editing. Here, the authors report a strategy for identifying zinc finger nucleases that discriminate between desired targets and provide genome-wide specificity.

The Soft Touch: Low-Affinity Transcription Factor Binding Sites in Development and Evolution

Transcription factor proteins regulate gene expression by binding to specific DNA regions. Most studies of transcription factor binding sites have focused on the highest affinity sites for each factor. There is abundant evidence, however, that binding sites with a range of affinities, including very low affinities, are critical to gene regulation. Here, we present the theoretical and experimental evidence for the importance of low-affinity sites in gene regulation and development. We also discuss the implications of the widespread use of low-affinity sites in eukaryotic genomes for robustness, precision, specificity, and evolution of gene regulation.

Gene Regulatory Mechanisms Underlying the Spatial and Temporal Regulation of Target-Dependent Gene Expression in Drosophila Neurons

Neuronal differentiation often requires target-derived signals from the cells they innervate. These signals typically activate neural subtype-specific genes, but the gene regulatory mechanisms remain largely unknown. Highly restricted expression of the FMRFa neuropeptide in Drosophila Tv4 neurons requires target-derived BMP signaling and a transcription factor code that includes Apterous. Using integrase transgenesis of enhancer reporters, we functionally dissected the Tv4-enhancer of FMRFa within its native cellular context. We identified two essential but discrete cis-elements, a BMP-response element (BMP-RE) that binds BMP-activated pMad, and a homeodomain-response element (HD-RE) that binds Apterous. These cis-elements have low activity and must be combined for Tv4-enhancer activity. Such combinatorial activity is often a mechanism for restricting expression to the intersection of cis-element spatiotemporal activities. However, concatemers of the HD-RE and BMP-RE cis-elements were found to independently generate the same spatiotemporal expression as the Tv4-enhancer. Thus, the Tv4-enhancer atypically combines two low-activity cis-elements that confer the same output from distinct inputs. The activation of target-dependent genes is assumed to 'wait' for target contact. We tested this directly, and unexpectedly found that premature BMP activity could not induce early FMRFa expression; also, we show that the BMP-insensitive HD-RE cis-element is activated at the time of target contact. This led us to uncover a role for the nuclear receptor, seven up (svp), as a repressor of FMRFa induction prior to target contact. Svp is normally downregulated immediately prior to target contact, and we found that maintaining Svp expression prevents cis-element activation, whereas reducing svp gene dosage prematurely activates cis-element activity. We conclude that the target-dependent FMRFa gene is repressed prior to target contact, and that target-derived BMP signaling directly activates FMRFa gene expression through an atypical gene regulatory mechanism.

Affinity regression predicts the recognition code of nucleic acid-binding proteins

Predicting the affinity profiles of nucleic acid-binding proteins directly from the protein sequence is a major unsolved problem. We present a statistical approach for learning the recognition code of a family of transcription factors (TFs) or RNA-binding proteins (RBPs) from high-throughput binding assays. Our method, called affinity regression, trains on protein binding microarray (PBM) or RNA compete experiments to learn an interaction model between proteins and nucleic acids, using only protein domain and probe sequences as inputs. By training on mouse homeodomain PBM profiles, our model correctly identifies residues that confer DNA-binding specificity and accurately predicts binding motifs for an independent set of divergent homeodomains. Similarly, learning from RNA compete profiles for diverse RBPs, our model can predict the binding affinities of held-out proteins and identify key RNA-binding residues. More broadly, we envision applying our method to model and predict biological interactions in any setting where there is a high-throughput ‘affinity’ readout.

Fine-tuning the onset of myogenesis by homeobox proteins that interact with the Myf5 limb enhancer

Skeletal myogenesis in vertebrates is initiated at different sites of skeletal muscle formation during development, by activation of specific control elements of the myogenic regulatory genes. In the mouse embryo, Myf5 is the first myogenic determination gene to be expressed and its spatiotemporal regulation requires multiple enhancer sequences, extending over 120 kb upstream of the Mrf4-Myf5 locus. An enhancer, located at −57/−58 kb from Myf5, is responsible for its activation in myogenic cells derived from the hypaxial domain of the somite, that will form limb muscles. Pax3 and Six1/4 transcription factors are essential activators of this enhancer, acting on a 145-bp core element. Myogenic progenitor cells that will form the future muscle masses of the limbs express the factors necessary for Myf5 activation when they delaminate from the hypaxial dermomyotome and migrate into the forelimb bud, however they do not activate Myf5 and the myogenic programme until they have populated the prospective muscle masses. We show that Msx1 and Meox2 homeodomain-containing transcription factors bind in vitro and in vivo to specific sites in the 145-bp element, and are implicated in fine-tuning activation of Myf5 in the forelimb. Msx1, when bound between Pax and Six sites, prevents the binding of these key activators, thus inhibiting transcription of Myf5 and consequent premature myogenic differentiation. Meox2 is required for Myf5 activation at the onset of myogenesis via direct binding to other homeodomain sites in this sequence. Thus, these homeodomain factors, acting in addition to Pax3 and Six1/4, fine-tune the entry of progenitor cells into myogenesis at early stages of forelimb development.

Summary: Homeodomain factors Msx1 and Meox2, acting in addition to Pax3 and Six1/4, fine-tune the entry of progenitor cells into myogenesis at early stages of forelimb development via modulation of limb enhancer gene Myf5.

Transcriptional Dynamics Driving MAMP-Triggered Immunity and Pathogen Effector-Mediated Immunosuppression in Arabidopsis Leaves Following Infection with Pseudomonas syringae pv tomato DC3000

Transcriptional reprogramming is integral to effective plant defense. Pathogen effectors act transcriptionally and posttranscriptionally to suppress defense responses. A major challenge to understanding disease and defense responses is discriminating between transcriptional reprogramming associated with microbial-associated molecular pattern (MAMP)-triggered immunity (MTI) and that orchestrated by effectors. A high-resolution time course of genome-wide expression changes following challenge with Pseudomonas syringae pv tomato DC3000 and the nonpathogenic mutant strain DC3000hrpA- allowed us to establish causal links between the activities of pathogen effectors and suppression of MTI and infer with high confidence a range of processes specifically targeted by effectors. Analysis of this information-rich data set with a range of computational tools provided insights into the earliest transcriptional events triggered by effector delivery, regulatory mechanisms recruited, and biological processes targeted. We show that the majority of genes contributing to disease or defense are induced within 6 h postinfection, significantly before pathogen multiplication. Suppression of chloroplast-associated genes is a rapid MAMP-triggered defense response, and suppression of genes involved in chromatin assembly and induction of ubiquitin-related genes coincide with pathogen-induced abscisic acid accumulation. Specific combinations of promoter motifs are engaged in fine-tuning the MTI response and active transcriptional suppression at specific promoter configurations by P. syringae.

Homeodomain proteins: an update

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.

The HOX-Apoptosis Regulatory Interplay in Development and Disease

Apoptosis is a cellular suicide program, which is on the one hand used to remove superfluous cells thereby promoting tissue or organ morphogenesis. On the other hand, the programmed killing of cells is also critical when potentially harmful cells emerge in a developing or adult organism thereby endangering survival. Due to its critical role apoptosis is tightly controlled, however so far, its regulation on the transcriptional level is less studied and understood. Hox genes, a highly conserved gene family encoding homeodomain transcription factors, have crucial roles in development. One of their prominent functions is to shape animal body plans by eliciting different developmental programs along the anterior-posterior axis. To this end, Hox proteins transcriptionally regulate numerous processes in a coordinated manner, including cell-type specification, differentiation, motility, proliferation as well as apoptosis. In this review, we will focus on how Hox proteins control organismal morphology and function by regulating the apoptotic machinery. We will first focus on well-established paradigms of Hox-apoptosis interactions and summarize how Hox transcription factors control morphological outputs and differentially shape tissues along the anterior-posterior axis by fine-tuning apoptosis in a healthy organism. We will then discuss the consequences when this interaction is disturbed and will conclude with some ideas and concepts emerging from these studies.

The achaete-scute complex in Diptera: patterns of noncoding sequence evolution

The achaete‐scute complex (AS‐C) has been a useful paradigm for the study of pattern formation and its evolution. achaete‐scute genes have duplicated and evolved distinct expression patterns during the evolution of cyclorraphous Diptera. Are the expression patterns in different species driven by conserved regulatory elements? If so, when did such regulatory elements arise? Here, we have sequenced most of the AS‐C of the fly Calliphora vicina (including the genes achaete, scute and lethal of scute) to compare noncoding sequences with known cis‐regulatory sequences in Drosophila. The organization of the complex is conserved with respect to Drosophila species. There are numerous small stretches of conserved noncoding sequence that, in spite of high sequence turnover, display binding sites for known transcription factors. Synteny of the blocks of conserved noncoding sequences is maintained suggesting not only conservation of the position of regulatory elements but also an origin prior to the divergence between these two species. We propose that some of these enhancers originated by duplication with their target genes.

The Homeodomain Iroquois Proteins Control Cell Cycle Progression and Regulate the Size of Developmental Fields

During development, proper differentiation and final organ size rely on the control of territorial specification and cell proliferation. Although many regulators of these processes have been identified, how both are coordinated remains largely unknown. The homeodomain Iroquois/Irx proteins play a key, evolutionarily conserved, role in territorial specification. Here we show that in the imaginal discs, reduced function of Iroquois genes promotes cell proliferation by accelerating the G1 to S transition. Conversely, their increased expression causes cell-cycle arrest, down-regulating the activity of the Cyclin E/Cdk2 complex. We demonstrate that physical interaction of the Iroquois protein Caupolican with Cyclin E-containing protein complexes, through its IRO box and Cyclin-binding domains, underlies its activity in cell-cycle control. Thus, Drosophila Iroquois proteins are able to regulate cell-autonomously the growth of the territories they specify. Moreover, our results provide a molecular mechanism for a role of Iroquois/Irx genes as tumour suppressors.

The correct development of body organs, with their characteristic size and shape, requires the coordination of cell division and cell differentiation. Here we show that the Iroquois proteins (Irx in vertebrates) slow down cell division in the Drosophila imaginal discs, in addition to their well-known role in cell fate and territorial specification. In humans, inactivating mutations at the Irx genes are associated to several types of cancer, thus allowing their classification as tumour suppressor genes. We have observed that Drosophila Iroquois genes similarly behave as tumour suppressor genes. Iroquois proteins belong to a family of homeodomain-containing transcriptional regulators. However, our results indicate that they control cell division by a transcription independent mechanism based on their physical interaction with Cyclin E containing complexes, a key player in cell-cycle progression. We have identified two evolutionary conserved domains of Iroquois proteins, different from the homeodomain, involved in that interaction. This new function of Iroquois proteins places them in a key position to coordinate growth and differentiation during normal development. Our results further suggest a molecular mechanism for their role in tumour suppression. Future studies of Irx genes should help to determine if a similar mechanism could operate to help cancer progression when Irx activity is compromised.

The Evolutionarily Conserved LIM Homeodomain Protein LIM-4/LHX6 Specifies the Terminal Identity of a Cholinergic and Peptidergic C. elegans Sensory/Inter/Motor Neuron-Type

The expression of specific transcription factors determines the differentiated features of postmitotic neurons. However, the mechanism by which specific molecules determine neuronal cell fate and the extent to which the functions of transcription factors are conserved in evolution are not fully understood. In C. elegans, the cholinergic and peptidergic SMB sensory/inter/motor neurons innervate muscle quadrants in the head and control the amplitude of sinusoidal movement. Here we show that the LIM homeobox protein LIM-4 determines neuronal characteristics of the SMB neurons. In lim-4 mutant animals, expression of terminal differentiation genes, such as the cholinergic gene battery and the flp-12 neuropeptide gene, is completely abolished and thus the function of the SMB neurons is compromised. LIM-4 activity promotes SMB identity by directly regulating the expression of the SMB marker genes via a distinct cis-regulatory motif. Two human LIM-4 orthologs, LHX6 and LHX8, functionally substitute for LIM-4 in C. elegans. Furthermore, C. elegans LIM-4 or human LHX6 can induce cholinergic and peptidergic characteristics in the human neuronal cell lines. Our results indicate that the evolutionarily conserved LIM-4/LHX6 homeodomain proteins function in generation of precise neuronal subtypes.

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.

A widespread role of the motif environment in transcription factor binding across diverse protein families

Transcriptional regulation requires the binding of transcription factors (TFs) to short sequence-specific DNA motifs, usually located at the gene regulatory regions. Interestingly, based on a vast amount of data accumulated from genomic assays, it has been shown that only a small fraction of all potential binding sites containing the consensus motif of a given TF actually bind the protein. Recent in vitro binding assays, which exclude the effects of the cellular environment, also demonstrate selective TF binding. An intriguing conjecture is that the surroundings of cognate binding sites have unique characteristics that distinguish them from other sequences containing a similar motif that are not bound by the TF. To test this hypothesis, we conducted a comprehensive analysis of the sequence and DNA shape features surrounding the core-binding sites of 239 and 56 TFs extracted from in vitro HT-SELEX binding assays and in vivo ChIP-seq data, respectively. Comparing the nucleotide content of the regions around the TF-bound sites to the counterpart unbound regions containing the same consensus motifs revealed significant differences that extend far beyond the core-binding site. Specifically, the environment of the bound motifs demonstrated unique sequence compositions, DNA shape features, and overall high similarity to the core-binding motif. Notably, the regions around the binding sites of TFs that belong to the same TF families exhibited similar features, with high agreement between the in vitro and in vivo data sets. We propose that these unique features assist in guiding TFs to their cognate binding sites.

Protein-DNA binding in high-resolution

Recent advances in experimental and computational methodologies are enabling ultra-high resolution genome-wide profiles of protein-DNA binding events. For example, the ChIP-exo protocol precisely characterizes protein-DNA cross-linking patterns by combining chromatin immunoprecipitation (ChIP) with 5' → 3' exonuclease digestion. Similarly, deeply sequenced chromatin accessibility assays (e.g. DNase-seq and ATAC-seq) enable the detection of protected footprints at protein-DNA binding sites. With these techniques and others, we have the potential to characterize the individual nucleotides that interact with transcription factors, nucleosomes, RNA polymerases and other regulatory proteins in a particular cellular context. In this review, we explain the experimental assays and computational analysis methods that enable high-resolution profiling of protein-DNA binding events. We discuss the challenges and opportunities associated with such approaches.

Role of HOXA9 in leukemia: dysregulation, cofactors and essential targets

HOXA9 is a homeodomain-containing transcription factor that has an important role in hematopoietic stem cell expansion and is commonly deregulated in acute leukemias. A variety of upstream genetic alterations in acute myeloid leukemia lead to overexpression of HOXA9, which is a strong predictor of poor prognosis. In many cases, HOXA9 has been shown to be necessary for maintaining leukemic transformation; however, the molecular mechanisms through which it promotes leukemogenesis remain elusive. Recent work has established that HOXA9 regulates downstream gene expression through binding at promoter distal enhancers along with a subset of cell-specific cofactor and collaborator proteins. Increasing efforts are being made to identify both the critical cofactors and target genes required for maintaining transformation in HOXA9-overexpressing leukemias. With continued advances in understanding HOXA9-mediated transformation, there is a wealth of opportunity for developing novel therapeutics that would be applicable for greater than 50% of AML with overexpression of HOXA9.

A new hereditary congenital facial palsy case supports arg5 in HOX-DNA binding domain as possible hot spot for mutations

Moebius syndrome (MBS) is a rare congenital disorder characterized by rhombencephalic mal development, mainly presenting with facial palsy with limited gaze abduction. Most cases are sporadic, possibly caused by a combination of environmental and genetic factors; however, no proven specific associations have been yet established. Hereditary congenital facial palsy (HCFP) is an autosomal dominant congenital dysinnervation syndrome, recognizable by the isolated dysfunction of the seventh cranial nerve. Mutant mice for Hoxb1 were reported to present with facial weakness, resembling MBS. Recently a homozygous mutation altering arg5 residue of HOXB1 homeodomain into cys5 was identified in two families with HCFP. We screened 95 sporadic patients diagnosed as MBS or HCFP for mutations in HOXB1. A novel homozygous alteration was identified in one HCFP case, affecting the same residue, resulting to his5. In silico protein analysis predicted stronger HOXB1-DNA binding properties for his5 than cys5 that resulted to milder phenotype. It should be noted that, inclusive of the previous report, only two mutations revealed in HOXB1 associated with HCFP involved the same amino acid arg5 in HOXB1 residing in HOXB1-DNA-PBX1 ternary complex.

Differential Nucleosome Occupancies across Oct4-Sox2 Binding Sites in Murine Embryonic Stem Cells

The binding sequence for any transcription factor can be found millions of times within a genome, yet only a small fraction of these sequences encode functional transcription factor binding sites. One of the reasons for this dichotomy is that many other factors, such as nucleosomes, compete for binding. To study how the competition between nucleosomes and transcription factors helps determine a functional transcription factor site from a predicted transcription factor site, we compared experimentally-generated in vitro nucleosome occupancy with in vivo nucleosome occupancy and transcription factor binding in murine embryonic stem cells. Using a solution hybridization enrichment technique, we generated a high-resolution nucleosome map from targeted regions of the genome containing predicted sites and functional sites of Oct4/Sox2 regulation. We found that at Pax6 and Nes, which are bivalently poised in stem cells, functional Oct4 and Sox2 sites show high amounts of in vivo nucleosome displacement compared to in vitro. Oct4 and Sox2, which are active, show no significant displacement of in vivo nucleosomes at functional sites, similar to nonfunctional Oct4/Sox2 binding. This study highlights a complex interplay between Oct4 and Sox2 transcription factors and nucleosomes among different target genes, which may result in distinct patterns of stem cell gene regulation.

Toward a new twist in Hox and TALE DNA-binding specificity

Hox proteins gain specificity by interacting with TALE-class cofactors. In a recent issue of Cell and in this issue of Developmental Cell, Crocker et al. (2015) and Amin et al. (2015), respectively, demonstrate that non-canonical Hox/TALE binding sequences play a major role in the regionalized regulation of target gene expression in vivo.

Inhibitory activities of short linear motifs underlie Hox interactome specificity in vivo

Hox proteins are well-established developmental regulators that coordinate cell fate and morphogenesis throughout embryogenesis. In contrast, our knowledge of their specific molecular modes of action is limited to the interaction with few cofactors. Here, we show that Hox proteins are able to interact with a wide range of transcription factors in the live Drosophila embryo. In this context, specificity relies on a versatile usage of conserved short linear motifs (SLiMs), which, surprisingly, often restrains the interaction potential of Hox proteins. This novel buffering activity of SLiMs was observed in different tissues and found in Hox proteins from cnidarian to mouse species. Although these interactions remain to be analysed in the context of endogenous Hox regulatory activities, our observations challenge the traditional role assigned to SLiMs and provide an alternative concept to explain how Hox interactome specificity could be achieved during the embryonic development.

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

In all animals, it is important that cells are correctly organised into tissues and organs. This organisation starts in the embryo, and cells are instructed to perform different roles depending on their position within the body.

A family of proteins called the Hox proteins coordinates the organisation of the cells in the animal embryo by binding to and controlling the expression of specific genes. To properly control their target genes, Hox proteins need to interact with other proteins called transcription factors that can also bind to the genes. However, only a few of these transcription factors have been identified so far, and it is not clear how Hox proteins are able to interact with them.

Here, Baëza, Viala, Heim et al. identified several more transcription factors that can bind to the Hox proteins in fruit fly embryos. The experiments show that Hox proteins are able to bind to many transcription factors that are very different from each other. Baëza, Viala, Heim et al. also show that two short sections within the Hox proteins known as short linear motifs are important for controlling these interactions. A fly Hox protein that was missing these motifs was able to interact with new transcription factors. This inhibitory role was found in Hox proteins from mice and sea anemones, suggesting that these motifs may play the same role in all animals.

Baëza, Viala, Heim et al.'s findings challenge the traditional view of the role of the short linear motifs in interactions between proteins. Also, the findings provide an alternative explanation for how the Hox proteins are only able to interact with particular transcription factors in animal embryos. The next step will be to find out whether the inhibitory role of short linear motifs could more generally apply to many other protein families.

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

The evolutionary origination and diversification of a dimorphic gene regulatory network through parallel innovations in cis and trans

The origination and diversification of morphological characteristics represents a key problem in understanding the evolution of development. Morphological traits result from gene regulatory networks (GRNs) that form a web of transcription factors, which regulate multiple cis-regulatory element (CRE) sequences to control the coordinated expression of differentiation genes. The formation and modification of GRNs must ultimately be understood at the level of individual regulatory linkages (i.e., transcription factor binding sites within CREs) that constitute the network. Here, we investigate how elements within a network originated and diversified to generate a broad range of abdominal pigmentation phenotypes among Sophophora fruit flies. Our data indicates that the coordinated expression of two melanin synthesis enzymes, Yellow and Tan, recently evolved through novel CRE activities that respond to the spatial patterning inputs of Hox proteins and the sex-specific input of Bric-à-brac transcription factors. Once established, it seems that these newly evolved activities were repeatedly modified by evolutionary changes in the network’s trans-regulators to generate large-scale changes in pigment pattern. By elucidating how yellow and tan are connected to the web of abdominal trans-regulators, we discovered that the yellow and tan abdominal CREs are composed of distinct regulatory inputs that exhibit contrasting responses to the same Hox proteins and Hox cofactors. These results provide an example in which CRE origination underlies a recently evolved novel trait, and highlights how coordinated expression patterns can evolve in parallel through the generation of unique regulatory linkages.

The genomic content of regulatory genes such as transcription factors is surprisingly conserved between diverse animal species, raising the paradox of how new traits emerge, and are subsequently modified and lost. In this study we make a connection between the developmental basis for the formation of a fruit fly trait and the evolutionary basis for that trait’s origin, diversification, and loss. We show how the origin of a novel pigmentation trait is associated with the evolution of two regulatory sequences that control the co-expression of two key pigmentation genes. These sequences interact in unique ways with evolutionarily conserved Hox transcription factors to drive gene co-expression. Once these unique connections evolved, the alteration of this trait appears to have proceeded through changes to regulatory genes rather than regulatory sequences of the pigmentation genes. Thus, our findings support a scenario where regulatory sequence evolution provided new functions to old transcription factors, how co-expression can emerge from different utilizations of the same transcription factors, and that trait diversity was surprisingly shaped by changes in some manner to the deeply conserved regulatory genes.

Structure-based discovery of NANOG variant with enhanced properties to promote self-renewal and reprogramming of pluripotent stem cells

NANOG (from Irish mythology Tír na nÓg) transcription factor plays a central role in maintaining pluripotency, cooperating with OCT4 (also known as POU5F1 or OCT3/4), SOX2, and other pluripotency factors. Although the physiological roles of the NANOG protein have been extensively explored, biochemical and biophysical properties in relation to its structural analysis are poorly understood. Here we determined the crystal structure of the human NANOG homeodomain (hNANOG HD) bound to an OCT4 promoter DNA, which revealed amino acid residues involved in DNA recognition that are likely to be functionally important. We generated a series of hNANOG HD alanine substitution mutants based on the protein-DNA interaction and evolutionary conservation and determined their biological activities. Some mutant proteins were less stable, resulting in loss or decreased affinity for DNA binding. Overexpression of the orthologous mouse NANOG (mNANOG) mutants failed to maintain self-renewal of mouse embryonic stem cells without leukemia inhibitory factor. These results suggest that these residues are critical for NANOG transcriptional activity. Interestingly, one mutant, hNANOG L122A, conversely enhanced protein stability and DNA-binding affinity. The mNANOG L122A, when overexpressed in mouse embryonic stem cells, maintained their expression of self-renewal markers even when retinoic acid was added to forcibly drive differentiation. When overexpressed in epiblast stem cells or human induced pluripotent stem cells, the L122A mutants enhanced reprogramming into ground-state pluripotency. These findings demonstrate that structural and biophysical information on key transcriptional factors provides insights into the manipulation of stem cell behaviors and a framework for rational protein engineering.

Hoxa2 selectively enhances Meis binding to change a branchial arch ground state

Hox transcription factors (TFs) are essential for vertebrate development, but how these evolutionary conserved proteins function in vivo remains unclear. Because Hox proteins have notoriously low binding specificity, they are believed to bind with cofactors, mainly homeodomain TFs Pbx and Meis, to select their specific targets. We mapped binding of Meis, Pbx, and Hoxa2 in the branchial arches, a series of segments in the developing vertebrate head. Meis occupancy is largely similar in Hox-positive and -negative arches. Hoxa2, which specifies second arch (IIBA) identity, recognizes a subset of Meis prebound sites that contain Hox motifs. Importantly, at these sites Meis binding is strongly increased. This enhanced Meis binding coincides with active enhancers, which are linked to genes highly expressed in the IIBA and regulated by Hoxa2. These findings show that Hoxa2 operates as a tissue-specific cofactor, enhancing Meis binding to specific sites that provide the IIBA with its anatomical identity.

•

Meis provides a ground state that is common to all the branchial arches

•

Hoxa2 recognizes Meis prebound sites in the second arch that contain Hox motifs

•

Hoxa2 enhances Meis binding, which coincides with active enhancers, at these sites

•

Hoxa2 modulates the ground-state binding of Meis to instruct second arch identity

Low affinity binding site clusters confer hox specificity and regulatory robustness

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.

Understanding DNA-binding specificity by bacteria hybrid selection

Understanding how sequence-specific protein-DNA interactions direct cellular function is of great interest to the research community. High-throughput methods have been developed to determine DNA-binding specificities; one such technique, the bacterial one-hybrid (B1H) system, confers advantages including ease of use, sensitivity and throughput. In this review, we describe the evolution of the B1H system as a tool capable of screening large DNA libraries to investigate protein-DNA interactions of interest. We discuss how DNA-binding specificities produced by the B1H system have been used to predict regulatory targets. Additionally, we examine how this approach has been applied to characterize two common DNA-binding domain families-homeodomains and Cys2His2 zinc fingers-both in organism-wide studies and with synthetic approaches. In the case of the former, the B1H system has produced large catalogs of protein specificity and nuanced information about previously recovered DNA targets, thereby improving our understanding of these proteins' functions in vivo and increasing our capacity to predict similar interactions in other species. In the latter, synthetic screens of the same DNA-binding domains have further refined our models of specificity, through analyzing comprehensive libraries to uncover all proteins able to bind a complete set of targets, and, for instance, exploring how context-in the form of domain position within the parent protein-may affect specificity. Finally, we recognize the limitations of the B1H system and discuss its potential for use in the production of designer proteins and in studies of protein-protein interactions.

Structure-based modeling of protein: DNA specificity

Protein:DNA interactions are essential to a range of processes that maintain and express the information encoded in the genome. Structural modeling is an approach that aims to understand these interactions at the physicochemical level. It has been proposed that structural modeling can lead to deeper understanding of the mechanisms of protein:DNA interactions, and that progress in this field can not only help to rationalize the observed specificities of DNA-binding proteins but also to allow researchers to engineer novel DNA site specificities. In this review we discuss recent developments in the structural description of protein:DNA interactions and specificity, as well as the challenges facing the field in the future.

Determination and inference of eukaryotic transcription factor sequence specificity

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Divergent mechanisms regulate conserved cardiopharyngeal development and gene expression in distantly related ascidians

Ascidians present a striking dichotomy between conserved phenotypes and divergent genomes: embryonic cell lineages and gene expression patterns are conserved between distantly related species. Much research has focused on Ciona or Halocynthia spp. but development in other ascidians remains poorly characterized. In this study, we surveyed the multipotent myogenic B7.5 lineage in Molgula spp. Comparisons to the homologous lineage in Ciona revealed identical cell division and fate specification events that result in segregation of larval, cardiac, and pharyngeal muscle progenitors. Moreover, the expression patterns of key regulators are conserved, but cross-species transgenic assays uncovered incompatibility, or ‘unintelligibility’, of orthologous cis-regulatory sequences between Molgula and Ciona. These sequences drive identical expression patterns that are not recapitulated in cross-species assays. We show that this unintelligibility is likely due to changes in both cis- and trans-acting elements, hinting at widespread and frequent turnover of regulatory mechanisms underlying otherwise conserved aspects of ascidian embryogenesis.

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

When two species have features that look similar, this may be because the features arise by the same processes during development. Other features may look similar yet develop by different mechanisms. ‘Developmental system drift’ refers to the process where a physical feature remains unaltered during evolution, but the underlying pathway that controls its development is changed. However, to date, there have been only a few experimental studies that support this idea.

Ascidians—also commonly known as sea squirts—are vase-like marine creatures, which start off as tadpole-like larvae that swim around until they find a place to settle down and attach themselves. Once attached, the sea squirts lose the ability to swim and start feeding, typically by filtering material out of the seawater. Sea squirts and their close relatives are the invertebrates (animals without backbones) that are most closely related to all vertebrates (animals with backbones), including humans. Furthermore, although different species of sea squirt have almost identical embryos, their genomes are very different.

Stolfi et al. have now studied whether developmental system drift may have occurred during the evolution of ascidians, by analyzing different species of sea squirt named Molgula and Ciona. Stolfi et al. compared the genomes of Molgula and Ciona and studied the expression of genes in the cells that give rise to the heart and the muscles of the head. As an embryo develops, specific genes are switched on or off, and these patterns of gene activation were broadly identical in the two species of sea squirt examined.

Enhancers are sequences of DNA that control when and how a gene is switched on. Given the similarities between the development of heart and head muscle cells in the different sea squirts, Stolfi et al. looked to see if the mechanisms of gene expression, and therefore the enhancers, were also conserved. Unexpectedly, this was not the case. When enhancers from Molgula were introduced into Ciona (and vice versa), these sequences were unable to switch on gene expression—thus enhancers from one sea squirt species could not function in the other.

Stolfi et al. conclude that the developmental systems may have drifted considerably during evolution of the sea squirts, in spite of their nearly identical embryos. This reinforces the view that different paths can lead to the formation of similar physical features.

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

High-resolution specificity from DNA sequencing highlights alternative modes of Lac repressor binding

Knowing the specificity of transcription factors is critical to understanding regulatory networks in cells. The lac repressor-operator system has been studied for many years, but not with high-throughput methods capable of determining specificity comprehensively. Details of its binding interaction and its selection of an asymmetric binding site have been controversial. We employed a new method to accurately determine relative binding affinities to thousands of sequences simultaneously, requiring only sequencing of bound and unbound fractions. An analysis of 2560 different DNA sequence variants, including both base changes and variations in operator length, provides a detailed view of lac repressor sequence specificity. We find that the protein can bind with nearly equal affinities to operators of three different lengths, but the sequence preference changes depending on the length, demonstrating alternative modes of interaction between the protein and DNA. The wild-type operator has an odd length, causing the two monomers to bind in alternative modes, making the asymmetric operator the preferred binding site. We tested two other members of the LacI/GalR protein family and find that neither can bind with high affinity to sites with alternative lengths or shows evidence of alternative binding modes. A further comparison with known and predicted motifs suggests that the lac repressor may be unique in this ability and that this may contribute to its selection.

The role of Six1 in the genesis of muscle cell and skeletal muscle development

The sine oculis homeobox 1 (Six1) gene encodes an evolutionarily conserved transcription factor. In the past two decades, much research has indicated that Six1 is a powerful regulator participating in skeletal muscle development. In this review, we summarized the discovery and structural characteristics of Six1 gene, and discussed the functional roles and molecular mechanisms of Six1 in myogenesis and in the formation of skeletal muscle fibers. Finally, we proposed areas of future interest for understanding Six1 gene function.

Absence of a simple code: how transcription factors read the genome

Transcription factors (TFs) influence cell fate by interpreting the regulatory DNA within a genome. TFs recognize DNA in a specific manner; the mechanisms underlying this specificity have been identified for many TFs based on 3D structures of protein-DNA complexes. More recently, structural views have been complemented with data from high-throughput in vitro and in vivo explorations of the DNA-binding preferences of many TFs. Together, these approaches have greatly expanded our understanding of TF-DNA interactions. However, the mechanisms by which TFs select in vivo binding sites and alter gene expression remain unclear. Recent work has highlighted the many variables that influence TF-DNA binding, while demonstrating that a biophysical understanding of these many factors will be central to understanding TF function.

Occupancy of tissue-specific cis-regulatory modules by Otx2 and TLE/Groucho for embryonic head specification

Head specification by the head-selector gene, orthodenticle (otx), is highly conserved among bilaterian lineages. However, the molecular mechanisms by which Otx and other transcription factors (TFs) interact with the genome to direct head formation are largely unknown. Here we employ ChIP-seq and RNA-seq approaches in Xenopus tropicalis gastrulae and find that occupancy of the corepressor, TLE/Groucho, is a better indicator of tissue-specific cis-regulatory modules (CRMs) than the coactivator p300, during early embryonic stages. On the basis of TLE binding and comprehensive CRM profiling, we define two distinct types of Otx2- and TLE-occupied CRMs. Using these devices, Otx2 and other head organizer TFs (for example, Lim1/Lhx1 (activator) or Goosecoid (repressor)) are able to upregulate or downregulate a large battery of target genes in the head organizer. An underlying principle is that Otx marks target genes for head specification to be regulated positively or negatively by partner TFs through specific types of CRMs.

Common transcriptional mechanisms for visual photoreceptor cell differentiation among Pancrustaceans

A hallmark of visual rhabdomeric photoreceptors is the expression of a rhabdomeric opsin and uniquely associated phototransduction molecules, which are incorporated into a specialized expanded apical membrane, the rhabdomere. Given the extensive utilization of rhabdomeric photoreceptors in the eyes of protostomes, here we address whether a common transcriptional mechanism exists for the differentiation of rhabdomeric photoreceptors. In Drosophila, the transcription factors Pph13 and Orthodenticle (Otd) direct both aspects of differentiation: rhabdomeric opsin transcription and rhabdomere morphogenesis. We demonstrate that the orthologs of both proteins are expressed in the visual systems of the distantly related arthropod species Tribolium castaneum and Daphnia magna and that their functional roles are similar in these species. In particular, we establish that the Pph13 homologs have the ability to bind a subset of Rhodopsin core sequence I sites and that these sites are present in key phototransduction genes of both Tribolium and Daphnia. Furthermore, Pph13 and Otd orthologs are capable of executing deeply conserved functions of photoreceptor differentiation as evidenced by the ability to rescue their respective Drosophila mutant phenotypes. Pph13 homologs are equivalent in their ability to direct both rhabdomere morphogenesis and opsin expression within Drosophila, whereas Otd paralogs demonstrate differential abilities to regulate photoreceptor differentiation. Finally, loss-of-function analyses in Tribolium confirm the conserved requirement of Pph13 and Otd in regulating both rhabdomeric opsin transcription and rhabdomere morphogenesis. Taken together, our data identify components of a regulatory framework for rhabdomeric photoreceptor differentiation in Pancrustaceans, providing a foundation for defining ancestral regulatory modules of rhabdomeric photoreceptor differentiation.

Visual systems are populated by one of two fundamental types of photoreceptors, ciliary and rhabdomeric. Each photoreceptor type is defined by the opsin molecule expressed and the final morphological form adapted to house the phototransduction machinery. Here we address whether a common transcriptional mechanisms exists for the differentiation of rhabdomeric photoreceptors. We demonstrate that orthologs of two Drosophila (fruit fly) transcription factors, Pph13 and Orthodenticle, are expressed in photoreceptors of Pancrustaceans, Tribolium (red flour beetle) and Daphnia (water flea), and are capable of executing conserved functions of rhabdomeric photoreceptor differentiation. In particular, Tribolium and Daphnia orthologs are capable of substituting and rescuing the photoreceptor differentiation defects observed in their corresponding Drosophila mutants. Furthermore, loss of function analysis in Tribolium of both Pph13 and orthodenticle genes demonstrate they regulate opsin transcription and morphogenesis of the photoreceptor apical membrane. Our data illuminate a framework for rhabdomeric photoreceptor differentiation and provide the foundation for defining the ancestral regulatory modules for rhabdomeric differentiation and potential modifications that underlie the functional diversity observed in rhabdomeric photoreceptors.

The Cryptosporidium parvum ApiAP2 gene family: insights into the evolution of apicomplexan AP2 regulatory systems

We provide the first comprehensive analysis of any transcription factor family in Cryptosporidium, a basal-branching apicomplexan that is the second leading cause of infant diarrhea globally. AP2 domain-containing proteins have evolved to be the major regulatory family in the phylum to the exclusion of canonical regulators. We show that apicomplexan and perkinsid AP2 domains cluster distinctly from other chromalveolate AP2s. Protein-binding specificity assays of C. parvum AP2 domains combined with motif conservation upstream of co-regulated gene clusters allowed the construction of putative AP2 regulons across the in vitro life cycle. Orthologous Apicomplexan AP2 (ApiAP2) expression has been rearranged relative to the malaria parasite P. falciparum, suggesting ApiAP2 network rewiring during evolution. C. hominis orthologs of putative C. parvum ApiAP2 proteins and target genes show greater than average variation. C. parvum AP2 domains display reduced binding diversity relative to P. falciparum, with multiple domains binding the 5'-TGCAT-3', 5'-CACACA-3' and G-box motifs (5'-G[T/C]GGGG-3'). Many overrepresented motifs in C. parvum upstream regions are not AP2 binding motifs. We propose that C. parvum is less reliant on ApiAP2 regulators in part because it utilizes E2F/DP1 transcription factors. C. parvum may provide clues to the ancestral state of apicomplexan transcriptional regulation, pre-AP2 domination.

The Cavβ1a subunit regulates gene expression and suppresses myogenin in muscle progenitor cells

Ca_vβ_1a acts as a voltage-gated calcium channel-independent regulator of gene expression in muscle progenitor cells and is required for their normal expansion during myogenic development.

Voltage-gated calcium channel (Ca_v) β subunits are auxiliary subunits to Ca_vs. Recent reports show Ca_vβ subunits may enter the nucleus and suggest a role in transcriptional regulation, but the physiological relevance of this localization remains unclear. We sought to define the nuclear function of Ca_vβ in muscle progenitor cells (MPCs). We found that Ca_vβ_1a is expressed in proliferating MPCs, before expression of the calcium conducting subunit Ca_v1.1, and enters the nucleus. Loss of Ca_vβ_1a expression impaired MPC expansion in vitro and in vivo and caused widespread changes in global gene expression, including up-regulation of myogenin. Additionally, we found that Ca_vβ_1a localizes to the promoter region of a number of genes, preferentially at noncanonical (NC) E-box sites. Ca_vβ_1a binds to a region of the Myog promoter containing an NC E-box, suggesting a mechanism for inhibition of myogenin gene expression. This work indicates that Ca_vβ_1a acts as a Ca_v-independent regulator of gene expression in MPCs, and is required for their normal expansion during myogenic development.

A DNA-binding-site landscape and regulatory network analysis for NAC transcription factors in Arabidopsis thaliana

Target gene identification for transcription factors is a prerequisite for the systems wide understanding of organismal behaviour. NAM-ATAF1/2-CUC2 (NAC) transcription factors are amongst the largest transcription factor families in plants, yet limited data exist from unbiased approaches to resolve the DNA-binding preferences of individual members. Here, we present a TF-target gene identification workflow based on the integration of novel protein binding microarray data with gene expression and multi-species promoter sequence conservation to identify the DNA-binding specificities and the gene regulatory networks of 12 NAC transcription factors. Our data offer specific single-base resolution fingerprints for most TFs studied and indicate that NAC DNA-binding specificities might be predicted from their DNA-binding domain's sequence. The developed methodology, including the application of complementary functional genomics filters, makes it possible to translate, for each TF, protein binding microarray data into a set of high-quality target genes. With this approach, we confirm NAC target genes reported from independent in vivo analyses. We emphasize that candidate target gene sets together with the workflow associated with functional modules offer a strong resource to unravel the regulatory potential of NAC genes and that this workflow could be used to study other families of transcription factors.

Ecdysone response elements in the distal promoter of the Bombyx Broad-Complex gene, BmBR-C

The Bombyx mori silkworm's homologue of the Broad-Complex gene (BmBR-C) is transcribed from two promoters: a distal promoter (Pdist) and a proximal promoter (Pprox). As determined by a luciferase assay, the transcriptional activity of Pdist, but not Pprox, was activated by ecdysone. Further analyses using reporters driven by sequential deletion Pdist mutants indicated that two regions, ecdysone responsive element (EcRE)-D and EcRE-P, -4950 bp and -3480 bp upstream from the distal transcription start site, respectively, were important in the responsiveness of Pdist to 20-hydroxyecdysone (20E); however, no significant sequence similarities were found between the canonical EcRE and the EcRE-D or EcRE-P regions. Electrophoretic mobility shift assays showed that both the EcRE-D and -P sequences specifically bound to Bombyx protein(s). Sequence analyses and competition assays suggested that the protein(s) bound to EcRE-P might include components other than the ecdysone receptor (EcR), suggesting that BmBR-C transcription was indirectly activated by ecdysone through the EcRE-P. Remarkably, protein binding to the mid-region of the EcRE-D, EcRE-Db, was competitively inhibited by an oligonucleotide containing the Drosophila hsp27 EcRE sequence. Furthermore, an anti-EcR antibody interfered with the formation of the protein-EcRE-Db complex. These results indicated that a functional Bombyx ecdysone receptor binds to EcRE-D and activates the expression of BmBR-C.

Genomics meets proteomics: identifying the culprits in disease

Genome-wide association studies (GWAS) revealed genomic risk loci that potentially have an impact on disease and phenotypic traits. This extensive resource holds great promise in providing novel directions for personalized medicine, including disease risk prediction, prevention and targeted medication. One of the major challenges that researchers face on the path between the initial identification of an association and precision treatment of patients is the comprehension of the biological mechanisms that underlie these associations. Currently, the focus to solve these questions lies on the integrative analysis of system-wide data on global genome variation, gene expression, transcription factor binding, epigenetic profiles and chromatin conformation. The generation of this data mainly relies on next-generation sequencing. However, due to multiple recent developments, mass spectrometry-based proteomics now offers additional, by the GWAS field so far hardly recognized possibilities for the identification of functional genome variants and, in particular, for the identification and characterization of (differentially) bound protein complexes as well as physiological target genes. In this review, we introduce these proteomics advances and suggest how they might be integrated in post-GWAS workflows. We argue that the combination of highly complementary techniques is powerful and can provide an unbiased, detailed picture of GWAS loci and their mechanistic involvement in disease.

SoxNeuro orchestrates central nervous system specification and differentiation in Drosophila and is only partially redundant with Dichaete

Background

Sox proteins encompass an evolutionarily conserved family of transcription factors with critical roles in animal development and stem cell biology. In common with vertebrates, the Drosophila group B proteins SoxNeuro and Dichaete are involved in central nervous system development, where they play both similar and unique roles in gene regulation. Sox genes show extensive functional redundancy across metazoans, but the molecular basis underpinning functional compensation mechanisms at the genomic level are currently unknown.

Results

Using a combination of genome-wide binding analysis and gene expression profiling, we show that SoxNeuro directs embryonic neural development from the early specification of neuroblasts through to the terminal differentiation of neurons and glia. To address the issue of functional redundancy and compensation at a genomic level, we compare SoxNeuro and Dichaete binding, identifying common and independent binding events in wild-type conditions, as well as instances of compensation and loss of binding in mutant backgrounds.

Conclusions

We find that early aspects of group B Sox functions in the central nervous system, such as stem cell maintenance and dorsoventral patterning, are highly conserved. However, in contrast to vertebrates, we find that Drosophila group B1 proteins also play prominent roles during later aspects of neural morphogenesis. Our analysis of the functional relationship between SoxNeuro and Dichaete uncovers evidence for redundant and independent functions for each protein, along with unexpected examples of compensation and interdependency, thus providing new insights into the general issue of transcription factor functional redundancy.

The homeobox transcription factor Irxl1 negatively regulates MyoD expression and myoblast differentiation

Irxl1/Mkx (Iroquois homeobox-like 1/Mohawk) encodes a member of the TALE subfamily of homeodomain proteins. It is expressed in multiple mesoderm-derived tissues and has recently been shown to regulate tendon differentiation during mouse embryonic development. Previously we showed that knockdown of Irxl1 in zebrafish caused a deficit in neural crest cells which consequently resulted in deformation of craniofacial muscles and arch cartilages. Here, we further demonstrate that loss of Irxl1 function results in deformed somites with disordered muscle fibers and myotendinous junctions. Because expression of myoD is increased in the somites of Irxl1 knockdown morphants, we test whether Irxl1 negatively regulates myoD expression. When stable C2C12 myoblasts overexpressing Irxl1/Mkx were induced to differentiate, myotube formation was inhibited and protein levels of myoD and myosin heavy chain were decreased accordingly. A series of deletion constructs of myoD promoter fragments were tested by luciferase reporter assays, which identified a promoter fragment that is necessary and sufficient for Irxl1-mediated repression. Direct interaction of Irxl1 and myoD promoter was subsequently elucidated by yeast one-hybrid assays, electrophoretic mobility shift assays and chromatin immunoprecipitation analysis. Furthermore, mouse Mkx also binds to and represses myoD promoter. These results indicate that Irxl1/Mkx can repress myoD expression through direct binding to its promoter and may thus play a negative regulatory role in muscle differentiation.

Impacts of the ubiquitous factor Zelda on Bicoid-dependent DNA binding and transcription in Drosophila

The Drosophila transcription factor Bicoid (Bcd) binds thousands of genomic sites during early embryogenesis, but it is unclear how many of these binding events are functionally important. Here, Small and colleagues test the role of the maternal factor Zelda (Zld) in Bcd-mediated binding and transcription. Embryos lacking Zld show enhanced Bcd binding to a subset of genomic locations, causing early activation of target genes normally silent until later stages. This study demonstrates a critical role for Zld in controlling Bcd binding and target gene activation in the early embryo.

In vivo cross-linking studies suggest that the Drosophila transcription factor Bicoid (Bcd) binds to several thousand sites during early embryogenesis, but it is not clear how many of these binding events are functionally important. In contrast, reporter gene studies have identified >60 Bcd-dependent enhancers, all of which contain clusters of the consensus binding sequence TAATCC. These studies also identified clusters of TAATCC motifs (inactive fragments) that failed to drive Bcd-dependent activation. In general, active fragments showed higher levels of Bcd binding in vivo and were enriched in predicted binding sites for the ubiquitous maternal protein Zelda (Zld). Here we tested the role of Zld in Bcd-mediated binding and transcription. Removal of Zld function and mutations in Zld sites caused significant reductions in Bcd binding to known enhancers and variable effects on the activation and spatial positioning of Bcd-dependent expression patterns. Also, insertion of Zld sites converted one of six inactive fragments into a Bcd-responsive enhancer. Genome-wide binding experiments in zld mutants showed variable effects on Bcd-binding peaks, ranging from strong reductions to significantly enhanced levels of binding. Increases in Bcd binding caused the precocious Bcd-dependent activation of genes that are normally not expressed in early embryos, suggesting that Zld controls the genome-wide binding profile of Bcd at the qualitative level and is critical for selecting target genes for activation in the early embryo. These results underscore the importance of combinatorial binding in enhancer function and provide data that will help predict regulatory activities based on DNA sequence.

Mild nasal clefting may be predictive for ALX4 heterozygotes

Heterozygous loss-of-function mutations in ALX4 are responsible for enlarged parietal foramina, whereas patients with biallelic ALX4 mutations display a phenotypic spectrum of clinical findings, from mild to severe alopecia, cranium bifidum, hypertelorism, microphthalmia, with alar clefting being the pivotal sign in all affecteds. We report on four affected individuals in a three-generation family, displaying a phenotypic spectrum ranging from mild nasal clefting and broad columella to subtle changes in nasal configuration in addition to parietal foramina, caused by a novel ALX4 mutation (c.646C>G, p.Arg216Gly). This is the second report of a family showing vertical transmission of a dominant ALX4 mutation with facial involvement in addition to parietal foramina, mimicking mild recessive ALX4 phenotype. We discuss possible pathological mechanisms that may have lead to phenotypic variation in the family and challenges in genetic counseling.

Structural basis for DNA binding specificity by the auxin-dependent ARF transcription factors

Auxin regulates numerous plant developmental processes by controlling gene expression via a family of functionally distinct DNA-binding auxin response factors (ARFs), yet the mechanistic basis for generating specificity in auxin response is unknown. Here, we address this question by solving high-resolution crystal structures of the pivotal Arabidopsis developmental regulator ARF5/MONOPTEROS (MP), its divergent paralog ARF1, and a complex of ARF1 and a generic auxin response DNA element (AuxRE). We show that ARF DNA-binding domains also homodimerize to generate cooperative DNA binding, which is critical for in vivo ARF5/MP function. Strikingly, DNA-contacting residues are conserved between ARFs, and we discover that monomers have the same intrinsic specificity. ARF1 and ARF5 homodimers, however, differ in spacing tolerated between binding sites. Our data identify the DNA-binding domain as an ARF dimerization domain, suggest that ARF dimers bind complex sites as molecular calipers with ARF-specific spacing preference, and provide an atomic-scale mechanistic model for specificity in auxin response.

An improved predictive recognition model for Cys(2)-His(2) zinc finger proteins

Cys(2)-His(2) zinc finger proteins (ZFPs) are the largest family of transcription factors in higher metazoans. They also represent the most diverse family with regards to the composition of their recognition sequences. Although there are a number of ZFPs with characterized DNA-binding preferences, the specificity of the vast majority of ZFPs is unknown and cannot be directly inferred by homology due to the diversity of recognition residues present within individual fingers. Given the large number of unique zinc fingers and assemblies present across eukaryotes, a comprehensive predictive recognition model that could accurately estimate the DNA-binding specificity of any ZFP based on its amino acid sequence would have great utility. Toward this goal, we have used the DNA-binding specificities of 678 two-finger modules from both natural and artificial sources to construct a random forest-based predictive model for ZFP recognition. We find that our recognition model outperforms previously described determinant-based recognition models for ZFPs, and can successfully estimate the specificity of naturally occurring ZFPs with previously defined specificities.

Deciphering RNA regulatory elements in trypanosomatids: one piece at a time or genome-wide?

Morphological and metabolic changes in the life cycle of Trypanosoma brucei are accomplished by precise regulation of hundreds of genes. In the absence of transcriptional control, RNA-binding proteins (RBPs) shape the structure of gene regulatory maps in this organism, but our knowledge about their target RNAs, binding sites, and mechanisms of action is far from complete. Although recent technological advances have revolutionized the RBP-based approaches, the main framework for the RNA regulatory element (RRE)-based approaches has not changed over the last two decades in T. brucei. In this Opinion, after highlighting the current challenges in RRE inference, we explain some genome-wide solutions that can significantly boost our current understanding about gene regulatory networks in T. brucei.

Transcriptional enhancers: from properties to genome-wide predictions

Cellular development, morphology and function are governed by precise patterns of gene expression. These are established by the coordinated action of genomic regulatory elements known as enhancers or cis-regulatory modules. More than 30 years after the initial discovery of enhancers, many of their properties have been elucidated; however, despite major efforts, we only have an incomplete picture of enhancers in animal genomes. In this Review, we discuss how properties of enhancer sequences and chromatin are used to predict enhancers in genome-wide studies. We also cover recently developed high-throughput methods that allow the direct testing and identification of enhancers on the basis of their activity. Finally, we discuss recent technological advances and current challenges in the field of regulatory genomics.