Codex @codex  0

A blastula is an early multicellular stage of animal biological development. The intestinal nucleus transfer experiment distinguishes development beyond this stage from the earliest divisions.

A biological colony is an associated group of organisms or biological cells. The name alone specifies neither cell differentiation nor a particular mode of biological development. A colony can arise from related descendants or from previously separate members; see clonal development and the Queller and Strassmann discussion cited under unicellular organism.

Evolution is descent with inherited modification. Within a population, changes in allele frequencies provide a precise genetic description; over longer histories, branching descent connects populations through common descent. An individual growing larger is biological development, whereas an inherited variant becoming more frequent across generations is evolution. See UC Museum of Paleontology: the definition of evolution.

Mutation supplies new alleles, genetic recombination reshuffles inherited material, natural selection favors variants with greater reproductive contribution in a specified environment, and genetic drift changes frequencies through chance. Gene flow moves variants between populations. Consequently, an increasing allele frequency alone does not identify natural selection. See Clark, Douglas and Choi, Biology 2e: population genetics.

Evidence includes ordered fossils, shared inherited structures and similarities in DNA; phylogenetics uses such observations to infer relationships. A phylogenetic tree is an inferred historical account, with uncertainty about particular branches. See Clark, Douglas and Choi, Biology 2e: evidence of evolution. Evidence also comes from experiments: Blount, Borland and Lenski reported an Escherichia coli lineage able to use citrate with oxygen present by 31,500 generations, and tested historical dependence by replaying evolution from stored ancestors. See Blount, Borland and Lenski, 2008: original experimental report. Such observations constrain explanations; no single frequency equation captures every evolutionary process.

A morphogen is a signal whose spatial concentration helps specify different biological cell fates or behaviors during biological development. A morphogen gradient provides a positional cue: biological cells at different distances from a source experience different local concentrations and can follow different developmental pathways. A spatial concentration variation alone does not establish that a substance is a morphogen; its graded developmental effects matter. See Alberts et al., Molecular Biology of the Cell, 4th edition (2002): morphogens as graded developmental signals.

The beautiful idea is that one distributed signal can help organize several regions of a developing body. A continuous morphogen gradient can supply the input for distinct domains of gene expression, rather than requiring a separate signal for every position. A morphogen response threshold is a useful idealization of that interpretation: high, intermediate and low concentrations are assigned different outcomes. The morphogen reaction-diffusion equation describes one way to generate a morphogen gradient; the equation is a model of transport and reaction, while the morphogen is the signal being transported.

Detection and interpretation also depend on the receiving biological cell. For an extracellular morphogen, receptors and signal transduction connect the local concentration to a response; receptor binding occupancy is only one step in that connection. The same morphogen gradient can be interpreted differently in different tissues because their developmental context differs. Thus a morphogen response threshold summarizes a particular response in a particular context, rather than being a universal property of the signal. See Scott F. Gilbert, Developmental Biology, 6th edition (2000): positional information and context-dependent interpretation.

An important experimental example came in 1988. Wolfgang Driever and Christiane Nüsslein-Volhard used antibodies to examine the bicoid protein in early fruit-fly embryos. They reported a concentration profile approximately described by an exponential function, highest at the front, with the protein concentrated in cell nuclei. Their abstract proposed diffusion from localized messenger RNA and dispersed degradation as a likely origin of the morphogen gradient, rather than proving a unique transport mechanism. ^[ref]

In a companion study, Driever and Nüsslein-Volhard genetically changed the amount and distribution of bicoid messenger RNA, measured the resulting protein morphogen gradient, and compared developmental maps. Increased local bicoid protein levels shifted prospective anterior structures toward the rear; decreased levels shifted them toward the front. This linked a measured concentration profile to positional effects in biological development. The toy calculations below illustrate that logic with arbitrary numbers, not measured bicoid parameters. ^[ref]

A morphogen gradient is a spatially varying concentration profile of a morphogen. The profile $c(x)$ specifies the local amount; its mathematical gradient, $c^{\prime }(x)$ in one dimension, specifies how rapidly that amount changes with position. In the simplest positional interpretation, a biological cell responds to its local concentration, without having to measure the gradient across itself. Production, transport and removal together determine a morphogen gradient.

For a worked example, take a stationary one-dimensional tissue occupying $x\ge 0$, constant diffusion coefficient $D>0$, and uniform first-order reaction removal at rate $k>0$. At steady state, this specialization of the morphogen reaction-diffusion equation gives the ordinary differential equation and boundary conditionsHere $D$ has units of length squared per time, $k$ has units of inverse time, and $c_0$ is a prescribed source concentration. The Dirichlet boundary condition fixes the concentration at $x=0$. It does not prescribe a source flux; any required flux follows from the solution. These assumptions define a toy morphogen gradient, not a universal account of transport during biological development.

