The homeobox genes are a large and diverse group of genes, many of which play important roles in the embryonic development of animals. Increasingly, homeobox genes are being compared between genomes in an attempt to understand the evolution of animal development. Despite their importance, the full diversity of human homeobox genes has not previously been described.

Although homeobox genes have been the subject of many studies, little is known about the main amino acid changes that occurred early in the evolution of genes belonging to different classes.

Homeobox genes constitute a large family of genes widely studied because of their role in the establishment of the body pattern. However, they are also involved in many other events during development and adulthood. The main olfactory epithelium (MOE) is an excellent model to study neurogenesis in the adult nervous system. Analyses of homeobox genes during development show that some of these genes are involved in the formation and establishment of cell diversity in the MOE. Moreover, the mechanisms of expression of odorant receptors (ORs) constitute one of the biggest enigmas in the field. Analyses of OR promoters revealed the presence of homeodomain binding sites in their sequences. Here we characterize the expression patterns of a set of 49 homeobox genes in the MOE with in situ hybridization. We found that seven of them (Dlx3, Dlx5, Dlx6, Msx1, Meis1, Isl1, and Pitx1) are zonally expressed. The homeobox gene Emx1 is expressed in three guanylate cyclase(+) populations, two located in the MOE and the third one in an olfactory subsystem known as Grüneberg ganglion located at the entrance of the nasal cavity. The homeobox gene Tshz1 is expressed in a unique patchy pattern across the MOE. Our findings provide new insights to guide functional studies that aim to understand the complexity of transcription factor expression and gene regulation in the MOE. J. Comp. Neurol. 524:2713-2739, 2016. © 2016 Wiley Periodicals, Inc.

Linkage studies in autism have identified susceptibility loci on chromosomes 2q and 7q, regions containing the DLX1/2 and DLX5/6 bigene clusters. The DLX genes encode homeodomain transcription factors that control craniofacial patterning and differentiation and survival of forebrain inhibitory neurons. We investigated the role that sequence variants in DLX genes play in autism by in-depth resequencing of these genes in 161 autism probands from the AGRE collection.

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

The full complement of homeobox transcription factor sequences, including genes and pseudogenes, was determined from the analysis of 10 complete genomes from flowering plants, moss, Selaginella, unicellular green algae, and red algae. Our exhaustive genome-wide searches resulted in the discovery in each class of a greater number of homeobox genes than previously reported. All homeobox genes can be unambiguously classified by sequence evolutionary analysis into 14 distinct classes also characterized by conserved intron-exon structure and by unique codomain architectures. We identified many new genes belonging to previously defined classes (HD-ZIP I to IV, BEL, KNOX, PLINC, WOX). Other newly identified genes allowed us to characterize PHD, DDT, NDX, and LD genes as members of four new evolutionary classes and to define two additional classes, which we named SAWADEE and PINTOX. Our comprehensive analysis allowed us to identify several newly characterized conserved motifs, including novel zinc finger motifs in SAWADEE and DDT. Members of the BEL and KNOX classes were found in Chlorobionta (green plants) and in Rhodophyta. We found representatives of the DDT, WOX, and PINTOX classes only in green plants, including unicellular green algae, moss, and vascular plants. All 14 homeobox gene classes were represented in flowering plants, Selaginella, and moss, suggesting that they had already differentiated in the last common ancestor of moss and vascular plants.

All bilaterian animals share a general genetic framework that controls the formation of their body structures, although their forms are highly diversified. The Hox genes that encode transcription factors play a central role in this framework. All Hox proteins contain a highly conserved homeodomain encoded by the homeobox motif, but the other regions are generally assumed to be less conserved. In this study, we used comparative genomic methods to infer possible functional elements in the coding regions of mammalian Hox genes.

Six genes are vertebrate homologues of the homeobox-containing gene sine oculis, which plays an essential role in controlling Drosophila compound eye development. Here we report the identification and expression patterns of all three subfamilies of Xenopus Six genes. Two Six2 subfamily genes (Six1, Six2) showed very similar expression patterns in cranial ganglia, otic placodes and the eyes. Non-neural expression of Six1 and Six2 was observed with mesodermal head mesenchyme, somites and their derivatives, the muscle anlagen of the embryonic trunk. In addition, Six2 expression was also found with mesenchyme associated with the developing stomach and pronephros. Expression of Six3 subfamily genes (Six3.1, Six3.2, Six6.1, and Six6.2) was restricted to the developing head, where expression was especially observed in derivatives of the forebrain (eyes, optic stalks, the hypothalamus and pituitary gland). Interestingly, expression of all Six3 subfamily members but Six6.2 was also found with the pineal gland primordium and the tegmentum. Expression of Six4 subfamily genes (Six4.1, Six4.2) was present in the developing visceral arches, placodal derivatives (otic vesicle, olfactory system), head mesenchyme and the eye. The observed dynamic expression patterns are largely conserved between lower and higher vertebrates and imply important roles of Six family genes not only in eye formation and myogenesis, but also in the development of the gut, the kidney and of placode-derived structures.

