A biophysical mechanism may control the collinearity of Hoxdgenes during the early phase of limb development
© Henry Stewart Publications 2009
Received: 17 November 2008
Accepted: 17 November 2008
Published: 1 April 2009
A biophysical model has been proposed which deals with the observed collinearity of Hox gene expressions in developing vertebrate limbs. It is assumed that physical forces gradually dislocate the genes of the Hoxd cluster from inside the chromosome territory into the interchromosome domain, where the genes are activated. In particular, the action of Coulomb electric forces has been estimated in detail. Genetic engineering experiments (deletions, duplications and transpositions) were recently reported for Hoxd expression during limb development. Here, we analyse these results and show that the biophysical model explains them successfully.
The collinearity of Hox gene expressions is a surprising feature observed during animal development along the primary and secondary body axes [1, 2]. These gene expressions follow in space and time the sequence of gene ordering from 3' to 5' in the Hox cluster [2, 3]. In the Hoxd cluster, the genes are lined up in the 3' to 5' direction and their expression domains (Hoxd1, Hoxd2, Hoxd3, ...) along the anteroposterior axis follow a peculiar pattern: the anterior border of expression for each gene is shifted in the anterior-posterior direction -- for example, the anterior border of Hoxd1 expression is more anterior than the anterior border of Hoxd2 expression etc. Thus, a pattern of partially overlapping expression domains along the head to tail axis is formed ('spatial collinearity') . Furthermore, the genes of the Hoxd cluster are activated one after the other in the same 3' to 5' order ('temporal collinearity') [2, 3]. Tarchini and Duboule  have put forward a model in which, during early limb budding (first wave), two regulatory influences control the activation of Hoxd genes: one mechanism is located at the telomeric (3')side of the cluster (early limb control regulation; ELCR) and the other relies on sequences centromeric to the cluster (POST). Thus, in order to describe the data in terms of biomolecular processes only, two opposing collinearity mechanisms have to be introduced .
The biophysical model -- Comparison with experiments
The biophysical model does not merely describe Hox gene collinearity. Its aim is to explain it and to provide the mechanism that causes the modes of Hox gene expression by integrating and coordinat-ing actions at both the multicellular and subcellular levels. Of course, the control regions, such as ELCR and any other regulatory elements that have yet to be discovered, are essential tools of the biochemical machinery that participate in the realisation of Hox gene expression. We believe, however, that this essential biomolecular pathway is triggered and controlled by a more fundamental, single and unidirectional mechanism, such as the one presented here. The biophysical model reproduces spatial, temporal and quantitative collinearities well .
Gene deletions, duplications and transpositions
Duboule and collaborators [4, 6, 7] have performed detailed and extensive experiments of Hoxd gene deletions and duplications in order to explore the features of collinearity during vertebrate limb development. Their findings for the early phase of limb development  are highly illuminating, and it is worth comparing them with the biophysical model predictions.
Suppose the subcluster of the 5' Hoxd genes has an evenly distributed negative charge N (N = 10 in arbitrary units) and that every gene region (Hoxd 9, ... Hoxd13) carries a charge N/5 = 2. Before activation of the 5' Hoxd subcluster, all genes are embedded inside the CT. In this case, we can set P = 0 and the Coulomb force is 0. For activation of the first gene (Hoxd 9), we set P = 2 (in arbitrary units) and P = 10 when all genes (Hoxd9, ..., Hoxd13) are extruded in the ICD.
In this truncated Coulomb equation, we neglect the dependence of F on the distance between positive and negative charges, and the Coulomb's constant is set at 1 (in appropriate units). The Coulomb force causes gene fibre decondensation and extrusion of length L. We assume, in a first approximation, that this extrusion length L is proportional to F in the spirit of Hooke's law of elastic spring expansion. Deletion of one Hoxd gene region causes a reduction in the negative charge of the cluster; therefore, the Coulomb force will decrease and consequently the extrusion length of the Hoxd fibre will be shortened. The opposite will occur with a duplication of one Hoxd gene region in the cluster: the larger the duplication, the longer the gene fibre extrusion. In the supplementary material, we present some indicative experiments of gene deletions, duplications and transpositions  and compare them with the biophysical model predictions. The agreement is astonishing.
In an analysis of the Hoxd locus in mouse development, activation along the primary anteroposterior axis is accompanied by fibre extrusion from the CT . In the limb bud, at the same stage (E 9.5), however, no looping out is detected, even though the locus is still decondensed. A further confirmation of this difference is needed, taking into account that at that stage Hoxd expressions (eg Hoxd4) in the mouse limb bud are much weaker than in the tailbud .
The biophysical model explains how collinearity phenomena occur during the early stages of limb development. At later stages, when all Hoxd genes are expressed, secondary regulations and gene interactions modulate their expressions and the complex landscape of digits starts emerging. Recently, Montavon et al.  analysed experimentally and theoretically the features of quantitative collinearity in developing digits. They used some of their essential findings to formulate a model for the differential importance of gene rank -- the relative position of Hoxd genes with respect to the 5' extremity of the cluster. For instance, they observed that when Hoxd13 (the most strongly expressed gene) was deleted, the gene that substituted for it in the last 5' position was expressed with an increased intensity. It is instructive to view this result in the framework of the biophysical model: as soon as Hoxd13 is activated, the whole gene fibre is extruded and Hoxd13 is located closest to the surface of the CT, and its expression is therefore strongest. (This occurs in the posterior domain of the limb bud.) When Hoxd13 is deleted, the extruding Coulomb force F decreases and the extruded fibre retreats toward the interior of the CT. The consecutive gene (eg Hoxd12) is dislocated towards the CT surface and consequently its expression is strengthened. This is the essence of the observation by Montavon et al.  and it indicates, on the one hand, that this observation constitutes additional evidence in favour of the biophysical model. On the other hand, it shows that our interpretation of quantitative collinearity can be extended to all stages of limb development. Finally, the above authors have managed to establish a quasi-linear relationship between the total transcriptional activity at the locus and the number of transcription units -- e.g. active genes in developing digits (Figure 4b in reference ). This remarkable correlation supports our hypothesis of Coulomb force variations and the corresponding proportional changes in the extruded fibre lengths (reminiscent of Hooke's linear law of elastic expansions).
For cells expressing only one gene (Hoxd 9):
P = 2, N = 10, F =2 × 10 = 20, L = 1 (Figure 1a), Domain D9
For cells expressing two genes (Hoxd9 and Hoxd10):
P = 4, N = 10, F =4 × 10 = 40, L = 2 (Figure 1b), Domain D10 < Domain D9
For cells expressing three genes (Hoxd9, Hoxd10 and HoxdH):
P = 6, N = 10, F = 6 × 10 = 60, L = 3, Domain D11 < Domain D10 < Domain D9...etc.
In experiments where one gene of the subcluster is deleted, the negative charge of the cluster will decrease to N = 8 and the length of the extruded fibre will be shorter. Two gene deletions will produce an even bigger reduction in the length of the extruded fibre. By contrast, when one gene is duplicated, N will increase to 12 and the extruded fibre will be longer etc.
In the following, we estimate how deletions and duplications modify gene expressions, and particularly in domains D9, D10 and D11.
(a) Posterior gene deletions
For the complete extrusion of Hoxd9, a stronger force is needed and this is achieved by a posterior shift to cells where the morphogen value is higher, as is the corresponding P value: P (++ → +++).
If more than one posterior gene is deleted, the posterior shift of Hoxd9 activation will increase. Tarchini and Duboule  designated this as an unexpected observation. For the biophysical model, this posteriorisation is an obvious outcome.
(b) Anterior gene deletions
For activation of Hoxd10, F can be reduced to F1, and this occurs at more anterior cell positions with alower P value: P (++++ → +++). An ectopic expression of Hoxd10 in the anterior region of the limb bud will be observed, associated with a premature expression of the gene. This was verified by Tarchini and Duboule .
(c) Posterior gene duplications
Compared to gene deletions, gene duplications have an opposite effect on Hoxd gene expressions. Consider, again, cells in which Hoxd 9 and Hoxd10 are activated, where the pulling force is (at least) F2 = 40. If posterior gene Hoxd11 is duplicated (Figure 2Sa (Figure 2)), the charge of the cluster will increase to N = 12.
The pulling force on the fibre will become: F = 4 × 12 = 48, which is higher than F2. Activation of Hoxd10, however, can start with F2; therefore, F can be reduced to this value. This is achieved by an anterior shift of D10; in this case, the morphogen decreases, and so does P (++++ → +++) and F → F2. The result is an anteriorisation and premature expression of Hoxd10.
Let us examine Hoxd10 expression at a position and time without Hoxd11 duplication. At this position, after Hoxd11 duplication, Hoxd10 will be abnormally translocated further inside the ICD, away from the CT border. A downregulation of Hoxd10 expression is expected at this position (quantitative collinearity). This was verified by Kmita et al. 
(d) Anterior gene duplications
Consider the case of two activated genes (Hoxd9 and Hoxd10), where F2 = 40. If anterior gene Hoxd 9 is duplicated, the pulling force on the fibre increases from F2 to F = 48. For activation of Hoxd10, however, this increase is not sufficient, since now three genes must be extruded. The corresponding pulling force for three genes is F3 = 60 (Figure S2b (Figure 3)). For a stronger force F3, P must take a higher value: P (+ + + → + + + + +). Therefore, after an anterior gene duplication, the expression of Hoxd10 will be shifted posteriorly, in agreement with the findings of Tarchini and Duboule .
These calculations can be applied to extended gene deletions and duplications incorporating several genes from the entire Hoxd cluster. Following the above rules, the estimation of anteriorisations or posteriorisations is straightforward, and the results agree with the findings of Duboule and co-workers. The conclusion is that Hoxd gene deletions and duplications strongly support the present electric force mechanism of the biophysical model.
Hoxbl transposition in Hoxd
Recently, the Hoxb1 transgene was transposed at the 5' end of the Hoxd cluster . The results from this transposition are in agreement with the biophysical model. For example, in the wild-type mouse embryo at stage E7.5, Hoxd13 is not expressed in the primitive streak tissue. By contrast, Hoxd13 is activated when the Hoxb1-LacZ reporter is inserted at the 5' end of the Hoxd cluster. This finding is explained by the biophysical model, since, in analogy with the Hoxd regions, the Hoxb1 transgene carries a negative charge and the total N increases in the transgene embryo (as analysed in detail above). The pulling force on the gene fibre therefore increases, so that Hoxd13 loops out of the CT and its activation is possible. In the wild-type limb bud at E9.5, the Hoxd region, although decondensed, does not extrude out of the CT . In the Hoxb1-LacZ embryo, however, Hoxd does loop out and this can be attributed, again, to the increased force of attraction that pulls the gene fibre out of the CT.
- Lewis EB: A gene complex controlling segmentation in Drosophila. Nature. 1978, 276: 565-570. 10.1038/276565a0.View ArticlePubMedGoogle Scholar
- Nolte C, Krumlauf R: Expression of Hox genes in the nervous system of vertebrates. HOX Gene Expression. Edited by: Papageorgiou, S. 2007, Landes Bioscience and Springer Science+Business Media, New York, NY USA.Google Scholar
- Nelson CE, Morgan BA, Burke AC, et al: Analysis of Hox gene expression in the chick limb bud. Development. 1996, 122: 1449-1466.PubMedGoogle Scholar
- Tarchini B, Duboule D: Control of Hoxd genes' collinearity during early limb development'. Developmental Cell. 2006, 10: 93-103. 10.1016/j.devcel.2005.11.014.View ArticlePubMedGoogle Scholar
- Papageorgiou S: Pulling forces acting on Hox gene clusters cause expression collinearity. Int J Dev Biol. 2006, 50: 301-308. 10.1387/ijdb.052034sp.View ArticlePubMedGoogle Scholar
- Kmita M, Fraudeau N, Herault Y, Duboule D: Serial deletions and duplications suggest a mechanism for the collinearity of Hoxd genes in limbs. Nature. 2002, 420: 145-150. 10.1038/nature01189.View ArticlePubMedGoogle Scholar
- Montavon T Le, Garrec J-F, Kerszberg M, Duboule D: Modeling Hox gene regulation in digits reverse collinearity and the molecular origin ofthumbness. Genes Dev. 2008, 22: 346-359. 10.1101/gad.1631708.PubMed CentralView ArticlePubMedGoogle Scholar
- Morey C, Da Silva NR, Kmita M, et al: Ectopic nuclear reorganization driven by a Hoxb1 transgene transposed into Hoxd. J Cell Sci. 2008, 121: 571-577. 10.1242/jcs.023234.PubMed CentralView ArticlePubMedGoogle Scholar
- Morey C, Da Silva NR, Perry IP, Bickmore WA: Nuclear reorganization and chromatin decondensation are conserved, but distinct, mechanisms linked to Hox gene activation. Development. 2007, 134: 909-919. 10.1242/dev.02779.View ArticlePubMedGoogle Scholar