Mechanism and Bicoid-dependent control of hairy stripe 7 expression in the posterior region of the Drosophila embryo

Global repression by tailless during segmentation

The orphan nuclear receptor Tailless (Tll) exhibits conserved roles in brain formation and maintenance that are shared, for example, with vertebrate orthologous forms (Tlx). However, the early expression of tll in two gap domains in the segmentation cascade of Drosophila is unusual even for most other insects. Here we investigate tll regulation on pair-rule stripes. With ectopic misexpression of tll we detected unexpected repression of almost all pair-rule stripes of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz). Examining Tll embryonic ChIP-chip data with regions mapped as Cis-Regulatory Modules (CRMs) of pair-rule stripes we verified Tll interactions to these regions. With the ChIP-chip data we also verified Tll interactions to the CRMs of gap domains and in the misexpression assay, Tll-mediated repression on Kruppel (Kr), kni (kni) and giant (gt) according to their differential sensitivity to Tll. These results with gap genes confirmed previous data from the literature and argue against indirect repression roles of Tll in the striped pattern. Moreover, the prediction of Tll binding sites in the CRMs of eve stripes and the mathematical modeling of their removal using an experimentally validated theoretical framework shows effects on eve stripes compatible with the absence of a repressor binding to the CRMs. In addition, modeling increased tll levels in the embryo results in the differential repression of eve stripes, agreeing well with the results of the misexpression assay. In genetic assays we investigated eve 5, that is strongly repressed by the ectopic domain and representative of more central stripes not previously implied to be under direct regulation of tll. While this stripe is little affected in tll-, its posterior border is expanded in gt- but detected with even greater expansion in gt-;tll-. We end up by discussing tll with key roles in combinatorial repression mechanisms to contain the expression of medial patterns of the segmentation cascade in the extremities of the embryo.

A timer gene network is spatially regulated by the terminal system in the Drosophila embryo

In insect embryos, anteroposterior patterning is coordinated by the sequential expression of the 'timer' genes caudal, Dichaete, and odd-paired, whose expression dynamics correlate with the mode of segmentation. In Drosophila, the timer genes are expressed broadly across much of the blastoderm, which segments simultaneously, but their expression is delayed in a small 'tail' region, just anterior to the hindgut, which segments during germband extension. Specification of the tail and the hindgut depends on the terminal gap gene tailless, but beyond this the regulation of the timer genes is poorly understood. We used a combination of multiplexed imaging, mutant analysis, and gene network modelling to resolve the regulation of the timer genes, identifying 11 new regulatory interactions and clarifying the mechanism of posterior terminal patterning. We propose that a dynamic Tailless expression gradient modulates the intrinsic dynamics of a timer gene cross-regulatory module, delineating the tail region and delaying its developmental maturation.

Evidence for the temporal regulation of insect segmentation by a conserved sequence of transcription factors

Long-germ insects, such as the fruit fly Drosophila melanogaster, pattern their segments simultaneously, whereas short-germ insects, such as the beetle Tribolium castaneum, pattern their segments sequentially, from anterior to posterior. While the two modes of segmentation at first appear quite distinct, much of this difference might simply reflect developmental heterochrony. We now show here that, in both Drosophila and Tribolium, segment patterning occurs within a common framework of sequential Caudal, Dichaete, and Odd-paired expression. In Drosophila these transcription factors are expressed like simple timers within the blastoderm, while in Tribolium they form wavefronts that sweep from anterior to posterior across the germband. In Drosophila, all three are known to regulate pair-rule gene expression and influence the temporal progression of segmentation. We propose that these regulatory roles are conserved in short-germ embryos, and that therefore the changing expression profiles of these genes across insects provide a mechanistic explanation for observed differences in the timing of segmentation. In support of this hypothesis we demonstrate that Odd-paired is essential for segmentation in Tribolium, contrary to previous reports.

Dense Bicoid hubs accentuate binding along the morphogen gradient

Morphogen gradients direct the spatial patterning of developing embryos; however, the mechanisms by which these gradients are interpreted remain elusive. Here we used lattice light-sheet microscopy to perform in vivo single-molecule imaging in early Drosophila melanogaster embryos of the transcription factor Bicoid that forms a gradient and initiates patterning along the anteroposterior axis. In contrast to canonical models, we observed that Bicoid binds to DNA with a rapid off rate throughout the embryo such that its average occupancy at target loci is on-rate-dependent. We further observed Bicoid forming transient "hubs" of locally high density that facilitate binding as factor levels drop, including in the posterior, where we observed Bicoid binding despite vanishingly low protein levels. We propose that localized modulation of transcription factor on rates via clustering provides a general mechanism to facilitate binding to low-affinity targets and that this may be a prevalent feature of other developmental transcription factors.

Genome-wide measurement of spatial expression in patterning mutants of Drosophila melanogaster

Patterning in the Drosophila melanogaster embryo is affected by multiple maternal factors, but the effect of these factors on spatial gene expression has not been systematically analyzed. Here we characterize the effect of the maternal factors Zelda, Hunchback and Bicoid by cryosectioning wildtype and mutant blastoderm stage embryos and sequencing mRNA from each slice. The resulting atlas of spatial gene expression highlights the intersecting roles of these factors in regulating spatial patterns, and serves as a resource for researchers studying spatial patterning in the early embryo. We identify a large number of genes with both expected and unexpected patterning changes, and through integrated analysis of transcription factor binding data identify common themes in genes with complex dependence on these transcription factors.