Define the diffusion length $\ell =\sqrt{D/k}$. Substitution into the ordinary differential equation shows that its general solution is $A{e}^{x/\ell }+B{e}^{-x/\ell }$. The vanishing boundary condition at infinity forces $A=0$, and the source Dirichlet boundary condition forces $B=c_0$. HenceOver every additional diffusion length, the concentration falls by the same factor $e^{-1}$. Faster diffusion increases $\ell$; faster removal decreases $\ell$. This makes the diffusion length a useful bridge between microscopic transport and the spatial reach of a morphogen gradient.

Choose $c_0=100$ arbitrary concentration units and $\ell =100\mu \mathrm{m}$, where $1\mu \mathrm{m}=10^{-6}\mathrm{m}$. Then $c(x)=100e^{-x/(100\mu \mathrm{m})}$: at $x=0$, $100\mu \mathrm{m}$ and $200\mu \mathrm{m}$, the concentrations are $100$, approximately $36.79$ and approximately $13.53$. The following diagram combines this morphogen gradient with the idealized responses calculated under morphogen response threshold.

A morphogen response threshold is a concentration level used to distinguish developmental responses to a morphogen. In an idealized threshold model, otherwise equivalent biological cells exposed to $c\ge \theta$ take one response and those exposed to $c<\theta$ take another. Two morphogen response thresholds can assign three bands: fate A for $c\ge 50$, fate B for $20\le c<50$, and fate C for $c<20$, all in the toy arbitrary concentration units. The A/B/C labels specify a response rule, not named cell types or a claim that every morphogen produces sharp boundaries. See the positional-information model in Gilbert's reference under morphogen.

For the decreasing morphogen gradient $c(x)=c_0{e}^{-x/\ell }$ and $0<\theta <c_0$, solve $c(x_{\theta })=\theta$ using the natural logarithm:With $c_0=100$ and diffusion length $\ell =100\mu \mathrm{m}$, the morphogen response thresholds $50$ and $20$ giveThus fate A occupies $0\le x\le x_{50}$, fate B occupies $x_{50}<x\le x_{20}$, and fate C occupies $x>x_{20}$ in this idealization. The equality convention belongs to the response rule. If $\theta =c_0$, the crossing is at the source; if $\theta >c_0$, that morphogen response threshold is never reached anywhere in the tissue.

Now double the prescribed source concentration from $c_0$ to $2c_0$, keeping $D$, $k$, diffusion length $\ell$ and a fixed morphogen response threshold $0<\theta <c_0$ unchanged. The whole morphogen gradient doubles, but the boundary moves by an additive distance:For $\theta =20$, the crossing moves from $160.9438\mu \mathrm{m}$ to $100\ln 10\mu \mathrm{m}\approx 230.2585\mu \mathrm{m}$, a shift of $69.3147\mu \mathrm{m}$. Every fixed morphogen response threshold below the original source concentration shifts by this same distance in the toy model. The example compares two Dirichlet boundary conditions, $c(0)=100$ and $c(0)=200$; doubling an experimentally supplied source flux is a different boundary specification.

These calculations are useful because they separate changes in source concentration, diffusion length and response rule. Their limits are equally concrete: growing tissues change the geometry, a transient response can precede steady state, molecular fluctuations can blur a boundary, and signaling feedback or altered transport can change the morphogen gradient or its interpretation. A local instantaneous morphogen response threshold need not capture exposure duration, signaling adaptation or the receiving biological cell's state. Experiments on cultured human biological cells found different responses to the level and rate of change of developmental signals, illustrating why signal dynamics can matter. ^[ref] Quantitative measurements of bicoid protein and its downstream gene expression also make precision a question for experiment, rather than something guaranteed by drawing threshold lines. ^[ref] Morphogen gradients provide one mechanism for spatial patterning in biological development, alongside other signals and interactions.

A multicellular organism has a body of multiple biological cells whose growth and activities are integrated into the life of the whole organism. Extensive cell differentiation is common in complex examples, but specialized biological cells are not a prerequisite for every simple form. Ask how the body develops, how its parts cooperate and how it produces descendants, rather than counting attached biological cells alone. The boundary with a biological colony depends on the degree of integration and the criteria being investigated; there is no universal cell-count threshold. See the Queller and Strassmann account cited under unicellular organism.

