Combinations of closely situated cis-acting elements determine tissue-specific patterns and anterior extent of early Hoxc8 expression

Identification of a Chondrocyte-Specific Enhancer in the Hoxc8 Gene

Hox genes encode transcription factors whose roles in patterning animal body plans during embryonic development are well-documented. Multiple studies demonstrate that Hox genes continue to act in adult cells, in normal differentiation, in regenerative processes, and, with abnormal expression, in diverse types of cancers. However, surprisingly little is known about the regulatory mechanisms that govern Hox gene expression in specific cell types, as they differentiate during late embryonic development, and in the adult organism. The murine Hoxc8 gene determines the identity of multiple skeletal elements in the lower thoracic and lumbar region and continues to play a role in the proliferation and differentiation of cells in cartilage as the skeleton matures. This study was undertaken to identify regulatory elements in the Hoxc8 gene that control transcriptional activity, specifically in cartilage-producing chondrocytes. We report that an enhancer comprising two 416 and 224 bps long interacting DNA elements produces reporter gene activity when assayed on a heterologous transcriptional promoter in transgenic mice. This enhancer is distinct in spatial, temporal, and molecular regulation from previously identified regulatory sequences in the Hoxc8 gene that control its expression in early development. The identification of a tissue-specific Hox gene regulatory element now allows mechanistic investigations into Hox transcription factor expression and function in differentiating cell types and adult tissues and to specifically target these cells during repair processes and regeneration.

HOXC8 regulates self-renewal, differentiation and transformation of breast cancer stem cells

BACKGROUND

Homeobox genes are master regulators of cell fate during embryonic development and their expression is altered in cancer. By regulating the balance between cell proliferation and differentiation, they maintain homeostasis of normal tissues. Here, we screened the expression of homeobox genes in mammary stem cells to establish their role in stem cells transformation in breast cancer.

METHODS

Using a Homeobox Genes PCR array, we screened 83 homeobox genes in normal cancer breast stem/progenitor cells isolated by flow cytometry. The candidate gene HOXC8 epigenetic regulation was studied by DNA methylation and miRNA expression analyses. Self-renewal and differentiation of HOXC8-overexpressing or knockdown cells were assessed by flow cytometry and mammosphere, 3D matrigel and soft agar assays. Clinical relevance of in vitro findings were validated by bioinformatics analysis of patient datasets from TCGA and METABRIC studies.

RESULTS

In this study we demonstrate altered expression of homeobox genes in breast cancer stem/progenitor cells. HOXC8 was consistently downregulated in stem/progenitor cells of all breast molecular subtypes, thus representing an interesting tumour suppressor candidate. We show that downregulated expression of HOXC8 is associated with DNA methylation at the gene promoter and expression of miR196 family members. Functional studies demonstrated that HOXC8 gain of function induces a decrease in the CD44⁺/CD24^-/low cancer stem cell population and proportion of chemoresistant cells, with a concomitant increase in CD24⁺ differentiated cells. Increased HOXC8 levels also decrease the ability of cancer cells to form mammospheres and to grow in anchorage-independent conditions. Furthermore, loss of HOXC8 in non-tumorigenic mammary epithelial cells expands the cancer stem/progenitor cells pool, increases stem cell self-renewal, prevents differentiation induced by retinoic acid and induces a transformed phenotype.

CONCLUSIONS

Taken together, our study points to an important role of homeobox genes in breast cancer stem/progenitor cell function and establishes HOXC8 as a suppressor of stemness and transformation in the mammary gland lineage.

Cdx is crucial for the timing mechanism driving colinear Hox activation and defines a trunk segment in the Hox cluster topology

Cdx and Hox transcription factors are important regulators of axial patterning and are required for tissue generation along the vertebrate body axis. Cdx genes have been demonstrated to act upstream of Hox genes in midgestation embryos. Here, we investigate the role of Cdx transcription factors in the gradual colinear activation of the Hox clusters. We found that Hox temporally colinear expression is severely affected in epiblast stem cells derived from Cdx null embryos. We demonstrate that after initiation of 3' Hox gene transcription, Cdx activity is crucial for H3K27ac deposition and for accessibility of cis-regulatory elements around the central - or 'trunk' - Hox genes. We thereby identify a Cdx-responsive segment of HoxA, immediately 5' to the recently defined regulatory domain orchestrating initial transcription of the first Hox gene. We propose that this partition of HoxA into a Wnt-driven 3' part and the newly found Cdx-dependent middle segment of the cluster, forms a structural fundament of Hox colinearity of expression. Subsequently to initial Wnt-induced activation of 3' Hox genes, Cdx transcription factors would act as crucial effectors for activating central Hox genes, until the last gene of the cluster arrests the process.

Embigin, regulated by HOXC8, plays a suppressive role in breast tumorigenesis

The transmembrane glycoprotein embigin (EMB) belongs to the immunoglobulin superfamily (IgSF) and a number of IgSF members have been identified as biomarkers for cancer progression. In this study, we show that embigin is transcriptionally regulated by Homeobox C8 (HOXC8) in breast cancer cells and embigin expression suppresses breast tumorigenesis. With aid of Western blot, luciferase reporter gene assay and chromatin immunoprecipitation, we reveal that HOXC8 binds to the EMB promoter at the region of nucleotides -2303 to -2315 and acts as a transcription inhibitor to suppress embigin expression. Depletion of embigin leads to increase in proliferation, anchorage-independent growth and migration of breast cancer cells, and the inhibitory effects mediated by HOXC8 knockdown on breast tumorigenesis can be largely rescued by depletion of embigin expression in breast cancer cells, suggesting that HOXC8 regulates breast tumorigenesis, at least partly, through regulating embigin expression. Moreover, we show that loss of embigin promotes proliferation, anchorage-independent growth, and migration ability of normal mammary epithelial MCF10A cells. The analyses of publically available human breast tumor microarray gene expression database show that low embigin levels correlate with short survival of breast tumor patients, particularly with basal-like tumor patients, and embigin expression is low specifically in patients with basal-like, ER-/HER2- tumors. Taken together, our study demonstrates that low/loss of embigin plays an important role in the progression of breast tumors.

Patterns and evolution of ACGT repeat cis-element landscape across four plant genomes

BACKGROUND

