Introduction

Polysaccharides represent the largest fraction of biomass on Earth1,2 and are constantly degraded and remineralized by microorganisms. Polysaccharides, also known as glycans, are long chains of monosaccharide units produced by cells for structural support (e.g., cellulose, chitin, and alginate) or energy storage (e.g., starch or glycogen)1. Heterotrophic microbes obtain nutrients and energy from the polysaccharide breakdown products. They often use exoenzymes, either secreted or anchored in the cell membrane2,3, to cleave these large polymers into smaller units that can be taken up by cells. The formation of dense cell groups is observed during growth on polysaccharides by diverse bacteria, including the soil- and gut-dwelling Bacillus subtilis3, the oligotrophic fresh-water Caulobacter crescentus4, and several representatives of the copiotrophic5,6 marine Vibrio spp.7,8. It has been suggested that “cooperative” growth3, where dense cell groups reduce diffusional loss of valuable degradation products and secreted exoenzymes, represents a general principle in polysaccharide degradation4,8. Yet, how bacteria subsequently disperse from exhausted polysaccharide sources and navigate towards new polysaccharide hotspots remains poorly understood.

To study the cellular mechanisms that govern the switch between polysaccharide degradation and dispersal toward new polysaccharide hotspots, we worked with the ubiquitous marine Gammaproteobacterium Vibrio cyclitrophicus ZF270, a degrader of the prevalent marine polysaccharide alginate6,9. Alginate is a linear polysaccharide that is produced by brown algae as a cell wall component, as well as by certain bacteria. It can be cleaved into oligomers by endo-acting alginate lyases or into monomers by exo-acting alginate lyases (alginate lyases reviewed here10), namely β-D-mannuronic acid and α-L-guluronic acid11. V. cyclitrophicus ZF270 is found predominantly in the large-particle fractions of coastal water12, can attach and form biofilms on the surfaces of multiple polysaccharides including alginate13, and was shown to secrete alginate lyases during alginate degradation8. In this study, we used microfluidics coupled to automated time-lapse imaging to quantify the growth dynamics, group formation, and motility of V. cyclitrophicus ZF270 at the single-cell level under constant supply of either polymeric alginate or alginate degradation product in the form of digested alginate, as well as upon a transition from alginate to digested alginate. Furthermore, we used RNA-sequencing to compare the gene expression of V. cyclitrophicus ZF270 grown on alginate and digested alginate. We found striking responses to the form of alginate in growth rate, group formation, motility and chemotaxis, as well as in the expression of corresponding genes. Overall, our work provides insights into the metabolic and cellular regulation that allows cells to forage in heterogeneous nutrient-scapes through degradation-dispersal cycles.

Results

Extracellular break down of alginate delays population growth

To probe the phenotypic and metabolic regulation of bacterial cells during the progressive degradation of a polysaccharide source, we developed an experimental system consisting of the marine bacterium V. cyclitrophicus ZF270 growing on alginate, a highly abundant polysaccharide in the ocean. To mimic the local nutrient environment during the colonization of a new polysaccharide source, we used 0.1% weight per volume (w/v) algae-derived alginate in its soluble form (in the following also simply referred to as “alginate”). Since commercially available alginate breakdown products of specific sizes are limited and expensive, with only one monomeric component available at a considerable cost, we simulate an advanced stage of polysaccharide degradation by supplying cells with digested alginate (0.1% (w/v)). This digested alginate is prepared by treating alginate with a readily available endo-acting alginate lyase (see Materials and Methods). This mimics an environment where degradation products like monomers and oligomers become abundant through the action of extracellular alginate lyases. The commercially available alginate lyase has been shown to produce alginate oligomers of progressively smaller size over extended digestion periods14. We used Liquid Chromatography-Mass Spectrometry (LC-MS) to analyze the composition of the digested alginate, which had been subject to 48 hours of digestion. We found that digested alginate contained more monosaccharides than untreated alginate, which contained more alginate molecules of higher molecular weight (Fig. S1).

Using this system, we first set out to investigate the growth dynamics of V. cyclitrophicus ZF270 in well-mixed batch cultures containing either alginate or digested alginate as a sole carbon source (Fig. S2A). Measurements of optical density showed that the onset of growth on alginate was delayed by about 7.5 hours compared to growth on digested alginate (Fig. S2B, D and E). This is consistent with previous observations that in well-mixed environments the lag time of bacteria growing on alginate can be reduced by the external supplementation of alginate lyases8. We also observed that the optical density at stationary phase was higher when cells were grown on alginate (Fig. S2B and C). However, colony counts did not show a significant difference in cell numbers (Fig. S3), suggesting that the increased optical density may stem from aggregation of cells in the alginate medium, as observed for other Vibrio species7. Overall, these findings indicate that growth on polysaccharides such as alginate in well-mixed cultures is initially limited by the extracellular breakdown of the polysaccharide, though similar cell numbers were reached eventually.

Large cell groups form on alginate but not on digested alginate

To better understand how cells react to the changing degree of depolymerization of a polysaccharide source during degradation, we investigated the growth of V. cyclitrophicus ZF270 at the level of single cells on alginate and digested alginate. For this purpose, we grew the cells in microfluidic growth chambers, as depicted in Figure 1A and described in detail by Dal Co et al.15. In brief, the microfluidic chips are made of an inert polymer (polydimethylsiloxane) bound to a glass coverslip. The PDMS layer contains flow channels through which the culture medium is pumped continuously. Each channel is connected to several growth chambers that are laterally positioned. The dimensions of these growth chambers (height: 0.85 µm, length: 60 µm, width: 90-120 µm) allow cells to freely move and grow as monolayers. The culture medium, containing either alginate or digested alginate in their soluble form, is constantly pumped through the flow channel and enters the growth chambers primarily through diffusion15,16,4,17,8. Therefore, the number of cells and their positioning within microfluidic chambers is determined by the cellular growth rate as well as by cell movement4. This setup combined with time-lapse microscopy allowed us to follow the development of cell communities over time. We found that over 24 hours dense groups of more than 1000 cells formed on alginate, filling the entire microfluidic chamber (Fig. 1B). In contrast, cells supplied with digested alginate grew in smaller groups that never exceeded 100 cells per chamber (Fig. 1C and D). Reconstructed cell lineages revealed that the large cell groups on alginate formed because cells often did not disperse after division and thereby formed dense cell groups originating from a single cell lineage (Fig. 1B). The lower cell density on digested alginate could be caused by slower growth or by more cells leaving the chambers. As the maximum growth rate both in bulk and on single cell level is similar in alginate and digested alginate (Fig. S2E and F), it is more likely that the lower cell density in digested alginate is caused by cell dispersal. To test the role of increased viscosity of polymeric alginate in causing the increased aggregation of cells, we measured the viscosity of 0.1% (w/v) alginate or digested alginate dissolved in TR media. For alginate, the viscosity was 1.03±0.01 mPa·s (mean and standard deviation of three technical replicates) whereas the viscosity of digested alginate in TR media was found to be 0.74±0.01 mPa·s. Both these values are relatively close to the viscosity of water at this temperature (0.89 mPa·s18) and, while they may affect swimming behavior19, they are insufficient to physically restrain cell movement20. Overall, our observations suggest that cells can modulate their propensity to form groups depending on the state of polysaccharide degradation in their local environment. Similar observations were made with a different model system (Caulobacter crescentus growing on the polysaccharide xylan)16, indicating that group formation on polymeric nutrient sources may be a general mechanism of bacteria that degrade polysaccharides extracellularly.