Maternal co-ordinate gene regulation and axis polarity in the scuttle fly Megaselia abdita

Axis specification and segment determination in dipteran insects are an excellent model system for comparative analyses of gene network evolution. Antero-posterior polarity of the embryo is established through systems of maternal morphogen gradients. In Drosophila melanogaster, the anterior system acts through opposing gradients of Bicoid (Bcd) and Caudal (Cad), while the posterior system involves Nanos (Nos) and Hunchback (Hb) protein. These systems act redundantly. Both Bcd and Hb need to be eliminated to cause a complete loss of polarity resulting in mirror-duplicated abdomens, so-called bicaudal phenotypes. In contrast, knock-down of bcd alone is sufficient to induce double abdomens in non-drosophilid cyclorrhaphan dipterans such as the hoverfly Episyrphus balteatus or the scuttle fly Megaselia abdita. We investigate conserved and divergent aspects of axis specification in the cyclorrhaphan lineage through a detailed study of the establishment and regulatory effect of maternal gradients in M. abdita. Our results show that the function of the anterior maternal system is highly conserved in this species, despite the loss of maternal cad expression. In contrast, hb does not activate gap genes in this species. The absence of this activatory role provides a precise genetic explanation for the loss of polarity upon bcd knock-down in M. abdita, and suggests a general scenario in which the posterior maternal system is increasingly replaced by the anterior one during the evolution of the cyclorrhaphan dipteran lineage.

The basic head-to-tail polarity of an animal is established very early in development. In dipteran insects (flies, midges, and mosquitoes), polarity is established with the help of so-called morphogen gradients. Morphogens are regulatory proteins that are distributed as a concentration gradient, often involving diffusion from a localised source. This graded distribution then leads to the concentration-dependent activation of different target genes along the embryo’s axis. We examine this process, which differs to a surprising extent between dipteran species, in the scuttle fly Megaselia abdita, and compare our results to the model organism Drosophila melanogaster. In this way, we not only gain insights into how the mechanisms that establish polarity function differently in different species, but also how the system has evolved since these two flies shared a common ancestor. Specifically, we pin down the main difference between Drosophila and Megaselia in the altered function of the maternal Hunchback morphogen gradient, which activates target genes in the former, but not the latter species, where it has been completely replaced by the Bicoid morphogen during evolution.

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.

Quantitative modeling of a gene's expression from its intergenic sequence

Modeling a gene's expression from its intergenic locus and trans-regulatory context is a fundamental goal in computational biology. Owing to the distributed nature of cis-regulatory information and the poorly understood mechanisms that integrate such information, gene locus modeling is a more challenging task than modeling individual enhancers. Here we report the first quantitative model of a gene's expression pattern as a function of its locus. We model the expression readout of a locus in two tiers: 1) combinatorial regulation by transcription factors bound to each enhancer is predicted by a thermodynamics-based model and 2) independent contributions from multiple enhancers are linearly combined to fit the gene expression pattern. The model does not require any prior knowledge about enhancers contributing toward a gene's expression. We demonstrate that the model captures the complex multi-domain expression patterns of anterior-posterior patterning genes in the early Drosophila embryo. Altogether, we model the expression patterns of 27 genes; these include several gap genes, pair-rule genes, and anterior, posterior, trunk, and terminal genes. We find that the model-selected enhancers for each gene overlap strongly with its experimentally characterized enhancers. Our findings also suggest the presence of sequence-segments in the locus that would contribute ectopic expression patterns and hence were "shut down" by the model. We applied our model to identify the transcription factors responsible for forming the stripe boundaries of the studied genes. The resulting network of regulatory interactions exhibits a high level of agreement with known regulatory influences on the target genes. Finally, we analyzed whether and why our assumption of enhancer independence was necessary for the genes we studied. We found a deterioration of expression when binding sites in one enhancer were allowed to influence the readout of another enhancer. Thus, interference between enhancer activities was a possible factor necessitating enhancer independence in our model.

Adaptation of the length scale and amplitude of the Bicoid gradient profile to achieve robust patterning in abnormally large Drosophila melanogaster embryos

The formation of patterns that are proportional to the size of the embryo is an intriguing but poorly understood feature of development. Molecular mechanisms controlling such proportionality, or scaling, can be probed through quantitative interrogations of the properties of morphogen gradients that instruct patterning. Recent studies of the Drosophila morphogen gradient Bicoid (Bcd), which is required for anterior-posterior (AP) patterning in the early embryo, have uncovered two distinct ways of scaling. Whereas between-species scaling is achieved by adjusting the exponential shape characteristic of the Bcd gradient profile, namely, its length scale or length constant (λ), within-species scaling is achieved through adjusting the profile's amplitude, namely, the Bcd concentration at the anterior (B0). Here, we report a case in which Drosophila melanogaster embryos exhibit Bcd gradient properties uncharacteristic of their size. The embryos under investigation were from a pair of inbred lines that had been artificially selected for egg size extremes. We show that B0 in the large embryos is uncharacteristically low but λ is abnormally extended. Although the large embryos have more total bcd mRNA than their smaller counterparts, as expected, its distribution is unusually broad. We show that the large and small embryos develop gene expression patterns exhibiting boundaries that are proportional to their respective lengths. Our results suggest that the large-egg inbred line has acquired compensating properties that counteract the extreme length of the embryos to maintain Bcd gradient properties necessary for robust patterning. Our study documents, for the first time to our knowledge, a case of within-species Bcd scaling achieved through adjusting the gradient profile's exponential shape characteristic, illustrating at a molecular level how a developmental system can follow distinct operational paths towards the goal of robust and scaled patterning.