Recently, we have presented a scheme, termed "NKL-code", which describes physiological expression patterns of NKL homeobox genes in early hematopoiesis and in lymphopoiesis including main stages of T-, B- and NK-cell development. Aberrant activity of these genes underlies the generation of hematological malignancies notably T-cell leukemia. Here, we searched for deregulated NKL homeobox genes in main entities of T-cell lymphomas comprising angioimmunoblastic T-cell lymphoma (AITL), anaplastic large cell lymphoma (ALCL), adult T-cell leukemia/lymphoma (ATLL), hepatosplenic T-cell lymphoma (HSTL), NK/T-cell lymphoma (NKTL) and peripheral T-cell lymphoma (PTCL). Our data revealed altogether 19 aberrantly overexpressed genes in these types, demonstrating deregulated NKL homeobox genes involvement in T-cell lymphomas as well. For detailed analysis we focused on NKL homeobox gene MSX1 which is normally expressed in NK-cells. MSX1 was overexpressed in subsets of HSTL patients and HSTL-derived sister cell lines DERL-2 and DERL-7 which served as models to characterize mechanisms of deregulation. We performed karyotyping, genomic and expression profiling, and whole genome sequencing to reveal mutated and deregulated gene candidates, including the fusion gene CD53-PDGFRB. Subsequent knockdown experiments allowed the reconstruction of an aberrant network involved in MSX1 deregulation, including chromatin factors AUTS2 and mutated histone HIST1H3B(K27M). The gene encoding AUTS2 is located at chromosome 7q11 and may represent a basic target of the HSTL hallmark aberration i(7q). Taken together, our findings highlight an oncogenic role for deregulated NKL homeobox genes in T-cell lymphoma and identify MSX1 as a novel player in HSTL, implicated in aberrant NK- and T-cell differentiation.

Dlx5 and Dlx6 encode distal-less homeodomain transcription factors that are present in the genome as a linked pair at a single locus. Dlx5 and Dlx6 have redundant roles in craniofacial, skeletal, and uterine development. Previously, we performed a transcriptome comparison for anti-Müllerian hormone (AMH)-induced genes expressed in the Müllerian duct mesenchyme of male and female mouse embryos. In that study, we found that Dlx5 transcripts were nearly seven-fold higher in males compared to females and Dlx6 transcripts were found only in males, suggesting they may be AMH-induced genes. Therefore, we investigated the role of Dlx5 and Dlx6 during AMH-induced Müllerian duct regression. We found that Dlx5 was detected in the male Müllerian duct mesenchyme from E14.5 to E16.5. In contrast, in female embryos Dlx5 was detected in the Müllerian duct epithelium. Dlx6 expression in Müllerian duct mesenchyme was restricted to males. Dlx6 expression was not detected in female Müllerian duct mesenchyme or epithelium. Genetic experiments showed that AMH signaling is necessary for Dlx5 and Dlx6 expression. Müllerian duct regression was variable in Dlx5 homozygous mutant males at E16.5, ranging from regression like controls to a block in Müllerian duct regression. In E16.5 Dlx6 homozygous mutants, Müllerian duct tissue persisted primarily in the region adjacent to the testes. In Dlx5-6 double homozygous mutant males Müllerian duct regression was also found to be incomplete but more severe than either single mutant. These studies suggest that Dlx5 and Dlx6 act redundantly to mediate AMH-induced Müllerian duct regression during male differentiation.

We performed a genome-wide brain expression study to reveal the underpinnings of diseases linked to a repeat expansion in chromosome 9 open reading frame 72 (C9ORF72).

The homeobox encodes a DNA-binding domain found in transcription factors regulating key developmental processes. The most notable examples of homeobox containing genes are the Hox genes, arranged on chromosomes in the same order as their expression domains along the body axis. The mechanisms responsible for the synchronous regulation of Hox genes and the molecular function of their colinearity remain unknown. Here we report the discovery of a conserved structural signature of the 180-base pair DNA fragment comprising the homeobox. We demonstrate that the homeobox DNA has a characteristic 3-base-pair periodicity in the hydroxyl radical cleavage pattern. This periodic pattern is significant in most of the 39 mammalian Hox genes and in other homeobox-containing transcription factors. The signature is present in segmented bilaterian animals as evolutionarily distant as humans and flies. It remains conserved despite the fact that it would be disrupted by synonymous mutations, which raises the possibility of evolutionary selective pressure acting on the structure of the coding DNA. The homeobox coding DNA may therefore have a secondary function, possibly as a regulatory element. The existence of such element may have important consequences for understanding how these genes are regulated.