Large cell groups form on alginate but not on digested alginate.

(A) Schematic representation of the setup of the microfluidic experiments. (B and C) Representative images at different time points of V. cyclitrophicus ZF270 cells growing in microfluidic chambers, described in detail by Dal Co et al.15, with (B) alginate medium or (C) digested alginate medium, both in their soluble form (not visible). Cells are false-colored according to their lineage identities based on cell segmentation and tracking over 24 hours. Cells without identified progenitors are colored in dark blue. See Supplemental Videos S1 (alginate) and S2 (digested alginate) for time-lapse videos. (D) Cell numbers within microfluidic chambers supplied with alginate (orange) are substantially higher than cell numbers within microfluidic chambers supplied with digested alginate (blue) (Logistic growth regression for alginate: R2 = 0.99, maximal number of cells = 1217-1564, k = 0.24-0.38 h-1; for digested alginate: R2 = 0.86-0.97, maximal number of cells = 47-100, k = 0.07-0.4 h-1). Circles indicate the number of cells present at a given time point in each chamber (nchambers = 7). Data for chambers with alginate originate from D’Souza et al.8. Lines are fits of a logistic growth regression line for each condition.

Transition from alginate to digested alginate triggers density-dependent dispersal of cells

To investigate the transition from the large cell groups formed on alginate to the small groups formed on digested alginate, we subjected V. cyclitrophicus ZF270 cells grown on alginate to a switch to digested alginate (Fig. 2A). Following the limited cell motility on alginate, this switch led to a rapid decrease in cell density within the growth chambers (Fig. 2B), presumably caused by cell dispersal. As we previously reported for other Vibrionaceae isolates, the growth rate of the cells on alginate was dependent on the local cell density: Initially the growth rate increased with cell density but then decreased at high cell densities, indicating that cell groups can benefit from the sharing of breakdown products generated by each other’s exoenzymes, but also increasingly compete for nutrients8. This led us to investigate whether cells in larger groups, potentially experiencing stronger nutrient competition, might have a higher propensity to disperse after a switch to digested alginate than cells in smaller groups. We indeed found that the nutrient switch caused a few or no cells to disperse from small cell groups (Fig. 2B), whereas a large fraction of cells from large cell groups dispersed (Fig. 2C). In fact, the e fraction of cells that dispersed upon imposition of the nutrient switch showed a strong positive relationship with the number of cells present, meaning that cells in chambers with many cells were more likely to disperse than cells in chambers with fewer cells (Fig. 2C). Thus, during the transition from polysaccharides to degradation products, we found that the dispersal rate of cells depends on the size of the cell groups, likely through increased motility of cells in large groups.

Transition from alginate to digested alginate triggers density-dependent dispersal of cells.

(A) Representative time-lapse images of V. cyclitrophicus ZF270 cells (phase contrast microscopy) in microfluidic growth chambers that were initially exposed to alginate and then switched to digested alginate. (B) Number of cells in different chambers over time, each chamber indicated by a unique color (n = 8). The carbon source is indicated by the colored background (orange: alginate; blue: digested alginate). See Supplemental Video 3 for a time-lapse video. (C) Positive relationship between the number of cells in the microfluidic growth chamber at the time of the switch and the fraction of cells that disperse after the switch. Each circle represents one growth chamber with colors corresponding to (B), and the line depicts a linear regression fit (R2 = 0.92, slope = 0.12).

Cells growing on digested alginate are more motile and polymeric alginate acts as chemoattractant

To investigate whether increased cell motility of V. cyclitrophicus ZF270 underlies the dispersal of cells observed after a shift to digested alginate, as also previously observed for C. crescentus cells on the monosaccharide xylose4,16, we quantified the motility of single cells supplied with either alginate or digested alginate. We found that the single-cell swimming speed as well as the swimming distance were significantly larger for cells supplied with digested alginate compared to cells supplied with alginate, and that a larger fraction of cells was motile (Fig. 3). This confirmed increased motility as a cellular response to the exposure to degradation products in the form of digested alginate. Overall, these findings suggest that once the breakdown of a polysaccharide source makes breakdown products available in the local environment, a fraction of cells becomes motile and disperses.

To understand whether alginate polymers or alginate breakdown products act as chemoattractants on motile cells, we measured the chemotactic strength towards alginate and digested alginate using the In Situ Chemotaxis Assay (ISCA)21,22. Interestingly, we found V. cyclitrophicus ZF270 to significantly chemotax toward alginate (chemotactic index Ic > 1) but not significantly toward digested alginate (Fig. 3F). This suggested that the increase in motility is accompanied by chemotaxis toward alginate polymers.

Cells are more motile on digested alginate than alginate and show chemotaxis towards alginate.

(A and B) Spatial trajectories of cells supplied with (A) alginate or (B) digested alginate in representative microfluidic growth chambers are shown. Black points mark the starting point of each trajectory, pink points mark the end point of each trajectory, and colored lines mark the trajectories of cells. (C) Distributions of the mean single-cell swimming speeds (Nested t-test, p-value < 0.0007, t = 4.803, df = 10, ncells = 86 vs 375 in nchambers = 5) are shown. (D) Distributions of cell displacement over the course of a trajectory (Nested t-test test, p-value < 0.0131, t = 4.39, df = 10, ncells = 86 vs 375, and nchambers = 5) are shown. In (C) and (D) the red horizontal lines indicate the mean while black lines depict the 25th and 75th quartiles of the distribution. (E) The mean fraction of motile cells in each chamber, where motile cells are defined as cells with displacement greater than 1 µm (Mann-Whitney test on the means of five growth chambers, p-value = 0.034). In C, D, and E, each chamber was considered as an independent replicate. (F) Chemotactic index (IC) quantified by In Situ Chemotaxis Assay (ISCA) (Tukey multiple comparisons of means, 95% family-wise confidence levels as error bars, p-value < 0.05, n = 3). Asterisks indicate statistically significant differences. See Supplemental Videos S4 and S5 for time-lapse videos of swimming cells.