Transcription factor binding is regulated by several interactions, primarily involving cis-element binding. These binding sites maintain specificity by means of their sequence, and other additional factors such as inter-motif distance and spacer specificity. The ACGT core sequence has been established as a functionally important cis-element which frequently regulates gene expression in synergy with other cis-elements. In this study, we used two monocotyledonous - Oryza sativa and Sorghum bicolor, and two dicotyledonous species - Arabidopsis thaliana and Glycine max to analyze the conservation of co-occurring ACGT core elements in plant promoters with respect to spacer distance between them. Using data generated from Arabidopsis thaliana and Oryza sativa, we also identified conserved regions across all spacers and possible conditions regulating gene promoters with multiple ACGT cis-elements.

RESULTS

Our data indicated specific predominant spacer lengths between co-occurring ACGT elements, but these lengths were not universally conserved across all species under analysis. However, the frequency distribution indicated local regions of high correlation among monocots and dicots. Sequence specificity data clearly revealed a preference for G at the first and C at the terminal position of a spacer sequence, suggesting that the G-box motif is the most prevalent for the ACGT class of promoters. Using gene expression databases, we also observed trends suggesting that co-occurring ACGT elements are responsible for gene regulation in response to exogenous stress. Conservation in patterns of ACGT (N) ACGT among orthologous genes also indicated the possibility that emergence of functional significance across species was a result of parallel evolution of these cis-elements.

CONCLUSIONS

Although the importance of ACGT elements has been acknowledged for several plant species, ours is the first study that attempts to compare their occurrence across four species and analyze conservation among them. The apparent preference for particular spacer distances suggest that these motifs might be implicated in important physiological functions which are yet to be identified. Variations in correlation patterns among monocots and dicots might arise out of differences in transcriptional regulation in the two classes. In accordance with literature, we established the involvement of co-occurring ACGT elements in stress responses and showed how this regulation differs with variation in the ACGT (N) ACGT motif. We believe that our study will be an essential resource in determining optimum spacer length and spacer sequence between ACGT elements for promoter design in future.

Homeobox C9 is not potentially related to congenital heart disease in Chinese patients

BACKGROUND

Congenital heart disease (CHD) is one of the most common human birth defects. The etiology and pathogenesis of CHD are complex and involve several genes as well as multiple changes in signaling pathways. The aim of this study was to identify potential pathological mutations in the Homeobox C9 (Hoxc9) gene in 350 Chinese children with CHD to further understand the etiology of CHD.

METHOD

Sequence analysis of the Hoxc9 gene in 350 nonsyndromic patients with CHD Result: We did not identify any nonsynonymous variants in the coding regions of Hoxc9 in the patients with CHD. We found one synonymous variant c.C564T (p. his188his) in one ventricular septal defect patient. We also identified four previously reported polymorphisms (rs56368105, rs12817092, rs34079606, and rs2241820) in CHD.

CONCLUSIONS

We did not find any diagnostic alterations in the coding regions of Hoxc9 among the patients with CHD. Nevertheless, to our knowledge, this is the first study of Hoxc9 in nonsyndromic CHD and has expanded our overall knowledge of the etiology of this disease.

Shadow enhancers: frequently asked questions about distributed cis-regulatory information and enhancer redundancy

This paper, in the form of a frequently asked questions page (FAQ), addresses outstanding questions about "shadow enhancers", quasi-redundant cis-regulatory elements, and their proposed roles in transcriptional control. Questions include: What exactly are shadow enhancers? How many genes have shadow/redundant/distributed enhancers? How redundant are these elements? What is the function of distributed enhancers? How modular are enhancers? Is it useful to study a single enhancer in isolation? In addition, a revised definition of "shadow enhancers" is proposed, and possible mechanisms of shadow enhancer function and evolution are discussed.

Hoxc8 downregulates Mgl1 tumor suppressor gene expression and reduces its concomitant function on cell adhesion

Hoxc8 is a homeobox gene family member, which is essential for growth and differentiation. Mgl1, a mouse homologue of the Drosophila tumor suppressor gene lgl, was previously identified as a possible target of Hoxc8. However, the biological effects and underlying molecular mechanism of Hoxc8 regulation on Mgl1 has not been fully established. The endogenous expression patterns of Hoxc8 were inversely correlated with those of Mgl1 in different types of cells and tissues. Here we showed that Hoxc8 overexpression downregulated the Mgl1 mRNA expression. Characterization of the ~2 kb Mgl1 promoter region revealed that the upstream sequence contains several putative Hox core binding sites and chromatin immunoprecipitation assay confirmed that Hoxc8 directly binds to the 5' upstream region of Mgl1. The promoter activity of this region was diminished by Hoxc8 expression but resumed by knockdown of Hoxc8 using siRNA against Hoxc8. Functional study of Mgl1 in C3H10T1/2 cells revealed a significant reduction in cell adhesion upon expression of Hoxc8. Taken together, our data suggest that Hoxc8 downregulates Mgl1 expression via direct binding to the promoter region, which in turn reduces cell adhesion and concomitant cell migration.

Comparison of diverged Hoxc8 early enhancer activities reveals modification of regulatory interactions at conserved cis-acting elements

The Hoxc8 early enhancer that controls the initiation and establishment of Hoxc8 expression in the developing mouse embryo is found in different vertebrate lineages including mammals, birds and fish. Mouse and Fugu Hoxc8 early enhancers (200 bp) have diverged in the composition of elements located towards the 3' region. However, they share cis-acting elements A-E located in the 5' region. Mutations at these elements in the context of the mouse Hoxc8 early enhancer affect reporter gene expression in the posterior neural tube, somites and lateral plate mesoderm of day 9.5 mouse embryos. Here, we demonstrate that mutations introduced at the same elements but in the context of the Fugu Hoxc8 early enhancer had different consequences on the reporter gene expression in transgenic mouse embryos. Furthermore, in contrast to the mouse enhancer the Fugu enhancer does not utilize elements D and E in achieving posterior neural tube and somite expression. These results suggest that the diverged sequences prevent regulatory interactions at conserved cis-acting elements. We propose that divergent sequences modify regulatory interactions at conserved elements by providing a "contextual change". Our finding that the enhancer elements do not act in a unitary fashion but function in the context of the surrounding sequence brings a new dimension to the study of cis-regulatory evolution.

Mouse naked cuticle 2 (mNkd2) as a direct transcriptional target of Hoxc8 in vivo

