Mutual repression enhances the steepness and precision of gene expression boundaries

Embryonic development is driven by spatial patterns of gene expression that determine the fate of each cell in the embryo. While gene expression is often highly erratic, embryonic development is usually exceedingly precise. In particular, gene expression boundaries are robust not only against intra-embryonic fluctuations such as noise in gene expression and protein diffusion, but also against embryo-to-embryo variations in the morphogen gradients, which provide positional information to the differentiating cells. How development is robust against intra- and inter-embryonic variations is not understood. A common motif in the gene regulation networks that control embryonic development is mutual repression between pairs of genes. To assess the role of mutual repression in the robust formation of gene expression patterns, we have performed large-scale stochastic simulations of a minimal model of two mutually repressing gap genes in Drosophila, hunchback (hb) and knirps (kni). Our model includes not only mutual repression between hb and kni, but also the stochastic and cooperative activation of hb by the anterior morphogen Bicoid (Bcd) and of kni by the posterior morphogen Caudal (Cad), as well as the diffusion of Hb and Kni between neighboring nuclei. Our analysis reveals that mutual repression can markedly increase the steepness and precision of the gap gene expression boundaries. In contrast to other mechanisms such as spatial averaging and cooperative gene activation, mutual repression thus allows for gene-expression boundaries that are both steep and precise. Moreover, mutual repression dramatically enhances their robustness against embryo-to-embryo variations in the morphogen levels. Finally, our simulations reveal that diffusion of the gap proteins plays a critical role not only in reducing the width of the gap gene expression boundaries via the mechanism of spatial averaging, but also in repairing patterning errors that could arise because of the bistability induced by mutual repression.     

The role of Bicoid cooperative binding in the patterning of sharp borders in Drosophila melanogaster

In Drosophila embryonic development, the Bicoid (Bcd) protein establishes positional information of downstream developmental genes like hunchback (hb), which has a strong anterior expression and a sharp on-off boundary in the mid-embryo. The role of Bcd cooperative binding in the positioning of the Hb pattern has been previously demonstrated. However, there are discrepancies in the reported results about the role of this mechanism in the sharp Hb border. Here, we determined the Hill coefficient (n(H)) required for Bcd to generate the sharp border of Hb in wild-type (WT) embryos. We found that an n(H) of approximately 6.3 (s.d. 1.4) and 10.8 (s.d. 4.0) is required to account for Hb sharpness at early and late cycle 14A, respectively. Additional mechanisms are possibly required because the high n(H) is likely unachievable for Bcd binding to the hb promoter. To test this idea, we determined the n(H) required to pattern the Hb profile of 15 embryos expressing an hb(14F) allele that is defective in self-activation and found n(H) to be 3.0 (s.d. 1.0). This result indicates that in WT embryos, the hb self-activation is important for Hb sharpness. Corroborating our results, we also found a progressive increase in the required value of n(H) spanning from 4.0 to 9.2 by determining this coefficient from averaged profiles of eight temporal classes at cycle 14A (T1 to T8). Our results indicate that there is a transition in the mechanisms responsible for the sharp Hb border during cycle 14A: in early stages of this cycle, Bcd cooperative binding is primarily responsible for Hb sharpness; in late cycle 14A, hb self-activation becomes the dominant mechanism.     

Differential affinity and cooperativity functions of the amino-terminal 70 residues of lambda integrase

The site-specific recombinase (Int) of bacteriophage lambda is a heterobivalent DNA-binding protein that binds two different classes of DNA-binding sites within its recombination target sites. The several functions of Int are apportioned between a large carboxy-terminal domain that cleaves and ligates DNA at each of its four "core-type" DNA-binding sites and a small amino-terminal domain, whose primary function is binding to each of its five "arm-type" DNA sites, which are distant from the core region. Int bridges between the two classes of binding sites are facilitated by accessory DNA-bending proteins that along with Int comprise higher-order recombinogenic complexes. We show here that although the 64 amino-terminal residues of Int bind efficiently to a single arm site, this protein cannot form doubly bound complexes on adjacent arm sites. However, 1-70 Int does show the same cooperative binding to adjacent arm sites as the full length protein. We also found that 1-70 Int specifies cooperative interactions with the accessory protein Xis when the two are bound to their adjacent cognate sites P2 and X1, respectively. To complement the finding that these two amino-terminal domain functions (along with arm DNA binding) are all specified by residues 1-70, we determined that Thr75 is the first residue of the minimal carboxy-terminal domain, thereby identifying a specific interdomain linker region. We have measured the affinity constants for Int binding to each of the five arm sites and the cooperativity factors for Int binding to the two pairs of adjacent arm sites, and we have identified several DNA structural features that contribute to the observed patterns of Int binding to arm sites. Taken together, the results highlight several interesting features of arm DNA binding that invite speculation about additional levels of complexity in the regulation of lambda site-specific recombination.     

Mutations in the amino-terminal domain of lambda-integrase have differential effects on integrative and excisive recombination

Lambda integrase (Int) forms higher-order protein-DNA complexes necessary for site-specific recombination. The carboxy-terminal domain of Int (75-356) is responsible for catalysis at specific core-type binding sites whereas the amino-terminal domain (1-70) is responsible for cooperative arm-type DNA binding. Alanine scanning mutagenesis of residues 64-70, within full-length integrase, has revealed differential effects on cooperative arm binding interactions that are required for integrative and excisive recombination. Interestingly, while these residues are required for cooperative arm-type binding on both P'1,2 and P'2,3 substrates, cooperative binding at the arm-type sites P'2,3 was more severely compromised than binding at arm-type sites P'1,2 for L64A. Concomitantly, L64A had a much stronger effect on integrative than on excisive recombination. The arm-binding properties of Int appear to be intrinsic to the amino-terminal domain because the phenotype of L64A was the same in an amino-terminal fragment (Int 1-75) as it was in the full-length protein.     