Altered gene expression in central carbon metabolism, enzyme production, secretion and transporters, motility, and quorum sensing underlies the late-stage alginate degradation and cell dispersal

Next, we sought to elucidate the molecular mechanisms underlying the observed phenomenological disparities between cells cultivated on alginate and digested alginate. Due to the challenge of generating knock-out mutants in natural isolates, we used transcriptomics to investigate the differentially expressed genes of V. cyclitrophicus ZF270 under these respective conditions. To obtain a high-quality reference genome of V. cyclitrophicus ZF270, we sequenced and assembled a new reference genome using combined short and long read sequencing (BioProject PRJNA991487). We then grew cultures on either alginate or digested alginate until mid-exponential phase. To understand which cellular functions were affected by the expression changes, we first performed differential gene expression analysis using the software DESeq2. Here, genes exhibiting a log2 fold expression change greater than 0.5 or smaller than -0.5 between the two conditions, with a Benjamini-Hochberg(BH)-adjusted p-value below 0.01, were considered to be differentially expressed (Table S1). Next, we investigated which KEGG categories were enriched in either genes with increased or decreased expression via Gene Set Enrichment Analysis (GSEA)23 (Table S2) and found nine categories significantly enriched in genes with increased gene expression on digested alginate, and three categories significantly enriched in genes with decreased gene expression (Fig. 4A).

Cells growing on digested alginate showed expression changes in parts of the central metabolism and translation, compared to cells growing on alginate. Specifically, two pathways of the central metabolism were significantly enriched in genes with lower expression on digested alginate (”Valine, leucine and isoleucine biosynthesis” and “Propanoate metabolism” with a negative normalized enrichment score (NES), Fig. 4A, S4A and B, Table S3 and S4). We also found the gene set that encodes ribosomal proteins significantly enriched with genes highly expressed on digested alginate (positive NES, Fig. 4A and S4C, Table S5). This implies that cells invest proportionally more of their transcriptome into the production of new proteins, a sign of faster growth24,25. Both observations likely relate to the different growth dynamics of V. cyclitrophicus ZF270 on digested alginate compared to alginate (Fig. S2A), where cells in digested alginate medium reached their maximal growth rate 7.5 hours earlier and thus showed a shorter lag time (Fig. S2D). As a consequence,, the growth rate at the time of RNA extraction (mid-to-late exponential phase) may have differed, even though the maximum growth rate of cells grown in alginate medium and digested alginate medium were not found to be significantly different (Fig. S2E).

The expression of transporters and secretion systems was generally increased in cells growing on digested alginate, compared to cells growing on alginate. This includes genes with a function in “Bacterial secretion systems”, namely secretion systems I to VI, which mediate protein export through the inner and outer membranes of Gram-negative bacteria (Fig. 4A and S4D, Table S6). Notably, eleven of the thirteen General Secretion Pathway (GSP) genes (part of the type II secretion system) showed 1.1 to 3.4-fold increased expression levels on digested alginate (p-values < 0.01, Fig. 4A, Table S8). The GSP is known to facilitate secretion of various extracellular enzymes like chitinases, proteases, and lipases26,27, therefore, the positive enrichment on digested alginate may be linked to the export of extracellular alginate lyases. The KEGG BRITE category “Secretion system” was also enriched on digested alginate (Fig. 4A and S4D, Table S7). It contains additional genes involved in protein export, implicating them as interesting subjects to further research in their role in protein secretion and cell attachment and detachment during degradation-dispersal cycles of polysaccharide-degrading bacteria. Also, the gene set of transporters and in particular ABC transporters were positively enriched (Fig. 4A, Table S10 and S11). The latter showed enrichment especially in genes related to saccharide, iron, zinc, and phosphate transport (Fig. S4E), which suggests that cells growing on digested alginate invest proportionally more of their transcriptome into uptake of not only saccharides but also essential nutrients like iron, zinc, and phosphate, which may become growth-limiting when degradation products are abundantly available.

Both motility and quorum sensing genes increased in expression in cells growing on digested alginate. The positive enrichment of the gene set containing bacterial motility proteins aligned with our expectations based on the increase in motile cells that we observed in Figure 3E (Fig. 4A, Table S12). The set of quorum sensing genes was also positively enriched in cells growing on digested alginate (Fig. 4A and S4F, Table S13). This role in dispersal is in agreement with a previous study that showed induction of the quorum sensing master regulator in V. cholerae cells during dispersal from biofilms on a similar time scale as here (less than an hour)28. Quorum sensing is known to control biofilm formation in the well studied model system Vibrio cholerae29,30 and may also orchestrate the density-dependent dispersal in the presence of degradation products observed in our study, though the particular signaling cues remain to be uncovered. The strong cellular response to the degree of alginate depolymerization was also emphasized by the finding that transcription factor genes were enriched among genes with decreased expression on digested alginate (Fig. 4, Table S14). The set of genes associated to KEGG’s beta-lactam resistance category was enriched in genes with increased expression on digested alginate. However, most genes in this gene set are also associated to the KEGG pathways “Transporters”, “Quorum sensing”, “Chromosome and associated genes”, and “Peptidoglycan biosynthesis” and thus the enrichment likely reflects expression changes in these categories and may relate to the difference in growth dynamics (Fig. 4A and S4G, Table S15). Overall, the observed gene expression changes provide insights into the molecular mechanisms underlying the substantial adaptations of V. cyclitrophicus ZF270 cells to the polymeric and digested form of alginate.

Twelve functional gene sets are enriched in genes with increased or decreased expression in cells grown on digested alginate.

Gene set enrichment analysis (GSEA) with (A) all KEGG pathways and KEGG BRITE categories as gene sets or with (B) a custom alginate utilization, flagellar assembly, and flagellum-driven chemotaxis gene set was performed comparing the gene counts of the transcriptome of V. cyclitrophicus ZF270 cultures grown on digested alginate and alginate (six replicates each). Gene sets with a positive enrichment score were enriched with genes with higher expression in cells grown on digested alginate relative to cells grown on alginate (BH-adjusted p-value < 0.05), whereas gene sets with negative enrichment scores were significantly enriched with genes with decreased expression on digested alginate. The number in brackets indicates the number of genes with unique K number per gene set (A) and the number of genes per gene set (B) within the V. cyclitrophicus ZF270 genome.

Growth on digested alginate is associated with increased expression of alginate catabolism, flagellar motility and chemotaxis genes

