Long-range formation of the Bicoid gradient requires multiple dynamic modes that spatially vary across the embryo

Morphogen gradients provide essential positional information to gene networks through their spatially heterogeneous distribution, yet how they form is still hotly contested, with multiple models proposed for different systems. Here, we focus on the transcription factor Bicoid (Bcd), a morphogen that forms an exponential gradient across the anterior-posterior (AP) axis of the early Drosophila embryo. Using fluorescence correlation spectroscopy we find there are spatial differences in Bcd diffusivity along the AP axis, with Bcd diffusing more rapidly in the posterior. We establish that such spatially varying differences in Bcd dynamics are sufficient to explain how Bcd can have a steep exponential gradient in the anterior half of the embryo and yet still have an observable fraction of Bcd near the posterior pole. In the nucleus, we demonstrate that Bcd dynamics are impacted by binding to DNA. Addition of the Bcd homeodomain to eGFP::NLS qualitatively replicates the Bcd concentration profile, suggesting this domain regulates Bcd dynamics. Our results reveal how a long-range gradient can form while retaining a steep profile through much of its range.

Revisiting bicoid function: complete inactivation reveals an additional fundamental role in Drosophila egg geometry specification

INTRODUCTION

The bicoid (bcd) gene in Drosophila has served as a paradigm for a morphogen in textbooks for decades. Discovered in 1986 as a mutation affecting anterior development in the embryo, its expression pattern as a protein gradient later confirmed the prediction from transplantation experiments. These experiments suggested that the protein fulfills the criteria of a true morphogen, with the existence of a homeodomain crucial for activation of genes along the anterior-posterior axis, based on the concentration of the morphogen. The bcd gene undergoes alternative splicing, resulting in, among other isoforms, a small and often neglected isoform with low abundance, which lacks the homeodomain, termed small bicoid (smbcd). Most importantly, all known classical strong bcd alleles used in the past to determine bcd function apparently do not affect the function of this isoform.

RESULTS

To overcome the uncertainty regarding which isoform regulates what, I removed the bcd locus entirely using CRISPR technology. bcdCRISPR eggs exhibited a short and round appearance. The phenotype could be ascribed to smbcd because all bcd alleles affecting the function of the major transcript, termed large bicoid (lgbcd) showed normally sized eggs. Several patterning genes for the embryo showed expression in the oocyte, and their expression patterns were altered in bcdCRISPR oocytes. In bcdCRISPR embryos, all downstream segmentation genes showed altered expression patterns, consistent with the expression patterns in "classical" alleles; however, due to the altered egg geometry resulting in fewer blastoderm nuclei, additional constraints came into play, further affecting their expression patterns.

CONCLUSIONS

This study unveils a novel and fundamental role of bcd in shaping the egg's geometry. This discovery demands a comprehensive revision of our understanding of this important patterning gene and prompts a reevaluation of past experiments conducted under the assumption that bcd mutants were bcdnull-mutants.

Shaping the scaling characteristics of gap gene expression patterns in Drosophila

How patterns are formed to scale with tissue size remains an unresolved problem. Here we investigate embryonic patterns of gap gene expression along the anterior-posterior (AP) axis in Drosophila. We use embryos that greatly differ in length and, importantly, possess distinct length-scaling characteristics of the Bicoid (Bcd) gradient. We systematically analyze the dynamic movements of gap gene expression boundaries in relation to both embryo length and Bcd input as a function of time. We document the process through which such dynamic movements drive both an emergence of a global scaling landscape and evolution of boundary-specific scaling characteristics. We show that, despite initial differences in pattern scaling characteristics that mimic those of Bcd in the anterior, such characteristics of final patterns converge. Our study thus partitions the contributions of Bcd input and regulatory dynamics inherent to the AP patterning network in shaping embryonic pattern's scaling characteristics.

On the importance of protein diffusion in biological systems: The example of the Bicoid morphogen gradient