Lack of tailless leads to an increase in expression variability in Drosophila embryos

Developmental processes are robust, or canalised: dynamic patterns of gene expression across space and time are regulated reliably and precisely in the presence of genetic and environmental perturbations. It remains unclear whether canalisation relies on specific regulatory factors (such as heat-shock proteins), or whether it is based on more general redundancy and distributed robustness at the network level. The latter explanation implies that mutations in many regulatory factors should exhibit loss of canalisation. Here, we present a quantitative characterisation of segmentation gene expression patterns in mutants of the terminal gap gene tailless (tll) in Drosophila melanogaster. Our analysis provides new insights into the dynamic mechanisms underlying gap gene regulation, and reveals significantly increased variability of gene expression in the mutant compared to the wild-type background. We show that both position and timing of posterior segmentation gene expression domains vary strongly from embryo-to-embryo in tll mutants. This variability must be caused by a vulnerability in the regulatory system which is hidden or buffered in the wild-type, but becomes uncovered by the deletion of tll. Our analysis provides evidence that loss of canalisation in mutants could be more widespread than previously thought.

Cdx1 refines positional identity of the vertebrate hindbrain by directly repressing Mafb expression

An interplay of transcription factors interprets signalling pathways to define anteroposterior positions along the vertebrate axis. In the hindbrain, these transcription factors prompt the position-appropriate appearance of seven to eight segmental structures, known as rhombomeres (r1-r8). The evolutionarily conserved Cdx caudal-type homeodomain transcription factors help specify the vertebrate trunk and tail but have not been shown to directly regulate hindbrain patterning genes. Mafb (Kreisler, Krml1, valentino), a basic domain leucine zipper transcription factor, is required for development of r5 and r6 and is the first gene to show restricted expression within these two segments. The homeodomain protein vHnf1 (Hnf1b) directly activates Mafb expression. vHnf1 and Mafb share an anterior expression limit at the r4/r5 boundary but vHnf1 expression extends beyond the posterior limit of Mafb and, therefore, cannot establish the posterior Mafb expression boundary. Upon identifying regulatory sequences responsible for posterior Mafb repression, we have used in situ hybridization, immunofluorescence and chromatin immunoprecipitation (ChIP) analyses to determine that Cdx1 directly inhibits early Mafb expression in the neural tube posterior of the r6/r7 boundary, which is the anteriormost boundary of Cdx1 expression in the hindbrain. Cdx1 dependent repression of Mafb is transient. After the 10-somite stage, another mechanism acts to restrict Mafb expression in its normal r5 and r6 domain, even in the absence of Cdx1. Our findings identify Mafb as one of the earliest direct targets of Cdx1 and show that Cdx1 plays a direct role in early hindbrain patterning. Thus, just as Cdx2 and Cdx4 govern the trunk-to-tail transition, Cdx1 may regulate the hindbrain-to-spinal cord transition.

Non-additive interactions involving two distinct elements mediate sloppy-paired regulation by pair-rule transcription factors

The relatively simple combinatorial rules responsible for establishing the initial metameric expression of sloppy-paired-1 (slp1) in the Drosophila blastoderm embryo make this system an attractive model for investigating the mechanism of regulation by pair-rule transcription factors. This investigation of slp1 cis-regulatory architecture identifies two distinct elements, a proximal early stripe element (PESE) and a distal early stripe element (DESE) located from -3.1kb to -2.5kb and from -8.1kb to -7.1kb upstream of the slp1 promoter, respectively, that mediate this early regulation. The proximal element expresses only even-numbered stripes and mediates repression by Even-skipped (Eve) as well as by the combination of Runt and Fushi-tarazu (Ftz). A 272 basepair sub-element of PESE retains an Eve-dependent repression, but is expressed throughout the even-numbered parasegments due to the loss of repression by Runt and Ftz. In contrast, the distal element expresses both odd and even-numbered stripes and also drives inappropriate expression in the anterior half of the odd-numbered parasegments due to an inability to respond to repression by Eve. Importantly, a composite reporter gene containing both early stripe elements recapitulates pair-rule gene-dependent regulation in a manner beyond what is expected from combining their individual patterns. These results indicate that interactions involving distinct cis-elements contribute to the proper integration of pair-rule regulatory information. A model fully accounting for these results proposes that metameric slp1 expression is achieved through the Runt-dependent regulation of interactions between these two pair-rule response elements and the slp1 promoter.

Association of orthodenticle with natural variation for early embryonic patterning in Drosophila melanogaster

