Analysis of the murine Hoxa-9 cDNA: an alternatively spliced transcript encodes a truncated protein lacking the homeodomain

HOXA9 transcription factor is a double-edged sword: from development to cancer progression

The HOXA9 transcription factor serves as a molecular orchestrator in cancer stemness, epithelial-mesenchymal transition (EMT), metastasis, and generation of the tumor microenvironment in hematological and solid malignancies. However, the multiple modes of regulation, multifaceted functions, and context-dependent interactions responsible for the dual role of HOXA9 as an oncogene or tumor suppressor in cancer remain obscure. Hence, unravelling its molecular complexities, binding partners, and interacting signaling molecules enables us to comprehend HOXA9-mediated transcriptional programs and molecular crosstalk. However, it is imperative to understand its central role in fundamental biological processes such as embryogenesis, foetus implantation, hematopoiesis, endothelial cell proliferation, and tissue homeostasis before designing targeted therapies. Indeed, it presents an enormous challenge for clinicians to selectively target its oncogenic functions or restore tumor-suppressive role without altering normal cellular functions. In addition to its implications in cancer, the present review also focuses on the clinical applications of HOXA9 in recurrence and drug resistance, which may provide a broader understanding beyond oncology, open new avenues for clinicians for accurate diagnoses, and develop personalized treatment strategies. Furthermore, we have also discussed the existing therapeutic options and accompanying challenges in HOXA9-targeted therapies in different cancer types.

False-positive IRESes from Hoxa9 and other genes resulting from errors in mammalian 5' UTR annotations

Hyperconserved genomic sequences have great promise for understanding core biological processes. It has been recently proposed that scores of hyperconserved 5' untranslated regions (UTRs), also known as transcript leaders (hTLs), encode internal ribosome entry sites (IRESes) that drive cap-independent translation, in part, via interactions with ribosome expansion segments. However, the direct functional significance of such interactions has not yet been definitively demonstrated. We provide evidence that the putative IRESes previously reported in Hox gene hTLs are rarely included in transcript leaders. Instead, these regions function independently as transcriptional promoters. In addition, we find the proposed RNA structure of the putative Hoxa9 IRES is not conserved. Instead, sequences previously shown to be essential for putative IRES activity encode a hyperconserved transcription factor binding site (E-box) that contributes to its promoter activity and is bound by several transcription factors, including USF1 and USF2. Similar E-box sequences enhance the promoter activities of other putative Hoxa gene IRESes. Moreover, we provide evidence that the vast majority of hTLs with putative IRES activity overlap transcriptional promoters, enhancers, and 3' splice sites that are most likely responsible for their reported IRES activities. These results argue strongly against recently reported widespread IRES-like activities from hTLs and contradict proposed interactions between ribosomal expansion segment ES9S and putative IRESes. Furthermore, our work underscores the importance of accurate transcript annotations, controls in bicistronic reporter assays, and the power of synthesizing publicly available data from multiple sources.

Evolutionarily conserved inhibitory uORFs sensitize Hox mRNA translation to start codon selection stringency

Translation start site selection in eukaryotes is influenced by context nucleotides flanking the AUG codon and by levels of the eukaryotic translation initiation factors eIF1 and eIF5. In a search of mammalian genes, we identified five homeobox (Hox) gene paralogs initiated by AUG codons in conserved suboptimal context as well as 13 Hox genes that contain evolutionarily conserved upstream open reading frames (uORFs) that initiate at AUG codons in poor sequence context. An analysis of published cap analysis of gene expression sequencing (CAGE-seq) data and generated CAGE-seq data for messenger RNAs (mRNAs) from mouse somites revealed that the 5' leaders of Hox mRNAs of interest contain conserved uORFs, are generally much shorter than reported, and lack previously proposed internal ribosome entry site elements. We show that the conserved uORFs inhibit Hox reporter expression and that altering the stringency of start codon selection by overexpressing eIF1 or eIF5 modulates the expression of Hox reporters. We also show that modifying ribosome homeostasis by depleting a large ribosomal subunit protein or treating cells with sublethal concentrations of puromycin leads to lower stringency of start codon selection. Thus, altering global translation can confer gene-specific effects through altered start codon selection stringency.

Identification of a Modified HOXB9 mRNA in Breast Cancer