Morphogens are proteins that form concentration gradients in embryos and developing tissues, where they act as postal codes, providing cells with positional information and allowing them to behave accordingly. Bicoid was the first discovered morphogen, and remains one of the most studied. It regulates segmentation in flies, forming a striking exponential gradient along the anterior-posterior axis of early Drosophila embryos, and activating the transcription of multiple target genes in a concentration-dependent manner. In this review, the work done by us and by others to characterize the mobility of Bicoid in D. melanogaster embryos is presented. The central role played by the diffusion of Bicoid in both the establishment of the gradient and the activation of target genes is discussed, and placed in the context of the need for these processes to be all at once rapid, precise and robust. The Bicoid system, and morphogen gradients in general, remain amongst the most amazing examples of the coexistence, often observed in living systems, of small-scale disorder and large-scale spatial order. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Concentration dependent chromatin states induced by the bicoid morphogen gradient

In Drosophila, graded expression of the maternal transcription factor Bicoid (Bcd) provides positional information to activate target genes at different positions along the anterior-posterior axis. We have measured the genome-wide binding profile of Bcd using ChIP-seq in embryos expressing single, uniform levels of Bcd protein, and grouped Bcd-bound targets into four classes based on occupancy at different concentrations. By measuring the biochemical affinity of target enhancers in these classes in vitro and genome-wide chromatin accessibility by ATAC-seq, we found that the occupancy of target sequences by Bcd is not primarily determined by Bcd binding sites, but by chromatin context. Bcd drives an open chromatin state at a subset of its targets. Our data support a model where Bcd influences chromatin structure to gain access to concentration-sensitive targets at high concentrations, while concentration-insensitive targets are found in more accessible chromatin and are bound at low concentrations. This may be a common property of developmental transcription factors that must gain early access to their target enhancers while the chromatin state of the genome is being remodeled during large-scale transitions in the gene regulatory landscape.

Improved readout precision of the Bicoid morphogen gradient by early decoding

Transcription factors (TFs) bind to specific DNA sequences to induce or repress gene expression. Expression levels can be tuned by changing TF concentrations, but the precision of such tuning is limited, since the fraction of time a TF occupies its binding site is subject to stochastic fluctuations. Bicoid (Bcd) is a TF that patterns the early Drosophila embryo by establishing an anterior-to-posterior concentration gradient and activating specific gene targets ("gap genes") in a concentration-dependent manner. Recently, the Bcd gradient and its in-vivo diffusion were quantified in live embryos, raising a quandary: the precision by which the Bcd target genes are defined (one-cell resolution) appeared to exceed the physical limits set by the stochastic binding of Bcd to DNA. We hypothesize that early readout of Bcd could account for the observed precision. Specifically, we consider the possibility that gap genes begin to be expressed earlier than typically measured experimentally, at a time when the distance between the nuclei is large. At this time, the difference in Bcd concentration between adjacent nuclei is large, enabling better tolerance for measurement imprecision. We show that such early decoding can indeed increase the accuracy of gap-gene expression, and that the initial pattern can be stabilized during subsequent divisions.

An unusual class of PITX2 mutations in Axenfeld-Rieger syndrome

BACKGROUND

Mutations in the PITX2 homeobox gene are known to contribute to Axenfeld-Rieger syndrome (ARS), an autosomal-dominant developmental disorder. Although most mutations are in the homeodomain and result in a loss of function, there is a growing subset in the C-terminal domain that has not yet been characterized. These mutations are of particular interest because the C-terminus has both inhibitory and stimulatory activities.

METHODS

In this study we used a combination of in vitro DNA binding and transfection reporter assays to investigate the fundamental issue of whether C-terminal mutations result in gain or loss of function at a cellular level.

RESULTS

We report a new frameshift mutation in the PITX2 allele that predicts a truncated protein lacking most of the C-terminal domain (D122FS). This newly reported mutant and another ARS C-terminal mutant (W133Stop) both have greater binding than wild-type to the bicoid element. Of interest, the mutants yielded approximately 5-fold greater activation of the prolactin promoter in CHO cells, even though the truncated proteins were expressed at lower levels than the wild-type protein. The truncated proteins also had greater than wild-type activity in 2 other cell lines, including the LS8 oral epithelial line that expresses the endogenous Pitx2 gene.

CONCLUSIONS