Mouse naked cuticle 2 (mNkd2), the mammalian homolog of the Drosophila segment polarity gene naked cuticle (nkd), encodes an EF hand protein that regulates early Wg activity by acting as an inducible antagonist. The transcription factor, Hoxc8, a member of the homeobox gene family, is vital for growth and differentiation. Chromatin immunoprecipitation (ChIP) assay, an electrophoretic mobility shift assay (EMSA), and a reporter assay demonstrated that endogenous Hoxc8 protein binds directly to the enhancer region of the mNkd2 gene, implying a Hoxc8-dependent transcriptional activity. Introduction of exogenous Hoxc8 into NIH3T3 cell lines lacking wild-type Hoxc8 dramatically reduced expression of mNkd2 mRNA. If, as the results suggest, mNkd2 is a direct target of Hoxc8, it represents a novel mechanism by which Hoxc8 might cross-talk with the Wnt signaling pathway by regulating mNkd2.

Identification of regulatory elements by gene family footprinting and in vivo analysis

Gene families of recently duplicated but subsequently diverged genes provide an unique opportunity for comparative analysis of regulatory elements. We have studied the human SPRR gene family of small proline rich proteins involved in barrier function of stratified squamous epithelia. These genes are all expressed in normal human keratinocytes, but respond differently to environmental insults. Comparisons of the functional promoter regions allows the rapid identification of both conserved and of novel regulatory elements that appeared after gene duplication. Competitive electrophoretic mobility shift assays can be used to confirm their presence. Here we show the power of gene family footprinting by the identification of two novel elements in the SPRR3 promoter, not present in SPRR1A and SPRR2A. One of these elements binds a protein similar to GAAP-1, a pro-apoptotic activator of IRF-1 and p53. In vivo analysis shows that this element functions as an inhibitor of SPRR3 transcription. The second novel element functions as an activator of promoter activity and is characterized by its A/T rich sequence. The latter interacting protein indeed binds through contacts in the minor groove, and strikingly, depends on the presence of calcium for DNA interaction.

Evidence for positive and negative regulation of the mouseCdx2 gene

The caudal family of transcription factors specifies posterior structures during mouse development. We describe the cis-regulatory regions that control mouse Cdx2 expression in the posterior neural tube and mesoderm. An 11.4 kb genomic fragment directs reporter gene expression in a pattern reflecting endogenous Cdx2 expression. A crucial enhancer is located in a 1 kb fragment upstream of the Cdx2 transcriptional start site. This enhancer by itself directs reporter gene expression to more anterior levels in the neural tube compared to the endogenous Cdx2 expression, suggesting the presence of negative regulatory elements outside the 1 kb fragment. A second enhancer, located in the first intron directs robust expression to the posterior two-thirds of the developing embryo in a pattern that is ectopic to Cdx2 expression. The intronic enhancer activity is silenced in the context of the larger 11.4 kb reporter construct. Intron 1 contains two independent enhancers that specifically direct expression to mesoderm (MSE) and neural tube (NSE). Phylogenetic comparison of vertebrate Cdx2 sequences indicates several conserved regions of sequences within the three-enhancer regions. A transcription factor database search suggests potential binding sites for factors involved in FGF and Wnt signaling pathways.

Comparative genomics in vertebrate evolution and development

The vast quantities of publicly available DNA sequencing data and genome resources are enabling biologists to investigate age-old problems in biology that were not addressable previously. In this review, we discuss how comparative genomics is practiced and how the data can be used to make biological inferences with respect to vertebrate evolution and development. Examples are taken from the well-known HOX clusters, which are always a high-priority target for genomic analyses due to their inferred role in the evolution of metazoans. In addition, we briefly discuss the application of genomic approaches to problems in comparative endocrinology.

Identification of a Hoxc8-regulated transcriptional network in mouse embryo fibroblast cells

The transcription factor, Hoxc8, is a member of the homeobox gene family that is vital for growth and differentiation. Previously, we identified 34 genes whose expression levels were changed at least 2-fold by forced expression of Hoxc8 in C57BL/6J mouse embryo fibroblast cells using a mouse 16,463-gene oligonucleotide microarray. In the present study, we used the combined power of microarray profiling, global Hoxc8 DNA-binding site analysis, and high-throughput chromatin immunoprecipitation assays to identify direct and biologically relevant targets of Hoxc8 in vivo. Here we show that 19 of the 34 responsive genes contain Hoxc8 consensus DNA-binding sequence(s) in their regulatory regions. Chromatin immunoprecipitation analysis indicated that Hoxc8-DNA interaction was detected in five of the 19 candidate genes. All of these five target genes have been implicated in oncogenesis, cell adhesion, proliferation, and apoptosis. Overall, the genes described here should aid in the understanding of global regulatory networks of Hox genes and to provide valuable insight into the molecular basis of Hoxc8 in development and carcinogenesis.

Hoxc8 early enhancer of the Indonesian coelacanth, Latimeria menadoensis

Hoxc8 early enhancer controls the initiation and establishment phase of Hoxc8 expression in the mouse. Comparative studies indicate the presence of Hoxc8 early enhancer sequences in different vertebrate clades including mammals, birds and fish. Previous studies have shown differences between teleost and mammalian Hoxc8 early enhancers with respect to sequence and organization of protein binding elements. This raises the question of when the Hoxc8 early enhancer arose and how it has become modified in different vertebrate lineages. Here, we describe Hoxc8 early enhancer from the Indonesian coelacanth, Latimeria menadoensis. Coelacanths are the only extant lobefinned fish whose genome is tractable to genome analysis. The Latimeria Hoxc8 early enhancer sequence more closely resembles that of the mouse than that of Fugu or zebrafish. When assayed for enhancer activity by reporter gene analysis in transgenic mouse embryos, Latimeria Hoxc8 early enhancer directs expression to the posterior neural tube and mesoderm similar to that of the mouse enhancer. These observations support a close relationship between coelacanths and tetrapods and place the origin of a common Hoxc8 early enhancer sequence within the sarcopterygian lineage. The divergence of teleost (actinopterygii) Hoxc8 early enhancer may reflect a case of relaxed selection or other forms of instability induced by genome duplication events.

The identification of Hoxc8 target genes