Lin-11, Isl-1 and Mec-3 domains (LIM) homeobox genes are among the most important sub-families of homeobox genes. These genes are thought to play an important role in cancer. In this study, the protein expression of these genes was examined in urothelial carcinoma of the bladder. The expression pattern of Islet-1 (ISL1) and LIM homeobox 5 (LHX5) across different cancer stages and grades, as well as the association between the protein expression of these genes and patient demographics and clinicopathological features, were examined.

We present a loss-of-function study using antisense morpholino (MO) reagents for the organizer-specific gene Goosecoid (Gsc) and the ventral genes Vent1 and Vent2. Unlike in the mouse Gsc is required in Xenopus for mesodermal patterning during gastrulation, causing phenotypes ranging from reduction of head structures-including cyclopia and holoprosencephaly-to expansion of ventral tissues in MO-injected embryos. The overexpression effects of Gsc mRNA require the expression of the BMP antagonist Chordin, a downstream target of Gsc. Combined Vent1 and Vent2 MOs strongly dorsalized the embryo. Unexpectedly, simultaneous depletion of all three genes led to a rescue of almost normal development in a variety of embryological assays. Thus, the phenotypic effects of depleting Gsc or Vent1/2 are caused by the transcriptional upregulation of their opposing counterparts. A principal function of Gsc and Vent1/2 homeobox genes might be to mediate a self-adjusting mechanism that restores the basic body plan when deviations from the norm occur, rather than generating individual cell types. The results may shed light on the molecular mechanisms of genetic redundancy.

In Hodgkin lymphoma (HL) we recently reported that deregulated homeobox gene MSX1 mediates repression of the B-cell specific transcription factor ZHX2. In this study we investigated regulation of MSX1 in this B-cell malignancy. Accordingly, we analyzed expression and function of OTX homeobox genes which activate MSX1 transcription during embryonal development in the neural plate border region. Our data demonstrate that OTX1 and OTX2 are aberrantly expressed in both HL patients and cell lines. Moreover, both OTX loci are targeted by genomic gains in overexpressing cell lines. Comparative expression profiling and subsequent pathway modulations in HL cell lines indicated that aberrantly enhanced FGF2-signalling activates the expression of OTX2. Downstream analyses of OTX2 demonstrated transcriptional activation of genes encoding transcription factors MSX1, FOXC1 and ZHX1. Interestingly, examination of the physiological expression profile of ZHX1 in normal hematopoietic cells revealed elevated levels in T-cells and reduced expression in B-cells, indicating a discriminatory role in lymphopoiesis. Furthermore, two OTX-negative HL cell lines overexpressed ZHX1 in correlation with genomic amplification of its locus at chromosomal band 8q24, supporting the oncogenic potential of this gene in HL. Taken together, our data demonstrate that deregulated homeobox genes MSX1 and OTX2 respectively impact transcriptional inhibition of (B-cell specific) ZHX2 and activation of (T-cell specific) ZHX1. Thus, we show how reactivation of a specific embryonal gene regulatory network promotes disturbed B-cell differentiation in HL.

The EMX (Empty Spiracles Homeobox) genes EMX1 and EMX2 are two homeodomain gene members of the EMX family of transcription factors involved in the regulation of various biological processes, such as cell proliferation, migration, and differentiation, during brain development and neural crest migration. They play a role in the specification of positional identity, the proliferation of neural stem cells, and the differentiation of certain neuronal cell phenotypes. In general, they act as transcription factors in early embryogenesis and neuroembryogenesis from metazoans to higher vertebrates. The EMX1 and EMX2's potential as tumor suppressor genes has been suggested in some cancers. Our work showed that EMX1/EMX2 act as tumor suppressors in sarcomas by repressing the activity of stem cell regulatory genes (OCT4, SOX2, KLF4, MYC, NANOG, NES, and PROM1). EMX protein downregulation, therefore, induced the malignance and stemness of cells both in vitro and in vivo. In murine knockout (KO) models lacking Emx genes, 3MC-induced sarcomas were more aggressive and infiltrative, had a greater capacity for tumor self-renewal, and had higher stem cell gene expression and nestin expression than those in wild-type models. These results showing that EMX genes acted as stemness regulators were reproduced in different subtypes of sarcoma. Therefore, it is possible that the EMX genes could have a generalized behavior regulating proliferation of neural crest-derived progenitors. Together, these results indicate that the EMX1 and EMX2 genes negatively regulate these tumor-altering populations or cancer stem cells, acting as tumor suppressors in sarcoma.

