Hox Gene Collinearity Research Articles

Hox Gene Collinearity with Pulling Physical Forces Creates a Hox Gene Clustering in Embryos of Vertebrates and Invertebrates: Complete or Split Clusters

Hox gene clusters are crucial in embryogenesis. It was observed that some Hox genes are located in order along the telomeric to centromeric direction of the DNA sequence: Hox1, Hox2, Hox3…. These genes are expressed in the same order in the ontogenetic units of the Drosophila embryo along the anterior–posterior axis. The two entities (genome and embryo) differ significantly in linear size and in-between distance. This strange phenomenon was named spatial collinearity (SP). Later, it was observed that, particularly in the vertebrates, a temporal collinearity (TC) coexists: first Hox1 is expressed, later Hox2, and later on Hox3…. According to a biophysical model (BM), pulling forces act at the anterior end of the cluster while a cluster fastening applies at the posterior end. Hox clusters are irreversibly elongated along the force direction. During evolution, the elongated Hox clusters are broken at variable lengths, thus split clusters may be created. An empirical rule was formulated, distinguishing development due to a complete Hox cluster from development due to split Hox clusters. BM can explain this empirical rule. In a spontaneous mutation, where the cluster fastening is dismantled, a weak pulling force automatically shifts the cluster inside the Hox activation domain. This cluster translocation can probably explain the absence of temporal collinearity in Drosophila.

Physical Laws Shape Up HOX Gene Collinearity.

Hox gene collinearity (HGC) is a multi-scalar property of many animal phyla particularly important in embryogenesis. It relates entities and events occurring in Hox clusters inside the chromosome DNA and in embryonic tissues. These two entities differ in linear size by more than four orders of magnitude. HGC is observed as spatial collinearity (SC), where the Hox genes are located in the order (Hox1, Hox2, Hox3 …) along the 3′ to 5′ direction of DNA in the genome and a corresponding sequence of ontogenetic units (E1, E2, E3, …) located along the Anterior—Posterior axis of the embryo. Expression of Hox1 occurs in E1, Hox2 in E2, Hox3 in E3, etc. Besides SC, a temporal collinearity (TC) has been also observed in many vertebrates. According to TC, first Hox1 is expressed in E1; later, Hox2 is expressed in E2, followed by Hox3 in E3, etc. Lately, doubt has been raised about whether TC really exists. A biophysical model (BM) was formulated and tested during the last 20 years. According to BM, physical forces are created which pull the Hox genes one after the other, driving them to a transcription factory domain where they are transcribed. The existing experimental data support this BM description. Symmetry is a physical–mathematical property of matter that was explored in depth by Noether who formulated a ground-breaking theory (NT) that applies to all sizes of matter. NT may be applied to biology in order to explain the origin of HGC in animals developing not only along the A/P axis, but also to animals with circular symmetry.

Hox Gene Collinearity May Be Related to Noether Theory on Symmetry and Its Linked Conserved Quantity

Hox Gene Collinearity (HGC) is a fundamental property that controls the development of many animal species, including vertebrates. In the Hox gene clusters, the genes are located in a sequential order Hox1, Hox2, Hox3, etc., along the 3’ to 5’ direction of the cluster in the chromosome. During Hox cluster activation, the Hox genes are expressed sequentially in the ontogenetic units D1, D2, D3, etc., along the anterior–posterior axis (A-P) of the early embryo. This collinearity, first observed by E.B. Lewis, is surprising because the spatial collinearity of these structures (Hox clusters and embryos) correlates entities that differ by about four orders of magnitude. Biomolecular mechanisms alone cannot explain such correlations. Long-range physical interactions, such as diffusion or electric attractions, should be involved. A biophysical model (BM) was formulated, which, in alignment with the biomolecular processes, successfully describes the existing vertebrate genetic engineering data. One hundred years ago, Emmy Noether made a fundamental discovery in mathematics and physics. She proved, rigorously, that a physical system obeying a symmetry law (e.g., rotations or self-similarity) is followed by a conserved physical quantity. It is argued here that HGC obeys a ‘primitive’ self-similarity symmetry. In this case, the associated primitive conserved quantity is the irreversibly increasing ‘ratchet’-like Hoxgene ordering where some genes may be missing. The genes of a vertebrate Hox clusterare located along a finite straight line. The same order follows the ontogenetic unitsof the vertebrate embryo. Therefore, HGC is a manifestation of a primitive Noether Theory (NT). NT may be applied to other than the vertebrate case, for instance, to animals with a circular topological symmetry. For example, the observed abnormal Hox gene ordering of the echinoderm Hox clusters may be reproduced by a double-strand break of the circular Hox gene ordering and its subsequent incorporation in the flanking chromosome.