First identified as a developmental gene, HOXB9 is also known to be involved in tumor biological processes, and its aberrant expression correlates with poor prognosis of various cancers. In this study, we isolated a homeodomain-less, novel HOXB9 variant (HOXB9v) from human breast cancer cell line-derived mRNA. We confirmed that the novel variant was produced from variationless HOXB9 genomic DNA. RT-PCR of mRNA isolated from clinical samples and reanalysis of publicly available RNA-seq data proved that the new transcript is frequently expressed in human breast cancer. Exogenous HOXB9v expression significantly enhanced the proliferation of breast cancer cells, and gene ontology analysis indicated that apoptotic signaling was suppressed in these cells. Considering that HOXB9v lacks key domains of homeobox proteins, its behavior could be completely different from that of the previously described variationless HOXB9. Because none of the previous studies on HOXB9 have considered the presence of HOXB9v, further research analyzing the two transcripts individually is warranted to re-evaluate the true role of HOXB9 in cancer.

Mapping Whole-Transcriptome Splicing in Mouse Hematopoietic Stem Cells

Hematopoietic stem cells (HSCs) are rare cells that generate all the various types of blood and immune cells. High-quality transcriptome data have enabled the identification of significant genes for HSCs. However, most genes are expressed in various forms by alternative splicing (AS), extending transcriptome complexity. Here, we delineate AS to determine which isoforms are expressed in mouse HSCs. Our analysis of microarray and RNA-sequencing data includes differential expression of splicing factors that may regulate AS, and a complete map of splicing isoforms. Multiple types of isoforms for known HSC genes and unannotated splicing that may alter gene function are presented. Transcriptome-wide identification of genes and their respective isoforms in mouse HSCs will open another dimension for adult stem cells.

The leukemogenicity of Hoxa9 depends on alternative splicing

Although the transforming potential of Hox genes is known for a long time, it is not precisely understood to which extent splicing is important for the leukemogenicity of this gene family. To test this for Hoxa9, we compared the leukemogenic potential of the wild-type Hoxa9, which undergoes natural splicing, with a full-length Hoxa9 construct, which was engineered to prevent natural splicing (Hoxa9FLim). Inability to undergo splicing significantly reduced in vivo leukemogenicity compared to Hoxa9-wild-typed. Importantly, Hoxa9FLim could compensate for the reduced oncogenicity by collaborating with the natural splice variant Hoxa9T, as co-expression of Hoxa9T and Hoxa9FLim induced acute myeloid leukemia (AML) after a comparable latency time as wild-type Hoxa9. Hoxa9T on its own induced AML after a similar latency as Hoxa9FLim, despite its inability to bind DNA. These data assign splicing a central task in Hox gene mediated leukemogenesis and suggest an important role of homeodomain-less splice variants in hematological neoplasms.

The regulation of Hox gene expression during animal development

Hox genes encode a family of transcriptional regulators that elicit distinct developmental programmes along the head-to-tail axis of animals. The specific regional functions of individual Hox genes largely reflect their restricted expression patterns, the disruption of which can lead to developmental defects and disease. Here, we examine the spectrum of molecular mechanisms controlling Hox gene expression in model vertebrates and invertebrates and find that a diverse range of mechanisms, including nuclear dynamics, RNA processing, microRNA and translational regulation, all concur to control Hox gene outputs. We propose that this complex multi-tiered regulation might contribute to the robustness of Hox expression during development.

Two isoforms of HOXA9 function differently but work synergistically in human MLL-rearranged leukemia

HOXA9 plays a critical role in both normal hematopoiesis and leukemogenesis, particularly in the development and maintenance of mixed lineage leukemia (MLL)-rearranged leukemia. Through reverse transcription-polymerase chain reaction (RT-PCR) analysis of HOXA9 transcripts in human leukemia and normal bone marrow samples, we identified a truncated isoform of HOXA9, namely HOXA9T, and found that both HOXA9T and canonical HOXA9 were highly expressed in leukemia cell lines bearing MLL rearrangements, relative to human normal bone marrow cells or other subtypes of leukemia cells. A frameshift in HOXA9T in exon I causes a premature stop codon upstream of the PBX-binding domain and the homeodomain, which leads to the generation of a non-homeodomain-containing protein. Unlike the canonical HOXA9, HOXA9T alone cannot transform normal bone marrow progenitor cells. Moreover, HOXA9T cannot cooperate with MEIS1 to transform cells, despite the presence of a MEIS1-binding domain. Remarkably, although the truncated isoforms of many proteins function as dominant-negative competitors or inhibitors of their full-length counterparts, this is not the case for HOXA9T; instead, HOXA9T synergized with HOXA9 in transforming mouse normal bone marrow progenitor cells through promoting self-renewal and proliferation of the cells. Collectively, our data indicate that both truncated and full-length forms of HOXA9 are highly expressed in human MLL-rearranged leukemia, and the truncated isoform of HOXA9 might also play an oncogenic role by cooperating with canonical HOXA9 in cell transformation and leukemogenesis.

