(A) In all simulated and experimental communities (see ‘Materials and methods’), two populations of cells, marked in red and green, were initially randomly distributed on a surface unless otherwise …
Parameter values used in the fitness model.
(A) To incorporate realistic assumptions about cell rearrangement in three-dimensional communities, we monitored the growth of a single fluorescent yeast cell into a microcolony on top of solid …
(A) The fitness model demonstrated that interactions conferring large benefit to at least one partner could lead to ratio convergence. Simulation parameters are listed in Figure 1—source data 1. (B) …
The fitness model shows that if two partner populations inhibit each other sufficiently strongly compared to their basal fitness, population ratios may diverge rapidly.
Competitive communities of strains with equal fitness (equal-fitness competition, abbreviated as ‘Eq-fitness Comp.’) showed population segregation as suggested by static late-stage top-views (A, …
Definitions and values of parameters used in the diffusion model.
Engineered yeast communities consisted of two non-mating fluorescent strains of yeast, G and R, engaged in a metabolic interaction. (A) In competitive communities, two prototrophic strains engaged …
In the diffusion model, community grew on the top surface of an agarose column (light blue). Cell grids (black borders) represented cells at different states: live R (light red), live G (light …
(A) Starting from densities and ratios spanning orders of magnitude, experimental obligatory cooperative yeast communities developed a consistent characteristic vertical patch size of ∼10 to 20 …
Vertical cross-sections in (A) and (B) were from different growth stages of obligatory cooperative yeast communities starting at total 500 cells/mm2 and 1:1 population ratio. Vertical cross-section …
(A) In engineered yeast communities, even though both obligatory cooperative and non-cooperative communities directly above the high-density inoculation spot showed high population intermixing …
Parameter values used in the fitness model in Figure 3B.
(A) Partner intermixing is insensitive to diffusion constants. In the diffusion model, obligatory cooperative communities with diffusion constants ranging from 20 to 360 µm2/s (corresponding to the …
(A) In the absence of sulfate and hydrogen, the bacterium Desulfovibrio vulgaris (Dv) and the archaeon Methanococcus maripaludis (Mm) cooperate through redox coupling. Dv ferments lactate and …
(A-C) Examples of networks in which cooperative pairs intermixed (A), non-cooperative pairs intermixed (B), or cooperative pairs did not intermix (C) are shown. In the schematic network diagrams, …
Yeast cells bud to the sides when there is available space and bud upward when sufficiently confined (corresponding to Figure 1—figure supplement 1A). To infer the process of cell rearrangement in …
Top views of an equal-fitness competitive community suggest population segregation (corresponding to Figure 2A, left). Competitive communities of strains with equal fitness showed population …
In unequal-fitness competition, the fitter population gradually covers the less fit population (corresponding to Figure 2A, right). Here, G is fitter than R. The community started from a uniform …
Top views of an obligatory commensal community suggest population layering (corresponding to Figure 2B). For a detailed explanation of the growth kinetics of the community (), please refer to the …
Top views of an obligatory cooperative community suggest populations intermixing (corresponding to Figure 2C). For a detailed explanation of the growth kinetics of the community (), please refer …
Summary of steady-state occupancy, conditions to achieve steady-state, and the stability of steady-state for six types of ecological interactions.
Interaction | Steady-state occupancy | Steady-state condition | Stability |
---|---|---|---|
G[∼∼]R | Any | Only when | Unstable |
G[↑ ↑]R | Globally stable | ||
G[∼↓]R | Unstable | ||
G[∼↑]R | Globally stable | ||
G[↓↓]R | Unstable | ||
G[↓↑]R | Locally stable (when ; see Figure 1—figure supplement 2) |
MATLAB code for simulating cooperative communities in the fitness model.
Other ecological interactions can be simulated by changing the values of the interaction parameters r12 and r21.
MATLAB code for simulating cooperative communities in the diffusion model.
In this version, the diffusion constant inside the community is assumed to be similar to the diffusion constant of agarose (Figure 2—figure supplement 2A).
MATLAB code for simulating competitive communities in the diffusion model.
In this version, the diffusion constant inside the community is assumed to be smaller than the diffusion constant of agarose (Figure 2—figure supplement 2B).
MATLAB function for in-plane rearrangement of cells in both the fitness and the diffusion models.