Although it is well established that cis-acting regulatory variation contributes to morphological evolution between species, few concrete examples of polymorphism affecting developmental patterning within species have been demonstrated. Early embryogenesis in Drosophila is initiated by a gradient of Bicoid morphogen activity that results in differential expression of multiple target genes. In a screen for genetic variation affecting this process, we surveyed 96 wild-type lines of Drosophila melanogaster for polymorphisms in binding sites within 16 Bicoid cis-regulatory response elements. One common polymorphism in the orthodenticle (otd) early head enhancer is associated with a complex series of indels/substitutions that define two distinct haplotypes. The middle region of this enhancer exhibits an unusual pattern of nucleotide diversity that does not easily fit into standard models of selection and demography. Population Gene Expression Maps, generated by extracting binary expression profiles from normalized embryo images, revealed a ventral reduction of otd transcript abundance in one of the haplotypes that was recapitulated in expression of transgenic constructs containing the two alleles. We thus demonstrate that even a process as robust as early developmental patterning is affected by standing genetic variation, intriguingly involving otd, whose morphogenetic function bicoid is thought to have displaced during dipteran evolution.

Predicting expression patterns from regulatory sequence in Drosophila segmentation

The establishment of complex expression patterns at precise times and locations is key to metazoan development, yet a mechanistic understanding of the underlying transcription control networks is still missing. Here we describe a novel thermodynamic model that computes expression patterns as a function of cis-regulatory sequence and of the binding-site preferences and expression of participating transcription factors. We apply this model to the segmentation gene network of Drosophila melanogaster and find that it predicts expression patterns of cis-regulatory modules with remarkable accuracy, demonstrating that positional information is encoded in the regulatory sequence and input factor distribution. Our analysis reveals that both strong and weaker binding sites contribute, leading to high occupancy of the module DNA, and conferring robustness against mutation; short-range homotypic clustering of weaker sites facilitates cooperative binding, which is necessary to sharpen the patterns. Our computational framework is generally applicable to most protein-DNA interaction systems.

cis-Decoder discovers constellations of conserved DNA sequences shared among tissue-specific enhancers

: The use of cis-Decoder, a new tool for discovery of conserved sequence elements that are shared between similarly regulating enhancers, suggests that enhancers use overlapping repertoires of highly conserved core elements.

A systematic approach is described for analysis of evolutionarily conserved cis-regulatory DNA using cis-Decoder, a tool for discovery of conserved sequence elements that are shared between similarly regulated enhancers. Analysis of 2,086 conserved sequence blocks (CSBs), identified from 135 characterized enhancers, reveals most CSBs consist of shorter overlapping/adjacent elements that are either enhancer type-specific or common to enhancers with divergent regulatory behaviors. Our findings suggest that enhancers employ overlapping repertoires of highly conserved core elements.

AcaudalmRNA gradient controls posterior development in the waspNasonia

One of the earliest steps of embryonic development is the establishment of polarity along the anteroposterior axis. Extensive studies of Drosophila embryonic development have elucidated mechanisms for establishing polarity, while studies with other model systems have found that many of these molecular components are conserved through evolution. One exception is Bicoid, the master organizer of anterior development in Drosophila and higher dipterans, which is not conserved. Thus, the study of anteroposterior patterning in insects that lack Bicoid can provide insight into the evolution of the diversity of body plan patterning networks. To this end, we have established the long germ parasitic wasp Nasonia vitripennis as a model for comparative studies with Drosophila.Here we report that, in Nasonia, a gradient of localized caudal mRNA directs posterior patterning, whereas, in Drosophila, the gradient of maternal Caudal protein is established through translational repression by Bicoid of homogeneous caudalmRNA. Loss of caudal function in Nasonia results in severe segmentation defects. We show that Nasonia caudal is an activator of gap gene expression that acts far towards the anterior of the embryo, placing it atop a cascade of early patterning. By contrast, activation of gap genes in flies relies on redundant functions of Bicoid and Caudal, leading to a lack of dramatic action on gap gene expression: caudal instead plays a limited role as an activator of pair-rule gene expression. These studies,together with studies in short germ insects, suggest that caudal is an ancestral master organizer of patterning, and that its role has been reduced in higher dipterans such as Drosophila.

The tailless nuclear receptor acts as a dedicated repressor in the early Drosophila embryo

Tailless is an orphan nuclear receptor that controls terminal body patterning in Drosophila. Genetic analyses have revealed both positive and negative regulatory interactions of Tailless with various target genes, leading to the idea that, like many other nuclear receptors, Tailless mediates both activation and repression of transcription. In this paper, we have examined the consequences of converting Tailless into an obligate repressor and compared the activities of the resulting protein with those of wild-type Tailless. We find that this repressor form of Tailless behaves like the intact protein in gain- and loss-of-function experiments, being sufficient to support normal embryonic development and establish accurate patterns of gene expression even for positive Tailless targets such as hunchback and brachyenteron. This suggests that Tailless functions exclusively as a transcriptional repressor in the embryo and that the observed positive interactions of Tailless with specific targets are secondary effects involving repression of repressors. We provide evidence that knirps is one such repressor gene acting between Tailless and its indirect positive targets. Finally, our results indicate that Tailless exerts an active mechanism of repression via its ligand-binding domain and that this activity is largely independent of the activation function 2 (AF2) motif characteristic of most nuclear receptors.

The interpretation of morphogen gradients