Interaction of nuclear factor EF-1A with the polyomavirus enhancer region.

Enhancer factor 1A (EF-1A) is a mammalian nuclear protein that previously was shown to bind cooperatively to the repeated core enhancer element I sequence in the adenovirus E1A enhancer region. We now have characterized three binding sites for EF-1A in the polyomavirus A2 (Py) enhancer region. Site 1 resides in the Py A enhancer domain, and sites 2 and 3 reside in the Py B enhancer domain. EF-1A binding to Py site 1 is independent of cooperation with other EF-1A sites or the adjacent binding sites for PEA-1 and PEA-2, two murine nuclear factors that bind in the Py A enhancer domain. EF-1A binding to Py sites 2 and 3, in contrast, is cooperative, similar to the situation previously observed with binding sites in the adenovirus E1A enhancer region. In a transient replication assay, EF-1A site 1 functions synergistically with the PEA-1 and PEA-2 sites in the A enhancer domain to enhance Py replication. The functional cooperativity observed with the EF-1A, PEA-1, and PEA-2 sites in vivo does not reflect cooperative DNA binding interactions, as detected in vitro. Py EF-1A site 1 alone is capable of weakly stimulating Py replication. EF-1A site 1 overlaps with the binding sites for the murine nuclear protein PEA-3 and the ets family of oncoproteins.     

Molecular mechanisms of pattern formation by the BRE enhancer of the Ubx gene.

The core activity of the Ubx gene enhancer BRE (bx region enhancer) is encoded within a 500 bp module. bx DNA outside this active module increases the level of expression, expands the expression into ventro-lateral ectoderm and partially stabilizes the late expression pattern. The products of the gap genes hb and tll and of the pair-rule gene ftz bind to the 500 bp BRE module and control directly its initial pattern of expression. ftz enhances expression in even-numbered parasegments within the correct spatial domain whose boundaries are set by hb and tll. In addition, en and twi products activate the enhancer, probably directly. en broadens the parasegmental stripe while twi cooperates with ftz to enhance expression in the mesoderm. Binding sites for the five regulators are closely clustered, often overlapping extensively with one another. In vitro, hb blocks the binding of ftz and can also displace ftz protein pre-bound to an overlapping site, suggesting that competitive binding and/or interference by hb sets the initial boundaries of the domain of expression. Our results also suggest that this interaction is short-range and the long distance interactions among different enhancers may depend on each enhancer's ability to complex with the promoter.     

The bx region enhancer, a distant cis-control element of the Drosophila Ubx gene and its regulation by hunchback and other segmentation genes.

The Drosophila homeotic gene Ultrabithorax (Ubx) is regulated by complex mechanisms that specify the spatial domain, the timing and the activity of the gene in individual tissues and in individual cells. In early embryonic development, Ubx expression is controlled by segmentation genes turned on earlier in the developmental hierarchy. Correct Ubx expression depends on multiple regulatory sequences located outside the basal promoter. Here we report that a 500 bp DNA fragment from the bx region of the Ubx unit, approximately 30 kb away from the promoter, contains one of the distant regulatory elements (bx region enhancer, BRE). During early embryogenesis, this enhancer element activates the Ubx promoter in parasegments (PS) 6, 8, 10, and 12 and represses it in the anterior half of the embryo. The repressor of the anterior Ubx expression is the gap gene hunchback (hb). We show that the hb protein binds to the BRE element and that such binding is essential for hb repression in vivo, hb protein also binds to DNA fragments from abx and bxd, two other regulatory regions of the Ubx gene. We conclude that hb represses Ubx expression directly by binding to BRE and probably other Ubx regulatory elements. In addition, the BRE pattern requires input from other segmentation genes, among them tailless and fushi tarazu but not Krüppel and knirps.     

CTCF-dependent enhancer blockers at the upstream region of the chicken alpha-globin gene domain

The eukaryotic genome is partitioned into chromatin domains containing coding and intergenic regions. Insulators have been suggested to play a role in establishing and maintaining chromatin domains. Here we describe the identification and characterization of two separable enhancer blocking elements located in the 5′ flanking region of the chicken α-globin domain, 11–16 kb upstream of the embryonic α-type π gene in a DNA fragment harboring a MAR (matrix attachment region) element and three DNase I hypersensitive sites (HSs). The most upstream enhancer blocking element co-localizes with the MAR element and an erythroid-specific HS. The second enhancer blocking element roughly co-localizes with a constitutive HS. The third erythroid-specific HS present within the DNA fragment studied harbors a silencing, but not an enhancer blocking, activity. The 11 zinc-finger CCCTC-binding factor (CTCF), which plays an essential role in enhancer blocking activity in many previously characterized vertebrate insulators, is found to bind the two α-globin enhancer blocking elements. Detailed analysis has demonstrated that mutation of the CTCF binding site within the most upstream enhancer blocking element abolishes the enhancer blocking activity. The results are discussed with respect to special features of the tissue-specific α-globin gene domain located in a permanently open chromatin area.