Strains belonging to the same species display distinct growth dynamics on the marine polysaccharide alginate
We first quantified the growth dynamics of the 12 Vibrionaceae strains (Supplementary Table 1) on alginate in well-mixed batch cultures. Growth of populations was initiated at approximately the same inoculum density (105 colony forming units (c.f.u.) ml−1). We tracked the growth dynamics by measuring the optical density at 600 nm and compared the maximum population size reached over the course of 36 h (Fig. 1 and S1). We found significant differences in the maximal optical density achieved by different strains within each species (Fig. 1 and S1). In V. splendidus, strains 12B01 and FF6 reached a lower maximum population size compared to strains 1S124 and 13B01 (Fig. 1 and S1A). In V. cyclitrophicus, strain ZF270 reached a lower maximum population size compared to strains 1F175, 1F111, and ZF28 (Fig. 1 and S1A). Similarly, in V. sp. F13, strain 9ZC77 reached a lower maximum population size than strains 9CS106, 9ZC13, and ZF57 (Fig. 1 and S1A). These findings suggest that some strains are limited in their growth abilities in well-mixed environments, perhaps as a consequence of differences in the amount and activity of enzymes they release (Supplementary Table 1).
Maximum optical density (measured at 600 nm) achieved by populations of strains belonging to Vibrio splendidus, Vibrio cyclitrophicus, and Vibrio sp. F13 during the course of a 36 h growth cycle on the same concentration (0.1% weight/volume) of the polysaccharide alginate. Points and error bars indicate the mean of measurements across populations within each ecotype (npopulations = 3) and the 95% confidence interval (CI), respectively. Different letters indicate statistically significant differences between strains within one species (One-way ANOVA and Dunnett’s post-hoc test; V. splendidus: p < 0.0001, F = 325.9, R2 = 0.99; V. cyclitrophicus: p = 0.0099, F = 7.60, R2 = 0.74; V. sp. F13: p < 0.0001, F = 245, R2 = 0.98). For growth curves of individual strains, see Fig. S1.
Differences in enzyme secretion are associated with differences in the growth dynamics of alginate degraders
We hypothesized that reduced secretion of alginate lyases is responsible for limiting the growth of strains that achieve low optical densities in batch culture (Fig. 1). To test our prediction, we measured the level of secretions of alginate degrading enzymes by the different strains using a plate assay. In this assay, alginate reacts with iodine to produce a violet compound. The breakdown of alginate by secreted enzymes that diffuse away from the colony creates simpler oligomers or monomers that yield colorless halos when stained with iodine (Fig. S2A). Therefore, the level of enzyme secretion and activity can be quantified by measuring the size of the alginate-free halo generated by cell colonies as they grow on an alginate-containing plate (see Methods and Fig. S2A) [21]. Strains that do not produce halos but can grow on alginate likely possess enzymes that are bound to the membrane or cell-surface. While this assay provides a quantitative measure of the activities of secreted enzymes, it cannot be used to determine the amount of enzymes secreted, since a lower quantity of highly active enzymes and a higher quantity of less-active enzymes could give the same clearing radius.
Enzyme secretion as measured by halo diameter varied among strains within each species. Among V. splendidus strains, 12B01 and FF6 had the lowest halo diameters (Fig. S2B), while among V. cyclitrophicus strains, ZF28 had a smaller halo diameter than the other strains (Fig. S2B). In contrast to the other V. sp. F13 strains, 9ZC77 did not produce a detectable halo (Fig. S2B). This observation suggests that 9ZC77 either does not produce substantial enzymatic secretions or the alginate lyases are primarily membrane- or cell-surface bound. Comparison of the results from the plate assays (Fig. 2A, Fig. S2) and the batch cultures (Fig. 1) reveals that the level of secretion of alginate lyases by the different strains is significantly correlated with the maximum population size they achieved when growing on alginate under well-mixed conditions (Fig. 2A). This correlation supports our hypothesis that growth in batch culture is limited by the ability to produce secreted alginate lyases.
A Correlation between maximum optical density in batch culture (Fig. 1) and halo diameters in the plate assays across the 12 strains. The line indicates the fit of a linear regression (slope = 0.05, R2 = 0.67, P < 0.0001; Spearman’s correlation: r = 0.90, P = 0.0002). B Supplementation with external alginate lyase reduces the time to reach exponential growth phase of strains that secrete low levels of alginate lyase to a greater extent compared to high secretors of alginate lyase. Circles represent the mean of the populations of each ecotype (npopulations = 3) and error bars indicate 95% confidence intervals. Letters indicate statistically significant differences between strains within each species (One-way ANOVA with Dunnett’s posthoc test; V. splendidus: F = 114.0, p < 0.0001, R2 = 0.97; V. cyclitrophicus: F = 2.36, p = 0.14, R2 = 0.46; V. sp. F13: F = 47.27, p < 0.0001, R2 = 0.94). Open circles indicate non-significant differences from 0, indicating no reduction in the time to reach exponential phase; filled circles indicate significant differences from 0 (One sample t-test; 12B01: t = 24.4, df = 2, p = 0.001; FF6: t = 16.1, df = 2, p = 0.003; 1S124: t = 53, df = 2, p = 0.004; 13B01: t = 3.2, df = 2, p = 0.08; 1F175: t = 3.8, df = 2, p = 0.06; 1F111: t = 1.9, df = 2, p = 0.19; ZF28: t = 19, df = 2, p = 0.002; ZF270: t = 3.8, df = 2, p = 0.006; 9CS106: t = 23, df = 2, p = 0.0019; 9ZC13: t = 0.4, df = 2, p = 0.67; 9ZC77: t = 19.6, df = 2, p = 0.002; ZF57: t = 2.4, df = 2, p = 0.13). See Fig. S3 for growth dynamics of each ecotype on alginate as well as correlations between enzyme secretion and growth kinetics.
To further test our prediction that alginate lyase production capability constrains growth, we determined whether externally supplementing alginate lyases could alleviate the growth limitation of the strains that are low enzyme secretors. We grew the 12 strains individually as batch cultures in alginate medium that was supplemented with a commercial mixture of alginate lyases and measured their growth dynamics. The results show that alginate lyase supplementation did not alter the maximum population size reached by the low secreting strains (Fig. S3A). Commercial alginate lyases likely produce a mixture of oligomers and not just monomers. Therefore, cells of lower-growing alginate degraders could still be limited by the availability of monomeric hydrolysis products. We then calculated the time to reach exponential growth phase of each ecotype when growing in the presence of external alginate lyases and compared this with the time taken in the absence of this external supplementation. We found that the strains with a low level of enzyme secretion and activity, such as 12B01, FF6 and 9ZC77, showed the largest reduction in the time required to reach exponential phase upon supplementation with external alginate lyases (Fig. 2B and S3B, C). However, strains 12B01 and FF6 still have very low optical densities in these culture conditions, even in the presence of externally added enzymes, suggesting factors other than enzyme production limit the growth of these strains.
Strains with relatively higher enzyme secretions such as 9CS106, ZF57 and 9ZC13 showed varying levels of benefits from enzyme supplementation. It is known that such differences can emerge due to variations in enzyme catalytic properties, transport efficiencies and physiological properties of the different strains [5, 22]. In contrast to the effect of enzyme supplementation for V. splendidus and V. sp. F13, the effect of enzyme supplementation was more limited in V. cyclitrophicus (Fig. 2A), the species in which differences in enzyme activity were most limited (Figs. S2B, S3A, D). Taken together, our results suggest that strains that secrete low levels of alginate lyases benefit the most from external alginate lyase supplementation, indicating that the growth rate of these strains is limited by the ability to produce alginate lyases.
Strains with low levels of enzyme secretion display higher levels of aggregation in environments that favor spatial structure
Cell aggregation can reduce the growth limitations observed on polysaccharides in absence of supplementation and based on native enzymatic ability of strains [16]. Specifically, we expect that by reducing dispersal and forming aggregates, cells could reduce the diffusional loss of alginate lyases or breakdown products and thereby could alleviate growth limitation of strains with low secretion levels (we use the term “aggregate” here without implying that cells are actively moving towards each other; cell aggregates can also be formed if cells do not disperse after cell division, as observed previously [16]. We tested this prediction on a subset of the original 12 Vibrionaceae strains. From each species, we selected a representative ecotype producing high levels of alginate lyase (V.splendidus: 13B01; V. cyclitrophicus: 1F111; and V. sp. F13: 9ZC13) and one producing low levels of alginate lyase (V.splendidus: 12B01; V. cyclitrophicus: ZF270; and V. sp. F13: 9ZC77). To test if strains aggregate when growing on the polysaccharide alginate, we grew cells of the different strains within microfluidic growth chambers. Within these chambers cells receive a constant supply of alginate and divide, either forming aggregates by growing as monolayers or dispersing after division. This approach enabled the quantification of growth and movement using time-lapse microscopy coupled to automated image analysis.
We found that aggregative behaviors were common in all the strains across the three species tested (Fig. 3A–F). By following the growth dynamics of individual cells using cell segmentation and tracking, we were able to map lineages of all cells present in each growth chamber based on division events. This analysis revealed that cells emerging from the same founder cell tended to remain close to each other (Fig. 3A–F, Supplementary Video 1–6). However, for all three species, the low-secreting strains, (Fig. 3G–I), displayed a higher degree of aggregation compared to the corresponding high-secreting strains, 13B01 (Fig. 3G), 1F111 (Fig. 3H), and 9ZC13 (Fig. 3I). The maximum number of cells within aggregates in a chamber was on average 5.54, 3.21, and 1.89 times higher for the low-secreting strains of V. splendidus, V. cyclitrophicus, and V. sp. F13, respectively, compared to the corresponding high-secreting strains (Fig. S4). These findings indicate that cells of low-secreting strains, at least in the context of the three species-pairs tested, form larger aggregates while growing on alginate.
Representative images of microfluidic growth chambers containing cells of (A) Vibrio splendidus 13B01, (B) Vibrio splendidus 12B01, (C) Vibrio cyclitrophicus 1F111, (D) Vibrio cyclitrophicus ZF270, (E) Vibrio sp. F13 9ZC13, and (F) Vibrio sp. F13 9ZC77. Cells that originate from the same founder cell are shown in the same color (based on cell tracking through divisions). The three high-secreting strains (13B01, 1F111 and 9ZC13) are shown with orange panel outlines and the three low alginate lyase secreting strains (12B01, ZF270 and 9ZC77) are shown with blue panel outlines. G–I The increase in the number of resident cells in growth chambers over time for strains of (G) V. splendidus, (H) V. cyclitrophicus and (I) Vibrio sp. F13. The levels of enzymatic secretions for the strains (from Fig. S2B) are shown in the inset of panels G–I. The number of cells in a chamber depends on both the rate at which cells divide and the fraction of cells that leaves the chamber. Circles indicate the number of cells present at a given time point in each chamber (nchambers = 7), with a logistic growth regression line for each chamber: 13B01, R2 = 0.97–0.99, max cells (mean ± standard deviation) = 156.90 ± 83.21, k = 0.09–0.20 h−1; 12B01, R2 = 0.99 in all cases, max cells = 858.70 ± 67.78, k = 0.18–0.20 h−1; 1F111, R2 = 0.89–0.99, max cells = 410.40 ± 106.40, k = 0.10–0.53 h−1; ZF270, R2 = 0.86–0.99, max cells = 1318 ± 198.70, k = 0.10–0.30 h−1; 9ZC13, R2 = 0.91–0.99, max cells = 545.40 ± 50.69, k = 0.08–0.21 h−1; 9ZC77, R2 = 0.98–0.99, max cells = 1033 ± 72.12, k = 0.36–0.40 h−1). See supplemental videos S1–S6 for timelapse image series of cells.
Density dependence is stronger in strains with low alginate lyase secretions
Given that low secretion by strains is associated with greater aggregation, we then quantified if the differential aggregation behaviors of strains influence growth abilities of individual cells on the polysaccharide alginate. For this, we measured the growth rate of individual cells within microfluidic chambers using automated image analysis (Fig. S5A–C). We then analyzed the relationship between the median growth rate of all the cells present in a chamber during a given time interval and the total number of cells that were present in the chamber during that interval. For cells of all six strains across three species, the median growth rate increases with the increase in the number of resident cells in a chamber up to a point, and then decreases as the number of cells increases further (Fig. 4). However, the number of cells required to achieve the maximal growth rate differed between the strains. For two of the three species tested, the maximal growth rate of the low-secreting ecotype occurred at a significantly higher cell density than that of the high-secreting ecotype (Fig. 4A, B, E, F). The comparison for V. cyclitrophicus did not follow this trend, with no significant difference between the two V. cyclitrophicus strains in the cell density associated with the maximal growth rate (Fig. 4C, D), likely because the two strains are much more similar in their level of secretion in comparison with the other species (Fig. 2A). The increased cell density due to aggregation within the microfluidic chambers allows cells to overcome growth limitations that strains face in well-mixed environments (Fig. S6A). We find that high-secreting strains have similar or increased final growth rates compared to the low-secreting strains when growing in microfluidic chambers. This indicates that reduced group sizes of high-secreting strains (Fig. 3G–I) arise as a result of increased dispersal rather than reduced growth rates. Well-mixed environments, as in the batch-culture experiments, limit the growth rate of the two strains with the lowest levels of enzyme secretion, 12B01 and 9ZC77, relative to their high-secreting counterparts, 13B01 and 9ZC13, respectively (Fig. S6B). In contrast to growth in well-mixed environments, growing in spatially structured environments, such as within the microfluidic growth chambers, 12B01 and 9ZC77 are able to improve growth substantially relative to 13B01 and 9ZC13, respectively (Fig. S6B). The growth of 12B01 in microfluidic chambers is especially interesting since the growth of this strain was negligible in well-mixed environments (Fig. 1) even though it produced enzymes (Fig. S2B). These observations are indicative of other yet unknown biological factors that likely influence growth on polymers. In addition, 9ZC77 cells had no noticeable enzyme secretions (Fig. S2B) but displayed density dependent growth dynamics. This finding suggests that strains with cell-bound enzymes have a disadvantage in well-mixed environments but can benefit from surface-association. The growth rate of V. cyclitrophicus ecotype ZF270 relative to that of ecotype 1F111 is similar in both well-mixed and microfluidic growth environments (Fig. S6B). These findings indicate that surface association and aggregate formation allow strains, especially those that secrete very low levels of alginate lyases, to increase cell density and enables cells to benefit from each other’s degradation activities, thereby allowing intercellular synergy amongst cells within aggregates in spatially structured environments. The increased cell density and surface-association thus enables cells in spatially structured environments to overcome growth limitations that arise when growing on polysaccharides in well-mixed environments.
A–C Median single-cell growth rate (h−1) of (A) Vibrio splendidus, (B) Vibrio cyclitrophicus and (D) Vibrio sp. F13 cells as a function of the number of cells growing together in a microfluidic chamber. Cells were binned into 2-hour intervals based on their birth times over a 28 h experiment. This approach also yielded the number of cells present during growth in each time interval. The relationship between median growth rate of all cells present in a 2 h interval in a given chamber and the number of cells in the chamber was determined using a non-linear regression model for each chamber. The levels of enzymatic secretions for the strains (from Fig. S2B) are shown in the inset of panels. High-secretion strains are indicated in orange whereas relatively low-secretion strains are indicated in blue. Circles represent data for a single bin from one chamber, and lines indicate the trajectory of growth rates for each chamber (nchambers = 7; 13B01, R2 = 0.52–0.76; 12B01, R2 = 0.66–0.95; 1F11, R2 = 0.88–0.91; ZF270, R2 = 0.71–0.93; 9ZC13, R2 = 0.73–0.95; 9ZC77, R2 = 0.57–0.80). D–F Number of cells required to reach half of the maximum growth rate for strains of (D) V. splendidus, (E) V. cyclitropicus, and (F) V. sp. F13. The low-secreting V. splendidus ecotype 12B01 required a five-fold higher number of cells to achieve maximal growth rates compared to the high-secreting 13B01 cells (D; 13B01 median 8.99 cells and 12B01 median 45.25 cells). In the case of V. cyclitrophicus, both strains required a similar number of cells to reach maximal growth rates (E; 1F111 median 25.48 cells and ZF270 median 41.34 cells). For V. sp. F13, the low-secreting ecotype 9ZC77 required an eight-fold higher cell number to achieve maximal growth rates compared to the high-secreting ecotype 9ZC13 (F; 9ZC13 median 23.90 cells and 9ZC77 median 210.40 cells). Asterisks indicate statistically significant differences (P < 0.001) between strains (Mann–Whitney tests, nchambers = 7; V. splendidus, p = 0.008, Hodges–Lehmann difference = 36.75; V. cyclitrophicus p = 0.15, Hodges–Lehmann difference = 19.85; V. sp. F13 p = 0.003, Hodges–Lehmann difference = 188.5). Box plots extend from the 25th to 75th percentiles and whiskers indicate the 10th and 90th percentiles of median growth rates.
Once the growth rate has reached its maximum value, the growth rate then plateaus for V. splendidus strains or begins to decrease for V. cyclitrophicus and V. sp. F13 strains with further increase in the cell numbers within the growth chambers (Fig. 4A, C, E). We measured the peak cell densities reached within the growth chambers (13B01: 250.30 ± 129 (mean cell number across replicate chambers ± standard deviation), 12B01: 1439 ± 98.95, 1F111: 255 ± 42.85, ZF270: 418.60 ± 34.11, 9ZC13: 228.80 ± 39.61 and 9ZC77: 358.5 ± 109.10) and found that the benefit of intercellular interactions is balanced (for V. splendidus) or outweighed (for V. cyclitrophicus and V. sp. F13) by the negative effects of competitive interactions between cells within aggregates.
Since maximum growth rates are dependent on cell densities, we probed if an increased initial inoculum density would allow low-secreting strains to increase growth rates in well-mixed batch culture environments. In line with our findings in microfluidic chambers, low-secreting strains 12B01 and 9ZC77 benefit from a higher inoculum density (Fig. S7A–C). An increased inoculum density reduces the growth rate of high-secreting strains in well-mixed environments (Fig. S7A–C). These findings lend further credence to our idea that high cell densities due to aggregate formation can benefit low enzyme secreting strains.
A mathematical model offers a mechanistic explanation for the role of cell density and enzyme secretion on polymer degradation and growth rates
To obtain a mechanistic understanding of the roles of cell density and enzyme secretion rate on growth rates of strains, we developed an individual-based model of the growth rate of cells on polysaccharides taking into account cell density, enzymatic breakdown and loss of breakdown products. Our goal was to describe a system composed of spatially structured populations in nature similar to cells growing as a monolayer in microfluidic growth chambers. We modeled a 2D system (Supplementary Information and Table S2) where cells occupy sites on a 50 × 50 grid (Fig. 5A). We assumed diffusion as the mode of transport within the 2D chamber for polymeric substrates, secreted enzymes, and breakdown products. The loss of breakdown products (oligomers) and enzymes from the opening of the chamber by flow is simulated by assuming constant zero concentrations as the boundary condition. Individual cells are initially uniformly distributed on the 2D lattice, where they secrete enzymes that breakdown polysaccharides, consume oligomers, and grow in response to local oligomer concentrations. Monod-type kinetics were used to model substrate uptake and cell growth (see methods for details). The rate of enzyme secretion, SE, was assumed to follow first-order kinetics as a function of cell biomass, SE = KenzB, as previously observed experimentally, where Kenz is the enzymatic activity and an intrinsic cellular property and B is the individual cell biomass.
A Schematic representation of the mathematical model that simulates microfluidic chambers as a 2D lattice with 50 × 50 µm grid. Individual cells occupy single grid cells where they act as sink terms for the oligomer diffusion–reaction model. A Dirichlet boundary condition is used at the chamber periphery to simulate the effect of flow, which allows for constant concentrations of polysaccharide (0.1 mg/L) and oligomer (0 mg/L) at the boundary. Cells are uniformly distributed in the chamber. B Average growth rates of individual cells 20 h after the start of the simulations as a function of enzymatic activity, Kenz, and cell density, ρ. The growth rate was normalized to the maximum values obtained across the full set of conditions. C Oligomer profiles and single-cell growth rates for simulated chambers after 20 h. Two scenarios are shown, taken from simulations in Fig. 5B (I, ρ = 0.3 µm−2 and Kenz = 0.7 hr−1; II, ρ = 0.3 µm−2 and Kenz = 0.1 h−1). The values of growth rates and oligomer concentrations are normalized to the maximum values observed in the simulations to allow comparison. D Depolymerization and E growth rates of polysaccharides as a function of the diffusion coefficient of oligomers, polysaccharide, and enzymes. The diffusion coefficient was normalized by the diffusion coefficient in free water. Simulation parameters are described in detail in Table S2.
We asked whether our model predicts similar density-dependent growth dynamics to those observed in the experimental observations. We simulated bacterial growth in the 2D chambers and calculated the mean growth rates over a wide range of initial cell densities and enzymatic activities. Consistent with our experimental observations, the model predicted that the growth rate of cells with lower enzymatic activity is higher at an increased cell density (Fig. 5B). In comparison, cells with a higher enzymatic activity through increased secretions have higher growth rates at low densities (Fig. 5B). To understand the biophysical mechanism underlying these effects, we examined the spatial profile of oligomer concentrations simulated by the model and found that lower enzymatic activity leads to the accumulation of oligomers and supports higher growth rates for cells regardless of their spatial arrangement. However, in higher enzymatic activity, oligomers arising as result of polysaccharide degradation are increasingly lost to diffusion without being available for bacterial cells to uptake (Fig. 5C). Our results indicate that the rate of oligomer loss to diffusion and uptake plays a dominant role in shaping the nature of intercellular interactions in polysaccharide-degrading populations. We thus performed simulations where the diffusion coefficients of polymers and oligomers were systematically varied relative to diffusion in free water to yield low rates of diffusion such as in biofilm or cell aggregations and high rates of diffusion similar to boundary conditions in a quiescent region in contact with flow. The model suggests that a scenario with a diffusion coefficient that is an order of magnitude lower than that in bulk water results in the highest average growth rate for the bacterial population (Fig. 5DE). A higher diffusion rate that resembles flow, while increasing the polymer degradation rate (Fig. 5D), also enhances the loss of oligomers through diffusion, making them unavailable to bacterial cells and thus reducing the growth rate. A similar effect emerges at higher cell densities, where faster polymer degradation by cells at the opening of the growth chambers can result in loss of oligomers, thereby making them unavailable to other cells and leads to negative density dependent growth effects (Fig. 5C, E). In absence of experimental measurements, constant physiological properties for maximum growth rate, half-saturation constant and yield of converting oligomers to biomass were assumed for all bacterial cell types. Recent studies had shown variabilities in half-saturation constant across bacterial strains [5] and further studies are required to investigate the role of these physiological differences on bacterial collective behavior.
Our results indicate that enzymatic secretion capabilities can influence growth dynamics and intercellular interactions within bacterial populations growing on polysaccharides. Differences in enzyme secretion levels can arise as a result of strains harboring functionally distinct alginate lyases or as a result of distinct localization of alginate lyases [18, 22,23,24]. For instance, as previously explained V. splendidus 12B01 has alginate lyases that are largely intra-cellular or membrane-bound whereas alginate lyases of V. splendidus 13B01 are largely secreted into the extracellular environment [18, 25]. In addition, strains can harbor distinct classes of alginate lyases and the breakdown products generated by these enzymes can differ [26]. Our findings suggest that cells that are genetically predisposed to produce low levels of enzymes or have predominantly membrane-bound alginate lyases can overcome growth disadvantages that exist in well-mixed environments by aggregating. These collective behaviors allow cells to generate 5–10-fold higher local cell densities enabling growth at maximal rates. Differences in growth patterns between rumen-derived strains of Bacteroides thetaiotaomicron on the polysaccharide yeast α-mannan have also been shown to emerge as a consequence of variation in substrate uptake affinities [5]. Our finding that supplementation of alginate lyase, which should increase the pool of oligo- or monomeric hydrolysis products, did not increase maximum growth yields of the low alginate lyase secreting strains suggests that uptake capabilities of strains could also play a part in the observed growth constraints on polymers.
Polysaccharides like alginate or chitin exist predominantly as particulates in natural ecosystems [27] and offer surfaces that enable bacterial cells to form dense colonies. It is known that surface association can induce the expression of polysaccharide utilization genes in bacterial populations [28]. In addition, aggregation behaviors increase local cell densities and enable intercellular interactions that can allow individuals to benefit from each other’s metabolic activities [16, 29]. Strains that have increased enzymatic secretions can breakdown polysaccharides by adopting dispersed growth modes, that can facilitate migration of cells to new polysaccharide particles in nature. Additionally, in communities where multiple strains coexist [7, 8, 20, 30], the degradative activity of high enzyme secretors can create public goods that can positively influence the growth of not only the low enzyme producers but also strains that lack the ability to degrade polysaccharides and rely on cross-feeding of monomers from the polymer degraders [18, 31, 32].
Dispersal and collective behaviors have been shown to evolve together in multiple organisms [33, 34]. Differences in mobility can dramatically influence the likelihood of individuals in a population engaging in cooperation [33]. When the cost of aggregation behaviors is high, cooperation between individuals of a population is favored when individuals disperse less [34]. Less costly behaviors are relatively less sensitive to dispersal of individuals [34]. While aggregation behaviors do provide a benefit to individuals that require high densities, increased dispersal can allow individuals that produce higher levels of public goods to avoid invasion by consumer cells which do not engage in public good production [35].
The cellular processes that link aggregation, motility and polymer metabolism in marine bacterial populations are yet unclear. The formation of bacterial aggregates and production of polysaccharide-degrading enzymes both require substantial cellular resources and thus likely impose substantial metabolic costs on cells [10, 36]. This will result in cells facing trade-offs between allocating resources to enzyme production versus aggregation [14], both of which enable cells to grow on polymers. A potential consequence is that marine polysaccharide degraders are constrained by trade-offs such that they forego enzymatic capabilities in favor of aggregation or alternately acquire genes encoding enzymes that can either be produced in large amounts and also be secreted instead of aggregate formation capabilities. Aggregation of bacterial cells is known to be coordinated through quorum sensing [37, 38]. However, the molecular mechanisms that regulate collective behaviors in response to the nutrient composition in the natural environment are yet unclear. Our work is a first step in this direction and provides a direct link between the dynamics of enzyme secretion, growth physiology and collective behaviors of polysaccharide degrading bacterial populations.
Source: Ecology - nature.com
