Novel interactions between vertebrate Hox genes

Understanding why metazoan Hox/HOM-C genes are expressed in spatiotemporal sequences showing colinearity with their genomic sequence is a central challenge in developmental biology. Here, we studied the consequences of ectopically expressing Hox genes to investigate whether Hox-Hox interactions might help to order gene expression during very early vertebrate embryogenesis. Our study revealed conserved autoregulatory loops for the Hox4 and Hox7 paralogue groups, detected following ectopic expression Hoxb-4 or HOXD4, and Hoxa-7, respectively. We also detected specific induction of 5' posterior Hox genes; Hoxb-5 to Hoxb-9, following ectopic expression of Hoxb-4/HOXD4; Hoxb-8 and Hoxb-9 following ectopic expression of Hoxa-7. Additionally, we observed specific repression of 3' anterior genes, following ectopic expression of Hox4 and Hox7 paralogues. We found that induction of Hoxb-4 and Hoxb-5 by Hoxb-4 can be direct, whereas induction of Hoxb-7 is indirect, suggesting the possibility of an activating cascade. Finally, we found that activation of Hoxb-4 itself and of posterior Hox genes by Hoxb-4 can be both non-cell-autonomous, as well as direct. We believe that our findings could be important for understanding how a highly ordered Hox expression sequence is set up in the early vertebrate embryo.     

Nucleotides flanking a conserved TAAT core dictate the DNA binding specificity of three murine homeodomain proteins.

Murine homeobox genes play a fundamental role in directing embryogenesis by controlling gene expression during development. The homeobox encodes a DNA binding domain (the homeodomain) which presumably mediates interactions of homeodomain proteins with specific DNA sites in the control regions of target genes. However, the bases for these selective DNA-protein interactions are not well defined. In this report, we have characterized the DNA binding specificities of three murine homeodomain proteins, Hox 7.1, Hox 1.5, and En-1. We have identified optimal DNA binding sites for each of these proteins by using a random oligonucleotide selection strategy. Comparison of the sequences of the selected binding sites predicted a common consensus site that contained the motif (C/G)TAATTG. The TAAT core was essential for DNA binding activity, and the nucleotides flanking this core directed binding specificity. Whereas variations in the nucleotides flanking the 5' side of the TAAT core produced modest alterations in binding activity for all three proteins, perturbations of the nucleotides directly 3' of the core distinguished the binding specificity of Hox 1.5 from those of Hox 7.1 and En-1. These differences in binding activity reflected differences in the dissociation rates rather than the equilibrium constants of the protein-DNA complexes. Differences in DNA binding specificities observed in vitro may contribute to selective interactions of homeodomain proteins with potential binding sites in the control regions of target genes.     

The optimal binding sequence of the Hox11 protein contains a predicted recognition core motif.

HOX11 is a homeobox-containing oncogene of specific T-cell leukemias. We determined the DNA binding specificity of the Hox11 protein by using a novel technique of random oligonucleotide selection developed in this study. The optimal Hox11 binding sequence, GGCGGTAAGTGG, contained a core TAAGTG motif that is consistent with a prediction based on the residues at specific positions that potentially make DNA base contacts and models of homeodomain-DNA interaction proposed from studies with other homeodomains. The specific interaction between Hox11 and the selected optimal binding sequence was further confirmed by band-shift and DNA competition assays. Given that the Hox11 homeodomain shares low homology with other well studied homeodomains, the presence of a predictable recognition core motif in its optimal binding sequence supports the notion that different homeodomains interact with DNA in a similar manner, through highly conserved residues at specific positions that allow contact with DNA.     

Hox homeodomain proteins exhibit selective complex stabilities with Pbx and DNA.