The results indicate that the PITX2 C-terminal domain has inhibitory activity and support the notion that ARS may also be caused by gain-of-function mutations.

The solution structure of the native K50 Bicoid homeodomain bound to the consensus TAATCC DNA-binding site

The solution structure of the homeodomain of the Drosophila morphogenic protein Bicoid (Bcd) complexed with a TAATCC DNA site is described. Bicoid is the only known protein that uses a homeodomain to regulate translation, as well as transcription, by binding to both RNA and DNA during early Drosophila development; in addition, the Bcd homeodomain can recognize an array of different DNA sites. The dual functionality and broad recognition capabilities signify that the Bcd homeodomain may possess unique structural/dynamic properties. Bicoid is the founding member of the K50 class of homeodomain proteins, containing a lysine residue at the critical 50th position (K50) of the homeodomain sequence, a residue required for DNA and RNA recognition; Bcd also has an arginine residue at the 54th position (R54), which is essential for RNA recognition. Bcd is the only known homeodomain with the K50/R54 combination of residues. The Bcd structure indicates that this homeodomain conforms to the conserved topology of the homeodomain motif, but exhibits a significant variation from other homeodomain structures at the end of helix 1. A key result is the observation that the side-chains of the DNA-contacting residues K50, N51 and R54 all show strong signs of flexibility in the protein-DNA interface. This finding is supportive of the adaptive-recognition theory of protein-DNA interactions.

A kinetic mechanism for Drosophila bicoid cooperative binding

The Bicoid (Bcd) protein is a concentration-dependent transcriptional activator in the embryo of Drosophila melanogaster. Bcd regulates the expression of the maternal and zygotic gene hunchback (hb) that shows a step-like-function expression pattern, in the anterior half of the egg. The regulatory region of hb contains six major binding sites for the Bcd protein, named A1, A2, A3 (strong sites), and X1, X2, X3 (weak sites). Cooperativity between Bcd molecules binding to the hb enhancer element has been characterized as an important mechanism for the step-like shape of hb anterior expression domain. The objective of the present report is to analyse the mechanism of this cooperative binding based on a reaction network model. Using this method we have analysed experimental results from the literature describing how the Bcd protein binds to hb enhancer elements containing the A1 or X1 site alone or these two sites together at wild type distance. This approach allows us to estimate the kinetic constants of protein-protein and protein-DNA interactions. Moreover our results suggest that binding of a Bcd dimer to the hb enhancer element is more stable than binding of a monomer. We propose a cooperative kinetic mechanism for binding of Bcd to the hb enhancer element: First, a monomer binds to the site with a relatively low affinity; after that, another monomer binds to the first one with higher affinity, generating a dimer bound to the site. This yet unreported monomer-monomer cooperative mechanism takes place for occupancy of either one-site or two-site enhancer elements. In addition, we find cooperativity between neighbor sites, as previously reported in the literature.

How to get ahead: the origin, evolution and function ofbicoid

In Drosophila, a Bcd protein gradient orchestrates patterning along the anteroposterior embryonic axis. However, studies of basal flies and other insects have revealed that bcd is a derived Hox3 gene found only in higher dipterans. To understand how bcd acquired its role in flies and how anteroposterior patterning mechanisms have evolved, I first review key features of bcd function in Drosophila: anterior localization and transcriptional and translation control of gene expression. I then discuss investigations of bcd in other higher dipterans that have provided insight into the evolution of regulatory interactions and the Bcd gradient. Finally, I review studies of Drosophila and other insects that address the evolution of bcd function and integration of bcd into ancestral regulatory mechanisms. I suggest further comparative studies may allow us to identify the intermediate steps in bcd evolution. This will make bcd a paradigm for the origin and evolution of genes and regulatory networks.

A composite motif of the Drosophila morphogenetic protein bicoid critical to transcription control

