1. Introduction
In 1978 E.B. Lewis discovered a fundamental
property of Developmental Biological: Hox Gene Collinearity (HGC).
HGC is the long- range interaction of entities where the archetypical example
is the Coulomb force relating material objects at different geometric ranges.
In particular, HGC correlates genes (entities inside the cell nucleus) with
embryonic units. These entities are quite apart from each other.
Hox genes play an important role in the
development of most animals and plants. Some Hox genes form clusters which are
crucial for the Embryogenesis of Metazoa. The importance of this clustering was
noticed by Lewis who studied the genetics of Drosophila [
1]. He observed that some genes of the genome
(later coined Hox genes) were located in order along the telomeric to
centromeric direction and were denoted (Hox1, Hox2, Hox3…). Lewis noticed that
the genes of these clusters were expressed in the same order along the
anterior-posterior axis of the Drosophila embryo (
Figure 1) [
1].
This is an astonishing event since this correlation occurs between extremely
distant domains - in one hand the gene sequence in the cell nucleus and in the
other the Drosophila embryo. Biomolecular interactions
alone cannot create such correlations [
2].
This surprising phenomenon was named
Spatial Collinearity (SC). Some
years later, another collinearity was observed particularly in the Vertebrates:
Temporal Collinearity (TC). According to TC, the first Hox gene (Hox1)
of the Hox cluster starts being expressed. Later, Hox2 is expressed and even
later Hox3 followed until all Hox genes are expressed following the sequence
Hox1, Hox2, Hox3… [
3].
In order to explain these phenomena, I proposed in
2001 a
Biophysical Model (BM) according to which, pulling physical
forces could explain the experimental findings [
2,
4] (
Figure 1a, b). Several predictions of BM
were compared to many collected data and the comparison confirmed BM [
5,
6].
A simple heuristic expression for these pulling
forces F was proposed [
7,
8].
In Eq.(1), F is a simplification of the pulling Coulomb force since the distance between the electric charges (N and P) is missing. This ‘quasi Coulomb force’ F is applied at the telomeric end of the Hox cluster. It turns out that this arbitrary omission of the Coulomb range reflects a deep connection of these ‘quasi- Coulomb’ forces to the fundamental phenomenon of
Symmetry (see
Section 2.1 and
Section 2.2). Ν and P stand for the negative and positive electric charges acting on a Hox cluster. In the above heuristic formulation, the Hox cluster consists of a deployed finite sequence of Hox genes along the telomeric to centromeric ends of the cluster (Hox1, Hox2, Hox3…). The numbers assign the
gene order in the cluster. These numbers determine the
order membership to the
Paralogy Group (PG). (Here is followed Duboule’s definition of PG [
9]).
As mentioned above, (N) represents the microscopic contribution to F and it is a real entity – the negative electric charge of the DNA sequence in the cell nucleus. (P) represents a positively charged molecular structure located opposite the telomeric end of the Hox cluster (
Figure 1). Contrary to N, P is a fictitious entity as yet, standing for the embryonic-macroscopic contribution to F. Note that the known morphogens of the present time like Sonic Hedgehog, Fibroblast Growth Factors, Retinoic Acid and the plethora of other morphogenetic factors were fictitious fifty years ago. The existence of P does not contradict any First Principle so it is legitimate to anticipate its existence as advocated in [
4]. F pulls the Hox genes sequentially out of the cluster (
Figure 1). As mentioned already, Eq.(1) is a heuristic expression that was successfully tested in several experiments [
5,
6,
7,
8].
Hox genes control the normal development of animals (wild type). Spontaneous mutation of these genes cause severe malformations (Homeosis), consisting of parts of the animal growing in the wrong location of the body. In Homeosis PG ordering is violated.
About twenty years ago, an important advancement was achieved concerning the transfer of specific molecules from outside the cell into the inner domain of the cell nucleus [
10,
11,
12,
13]. For example, it was noticed that significant amounts of Activin are gathered outside the cell nucleus. Controlled amounts of this activin were
transduced inside the nucleus causing specific modifications on the genome. It is assumed that BM combined with the action of transduction technology can affect Hox gene expression. This possibility is incorporated in the present Review.
BM is based on the hypothesis that pulling forces are applied at the telomeric end of the Hox cluster. This hypothesis was elaborated in detail and it was concluded that the cluster is elongated along the direction of the force [
14]. This BM prediction wash later confirmed [
15,
16,
17]. In some cases the measured elongation of the activated Hox cluster was five times longer than the length of the inactive Hox cluster [
16]. When Hox cluster activation is initiated, a weak force (F1) pulls the first gene of the cluster (Hox1) out of its niche toward the interchromosome domain (ICD) (
Figure 1(a), (b)).
Particularly Hox1 is directed towards the transcription factory domain (TFD) where Hox gene activation (expression) is possible [
18,
19]. The pulling forces increase irreversibly: i) increasing along the gene location from Hox1 towards Hox13 (Hox1 → Hox13), and ii) following the time course t1 → t3. Under force F2 (F2 > F1) Hox1 and Hox2 are extruded from their niche. This process continues until all Hox genes of a
complete Hox cluster are transferred in the TFD. In the mechanical analogue, for the efficient function of an elongated elastic spring, besides the pulling force at one of its ends, a proper fastening should be applied at its other end. Accordingly, the Hox cluster should be fastened at the centromeric end of the cluster (
Figure 1(b)).
The Vertebrate Hox clusters comprise four homologue clusters
(HoxA, HoxB, HoxC, and HoxD) as shown in
Figure 1(c) [
11]. Each homologue cluster is included in a separate chromosome. In these homologue clusters the PG identity is conserved. However, in the course of Evolution, modifications of the mouse genes are possible up to the point of gene deletion. This ordered mouse Hox clusters remind of a ‘ratchet’ allowing motion in one irreversible direction only [
8]. Note that some ‘teeth’ of the ratchet may be missing (corresponding to Hox genes deleted in the course of Evolution).
Besides Vertebrates, the contemporary cephalochordate Amphioxus is a descendant of the ancestor Amphioxus. This ancestor Amphioxus was the ancestor of both Drosophila and vertebrates [
9]. Amphioxus lived after the Cambrian period of evolutionary explosion 500 million years ago (Mya). Vertebrates and Drosophila appeared a few Mya later. Amphioxus has 14 Hox genes whereas vertebrates and Drosophila have 13 (Hox14 is missing).
3. Complete vs Split Hox Clusters
Hox Gene Collinearity (both Spatial and Temporal) has been unequivocally confirmed in the Vertebrates. However, in recent years it was found that this is not true in many other animal species, particularly in invertebrates. For instance, it was observed that Hox collinearity is violated in the lophotrochozoa and this violation was associated with the brachiopods whose Hox cluster is broken [
26,
27]. In brachiopods both spatial and temporal collinearities are violated, while lophotrochozoan morphological novelties result from Hox Collinearity violation [
26]. It was argued above that for the insertion of a circularly organized Hox cluster in the flanking genome, a break (split) of the cluster is necessary. It is clear that Hox cluster splitting is a necessary step for evolutionary novelties.
It has been emphasized that tight Hox clustering is lost during Evolution [
28,
29,
30,
31]. More specifically D. Ferrier and P. Holland observed that Hox clusters are constrained by TC in their gene order [
28,
29]. Moreover,
complete Hox clusters are associated with the spectrum of Hox gene expressions along the
whole Anterior – Posterior embryonic axis. Otherwise, if the gene expression does not extend along the whole A-P axis, the Hox gene cluster is split. In this cluster splitting tendency, the Hox clusters ‘may fall apart when TC is no more needed for [
28,
29]. Similar arguments were put forward by several other authors before and after the above observation [
29,
30].
Drosophila has a typically split Hox cluster. The Drosophila Hox cluster has 13 genes consisting of two subclusters ANT
-C and
BX-C depicted in
Figure 3 [
29].
Drosophila, together with the vertebrate Hox cluster, originates from a large ancestral Hox cluster. Cloning had later identified Amphioxus as the common ancestor of insects and vertebrates [
28], and a one-to-one correspondence between the Amphioxus Hox genes and the Drosophila Hox genes was confirmed [
28,
29]. However in this correspondence some Drosophila Hox genes of the ANT-C subcluster developed novel evolutionary non-Hox functions. For instance, the Drosophila complex of Hox genes (zen1, zen2, bcd) corresponds to the ancestral Amphioxus Hox3 gene. Some Central genes have evolved from tandem evolutionary duplications [
29]. (See in [
29]
Figure 1). The BX-C subcluster consists of the 3 last genes (Ubx, Abd-A, Abd-B) of the Drosophila Hox cluster. The summarizing conclusion from the above analysis is that TC is responsible for a
complete Hox cluster. If this is not possible (or not needed) the Hox cluster is split [
28]. In any case SC is a necessity for a Hox cluster.
As stressed above, if TC extends to a fraction only of the Hox cluster range, the cluster is expected to split. In a way, Hox expression in a range between Anterior and Posterior ends (
‘space’) is translated into
‘time)’- where TC coordinates a ‘Hox clock’- more exactly a ‘Hox timer’ since the time course is irreversible [
30].
5. Complex Expression Patterns
According to BM, it is expected that complex patterns can be created by splitting the early (simply connected) Hox gene expression into expression domains separated by a ‘ditch’ zone (cf.
Figure 4a) [
37,
38]. This splitting was already confirmed in 2013 [
37]. Following this line of thought [
37,
38], BM predicts that a small DNA strip containing (Hox10 and Hox11) has a strange expression behavior in time (
Figure 4). These genes can be pushed in and out of the Hox cluster activation domain which is depicted in the dark blue `circle (
Figure 4b). In the left graph, Hox10, Hox11, Hox12 are activated. In the middle graph, (Hox10, Hox11) are pushed out of the activation domain. In the right graph, (Hox10, Hox11) reenter later in the activation domain following the increased force of BM in the time course [
37,
38]
.
Recently it was unexpectedly observed in limb digit condensations a local expression disappearance in the interdigital area as shown in Figures
{ 2 (m, o) and
2 (n, p)} of [
39]. These figures are reminiscent of the (Hox10,Hox11) disappearance in
Figure 4b of the present Reviewed. The above reminiscence is a further confirmation of the BM pulling forces.
It is interesting that the theoretical prediction ‘Biophysics
precedes Biochemistry’ [
38] was experimentally confirmed soon after ‘…structural organization of HoxD cluster may
predate transcriptional activation’ [
16].