Integration of RNA processing and expression level control modulates the function of the Drosophila Hox gene Ultrabithorax during adult development

Although most metazoan genes undergo alternative splicing, the functional relevance of the majority of alternative splicing products is still unknown. Here we explore this problem in the Drosophila Hox gene Ultrabithorax (Ubx). Ubx produces a family of six protein isoforms through alternative splicing. To investigate the functional specificity of the Ubx isoforms, we studied their role during the formation of the Drosophila halteres, small dorsal appendages that are essential for normal flight. Our work shows that isoform Ia, which is encoded by all Ubx exons, is more efficient than isoform IVa, which lacks the amino acids coded by two small exons, in controlling haltere development and regulating Ubx downstream targets. However, our experiments also demonstrate that the functional differences among the Ubx isoforms can be compensated for by increasing the expression levels of the less efficient form. The analysis of the DNA-binding profiles of Ubx isoforms to a natural Ubx target, spalt, shows no major differences in isoform DNA-binding activities, suggesting that alternative splicing might primarily affect the regulatory capacity of the isoforms rather than their DNA-binding patterns. Our results suggest that to obtain distinct functional outputs during normal development genes must integrate the generation of qualitative differences by alternative splicing to quantitative processes affecting isoform protein expression levels.

Conservation of the TGFbeta/Labial homeobox signaling loop in endoderm-derived cells between Drosophila and mammals

BACKGROUND/AIMS

Midgut formation in Drosophila melanogaster is dependent upon the integrity of a signaling loop in the endoderm which requires the TGFbeta-related peptide, Decapentaplegic, and the Hox transcription factor, Labial. Interestingly, although Labial-like homeobox genes are present in mammals, their participation in endoderm morphogenesis is not clearly understood.

METHODS

We report the cloning, expression, localization, TGFbeta inducibility, and biochemical properties of the mammalian Labial-like homeobox, HoxA1, in exocrine pancreatic cells that are embryologically derived from the gut endoderm.

RESULTS

HoxA1 is expressed in pancreatic cell populations as two alternatively spliced messages, encoding proteins that share their N-terminal domain, but either lack or include the homeobox at the C-terminus. Transcriptional regulatory assays demonstrate that the shared N-terminal domain behaves as a strong transcriptional activator in exocrine pancreatic cells. HoxA1 is an early response gene for TGFbeta(1) in pancreatic epithelial cell populations and HoxA1 protein co-localizes with TGFbeta(1) receptors in the embryonic pancreatic epithelium at a time when exocrine pancreatic morphogenesis occurs (days E16 and E17).

CONCLUSIONS

These results report a role for HoxA1 in linking TGFbeta-mediated signaling to gene expression in pancreatic epithelial cell populations, thus suggesting a high degree of conservation for a TGFbeta/labial signaling loop in endoderm-derived cells between Drosophila and mammals. and IAP.

The truncated Hoxa1 protein interacts with Hoxa1 and Pbx1 in stem cells

Hox genes contain a homeobox encoding a 60-amino acid DNA binding sequence. The Hoxa1 gene (Hox1.6, ERA1) encodes two alternatively spliced mRNAs that encode distinct proteins, one with the homeodomain (Hoxa1-993), and another protein lacking this domain (Hoxa1-399). The functions of Hoxa1-399 are unknown. We detected Hoxa1-993 and Hoxa1-399 by immunoprecipitation using Hoxa1 antibodies. To assess whether Hoxa1-399 functions in cellular differentiation, we analyzed Hoxb1, a Hoxa1 target gene. Hoxa1-993 and its cofactor, Pbx1, bind to the Hoxb1 SOct-R3 promoter to transcriptionally activate a luciferase reporter. Results from F9 stem cells that stably express ectopic Hoxa1-399 (the F9-399 line) show that Hoxa1-399 reduces this transcriptional activation. Gel shift assays demonstrate that Hoxa1-399 reduces Hoxa1-993/Pbx1 binding to the Hoxb1 SOct-R3 region. GST pull-down experiments suggest that Hoxa1-399, Hoxa1-993, and Pbx1 form a trimer. However, the F9-399 line exhibits no differences in RA-induced proliferation arrest or endogenous Hoxb1, Pbx1, Hoxa5, Cyp26a1, GATA4, or Meis mRNA levels when compared to F9 wild-type.