Morphogens act as graded positional cues that control cell fate specification in many developing tissues. This concept, in which a signalling gradient regulates differential gene expression in a concentration-dependent manner, provides a basis for understanding many patterning processes. It also raises several mechanistic issues, such as how responding cells perceive and interpret the concentration-dependent information provided by a morphogen to generate precise patterns of gene expression and cell differentiation in developing tissues. Here, we review recent work on the molecular features of morphogen signalling that facilitate the interpretation of graded signals and attempt to identify some emerging common principles.

Localized maternal orthodenticle patterns anterior and posterior in the long germ wasp Nasonia

The Bicoid (Bcd) gradient in Drosophila has long been a model for the action of a morphogen in establishing embryonic polarity. However, it is now clear that bcd is a unique feature of higher Diptera. An evolutionarily ancient gene, orthodenticle (otd), has a bcd-like role in the beetle Tribolium. Unlike the Bcd gradient, which arises by diffusion of protein from an anteriorly localized messenger RNA, the Tribolium Otd gradient forms by translational repression of otd mRNA by a posteriorly localized factor. These differences in gradient formation are correlated with differences in modes of embryonic patterning. Drosophila uses long germ embryogenesis, where the embryo derives from the entire anterior-posterior axis, and all segments are patterned at the blastoderm stage, before gastrulation. In contrast, Tribolium undergoes short germ embryogenesis: the embryo arises from cells in the posterior of the egg, and only anterior segments are patterned at the blastoderm stage, with the remaining segments arising after gastrulation from a growth zone. Here we describe the role of otd in the long germband embryo of the wasp Nasonia vitripennis. We show that Nasonia otd maternal mRNA is localized at both poles of the embryo, and resulting protein gradients pattern both poles. Thus, localized Nasonia otd has two major roles that allow long germ development. It activates anterior targets at the anterior of the egg in a manner reminiscent of the Bcd gradient, and it is required for pre-gastrulation expression of posterior gap genes.

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.

cis-regulatory logic of short-range transcriptional repression in Drosophila melanogaster

Bioinformatics analysis of transcriptional control is guided by knowledge of the characteristics of cis-regulatory regions or enhancers. Features such as clustering of binding sites and co-occurrence of binding sites have aided enhancer identification, but quantitative predictions of enhancer function are not yet generally feasible. To facilitate the analysis of regulatory sequences in Drosophila melanogaster, we identified quantitative parameters that affect the activity of short-range transcriptional repressors, proteins that play key roles in development. In addition to the previously noted distance dependence, repression is strongly influenced by the stoichiometry, affinity, spacing, and arrangement of activator binding sites. Repression is insensitive to the type of activation domain, suggesting that short-range repression may primarily affect activators at the level of DNA binding. The activity of several short-range, but not long-range, repressors is circumscribed by the same quantitative parameters. This cis-regulatory "grammar" may aid the identification of enhancers regulated by short-range repressors and facilitate bioinformatic prediction of the functional output of transcriptional regulatory sequences.

The role of binding site cluster strength in Bicoid-dependent patterning in Drosophila

The maternal morphogen Bicoid (Bcd) is distributed in an embryonic gradient that is critical for patterning the anterior-posterior (AP) body plan in Drosophila. Previous work identified several target genes that respond directly to Bcd-dependent activation. Positioning of these targets along the AP axis is thought to be controlled by cis-regulatory modules (CRMs) that contain clusters of Bcd-binding sites of different "strengths." Here we use a combination of Bcd-site cluster analysis and evolutionary conservation to predict Bcd-dependent CRMs. We tested 14 predicted CRMs by in vivo reporter gene assays; 11 show Bcd-dependent activation, which brings the total number of known Bcd target elements to 21. Some CRMs drive expression patterns that are restricted to the most anterior part of the embryo, whereas others extend into middle and posterior regions. However, we do not detect a strong correlation between AP position of target gene expression and the strength of Bcd site clusters alone. Rather, we find that binding sites for other activators, including Hunchback and Caudal correlate with CRM expression in middle and posterior body regions. Also, many Bcd-dependent CRMs contain clusters of sites for the gap protein Kruppel, which may limit the posterior extent of activation by the Bcd gradient. We propose that the key design principle in AP patterning is the differential integration of positive and negative transcriptional information at the level of individual CRMs for each target gene.

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.

Transcriptional control in the segmentation gene network of Drosophila

The segmentation gene network of Drosophila consists of maternal and zygotic factors that generate, by transcriptional (cross-) regulation, expression patterns of increasing complexity along the anterior-posterior axis of the embryo. Using known binding site information for maternal and zygotic gap transcription factors, the computer algorithm Ahab recovers known segmentation control elements (modules) with excellent success and predicts many novel modules within the network and genome-wide. We show that novel module predictions are highly enriched in the network and typically clustered proximal to the promoter, not only upstream, but also in intronic space and downstream. When placed upstream of a reporter gene, they consistently drive patterned blastoderm expression, in most cases faithfully producing one or more pattern elements of the endogenous gene. Moreover, we demonstrate for the entire set of known and newly validated modules that Ahab's prediction of binding sites correlates well with the expression patterns produced by the modules, revealing basic rules governing their composition. Specifically, we show that maternal factors consistently act as activators and that gap factors act as repressors, except for the bimodal factor Hunchback. Our data suggest a simple context-dependent rule for its switch from repressive to activating function. Overall, the composition of modules appears well fitted to the spatiotemporal distribution of their positive and negative input factors. Finally, by comparing Ahab predictions with different categories of transcription factor input, we confirm the global regulatory structure of the segmentation gene network, but find odd skipped behaving like a primary pair-rule gene. The study expands our knowledge of the segmentation gene network by increasing the number of experimentally tested modules by 50%. For the first time, the entire set of validated modules is analyzed for binding site composition under a uniform set of criteria, permitting the definition of basic composition rules. The study demonstrates that computational methods are a powerful complement to experimental approaches in the analysis of transcription networks.