Bicoid is a molecular morphogen-controlling embryonic patterning in Drosophila. It is a homeodomain-containing protein that activates specific target genes during early embryogenesis. Our recent studies have identified a domain of Bcd located outside its homeodomain and referred to as a self-inhibitory domain that can dramatically repress its own ability to activate transcription. Here we present evidence that the self-inhibitory function is evolutionarily conserved. A systematic analysis of this domain reveals a composite 10-amino acid motif with interdigitating residues that regulate Bcd activity in opposite manners. Mutations within the Bcd motif can exert their respective effects when the self-inhibitory domain is grafted to an entirely heterologous activator, but they do not affect DNA binding in vitro or subcellular localization of Bcd in cells. We further show that the self-inhibitory domain of Bcd can interact with Sin3A, a component of the histone deacetylase co-repressor complex. Our study suggests that the activity of Bcd is intricately controlled by multiple mechanisms involving the actions of co-repressor proteins.

Enhancer sequences influence the role of the amino-terminal domain of bicoid in transcription

Bicoid (Bcd) is a Drosophila melanogaster morphogenetic gradient that controls embryonic patterning by activating target gene expression in a concentration-dependent manner. In this study we describe experiments to determine how different enhancers respond to Bcd distinctively, focusing on two natural Bcd-responsive enhancer elements, hunchback (hb) and knirps (kni). Our results show that, on the hb enhancer element, the amino-terminal domain of Bcd (residues 1 to 91) plays primarily an inhibitory role, whereas on the kni enhancer element this same Bcd domain plays a positive role at low protein concentrations. We further demonstrate that while the amino-terminal domain is largely dispensable for cooperative binding to the hb enhancer element, it is preferentially required for cooperative binding to the kni enhancer element. Alteration of the arrangement of Bcd binding sites in the kni enhancer element reduces the role of the amino-terminal domain in cooperative DNA binding but increases the effectiveness of the self-inhibitory function. In addition, elimination of symmetric pairs of Bcd binding sites in the kni enhancer element reduces both DNA binding and activation by Bcd. We propose that the amino-terminal domain of Bcd is an enhancer-specific switch that contributes to the protein's ability to activate different target genes in distinct manners.

Dominant negative dimerization of a mutant homeodomain protein in Axenfeld-Rieger syndrome

Axenfeld-Rieger syndrome is an autosomal-dominant disorder caused by mutations in the PITX2 homeodomain protein. We have studied the mechanism underlying the dominant negative K88E mutation, which occurs at position 50 of the homeodomain. By using yeast two-hybrid and in vitro pulldown assays, we have documented that PITX2a can form homodimers in the absence of DNA. Moreover, the K88E mutant had even stronger dimerization ability, primarily due to interactions involving the C-terminal region. Dimerization allowed cooperative binding of wild-type (WT) PITX2a to DNA containing tandem bicoid sites in a head-to-tail orientation (Hill coefficient, 1.73). In contrast, the WT-K88E heterodimer bound the tandem sites with greatly reduced cooperativity and decreased transactivation activity. To further explore the role of position 50 in PITX2a dimerization, we introduced a charge-conservative mutation of lysine to arginine (K88R). The K88R protein had greatly reduced binding to a TAATCC element and did not specifically bind any other TAATNN motif. Like K88E, K88R formed relatively stronger dimers with WT. As predicted by our model, the K88R protein acted in a dominant negative manner to suppress WT PITX2a activity. These results suggest that the position 50 residue in the PITX2 homeodomain plays an important role in both DNA binding and dimerization activities.

The activity of the Drosophila morphogenetic protein Bicoid is inhibited by a domain located outside its homeodomain

The Drosophila morphogenetic protein Bicoid (Bcd) is a homeodomain-containing activator that stimulates the expression of target genes during early embryonic development. We demonstrate that a small domain of Bcd located immediately N-terminally of the homeodomain represses its own activity in Drosophila cells. This domain, referred to as a self-inhibitory domain, works as an independent module that does not rely on any other sequences of Bcd and can repress the activity of heterologous activators. We further show that this domain of Bcd does not affect its properties of DNA binding or subcellular distribution. A Bcd derivative with point mutations in the self-inhibitory domain severely affects pattern formation and target gene expression in Drosophila embryos. We also provide evidence to suggest that the action of the self-inhibitory domain requires a Drosophila co-factor(s), other than CtBP or dSAP18. Our results suggest that proper action of Bcd as a transcriptional activator and molecular morphogen during embryonic development is dependent on the downregulation of its own activity through an interaction with a novel co-repressor(s) or complex(es).