MicroRNA-126 regulates HOXA9 by binding to the homeobox

The PicTar program predicted that microRNA-126 (miR-126), miR-145, and let-7s target highly conserved sites within the Hoxa9 homeobox. There are increased nucleotide constraints in the three microRNA seed sites among Hoxa9 genes beyond that required to maintain protein identity, suggesting additional functional conservation. In preliminary experiments, forced expression of these microRNAs in Hoxa9-immortalized bone marrow cells downregulated the HOXA9 protein and caused loss of biological activity. The microRNAs were shown to target their predicted sites within the homeobox. miR-126 and Hoxa9 mRNA are coexpressed in hematopoietic stem cells and downregulated in parallel during progenitor cell differentiation; however, miR-145 is barely detectable in hematopoietic cells, and let-7s are highly expressed in bone marrow progenitors, suggesting that miR-126 may function in normal hematopoietic cells to modulate HOXA9 protein. In support of this hypothesis, expression of miR-126 alone in MLL-ENL-immortalized bone marrow cells decreased endogenous HOXA9 protein, while inhibition of endogenous miR-126 increased expression of HOXA9 in F9 cells.

Transcriptional complexity of the HOXA9 locus

Evolutionarily conserved HOX genes play an important role during development and hematopoiesis. HOX protein products are transcription factors whose precise mechanism of action is still poorly understood. Regulation of HOX gene expression has been the topic of various studies. While alternative splicing and alternative promoter usage have been known to increase the number of transcripts across the HOX clusters, more recently high-throughput analyses have identified a number of new coding and noncoding RNA molecules whose function is not known. Here we review the transcriptome of the most studied HOX locus, HOXA9. Strict control of HOXA9 expression has been shown to play a critical role in hematopoiesis while aberrant expression has been shown to be important to the development of leukemia. However, it is still unclear how various transcripts from this locus are regulated and what specific role(s) each one of them plays.

Role of HOXA7 to HOXA13 and PBX1 genes in various forms of MRKH syndrome (congenital absence of uterus and vagina)

The Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome refers to the congenital absence or severe hypoplasia of the female genital tract, often described as uterovaginal aplasia which is the prime feature of the syndrome. It is the second cause of primary amenorrhea after gonadal dysgenesis and occurs in approximately 1 in 4500 women. Aetiology of this syndrome remains poorly understood. Frequent association of other malformations with the MRKH syndrome, involving kidneys, skeleton and ears, suggests the involvement of major developmental genes such as those of the HOX family. Indeed mammalian HOX genes are well known for their crucial role during embryogenesis, particularly in axial skeleton, hindbrain and limb development. More recently, their involvement in organogenesis has been demonstrated notably during urogenital differentiation. Although null mutations of HOX genes in animal models do not lead to MRKH-like phenotypes, dominant mutations in their coding sequences or aberrant expression due to mutated regulatory regions could well account for it. Sequence analysis of coding regions of HOX candidate genes and of PBX1, a likely HOX cofactor during Müllerian duct differentiation and kidney morphogenesis, did not reveal any mutation in patients showing various forms of MRKH syndrome. This tends to show that HOX genes are not involved in MRKH syndrome. However it does not exclude that other mechanisms leading to HOX dysfunction may account for the syndrome.

Leukemogenic MLL fusion proteins bind across a broad region of the Hox a9 locus, promoting transcription and multiple histone modifications