Why a Constant Number of Vertebrae? Digital Control of Segmental Identity during Vertebrate Development: The Somite Cycle Controls a Digital, Chromatin-Based Counter That Defines Segmental Identity and Body Plans in Vertebrate Animals.

It is not understood how the numbers and identities of vertebrae are controlled during mammalian development. The remarkable robustness and conservation of segmental numbers may suggest the digital nature of the underlying process. The study proposes a mechanism that allows cells to obtain and store the segmental information in digital form, and to produce a pattern of chromatin accessibility that in turn regulates Hox gene expression specific to the metameric segment. The model requires that a regulatory element be present such that the number of occurrences of the motif between two consecutive Hox genes equals the number of segments under the control of the anterior gene. This is true for the recently discovered hydroxyl radical cleavage 3bp-periodic (HRC3) motif, associated with histone modifications and developmental genes. The finding not only allows the correct prediction of the numbers of segments using only sequence information, but also resolves the 40-year-old enigma of the function of temporal and spatial collinearity of Hox genes. The logic of the mechanism is illustrated in the attached animated video. How different aspects of the proposed mechanism can be tested experimentally is also discussed.

De novo transcription of multiple Hox cluster genes takes place simultaneously in early Xenopus tropicalis embryos.

ABSTRACThox genes are found as clusters in the genome in most bilaterians. The order of genes in the cluster is supposed to be correlated with the site of expression along the anterior-posterior body axis and the timing of expression during development, and these correlations are called spatial and temporal collinearity, respectively. Here we studied the expression dynamics of all hox genes of the diploid species Xenopus tropicalis in four Hox clusters (A–D) by analyzing high-temporal-resolution RNA-seq databases and the results showed that temporal collinearity is not supported, which is consistent with our previous data from allotetraploid Xenopus laevis. Because the temporal collinearity hypothesis implicitly assumes the collinear order of gene activation, not mRNA accumulation, we determined for the first time the timing of when new transcripts of hox genes are produced, by detecting pre-spliced RNA in whole embryos with reverse transcription and quantitative PCR (RT-qPCR) for all hoxa genes as well as several selected hoxb, hoxc and hoxd genes. Our analyses showed that, coinciding with the RNA-seq results, hoxa genes started to be transcribed in a non-sequential order, and found that multiple genes start expression almost simultaneously or more posterior genes could be expressed earlier than anterior ones. This tendency was also found in hoxb and hoxc genes. These results suggest that temporal collinearity of hox genes is not held during early development of Xenopus.

Retinoic acid and CTCF play key roles in inducing the collinear expression of the Hoxa cluster.

During the development of an embryo, the initiation of the collinear expression of Hox genes is essential for the proper formation of the anteroposterior body axis. Retinoic acid (RA), a natural derivative of vitamin A, plays a role in vertebrate development by regulating Hox gene expression. CCCTC-binding factor (CTCF), an insulator protein that controls gene transcription, also regulates the expression of Hox genes by binding to the CTCF-binding sites (CBSs). It has been reported that upon RA signaling, retinoic acid response elements (RAREs) located in the Hox clusters become occupied. Interestingly, RAREs exist in close proximity with CBSs, and therefore when RA is bound, CTCF cannot bind. Without CTCF and its insulator activities, the repressive domain in the chromatin becomes open for gene transcription. Here, we examine the relationship between RA and CTCF during the RA-induced expression of the Hoxa cluster genes, using F9 murine embryonic teratocarcinoma cells as a model system. We treated F9 cells with RA for different time, confirmed the collinear expression of Hoxa genes, and validated CTCF-binding in F9 cells as well as in CTCF-overexpressing F9 cells, in the presence of RA. The present study suggests that RA and CTCF pose antagonistic effects on each other during vertebrate development to attain Hox gene collinearity.

Physical Forces May Cause the HoxD Gene Cluster Elongation.