Genome-wide mapping of in vivo targets of the Drosophila transcription factor Kruppel

Krüppel (Kr), a member of the gap class of Drosophila segmentation genes, encodes a DNA binding zinc finger-type transcription factor. In addition to its segmentation function at the blastoderm stage, Krüppel also plays a critical role in organ formation during later stages of embryogenesis. To systematically identify in vivo target genes of Krüppel, we isolated DNA fragments from the Krüppel-associated portion of chromatin and used them to find and map Krüppel-dependent cis-acting regulatory sites in the Drosophila genome. We show that Krüppel binding sites are not enriched in Krüppel-associated chromatin and that the clustering of Krüppel binding sites, as found in the cis-acting elements of Krüppel-dependent segmentation genes used for in silico searches of Krüppel target genes, is not a prerequisite for the in vivo binding of Krüppel to its regulatory elements. Results obtained with the newly identified target gene ken and barbie (ken) indicate that Krüppel represses transcription and thereby restricts the spatial expression pattern of ken during blastoderm and gastrulation.

Information display by transcriptional enhancers

Transcriptional enhancers integrate positional and temporal information to regulate the complex expression of developmentally controlled genes. Current models suggest that enhancers act as computational devices, receiving multiple inputs from activators and repressors and resolving them into a single positive or a negative signal that is transmitted to the basal transcriptional machinery. We show that a simple, compact enhancer is capable of representing both repressed and activated states at the same time and in the same nucleus. This finding suggests that closely apposed factor binding sites, situated within compact cis-elements, can be independently interpreted by the transcriptional machinery, possibly through successive enhancer-promoter interactions. These results provide clear evidence that the computational functions usually ascribed to the enhancer itself are actually shared with the basal machinery. In contrast to the autonomous computer model of enhancer function, an information-display or 'billboard' model of enhancer activity may better describe many developmentally regulated transcriptional enhancers.

CtBP-independent repression in the Drosophila embryo

There are three mechanisms of transcriptional repression in eukaryotes. The first is quenching, whereby repressors and activators co-occupy closely linked sites and then the repressor inhibits adjacent activators. The second is direct repression, in which repressors block the function of the core transcription complex. The third is competition, in which repressors compete with activators for a common DNA-binding site. Previous studies have shown that the Drosophila melanogaster CtBP corepressor (dCtBP) is essential for the quenching activity of three short-range sequence-specific repressors in the early Drosophila embryo: Krüppel, Knirps, and Snail. Here we demonstrate that dCtBP is dispensable for target enhancers that contain overlapping activator and repressor binding sites. However, it is essential when Krüppel and Knirps repressor sites do not overlap activator sites but are instead located adjacent to either activators or the core promoter. These findings provide evidence that competition is distinct from quenching and direct repression. Quenching and direct repression depend on dCtBP, whereas competition does not.

Ingeneue: a versatile tool for reconstituting genetic networks, with examples from the segment polarity network

Here we describe a software tool for synthesizing molecular genetic data into models of genetic networks. Our software program Ingeneue, written in Java, lets the user quickly turn a map of a genetic network into a dynamical model consisting of a set of ordinary differential equations. We developed Ingeneue as part of an ongoing effort to explore the design and evolvability of genetic networks. Ingeneue has three principal advantages over other available mathematical software: it automates instantiation of the same network model in each cell in a 2-D sheet of cells; it constructs model equations from pre-made building blocks corresponding to common biochemical processes; and it automates searches through parameter space, sensitivity analyses, and other common tasks. Here we discuss the structure of the software and some of the issues we have dealt with. We conclude with some examples of results we have achieved with Ingeneue for the Drosophila segment polarity network.

An anterior function for the Drosophila posterior determinant Pumilio

Bicoid is a key determinant of anterior Drosophila development. We demonstrate that the prototypical Puf protein Pumilio temporally regulates bicoid (bcd) mRNA translation via evolutionarily conserved Nanos response elements (NRE) in its 3′UTR. Disruption of Pumilio-bcd mRNA interaction by either Pumilio or bcd NRE mutations caused delayed bcd mRNA deadenylation and stabilization, resulting in protracted Bicoid protein expression during embryogenesis. Phenotypically, embryos from transgenic mothers that harbor bcd NRE mutations exhibited dominant anterior patterning defects and we discovered similar head defects in embryos from pum– mothers. Hence, Pumilio is required for normal anterior development. Since bcd mRNA resides outside the posterior gradient of the canonical partner of Pumilio, Nanos, our data suggest that Pumilio can recruit different partners to specifically regulate distinct mRNAs.

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).

Extraction of functional binding sites from unique regulatory regions: the Drosophila early developmental enhancers