Hox genes encode transcription factors that control spatial patterning during embryogenesis. To date, downstream targets of Hox genes have proven difficult to identify. Here, we describe studies designed to identify target genes under the control of the murine transcription factor Hoxc8. We used a mouse 16,463 gene oligonucleotide microarray to identify mRNAs whose expression was altered by the overexpression of Hoxc8 in C57BL/6J mouse embryo fibroblasts (MEF) in cell culture (in vitro). We identified a total of 34 genes whose expression was changed by 2-fold or greater: 16 genes were up-regulated, and 18 genes were down-regulated. The majority of genes encoded proteins involved in critical biological processes, such as cell adhesion, migration, metabolism, apoptosis, and tumorigenesis. Two genes showed high levels of regulation: (i) secreted phosphoprotein 1 (Spp1), also known as osteopontin (OPN), was down-regulated 4.8-fold, and (ii) frizzled homolog 2 (Drosophila) (Fzd2) was up-regulated 4.4-fold. Chromatin immunoprecipitation (ChIP) analysis confirmed the direct interaction between the OPN promoter and Hoxc8 protein in vivo, supporting the view that OPN is a direct transcriptional target of Hoxc8.

Comparative cis-regulatory analyses identify new elements of the mouse Hoxc8 early enhancer

The Hoxc8 early enhancer is a 200 bp region that controls the early phase of Hoxc8 expression during mouse embryonic development. This enhancer defines the domain of Hoxc8 expression in the neural tube and mesoderm of the posterior regions of the developing embryo. Five distinct cis-acting elements, A-E, were previously shown to govern early phase Hoxc8 expression. Significant divergence between mammalian and fish Hoxc8 early enhancer sequences and activities suggested additional cis-acting elements. Here we describe four additional cis-acting elements (F-I) within the 200 bp Hoxc8 early enhancer region identified by comparative regulatory analysis and transgene-mutation studies. These elements affect posterior neural tube and mesoderm expression of the reporter gene, either singly or in combination. Surprisingly, these new elements are missing from the zebrafish and Fugu Hoxc8 early enhancer sequences. Considering that fish enhancers direct robust reporter expression in transgenic mouse embryos, it is tempting to postulate that fish and mammalian Hoxc8 early enhancers utilize different sets of elements to direct Hoxc8 early expression. These observations reveal a remarkable plasticity in the Hoxc8 early enhancer, suggesting different modes of initiation and establishment of Hoxc8 expression in different species. We postulate that extensive restructuring and remodeling of Hox cis-regulatory regions occurring in different taxa lead to relatively different Hox expression patterns, which in turn may act as a driving force in generating diverse axial morphologies.

The anthropoid postcranial axial skeleton: Comments on development, variation, and evolution

Within-species phenotypic variation is the raw material on which natural selection acts to shape evolutionary change, and understanding more about the developmental genetics of intraspecific as well as interspecific phenotypic variation is an important component of the Evo-Devo agenda. The axial skeleton is a useful system to analyze from such a perspective. Its development is increasingly well understood, and between-species differences in functionally important developmental parameters are well documented. I present data on intraspecific variation in the axial postcranial skeleton of some Primates, including hominoids (apes and humans). Hominoid species are particularly valuable, because counts of total numbers of vertebrae, and hence original somite numbers, are available for large samples. Evolutionary changes in the axial skeleton of various primate lineages, including bipedal humans, are reviewed, and hypotheses presented to explain the changes in terms of developmental genetics. Further relevant experiments on model organisms are suggested in order to explore more fully the differences in developmental processes between primate species, and hence to test these hypotheses.

Divergence of Hoxc8 early enhancer parallels diverged axial morphologies between mammals and fishes

There is considerable interest in understanding how cis-regulatory modifications drive morphological changes across species. Because developmental regulatory genes, including Hox genes, are remarkably conserved, their noncoding regulatory regions are likely sources for variations. Modifications of Hox cis-regulatory elements have potential to alter Hox gene expression and, hence, axial morphologies. In vertebrates, differences in the axial levels of Hox gene expression correlate with differences in the number and relative position of thoracic vertebrae. Variation in cis-regulatory elements of Hox genes can be identified by comparative sequence and reporter gene analyses in transgenic mouse embryos. Using these approaches, we show a remarkable divergence of the Hoxc8 early enhancers between mammals and fishes representing diverse axial morphologies. Extensive restructuring of the Hoxc8 early enhancer including nucleotide substitutions, inversion, and divergence result in distinct patterns of reporter gene expression along the embryonic axis. Our results provide an evolutionary perspective on how the enhancer elements are engineered and support the hypothesis that remodeling of Hox regulatory elements in different species has played a significant role in generating morphological diversity.

Enhancer timing of Hox gene expression: deletion of the endogenous Hoxc8 early enhancer

The proper expression of Hox genes is necessary for the accurate patterning of the body plan. The elucidation of the developmental genetic basis of transcriptional regulation of Hox genes by the study of their cis-regulatory elements provides crucial information regarding the establishment of axial specification. In this report, we investigate the role of the early enhancer (EE) of the murine Hoxc8 gene to better understand its role in pattern formation. Previous reports show that knockouts of the endogenous Hoxc8 coding region result in a combination of neural, behavioral and skeletal phenotypes. In this report, we limit ourselves to a consideration of the skeletal abnormalities. Early reports from our laboratory based on exogenous transgenic reporter constructs implicate a 200 bp non-coding element 3 kb upstream of the Hoxc8 promoter as a crucial enhancer that regulates the transcription of Hoxc8. In the present work, we have deleted this regulatory region from the endogenous genome using embryonic stem cell technology. Our results show that the deletion of the EE results in a significant delay in the temporal expression of Hoxc8. We also show that the deletion of the EE does not eliminate the expression of the Hoxc8 protein, but delays the attainment of control levels of expression and anterior and posterior boundaries of expression on the AP axis. The temporal delay in Hoxc8 expression is sufficient to produce phenocopies of many of the axial skeletal defects associated with the complete absence of Hoxc8 gene product as previously reported for the Hoxc8-null mutation. Our results are consistent with emerging evidence that the precise temporal expression of Hox genes is crucial for the establishment of regional identities. The fact that the EE deletion does not eliminate Hoxc8 expression indicates the existence of a Hoxc8 transcriptional regulatory apparatus independent to some degree of the Hoxc8 EE. In a comparison of our results with those reported previously by others investigating temporal control of Hox gene expression, we have discovered a structural similarity between the Hoxc8 EE reported here and a transcriptional control element located in the Hoxd11 region. We speculate that a distributed system of expression timing control may exist that is similar the one we propose for Hoxc8. Last, our data is consistent with the position that disparate regulatory pathways are responsible for the expression of Hoxc8 in the organogenesis of somites, neural tube and limb bud.