Chromosome translocations involving the mixed lineage leukemia gene MLL are associated with aggressive acute leukemias in both children and adults. Leukemogenic MLL fusion proteins delete the MLL SET domain Lys(4) methyltransferase activity and fuse MLL to 1 of >40 different translocation partners. Some MLL fusion proteins involve nuclear proteins that are transcriptional activators, whereas others have transcriptional activating activity but instead dimerize the truncated MLL molecule. Both types of MLL fusion proteins enforce persistent expression of Hox a9 and Meis1, which is pivotal for leukemogenesis through mechanisms that remain obscure. Here, we show that nuclear and dimerizable forms of MLL bind with a similar pattern to the Hox a9 locus that overlaps the distribution of wild-type MLL and deregulate transcription of three isoforms of Hox a9. Induction of MLL fusion protein activity is associated with increased levels of histone acetylation and Lys(4) methylation at Hox target genes. In addition, the MLL-ENL-ER protein, but not dimerized MLL, also induces dimethylation of histone H3 at Lys(79), suggesting alternative mechanisms for transcriptional activation.

MLL associates specifically with a subset of transcriptionally active target genes

MLL (mixed-lineage leukemia) is a histone H3 Lys-4 specific methyltransferase that is a positive regulator of Hox expression. MLL rearrangements and amplification are common in acute lymphoid and myeloid leukemias and myelodysplastic disorders and are associated with abnormal up-regulation of Hox gene expression. Although MLL is expressed throughout hematopoiesis, Hox gene expression is sharply down-regulated during differentiation, suggesting that either the activity of MLL or its association with target promoters must be regulated. Here we show that MLL associates with actively transcribed genes but does not remain bound after transcriptional down-regulation. Surprisingly, MLL is associated not only with promoter regions but also is distributed across the entire coding regions of genes. MLL interacts with RNA polymerase II (pol II) and colocalizes with RNA pol II at a subset of actively transcribed target in vivo. Loss of function Mll results in defects in RNA pol II distribution. Together the results suggest that an intimate association between MLL and RNA pol II occurs at MLL target genes in vivo that is required for normal initiation and/or transcriptional elongation.

Physical association and coordinate function of the H3 K4 methyltransferase MLL1 and the H4 K16 acetyltransferase MOF

A stable complex containing MLL1 and MOF has been immunoaffinity purified from a human cell line that stably expresses an epitope-tagged WDR5 subunit. Stable interactions between MLL1 and MOF were confirmed by reciprocal immunoprecipitation, cosedimentation, and cotransfection analyses, and interaction sites were mapped to MLL1 C-terminal and MOF zinc finger domains. The purified complex has a robust MLL1-mediated histone methyltransferase activity that can effect mono-, di-, and trimethylation of H3 K4 and a MOF-mediated histone acetyltransferase activity that is specific for H4 K16. Importantly, both activities are required for optimal transcription activation on a chromatin template in vitro and on an endogenous MLL1 target gene, Hox a9, in vivo. These results indicate an activator-based mechanism for joint MLL1 and MOF recruitment and targeted methylation and acetylation and provide a molecular explanation for the closely correlated distribution of H3 K4 methylation and H4 K16 acetylation on active genes.

Evolution of microRNAs located withinHox gene clusters

MicroRNAs (miRNAs) form an abundant class of non-coding RNA genes that have an important function in post-transcriptional gene regulation and in particular modulate the expression of developmentally important transcription factors including Hox genes. Two families of microRNAs are genomically located in intergenic regions in the Hox clusters of vertebrates. Here we describe their evolution in detail. We show that the micro RNAs closely follow the patterns of protein evolution in the Hox clusters, which is characterized by cluster duplications followed by differential gene loss.

A conserved non-homeodomain Hoxa9 isoform interacting with CBP is co-expressed with the 'typical' Hoxa9 protein during embryogenesis

Various Hox genes are known to produce alternative transcripts encoding different isoforms whose physiological relevance during development is not yet understood. In this work, we analysed two different Hoxa9 mRNAs encoding a full-length protein (Hoxa9) or a protein lacking the homeodomain (Hoxa9T). First, we demonstrated that these transcripts are conserved from birds to mammals. We then showed that both transcripts are present throughout embryogenesis and that Hoxa9T transcript is particularly abundant in embryonic genital tract, kidney, forelimb and tail. We further found that both isoforms are able to interact with CBP, suggesting a competition between Hoxa9 and Hoxa9T with this protein.

Disruption of a novel imprinted zinc-finger gene, ZNF215, in Beckwith-Wiedemann syndrome

