1. Introduction
The Golden Ratio
was recognized at least as early as 500 BCE by Phidias, after whom the symbol
remains named [
1]. Naturalists over the centuries have commented on its presence in plants, mollusks, and vertebrates, it has been depicted in the arts, and it has been the subject of teleological conjecture [
2,
3,
4]. Multicellular organisms evolved from unicellular organisms in the Precambrian period 580 million years ago, and appear to have done so on more than one occasion [
5]. Fibonacci forms are found among these organisms [
6], and are found among primitive plants such as algae [
7]. Nevertheless, the molecular or cellular mechanisms for its presence in multicellular organisms remain unknown. It plays a role in plant phyllotaxis and perhaps in cellular packing, but a broader adaptive value, if any, remains unknown [
8,
9,
10,
11]. We investigate Golden Ratio mathematics and its implications regarding the ontogeny of multicellular organisms.
The Golden Ratio is real-valued while cell population counts are whole-numbered. The Fibonacci numbers are the particular whole numbers that obey the recurrence relation
, with
and
, giving the sequence
. Binet’s formula offers a way to connect the Golden Ratio to the Fibonacci numbers,
Like the Golden Ratio, forms of the Fibonacci numbers appear in ancient culture [
12,
13].
In this paper, we shift the Fibonacci index by one, with
. With
and
, this gives the sequence
. Under this convention, if we start with one cell and assume the cell population grows as the Fibonacci numbers, then
gives the population cell count at the
generation. Furthermore, this convention offers a combinatorial interpretation of
, giving the number of ways that the number
n can be composed as sums of 1 and 2 [
14]. This is seen in
Section 7.
Single-celled and multicellular organisms both perform replication. Some single-celled organisms exhibit adhesion. For example, the bacterium Streptococcus pneumonia may adhere into pairs and chains, and Staphylococcus aureus may adhere into clusters. We hypothesize that multicellular organisms differ from single-celled organisms by having replication and adhesion be mutually regulated phenomena. An example in multicellular organisms is contact inhibition, where the adhesion between two cells yields the coordinated replication of only one of them [
15]. The inter-regulation of adhesion and replication enables more extensive structures having cellular specialization amidst topological variations.
1.1. Replication
A cell in a multicellular organism has a cell cycle divided broadly into two intervals, interphase and mitosis, as shown in
Figure 1. Interphase includes the time interval where a cell makes its specialized contribution to an organism and is where it makes preparatory changes so that it may divide. The proportions of time spent in interphase and mitosis vary greatly, but generally, a cell spends much more time in interphase.
The first phase in interphase is (G as in gap). In a cell increases its supply of proteins, organelles (such as mitochondria, which catalyze the conversion of oxygen into energy), and it grows in size. It may transition to and from a non-growing phase .
Chromosomal DNA is replicated in the S interval of interphase. After DNA replication, each chromosome has a pair with an identical DNA sequence. is a growth phase after DNA replication.
Mitosis (M) is a brief interval of the cell cycle where the cell divides into two. It consists of prophase where the chromosomes become dense, metaphase where the chromosome pairs line up on the equatorial plane, anaphase where the chromosome pairs detach into sister chromatids, telophase where the sister chromatids move apart and form separate nuclei, and abscission where the membranes cleave to form 2 progeny cells.
Cytokinesis is the interval of mitosis between anaphase and telophase when the sister chromatids segregate. They are pulled apart by mitotic spindles composed of tubulin. Cytokinesis provides the proper distribution of genetic material and cytoplasm between the two progeny cells. It involves the formation of a contractile ring, composed of actin and myosin filaments, at the equatorial plane of the dividing cell. The contractile ring is oriented perpendicular to the spindles. The contractile ring contracts, leading to the formation of a cleavage furrow and eventually to abscission, the detachment of the two progeny cells. The relationships between the contractile ring and the spindle at abscission can encode polarities in the mitotic progeny [
16,
17,
18,
19,
20].
1.2. Adhesion
Adhesion is a persistent cell-to-cell connection with binding molecules that maintain physical proximity between the adhered cells. It may occur by direct cell-cell contact or be mediated by an intervening extracellular matrix. It may offer modes of communication or coordination between the cells.
There are broadly five types of molecular cell-cell adhesion by direct contact, illustrated in
Figure 2. A cell’s intracellular content is separated from the extracellular content by a bilipid membrane. With two of the adhesion types, gap junctions and tight junctions, lipid cell membranes between adjacent cells are intermingled. These adhesion types allow for the exchange of molecules between intracellular spaces. For the other three adhesion types, adherens junctions, desmosomes, and focal adhesions (having tunneling nanotubes) the adhesion is mediated by a multi-subunit protein structure that spans the extra to the intracellular space [
21]. The binding to the outer surface cell junctions generally induces metabolic changes within the cell. There is evidence that adhered cells can coordinate their cell cycles [
22,
23,
24]. Furthermore, the tunneling nanotubes of focal adhesions enable intercellular exchange of DNA and RNA [
21].
On the internal cell surface, the adhesion protein complex commonly binds to elements of the cytoskeleton. There are broadly three types of cytoskeleton protein. These are actin filaments, microtubules, and intermediate filaments such as vimentin [
25]. Of these, actin filaments and microtubules are capable of active contraction. They play roles in modifying cell shape and in cell migration. The cytoskeleton binds the inner surface of cell adhesion sites to intracellular organelles, including the cell nucleus. In mitosis, microtubules bind to chromosome pairs and pull them apart as a cell replicates. The cytoskeleton offers intracellular binding between the inner surface of cell adhesion sites [
26].
The five types of cell-cell adhesions have differences in their properties, but the property they have in common is that they physically bind adjacent cells together and may offer signaling channels. This is the property of cell adhesion that we focus on, and we use adhesion as an umbrella term, but we do not treat the five cell-cell adhesion types separately. The extracellular matrix may offer indirect forms of cell to cell communication between cells not sharing a membrane to membrane adhesion. For example, molecules of the extracellular matrix may be anisotropic, and this anisotropy may be reflected intracellularly in the cellular cytoskeleton [
27,
28,
29].
3. Candidate Molecular Mechanisms
An assignment at adhesion to mitose or not is consistent with laboratory evidence that the five categories of cell adhesion (
Figure 2) participate in gene and cell cycle regulation [
50,
51,
52,
53,
54,
55]. If these adhesion categories are able, in particular, to establish semaphore cell cycle coordination so that one of a pair undergoes mitosis, then Fibonacci population counts might be observed. If the paired cells are identical and the adhesion communication are symmetric, then it may be equal odds as to which cell of the adhered pair undergoes mitosis.
The molecular data are consistent with mutually reinforcing roles for mitosis and adhesion [
15]. This might be the case if there is polarity among mitotic progeny that influences the mitosis semaphore at adhesion. The polarity-dependent mitotic behavior could span generations. A candidate molecular mechanism might be a cytoskeletal polarity that is linked across the cell membrane to polarities of extracellular adhesion and the extracellular matrix [
27,
28,
29,
56,
57]. For example, the midbody, which forms along the cleavage line at telophase and can be asymmetrically inherited by a progeny [
58], is involved in cell proliferation [
59], and has attachment to cell surface adhesion structures [
54].
There are other opportunities for stereochemical variations between mitotic progeny with extracellular linkage. When chromosomes condense in the cell nucleus in prophase, mitotic spindles form between centrosomes at opposite poles of the cell and attach to the kinetochores of the chromosomes (
Figure 1). The spindles consist of microtubules which are motoric protein structures. They have a left chiral orientation as they link centrosome to centrosome at opposite cell poles. They attach to the kinetochores of chromosomes. Their chirality causes them to apply torque to chromosomes when they contract [
60]. The helical spindle attachment to chromosomal kinetochores opens opportunities for stereoisomeric distinction upon chromosomal separation [
61]
Figure 1). This opens an opportunity for enantiomeric chromosomal differences after metaphase that remain coordinated with extracellular adhesion patterns.
4. Simplicial Steps by Replication and Adhesion
The sheer accumulation of mass by replication and adhesion may not be the most efficient way for a multicellular organism to gain an adaptive edge. There may be instead advantages from qualitative changes in shape. We ask how replication and adhesion, founded in molecular and cell biology, can enable the topological transformation of organism shape.
A topological transformation of bioontology interest is the adoption of a tubular form. It enables a primitive organism to maintain a regulated internal environment with openings for exchanging nutrients and waste products with the external environment [
49]. A hollow morphology that hosts a separated internal environment arises in some of the earliest multicellular organisms in the Precambrian period 580 million years ago, such as those from the phylum Cnidaria [
62,
63]. Fibonacci forms are found among these tubular plants and animals [
6,
7].
Another topological form of biological interest is the torus. It is a surface of revolution of a circle about a coplanar axis, reminiscent of a “donut" [
64]. An internal torus form may offer a circulatory system. It may require a mechanical pump to impel cyclical fluid motion corresponding to a heart. There exist algorithms to construct three-dimensional objects such as a torus from 2-simplex meshes [
65]. A biological topic for investigation is whether there exist cellular and molecular mechanisms that might act like such algorithms to produce a torus. A torus provides greater control of the internal environment than a tube since the torus shape does not include a mouth-like opening. Instead, it relies on diffusion for chemical exchange with the environment. Among vertebrates, the circulatory system exchanges oxygen and carbon dioxide with the external environment by diffusion in the lungs.
We focus on the adoption by organisms of simple topological forms, the simplices. A topological k-simplex is a generalization of a directed multigraph having k vertices and a directed edge connecting every vertex pair, a face connecting every vertex triplet, and so on. A simplex reflects the most complex topological form representable by the fewest points. For a biological interpretation, we view cells as vertices and adhesions as edges, and we ask whether these roles might enable simplicial topological transformation. The study of simplices might offer insight into the genesis of complex biological forms from simpler ones by cellular replication and adhesion programs. A 0-simplex is a point, a 1-simplex is two points with an edge between them, a 2-simplex is a triangle having an outer boundary and a flat inner region, and so on.
In topology, there is an equivalence termed a homotopy between shapes that can be transformed into each other. If they can be continuously and invertibly transformed into each other, they are termed a homeomorphism. For example, the boundary of a 2-simplex is the maximally continuously contracted form of a circle, and a circle is a continuously contracted form of a tube. Transitively, the boundary of a 2-simplex is the maximally contracted form of a tube. We ask if there is an algebraic way to replicate the vertices of a n-simplex to produce an -simplex. If so, we ask if there exists a cellular and molecular apparatus for executing this algebraic procedure to produce a biological 2-simplex. If so, we further ask if there is a cellular and molecular apparatus for extending a biological 2-simplex toward a circle of enlarging cellular diameter. A desired property of the cellular and molecular apparatus is that it enables a self-similar automaton program where individual cells act similarly in response to similar stimuli.
We define a simplex in terms of the barycentric coordinate system.
Definition 1 (
n-Simplex)
. Given points in that are affinely independent, the n-simplex is the set of all points
The affine dependence criterion assures that the simplex properties are preserved by an affine transform, which can include a linear scaling, a rotation, or a translation.
Definition 2 (Vertex Replication). In an n-simplex with points , the replication of the point is given by a process:
Increase the space dimensions from to .
This increases the vector length of every point. Let be the first points. Let each be a copy of with a 0 placed in the last position.
Let be a copy of the point to be replicated, , with a 1 placed in the last position.
Here is an example progression from a 0-simplex to a 1-simplex to a 2-simplex by two vertex replications.
We find that a vertex replication in an
n-simplex yields an
-simplex.
Theorem 1 (Simplex Growth by Vertex Replication)
. A vertex replication in an n-simplex produces an -simplex. The new simplex is given by the set of points
Proof. By assumption, prior to replication, the points in are affinely independent. It suffices to show that the produced points in are affinely independent. This holds because the point produced by replication, , differs from all other points by having a 1 in its last vector position, whereas the other points have a 0 in that position. Therefore is affinely independent of all other points, which by assumption are affinely independent. □
In the contact inhibition model of Fibonacci cell count growth, an adhered cell pair assigns one to be a partner cell and the other a mitosis cell. Upon telophase, the mitosis cell forms a contractile ring. At cytokinesis, the contractile ring becomes the plane of abscission (
Figure 1). We hypothesize that if the pair polarity is perpendicular to the mitosis polarity, then upon abscission the polarity of the progeny cells will be perpendicular to the preceding pair polarity. This perpendicularity is equivalent to affine independence of the progeny adhesion direction to the preceding pair polarity. In addition, if the mitosis polarity of the mitosing cell relative to the partner cell is such that their adhesion spans the contractile ring, then both progeny cells will inherit adhesion to the partner cell, in addition to retaining adhesion to each other, as shown in the first three steps here
According to the simplex growth by vertex replication theorem, if the mitosis polarity is perpendicular to the pair polarity and if parental adhesions are inherited, then the original single cell produces a 2-simplex at the third generation.
The capacity for a mitosing cell to split and share adhesions among progeny cells to a partner cell as in Equation
6 implies particular properties of cell membranes and their embedded adhesion molecules. Rather than be an infinitesimally thin vertex between nodes as in graph theory, this implies an adhesion to be a bundle with a certain thickness that can be split. This appears to be consistent with the fluid mosaic biochemical theory where "cell membranes are viewed as two-dimensional solutions of oriented globular proteins and lipids" [
66]. Under this interpretation, the adhesion on the left hand side of Equation
6 is formed by an arrangement of oriented proteins shared by segments of the liquid membranes of the two partner cells on the left. If the liquid zone of the dividing cell spans the plane of abscission, then both progeny cells inherit an adhered liquid zone on their cell membrane. If the two progeny cells retain adhesion by oriented proteins between them as their liquid membrane zones split across the plan of abscission, the resulting three cells will have the triangular shape of a 2-simplex. These phenomena are consistent with the view of a cell membrane as a dynamic structure [
67].
For mitoses starting with the fourth generation, the partner cell’s adhesion does not span the contractile ring or the plane of abscission. The restriction of the adhesion to one side of the plane of abscission allows the 2-simplex to grow as a circle shape with every cell having a similar automaton program. The ensuing steps are not simplex progression because the orientations of the replication progeny are no longer affinely independent. After the 2-simplex the structure takes the form of a graph. An automaton program of self-similar replication and adhesion may be postulated for growth from a circle to a tube.
If the polarity of mitosis in the 1-simplex relative to the partner cell is such that adhesion does not span the contractile ring, then a 2-simplex will not form, as illustrated here.
A tubular morphology opens an opportunity for cellular specialization [
49]. A two-layered cellular boundary that segregates an internal environment appears in some of the earliest multicellular organisms, such as those from the phylum Cnidaria [
49,
62,
63,
68]. The outer layer is termed the epidermis, and the inner layer the gastrodermis [
63].
In vertebrates, tubular structures appear as the aero-digestive tract, the urinary tract, or endocrine/glandular secretion tracts—the hole in the 2-simplex functions as a mouth for nutrient entry into a regulated internal environment. There are internal structures that do not have a “mouth,” such as the circulatory system, the lymphatic system, and the cerebrospinal fluid system. These do not contract to a 2-simplex.
With classical power of two replication, two cells mitose simultaneously in the transition from the second to the third generation. The topological shape of the four-cell adhered product depends on each cell’s relation to the other’s contractile ring and abscission plane. If the mitotic polarity of each dividing cell is perpendicular to that of the other cell, and if each progeny retains adhesion to its sibling and inherits the adhesion of its parent, then the four-cell unit should take the shape of a tetrahedron. This shape is akin to a 3-simplex, which is the simplest contraction of a sphere. It has an internal environment without an opening to exchange nutrients and waste products with the external environment. If neither cell’s adhesions span the other’s contractile ring, then the four-cell result should take the shape of a row, as in Equation
8.
5. A Visual Representation for Algebraic Topology of Replication and Adhesion
While Fibonacci replication with adhesion across three generations to produce a 2-simplex offers the most direct path for an organism to form an internal environment, we wish to identify other paths toward hollow organs. For example, we may explore the more general ontogenic circumstance where asynchronously stem cells migrate into a vicinity, individually replicate, and then aggregate into cellular assemblies possibly having topological holes.
We offer a visual representation for symbolizing algebraically the replication and adhesion status of cellular aggregations, while depicting visually the opportunities for topological transformations. Plural adhesion slots per cell may be depicted. An arrow symbol gives the adhesion direction. The arrowhead is the = symbol, and adhesions on the left-hand side match to the right-hand side. Branches on either side of the arrowhead are interpreted as addition, +.
We select the arrowhead as the = symbol because the molecular apparatus of cellular adhesion and the extracellular matrix both appear to be anisotropic [
27,
69,
70]. While this implies that cell adhesion slots are anisotropic, this assumption is unnecessary. The visual representation could be adjusted accordingly, such as by using double-headed arrows.
In the following examples, we treat cells as replicating in Fibonacci patterns, but this restriction is unnecessary with this visual representation. Let us assume that a progenitor cell replicates in generations with cell counts given by the Fibonacci numbers . We assume that the two progeny cells of a replication adhere to each other. We refer to these as the horizontal cell adhesion slot for the convenience of representation on the printed page. We further assume that the individual cells in an aggregation may adhere to cells in another aggregation. We refer to these as the vertical cell adhesion slots for presentation on the printed page. We provide an over and under arrow notation to track adhesions between cells of different aggregations.
An example is
Figure 4 having cell aggregations of sizes
,
, and
. In
Figure 4, the cells have only one vertical receptor. Therefore the assembly of
Figure 4 cannot anneal with other Fibonacci-sized cell aggregations into larger structures.
The function of the horizontal and vertical cell adhesion slots may be served by the same molecular pattern, but we segregate their visual representation above and below the cell count symbol because they serve different roles. The precise positions on a cell surface are not the focus of analysis. Instead, we are interested in the qualitative topological properties and are less interested in this analysis of morphological forms that can be deformed continuously into each other as topological homotopies.
In this notation, for example, we might represent
. These three assemblies may adhere across the = sign. We introduce an over/under arrow notation to indicate the top or bottom inter-aggregation adhesion sites. The arrowhead may be read as the equal = sign. For the top vertical slots we may have for example
This may represent the biological process whereby the top inter-aggregation cell adhesion slots of cell aggregations of sizes
and
adhere to and occupy the top adhesion slots of a cell aggregation of size
. The vertical slots have sizes given by the Fibonacci number. This notation is specific to this application and is unrelated to the Fibonacci number of a graph [
71]. This triplet is the simplest Fibonacci adhesion event involving cell aggregations of different sizes.
We exclude from consideration reflexive arrangements such as
and
that interrupt the chaining.
With this restriction, the presence of two vertical slots implies a potential for infinite chaining. For example, a set of Fibonacci-sized cell aggregations
produced by kinetically tuned mitosis might assemble by overlapping triplet adhesion into a scaffold as
Such a structure may display a self-similar scaffold across the scales of size. Such a bio-fractal form may offer the biological efficiency of reuse of the same cell-to-cell adhesion molecules to maintain structural integrity across spatial scales. The same molecular apparatus that adheres one cell to another may adhere a lobule to a lobule, up to a limb to a trunk.
The Fibonacci numbers bring a rich set of combinatorial identities [
72]. These combinatorial identities may be visually encoded to depict adhesion patterns between cells and cellular aggregations. For example, some Fibonacci identities produce a larger Fibonacci number from a collection of smaller ones. An example is
If this represented a collection of Fibonacci-sized aggregations that adhered as
then as per the identity the assembly would have
cells. Accordingly, it may be contracted to the equivalent of a 0-simplex
When Fibonacci-sized cellular aggregations combine, the resulting cell count is given by the theory of Fibonacci compositions [
73].
7. Combinatorial Properties of Circular Forms
If several ring shapes of the same diameter were to stack up above and below the plane of the page and each cell were to have a further pair of unoccupied inter-aggregation adhesion slots oriented in the third dimension perpendicular to the plane of the page, then the rings could anneal into a tubular structure. We use the Fibonacci numbers to count the number of cells in a cellular aggregation, but the Fibonacci numbers offer a different combinatorial interpretation. For an open form of length
n cell diameters, the Fibonacci number
carries the combinatorial interpretation of counting the number of ways to fill those
n cell diameters with single cells and cell pairs. For a ring with circumference
n cell diameters, that count is given by the Lucas number
[
72,
74, p. 17]. For a given number of cell diameters, the Lucas numbers are larger than the Fibonacci numbers, as per the identity
recalling that in our convention
. This combinatorial relation is illustrated for a cellular 2-simplex in
Figure 5. As with the Fibonacci numbers, the ratio between adjacent Lucas numbers tends toward the Golden Ratio. Since the Lucas number is larger than the Fibonacci number for the same number of cells, a ring shape as a template offers a larger number of ways for single cells and pair cells to anneal to it than an open shape.
Figure 5.
Cellular adhesion interpretation of the tilings heuristic for Fibonacci and Lucas numbers. There is a tiling of 1 cell (gray) and one of two cells (yellow). The number of tilings for a linear cell arrangement of length n is given by the Fibonacci number . The number of tilings for a closed arrangement is given by the Lucas number where .
Figure 5.
Cellular adhesion interpretation of the tilings heuristic for Fibonacci and Lucas numbers. There is a tiling of 1 cell (gray) and one of two cells (yellow). The number of tilings for a linear cell arrangement of length n is given by the Fibonacci number . The number of tilings for a closed arrangement is given by the Lucas number where .
Figure 6.
A Scaffold Ring. All inter and intra-aggregation adhesion slots are occupied. The arrows depict the inter-aggregation adhesion slot bindings.
Figure 6.
A Scaffold Ring. All inter and intra-aggregation adhesion slots are occupied. The arrows depict the inter-aggregation adhesion slot bindings.
A circular assembly of cells appears to have greater entropy under a combinatorial Lucas tiling argument than an open assembly. This is further illustrated in Appendix A. If the same combinatorial forces act perpendicularly on the same cells, then a torus might have greater entropy than a tube. Indeed, the topological square of a circle
is a torus
. The shape resembles a circulatory system where pumped blood flows in circles. An idealized circulatory system with a two-chambered heart such as in gnathostome fish (Infraphylum Gnathostomata [
75]) is illustrated in
Figure 7.
Author Contributions
Conceptualization, W.E.B.H., H.A.L., K.K., Z.W., E.L., J.A., T.B.K. and L.T.; methodology, W.E.B.H., H.L.A. and L.T.; formal analysis, W.E.B.H., H.A.L. and L.T.; resources, W.E.B.H.; writing—original draft preparation, W.E.B.H.; writing—review and editing, W.E.B.H., H.A.L., E.L. and L.T.; visualization, W.E.B.H.; supervision, W.E.B.H. and L.T.; project administration, W.E.B.H. and L.T. All authors have read and agreed to the published version of the manuscript.