Developmental Biology and Human Evolution

▪ Abstract  Our understanding of developmental biology burgeoned during the last decade. This review summarizes recent advances, provides definitions and explanations of some basic principles, and does so in a way that will aid anthropologists in understanding their profound implications. Crucial concepts, such as developmental fields, selector and realizator genes, cell signaling mechanisms, and gene regulatory elements are briefly described and then integrated with the emergence of skeletal morphology. For the postcranium, a summary of events from limb bud formation, the appearance of anlagen, the expression of Hox genes, and the fundamentals of growth plate dynamics are briefly summarized. Of particular importance are revelations that bony morphology is largely determined by pattern formation, that growth foci such as physes and synovial joints appear to be regulated principally by positional information, and that variation in these fields is most likely determined by cis-regulatory elements acting on restricted numbers of anabolic genes downstream of selectors (such as Hox). The implications of these discoveries for the interpretation of both contemporary and ancient human skeletal morphology are profound. One of the most salient is that strain transduction now appears to play a much reduced role in shaping the human skeleton. Indeed, the entirety of “Wolff's Law” must now be reassessed in light of new knowledge about pattern formation. The review concludes with a brief discussion of some implications of these findings, including their impact on cladistics and homology, as well as on biomechanical and morphometric analyses of both ancient and modern human skeletal material.

Phylogenetic analysis of the mammalian Hoxc8 non-coding region

The non-coding intergenic regions of Hox genes are remarkably conserved among mammals. To determine the usefulness of this sequence for phylogenetic comparisons, we sequenced an 800-bp fragment of the Hoxc9-Hoxc8 intergenic region from several species belonging to different mammalian clades. Results obtained from the phylogenetic analysis are congruent with currently accepted mammalian phylogeny. Additionally, we found a TC mini satellite repeat polymorphism unique to felines. This polymorphism may serve as a useful marker to differentiate between mammalian species or as a genetic marker in feline matings. This study demonstrates usefulness of a comparative approach employing non-coding regions of Hox gene complexes.

Methylation of HoxA5 and HoxB5 and its relevance to expression during mouse development

Expression and function of homeobox genes (Hox genes) in development have been subject to extensive study in a variety of organisms including mammals, however practically nothing is known regarding the methylation patterns of these genes. Here we describe the methylation patterns of HoxA5 and HoxB5 in various tissues of fetal and adult mice and their relevance to expression. Both genes exhibit tissue specific methylation patterns that are established postnatally. This methylation appears to play a role in stabilizing the newly acquired silent state of the genes. In contrast to the postimplantation wave of de novo methylation that takes place across the mammalian genome, the methylation of the Hox genes represents a different time window for de novo methylation which might be characteristic of developmental genes. In the case of HoxA5 this postnatal de novo methylation can cover a domain of at least 25 kb that includes several genes of the HoxA cluster and the CpG islands within. Our observations suggest that the establishment of tissue specific methylation patterns of HoxA5 and HoxB5 and the relationship between these methylation patterns and activity are different from what had been known for non-developmental genes. This may reflect the specialized functions played by Hox genes in development.

The caudal homeodomain protein activates Drosophila E2F gene expression

The Drosophila caudal homeobox gene is required for definition of the anteroposterior axis and for gut development, and CDX1 and CDX2, human homologs, play crucial roles in the regulation of cell proliferation and differentiation in the intestine. Most studies have indicated tumor suppressor functions of Cdx2, with inhibition of proliferation, while the effects of Cdx1 are more controversial. The influence of Drosophila Caudal on cell proliferation is unknown. In this study, we found three potential Caudal binding sequences in the 5'-flanking region of the Drosophila E2F (DE2F) gene and showed by transient transfection assays that they are involved in Caudal transactivation of the dE2F gene promoter. Analyses with transgenic flies carrying an E2F-lacZ fusion gene, with and without mutation in the Caudal binding site, indicated that the Caudal binding sites are required for expression of dE2F in living flies. Caudal-induced E2F expression was also confirmed with a GAL4-UAS system in living flies. In addition, ectopic expression of Caudal with heat-shock promotion induced melanotic tumors in larvae. These results suggest that Caudal is involved in regulation of proliferation through transactivation of the E2F gene in Drosophila.

Hoxc-8 expression shows left-right asymmetry in the posterior lateral plate mesoderm

XHoxc-8 is the Xenopus homologue of the mouse Hoxc-8 gene, a homeodomain-containing transcription factor that is expressed in the posterior neural tube and adjacent tissues. Although XHoxc-8 has a very similar expression pattern to the Hoxc-8 gene in other species, it also displays distinct left/right asymmetry at later stages of development, being expressed in the posterior lateral mesoderm only on the left-hand side of the embryo.

Homeobox genes in gut development

Classical descriptions of gut development specify subdivision into foregut, midgut, and hindgut together with their derivatives. This is based on the anatomical localisation of the anterior and posterior intestinal portals separating the roof of the yolk sac from the foregut and hindgut diverticulae. When considering the molecular basis of intestinal differentiation, it is necessary to think in terms of the genes involved, and in this respect those containing the homeobox motif are important players in specifying the fate of both the endodermal and mesodermal components of the gut. In this review, evidence is considered for their role, with particular regard to the acquisition of positional information.

Spatially specific expression of Hoxb4 is dependent on the ubiquitous transcription factor NFY

Understanding how boundaries and domains of Hox gene expression are determined is critical to elucidating the means by which the embryo is patterned along the anteroposterior axis. We have performed a detailed analysis of the mouse Hoxb4 intron enhancer to identify upstream transcriptional regulators. In the context of an heterologous promoter, this enhancer can establish the appropriate anterior boundary of mesodermal expression but is unable to maintain it, showing that a specific interaction with its own promoter is important for maintenance. Enhancer function depends on a motif that contains overlapping binding sites for the transcription factors NFY and YY1. Specific mutations that either abolish or reduce NFY binding show that it is crucial for enhancer activity. The NFY/YY1 motif is reiterated in the Hoxb4 promoter and is known to be required for its activity. As these two factors are able to mediate opposing transcriptional effects by reorganizing the local chromatin environment, the relative levels of NFY and YY1 binding could represent a mechanism for balancing activation and repression of Hoxb4 through the same site.

