Metabolic constraints drive self-organization of specialized cell groups

How phenotypically distinct states in isogenic cell populations appear and stably co-exist remains unresolved. We find that within a mature, clonal yeast colony developing in low glucose, cells arrange into metabolically disparate cell groups. Using this system, we model and experimentally identify metabolic constraints sufficient to drive such self-assembly. Beginning in a uniformly gluconeogenic state, cells exhibiting a contrary, high pentose phosphate pathway activity state, spontaneously appear and proliferate, in a spatially constrained manner. Gluconeogenic cells in the colony produce and provide a resource, which we identify as trehalose. Above threshold concentrations of external trehalose, cells switch to the new metabolic state and proliferate. A self-organized system establishes, where cells in this new state are sustained by trehalose consumption, which thereby restrains other cells in the trehalose producing, gluconeogenic state. Our work suggests simple physico-chemical principles that determine how isogenic cells spontaneously self-organize into structured assemblies in complimentary, specialized states.


Introduction
During the course of development, groups of isogenic cells often form spatially organized, interdependent communities. The emergence of such phenotypically heterogeneous, spatially constrained sub-populations of cells is considered a requisite first step towards multicellularity. Here, clonal cells proliferate and differentiate into phenotypically distinct cells that stably coexist, and organize spatially with distinct patterns and shapes (Newman, 2016;Niklas, 2014). Through such collective behavior, groups of cells can maintain orientation, stay together, and specialize in different tasks through the division of labor, while remaining organized with intricate spatial arrangements (Ackermann, 2015;Newman, 2016). In both eukaryotic and prokaryotic microbes, such organization into structured, isogenic but phenotypically heterogeneous communities, is widely prevalent, and also reversible (Ackermann, 2015). Such phenotypic heterogeneity within groups of clonal cells enables several microbes to persist in fluctuating environments, thereby providing an adaptive benefit for the cell community (Wolf et al., 2005;Thattai and van Oudenaarden, 2004).
A well studied example of spatially organized, phenotypically heterogeneous groups of cells comes from the Dictyostelid social amoeba, which upon starvation transition from individual protists to collective cellular aggregates that go on to form slime-molds, or fruiting bodies (Bonner, 1949;Du et al., 2015;Kaiser, 1986). Indeed, most microbes show some such complex, heterogeneous cell behavior, for example in the extensive spatial organization within clonal bacterial biofilms and swarms (Kearns et al., 2004;Kolter, 2007), or in the individuality exhibited in Escherichia coli populations (Spudich and Koshland, 1976). Despite its popular perception as a unicellular microbe, natural isolates of the budding yeast, Saccharomyces cerevisiae, also form phenotypically heterogeneous, multicellular communities (Cáp et al., 2012;Koschwanez et al., 2011;Palková and Váchová, 2016;Ratcliff et al., 2012;Váchová and Palková, 2018;Veelders et al., 2010;Wloch-Salamon et al., 2017). However, despite striking descriptions on the nature and development of phenotypically heterogeneous states within groups of cells, the rules governing the emergence and maintenance of new phenotypic states within isogenic cell populations remain unclear.
Current studies emphasize genetic and epigenetic changes that are required to maintain phenotypic heterogeneity within a cell population (Ackermann, 2015;Sneppen et al., 2015). In particular, many studies emphasize stochastic gene expression changes that can drive phenotypic heterogeneity (Süel et al., 2007;Ackermann, 2015;Balázsi et al., 2011;Blake et al., 2003). Further, groups of cells can produce adhesion molecules to bring themselves together (Halfmann et al., 2012;Halme et al., 2004;Octavio et al., 2009;Váchová and Palková, 2018), or support possible codependencies (such as commensal or mutual dependencies on shared resources) within the populations (Ackermann, 2015). Such studies now provide insight into why such heterogeneous cell groups might exist, and what the evolutionary benefits might be. However, an underlying biochemical logic to explain how distinct, specialized cell states can emerge and persist in the first place is largely absent. This is particularly so for isogenic (and therefore putatively identical) groups of cells in seemingly uniform environments. In essence, are there simple chemical or physical constraints, derived from existing biochemical rules and limitations, that explain the emergence and maintenance of heterogeneous phenotypic states of groups of clonal cells in space and over time?
Contrastingly, a common theme occurs in nearly all described examples of phenotypically heterogeneous, isogenic groups of cells. This is a requirement of some 'metabolic stress' or nutrient limitation that is necessary for the emergence of phenotypic heterogeneity and spatial organization, typically in the form of metabolically inter-dependent cells (Ackermann, 2015;Campbell et al., 2016;Cáp et al., 2012;Johnson et al., 2012;Liu et al., 2015). This idea has been explored experimentally, where approaches that systems-engineer metabolic dependencies between non-isogenic cells can result in interdependent populations that constitute mixed communities (Campbell et al., 2016;Campbell et al., 2015;Embree et al., 2015;Wintermute and Silver, 2010). These findings suggest that biochemical constraints derived from metabolism may determine the nature of phenotypic heterogeneity, and the spatial organization of cells in distinct states within the population. Therefore, if we can understand what these biochemical constraints are, and discern how metabolic eLife digest Under certain conditions, single-celled microbes such as yeast and bacteria form communities of many cells. In some cases, the cells in these communities specialize to perform specific roles. By specializing, these cells may help the whole community to survive in difficult environments. These co-dependent communities have some similarities to how cells specialize and work together in larger living things -like animals or plants -that in some cases can contain trillions of cells.
Research has already identified the genes involved in creating communities from a population of identical cells. It is less clear how cells within these communities become specialized to different roles. The budding yeast Saccharomyces cerevisiae can help to reveal how genetic and environmental factors contribute to cell communities.
By growing yeast in conditions with a low level of glucose, Varahan et al. were able to form cell communities. The communities contained some specialized cells with a high level of activity in a biochemical system called the pentose phosphate pathway (PPP). This is unusual in low-glucose conditions. Further examination showed that many cells in the community produce a sugar called trehalose and, in parts of the community where trehalose levels are high, cells switch to the high PPP state and gain energy from processing trehalose.
These findings suggest that the availability of a specific nutrient (in this case, trehalose), which can be made by the cells themselves, is a sufficient signal to trigger specialization of cells. This shows how simple biochemistry can drive specialization and organization of cells. Certain infections are caused by cell communities called biofilms. These findings could also contribute to new approaches to preventing biofilms. This knowledge could in turn reveal how complex multi-cellular organisms evolved, and it may also be relevant to studies looking into the development of cancer.
states can be altered through these constraints, this may address how genetically identical cells can self-organize into distinct states.
In this study, using clonal yeast cells, we experimentally and theoretically show how metabolic constraints imposed on a population of isogenic cells can determine the production, accumulation,  Reconstructed bright-field images of a mature wild-type colony. Within the colony, a network of dark and bright regions is clearly visible, as classified based purely on optical density. We classify the cells in the dark region as dark cells, and in the peripheral light region as light cells. Scale bar: 2 mm. (C) Spatial distribution of mCherry fluorescence across a colony, indicating the activity of (i) a reporter for hexokinase (HXK1) activity, or (ii) a gluconeogenesis dependent reporter (PCK1), in two different colonies. The percentage of fluorescent cells (in isolated light and dark cells from the respective colonies) were also estimated by flow cytometry, and is shown as bar graphs. Scale bar: 2 mm. Also see  and utilization of a specific, shared resource. The selective utilization of this resource enables the spontaneous emergence and persistence of cells exhibiting a counter-intuitive metabolic state, with spatial organization. These metabolic constraints create inherent threshold effects, enabling some cells to switch to new metabolic states, while restraining other cells to the original state which produces the resource. This thereby drives the overall self-organization of cells into specialized, spatially ordered communities. Finally, this group of spatially organized, metabolically distinct cells confer a collective growth advantage to the community of cells, rationalizing why such spatial self-organization of cells into distinct metabolic states benefits the cell community.