For a more fine-grained understanding of the impact of the form of alginate on the alginate catabolism, we investigated specifically the expression of genes involved in alginate degradation, uptake, and catabolism. In the closed genome of V. cyclitrophicus ZF270, we identified genes that encode CAZymes responsible for alginate degradation, namely alginate lyase genes from the PL6, PL7, PL15, and PL17 family. We identified homologs of alginate transporters (porin kdgM, symporter toaA, toaB, and toaC) and metabolic enzymes that shunt into the Entner-Doudoroff pathway (DEHU reductase genes dehR, kdgK, and eda), based on the known genes in the alginate degradation pathway of Vibrionaceae31,32. We found that the expression of most of these genes increased significantly on digested alginate relative to alginate (Fig. 4B, Fig. 5A, S8A, Table S16). Surprisingly, also the expression of most alginate lyase genes increased on digested alginate, especially the expression of secreted alginate lyase genes. This indicates the production of “public” exoenzymes despite the abundance of monomeric and oligomeric degradation products in the digested alginate medium.

The increased motility of cells observed upon exposure to digested alginate (Fig. 3, Fig. 4A) led us to evaluate the expression of motility- and chemotaxis-associated genes across digested alginate and alginate treatments. As flagella are the main mode of motility in the genus Vibrio33,34, we focused on the expression of flagella-related genes. In the genome of V. cyclitrophicus ZF270 we found a gene cluster that encodes most genes involved in flagellar assembly and that was flanked by chemotaxis genes (hereon called flagellar locus, Fig. S8B, and Table S17). Overall, 21/34 flagellar locus genes were differentially expressed (log2 fold change > 0.5, BH-adjusted p-value < 0.01) and the majority (90%) of these differentially expressed genes showed increased expression on digested alginate (Fig. 4B). The flagellar genes flgA, fliC, and fliH and the chemotaxis gene cheW showed the strongest overexpression (5.4, 3.8, 2.0, and 2.6-fold, respectively). The expression of flagellar biosynthesis genes in Vibrionacaea occurs by a cascade of gene expression of four classes of genes (Class I - IV)35. We found the expression of the master regulator of the flagellar biosynthesis regulon, the Class I gene flrA, to be 1.6-fold increased on digested alginate (BH-adj. p-value 3e-19)36. The Class II regulatory genes flrBC, controlling Class III genes, and fliA, controlling Class IV genes, showed 1.5-fold and 1.3-fold increased expression on digested alginate (BH-adj. p-value 1e-12 for flrB, 2e-39 for flrC, 6e-15 for fliA) (Fig. 5B)34,37. These findings suggest that the increased phenotypic motility observed on digested alginate (Fig. 3) is related to the upregulation of flagellar biosynthesis genes. Additionally, the expression of the flagellum filament, encoded by fliC genes of the flagellar locus, was partially increased: flaD and flaC expression were increased by 3.8 and 1.6-fold, whereas the flaA gene was not significantly differentially expressed (BH-adj. p-value 7e-37, 2e-19, and 1e-02, respectively) (Fig. 5B, Table S6). This suggests that cells grown on digested alginate have a flaD-rich filament composition, which has been shown to alter the swimming and adhesion characteristics of bacterial cells38,39. Lastly, we found that most genes involved in chemotaxis and located in the flagellar locus are highly expressed in cells grown on digested alginate (Fig. 5C, Fig. 4B) and likely drive the chemotactic activity of V. cyclitrophicus ZF270 during dispersal. Overall, our findings elucidate that cellular responses upon exposure to degradation products manifest in increased expression of genes involved in extracellular alginate breakdown and alginate catabolism as well as flagellar assembly and chemotaxis.

Digested alginate increases expression of genes involved in alginate degradation, uptake and catabolism, as well as flagellar assembly and chemotaxis.

Genome-wide differential expression analysis where the log2 fold changes of gene expression on digested alginate compared to alginate is shown for (A) alginate lyases (PL6, PL7, PL15, PL17, scissors symbol), transporters (porin kdgM, symporter toaB, symporter toaC), and metabolic enzymes shunting into the Entner-Doudoroff pathway (DEHU reductase DehR, kdgK, eda), (B) genes of the flagellar locus associated with flagellar assembly and (C) adjacent chemotaxis genes. Genes displayed in (B) and (C) are part of the KEGG pathways “Bacterial motility proteins” and “Bacterial chemotaxis”. Differential expression analysis was performed to compute the Benjamini-Hochberg-adjusted Wald test p-value (”BH-adj. p-value”, text color and box outline color) and log2 fold change (box fill color) for each gene (box). For better visibility, genes that exhibited a log2 fold gene expression change greater than 1 (i.e., doubling of expression) or less than -1 (i.e., halving of expression) are designated maximum intensity of red or blue, respectively. Genes with BH-adj. p-value smaller than 0.01 were considered significantly differentially expressed. In (A), the location of the gene products was based on Figure 1 of Wargacki et al.31 with the exception of the alginate lyases (PL6, PL7, PL15, PL17) which were placed based on their signal peptides (S: extracellular, LS: membrane-embedded, none: cytosolic). In (B) and (C) the gene location and depiction were based on the KEGG pathway “Flagellar assembly” (map02040), “Bacterial chemotaxis” (map02030), and Figure 3 of Rajagopala et al.40. Genes without known cellular location were omitted here but displayed in the genomic architecture in Fig. S5. Arrow: activation; dashed arrow: modification; “flat” arrow: inhibition; OM: outer membrane; PM: periplasm; IM: inner membrane; PL: polysaccharide lyase family; kdgM: oligogalacturonate-specific outer membrane porin; toaABC: oligoalginate symporter; DEH: 4-deoxy-L-erythro-5-hexoseulose uronic acid; dehR: DEH reductase; KDG: 2-keto-3-deoxy-gluconate; kdgK: KDG kinase; KDPG: 2-keto-3-deoxy-6-phosphogluconate; eda: KDG-6-phosphate aldolase; GAP: glyceraldehyde 3-phosphate; ED: Entner-Doudoroff; ns: not significant, i.e. BH-adj. p-value > 0.01.

Discussion

On Earth organic carbon is mostly present in the form of polysaccharides1,2, which are often in a particulate state and form a heterogeneous resource landscape. Over the last years, the study of extracellular bacterial degradation of polysaccharides has revealed that bacterial growth on polysaccharides increases with increased cell density, enabling cells to benefit from the exoenzymes and extracellular degradation products of surrounding cells otherwise lost to diffusion (”cooperative growth”)4143. However, cells have been observed not only to aggregate on polysaccharide sources, but also to leave them before the source is depleted13,43. This cycle of biomass degradation and dispersal has long been discussed in the context of foraging e.g., 44,45,13,46,47, but the cellular mechanisms that drive the cell dispersal remain unclear.

