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

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.

Direct and Indirect Targeting of HOXA9 Transcription Factor in Acute Myeloid Leukemia

HOXA9 (Homeobox A9) is a homeotic transcription factor known for more than two decades to be associated with leukemia. The expression of HOXA9 homeoprotein is associated with anterior-posterior patterning during embryonic development, and its expression is then abolished in most adult cells, with the exception of hematopoietic progenitor cells. The oncogenic function of HOXA9 was first assessed in human acute myeloid leukemia (AML), particularly in the mixed-phenotype associated lineage leukemia (MPAL) subtype. HOXA9 expression in AML is associated with aggressiveness and a poor prognosis. Since then, HOXA9 has been involved in other hematopoietic malignancies and an increasing number of solid tumors. Despite this, HOXA9 was for a long time not targeted to treat cancer, mainly since, as a transcription factor, it belongs to a class of protein long considered to be an "undruggable" target; however, things have now evolved. The aim of the present review is to focus on the different aspects of HOXA9 targeting that could be achieved through multiple ways: (1) indirectly, through the inhibition of its expression, a strategy acting principally at the epigenetic level; or (2) directly, through the inhibition of its transcription factor function by acting at either the protein/protein interaction or the protein/DNA interaction interfaces.

Role of HOXA9 in leukemia: dysregulation, cofactors and essential targets

HOXA9 is a homeodomain-containing transcription factor that has an important role in hematopoietic stem cell expansion and is commonly deregulated in acute leukemias. A variety of upstream genetic alterations in acute myeloid leukemia lead to overexpression of HOXA9, which is a strong predictor of poor prognosis. In many cases, HOXA9 has been shown to be necessary for maintaining leukemic transformation; however, the molecular mechanisms through which it promotes leukemogenesis remain elusive. Recent work has established that HOXA9 regulates downstream gene expression through binding at promoter distal enhancers along with a subset of cell-specific cofactor and collaborator proteins. Increasing efforts are being made to identify both the critical cofactors and target genes required for maintaining transformation in HOXA9-overexpressing leukemias. With continued advances in understanding HOXA9-mediated transformation, there is a wealth of opportunity for developing novel therapeutics that would be applicable for greater than 50% of AML with overexpression of HOXA9.

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.

Phylogenetic patterns of emergence of new genes support a model of frequent de novo evolution

Background

New gene emergence is so far assumed to be mostly driven by duplication and divergence of existing genes. The possibility that entirely new genes could emerge out of the non-coding genomic background was long thought to be almost negligible. With the increasing availability of fully sequenced genomes across broad scales of phylogeny, it has become possible to systematically study the origin of new genes over time and thus revisit this question.

Results

We have used phylostratigraphy to assess trends of gene evolution across successive phylogenetic phases, using mostly the well-annotated mouse genome as a reference. We find several significant general trends and confirm them for three other vertebrate genomes (humans, zebrafish and stickleback). Younger genes are shorter, both with respect to gene length, as well as to open reading frame length. They contain also fewer exons and have fewer recognizable domains. Average exon length, on the other hand, does not change much over time. Only the most recently evolved genes have longer exons and they are often associated with active promotor regions, i.e. are part of bidirectional promotors. We have also revisited the possibility that de novo evolution of genes could occur even within existing genes, by making use of an alternative reading frame (overprinting). We find several cases among the annotated Ensembl ORFs, where the new reading frame has emerged at a higher phylostratigraphic level than the original one. We discuss some of these overprinted genes, which include also the Hoxa9 gene where an alternative reading frame covering the homeobox has emerged within the lineage leading to rodents and primates (Euarchontoglires).

Conclusions

We suggest that the overall trends of gene emergence are more compatible with a de novo evolution model for orphan genes than a general duplication-divergence model. Hence de novo evolution of genes appears to have occurred continuously throughout evolutionary time and should therefore be considered as a general mechanism for the emergence of new gene functions.

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.

Multiple promoters and alternative splicing: Hoxa5 transcriptional complexity in the mouse embryo

Background

The genomic organization of Hox clusters is fundamental for the precise spatio-temporal regulation and the function of each Hox gene, and hence for correct embryo patterning. Multiple overlapping transcriptional units exist at the Hoxa5 locus reflecting the complexity of Hox clustering: a major form of 1.8 kb corresponding to the two characterized exons of the gene and polyadenylated RNA species of 5.0, 9.5 and 11.0 kb. This transcriptional intricacy raises the question of the involvement of the larger transcripts in Hox function and regulation.