The early developmental enhancers of Drosophila melanogaster comprise one of the most sophisticated regulatory systems in higher eukaryotes. An elaborate code in their DNA sequence translates both maternal and early embryonic regulatory signals into spatial distribution of transcription factors. One of the most striking features of this code is the redundancy of binding sites for these transcription factors (BSTF). Using this redundancy, we explored the possibility of predicting functional binding sites in a single enhancer region without any prior consensus/matrix description or evolutionary sequence comparisons. We developed a conceptually simple algorithm, Scanseq, that employs an original statistical evaluation for identifying the most redundant motifs and locates the position of potential BSTF in a given regulatory region. To estimate the biological relevance of our predictions, we built thorough literature-based annotations for the best-known Drosophila developmental enhancers and we generated detailed distribution maps for the most robust binding sites. The high statistical correlation between the location of BSTF in these experiment-based maps and the location predicted in silico by Scanseq confirmed the relevance of our approach. We also discuss the definition of true binding sites and the possible biological principles that govern patterning of regulatory regions and the distribution of transcriptional signals.

Macro-evolution of the hairy enhancer in Drosophila species

It has been suggested that many of the changes in the developmental program might be in the cis-acting promoters and enhancer regions. Here I study the macro-evolutionary changes of an enhancer region for the early developmental gene hairy in Drosophila melanogaster, D. simulans, D. pseudoobscura, D. willistoni, D. nebulosa, D. hydei, and D. virilis. The enhancer region is characterized by small, highly conserved blocks interspersed among highly variable regions. Nevertheless, species phylogenies constructed by the enhancer sequences agree with the widely accepted phylogeny of these species. The evolution of the variable regions is consistent with a molecular clock, while the evolution of the conserved blocks is significantly different from a clock. In particular, the D. pseudoobscura lineage shows the highest degree of species-specific change consistent with changes in expression timing reported in an earlier study. It has been suggested that the variation in sequence length between highly conserved blocks may play a role in the coordination of regulatory processes, such as protein-protein interactions; thus, stabilizing selection has been suggested to act on the length variations. Here I develop a test for stabilizing selection on length variation and show that the hairy enhancer does not show statistically significant evidence for stabilizing selection. J. Exp. Zool. (Mol. Dev. Evol.) 291:175-185, 2001.

How gap genes make their domains: An analytical study based on data driven approximations

We consider a mathematical formulation of the problem of protein production during segment determination in the Drosophila blastoderm, together with some preliminary results of its analytical study. We reformulate the spatial difference equations as a set of nonlinear partial differential equations and obtain their dimensionless form in the continuum limit. Using previous results obtained by the gene circuit method, we find an asymptotic statement of the problem with a small parameter. Some results of the comparison method applied to the model are obtained, and exact stationary upper solutions are derived. They exhibit distinctive features of localized bell-shaped structures. (c) 2001 American Institute of Physics.

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.

Function of bicoid and hunchback homologs in the basal cyclorrhaphan fly Megaselia (Phoridae)

The Drosophila gene bicoid functions at the beginning of a gene cascade that specifies anterior structures in the embryo. Its transcripts are localized at the anterior pole of the oocyte, giving rise to a Bicoid protein gradient, which regulates the spatially restricted expression of target genes along the anterior-posterior axis of the embryo in a concentration-dependent manner. The morphogen function of Bicoid requires the coactivity of the zinc finger transcription factor Hunchback, which is expressed in a Bicoid-dependent fashion in the anterior half of the embryo. Whereas hunchback is conserved throughout insects, bicoid homologs are known only from cyclorrhaphan flies. Thus far, identification of hunchback and bicoid homologs rests only on sequence comparison. In this study, we used double-stranded RNA interference (RNAi) to address the function of bicoid and hunchback homologs in embryos of the lower cyclorrhaphan fly Megaselia abdita (Phoridae). Megaselia-hunchback RNAi causes hunchback-like phenotypes as observed in Drosophila, but Megaselia-bicoid RNAi causes phenotypes different from corresponding RNAi experiments in Drosophila and bicoid mutant embryos. Megaselia-bicoid is required not only for the head and thorax but also for the development of four abdominal segments. This difference between Megaselia and Drosophila suggests that the range of functional bicoid activity has been reduced in higher flies.

hairy stripe 7 element mediates activation and repression in response to different domains and levels of Krüppel in the Drosophila embryo

The Drosophila gap gene Krüppel (Kr) encodes a zinc finger-type transcription factor required for controlling the spatial expression of other segmentation genes during early blastoderm stage. Here we show that two independent and transferable repressor domains of Krüppel act to control expression of the pair-rule gene hairy, and that the minimal cis-acting element of hairy stripe7 (h7) mediates either Krüppel-dependent activation or repression in different regions of the blastoderm embryo. The C-terminal region of Krüppel which encompasses the predominant repressor domain is not essential for activation, but is required to fully suppress h7-mediated transcription in response to high levels of Krüppel activity. This domain contains an interaction motif for dCtBP, a homologue of the human co-repressor CtBP. dCtBP activity is, however, dispensable for Krüppel-mediated repression in the embryo since Krüppel-mediated repression functions in the absence of dCtBP. Possible modes of h7-mediated gene regulation in response to the different domains and levels of Krüppel are discussed.

The Cdx-1 and Cdx-2 homeobox genes in the intestine

