Multicellular Daisy Chains

• Multicellular daisy chains are sequential control architectures where gene-network outputs generate interaction networks that actuate the next gene expression stage.
• The approach employs dynamic graph formulations and experimental validations—using tools like fluorescent reporters and mutant analyses—to demonstrate stage-specific, robust self-organization.
• This framework not only explains E. coli self-organization but also paves the way for applications in synthetic biology, developmental biology, and bioengineering.

Multicellular daisy chains are a proposed control architecture for self-organizing cellular communities in which gene-network outputs at one multicellular stage generate a characteristic interaction network, and that interaction network then serves as the inducible input for the next stage of gene expression. In this view, a community is treated not simply as a set of cells responding to external nutrients or stresses, but as a dynamic graph whose nodes are cells and whose edges are changing physical or chemical interactions. The concept was advanced to explain how multicellular self-organization can produce robust expression of otherwise noisy genes, with Escherichia coli self-organization providing the motivating case: multicellular interactions act as inputs for key gene networks, while gene networks in turn generate the next multicellular interaction network (Allison, 15 Aug 2025).

1. Formal definition

A multicellular daisy chain is defined as a sequence of coupled mappings of the form gene-network output $\rightarrow$ multicellular interaction network and multicellular interaction network $\rightarrow$ next gene-network input. At stage $i$, the expression profile $g^{(i)}$ of key multicellular genes produces an interaction network $M^{(i)}$ among cells, and $M^{(i)}$ is then sensed by cells to induce the next expression profile $g^{(i+1)}$ (Allison, 15 Aug 2025).

With $g^{(i)} \in \mathbb{R}^k$ denoting the vector of expression levels of $k$ stage-$i$ genes, and $\rightarrow$0 the graph of multicellular interactions at stage $\rightarrow$1, the framework introduces two maps,

$\rightarrow$2

for gene outputs $\rightarrow$3 interaction network, and

$\rightarrow$4

for interaction network $\rightarrow$5 next gene inputs, such that

$\rightarrow$6

Chaining stages gives the composite

$\rightarrow$7

The central claim of the model is that each stage is activated by the internal, self-generated $\rightarrow$8 rather than by noisy external cues. On that basis, the succession of stages becomes inherently robust and self-referential. This suggests a form of developmental control in which a community’s later states are determined by interaction structures that it has already constructed.

2. Dynamical-graph formulation

The same idea can be written as a computational or dynamical-graph process. Letting $\rightarrow$9 denote the multicellular graph at generation $i$0, cells express gene set $i$1 in response to inputs from $i$2, and then create a new graph $i$3. Schematically,

$i$4

so that

$i$5

When $i$6 and $i$7 alternate in a tightly coupled sequence, graph evolution follows a daisy-chained dynamic in which each generation’s propagation rules derive from prior gene-network outputs (Allison, 15 Aug 2025).

This formulation places multicellular interaction networks at the center of control. The community is not merely a passive substrate on which intracellular regulation unfolds; rather, adjacency, confinement, tensile linkage, diffusional coupling, and related intercellular relations become stage-specific control variables. In this framework, the graph itself is both a product of gene expression and a determinant of subsequent gene expression. A plausible implication is that multicellular organization can function as a hierarchical signal-processing layer superimposed on intracellular regulatory circuitry.

3. Stagewise mechanism in E. coli

In Allison’s account of E. coli self-organization, multicellular daisy chains are instantiated through a succession of mechanical and chemical edges, including Ag43-mediated adhesion, fimbriation, flagellar tension, inner-cavity confinement, curli scaffolds, and diffusible signals. These edges function as either the output of one gene network or the input to the next (Allison, 15 Aug 2025).

Stage	Gene-network output	Interaction-network consequence / next input
1	flu-operon expression $i$8 produces Ag43 proteins	Adhesion edges between sister cells generate $i$9
2	$g^{(i)}$0 activates fim operon expression $g^{(i)}$1 via Rcs and QseBC	Type 1 fimbriae restructure cell-cell edges to $g^{(i)}$2
3	$g^{(i)}$3 induces the flagellar cascade $g^{(i)}$4, governed by FlhDC and YdiV	Asymmetric motility folds sister-cell pairs into rosettes and creates $g^{(i)}$5
4	$g^{(i)}$6 triggers curli and further fim expression $g^{(i)}$7	Polymers fill the cavity and adjust diffusion properties to yield $g^{(i)}$8
5	$g^{(i)}$9 activates rpoS transcription $M^{(i)}$0 and prepares the pga operon for PGA production $M^{(i)}$1	A terminal graph $M^{(i)}$2 is cemented

The stage sequence is described concretely. In stage 1, two-cell adhesion arises because flu-operon expression produces Ag43 proteins that create adhesion edges between sister cells. In stage 2, the tensile or redox shifts generated when sisters pull on each other activate fim operon expression through two-component sensors such as Rcs and QseBC. In stage 3, rigid fimbrial links and mechanical constraints induce the flagellar cascade, and asymmetric motility folds sister-cell pairs into rosettes with an internal tube-like cavity. In stage 4, confinement in that cavity concentrates chemical signals such as AI-3, while tension and osmolarity cues trigger curli and further fim expression. In stage 5, low nutrient flux and a high polymer barrier activate $M^{(i)}$3 transcription, insulating cells from future external variation and preparing the pga operon for PGA production (Allison, 15 Aug 2025).

Each arrow $M^{(i)}$4 and $M^{(i)}$5 is described as having been experimentally documented with fluorescent reporters, mechanical-tension measurements, sensor-mutant analyses, and microfluidic perturbations. Within the model, the developmental sequence is therefore not simply chronological; it is causally chained through interaction-network intermediates.