Clock regulatory elements control cyclic expression of Lunatic fringe during somitogenesis

Somitogenesis requires a segmentation clock and Notch signaling. Lunatic fringe (Lfng) expression in the presomitic mesoderm (PSM) cycles in the posterior PSM, is refined in the segmenting somite to the rostral compartment, and is required for segmentation. We identify distinct cis-acting regulatory elements for each aspect of Lfng expression. Fringe clock element 1 (FCE1) represents a conserved 110 bp region that is necessary to direct cyclic Lfng RNA expression in the posterior PSM. Mutational analysis of E boxes within FCE1 indicates a potential interplay of positive and negative transcriptional regulation by cyclically expressed bHLH proteins. A separable Lfng regulatory region directs expression to the prospective rostral aspect of the condensing somite. These independent Lfng regulatory cassettes advance a molecular framework for deciphering somite segmentation.

Evolution of functionally conserved enhancers can be accelerated in large populations: a population-genetic model

The evolution of cis-regulatory elements (or enhancers) appears to proceed at dramatically different rates in different taxa. Vertebrate enhancers are often very highly conserved in their sequences, and relative positions, across distantly related taxa. In contrast, functionally equivalent enhancers in closely related Drosophila species can differ greatly in their sequences and spatial organization. We present a population-genetic model to explain this difference. The model examines the dynamics of fixation of pairs of individually deleterious, but compensating, mutations. As expected, small populations are predicted to have a high rate of evolution, and the rate decreases with increasing population size. In contrast to previous models, however, this model predicts that the rate of evolution by pairs of compensatory mutations increases dramatically for population sizes above several thousand individuals, to the point of greatly exceeding the neutral rate. Application of this model predicts that species with moderate population sizes will have relatively conserved enhancers, whereas species with larger populations will be expected to evolve their enhancers at much higher rates. We propose that the different degree of conservation seen in vertebrate and Drosophila enhancers may be explained solely by differences in their population sizes and generation times.

An efficient cis-element discovery method using multiple sequence comparisons based on evolutionary relationships

The discovery of cis-element control motifs in noncoding DNA poses a difficult problem in genome analysis. Functional analysis by means of reporter constructs expressed in transgenic organisms is the most reliable method, but is by itself time-consuming and expensive. Searching noncoding DNA for known control motifs by sequence analysis is problematic, since protein binding motifs are short, in the range of 8-10 bp, and occur frequently by chance. Heretofore, the most reliable sequence analysis method has been the comparison of homologous sequence domains in related but moderately evolutionarily divergent species such as, for example, mouse and human. In such pairwise combinations, control regions are conserved because they serve a vital function and can be identified by their similar sequences. Single pairwise comparisons, however, allow the discovery of conserved sequence strings only at low resolution and without specific identity. We have investigated the possibility of using multiple sequence comparisons to correct these shortcomings. We applied this method to the Hoxc8 early enhancer region that has been previously analyzed in depth by functional methods and through its application successfully identified known protein binding cis-element motifs. Candidate protein binding sites could also be identified. This method, based on evolutionarily related sequence comparisons, should be quite useful as a prescreening step prior to functional analysis with corresponding savings in time and resources.

Hox genes and morphological identity: axial versus lateral patterning in the vertebrate mesoderm

The successful organization of the vertebrate body requires that local information in the embryo be translated into a functional, global pattern. Somite cells form the bulk of the musculoskeletal system. Heterotopic transplants of segmental plate along the axis from quail to chick were performed to test the correlation between autonomous morphological patterning and Hox gene expression in somite subpopulations. The data presented strengthen the correlation of Hox gene expression with axial specification and focus on the significance of Hox genes in specific derivatives of the somites. We have defined two anatomical compartments of the body based on the embryonic origin of the cells making up contributing structures: the dorsal compartment, formed from purely somitic cell populations; and the ventral compartment comprising cells from somites and lateral plate. The boundary between these anatomical compartments is termed the somitic frontier. Somitic tissue transplanted between axial levels retains both original Hox expression and morphological identity in the dorsal compartment. In contrast, migrating lateral somitic cells crossing the somitic frontier do not maintain donor Hox expression but apparently adopt the Hox expression of the lateral plate and participate in the morphology appropriate to the host level. Dorsal and ventral compartments, as defined here, have relevance for experimental manipulations that influence somite cell behavior. The correlation of Hox expression profiles and patterning behavior of cells in these two compartments supports the hypothesis of independent Hox codes in paraxial and lateral plate mesoderm.

Isolation and characterization of BEN, a member of the TFII-I family of DNA-binding proteins containing distinct helix-loop-helix domains

The transcriptional regulation of the Hoxc8 gene is controlled during early mouse embryogenesis by an enhanceosome-like control region, termed the early enhancer (EE), located 3 kb upstream from the Hoxc8 translation start site. The EE is involved in establishing the posterior expression pattern of Hoxc8 at embryonic day (E) 8.5-9. 0. Genetic and biochemical data have shown that nuclear factors interact with this region in a sequence-specific manner. We have used a yeast one-hybrid screen in a search for transcription factors that bind to EE motifs and have isolated a novel murine DNA-binding protein, termed BEN (binding factor for early enhancer). The ORF of BEN encodes a protein of 1072 amino acids and contains six helix-loop-helix domains, a hydrophobic leucine zipper-like motif, and a serine-rich repeat. The murine BEN gene is structurally similar to the human gene TFII-I in that both genes encode unique 95-amino acid long helix-loop/span-helix domains. The BEN gene produces several major transcripts (3.6, 4.4, and 5.9 kb) present in most adult tissues and shows discrete spatial and temporal domains of expression in areas of epithelial-mesenchymal interaction during mouse embryogenesis from E9.5 to E12.5. Several BEN-encoded polypeptides of different sizes ranging from 165 to 40 kDa were identified by Western blot analysis using BEN-specific polyclonal Abs. We propose, on the bases of sequence homology, that BEN is the mouse ortholog of the recently described human gene, WBSCR11, known also as GTF2IRD1, GTF3, Cream1, and MusTRD1. This gene is deleted hemizygously in individuals with Williams Syndrome, an autosomal dominant genetic condition characterized by complex physical, cognitive, and behavioral traits resulting from a perturbed developmental process.