Our work links these observations and connects them to the cellular mechanisms that underlie the degradation-dispersal cycles of bacterial degraders, which we see as basal drivers of the biogeochemical processing of polysaccharides in heterogeneous nutrient-scapes. When bacteria encounter a new source of biomass, their local environment likely contains few mono- or oligosaccharides but is rich in polysaccharides which usually require extracellular breakdown (Fig. 6, “Finding a new nutrient source”). General concepts of how bacteria recognize the presence of large biopolymers remain elusive, but it was proposed that “sentry” enzymes are constitutively expressed at a basal level to cleave mono- or oligosaccharides from polysaccharides, which cells can take up and which prime their metabolism for the degradation of the respective polysaccharide48,49. It is not known yet how widespread the concept of sentry enzymes may be, but the observation of a constitutively expressed PL7 and PL15 family alginate lyase gene in Z. galactanivorans48 is mirrored in our work by the constant expression of an extracellular PL7 family alginate lyase gene, which may act as sentry enzyme that helps to initiate alginate degradation when cells encounter alginate.

Growth on polysaccharides has been found to be dependent on the cell density, as increased cell density limits the loss of degradation products and exoenzymes by diffusion8 (Fig. 6, “Early-stage polysaccharide breakdown”). We found that cells grown on degradation products reach their maximal growth earlier and show increased expression of ribosomal biosynthesis, enzyme secretion, especially of secreted alginate lyases, transporters, quorum sensing and expression changes in the central carbon metabolism. The secretion of alginate lyases might seem surprising and wasteful in a monomer-rich environment. One reason for this observation may be that cells primarily rely on intracellular monosaccharide levels to trigger the upregulation of genes associated with polysaccharide degradation and catabolism, as has previously been observed for E. coli across various carbon sources50,51. In fact, the majority of carbon sources are sensed by prokaryotes through one-component sensors inside the cell50. In the one-component internal sensing scheme, the enzymes and transporters for the use of various carbon sources are expressed at basal levels, which leads to an increase in pathway intermediates upon nutrient availability. The pathway intermediates are sensed by an internal sensor, usually a transcription factor, and lead to the upregulation of transporter and enzyme expression50,51. This results in a positive feedback loop, which enables small changes in substrate abundance to trigger large transcriptional responses50,52. Thus, the presence of alginate breakdown products may likely result in increased expression of all components of the alginate degradation pathway, including the expression of degrading enzymes. As the gene expression analysis was performed on well-mixed cultures in culture medium containing alginate breakdown products, we therefore expect a strong stimulation of alginate catabolism. In a natural scenario, where cells disperse from a polysaccharide hotspot before its exhaustion, the expression of alginate catabolism genes may likely decrease again once the local concentration of breakdown products decreases. However, continued production of alginate lyases could also provide an advantage when encountering a new alginate source and continued production of alginate lyases may thus help cells to prepare for likely future environments. Further investigations of bacterial enzyme secretion in changing nutrient environments and at relevant spatial scales are required to improve our understanding of the regulation of enzyme secretion along nutrient gradients.

We show that cells respond to the exposure to degradation products with dispersal from dense cell groups by means of increased flagellum-driven swimming (Fig. 6, “Late-stage polysaccharide breakdown”), decreasing their local cell number. This finding matches with previous observations of cells leaving biopolymer particles before they are depleted13,43. A plausible explanation for this density-dependent dispersal is that cells in larger groups compete with each other for nutrients and space, while not profiting from cooperative degradation anymore due to the abundance of degradation products. Motility has also been shown to increase the encounter rate of cells with sources of nutrients53,54, suggesting motility as a strategy that allows cells to escape from the ensuing competition. Previous work in Caulobacter crescentus demonstrated that a flagellum knock-out mutant formed larger cell groups, resulting in reduced growth rates due to intercellular competition4,16. While it would be interesting to study non-motile mutants of V. cyclitrophicus ZF270, the non-tractability of natural isolates makes direct tests of molecular mechanisms difficult. Additionally, subjecting cells separately to the two monomeric units of alginate or oligomers of defined size could improve our understanding of the specific molecules that trigger motility, but this was experimentally not feasible.

Direct chemotaxis towards polysaccharides may facilitate the search for new polysaccharide sources after dispersal. We found that the presence of degradation products not only induces cell dispersal but also increases the expression of chemotaxis genes. Interestingly, we found that V. cyclitrophicus ZF270 cells show chemotaxis towards polymeric alginate but not digested alginate. This contrasts with previous findings for bacterial strains degrading the insoluble marine polysaccharide chitin, where chemotaxis was strongest towards chitin oligomers53, suggesting that oligomers may act as an environmental cue for polysaccharide nutrient sources55. However, recent work has shown that certain marine bacteria are attracted to the marine polysaccharide laminarin, and not laminarin oligomers56. Together with our results, this indicates that chemotaxis towards soluble polysaccharides may be mediated by the polysaccharide molecules themselves. The mechanism of this behavior is yet to be identified, but could be mediated by polysaccharide-binding proteins as have been found in Sphingomonas sp. A1 facilitating chemotaxis towards pectin57. Direct polysaccharide sensing adds complexity to chemosensing as polysaccharides cannot freely diffuse into the periplasm, which can lead to a trade-off between chemosensing and uptake58. Furthermore, most polysaccharides are not immediately metabolically accessible as they require degradation. But direct polysaccharide sensing can also provide certain benefits compared to using oligomers as sensory cues. First, it could enable bacterial strains to preferably navigate to polysaccharide nutrients sources that are relatively uncolonized and hence show little degradation activity. Second, strong chemotaxis towards degradation products could hinder a timely dispersal process as the dispersal then requires cells to travel against a strong attractant gradient formed by the degradation products. Overall, this strategy allows cells to alternate between degradation and dispersal to acquire carbon and energy in a heterogeneous world with nutrient hotspots44,5961.

Bacterial growth and regulation on patches of polysaccharides.

By integrating our results with previous studies on cooperative growth on the same system, as well as results on dispersal cycles in other systems, we highlight where the specific results of this work add to this framework (bold font). When cells encounter polymer sources, the colonization of the nutrient hotspot may be aided by the basal exoenzyme production of “sentry” enzymes (”Encounter of a new carbon source”). This phase is succeeded by group formation, which enables cells to benefit from exoenzymes of neighboring cells and diffusing degradation products (”Early-stage polysaccharide breakdown”). The following phase includes cooperative extracellular degradation of the polysaccharide source, further increasing the concentration of available degradation products. These degradation products trigger the overexpression of alginate degrading, importing, and catabolizing enzymes, ensuring swift polysaccharide degradation (”Late-stage polysaccharide breakdown”). The increased pool of breakdown products also cues flagellar swimming in a subpopulation of cells and increases the expression of chemotaxis genes. Polymeric alginate acts as chemoattractant towards new polysaccharide sources. Cells and molecules are not drawn to scale. Dark red pie symbols: intracellular and extracellular polysaccharide-degrading enzymes; orange shading: a polymeric carbon source; blue shading: monomeric or oligomeric degradation products.