Hox gene collinearity was discovered be Edward B. Lewis in 1978. It consists of the Hox1, Hox2, Hox3 ordering of the Hox genes in the chromosome from the telomeric to the centromeric side of the chromosome. Surprisingly, the spatial activation of the Hox genes in the ontogenetic units of the embryo follows the same ordering along the anterior-posterior embryonic axis. The chromosome microscale differs from the embryo macroscale by 3 to 4 orders of magnitude. The traditional biomolecular mechanisms are not adequate to comprise phenomena at so divergent spatial domains. A Biophysical Model of physical forces was proposed which can bridge the intermediate space and explain the results of genetic engineering experiments. Recent progress in constructing instruments and achieving high resolution imaging (e.g., 3D DNA FISH, STORM etc.) enable the assessment of the geometric structure of the chromatin during the different phases of Hox gene activation. It is found that the mouse HoxD gene cluster is elongated up to 5–6 times during Hox gene transcription. These unexpected findings agree with the BM predictions. It is now possible to measure several physical quantities inside the nucleus during Hox gene activation. New experiments are proposed to test further this model.

Possible rules for the ancestral origin of Hox gene collinearity

The Hox gene cluster is believed to have formed from a single ProtoHox gene by repeated cycles of the following events: tandem gene duplication, mutation to generate a new expression boundary along the embryonic axis, and acquisition of a new Hox patterning function. The Hox cluster in Bilateria evolved in compliance with the so-called collinearity rule. That is, the order of the genes along the chromosome corresponds with the order of their embryonic expression domains along the head-tail axis. Gaunt (2015) suggested that collinearity may have arisen as a mechanism to minimise the incidence of boundaries between active and inactive genes within the Hox cluster. We now attempt to clarify the model by presenting it in the form of three rules: 1) no two Hox genes may persist in the same cluster with the same anterior boundary of activity in the same tissue; 2) an inactive Hox gene must not be flanked by two active Hox genes; 3) an active Hox gene must not be flanked by two inactive genes. We provide evidence and illustrative computer simulations to show that these rules, which can apply only to partially overlapping patterns of Hox activity, may account for the ancestral origin of Hox gene collinearity.

Hox Gene Collinearity: From A-P Patterning to Radially Symmetric Animals.

Hox gene collinearity relates the gene order of the Hox cluster in the chromosome (telomeric to centromeric end) with the serial activation of these genes in the ontogenetic units along the Anterior-Posterior embryonic axis. Although this collinearity property is well respected in bilaterians (e.g. vertebrates), it is violated in other animals. The A-P axis is established in the early embryo of the sea urchin. Subsequently, rotational symmetry is superimposed when the vestibula larva is formed. In analogy to the linear A-P case, it is here hypothesized that the circular topology of the ontogenetic modules is associated to the architectural restructuring of the Hox loci where the two discrete ends of the Hox cluster approach each other so that an almost circular DNA contour is created. In the evolutionary process the circular mode undergoes double strand breaks and the generated cluster ends are attached to the open ends of the flanking chromosome. This event may lead to a novel gene ordering associated with an evolutionary innovation. For example, the loss of Hox4 is followed by the formation of a shorter gene circular arrangement. The opening of this contour at the missing Hox4 location and its connection to the chromosomal flanking ends leads to a new diversification namely the creation of the unusual gene order of the sea urchin Hox cluster.

The significance of Hox gene collinearity.

Arthropods and vertebrates inherited their Hox clusters from an ancestral cluster of at least six genes already present in their last common ancestor, Urbilateria. Clustering and a common transcriptional direction are both likely features of the way that the gene complex first arose in a process of tandem gene duplication. Spatial collinearity (correspondence between ordering of Hox genes along the chromosome and their expression patterns along the head-tail axis) has been conserved in many animal groups and is likely to have been already present in Urbilateria. It is not known why the Hox cluster evolved with spatial collinearity. Four models are discussed. These vary in the significance they place upon Hox chromatin structure, and also on whether they propose that collinearity is primarily concerned with establishment or maintenance of Hox expression. Published proposals to explain spatial collinearity, which invoke enhancer sharing, chromatin closing or chromatin opening, are either problematic or can offer only partial explanations. In an alternative proposal it is suggested here that spatial collinearity evolved principally to maximise physical segregation, and thereby minimise incidence of boundaries, between active and inactive genes within the Hox cluster. This is to minimise erroneous transfer of transcriptional activity, or inactivity, between adjacent Hox genes.

Evolutionary Constraints Favor a Biophysical Model Explaining Hox Gene Collinearity