Homeobox genes and gut development

The gut of vertebrates exhibits a common anteroposterior regional differentiation. The role of homeobox genes in establishing this pattern is inferred by their sites of expression. It is suggested that the primary source of positional information is in the endoderm, which subsequently establishes a 'dialogue' with the surrounding visceral layer of the lateral plate mesoderm. This results in the anatomical and physiological specialization of the adult gut.

Evolution of chordate hox gene clusters

In this article, we consider the role of the Hox genes in chordate and vertebrate evolution from the viewpoints of molecular and developmental evolution. Models of Hox cluster duplication are considered with emphasis on a threefold duplication model. We also show that cluster duplication is consistent with a semiconservative model of duplication, where following duplication, one daughter cluster remains unmodified, while the other diverges and assumes a new architecture and presumably new functions. Evidence is reviewed, suggesting that Hox gene enhancers have played an important role in body plan evolution. Finally, we contrast the invertebrates and vertebrates in terms of genome and Hox cluster duplication which are present in the latter, but not the former. We question whether gene duplication has been important in vertebrates for the introduction of novel features such as limbs, a urogenital system, and specialized neuromuscular interactions.

Comparative studies on mammalian Hoxc8 early enhancer sequence reveal a baleen whale-specific deletion of a cis-acting element

Variations in regulatory regions of developmental control genes have been implicated in the divergence of axial morphologies. To find potentially significant changes in cis-regulatory regions, we compared nucleotide sequences and activities of mammalian Hoxc8 early enhancers. The nucleotide sequence of the early enhancer region is extremely conserved among mammalian clades, with five previously described cis-acting elements, A-E, being invariant. However, a 4-bp deletion within element C of the Hoxc8 early enhancer sequence is observed in baleen whales. When assayed in transgenic mouse embryos, a baleen whale enhancer (unlike other mammalian enhancers) directs expression of the reporter gene to more posterior regions of the neural tube but fails to direct expression to posterior mesoderm. We suggest that regulation of Hoxc8 in baleen whales differs from other mammalian species and may be associated with variation in axial morphology.

Murine Hoxc-9 gene contains a structurally and functionally conserved enhancer

Reporter gene analysis of the Hoxc-9 genomic region in transgenic mice allowed us to identify a positional enhancer in the Hoxc-9 intron that drives expression in the posterior neural tube of midgestation mouse embryos in a Hoxc-9-related manner. Sequence comparison to the chicken Choxc-9 intron revealed the existence of two highly conserved sequence elements (CSEs) in a similar spatial arrangement. These structural similarities in the mammalian and avian lineage are mirrored by conserved function of the chicken Choxc-9 intron in transgenic mice. Deletion analysis of the two introns suggests that full activity of both enhancers depends on cooperation between the two CSEs located close to the respective 5' and 3' splice sites. Following the paradigm of phylogenetically conserved developmental control mechanisms, the Hoxc-9 intragenic enhancer was tested in Drosophila. Our data show that the mouse Hoxc-9 enhancer acts in a conserved fashion in transgenic flies, conferring posteriorly restricted reporter gene expression to the developing central nervous system in third instar larvae. This finding indicates that the Hoxc-9 intragenic enhancer is involved in transcriptional regulatory circuits conserved between vertebrates and arthropods.

Expression of the murine Hoxa4 gene requires both autoregulation and a conserved retinoic acid response element

Analysis of the regulatory regions of the Hox genes has revealed a complex array of positive and negative cis-acting elements that control the spatial and temporal pattern of expression of these genes during embryogenesis. In this study we show that normal expression of the murine Hoxa4 gene during development requires both autoregulatory and retinoic acid-dependent modes of regulation. When introduced into a Hoxa4 null background, expression of a lacZ reporter gene driven by the Hoxa4 regulatory region (Hoxa4/lacZ) is either abolished or significantly reduced in all tissues at E10. 5-E12.5. Thus, the observed autoregulation of the Drosophila Deformed gene is conserved in a mouse homolog in vivo, and is reflected in a widespread requirement for positive feedback to maintain Hoxa4 expression. We also identify three potential retinoic acid response elements in the Hoxa4 5′ flanking region, one of which is identical to a well-characterized element flanking the Hoxd4 gene. Administration of retinoic acid to Hoxa4/lacZ transgenic embryos resulted in stage-dependent ectopic expression of the reporter gene in the neural tube and hindbrain. When administered to Hoxa4 null embryos, however, persistent ectopic expression was not observed, suggesting that autoregulation is required for maintenance of the retinoic acid-induced expression. Finally, mutation of the consensus retinoic acid response element eliminated the response of the reporter gene to exogenous retinoic acid, and abolished all embryonic expression in untreated embryos, with the exception of the neural tube and prevertebrae. These data add to the evidence that Hox gene expression is regulated, in part, by endogenous retinoids and autoregulatory loops.

Modification of expression and cis-regulation of Hoxc8 in the evolution of diverged axial morphology

Differential Hox gene expression between vertebrate species has been implicated in the divergence of axial morphology. To examine this relationship, we have compared expression and transcriptional regulation of Hoxc8 in chicken and mouse. In both species, expression of Hoxc8 in the paraxial mesoderm and neural tube is associated with midthoracic and brachial identities, respectively. During embryogenesis, there is a temporal delay in the activation of Hoxc8 in chicken compared with mouse. As a result, chicken Hoxc8 expression in the paraxial mesoderm is at a posterior axial level, extending over a smaller domain compared with mouse Hoxc8 expression. This finding is consistent with a shorter thoracic region in chicken compared with mouse. In addition, the chicken Hoxc8 early enhancer, differing from its mouse counterpart in only a few specific nucleotides, directs a reporter gene expression to a more posterior domain in transgenic mouse embryos. These findings are consistent with the concept that the diversification of axial morphology has been achieved through changes in cis-regulation of developmental control genes.

A regulatory element of the empty spiracles homeobox gene is composed of three distinct conserved regions that bind regulatory proteins