Conclusion

The heterogeneous landscape of polysaccharide hotspots in natural systems requires bacteria to effectively break down polymeric carbohydrates as well as readily ensure dispersal to new nutrient hotspots. Our findings show that the degree of depolymerization of the polysaccharide influences this decision, altering the growth dynamics, metabolic activity, and motility of cells. Our study also contextualizes the surprising finding that foraging bacteria majorly leave polysaccharide particles before the last third of the particle is consumed43. Dispersal from a partially degraded carbon source may serve several purposes: (i) escaping competition that ensues within large cell groups, (ii) ensuring the spread of a part of the clonal population to new environments as bet hedging strategy47,62, here guided by chemotaxis towards new nutrient hotspots, (iii) preventing whole populations degrading a sinking marine particle or a deposited sediment particle to be buried in depth where nutrient hotspots become sparse43,63, and/or (iv) increase the genetic variation in bacterial populations47. However, dispersal may also occur when a nutrient source offers a surplus of carbon while other essential nutrients become limiting, as the increased expression of iron, zinc, and phosphate transporters in cells grown on digested alginate suggested. These findings emphasize that metabolic molecules can also act as triggers of dispersal, expanding upon the current perspective of dispersal in biofilms as a reaction to dispersal cues like NO, signaling molecules, nutrient starvation, and oxygen starvation64. The study of bacterial motility on increasingly complex biomass particles will reveal the role of the nutrient composition of the present nutrient hotspot on the bacterial decision-making. Overall, these new insights into the cellular mechanisms and regulation that drive degradation-dispersal cycles contribute to our understanding of the microbially driven remineralization of biomass, and factors that modulate this process. The open questions of how bacteria sense polysaccharides in their environment, which cell signaling pathways integrate the presence of degradation products in the cellular decision-making of degradation and dispersal, and to what extent cell populations coordinate this decision, present an exciting avenue of further research.

Materials and Methods

Bacterial strains, media and growth assays