Isolation of mutations that disrupt cooperative DNA binding by the Drosophila bicoid protein

Cooperative DNA binding is thought to contribute to the ability of the Drosophila melanogaster protein, Bicoid, to stimulate transcription of target genes in precise sub-domains within the embryo. As a first step toward testing this idea, we devised a genetic screen to isolate mutations in Bicoid that specifically disrupt cooperative interactions, but do not disrupt DNA recognition or transcription activation. The screen was carried out in Saccharomyces cerevisiae and 12 cooperativity mutants were identified. The mutations map across most of the Bicoid protein, with some located within the DNA-binding domain (homeodomain). Four homeodomain mutants were characterized in yeast and shown to activate a single-site reporter gene to levels comparable to that of wild-type, indicating that DNA binding per se is not affected. However, these mutants failed to show cooperative coupling between high and low-affinity sites, and showed reduced activation of a reporter gene carrying a natural Drosophila enhancer. Homology modeling indicated that none of the four mutations is in residues that contact DNA. Instead, these residues are likely to interact with other DNA-bound Bicoid monomers or other parts of the Bicoid protein. In vitro, the isolated homeodomains did not show strong cooperativity defects, supporting the idea that other regions of Bicoid are also important for cooperativity. This study describes the first systematic screen to identify cooperativity mutations in a eukaryotic DNA-binding protein.

Target selectivity of bicoid is dependent on nonconsensus site recognition and protein-protein interaction

We describe experiments to compare the activities of two Drosophila homeodomain proteins, Bicoid (Bcd) and an altered-specificity mutant of Fushi tarazu, Ftz(Q50K). Although the homeodomains of these proteins share a virtually indistinguishable ability to recognize a consensus Bcd site, only Bcd can activate transcription from natural enhancer elements when assayed in both yeast and Drosophila Schneider S2 cells. Our analysis of chimeric proteins suggests that both the homeodomain of Bcd and sequences outside the homeodomain contribute to its ability to recognize natural enhancer elements. We further show that, unlike the Bcd homeodomain, the Ftz(Q50K) homeodomain fails to recognize nonconsensus sites found in natural enhancer elements. The defect of a chimeric protein containing the homeodomain of Ftz(Q50K) in place of that of Bcd can be preferentially restored by converting the nonconsensus sites in natural enhancer elements to consensus sites. Our experiments suggest that the biological specificity of Bcd is determined by combinatorial contributions of two important mechanisms: the nonconsensus site recognition function conferred by the homeodomain and the cooperativity function conferred primarily by sequences outside the homeodomain. A systematic comparison of different assay methods and enhancer elements further suggests a fluid nature of the requirements for these two Bcd functions in target selection.

Reprogrammable recognition codes in bicoid homeodomain-DNA interaction

We describe experiments to determine how the homeodomain of the Drosophila morphogenetic protein Bicoid recognizes different types of DNA sequences found in natural enhancers. Our chemical footprint analyses reveal that the Bicoid homeodomain makes both shared and distinct contacts with a consensus site A1 (TAATCC) and a nonconsensus site X1 (TAAGCT). In particular, the guanine of X1 at position 4 (TAAGCT) is protected by Bicoid homeodomain. We provide further evidence suggesting that the unique arginine at position 54 (Arg 54) of the Bicoid homeodomain enables the protein to recognize X1 by specifically interacting with this position 4 guanine. We also describe experiments to analyze the contribution of artificially introduced Arg 54 to DNA recognition by other Bicoid-related homeodomains, including that from the human disease protein Pitx2. Our experiments demonstrate that the role of Arg 54 varies depending on the exact homeodomain framework and DNA sequences. Together, our results suggest that Bicoid and its related homeodomains utilize distinct recognition codes to interact with different DNA sequences, underscoring the need to study DNA recognition by Bicoid-class homeodomains in an individualized manner.

The homeodomain protein Barx2 contains activator and repressor domains and interacts with members of the CREB family