PAIRED (PRD)-like homeobox genes belong to a class of predicted transcription factor genes. Several of these PRD-like homeobox genes have been predicted in silico from genomic sequence but until recently had no evidence of transcript expression. We found recently that nine PRD-like homeobox genes, ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB, NOBOX, TPRX1 and TPRX2, were expressed in human preimplantation embryos. In the current study we characterized these PRD-like homeobox genes in depth and studied their functions as transcription factors. We cloned multiple transcript variants from human embryos and showed that the expression of these genes is specific to embryos and pluripotent stem cells. Overexpression of the genes in human embryonic stem cells confirmed their roles as transcription factors as either activators (CPHX1, CPHX2, ARGFX) or repressors (DPRX, DUXA, TPRX2) with distinct targets that could be explained by the amino acid sequence in homeodomain. Some PRD-like homeodomain transcription factors had high concordance of target genes and showed enrichment for both developmentally important gene sets and a 36 bp DNA recognition motif implicated in Embryo Genome Activation (EGA). Our data implicate a role for these previously uncharacterized PRD-like homeodomain proteins in the regulation of human embryo genome activation and preimplantation embryo development.

Homeodomain proteins are encoded by homeobox genes and regulate development and differentiation in many neuronal systems. The mouse vomeronasal organ (VNO) generates in situ mature chemosensory neurons from stem cells. The roles of homeodomain proteins in neuronal differentiation in the VNO are poorly understood. Here we have characterized the expression patterns of 28 homeobox genes in the VNO of C57BL/6 mice at postnatal stages using multicolor fluorescent in situ hybridization. We identified 11 homeobox genes (Dlx3, Dlx4, Emx2, Lhx2, Meis1, Pbx3, Pknox2, Pou6f1, Tshz2, Zhx1, Zhx3) that were expressed exclusively in neurons; 4 homeobox genes (Pax6, Six1, Tgif1, Zfhx3) that were expressed in all non-neuronal cell populations, with Pax6, Six1 and Tgif1 also expressed in some neuronal progenitors and precursors; 12 homeobox genes (Adnp, Cux1, Dlx5, Dlx6, Meis2, Pbx2, Pknox1, Pou2f1, Satb1, Tshz1, Tshz3, Zhx2) with expression in both neuronal and non-neuronal cell populations; and one homeobox gene (Hopx) that was exclusively expressed in the non-sensory epithelium. We studied further in detail the expression of Emx2, Lhx2, Meis1, and Meis2. We found that expression of Emx2 and Lhx2 initiated between neuronal progenitor and neuronal precursor stages. As far as the sensory neurons of the VNO are concerned, Meis1 and Meis2 were only expressed in the apical layer, together with Gnai2, but not in the basal layer.

Glioblastoma multiforme (GBM) exhibits considerable heterogeneity and associates with genome-wide alterations of the repressed chromatin marks DNA methylation and H3 lysine 27 trimethylation (H3K27me3). Tri-methylation on lysine 4 of histone H3 (H3K4me3) is an activating epigenetic mark that is enriched at promoter and promotes expression. It will be helpful in GBM diagnosis and treatment to identify the alteration of H3K4me3 between human GBM and GBM-surrounding tissues. Here, we performed an analysis using next-generation sequencing techniques to identify H3K4me3 modification in a case of GBM and the GBM-surrounding tissues. The results revealed a global decrease in H3K4me3 in GBM, especially at promoters and CpG islands. In GBM, homeobox genes gain H3K4me3, whereas the cell-cell adhesion-related cadherin genes lose H3K4me3. The products of the homeobox genes are highly connected with Ras-signalling and PI3K-Akt signalling pathways. Using The Cancer Genome Atlas (TCGA) data, we inferred the homeobox-regulated genes' expression is higher in 548 GBM cases than in 27 lower grade glioma cases giving that OLIG2 expression can be a reference. The results suggested that the H3K4me3 alteration is related to the formation and migration of GBM cells. We also found an extremely high reads count at epidermal growth factor receptor (EGFR) promoter, probably due to an amplification of copy number. Our analysis provides a case study about the change of H3K4me3 during shift to GBM.

A Xenopus gene whose expression can be activated by the organizer-specific homeobox genes goosecoid and Xnot2 was isolated by differential screening. The chordin gene encodes a novel protein of 941 amino acids that has a signal sequence and four Cys-rich domains. The expression of chordin starts in Spemann's organizer subsequent to that of goosecoid, and its induction by activin requires de novo protein synthesis. Microinjection of chordin mRNA induces twinned axes and can completely rescue axial development in ventralized embryos. This molecule is a potent dorsalizing factor that is expressed at the right time and in the right place to regulate cell-cell interactions in the organizing centers of head, trunk, and tail development.