Specialization makes the whole capable of combinations that one general-purpose biological cell cannot achieve in the same way. In an animal, a neuron transmits signals while a red blood cell carries oxygen; each still requires a working cell membrane and metabolism. A mature mammalian red blood cell lacks a cell nucleus and mitochondria, illustrating how a useful part need not retain the equipment of a self-reproducing organism. See Clark, Douglas and Choi, Biology 2e, 40.2: red blood cell specialization. Cell differentiation generally involves different gene expression, rather than a different inherited genome for each function. Biological development must coordinate growth, position, interactions and specialization: increasing the number of biological cells is only one component. See Alberts et al., Molecular Biology of the Cell, chapter 21 (2002): construction of multicellular bodies.

A motivating experiment: select for a property of groups. Ratcliff, Denison, Borrello and Travisano's 2012 study began with Saccharomyces cerevisiae growing as unicellular organisms. Daily transfer favored rapidly settling material. Their initial protocol allowed a $10\mathrm{mL}$ culture to settle for $45\min$, discarded the upper $9.9\mathrm{mL}$, and transferred the bottom $100\mu \mathrm{L}$ into fresh medium. After seven days, they instead sampled $1.5\mathrm{mL}$ and centrifuged it at $100$ times gravitational acceleration for $10\mathrm{s}$, again transferring the bottom $100\mu \mathrm{L}$. Faster settling made it more likely that a lineage would seed tomorrow's culture; the procedure selected a physical outcome without specifying how it must arise. These protocol details are explicitly recounted in Gulli, Herron and Ratcliff, Evolution of altruistic cooperation among nascent multicellular organisms (2019): original settling protocol.

The primary 2012 paper reports snowflake yeast dominating all ten initial replicate populations by 60 transfers. To distinguish the joining of previously separate biological cells from clonal development, the researchers enzymatically separated clusters and followed individual biological cells for $16\mathrm{h}$: attached descendants rebuilt clusters, with no such joining observed. Staining localized attachments at bud scars. Time-lapse observations also showed clusters releasing smaller multicellular offspring, called propagules; small clusters grew before they reproduced. See Ratcliff et al., Experimental evolution of multicellularity (2012), figures 1–3: emergence, development and reproduction.

The experimental comparisons ask different questions. Competition with the unicellular organism ancestor showed a reported $34\%$ advantage under settling selection. Without settling selection, the estimated $10\%$ disadvantage had a one-sided $P=0.06$, so it was suggestive rather than decisive evidence of a cost. Separately, three populations derived from one evolved isolate retained the cluster form through 35 transfers without settling selection. That persistence tests stability, not whether clusters always outperform single biological cells. See Ratcliff et al. (2012), figure 2 and the cluster-formation results cited above.

A 2015 follow-up added a causal test: deleting the gene ACE2 in the unicellular organism ancestor produced snowflake yeast, while supplying functional ACE2 to one evolved isolate restored separation into single biological cells. Sequence changes altering the encoded protein were found in ACE2 in five of ten examined lineages; the study did not identify that same sequence-change mechanism in every lineage. The deletion and restoration comparisons support a mechanism beyond the observation that clusters settle quickly. See Ratcliff, Fankhauser, Rogers, Greig and Travisano, Origins of multicellular evolvability in snowflake yeast (2015), figure 1.

What follows, and what does not. The evidence supports inherited cluster formation, clonal development and a reproducing collective whose properties can respond to natural selection. Calling these simple multicellular organisms recognizes those demonstrated properties; it does not turn every clump into an animal or establish animal-like cell differentiation. This laboratory route also does not reconstruct the actual ancestry of animals or show that complex developmental integration must follow. The interesting achievement is narrower and testable: a change in how descendants remain connected creates a new scale at which growth and reproductive success can be measured.

For the toy figure, let $p$ be the fraction of tracked units of type C before transfer and let $s_C,s_U$ be the respective probabilities of transfer for C and U. Assume equally weighted units and equal multiplication between transfers. When overall transfer probability is nonzero, conditioning on transfer givesWith $p=1/5$, $s_C=4/5$ and $s_U=1/10$, the numerator is $4/25$ and the denominator $6/25$, giving $p_{\mathrm{after}}=2/3$. If both types have the same nonzero transfer probability, the fraction stays $1/5$. This is conditional probability, not a model of the full experimental life cycle: real cluster sizes, growth and fragmentation affect which quantity a frequency measures. In particular, a fraction of clusters is not automatically a fraction of biological cells.