Eight of the nine homeobox genes of the Hoxb locus encode proteins which contain a conserved hexapeptide motif upstream from the homeodomain. All eight proteins (Hoxb-1-Hoxb-8) bind to a target oligonucleotide in the presence of Pbx1a under conditions where minimal or no binding is detected for the Hox or Pbx1a proteins alone. The stabilities of the Hox-Pbx1a-DNA complexes vary >100-fold, with the proteins from the middle of the locus (Hoxb-5 and Hoxb-6) forming very stable complexes, while Hoxb-4, Hoxb-7 and Hoxb-8 form complexes of intermediate stability and proteins at the 3'-side of the locus (Hoxb-1-Hoxb-3) form complexes which are very unstable. Although Hox-b proteins containing longer linker sequences between the hexapeptide and homeodomains formed unstable complexes, shortening the linker did not confer complex stability. Homeodomain swapping experiments revealed that this motif does not independently determine complex stability. Naturally occurring variations within the hexapeptides of specific Hox proteins also do not explain complex stability differences. However, two core amino acids (tryptophan and methionine) which are absolutely conserved within the hexapeptide domains appear to be required for complex formation. Removal of N- and C-terminal flanking regions did not influence complex stability and the members of paralog group 4 (Hoxa-4, b-4, c-4 and d-4), which share highly conserved hexapeptides, linkers and homeodomains but different flanking regions, form complexes of similar stability. These data suggest that the structural features of Hox proteins which determine Hox-Pbx1a-DNA complex stability reside within the precise structural relationships between the homeodomain, hexapeptide and linker regions.     

Evidence that Hoxa expression domains are evolutionarily transposed in spinal ganglia, and are established by forward spreading in paraxial mesoderm.

Transposition of anatomical structures along the anteroposterior axis has been a commonly used mechanism for changing body proportions during the course of evolutionary time. Earlier work (Gaunt, S.J., 1994. Conservation in the Hox code during morphological evolution. Int. J. Dev. Biol. 38, 549-552; Burke, A.C., Nelson, C.E., Morgan, B.A., Tabin, C., 1995. Hox genes and the evolution of vertebrate axial morphology. Development 121, 333-346) showed how transposition in mesodermal derivatives (vertebrae) could be attributed to transposition in the expression of Hox genes along the axial series of somites. We now show how transposition in the segmental arrangement of the spinal nerves can also be correlated with shifts in the expression domains of Hox genes. Specifically, we show how the expression domains of Hoxa-7, a-9 and a-10 in spinal ganglia correspond similarly in both mouse and chick with the positions of the brachial and lumbosacral plexuses, and that this is true even though the brachial plexus of chick is shifted posteriorly, relative to mouse, by seven segmental units. In spite of these marked species differences in the boundaries of Hoxa-7 expression, cis regulatory elements located up to 5 kb upstream of the chick Hoxa-7 gene showed much functional and structural conservation with those described in the mouse (Puschel, A.W., Balling, R., Gruss, P., 1991. Separate elements cause lineage restriction and specify boundaries of Hox-1.1 expression. Development 112, 279-287; Knittel, T., Kessel, M., Kim, M.H., Gruss, P., 1995. A conserved enhancer of the human and murine Hoxa-7 gene specifies the anterior boundary of expression during embryonal development. Development 121, 1077-1088). We also show that chick Hoxa-7 and a-10 expression domains spread forward into regions of somites that are initially negative for the expression of these genes. We discuss this as evidence that Hox expression in paraxial mesoderm spreads forward, as earlier found for neurectoderm and lateral plate mesoderm, in a process that occurs independently of cell movement.     

Characterization of Hoxa-10/Hoxa-11 transheterozygotes reveals functional redundancy and regulatory interactions

Hox genes show related sequences and overlapping expression domains that often reflect functional redundancy as well as a common evolutionary origin. To accurately define their functions, it has become necessary to compare phenotypes of mice with single and multiple Hox gene mutations. Here, we focus on two Abd-B-type genes, Hoxa-10 and Hoxa-11, which are coexpressed in developing vertebrae, limbs, and reproductive tracts. To assess possible functional redundancy between these two genes, Hoxa-10/Hoxa-11 transheterozygotes were produced by genetic intercrosses and analyzed. This compound mutation resulted in synergistic defects in transheterozygous limbs and reproductive tracts, but not in vertebrae. In the forelimb, distal radial/ulnar thickening and pisiform/triangular carpal fusion were observed in 35 and 21% of transheterozygotes, respectively, but were effectively absent in Hoxa-10 and Hoxa-11 +/- forelimbs. In the hindlimb, distal tibial/fibular thickening and loss of tibial/fibular fusion were observed in >80% of transheterozygotes but in no Hoxa-10 or Hoxa-11 +/- hindlimbs, and all transheterozygotes displayed reduced medial patellar sesamoids, compared to modest incidences in Hoxa-10 and Hoxa-11 +/- mutants. Furthermore, while the reproductive tracts of Hoxa-10 and Hoxa-11 single heterozygous mutants of both sexes were primarily unaffected, male transheterozygotes displayed cryptorchidism and abnormal tortuosity of the ductus deferens, and female transheterozygotes exhibited abnormal uterotubal junctions and narrowing of the uterus. In addition we observed that the targeted mutations of Hoxa-10 and Hoxa-11 each affected the expression of the other gene in the developing prevertebra and reproductive tracts. These results provide a measure of the functional redundancy of Hoxa-10 and Hoxa-11 and a deeper understanding of the phenotypes resulting in the single mutants and help elucidate the regulatory interactions between these two genes.     