Vibrio cyclitrophicus ZF270 (available through Culture Collection Of Switzerland; Accession number: 2043) cells were cultured in Marine Broth (DIFCO) and grown for 18 hours at 25 °C. Cells from these cultures were used for growth experiments in Tibbles Rawling (TR) salts minimal medium65,66 containing either 0.1% (weight/volume) algae-derived alginate (referred to as “alginate”) (Sigma Aldrich, CAS-number 9005-38-3) or 0.1% (weight/volume) digested alginate. At these concentrations, both alginate and digested alginate are soluble in the culture medium. The digested alginate was produced by enzymatically digesting 2% alginate with 1 unit ml-1 of alginate lyase (Sigma Aldrich, CAS-number 9024-15-1) at 37 °C for 48 hours. Carbon sources were prepared in nanopure water and filter sterilized using 0.40 µm Surfactant-Free Cellulose Acetate filters (Corning, USA). Well-mixed batch experiments in alginate or digested alginate medium were performed in 96-well plates (Greiner Bio) and growth dynamics were measured using a microwell plate reader (Biotek, USA). Plate reader assays were initiated as described previously67. Briefly, 1 ml from a culture grown for 18 h on Marine broth was centrifuged at 5000 g in 1.5 ml microfuge tubes for 5 minutes. The supernatant was discarded and the cell pellet was subjected to two rounds of washing with the basal TR salts medium. The cell-pellet was then resuspended in 1 ml of TR salts medium and 5 µl of this suspension inoculated into 195 µl TR medium with either carbon source (∼105 colony forming units (CFUs ml−1) in a 96-well plate (Greiner Bio). The optical density (600 nm) was then measured every 15 minutes for 40 hours. All measurements had six biological replicates.

Alginate oligosaccharide measurements

Oligosaccharide measurements were performed using liquid chromatography time of flight mass spectrometry (LC-QTOF-MS). Samples were prepared by diluting 1:20 in milliQ water and 5 µL of sample was injected per measurement. Chromatographic separation was performed using an Agilent 1290 stack, using an Agilent HILIC-Z column (2.7 µm particles, 2.1 x 50 mm). Mobile phase A contained 10% acetonitrile (Fisher Scientific) and 0.1% medronic acid (Agilent), and Mobile phase B contained 90% acetonitrile and 0.1% medronic acid. The separation was performed as follows: Mobile phase B 100% for 1 minute, gradient to 30% phase B over 3 minutes, 30% phase B for 30 seconds, and equilibration of 100% phase B for 5 minutes. The flow rate was 400 µL min−1 at 30 °C. Samples were measured using an Agilent 6520 mass spectrometer in negative mode, in 4 GHz high-resolution mode. Data analysis was performed in Agilent Quantitative Analysis software.

Microfluidics and time-lapse microscopy

Microfluidic experiments and microscopy were performed as described previously15,16,68. Cells were imaged within chambers of a PDMS (Sylgard-Dow) microfluidic chip that ranged in size from 60-120 × 60 × 0.56 μm (l × b × h). Within these chambers, cells can attach to the glass surface and experience the medium that diffuses through lateral flow channels. Imaging was performed using IX83 inverted microscope systems (Olympus, Japan) with automated stage controller (Marzhauser Wetzlar, Germany), shutter, and laser-based autofocus system (Olympus ZDC 2). Chambers were imaged in parallel on the same PDMS chip, and phase-contrast images of each position were taken every 8 or 10 minutes. The microscopy unit and PDMS chip were maintained at 25 °C using a cellVivo microscope incubation system (Pecon GmbH).

Viscosity of the alginate and digested alginate solution

We measured the viscosity of alginate solutions using shear rheology measurements. We use a 40 mm cone-plate geometry (4° cone) in a Netzsch Kinexus Pro+ rheometer. 1200 uL of sample was placed on the bottom plate, the gap was set at 150 um and the sample trimmed. We used a solvent trap to avoid sample evaporation during measurement. The temperature was set to 25°C using a Peltier element. We measure the dynamic viscosity over a range of shear rates = 0.1 – 100 s-1. We report the viscosity of each solution as the average viscosity measured over the shear rates 10 – 100 s-1, where the shear-dependence of the viscosity was low.

We measured the viscosity of 0.1% (w/V) alginate dissolved in TR media, which was 1.03 +/-0.01 mPa·s (reporting the mean and standard deviation of three technical replicates.). The viscosity of 0.1% digested alginate in TR media was found to be 0.74+/-0.01 mPa·s. This means that the viscosity of alginate in our microfluidic experiments is 36% higher than of digested alginate, but the viscosities are close to those expected of water (0.89 mPa·s at 25 degree celsius according to Berstad and colleagues (1988)18).

Motility assays

Cells were grown for 10 hours in Marine Broth (DIFCO) after which 10 µl of culture was used to inoculate culture tubes (Greiner) containing 5 ml of TR medium with either 0.1% alginate or 0.1% digested alginate. After 6 hours of growth at 25 °C, 2 µl of cell suspension was inoculated into microfluidic growth chambers. Cells within six replicate chambers were then imaged with the phase-contrast channel at a high frame rate (125 Hz, i.e., frames s-1) using the same microscopy setup described above.

Chemotaxis assays

To assess whether polymeric alginate and digested alginate attract Vibrio cyclitrophicus ZF270, we used the In Situ Chemotaxis Assay21,22 (ISCA), a microfluidic device consisting of a 5 × 5 array of microwells that can be individually loaded with solutions of different chemicals (110 µl each). Once the ISCA is deployed in an aqueous environment, the chemicals diffuse out of the wells through a small port, creating chemical gradients which will guide chemotactic bacteria inside the wells of the device21,22. Vibrio cyclitrophicus ZF210 was plated on Marine Agar (BD Difco) from a glycerol stock and grown for 16 h at 27 °C. A single colony was then incubated in 10% Marine Broth (BD Difco) in 0.22-μm filtered artificial seawater (Instant Ocean, Spectrum Brands) and grown overnight at 27 °C and 180 rpm. The culture was diluted down to 1 x 106 cells ml-1 in 0.22 μm filtered artificial seawater (Instant Ocean, Spectrum Brands) to perform the chemotaxis experiment. Both chemoattractants (alginate and digested alginate) were diluted in sterile seawater (35 g l-1; Instant Ocean, Spectrum Brands) at a final concentration of 0.1% and then filtered with a 0.2 μm filter (Millipore) to remove particles and potential contaminants. Within the ISCA, one full row of five wells was used per chemoattractant as technical replicates. The chemoattractants were injected in triplicate ISCA with a sterile 1 ml syringe (Codau) and needle (27 G, Henke Sass Wolf). A last row containing 0.2 μm-filtered seawater acted as negative control accounting for cells swimming in the device by random motility only. Experiments were conducted by incubating the ISCAs for 1 h in the diluted Vibrio cyclitrophicus ZF270 culture. Upon time completion, a sterile syringe and needle were used to retrieve the content of the wells and transferred to 1 ml microfuge tubes resulting in a pooling of a row of 5 wells containing the same sample. Sample staining was performed with SYBR Green I (Thermo Fisher) and the chemotactic response was quantified by counting cells using flow cytometry. The strength of the chemotactic response was determined by the mean chemotactic index (IC), defined as the ratio of the number of cells found in each chemoattractant to the number of cells in control wells containing filtered seawater (so that attraction corresponds to IC > 1).

Culturing and harvesting cells for transcriptomics

Cells were grown for 18 hours in Marine Broth (DIFCO) after which 1 ml of culture was centrifuged at 5000 g in 1.5 ml microfuge tubes for 5 minutes. The supernatant was discarded and the cell pellet was subjected to two rounds of washing with the basal TR salts medium. The cell pellet was then resuspended in 1 ml of TR salts medium and 250 µl of this suspension were used to inoculate 100 ml flasks (Schott-Duran) containing 10 ml of TR medium with either 0.1% alginate or 0.1% digested alginate. This was done in parallel for six flasks. Once cultures in the flasks reached mid-exponential phase (10 hours and 15 hours after inoculation for digested alginate and alginate, respectively) and had approximately the same OD (0.39 for digested alginate and 0.41 for alginate), 2 ml of cultures were harvested for RNA extraction. Samples were stabilized with the RNprotect reagent (Qiagen) and RNA was extracted using the RNeasy mini kit (Qiagen).

Sequencing and gene annotation of Vibrio cyclitrophicus ZF270

Long read sequencing using the Oxford Nanopore Platform (Long read DNA sequencing kit) and short read sequencing using the Illumina platform (Illumina DNA Prep kit and IDT 10bp UDI indices, and sequenced on an Illumina NextSeq 2000 producing 2x151bp reads) was performed by the Microbial Genome Sequencing Center, Pittsburgh, USA (MiGS), to create a new closed reference genome of V. cyclitrophicus ZF270 (BioProject PRJNA991487). Annotation of this genome was done with RASTtk (v2.0, Rapid Annotation using Subsystem Technology tool kit69). Additionally, KEGG Ontology identifiers (”K numbers”) were annotated with BlastKOALA (v2.2)70. Dedicated annotation of alginate lyase genes was performed by homology search for proteins belonging to the PL5, PL6, PL7, PL14, PL15, PL17, PL18, PL31, PL36, or PL39 family71 by dbCAN2 (v9.0)72. Enzymes for alginate transport and metabolism were identified by BLASTn-search73,74 of gene sequences of Vibrio splendidus 12B01, which were previously identified as minimum genetic prerequisites for alginate utilization and enabled alginate degradation when cloned into E. coli31.

Location prediction of alginate lyases

Signal peptides were annotated using SignalP (v.5.0)75, and LipoP (v.1.0)76. SignalP discriminated between 1) Sec/SPI: “standard” secretory signal peptides transported by the Sec translocon and cleaved by Signal Peptidase I (SPI), 2) Sec/SPII: lipoprotein signal peptides transported by the Sec translocon and cleaved by Signal Peptidase II (SPII), and 3) Tat/SPI: Tat signal peptides transported by the Tat translocon and cleaved by SPI. LipoP discriminates between 1) SPI: signal peptide, 2) SpII: lipoprotein signal peptide, and 3) TMH: n-terminal transmembrane helix. All predictions were in agreement, apart from one PL7 (gene 1136176.5.peg.4375) which was predicted by LipoP as cytoplasmic and by SignalP as equally likely cytoplasmic as containing a lipoprotein signal peptide.

Transcriptomic analysis: Sequencing, pre-processing, differential expression analysis, and functional analysis

Sequencing (12 M reads, 2 x 50 bp) of the isolated RNA was performed by MiGS after rRNA depletion using RiboZero Plus (Ilumina). cDNA libraries were prepared using an Illumina DNA Prep kit and IDT 10 bp UDI indices, and sequenced on an Illumina NextSeq 2000. Preprocessing of the raw reads was carried out as follows: Quality control was performed with FastQC (v0.11.9)77 and reads were trimmed with Trimmomatic (v0.38)78; the high-quality reads were mapped to the reference genome (described above) with Bowtie2 (v2.3.5.1)79; binarization, sorting, and indexing were done with Samtools (v1.10)80; gene counts were computed with the featureCount function of Subread (v2.0.1)81.