Barx1 and Barx2 are homeodomain proteins originally identified using regulatory elements of genes encoding certain cell adhesion molecules (CAMs). In the present study, we characterize regions of Barx2 that bind to regulatory elements of genes encoding three CAMs, L1, neuron-glia CAM (Ng-CAM), and neural CAM (N-CAM), and identify domains of Barx2 that regulate N-CAM transcription. The homeodomain of Barx2 was sufficient for binding to homeodomain binding sites (HBS) from all three CAM genes. The presence of a 17-amino acid Barx basic region resulted in a 2-fold decrease in binding to HBS sequences from the Ng-CAM and L1 genes, whereas it led to a 6.5-fold increase in binding to the HBS from the N-CAM promoter. Thus, the Barx basic region influences the strength and specificity of Barx2 binding to DNA. In co-transfection experiments, Barx2 repressed N-CAM promoter activity. A 24-residue N-terminal region of Barx2 was essential for repression. When this region was absent, Barx2 activated the N-CAM promoter. A 63-residue C-terminal domain was required for this activation. In GST pull-down experiments, Barx2 bound to proteins of the CREB family, CREB1 and ATF2. Overall, these findings provide a framework for understanding developmental and physiological contexts that influence repressor or activator functions of Barx2.

Chip interacts with diverse homeodomain proteins and potentiates bicoid activity in vivo

The Drosophila protein Chip potentiates activation by several enhancers and is required for embryonic segmentation. Chip and its mammalian homologs interact with and promote dimerization of nuclear LIM proteins. No known Drosophila LIM proteins, however, are required for segmentation, nor for expression of most genes known to be regulated by Chip. Here we show that Chip also interacts with diverse homeodomain proteins using residues distinct from those that interact with LIM proteins, and that Chip potentiates activity of one of these homeodomain proteins in Drosophila embryos and in yeast. These and other observations help explain the roles of Chip in segmentation and suggest a model to explain how Chip potentiates activation by diverse enhancers.

Drosophila Goosecoid requires a conserved heptapeptide for repression of paired-class homeoprotein activators

Goosecoid (Gsc) is a homeodomain protein expressed in the organizer region of vertebrate embryos. Its Drosophila homologue, D-Gsc, has been implicated in the formation of the Stomatogastric Nervous System. Although there are no apparent similarities between the phenotypes of mutations in the gsc gene in flies and mice, all known Gsc proteins can rescue dorsoanterior structures in ventralized Xenopus embryos. We describe how D-Gsc behaves as a transcriptional repressor in Drosophila cells, acting through specific palindromic HD binding sites (P3K). D-Gsc is a ‘passive repressor’ of activator homeoproteins binding to the same sites and an ‘active repressor’ of activators binding to distinct sites. In addition, D-Gsc is able to strongly repress transcription activated by Paired-class homeoproteins through P3K, via specific protein-protein interactions in what we define as ‘interactive repression’. This form of repression requires the short conserved GEH/eh-1 domain, also present in the Engrailed repressor. Although the GEH/eh-1 domain is necessary for rescue of UV-ventralized Xenopus embryos, it is dispensable for ectopic induction of Xlim-1 expression, demonstrating that this domain is not required for all Gsc functions in vivo. Interactive repression may represent specific interactions among Prd-class homeoproteins, several of which act early during development of invertebrate and vertebrate embryos.

Comparison of bicoid-dependent regulation of hunchback between Musca domestica and Drosophila melanogaster

Co-evolution between developmental regulatory elements is an important mechanism of evolution. This work compares the hunchback-bicoid interaction in the housefly Musca domestica with Drosophila melanogaster. The Musca HUNCHBACK protein is 66% conserved and partially rescues a hunchback mutant, yet the BICOID-dependent promoter (P2) of Musca hunchback is unexpectedly diverged from D. melanogaster. Introduced into D. melanogaster, this promoter drives a normal P2 pattern during the syncytial blastoderm stage but is expressed ectopically at the anterior pole of the embryo at later stages. We also report differences in the early expression of hunchback in Musca. We suggest that conservation of the morphogenetic function of bicoid in different sized embryos of higher diptera requires co-evolution of bicoid and its target binding sites.