The Hox gene collinearity enigma has often been approached using models based on biomolecular mechanisms. The biophysical model is an alternative approach based on the hypothesis that collinearity is caused by physical forces pulling the Hox genes from a territory where they are inactive to a distinct spatial domain where they are activated in a step by step manner. Such Hox gene translocations have recently been observed in support of the biophysical model. Genetic engineering experiments, performed on embryonic mice, gave rise to several unexpected mutant expressions that the biomolecular models cannot predict. On the contrary, the biophysical model offers convincing explanation. Evolutionary constraints consolidate the Hox clusters and as a result, denser and well organized clusters may create more efficient physical forces and a more emphatic manifestation of gene collinearity. This is demonstrated by stochastic modeling with white noise perturbing the expression of Hox genes. As study cases the genomes of mouse and amphioxus are used. The results support the working hypothesis that vertebrates have adopted their comparably more compact Hox clustering as a tool needed to develop more complex body structures. Several experiments are proposed in order to test further the physical forces hypothesis.

Comparison of Models for the Collinearity of Hox Genes in the Developmental Axes of Vertebrates

Hox gene clusters are very frequent in many animal genomes and their role in development is pivotal. Particularly in vertebrates, intensive efforts have established several properties of Hox clusters. The collinearity of Hox gene expressions (spatial, temporal and quantitative) is a common feature of the vertebrates. During the last decade, genetic engineering experiments have revealed some important facets of collinearity during limb and trunk development in mice. Two models have been proposed to explain all these properties. On one hand the ‘two-phases model’ makes use of the molecular regulatory mechanisms acting on the Hox genes. On the other hand, the’biophysical model’ is based on the signals transduced inside the cell nucleus and the generation of forces which apply on the cluster and lead to a coordinated activation of Hox genes. The two models differ fundamentally and a critical and detailed comparison is presented. Furthermore, experiments are proposed for which the two models provide divergent predictions. The outcome of these experiments will help to decide which of the two models is valid (if any).

Cofactor Binding Evokes Latent Differences in DNA Binding Specificity between Hox Proteins

Members of transcription factor families typically have similar DNA binding specificities yet execute unique functions in vivo. Transcription factors often bind DNA as multiprotein complexes, raising the possibility that complex formation might modify their DNA binding specificities. To test this hypothesis, we developed an experimental and computational platform, SELEX-seq, that can be used to determine the relative affinities to any DNA sequence for any transcription factor complex. Applying this method to all eight Drosophila Hox proteins, we show that they obtain novel recognition properties when they bind DNA with the dimeric cofactor Extradenticle-Homothorax (Exd). Exd-Hox specificities group into three main classes that obey Hox gene collinearity rules and DNA structure predictions suggest that anterior and posterior Hox proteins prefer DNA sequences with distinct minor groove topographies. Together, these data suggest that emergent DNA recognition properties revealed by interactions with cofactors contribute to transcription factor specificities in vivo.

Physical forces may cause Hox gene collinearity in the primary and secondary axes of the developing vertebrates

The features of spatial and temporal Hox gene collinearity along the anteroposterior and secondary axes of vertebrate development have been extensively studied. However, the understanding of these features remains problematic. Some genetic engineering experiments were performed and the consequent modifications of the Hoxd gene expressions in the vertebrate limb and trunk were presented. A two-phases model was proposed to describe the above results but still many data cannot be explained. In the present work a different mechanism is put forward in order to deal with the above experiments. This alternative mechanism (coined biophysical model), is based on the hypothesis that physical forces decondense and 'loop out' the chromatin fiber causing the observed Hox gene collinearity phenomena at the early stages of axonal development. The two models are compared in detail. The biophysical model adequately explains the data even in cases where the results are characterized as unexpected. Furthermore, the biophysical model predicts that the Hox gene expressions are entangled in space and time and this coupling is compatible with the data of the early developmental stages. Additional experiments are proposed for a direct test of this model.

Regulatory constraints in the evolution of the tetrapod limb anterior–posterior polarity

The anterior to posterior (A-P) polarity of the tetrapod limb is determined by the confined expression of Sonic hedgehog (Shh) at the posterior margin of developing early limb buds, under the control of HOX proteins encoded by gene members of both the HoxA and HoxD clusters. Here, we use a set of partial deletions to show that only the last four Hox paralogy groups can elicit this response: that is, precisely those genes whose expression is excluded from most anterior limb bud cells owing to their collinear transcriptional activation. We propose that the limb A-P polarity is produced as a collateral effect of Hox gene collinearity, a process highly constrained by its crucial importance during trunk development. In this view, the co-option of the trunk collinear mechanism, along with the emergence of limbs, imposed an A-P polarity to these structures as the most parsimonious solution. This in turn further contributed to stabilize the architecture and operational mode of this genetic system.