Results
Cells within S. cerevisiae colonies exhibit ordered metabolic specialization Using a well-studied S. cerevisiae isolate as a model (Reynolds and Fink, 2001), we established a simple system to study the formation of a clonal colony with irregular morphology. On 2% agar plates containing a complex rich medium with low glucose concentrations, S. cerevisiae forms rugose colonies with distinct architecture, after~5-6 days ( Figure 1A). Such colonies do not form in the typical, high (1-2%) glucose medium used for yeast growth ( Figure 1A). Thus, as previously well established (Granek and Magwene, 2010;Reynolds and Fink, 2001), glucose limitation (with other nutrients such as amino acids being non-limiting) drives this complex colony architecture formation. Currently, the description of such colonies is limited to this external rugose morphology, and does not describe the phenotypic states of cells and/or any spatial organization in the colony. With only such a description, as observed in Figure 1A, the mature colony surface has an internal circle and some radial streaks near the periphery. We carried out a more detailed observation of entire colonies under a microscopic bright-field (using a 4x lens). Here, we unexpectedly noticed what appeared to be distinct internal patterning, and apparent spatial organization of cells within the colony ( Figure 1B). As categorized purely based on these observed differences in visual optical density ('dark' or 'light'), regions between the colony center and periphery had optically dense (dark) networks spanning the circumference of the colony, interspersed with optically rare regions. In contrast, the periphery of the mature colony appeared entirely light ( Figure 1B). Based simply on these optical traits alone, we categorized cells present in these regions of the colony as dark cells and light cells ( Figure 1B). At this point, our description is visual and qualitative, and does not imply any other difference in the cells in either region. However, this visual description is both robust and simple, and hence we use this nomenclature for the remainder of this manuscript. Since these structured colonies form only in glucose-limited conditions, we hypothesized that dissecting the expected metabolic requirements during glucose limitation might reveal drivers of this internal organization within the colony. The expected metabolic requirements of cells growing in glucose limited conditions are as follows: first, all cells would be expected to have constitutively high expression of the high-affinity hexokinase (Hxk1p) (Lobo and Maitra, 1977;Rodríguez et al., 2001). Further, during glucose-limited growth, all cells are expected to carry out extensive gluconeogenesis, as the default metabolic state (Broach, 2012;Haarasilta and Oura, 1975;Yin et al., 2000). Indeed, we confirmed this second expectation by measuring the amounts of the gluconeogenic enzymes Pck1 (phosphoenolpyruvate carboxykinase) and Fbp1 (fructose-1,6-bisphosphatase), in short-term (4-5 hr) liquid cultures of log-phase cells growing in either high glucose medium, or in the same glucose-limited medium we used for colony growth. Expectedly, we observed very high amounts of these gluconeogenic enzymes in cells growing in glucose-limited medium (Figure 1-figure supplement 1A), reiterating that even in well-mixed glucose-limited, cells are in a strongly gluconeogenic state. In order to now examine the mature colony and dissecting expected metabolic requirements in these conditions, we first designed visual indicators for these metabolic hallmarks of yeast cell growth in low glucose. We engineered two different fluorescent reporters, one dependent on HXK1 expression (mCherry under the HXK1 gene promoter), and the second on PCK1 expression as an indicator of gluconeogenic activity (mCherry under the PCK1 gene promoter) (Figure 1-figure supplement 1B). Cells carrying these reporters were seeded to develop into colonies, and the expression levels of these reporters were monitored in the mature, rugose colony (5-6 days). Expectedly, the HXK1-promoter dependent reporter showed constitutive, high expression in all cells across the entire colony ( Figure 1C). Contrastingly, only the dark cells exhibited high gluconeogenesis reporter activity ( Figure 1C). Notably, the light cells entirely lacked detectable gluconeogenic reporter activity ( Figure 1C). To better quantify this phenomenon, cells were dissected out from dark or light regions respectively (under the light microscope, using a fine needle), and the percentage of fluorescent cells in each region was measured using flow cytometry. Based on flow cytometric readouts,~80% of the isolated dark cells showed strong fluorescence for the gluconeogenic reporter, while~97% of the light cells were non-fluorescent for gluconeogenic activity ( Figure 1C   observation was based solely on reporter activity, in order to more directly examine this observation, we estimated native protein amounts of enzymes associated with gluconeogenesis (Pck1, Fbp1, and Icl1-Isocitrate lyase from the glyoxylate shunt) in isolated light cells and dark cells. Only the dark cells showed expression of the gluconeogenic enzymes ( Figure 1D, Figure 1-figure supplement 2B). Finally, we measured steady-state amounts of trehalose and glycogen within dark and light cells, using these metabolites as unambiguous biochemical readouts of the end-point biochemical outputs of gluconeogenesis (Franç ois et al., 1991). We observed that the dark cells had substantially higher amounts of both trehalose and glycogen ( Figure 1E), indicating greater gluconeogenic activity in these cells. Collectively, these results strikingly reveal that intracellular gluconeogenic activity is spatially restricted to specific regions, resulting in a distinct pattern of metabolically specialized zones within the colony.

Cells organize into spatially restricted, contrary metabolic states within the colony
In the given nutrient conditions of low glucose, gluconeogenesis is an expected, constitutive metabolic process, essential for cells. This can therefore be considered as a necessary, permitted metabolic state in this condition. Paradoxically, in these mature colonies, gluconeogenic activity was spatially restricted to only within the dark cell region, with no discernible gluconeogenic activity in the cells located in the light region. This absence of gluconeogenic activity in these light cells, while concomitant with a constitutively high level of hexokinase activity, therefore poses a biochemical paradox. What might the metabolic state of these light cells be? To quickly address this using a crude but useful readout, we compared the ability of freshly isolated light and dark cells to proliferate in both gluconeogenic (low glucose), and non-gluconeogenic (high glucose) growth conditions. For simplicity, isolated light cells and dark cells were inoculated either into a medium where gluconeogenesis is essential (2% ethanol +glycerol as a sole carbon source), or in high (2%) glucose medium where cells rely on high glycolytic and pentose phosphate pathway (PPP) activity, and initial cell proliferation was monitored. Here, cells that had been growing in high glucose were used as a control. Expectedly, the dark cells grew robustly and reached significantly higher cell numbers (0D 600 ) compared to the light cells in the gluconeogenic condition ( Figure 2A). Conversely, light cells grew robustly when transferred to the high glucose medium, as compared to the dark cells ( Figure 2A). While this was an overly simple, and not definitive experiment, counter-intuitively, this result suggested that despite being in a low-glucose environment, the light cells were well suited for growth in high glucose, and therefore might be in a metabolic state suited for growth in glucose.
We therefore decided to more systematically investigate this phenomenon.
In the presence of glucose, yeast cells typically show high glycolytic and PPP activities, as part of the Crabtree (analogous to the Warburg) effect (Crabtree, 1929;De Deken, 1966; Figure 2B). Therefore, if the light cells in the colony were indeed behaving as though present in more glucosereplete conditions, they should exhibit high PPP activity. To test this, we first designed a fluorescent PPP-activity reporter (mCherry under the control of the transketolase 1 (TKL1) (Walfridsson et al., 1995) gene promoter,  analysis, using exogenously added 13 C Glucose, to compare flux of 13 C Glucose into the PPP metabolite ribulose-5-phosphate (R-5-P), in light and dark cells. The red circles represent 13 C labeled carbon atoms (n = 3). Also see Figure 2-figure supplement 1B. (E) Comparative metabolic-flux based analysis comparing 15 N incorporation into newly synthesized nucleotides, in dark and light cells. Also see S2, and Materials and methods. (F) Light cells and dark cells isolated from a 7 day old wild-type complex colony re-form indistinguishable mature colonies when re-seeded onto fresh agar plates, and allowed to develop for 7 days. Scale bar = 2 mm. Statistical significance was calculated using unpaired t test (*** indicates p<0.001, ** indicates p<0.01) and error bars represent standard deviation. The online version of this article includes the following figure supplement(s) for figure 2: shown as an overlaid quantitative heat-map histogram in Figure 2-figure supplement 1A. Next, we directly addressed the possibility of these light cells exhibiting relatively high PPP activity. For unambiguously testing this, we utilized a stable-isotope based metabolic flux experiment to assess the flux towards PPP in light and dark cells. Light and dark cells isolated from colonies were pulsed with 13 C-labeled glucose (for~5 min), metabolites extracted, and the incorporation of this carbon label into the late PPP intermediates ribulose-5-phosphate (R-5-P) and sedoheptulose-7-phosphate (S-7-P) was measured by liquid chromatography/mass spectrometry (LC/MS/MS). The relative amounts of these all-carbon labeled PPP intermediates were compared between the two cell types (light or dark). Notably, light cells incorporated significantly higher levels of 13 C labeled glucose into PPP metabolites compared to the dark cells ( Figure 2D, and Figure 2-figure supplement 1B, and see Table 1 for MS parameters), showing that the light cells are in a high PPP activity state. Finally, we assessed if other biochemical end-point outputs requiring high PPP activity/flux were also high in the light cells. High nucleotide synthesis is a canonical consequence of enhanced PPP activity   (Nelson and Cox, 2013). The carbon backbone (ribose-5-phosphate) of newly synthesized nucleotides is derived from the PPP, while the nitrogen backbone comes from amino acids (Nelson and Cox, 2013; Figure 2B, and see Table 1 for MS parameters). We devised another metabolic fluxbased experiment to assess de novo nucleotide biosynthesis in light and dark cells, as an end-point collective readout of high PPP activity coupled with amino acid utilization. Light and dark cells, isolated from colonies were pulsed with a 15 N-label (ammonium sulfate +aspartate), and incorporation of this label into nucleotides was measured by liquid chromatography/mass spectrometry (LC/MS/ MS). Light cells had higher flux into nucleotide biosynthesis, compared to the dark cells ( Figure 2E, and see Table 1 for MS parameters). Taken together, we surprisingly find that light cells exhibit multiple metabolic hallmarks of cells growing in glucose-replete conditions, including increased PPP activity, and increased nucleotide biosynthesis. Thus, in the spatially organized colony, the light cells and dark cells have contrary metabolic states. This is despite the expectation that the gluconeogenic state, exhibited by the dark cells, is the plausible metabolic state in the given growth conditions. Notably, the light cells or dark cells, when isolated and reseeded as a new colony, both develop into indistinguishable, complex colonies ( Figure 2F). This reiterates that these phenotypic differences between the light and dark cells are fully reversible, and do not require genetic changes. Collectively, these data reveal that cells within the colony organize into spatially separated, metabolically specialized regions. Within these regions, cells exhibit complimentary metabolic states. One of these states, where cells have high PPP activity, is counter-intuitive and cannot obviously be sustained given the external nutrient environment.

A mathematical model suggests constraints for the emergence and organization of cells in complimentary metabolic states
What determines the emergence and spatial organization of a group of cells, in these contrary metabolic states? Particularly, what can explain the emergence and proliferation of the light cells, which exhibit this counter-intuitive metabolic state, while the colony maintains a large subset of cells in the dark state? To address this, we built a coarse-grained mathematical model. This model incorporates  simple processes derived from our current experimental data, to simulate the formation of a colony of 'light' and 'dark' cells. The model was intentionally coarse-grained, since its purpose was only to find a minimal, biologically consistent combination of processes that is sufficient to produce the overall spatial structure and composition of cell states observed in the colonies. The intention behind the model was not to decipher all possible molecular details that explain this phenomenon. The model should only sufficiently account for both the emergence of light cells, as well as their spatial organization with dark cells. Such a model could therefore suggest constraints that determine the emergence of light cells, and the organization of the colony with the observed organization, which can then be experimentally tested.
While building this model, we included a range of processes that must be considered, based on our experimental data thus far ( Figure 3A). This includes (i) the dark cells switching to a light state, (ii) the production of some resource by dark cells, which may be shared/utilized by the cells, (iii) diffusion parameters for this resource, (iv) consumption of this resource, and (v) rates of cell division are included ( Figure 3A). Next, we constructed a two-dimensional square grid of 'locations' for groups of cells within the colony ( Figure 3B). Here, each location is either empty or occupied by a group of~100 cells (also see Materials and methods for full details). Note: we intentionally coarse-grain the grid (for computational simplicity, in order to simulate colony sizes comparable to real colonies) by approximating that the locations either consist of all light or all dark cells. This is a major simplification that was necessary. At each time step (12 min of real time), all the processes shown in Figure 3A are executed across the spatial grid using the outlined algorithm ( Figure 3C). In such an implemented algorithm, (i) all cells consume all available nutrients (present in saturating amounts), while free glucose concentrations are negligible, (ii) dark cells grow and divide in the given conditions, (iii) dark cells produce a resource/resources as a consequence of their existing metabolic (gluconeogenic) state, (iv) this resource diffuses around the grid and is freely available, (v) dark cells switch to the light state if sufficient resource is present at their location, and lastly, (vi) the resource when consumed can sustain the light state cells, which can expand if there is an empty location in the neighborhood. If the resource is not present in that location the light cells will switch back to dark cells. All processes occur at specified rates, allowing for stochasticity. Finally, this existence of a shared resource is surmised because, logically the emergence of light cells from dark can happen only if the local nutrient environment enables a switch to the new metabolic state.
In each simulation, empty grids are seeded with 1257 occupied locations, with 95-99% of the cells in the dark state. After~750 time steps (corresponding to~6 days) a simulated wild-type colony looks typically as shown in Figure 3D (also see Video 1). A range of resource amounts, growth rates and diffusion of the resource were included in control simulations (see Figure 3-figure supplement 1A-C). Strikingly, the simulated spatial organization ( Figure 3D, and Video 1) recapitulates most obvious features of a real colony ( Figure 3D). These are: at the edge of the initial circular inoculation of the colony is a ring of dark cells, the outermost part of the colony is made up of outcrops of light cells, and from this ring of dark cells emanate clusters of dark cells penetrating into the outcrops of light cells. This is despite the simplicity of the rules in the model, including its flattening into 2D. In the simulation, for the first 40-45 time steps, the colony remains small and predominantly dark, while the resource builds up. Then, dark cells start to switch to light. When this happens within the bulk of the colony, these light cells have restricted division due to spatial constraints. Around 100 to 150 time steps later, light cells emerge at the perimeter of the colony, and then rapidly divide and expand ( Figure 3D, Video 1). In order to test if the processes of Figure 3A are all required for this behavior, we examined three comprehensive control scenarios: (i) dark cells do not produce a resource (and therefore in this case for growth light cells depend only on amino acids or other pre-supplied resources in the medium), (ii) dark cells cannot switch to the light state, or (iii) light cells produce a resource that is needed by dark cells to grow (a straw-man scenario, since initially in the actual colony all cells were in a dark state, as shown earlier). None of these cases produces the wildtype spatial organization, over a wide range of parameter values ( Figure 3E, as well as a range of parameters explored in Figure 3-figure supplement 2A-D, and simulations in Videos 2-4).
Summarizing, this simple model successfully recaptures the general features of the spatial patterning and organization of real colonies. This includes the overall general architecture, and spatial organization of light and dark cells. Two simple take-home points emerge from this model, for such spatial distribution of cells in these two metabolic states, across the developing colony. First, the model requires that dark (gluconeogenic) cells will produce a resource that is needed by dark cells to switch to the light state. Second, a resource produced by the dark cells is required to sustain the light state. Collectively, in our model, these metabolic constraints are sufficient to determine the overall spatial organization of metabolically distinct, specialized cells.
Trehalose satisfies criteria to be the metabolic resource determining the emergence of light cells Does any gluconeogenic metabolite(s) determine the organization of these cells, consistent with Video 3. Dark cells do not switch to light cells. Simulation video showing changes in a colony where the dark cells don't switch to light cells but continue to produce resource on to the resource grid. In certain cases, there can be light cells at the very edge of the starting colony. This is because the composition of the colonies might be the same between simulations, but individual cell block locations are done at random. Such cells at the edge can utilize the shared resource and divide into empty space on the grid. https://elifesciences.org/articles/46735#video3 Video 4. Light cells share metabolic resource with dark cells (Wrong sharing). Simulation video showing changes in a colony where the dark cells don't share any resource for the light cells but the light cells provide amino acids for the dark cells to consume (wrong sharing). The dark cells have an abundance of amino acids to grow and divide. The final colonies predominantly comprise of dark cells. https://elifesciences.org/articles/46735#video4 these requirements suggested by experimental and modeled data? Such a metabolite must logically satisfy the following three criteria. First, this resource should be available in the extracellular environment (i.e. released by cells), second, cells must selectively be able to take up this resource, and third, the resource should be metabolized within cells to produce glucose/a glucose-like product capable of fueling a glycolytic and PPP-active state. Further, if this were indeed a 'controlling resource' that determined the emergence of light cells, preventing the uptake and utilization of this resource should prevent the emergence and proliferation of only the light cells, but leave the dark cells unaffected. In order to identify such a candidate metabolite, we considered all possible outputs of  gluconeogenesis: the storage carbohydrates/sugars glycogen and trehalose, the polysaccharides of the cell wall (chitin, mannans, glycans), and glycoproteins ( Figure 4A) (Jules et al., 2008;Kayikci and Nielsen, 2015). The large molecular size of glycogen, chitins, and complex glycosylated proteins, the lack of known cellular machinery for their uptake, and the difficulty in efficiently breaking them down make them all unlikely candidates to be the resource controlling the emergence of light cells. Contrastingly, trehalose has unique properties making it a plausible candidate. It is a small, non-reducing disaccharide composed of two glucose molecules. Trehalose has been observed in the extracellular environment in yeast (Parrou et al., 2005), and yeast can uptake trehalose through disaccharide transporters (Jules et al., 2008;Stambuk, 1998). Further, trehalose can be rapidly and specifically hydrolyzed to two glucose molecules, which can fuel glycolysis and re-entry into the cell division cycle (Laporte et al., 2011;Shi et al., 2010;Shi and Tu, 2013). These diverse data therefore presented trehalose as an excellent putative candidate metabolite that controlled the emergence of cells in the light state. To test this possibility, we first measured extracellular trehalose in colonies. Free trehalose was readily detectable in the extracellular environment of these colonies ( Figure 4B). To test if trehalose could be differentially transported into either light or dark cells, we first estimated amounts of a primary trehalose transporter, Mal11 (Stambuk, 1998) in these cells. Mal11 protein amounts were substantially higher in the light cells compared to the dark cells ( Figure 4C, and quantified in Figure 4-figure supplement 1A). To unambiguously, directly estimate trehalose uptake, we isolated light and dark cells from a mature colony, and exogenously added 13 C-trehalose. We then measured intracellular levels of labeled trehalose present in these cells, by extracting and estimating metabolite amounts (by LC/MS/MS) (see Table 1 for MS parameters). Notably, the light cells rapidly accumulated 13 C-trehalose ( Figure 4D), while the dark cells did not, suggesting robust, preferential uptake of extracellular trehalose. Finally, we estimated the ability of light and dark cells to break-down and utilize trehalose. For this, we first measured the expression of the predominant neutral trehalase in yeast (Nth1) (Jules et al., 2008), in the light and dark cells. Light cells had substantially higher Nth1 amounts than the dark cells ( Figure 4E, and quantified data shown in Figure 4-figure supplement 1B). We also measured enzymatic activity for Nth1 (in vitro, using cell lysates), and found that the light cells had~2 fold higher in vitro enzymatic activity, compared to the dark cells ( Figure 4F). Collectively, these data suggested that the light cells were uniquely able to preferentially take up more trehalose, break it down to glucose, in order to potentially utilize it to sustain a metabolic state with high PPP activity.

Trehalose uptake and utilization determines the existence of light cells
Since these data suggested that trehalose uptake and utilization would be preferentially high in the light cells, we directly tested this using a quantitative metabolic flux based approach. For this we used stable-isotope labeled trehalose, and measured trehalose uptake, breakdown and utilization. To the isolated light and dark cells, 13 C-labeled trehalose was externally provided, and intracellular metabolites extracted from the respective cells. The intracellular amounts of 13 C -labeled glycolytic and PPP intermediates were subsequently measured using LC/MS/MS ( Figure 5A and  Table 1 for MS parameters). 13 C -labeled glucose-6-phosphate (which enters both glycolysis and the PPP), the glycolytic intermediates glyceraldehyde-3-phosphate and 3phosphoglycerate, and the PPP intermediates 6-phosphogluconate, ribulose-5-phosphate and sedoheptulose-7-phosphate all rapidly accumulated exclusively in the light cells ( Figure 5A, Figure 5figure supplement 1A). Since the labeled carbon can come only from trehalose, these data indicate both the breakdown of trehalose to glucose, as well as the subsequent utilization of glucose for these pathways. Indeed, the labeled forms of these metabolites were only above the detection limit in dark cells ( Figure 5A and Figure 5-figure supplement 1A). Thus, these data demonstrate that external trehalose is preferentially taken up only by the light cells, and utilized to fuel the complimentary metabolic state of the light cells, with high glycolytic and PPP activity.
Finally, we tested if the sharing and differential utilization of trehalose determined both the emergence and the proliferation of light cells. An explicit prediction is made both in our model, and our hypothesis based on these experimental data. This is: preventing uptake and/or utilization of trehalose should prevent cells from switching to the light state. To test this prediction, we generated strains lacking NTH1 (which cannot break-down trehalose to glucose), and MAL11 (which will have reduced trehalose uptake), allowed them to develop into mature colonies, and compared the amounts of light cells in each. Compared to wild-type colonies, cells lacking the major trehalose uptake transporter (Dmal11) formed colonies with very few light cells ( Figure 5B). Note: while Mal11 shows a high affinity for trehalose, S. cerevisiae has multiple sugar transporters with reduced affinity for any disaccharide. Therefore, cells lacking MAL11 may take up trehalose with lower efficiency. In these cells, the ability to break-down trehalose remains intact. More importantly, in colonies of cells lacking trehalase (Dnth1), and which therefore cannot efficiently breakdown internal trehalose, had nearly no detectable light cells, based on brightfield microscope observations ( Figure 5B). This result was more quantitatively estimated in colonies of cells with these respective genetic backgrounds, using the expression of the fluorescent PPP reporter. Again, almost no PPP reporter activity was observed in the Dnth1 cell colonies, while very few cells with PPP reporter activity were seen in Dmal11 colonies ( Figure 5C). As controls, we ensured that there were no defects in the expression of the PPP reporters in cells from these genetic backgrounds. Correspondingly, we also quantified the percentage of dark, highly gluconeogenic cells (as determined using the gluconeogenesis reporter), in colonies from each of these genetic backgrounds. The percentage of gluconeogenic cells was proportionately higher in the Dmal11 (~73%), and Dnth1 (~80%) colonies compared to the wild-type colony (~65%) ( Figure 5D). Thus, controlling the uptake and utilization of the resource (trehalose) directly regulates the emergence of cells in the light state.
Collectively, these data demonstrate that trehalose is the shared gluconeogenic resource that determines the emergence, and persistence, of light cells within the structured colony.

A resource threshold effect controls cooperative switching of cells to the light state
Our experimental data showing the organization of dark and light cells was obtained from~5-6 day old, mature colonies. However, in our simulations of the temporal development of the colony, we observed that the dense network of dark (gluconeogenic) cells form first, followed by a very late appearance of light cells ( Figure 6A and Video 1). This late appearance of light cells in the simulations comes from an inherent threshold effect included within the model. Here, the external build-up of the shared resource made by the dark cells is required. At a sufficient built-up concentration, this resource will trigger the switching of some cells to light cells. Light cells in turn will consume the resource, reducing the available amounts, thereby preventing other cells from switching to this new state. This threshold-effect therefore predicts a delayed, rapid emergence of light cells, and also enables such a pattern of distinct cell groups to form. If this threshold requirement is removed in the simulation (for example when replaced by a rate of switching from dark to light that depends linearly on the amount of resource), the resultant colony remains small, and the organized pattern of cells in two states does not occur. This is shown in Figure 6A, and Video 5. This small colony size is largely due to low resource amounts to support the proliferation of the light cells, since there are insufficient dark cells remaining to produce the resource. This is also clearly seen in control simulations with a range of resource amounts, and linear switching, as shown in Figure 6-figure supplement 1A-1D. Contrastingly, in the model that successfully simulates the colony development, the externally available amount of the resource builds-up, reaches the threshold (where cells switch to the light state), and then rapidly decreases, if the light cells also consume the resource ( Figure 6B, upper panel).
This therefore prompted us to more closely examine the development of actual colonies over time, for these properties. We first estimated the relative amounts of extracellular, free trehalose in Video 5. Development of colony without a resource threshold. Simulation video showing changes in a colony where the dark cells switch to light by random chance (probability p=0.5). They don't need the resource levels to reach a certain threshold. Once they become light cells, they cannot switch back to dark. Due to lack of the produced resource, the colony doesn't grow discernibly.
https://elifesciences.org/articles/46735#video5 the colony over time. Notably, trehalose amounts steadily increased over 4 days, and subsequently rapidly decreased ( Figure 6B lower panel). This rapid decrease in trehalose after day four is despite a steady, continuing increase in the total number of cells in the colony (Figure 6-figure supplement 1E). We next monitored the development of colonies over time, to determine when the light cells emerge. Using just the bright-field image reconstruction of the colonies, during this time course, the intensity of dark cells steadily increased, and organized into the mesh-like network over 4 days ( Figure 6C). However, the light cells appeared only after~4 days, and rapidly increased in number ( Figure 6C). We more quantitatively estimated this, using strains expressing the  Figure 1A) and Dnth1 cells measured as a function of their ability to spread on a plate. Colony spreading was quantified by measuring the diameter of the colonies every day for 21 days (n = 3). (B) Cells in low glucose perform gluconeogenesis, as required in low glucose medium. As gluconeogenic reserves build up, trehalose builds up in the extracellular environment. At a threshold concentration of trehalose, some cells switch to a high glycolytic, PPP state. This state depends upon the utilization of trehalose to fuel it. This utilization of trehalose by the light cells results in decreased external trehalose to below a threshold. This in turn restrains the other, remaining cells in a gluconeogenic state, where they continue to produce trehalose. This gives rise to the final, self-organized community, with specialization of function and division of labor. Statistical significance was calculated using unpaired t test (** indicates p<0.01) and error bars represent standard deviation.
gluconeogenic-or the PPP-reporter ( Figure 6D). Notably, the increase in total fluorescence intensity due to the gluconeogenic-reporter in the colony (over time) was relatively linear over the first four days (r 2 = 0.99), increasing with the steady increase in the number of cells (Figure 6-figure supplement 1E). Contrastingly, the increase in the PPP reporter activity over the first five days was clearly non-linear and switch like, with very low signal intensity for the first three days, and then a rapid emergence of signal over days four and five ( Figure 6E). This indicated a cooperative, switch like emergence of, and increase in these light cells. A useful biophysical measure of cooperativity (more commonly used for protein-ligand binding characteristics) is the Hill coefficient. We adopted the Hill equation, using the amount of PPP reporter fluorescence (instead of ligand-receptor binding), to estimate cooperativity in the system. Over the first five days the increase in PPP-reporter activity showed a Hill coefficient greater than 1, indicating a positively cooperative switch of cells to the light state ( Figure 6E). This nicely correlates with the build-up, and rapid decrease in external trehalose around day 4 ( Figure 6B). These data also show that the peripheral location of the light cells cannot simply be due to possible greater access to glucose in the medium, since for the first~4 days, there are no cells in the periphery with high PPP reporter activity. Their emergence is indeed rapid, and switch-like.
In summary, data from model simulations and experiments show that initially the gluconeogenic cells increase in number, leading to release and build-up of the resource (trehalose) in the local environment. At this time there are no light cells in the colony. Above a threshold concentration of trehalose, some cells rapidly switch to light state with high PPP activity. The further expansion of these light cells correlates with rapid consumption of the extracellular trehalose that sustains this state. These data suggest a threshold effect, where the controlling resource, trehalose, needs to build up above a certain amount, in order for cells to switch to the contrary, high PPP activity state.
Cells in distinct metabolic states provide a collective growth advantage to the colony Finally, we wondered if such an emergence of light cells with high PPP activity might benefit the community of cells as a whole. To address this, we compared the long term colony expansion of wild-type cells, with colonies comprised of cells lacking the neutral trehalase (Dnth1). Cells in the Dnth1 colonies cannot utilize trehalose to produce glucose, and as shown earlier, will remain in a gluconeogenic state. Therefore, in these Dnth1 cell colonies light cells will be absent. However, these cells are still capable of normal gluconeogenesis (and trehalose production). Strikingly, we observed that as the respective colonies expanded over time (~21 days), the wild-type colonies spread over a significantly greater area on the plate, while the Dnth1 colonies were unable to expand as efficiently ( Figure 7A and Figure 7B). This shows that the emergence and proliferation of light cells are important for the expansion of the colony. Since the dark cells are required for the emergence and existence of the light cells, collectively, these data suggest how the community uses cells in distinct metabolic states to maximize growth and spatial expansion, possibly to forage for new nutrients.

Discussion
Collectively, we present a simple model proposing how cells in metabolically distinct states spontaneously emerge and spatially self-organize within a yeast colony, as summarized in Figure 7C. In low glucose conditions, cells begin in a uniform gluconeogenic state, which is the expected metabolic state in this nutrient condition. The gluconeogenic cells produce a resource (trehalose), that is now externally available. This resource builds up to above a threshold amount. At this threshold, some cells take up and consume trehalose, breaking it down to glucose. These cells spontaneously switch to the complimentary metabolic state, with high PPP and glycolytic activity (i.e. the light state) ( Figure 7C). These light cells can remain in this metabolic state only so as long as the resource (trehalose) is externally available. However, as trehalose is consumed by these cells, the available amount of external trehalose itself drops below the threshold. The surrounding dark cells therefore remain trapped in a gluconeogenic state, continuing to produce trehalose. Thereby, a predictable fraction of cells, constrained spatially, will remain in each metabolic state, resulting in specialized cell groups and division of metabolic labor. Thus, biochemically heterogeneous cell states can spontaneously emerge and spatially self-organize. An implicit concept emerging from this study is that of threshold amounts of a controlling or sentinel metabolite that regulates a switch to a new metabolic state. By definition, such a metabolite must be produced by cells present in a certain (original) metabolic state. But when this metabolite is utilized, it must have the ability to switch cells to an entirely distinct metabolic state. Further, the emergence and expansion of cells in the new state will be rapid and switch-like, resembling a bistable system (Pomerening, 2008). This idea of metabolites controlling cell states is an emerging area of interest (Cai and Tu, 2011;Krishna and Laxman, 2018), but has not been studied in the context of groups of cells organizing into distinct groups or metabolic states (and therefore different phenotypic properties).
We speculate what the advantages of such spatially organized, phenotypically distinct states within a group of clonal cells might be, in this example. Here, we observe that the organized community with cells in distinct states has clear advantages, in being able to spatially expand better ( Figure 7C). For sessile microbes such as yeast, this ability to forage for better nutrients is important for their survival. This might also convey other advantages, and uncovering those are obvious areas of future studies. Since the inherent properties of the cells in the distinct states are different, this raises the deeper possibility that these advantages come from physical and chemical properties of the cells, which arise from their distinct metabolic states. Regardless, our study substantially advances descriptions of yeast 'multicellularity' from simple dimorphism, aggregated cells, or three-dimensional colony forms (Cáp et al., 2012;Koschwanez et al., 2011;Palková and Váchová, 2016;Ratcliff et al., 2012;Váchová and Palková, 2018;Veelders et al., 2010;Wloch-Salamon et al., 2017), to self-organized, phenotypically heterogeneous cell states exhibiting division of labor and metabolic interdependence. Strikingly, the nature of spatial patterning allied with division of labor that we observe in yeast is reminiscent of true multicellular systems (Newman, 2016;Niklas, 2014). Also, the cell states in these yeast colonies can be considered commensal, since trehalose is a necessary output of gluconeogenesis, and therefore a default, biochemically non-limiting output in dark cells. Since trehalose controls the emergence and maintenance of light cells in the complimentary metabolic state, it thus can be considered a resource benefiting the light state. Thus simple, metabolism-derived constraints are sufficient to determine how contrary biochemical states can spontaneously emerge and be supported, in conjunction with spatial structure. Such organization of cells into specialized, labor-divided communities expands on the role of reaction-diffusion systems (particularly activator-depleted substrate schemes) in controlling cellular patterning (Gierer and Meinhardt, 1972;Kondo and Miura, 2010;Newman, 2016), with a metabolic resource threshold being central to the emergence and stabilization of a new phenotype (Cai and Tu, 2011;Krishna and Laxman, 2018). A deeper dissection of what such constraints can permit will therefore advance our general understanding of how specialized cell states can emerge and be stabilized.
Metabolic cross-feeding is best understood currently in multi-species microbial communities, where this has been inferred largely using inter-species genomic comparisons (Ackermann, 2015;D'Souza et al., 2018;Goldford et al., 2018;Tyson et al., 2004). Further, metabolic sharing has typically been demonstrated using synthetically engineered systems where mutual dependencies are created (Campbell et al., 2016;D'Souza et al., 2018;Mee et al., 2014;Pande et al., 2015;Wintermute and Silver, 2010). The spatial organizations of any such populations remain challenging to model. Biochemically identifying metabolites that are conclusively exchanged between cooperating cells remains difficult, and the significance of such putative metabolite exchange challenging to interpret (Ackermann, 2015;D'Souza et al., 2018). Finally, such studies have emphasized non-isogenic systems, where genetic changes stabilize different phenotypes, and auxotrophies define the nutrient sharing or cooperation (Ackermann, 2015). Contrastingly, here we directly identify a produced metabolic resource, and demonstrate how its availability and differential utilization can control the emergence of cells in opposing metabolic states, in a clonal population. We also explain how the spontaneous spatial organization into phenotypically distinct cell groups occurs. Thus, our study also goes beyond stochastic gene expression (Ackermann, 2015;Balázsi et al., 2011;Blake et al., 2003) to explain how phenotypic heterogeneity and specialization can emerge in clonal populations. By considering these metabolism-derived rules, and thereby manipulating available metabolic resources, we suggest how it can be viable to program the formation, structure or phenotypic composition of isogenic cell populations. Collectively, such simple physico-chemical constraints can advance our understanding of how isogenic cells can self-organize into specialized, labor-divided groups, as a first step towards multicellularity.

Yeast strains and growth media
The prototrophic sigma 1278b strain (referred to as wild-type or WT) was used in all experiments. Strains with gene deletions or chromosomally tagged proteins (at the C-terminus) were generated as described (Longtine et al., 1998). Strains used in this study are listed in Table 2. The growth medium used in this study is rich medium (1% yeast extract, 2% peptone and 2% glucose or 0.1% glucose).

Colony spotting assay
All strains were grown overnight at 30˚C in either rich medium or minimal medium. 5 microliters of the overnight cultures were spotted on rich medium (low glucose) (1% yeast extract, 2% peptone, 0.1% glucose and 2% agar). Plates were incubated at 30˚C for 7 days unless mentioned otherwise.

Colony imaging
For observing colony morphology, colonies were imaged using SZX-16 stereo microscope (Olympus) wherein the light source was above the colony. Bright-field imaging of 7 day old colonies were done using SZX-16 stereo microscope (Olympus) wherein the light source was below the colony. Epifluorescence microscopy imaging of 7 day old gluconeogenesis reporter colonies (pPCK1-mCherry), pentose phosphate pathway (PPP) reporter colonies (pTKL1-mCherry) and HXK1 reporter colonies (pHXK1-mCherry) were imaged using the red filter (excitation of 587 nm, emission of 610 nm) of SZX-16 stereo microscope (Olympus). Similar protocol was followed for imaging 1 day to 6 day old colonies.

Analysis of fluorescent cell populations in reporter strain colonies
Light cells and dark cells isolated from 7 day old wild-type colonies harboring either the gluconeogenesis reporter, PPP reporter or the HXK1 reporter were re-suspended in 1 ml of water. The percentage of fluorescent cells were determined by running the samples through a flow cytometer, and counting the total number of mCherry positive cells in a total of 1 million cells. Light cells and dark cells isolated from wild-type colonies without the fluorescent reporter were used as control.
Biochemical estimation of trehalose/glycogen levels Trehalose and glycogen from yeast samples were quantified as described previously, with minor modifications (Shi et al., 2010). 10 OD 600 of light cells and dark cells from 7 day old wild-type colonies (rich medium, 0.1% glucose) were collected. After re-suspension in water, 0.5 ml of cell suspension was transferred to four tubes (two tubes for glycogen assay and the other two tubes for trehalose assay). When sample collections were complete, cell samples (in 0.25 M sodium carbonate) were boiled at 95-98˚C for 4 hr, and then 0.15 ml of 1 M acetic acid and 0.6 ml of 0.2 M sodium acetate were added into each sample. Each sample was incubated overnight with 1 U/ml amyloglucosidase (Sigma-Aldrich) rotating at 57˚C for the glycogen assay, or 0.025 U/ml trehalase (Sigma-Aldrich) at 37˚C for the trehalose assay. Samples were then assayed for glucose using a glucose assay kit (Sigma-Aldrich). Glucose assays were done using a 96-well plate format. Samples were added into each well with appropriate dilution within the dynamic range of the assay (20-80 mg/ml glucose). The total volume of sample (with or without dilution) in each well was 40 microliters. The plate was pre-incubated at 37˚C for 5 min, and then 80 ml of the assay reagent from the kit was added into each well to start the colorimetric reaction. After 30 min of incubation at 37˚C, 80 microliters of 12 N sulfuric acid was added to stop the reaction. Absorbance at 540 nm was determined to assess the quantity of glucose liberated from either glycogen or trehalose. For measurement of extracellular trehalose measurement, single wild-type colony (1 day to 7 day old colony) was re-suspended in 100 microliters of water and centrifuged at 20000 g for 5 min. Supernatant was collected and buffered to a pH of 5.4 (optimal for trehalase activity) using sodium acetate buffer (pH 5.0). 0.025 U/ml trehalase (Sigma-Aldrich) was added and samples were incubated at 37˚C overnight. Glucose concentration was estimated as described earlier.

Neutral trehalase activity assay
Neutral trehalase activity assay was performed as described earlier with the following modifications (De Virgilio et al., 1991). Briefly, 1 OD 600 of light cells and dark cells isolated from 7 day old wildtype colonies (rich medium, 0.1% glucose) were washed twice with ice-cold water. For permeabilization, cells were re-suspended in tubes containing equal volume of 1% Triton-X in assay buffer (200 mM tricine buffer (Na + ) (pH 7.0)) and immediately frozen in liquid nitrogen. After thawing (1-4 min at 30˚C), the cells were centrifuged (2 min at 12000 g), washed twice with 1 ml of ice-cold assay buffer and immediately used for the trehalase assay. Trehalase assay was performed in 50 mM tricine buffer (Na + ) (pH 7.0), 0.1 M trehalose, 2 mM manganese chloride (MnCl 2 ) and the Triton X-100 permeabilized cells in a total volume of 400 microliters. After incubation for 30 min at 30˚C, the reaction was stopped in a boiling water bath for 3 min. Glucose concentration in the supernatant was determined using the glucose assay kit (Sigma-Aldrich).

Western blot analysis
Approximately ten OD 600 cells were collected from respective cultures, pelleted and flash frozen in liquid nitrogen until further use. The cells were re-suspended in 400 microliters of 10% trichloroacetic acid (TCA) and lysed by bead-beating three times: 30 s of beating and then 1 min of cooling on ice. The precipitates were collected by centrifugation, re-suspended in 400 microliters of SDS-glycerol buffer (7.3% SDS, 29.1% glycerol and 83.3 mM tris base) and heated at 100˚C for 10 min. The supernatant after centrifugation was treated as the crude extract. Protein concentrations from extracts were estimated using bicinchoninic acid assay (Thermo Scientific). Equal amounts of samples were resolved on 4% to 12% bis-tris gels (Invitrogen). Western blots were developed using the antibodies against the respective tags. We used the following primary antibody: 538 monoclonal FLAG M2 (Sigma-Aldrich). Horseradish peroxidase-conjugated secondary antibody (anti-mouse) was obtained from Sigma-Aldrich. For Western blotting, standard enhanced chemiluminescence reagents (GE Healthcare) were used.

N-and 13 C-based metabolite labeling experiments
For detecting 15 N label incorporation in nucleotides, 15 N Ammonium sulfate (Sigma-Aldrich) and 15 N Aspartate (Cambridge Isotope Laboratories) with all nitrogens labeled were used. For 13 C-labeling experiment, 13 C Trehalose with all carbons labeled (Cambridge Isotope Laboratories) was used. All the parent/product masses measured are enlisted in Table 1. For all the nucleotide measurements, release of the nitrogen base was monitored in positive polarity mode. For all sugar phosphates, the phosphate release was monitored in negative polarity mode. The HPLC and MS/MS protocol was similar to those explained above.

Building and implementing the model
The model simulation code is available via GitHub ref: https://github.com/vaibhhav/ yeastmetabcolony. These processes implement growth of cells, as well as production and consumption of amino acids and the shared metabolic resource. Subsequent to this, in each time step, we allow diffusion of the resource levels across the grid (the "external" level at the location, not the internal levels in cell blocks), using a numerical scheme called Forward Time Central Space (FTCS). Say that the value of the resource at time t and location (x,y) is given by U t x;y . The FTCS scheme updates the value simultaneously at all locations using the following formula: where DT is the time step and DL is the space step, or grid length, and D is the diffusion constant for the resource.

Model parameters
1. The parameters of the model are shown in Table 3. Time and length units are chosen such that each time step is one unit of time, and each grid length is one unit length. With these choice of units, the growth parameters for light and dark cells, respectively, are g l = 0.04, g d = 0.01. These were chosen so as to reflect the relative rates of diffusion and division. Light cells were observed to grow faster than the dark cells, so their respective growth parameters are set accordingly. 2. The switching threshold parameter (S = 3.0) was chosen to account for a delay in the switching of dark cells to another metabolic pathway via nutrient sensing as well as to give a reproducible facsimile of the experimental colonies. 3. The shared resource production value was chosen to be 7% (R = 0.07) of the minimum required to divide. In each time step, every block of dark cells adds this amount to the resource grid. This was chosen as a default value, which gave a reproducible facsimile of the experimental colonies. Other values were tried and their effect is seen in Figure 3-figure supplement 1C. 4. All cells consumed a small level of metabolites (the shared resource or amino acids) in each time step. This value was chosen to be 5% of the minimum required for division (C = 0.05). This gave division times that approximately matched the division times observed experimentally. 5. The switching probability (p=0.5) was chosen to add an element of stochasticity. So even if the threshold resource conditions are met, dark cells have a 50% chance to switch to light cells in that time step. 6. The choice of the diffusion constant (D = 0.24) is limited by the numerical stability of the FTCS scheme, which allows only a maximum value of D = 0.25. In real time and length units, this corresponds to a diffusion constant D eff of 8.7 Â 10 À13 m 2 /s. D eff is an order of magnitude smaller than the diffusion constant for sugars like glucose and sucrose in water (Roache, 1972). Since the agar used for the experiments is mostly water, the diffusion constants in water can be considered as a good reference point. Figure 3E, Figure 6A (bottom row) and Figure 6-figure supplement 1C-D showcase some of the final colonies generated by the simulations when the rules described above are varied. The following changes were made to the rules/parameters to generate these. 3E(i): No sharing: Set R = 0. 3E(ii): No switching from dark to light state: Set p=0. 3E(iii): 'Reverse' sharing: Set R = 0. When a cell block is light it adds R'=0.07 to the amino acid grid at the same location.

Variants of the wild-type model in different figures
6A(bottom row): No resource thresholding: Set S = 0. Figure 6-figure supplement 1C-D: Linear switching: Set S = 0. The probability of switching from dark to light state, p, is now a linear function of the locally available resource with a maximum value of 1.0. That is, p = max (m Ã U t x;y , 1.0), where m is a parameter that sets the slope of this linear relationship.

Data availability
Mathematical model data in this study are generated by computational simulations. All model parameters and equations are included in the text, and source code is included with this article.