Methodology/Principal Findings

We have undertaken the molecular characterization of the Hoxa5 larger transcripts. They initiate from two highly conserved distal promoters, one corresponding to the putative Hoxa6 promoter, and a second located nearby Hoxa7. Alternative splicing is also involved in the generation of the different transcripts. No functional polyadenylation sequence was found at the Hoxa6 locus and all larger transcripts use the polyadenylation site of the Hoxa5 gene. Some larger transcripts are potential Hoxa6/Hoxa5 bicistronic units. However, even though all transcripts could produce the genuine 270 a.a. HOXA5 protein, only the 1.8 kb form is translated into the protein, indicative of its essential role in Hoxa5 gene function. The Hoxa6 mutation disrupts the larger transcripts without major phenotypic impact on axial specification in their expression domain. However, Hoxa5-like skeletal anomalies are observed in Hoxa6 mutants and these defects can be explained by the loss of expression of the 1.8 kb transcript. Our data raise the possibility that the larger transcripts may be involved in Hoxa5 gene regulation.

Significance

Our observation that the Hoxa5 larger transcripts possess a developmentally-regulated expression combined to the increasing sum of data on the role of long noncoding RNAs in transcriptional regulation suggest that the Hoxa5 larger transcripts may participate in the control of Hox gene expression.

HOXA9 regulates BRCA1 expression to modulate human breast tumor phenotype

Breast cancer 1, early onset (BRCA1) expression is often reduced in sporadic breast tumors, even in the absence of BRCA1 genetic modifications, but the molecular basis for this is unknown. In this study, we identified homeobox A9 (HOXA9) as a gene frequently downregulated in human breast cancers and tumor cell lines and noted that reduced HOXA9 transcript levels associated with tumor aggression, metastasis, and patient mortality. Experiments revealed that loss of HOXA9 promoted mammary epithelial cell growth and survival and perturbed tissue morphogenesis. Restoring HOXA9 expression repressed growth and survival and inhibited the malignant phenotype of breast cancer cells in culture and in a xenograft mouse model. Molecular studies showed that HOXA9 restricted breast tumor behavior by directly modulating the expression of BRCA1. Indeed, ectopic expression of wild-type BRCA1 phenocopied the tumor suppressor function of HOXA9, and reducing BRCA1 levels or function inhibited the antitumor activity of HOXA9. Consistently, HOXA9 expression correlated with BRCA1 in clinical specimens and with tumor aggression in patients lacking estrogen receptor/progesterone receptor expression in their breast tissue. These findings indicate that HOXA9 restricts breast tumor aggression by modulating expression of the tumor suppressor gene BRCA1, which we believe provides an explanation for the loss of BRCA1 expression in sporadic breast tumors in the absence of BRCA1 genetic modifications.

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.

Fast and SNP-tolerant detection of complex variants and splicing in short reads

Motivation: Next-generation sequencing captures sequence differences in reads relative to a reference genome or transcriptome, including splicing events and complex variants involving multiple mismatches and long indels. We present computational methods for fast detection of complex variants and splicing in short reads, based on a successively constrained search process of merging and filtering position lists from a genomic index. Our methods are implemented in GSNAP (Genomic Short-read Nucleotide Alignment Program), which can align both single- and paired-end reads as short as 14 nt and of arbitrarily long length. It can detect short- and long-distance splicing, including interchromosomal splicing, in individual reads, using probabilistic models or a database of known splice sites. Our program also permits SNP-tolerant alignment to a reference space of all possible combinations of major and minor alleles, and can align reads from bisulfite-treated DNA for the study of methylation state.

Results: In comparison testing, GSNAP has speeds comparable to existing programs, especially in reads of ≥70 nt and is fastest in detecting complex variants with four or more mismatches or insertions of 1–9 nt and deletions of 1–30 nt. Although SNP tolerance does not increase alignment yield substantially, it affects alignment results in 7–8% of transcriptional reads, typically by revealing alternate genomic mappings for a read. Simulations of bisulfite-converted DNA show a decrease in identifying genomic positions uniquely in 6% of 36 nt reads and 3% of 70 nt reads.

Availability: Source code in C and utility programs in Perl are freely available for download as part of the GMAP package at http://share.gene.com/gmap.

Contact:twu@gene.com

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.