4. Robustness, canalization, and empirical support

The proposed robustness of multicellular daisy chains is attributed to hierarchical signal isolation. Because every stage’s inputs $M^{(i)}$6 are an internal function of the previous stage’s gene outputs, the community effectively re-encodes its environment in a multi-level filter. External fluctuations that do not disrupt the controlled edges of stage $M^{(i)}$7 are buffered from stage $M^{(i)}$8. The framework further emphasizes that many gene networks, including the curli master regulator CsgD and the Cpx/Rcs two-component systems, require multiple orthogonal signals such as mechanical plus chemical cues, thereby creating Boolean-like AND gates that distinguish the correct $M^{(i)}$9 from background noise (Allison, 15 Aug 2025).

This architecture is presented as enforcing stepwise canalization: only communities that correctly complete stage $M^{(i)}$0 can advance to stage $M^{(i)}$1. In that sense, developmental progression is linked to successful construction of the appropriate multicellular interaction network rather than solely to intracellular timing.

Several lines of evidence are cited.

• Synchronous activation during free self-organization: in mother-machine microchannels, many multicellular genes, including flagella, $M^{(i)}$2, and $M^{(i)}$3, exhibit stochastic or noisy transcription, but during free clonal self-organization these same genes switch on synchronously and nearly uniformly at their designated stages (Allison, 15 Aug 2025).
• Stage-specific developmental arrest in mutants: mutations that break individual links in the chain, including $M^{(i)}$4, $M^{(i)}$5, $M^{(i)}$6, $M^{(i)}$7, $M^{(i)}$8, and $M^{(i)}$9, arrest development precisely at the predicted stage, indicating the causal necessity of each input-output pair (Allison, 15 Aug 2025).
• Dynamic-graph simulations: mathematical simulations of dynamic graphs augmented by stage-specific logic functions $g^{(i+1)}$0 and $g^{(i+1)}$1 confirm that only a daisy-chained sequence can propagate large communities to a stable terminal graph without either collapsing through excess adhesion or dispersing through excess motility (Puri et al., 4 Mar 2025).

Taken together, these observations are used to argue that multicellular daisy chains provide predictability not by removing stochasticity at the level of individual genes, but by embedding gene-network transitions in a constrained sequence of self-generated multicellular states.

5. Terminological scope and relation to chain morphologies

The term “daisy chain” has a distinct usage in the statistical hydrodynamics of multicellular colonies. In work on flagellated swimming colonies, planar “daisy-chain” colonies refer to one-dimensional colonial morphologies rather than to stagewise control architectures. That literature studies transport and fluctuation properties of chainlike colonies using Ornstein–Uhlenbeck, or continuous-time AR(1), models for cycle-averaged flagellar force magnitude $g^{(i+1)}$2 and orientation $g^{(i+1)}$3, together with demographic variability in relaxed force, relaxed angle, and basal-attachment offset (Ashenafi et al., 2023).

For an $g^{(i+1)}$4-cell planar chain, the analysis provides closed-form expressions for mean swimming speed $g^{(i+1)}$5, translational diffusivity $g^{(i+1)}$6, orientational drift $g^{(i+1)}$7, rotational diffusivity $g^{(i+1)}$8, and orientation correlation time $g^{(i+1)}$9. It also derives large-$g^{(i)} \in \mathbb{R}^k$0 scalings and emphasizes the role of geometry, including the dependence of disorder effects on the spacing parameter through factors such as $g^{(i)} \in \mathbb{R}^k$1 (Ashenafi et al., 2023).

The distinction is substantive. In the hydrodynamic setting, a daisy chain is a geometric colony type whose dynamics follow from rigid-body low-Reynolds-number hydrodynamics, linear Stokes drag, and small-noise separation of time scales. In the multicellular-control setting, a daisy chain is a sequence of gene-network and interaction-network mappings that propagates development. A common misconception is to identify the latter with any chain-shaped multicellular arrangement. The available descriptions do not support that equivalence: multicellular daisy chains are defined by input-output coupling across stages, not by one particular morphology.

6. Broader implications and prospective extensions

Several implications are proposed for synthetic biology, developmental biology, and bioengineering. In synthetic biology, artificial daisy chains could be architected by linking synthetic adhesins, quorum-sensing circuits, two-component sensors, and polymer-export modules so that engineered cells build multilayered structures with predictable stage-by-stage gene induction. Recent work on synthetic cell-cell signaling by Toda et al. and synthetic adhesion molecules by Stevens et al. is described as potentially constituting the first modules of designer multicellular daisy chains (Allison, 15 Aug 2025).

In developmental biology, the framework recasts classical positional information into a graph-based mechanism in which cells read and rewrite their adjacency matrix as they differentiate. This perspective may help reverse-engineer complex metazoan lineages by identifying which cell-cell interactions serve as the critical inputs for each differentiation switch. In bioengineering and organoids, the proposal is that mimicking only the minimal subset of intermediate interactions $g^{(i)} \in \mathbb{R}^k$2, rather than an entire tissue environment, may be sufficient to coax stem cells down precise developmental trajectories with greater efficiency; early successes in spatially patterned organoid assembly by Hofer and Lutolf are described as hinting at this possibility (Allison, 15 Aug 2025).

The broader claim is that multicellular daisy chains provide a unifying model in which self-generated community structures function simultaneously as command signals and execution hardware for gene networks. The model therefore attributes developmental simplicity and evolutionary flexibility to the same architecture: a community’s future is written by its own past, yet individual links in the chain remain, in principle, rewritable. This suggests that multicellular control may, in some biological contexts, be simpler than the sum of the constituent intracellular parts.