The past years have witnessed an increasing number of reports relative to homeobox genes in endoderm-derived tissues. In this review, we focus on the caudal-related Cdx-1 and Cdx-2 homeobox genes to give an overview of the in vivo, in vitro, and ex vivo approaches that emphasize their primary role in intestinal development and in the control of intestinal cell proliferation, differentiation, and identity. The participation of these genes in colon tumorigenesis and their identification as important actors of the oncogenic process are also discussed.

Transducing positional information to the Hox genes: critical interaction of cdx gene products with position-sensitive regulatory elements

Studies of pattern formation in the vertebrate central nervous system indicate that anteroposterior positional information is generated in the embryo by signalling gradients of an as yet unknown nature. We searched for transcription factors that transduce this information to the Hox genes. Based on the assumption that the activity levels of such factors might vary with position along the anteroposterior axis, we devised an in vivo assay to detect responsiveness of cis-acting sequences to such differentially active factors. We used this assay to analyze a Hoxb8 regulatory element, and detected the most pronounced response in a short stretch of DNA containing a cluster of potential CDX binding sites. We show that differentially expressed DNA binding proteins are present in gastrulating embryos that bind to these sites in vitro, that cdx gene products are among these, and that binding site mutations that abolish binding of these proteins completely destroy the ability of the regulatory element to drive regionally restricted expression in the embryo. Finally, we show that ectopic expression of cdx gene products anteriorizes expression of reporter transgenes driven by this regulatory element, as well as that of the endogenous Hoxb8 gene, in a manner that is consistent with them being essential transducers of positional information. These data suggest that, in contrast to Drosophila Caudal, vertebrate cdx gene products transduce positional information directly to the Hox genes, acting through CDX binding sites in their enhancers. This may represent the ancestral mode of action of caudal homologues, which are involved in anteroposterior patterning in organisms with widely divergent body plans and modes of development.

Cooperative DNA-binding by Bicoid provides a mechanism for threshold-dependent gene activation in the Drosophila embryo

The Bicoid morphogen directs pattern formation along the anterior-posterior (A-P) axis of the Drosophila embryo. Bicoid is distributed in a concentration gradient that decreases exponentially from the anterior pole, however, it transcribes target genes such as hunchback in a step-function-like pattern; the expression domain is uniform and has a sharply defined posterior boundary. A 'gradient-affinity' model proposed to explain Bicoid action states that (i) cooperative gene activation by Bicoid generates the sharp on/off switch for target gene transcription and (ii) target genes with different affinities for Bicoid are expressed at different positions along the A-P axis. Using an in vivo yeast assay and in vitro methods, we show that Bicoid binds DNA with pairwise cooperativity; Bicoid bound to a strong site helps Bicoid bind to a weak site. These results support the first aspect of the model, providing a mechanism by which Bicoid generates sharp boundaries of gene expression. However, contrary to the second aspect of the model, we find no significant difference between the affinity of Bicoid for the anterior gene hunchback and the posterior gene knirps. We propose, instead, that the arrangement of Bicoids bound to the target gene presents a unique signature to the transcription machinery that, in combination with overall affinity, regulates the extent of gene transcription along the A-P axis.

Activation of posterior pair-rule stripe expression in response to maternal caudal and zygotic knirps activities

Drosophila pair-rule gene expression, in an array of seven evenly spaced stripes along the anterior-posterior axis of the blastoderm embryo, is controlled by distinct cis-acting stripe elements. In the anterior region, such elements mediate transcriptional activation in response to the maternal concentration gradient of the anterior determinant BICOID and repression by spatially distinct activities of zygotic gap genes. In the posterior region, activation of hairy stripe 6 has been shown to depend on the activity of the gap gene knirps, suggesting that posterior stripe expression is exclusively controlled by zygotic regulators. Here we show that the zygotic activation of hairy stripe 6 expression is preceded by activation in response to maternal caudal activity. Thus, transcriptional activation of posterior stripe expression is likely to be controlled by maternal and zygotic factors as has been observed for anterior stripes. The results suggest that activation and the expression level mediated by the hairy stripe 6-element depend on the number of activator binding sites, likely to involve additive rather than synergistic interactions. We found an identical transacting factor requirement for hairy stripe 6 and 7 expression. The arrangement of the corresponding binding sites for the common factors involved in the control of the two stripes share a high degree of similarity, but some of the factors exert opposite regulatory functions within the two enhancer elements.

Regulation of Drosophila spalt gene expression

The region-specific homeotic gene spalt is involved in the specification of terminal versus trunk structures during early Drosophila embryogenesis. Later in development spalt activity participates in specific processes during organogenesis and larval imaginal disc development. The multiple functions of spalt are reflected in distinct spatio-temporal expression patterns throughout development. Here we show that spalt cis-regulatory sequences for region-specific and organ-specific expression are clustered. Their organization may provide the structural basis for the diversification of expression pattern within the spalt/spalt related/spalt adjacent gene complex. We also examined the transacting factor requirement for the blastodermal spalt expression domains. They are under the genetic control of maternal and gap gene products and we show that these products are able to bind to corresponding spalt cis-acting sequences in vitro. The results suggest that the transacting factors, as defined by genetic studies, functionally interact with the spalt regulatory region. In addition, we provide evidence that a zygotic gene product of the terminal system, Tailless, cooperates with the maternal gene product Caudal and thereby activates gene expression in the terminal region of the embryo.