The genetics of Beckwith-Wiedemann syndrome (BWS) is complex and is thought to involve multiple genes. It is known that three regions on chromosome 11p15 (BWSCR1, BWSCR2, and BWSCR3) may play a role in the development of BWS. BWSCR2 is defined by two BWS breakpoints. Here we describe the cloning and sequence analysis of 73 kb containing BWSCR2. Within this region, we detected a novel zinc-finger gene, ZNF215. We show that two of its five alternatively spliced transcripts are disrupted by both BWSCR2 breakpoints. Parts of the 3' end of these splice forms are transcribed from the antisense strand of a second zinc-finger gene, ZNF214. We show that ZNF215 is imprinted in a tissue-specific manner.

Lhx9 and Lhx9alpha LIM-homeodomain factors: genomic structure, expression patterns, chromosomal localization, and phylogenetic analysis

Lhx9 is a LIM-homeodomain (LIM-hd) transcription factor expressed in the embryonic mouse brain. We report the isolation of Lhx9alpha, a cDNA encoding a truncated isoform of Lhx9 that lacks the recognition helix of the homeodomain and differs from Lhx9 cDNA in its 3'-coding and 3'-UTR sequences. Isolation of the Lhx9 gene showed that Lhx9 and Lhx9alpha are coded by six exons spanning 10 kb of genomic sequence and that Lhx9alpha is an isoform generated by alternative splicing of the fifth exon. Lhx9 was mapped to the subtelomeric region of chromosome 1. Further molecular analysis showed that Lhx9 is a new candidate gene for the unidentified dreher (dr) mutation in mouse. The comparison of genomic structure and molecular phylogenetic analysis led to the identification of six groups of LIM-hd proteins, a basis for further classification and knowledge of their evolutionary relationships. To investigate a possible role for Lhx9alpha, the expression patterns of Lhx9 and Lhx9alpha were compared during embryogenesis. Lhx9alpha was expressed at lower levels than Lhx9, with a similar but distinct pattern in the brain, especially in the neocortex. We suggest that Lhx9alpha could function as an endogenous dominant-negative form of Lhx9 during development, both to regulate in space and time the transcriptional effects of Lhx9 and to add a degree of refinement to the LIM-hd code.

Endothelial cells express a novel, tumor necrosis factor-alpha-regulated variant of HOXA9

The expression of the class 1 homeobox (HOX) family of "master control" transcription factors has been studied principally in embryogenesis and neoplasia in which HOX genes play a critical role in cell proliferation, migration, and differentiation. We wished to test whether HOX family members were also involved in a differentiation-like process occurring in normal, diploid adult cells, that is, cytokine-induced activation of endothelial cells (EC). Screening of a human EC cDNA library yielded several members of the A and B groups of HOX transcription factors. One clone represented a novel, alternatively spliced variant of the human HOXA9 gene containing a new exon and the expression of which was driven by a novel promoter. This variant termed HOXA9EC appeared restricted to cells of endothelial lineage, i.e. expressed by human EC from multiple sources, but not by fibroblasts, smooth muscle cells, or several transformed cell lines. HOXA9EC mRNA was rapidly down-regulated in EC in response to tumor necrosis factor-alpha due to an apparent reduction in transcriptional rate. Reporter construct studies showed that the 400 base pairs of genomic DNA directly 5' to the transcription initiation site of HOXA9EC contained the information required for both up-regulation in response to cotransfection with a HOXA9EC expression vector and tumor necrosis factor-alpha-dependent down-regulation of this gene. These results provide evidence of a novel HOX family member that may participate in either the suppression or the genesis of EC activation.

Deletions in HOXD13 segregate with an identical, novel foot malformation in two unrelated families

Synpolydactyly (SPD) is a dominantly inherited congenital limb malformation consisting of 3/4 syndactyly in the hands and 4/5 syndactyly in the feet, with digit duplication in the syndactylous web. The condition recently has been found to result from different-sized expansions of an amino-terminal polyalanine tract in HOXD13. We report a novel type of mutation in HOXD13, associated in some cases with features of classic SPD and in all cases with a novel foot phenotype. In two unrelated families, each with a different intragenic deletion in HOXD13, all mutation carriers have a rudimentary extra digit between the first and second metatarsals and often between the fourth and fifth metatarsals as well. This phenotype has not been reported in any mice with genetic modifications of the HoxD gene cluster. The two different deletions affect the first exon and the homeobox, respectively, in each case producing frameshifts followed by a long stretch of novel sequence and a premature stop codon. Although the affected genes may encode proteins that exert a dominant negative or novel effect, they are most likely to act as null alleles. Either possibility has interesting implications for the role of HOXD13 in human autopod development.