Primary structure, developmentally regulated expression and potential duplication of the zebrafish homeobox gene ZF-21.

We report the molecular cloning and characterization of a cDNA derived from a zebrafish gene (ZF-21) related to the mouse homeobox containing gene Hox2.1. Interesting information about the differential conservation of various domains was gained from comparisons between the putative protein sequences from ZF-21 (275 amino acids) and Hox2.1 (279 aa). A separate DNA binding domain including the ZF-21 homeodomain and 36 additional flanking residues is completely identical to the C-terminal part of Hox2.1. As a consequence, these two mouse and zebrafish proteins must have identical DNA binding properties. A lower level of sequence identity between the N-terminal coding regions of ZF-21 and Hox2.1 reduces the total protein homology to 81%. However, short stretches of perfect homology in these N-terminals suggests that the essential biochemical functions are the same. As expected for true homologues, the ZF-21 and Hox2.1 genes also share extensive similarities with respect to non-coding sequences and temporal expression during embryogenesis. The finding of a potential ZF-21 duplication is discussed in relation to functional and evolutionary aspects of vertebrate homeobox genes.     

Functional characterization of the Escherichia coli Fis-DNA binding sequence

The Escherichia coli protein Fis is remarkable for its ability to interact specifically with DNA sites of highly variable sequences. The mechanism of this sequence-flexible DNA recognition is not well understood. In a previous study, we examined the contributions of Fis residues to high-affinity binding at different DNA sequences using alanine-scanning mutagenesis and identified several key residues for Fis-DNA recognition. In this work, we investigated the contributions of the 15-bp core Fis binding sequence and its flanking regions to Fis-DNA interactions. Systematic base-pair replacements made in both half sites of a palindromic Fis binding sequence were examined for their effects on the relative Fis binding affinity. Missing contact assays were also used to examine the effects of base removal within the core binding site and its flanking regions on the Fis-DNA binding affinity. The results revealed that: (1) the − 7G and + 3Y bases in both DNA strands (relative to the central position of the core binding site) are major determinants for high-affinity binding; (2) the C⁵ methyl group of thymine, when present at the + 4 position, strongly hinders Fis binding; and (3) AT-rich sequences in the central and flanking DNA regions facilitate Fis-DNA interactions by altering the DNA structure and by increasing the local DNA flexibility. We infer that the degeneracy of specific Fis binding sites results from the numerous base-pair combinations that are possible at noncritical DNA positions (from − 6 to − 4, from − 2 to + 2, and from + 4 to + 6), with only moderate penalties on the binding affinity, the roughly similar contributions of − 3A or G and + 3T or C to the binding affinity, and the minimal requirement of three of the four critical base pairs to achieve considerably high binding affinities.     

矛尾鱼Hoxa—11基因克隆及5′端序列分析