Homeobox genes encode a class of highly evolutionarily conserved transcription factors that control embryonic development. The Drosophila melanogaster empty spiracles gene is the homolog of the two human homeobox genes EMX1 and EMX2. These genes are necessary for central nervous system development. We used a regulatory element of the empty spiracles gene to study the control of homeobox gene expression. The 1.2-kilobase (kb) cis-regulatory element located 3 kb 5' of the transcription start site of the empty spiracles gene was analyzed by evolutionary sequence comparisons, gel mobility shift assays, DNase footprinting, and the generation of transgenic flies. The corresponding element from a related species, Drosophila hydei, was cloned. Three discrete, approximately 100 base pair (bp) regions of sequence homology were identified. Each had two blocks of 10 to 40 bp of near perfect sequence identity. Fusion proteins were produced containing the Abdominal-B homeodomain or the empty spiracles homeodomain, known regulators of empty spiracles gene expression. Gel mobility shift assays showed that each of the three regions is bound by both proteins. DNase footprinting revealed closely linked empty spiracles and Abdominal-B binding sites. We then generated transgenic flies containing a reporter linked to individual conserved regions of the enhancer. Reporter expression was evident only outside of the usual empty spiracles expression domain. These elements are not sufficient alone; a combinatorial model is proposed. Conserved discrete areas within a homeobox gene regulatory element, which function as homeodomain protein transcription factor binding sites, are used in a combinatorial fashion to regulate these developmentally important genes.

Molecular evolution of Hox gene regulation: cloning and transgenic analysis of the lamprey HoxQ8 gene

The mammalian Hox clusters arose by duplication of a primordial cluster. The duplication of Hox clusters created redundancy within cognate groups, allowing for change in function over time. The lamprey, Petromyzon marinus, occupies an intermediate position within the chordates, both in terms of morphologic complexity and possibly cluster number. To determine the extent of divergence among Hox genes after duplication events within vertebrates, we analyzed Hox genes belonging to cognate group 8. Here we report characterization of the HoxQ8 gene, which shows conservation with mammalian genes in its amino-terminal, homeobox and hexapeptide sequences, and in the position of its splice sites. A beta-galactosidase reporter gene was introduced in the HoxQ8 genomic region by targeted recombinational cloning using a yeast-bacteria shuttle vector, pClasper. These reporter gene constructs were tested for their ability to direct region-specific expression patterns in transgenic mouse embryos. Lamprey enhancers direct expression to posterior neural tube but not to mesoderm, suggesting conservation of neuronal enhancers. In the presence of the mouse heat shock promoter, lamprey enhancers could also direct expression to the posterior mesoderm suggesting that there has been some divergence in promoter function. Our results suggest that comparative studies on Hox gene structure and analysis of regulatory elements may provide insights into changes concomitant with Hox cluster duplications in the chordates.

Interspecies exchange of a Hoxd enhancer in vivo induces premature transcription and anterior shift of the sacrum

The precise activation, in space and time, of vertebrate Hox genes is an essential requirement for normal morphogenesis. In order to assess for the functional potential of evolutionary conserved Hox regulatory sequences, a phylogenetically conserved bipartite regulatory element necessary for proper spatial and temporal activation of the Hoxd-11 gene was replaced by its fish counterpart in the HoxD complex of mice, using an ES cell-based targeted exchange. Fetuses carrying this replacement activated Hoxd-11 transcription prematurely, which led to a rostral shift of its expression boundary and a consequent anterior transposition of the sacrum. These results demonstrate the high phylogenetic conservation of regulatory mechanisms acting over vertebrate Hox complexes and suggest that minor time difference (heterochronies) in Hox gene activation may have contributed to important morphological variations in the course of evolution.

Deletion of a HoxD enhancer induces transcriptional heterochrony leading to transposition of the sacrum

A phylogenetically conserved transcriptional enhancer necessary for the activation of Hoxd-11 was deleted from the HoxD complex of mice by targeted mutagenesis. While genetic and expression analyses demonstrated the role of this regulatory element in the activation of Hoxd-11 during early somitogenesis, the function of this gene in developing limbs and the urogenital system was not affected, suggesting that Hox transcriptional controls are different in different axial structures. In the trunk of mutant embryos, transcriptional activation of Hoxd-11 and Hoxd-10 was severely delayed, but subsequently resumed with appropriate spatial distributions. The resulting caudal transposition of the sacrum indicates that proper vertebral specification requires a precise temporal control of Hox gene expression, in addition to spatial regulation. A slight time delay in expression (transcriptional heterochrony) cannot be compensated for at a later developmental stage, eventually leading to morphological alterations.

Activation of enhancer elements by the homeobox gene Cdx2 is cell line specific

Cdx2 is a caudal-related homeodomain transcription factor that is expressed in complex patterns during mouse development and at high levels in the intestinal epithelium of adult mice. Cdx2 activates transcription of intestinal gene promoters containing specific binding sites. Moreover, Cdx2 has been shown to induce intestinal differentiation in cell lines. In this study, we show that Cdx2 is able to bind to two well defined enhancer elements in the HoxC8 gene. We then demonstrate that Cdx2 is able to activate transcription of heterologous promoters when its DNA binding element is placed in an enhancer context. Furthermore, the ability to activate enhancer elements is cell-line dependent. When the Cdx2 activation domain was linked to the Gal4 DNA binding domain, the chimeric protein was able to activate Gal4 enhancer constructs in an intestinal cell line, but was unable to activate transcription in NIH3T3 cells. These data suggest that there are cell-specific factors that allow the Cdx2 activation domain to function in the activation of enhancer elements. We hypothesize that either a co-activator protein or differential phosphorylation of the activation domain may be the mechanism for intestinal cell line-specific function of Cdx2 and possibly in other tissues in early development.

The hardwiring of development: organization and function of genomic regulatory systems

The gene regulatory apparatus that directs development is encoded in the DNA, in the form of organized arrays of transcription factor target sites. Genes are regulated by interactions with multiple transcription factors and the target sites for the transcription factors required for the control of each gene constitute its cis-regulatory system. These systems are remarkably complex. Their hardwired internal organization enables them to behave as genomic information processing systems. Developmental gene regulatory networks consist of the cis-regulatory systems of all the relevant genes and the regulatory linkages amongst them. Though there is yet little explicit information, some general properties of genomic regulatory networks have become apparent. The key to understanding how genomic regulatory networks are organized, and how they work, lies in experimental analysis of cis-regulatory systems at all levels of the regulatory network.