Differential expression analysis was performed with DESeq2 (v1.30.1)82. In brief, DESeq2 normalizes the raw read counts with normalization factors (”size factors”) to account for differences in sequencing depth between samples. Subsequently, gene-wise dispersion estimates are computed for each gene separately using maximum likelihood, and then shrunk toward the values predicted by the dispersion-mean dependence curve to obtain final dispersion values. Finally, DESeq2 fits a negative binomial model to the read counts and performs significance testing using the Wald test. Here reported p-values result from the Wald test of read counts from the digested alginate condition compared to the alginate condition and were adjusted for multiple testing by Benjamini-Hochberg correction83 as implemented in the p.adjust function of base R (v4.1.2). The reported log2 fold changes indicate the log2(DESeq2-normalised reads in digested alginate condition / DESeq2-normalised reads in alginate condition) for each gene.

Visualization of gene maps was performed in R with the ggplot2 package (v3.4.0) and the extension gggenes (v0.4.1) by David Wilkins.

For systematic functional analysis we performed gene set enrichment analysis (GSEA)23 using the fgsea function of the fgsea package (v1.20.0) with minimal number of unique genes per gene set “minSize” = 5 and number of permutations “nPermSimple” = 1000000. In brief, GSEA takes the full gene list ranked by log2 fold change and annotated with K numbers as input and determines whether the member genes of any KEGG pathway are randomly distributed throughout the ranked gene list or whether they are primarily found at the top or bottom23. This is quantified by the enrichment score (ES), which corresponds to a weighted Kolmogorov-Smirnov-like statistic. The ES of each gene set is normalized to the mean enrichment of random samples of the same size to account for the size of the set, yielding the normalized enrichment score (NES). To estimate the significance level of the enrichment score, the p-value of the observed enrichment score is calculated relative to a null distribution that was computed from permuted data. The estimated significance level was adjusted to account for multiple hypothesis testing. As gene sets we chose all KEGG pathways and KEGG BRITE categories (as noted in Table S1 in column “KEGG_pathway”) within all genes of V. cyclitrophicus ZF270 annotated with a KEGG Ontology identifier (”K number”). Visualization of differential expression levels in KEGG pathways was performed with the R package pathview (v1.35.0)84. We also formed gene sets of the genes associated with alginate utilization (see “Sequencing and gene annotation of Vibrio cyclitrophicus ZF270”, Table S16), of the genes of the flagellar locus that map to the KEGG pathways “Bacterial motility proteins”, and of the genes of the flagellar locus that map to the KEGG pathways “Bacterial chemotaxis” (as noted in Table S17) within all genes of V. cyclitrophicus ZF270.

Image analysis

Cells within microscopy images were segmented and tracked using ilastik (v1.3) (”pixel classification workflow” and “tracking with probabilities workflow”). Phase contrast images were used for alignment, segmentation, tracking and linking. Images were cropped at the boundaries of each microfluidic chamber. The lineage identity of each single cell was assigned by ilastik’s tracking plugin and visualized by coloring the segmented cells, respectively. The growth rate of each cell was computed as the change of cell area over time, i.e., via a linear regression of the single-cell area over the time between consecutive cell divisions, based on ilastik’s segmentation and tracking output. Cells that were tracked over less than three frames were excluded. Measurement of swimming speeds and displacement of cells was performed using ilastik (v1.4), ImageJ (v2.3) and Trackmate (v7.5.2). Briefly, cell-segmentation (”pixel classification workflow” in ilastik) and tracking (”animal tracking workflow” in ilastik) were performed using the high frame rate phase contrast images in ilastik (v1.4). Cell trajectories and properties were then computed using the output of the ilastik workflow in Trackmate.

Dispersal analysis

For the analysis on the dispersal of cells (Fig. 2), we computed the cell number as the total number of cells within a microfluidic chamber. The change in the number of cells was computed by subtracting the number of cells before the medium switch (i.e., average number of cells between t = 1.9 to 2.1 h) from the number of cells after the medium switch (i.e., average number of cells between t = 3.9 to 5.5 h).

Datasets and statistical analysis

All batch experiments were replicated 3 to 6 times. Growth curves were analyzed in Python (v3.7) using the Amiga package (v1.1.0)85 and GraphPad Prism (v8, GraphPad Software). The microscopy dataset consisted of eight chambers each, corresponding to the eight replicates shown in Figure 1 and Figure 2. These were grouped into two biological replicates wherein each biological replicate was fed by media through a unique channel in a microfluidic chip. Cells with negative growth rates were excluded from the analysis after visual curation, as they represented artifacts, mistakes in segmentation or linking during the tracking process, or non-growing deformed cells. Each chamber was treated as an independent replicate. Comparisons were considered statistically significant when p < 0.05 or when the False Discovery Rate (FDR)-corrected q was smaller than 0.05. FDR corrections were applied when multiple t tests were performed for the same dataset. Measures of effect size are represented by the R2 or eta2 value. All statistical analyses were performed in GraphPad Prism v9.0 (GraphPad Software, USA), R v4.1.2, RStudio v1.1.463 (Posit, USA).

Acknowledgements

We thank past and present members of the Microbial Systems Ecology group for feedback. This research was supported by an ETH fellowship and a Marie Skłodowska-Curie Actions for People COFUND program fellowship (FEL-37-16-1) to GD; an ETH Career Seed Grant to GD (FEL-14 18-1), the Simons Foundation Collaboration on Principles of Microbial Ecosystems (PriME, #542379 and #542395) to MA and RS; and by ETH Zurich and Eawag.

Data Availability

Sequencing data is available on NCBI, BioProject PRJNA991487, upon publication. All further data and code is deposited on ERIC Open (https://opendata.eawag.ch/, https://doi.org/10.25678/0008MH), available upon publication. Code will also be made publicly accessible on https://github.com/Microbial-Systems-Ecology/publications upon publication.

Competing interests

The authors declare no competing interests.

Author Contributions

GD, AS, JS, and JK conceived the research along with MA. GD designed and performed the experiments. AS developed and implemented the algorithms to compute single-cell growth rates and visualize lineage identities from segmented and tracked cells. AS analyzed and visualized the RNA-seq data. JK conceived and implemented the cell dispersal analysis. EC conducted the ISCA experiment. SP carried out the mass spectrometry measurements. AS and GD analyzed the data with inputs from JK, JS, CM, MA, and RS. AS wrote the manuscript with inputs from GD, OS, MA, JK, JS, EC, SP, CM, and RS.