Hoxa-11基因调节鱼类鳍和四足动物肢的发育,在脊椎动物进化过程中起着重要的作用,利用人和鼠的Hoxa-11基因保守序列设计了两个兼并引物,通过PCR扩增到了矛尾鱼的Hoxa-11基因,经克隆和DNA序列分析,该片段为2065bp,包括绝大部分外显子Ⅰ,内含子和部分外显子Ⅱ,编码204个氨基酸,其氨基酸序列与人、鼠、鸡、蛙和斑马鱼的同源性分别为66.0%、67.6%、74.4%、72.8%和59.7%。外显子Ⅰ的长度从矛尾鱼到蛙、鸡、鼠和人呈现逐步上升趋势,人比矛尾鱼增长了16%,进一步分析,外显子Ⅰ可分为4个区域;两个高度保守区域,1个中度保守区域和1个可变区域,外显子Ⅰ的长度变化主要是由于可变区域内丙氨酸同聚物以及两侧富含甘氨酸和丝氨酸序列的累积。矛尾鱼只有1个由两个丙氨酸组成的同类物,蛙有1个由5个连续丙氨酸组成的同聚物,而鸡、鼠和人有3个丙氨酸同聚物,其中最大的同聚物由7个连续丙氨酸组成,而且在同聚物两侧出现了富含甘氨酸和丝氨酸序列。这表明可变区域可能与脊椎动物进化和鳍-肢转换过程中新功能的获得有关。同源异型盒所在的外显子Ⅱ区和剪接位点是高度保守的。内含子的长度变化较大,但在其内部也发现了两个高度保守的35bp和16bp的DNA片段,这两个片段在人、鼠、鸡、蛙和矛尾鱼中是完全相同的,这些序列的高度保守性提示其功能上的重要性。     

Identification and sequence analysis of the 5' end of the major chicken vitellogenin gene.

We have precisely determined the positions of the first three exons for the major chicken vitellogenin gene (VTG II) by a combination of S1 nuclease protection, primer extension and DNA sequencing experiments. In addition, we have determined the nucleotide sequences of the 5' flanking nuclease hypersensitive sites that we have previously shown are induced during the estrogen mediated activation of the VTG II gene in liver (1). One of these sites is found to be nearly identical to the enhancer core sequence of SV40. A computer assisted analysis of the DNA sequences upstream from the VTG II gene has revealed four short (7 to 9 base pair) sequence elements that are present in similar positions flanking the other major estrogen inducible gene for liver, very low density apolipoprotein II (apoVLDL II). For VTG II, these sequences are located between two of the induced nuclease hypersensitive sites that are liver specific. Sequences homologous to one element, located approximately 100 base pairs upstream from the mRNA cap sites of the VTG II and apoVLDL II genes, are also observed for three estrogen inducible genes that are expressed in the oviduct, although for each of these genes the sequence falls further upstream, between -220 and -200. We suggest that these conserved sequences may be important in mediating the tissue specific responses of these genes to estrogen.     

The interaction with DNA of wild-type and mutant fushi tarazu homeodomains.

The in vitro DNA binding properties of wild-type and mutant fushi tarazu homeodomains (ftz HD) have been analysed. The DNA binding properties of the ftz HD are very similar to those of the Antp HD. In interference experiments with mutant ftz HDs, close approaches between specific portions of the ftz HD peptide and specific regions of the binding site DNA were mapped. A methylation interference, G7 on the beta strand of BS2, is absent from the interference pattern with a mutant ftz HD [ftz (R43A) HD] in which the Arg43 at the second position of helix III (the recognition helix) is replaced by an Ala. This indicated that Arg43 of the ftz HD is in close proximity to the N7 of G7 of the beta strand of BS2 in the major groove. The methylation and ethylation interference patterns with the ftz (NTD) HD, in which the first six amino acids of the homeodomain were deleted, were extensively altered relative to the ftz HD patterns. Methylation of A11 and G12 of the alpha strand and ethylation of the phosphate of nucleotide A12 of the alpha strand no longer interfere with binding. This indicated that the first six amino acids of the homeodomain of ftz interact with A11 of the alpha strand in the minor groove, the phosphate of the nucleotide A13 on the alpha strand and G12 of the alpha strand in the adjacent major groove of BS2. In a binding study using a change of specificity mutation [ftz (Q50K) HD], in which the Gln50 at the ninth position of the third helix is exchanged for a Lys (as in the bicoid HD), and variant binding sites, we concluded that position 50 of the ftz HD and the ftz (Q50K) HD peptides interacts with base pairs at positions 6 and 7 of BS2. These three points of contact allowed us to propose a crude orientation of the ftz HD within the protein-DNA complex. We find that the ftz HD and the Antp HD peptides contact DNA in a similar way.