Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2002, p. 3996-4006, Vol. 68, No. 8
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.8.3996-4006.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
and Kam Tang1
Danish Institute for Fisheries Research, DK-2920 Charlottenlund,1 Marine Biological Laboratory, DK-3000 Helsingør, Denmark,3 Institute of Chemistry and Biology of the Marine Environment, University of Oldenburg, 26111 Oldenburg, Germany2
Received 4 March 2002/ Accepted 23 May 2002
|
|
|---|
|
|
|---|
In resolving the population dynamics of bacteria on aggregates one needs to consider the rate at which bacteria colonize the aggregates. Several studies have described the succession of microbes on aggregates (19, 26, 36, 41, 43), but few have examined the colonization process at a relevanti.e., short (minutes)time scale (15), and none have adequately considered the effect of the flow environment on colonization. A significant fraction of pelagic bacteria are motile (15, 18, 30), and simulation models have shown that motility and chemosensing capability are important for colonization, even of sinking aggregate (21, 23). Here we develop encounter models to predict colonization rates based on observed motility patterns of the bacteria and the hydrodynamic environments. We then compare model predictions with empirical measurements by using artificial aggregates.
|
|
|---|
![]() | (1) |
(i) Steady state.
The encounter rate kernel for scavenging, assuming a low Reynolds number flow (25), is:
![]() | (2a) |
![]() | (2b) |
![]() | (2c) |
![]() | (3) |
is the average run length (time) and
the mean value of the cosine of the angle between successive runs. The Sherwood number is the ratio between total (diffusive plus advective) transport of bacteria toward the aggregate and the transport due to diffusion alone (random walk motility). Sh is "1" in the absence of advection (i.e., stationary aggregates). For sinking aggregates, Sh can be approximated as follows (24):
![]() | (4) |
is the kinematic viscosity (
10-2 cm2 s-1). The Sh of typical marine snow aggregates in the size range of 0.1- to 1-cm radius varies from ca. 5 to 20 (24). A comparison of ß-values for the scavenging and motility colonization mechanisms demonstrates that colonization is dominated by bacterial motility, and we can thus ignore scavenging when considering motile bacteria.
Our experiments revealed that bacteria that have colonized an aggregate might detach (see Results). If we assume that a constant fraction per unit time (
) of the bacteria that have colonized the aggregate leave it again, then equation 1 can be modified to:
![]() | (5) |
![]() | (6a) |
![]() | (6b) |
(ii) Non-steady state.
For a sinking aggregate, transport of random walk bacteria toward the aggregate reaches steady state almost immediately (24). For bacteria swimming along straight lines, steady state is instantaneous. However, for a nonsinking aggregate colonized by random walk bacteria, steady state is approached slowly. The time-dependent encounter rate kernel for randomly walking bacteria is (33)
![]() | (7) |
10-6 to 10-5 cm2 s-1) and sphere sizes (a
0.2 cm) typical of our experiments. Obviously, steady state is not a valid assumption in still-water experiments. The cumulated number of bacteria that have colonized the particle at time t (Nt) under non-steady-state conditions is then
![]() | (8a) |
![]() | (8b) |
Model predictions.
We use equation 8a to describe the colonization process by randomly walking bacteria in still water where a steady state cannot be assumed. Equation 8a predicts that accumulation of bacteria on aggregate increases with the size of the aggregate (a) and diffusion coefficient of the bacteria (D), which in turn is a function of their swimming speed, run lengths, and turn angles (equation 3).
In all other cases (linear swimming, flow), we use equation 6a to describe the colonization process. It follows from equation 6a that the accumulation of bacteria on aggregate increases with the encounter rate kernel (ß), which is related to bacterium size, bacteria's motility patterns, aggregate sizes, and ambient flow rate (equations 2a, 2b, and 2c).
To test model predictions, we observed the swimming behavior of 10 bacterial strains (Table 1) and experimentally measured their rates of colonizing artificial aggregates (stationary and sinking). All strains were isolated from marine aggregates or diatoms collected in the Øresund, Denmark, and identified by means of 16S ribosomal DNA sequencing (H.-P. Grossart, unpublished data). We maintained bacterial cultures at exponential growth in marine broth (MB [MB2216; Difco]; 15 g liter-1). In most experiments, we aimed at an ambient concentration of 106 ml-1.
|
View this table: [in a new window] |
TABLE 1. Bacterial strains used in the experiments and the number and types of experiments conducted
|
Colonization experiments.
Bacterial colonization was examined by monitoring the accumulation of bacteria on artificial aggregates that were incubated in still or flowing seawater containing bacteria. Model aggregates were made from 2% agar as described by Cronenberg (12). Warm agar was dripped from a pipette into sterile seawater covered by a thin layer of paraffin oil. Almost-perfect agar spheres formed in the oil and subsequently sank into the seawater. Traces of paraffin oil on the spheres were rinsed off with sterilized seawater. The aggregate size (0.12- to 0.28-cm radius) was controlled by the width of the pipette mouth. In our standard experiments, we used spheres with a radius of 0.2 cm. In some experiments we added MB (15 g liter-1) or dimethylsulfoniopropionate (DMSP; 1 mM) to the agar. DMSP is a natural organic substrate leaking from aging algal cells and aggregates (27). In a few experiments we used diamond-cut glass pearls as model aggregates (0.15- to 0.38-cm radius).
A typical experiment consisted of simultaneously placing 24 aggregates in a bacterial suspension. Agar-aggregates were fixed on thin glass needles and suspended in the bacterial suspensions in either still water (2 liters) or in a flume (ca. 10 liters; flow velocities 0.1 to 0.3 cm s-1). Glass pearls were suspended on thin wires of stainless steel. Triplicate aggregates were collected typically at 0, 5, 10, 15, 20, 40, 80, and 160 min. Sampled aggregates were placed in a counting chamber made of an O-ring glued onto a microscopic slide, and a drop of formaldehyde and a drop of DAPI (4',6'-diamidino-2-phenylindole) were added. The chamber was then closed by a coverslip, and the bacteria were counted at x1,000 magnification under epifluorescence. Some bacteria may be lost when collecting aggregates and, hence, all abundance estimates are conservative. However, in one experiment we examined the effect of retrieving aggregates and found no measurable effect. Ambient concentrations of bacteria were measured at the onset and termination of each incubation experiment by filtering 1-ml suspensions onto 0.2-µm (pore-size) filters, staining with DAPI, and counting under epifluorescence.
We conducted 44 experiments in which we compared the different strains and examined the effects of flow, aggregate size, bacterial concentration, and the presence of chemical signals (MB and DMSP) in the aggregates on colonization pattern and rate (Table 1).
Observed colonization rates were fitted to either equation 8a (random walk bacteria in still water) or equation 6a (all other cases) by nonlinear regression routines to derive estimates of
and D (equation 8a) or
and ß (equation 6a). Unfortunately, the two parameters (
and D or
and ß) are strongly correlated, which results in large confidence limits. From the estimates of ß, we finally estimated D · Sh (equation 2c) or bacterial swimming speed, u (equation 2b). These parameters were compared against the u and D values derived from swimming behavior observations as independent validation of the models.
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Diffusion coefficients estimated from analyses of swimming tracks in 10 strains of bacteriaa
|
![]() View larger version (16K): [in a new window] |
FIG. 1. 2-D projections of swimming tracks of three different bacterial strains representative of the three motility types described in the text: nontumblers (strain HP15), rare tumblers (strain HP11), and frequent tumblers (strain HP46). The dots represent positions at 0.16-s (A and B) or 0.08-s (C) time intervals.
|
140° (
= -0.67). HP46 differed significantly from the rest, with a lower average turn angle (ca. 125°;
= -0.48). From the motility patterns we estimated the diffusion coefficients for the three motile groups (Table 2). Strains with long run lengths had diffusion coefficients on the order of 10-5 cm2 s-1, whereas strains that tumble frequently had diffusivities about an order of magnitude lower (ca. 10-6 cm2 s-1). No diffusion coefficient could be assigned the nontumbling species on the basis of equation 3.
Colonization pattern.
All experiments showed the same colonization pattern, i.e., a decelerating accumulation rate of bacteria (Fig. 2). We first interpreted this as a result of the non-steady-state condition in still water, but obviously equation 8b fits the observations very poorly since the accumulation of bacteria levels off much faster than this model can account for (Fig. 2A). In experiments where the aggregates were suspended in a flow we also found a saturating response (Fig. 2B). Because steady state is almost instantaneous in a flow regime, we had expected a linear increase here. Non-steady state, thus, cannot account for the saturating response in flow experiments either. Microscopic observations showed that bacteria that had attached to an aggregate might detach. We therefore conducted an experiment in which we incubated agar spheres for 60 min in a suspension of bacteria and then moved the spheres to sterile filtered seawater and followed the decrease of bacteria on the aggregate (Fig. 3). From the exponential decline we conclude that bacteria detach at a constant fraction per time; thus, encounter models need to incorporate the detachment rate (
in equation 6a and equation 8a; Fig. 2A).
![]() View larger version (23K): [in a new window] |
FIG. 2. Colonization of 0.2-cm radius agar spheres in still (A) and flowing (B) water by strain HP11 bacteria. Observations (solid symbols) are described by three different models: non-steady-state diffusion (equation 8b [solid line]), steady state with detachment (equation 6a [dotted line]), and non-steady-state diffusion with detachment (equation 8a [dashed line]). Error bars indicate the standard deviations of 10 measurements.
|
![]() View larger version (17K): [in a new window] |
FIG. 3. Colonization and detachment of bacteria (strain HP11) on 0.2-cm-radius agar spheres in still water. During the first 60 min, spheres were incubated at a bacterial concentration of 2.7 x 105 ml-1. At 60 min the spheres were transferred to sterile filtered seawater. Equation 8a was fitted to the data taken between 0 and 60 min, and an exponential model was fitted to the data taken from 60 min onward. Error bars indicate the standard deviations of 10 measurements.
|
![]() View larger version (29K): [in a new window] |
FIG. 4. Effect of ambient bacterial concentration on accumulation rate on 0.2-cm agar spheres of strain HP11 bacteria in still (A and B) and flowing water (0.33 cm s-1) (C and D). Panels A and C show absolute abundances of bacteria per aggregate, while panels B and D show abundances normalized with ambient concentration. Equation 8a has been fitted to the data (lines). Error bars indicate the standard deviations of 10 measurements.
|
![]() View larger version (24K): [in a new window] |
FIG. 5. Accumulation of different strains of bacteria on 0.2-cm agar spheres in still water. Equation 8a has been fitted to the data (lines). Error bars indicate the standard deviations of 10 measurements.
|
![]() View larger version (38K): [in a new window] |
FIG. 6. Average (± the standard error) diffusivities and detachment rates estimated from colonization rates of eight motile bacterial strains on 0.2-cm agar unenriched spheres in still water by fitting equation 8a to the observed accumulation of bacteria. The numbers within the columns indicate the numbers of experiments.
|
|
View this table: [in a new window] |
TABLE 3. Swimming velocities (u), diffusivities (D), and fractional detachment rates ( ) estimated from still-water colonization rates of 10 bacterial strainsa
|
Effect of aggregate size.
Our models predict that the total transport toward a spherical aggregate of bacteria with a random walk motility increases with aggregate radius, both in stationary and in sinking aggregates (cf. equations 2c and 7). At steady state, transport scales with aggregate radius (for constant Sh). At non-steady state, transport increases with radius to a power of 1 to 2 (equation 7). However, D and presumably
should be independent of aggregate size. These predictions are largely confirmed by observations (Fig. 7). Experiments with glass beads in still water followed the prediction most closely, with estimates of D and
varying by <10% for bead radii from 0.15 to 0.375 cm (Table 4). D · Sh for agar spheres in flow also varied very little, whereas experiments with agar spheres in still water yielded more variable estimates of D and
.
![]() View larger version (17K): [in a new window] |
FIG. 7. Effect of sphere size on the accumulation of strain HP11 bacteria in still and flowing water. Glass beads in still water (A), agar spheres in still water (B), and agar spheres in flowing water (0.27 cm s-1) (C). Equation 8a (still water) or equation 6a (flowing water) has been fitted to the data (lines). Error bars indicate the standard deviations of 10 measurements. Symbols: , small; , medium; , large.
|
|
View this table: [in a new window] |
TABLE 4. Estimates of D or D · Sh and (in parentheses) for experiments with various sphere sizesa
|
0.3 cm s-1) were ca. 5 (range, 1 to 9) and 20 times the estimated D in still water (Fig. 9), yielding estimated Sh values of 5 and 20, respectively. The expected Sh value is about 12 for both of these strains (equation 4).
![]() View larger version (14K): [in a new window] |
FIG. 8. Effect of water flow (0.1 or 0.3 cm s-1) on accumulation of strain HP11 (A) and strain HP25 (B) bacteria on 0.2-cm unenriched agar spheres. Equation 8a (still water) or equation 6a (flowing water) has been fitted to the data (lines). Error bars indicate the standard deviations of 10 measurements.
|
![]() View larger version (23K): [in a new window] |
FIG. 9. Average (± the standard error) estimates of diffusivity (D) or diffusivity multiplied by the Sherwood number (D · Sh) for strain HP11 and strain HP39 bacteria estimated from accumulation rates on 0.2-cm unenriched agar spheres by fitting equation 8a (still water) or equation 6a (flow) to the observed accumulation of bacteria.
|
1 µm (Table 1), the predicted encounter rate kernels are 0.5 and 1.5 x 10-8 ml s-1 for the low and high flow rates, respectively. The reason for this discrepancy is unclear.
Chemical signals (DMSP and MB).
The presence of marine broth in the agar substantially enhanced still-water colonization rates in the tumbling strains (HP33, HP39, and HP11), suggesting a chemosensory response of the bacteria (see Fig. 10A and B for examples). Among these three strains, only HP11 responded to the presence of DMSP in the agar (Fig. 10C). Nontumbling strain HP15 did not respond to the presence of MB in the agar (Fig. 10D)
![]() View larger version (25K): [in a new window] |
FIG. 10. Accumulation of strains HP11, HP39, and HP15 bacteria on enriched (open symbols, dotted lines) and unenriched (closed symbols, solid lines) 0.2-cm agar spheres in still water. Equation 8a has been fitted to the data (lines). In panels A, B, and D the agar spheres were enriched with 15 g of marine broth liter-1; in panel C the agar spheres were enriched with 1 mM DMSP. Error bars indicate the standard deviations of 10 measurements.
|
|
|
|---|
10-6 to 10-5 cm2 s-1as are the two estimates of swimming velocities for the nontumbling strains, u
10 to 30 µm s-1 (Tables 2 and 3). Also, the predicted effects of fluid flow (sinking), ambient bacterial concentration, and aggregate size on colonization are largely consistent with observations for the motile strains. However, the differences between strains in swimming behavior were only partly reflected in differences in colonization rates. For the group of bacterial strains with long tumbling intervals (HP1, HP5, and HP11), the correspondence between observed colonization rate and swimming behavior was good. However, for strains with short tumbling intervals (HP33, HP39, and HP46), D derived from colonization experiments (Table 3) was ca. 10-fold higher than that predicted from behavior observations and Gaussian random walk assumption (Table 2), i.e., behavioral observations would underestimate their colonization rates. One possible reason for this discrepancy is that Gaussian random walks may not appropriately describe frequently tumbling bacteria. Other types of random walks, Levy walks, provide a more efficient means of encountering food patches or particles and have been suggested as good descriptions of motilities of various organisms, including microbes (28, 42). Levy walks are characterized by bursts of tumbles interrupted by occasional very long runs, not unlike that seen for HP46 (Fig. 1).
The swimming velocities derived from colonization experiments for the nontumbling linear swimmers (strains HP4 and 15; Table 3) are less by a factor of 2 to 3 than those measured directly (Table 2). Thus, behavioral observations would overestimate the colonization rates. However, the efficiency of the linear swimming strategy is constrained by thermally driven Brownian rotation: the cells rotate continuously, causing constant and slight directional shifts (cf. Fig. 1A and B), and at sufficiently large spatial scales, rotational diffusion causes linear swimming to approach a diffusion process. The equivalent diffusion coefficients due to Brownian rotation for strains HP4 and HP15, computed from equations by Berg (8), are ca. 2 x 10-5 cm2 s-1, similar to those derived from the colonization experiments (Fig. 6). Therefore, Brownian rotation can account for the discrepancy between colonization measurements and behavioral observations.
Chemokinetic swimming behavior.
From the discussion above it follows that tumbling reduces the efficiency of encountering resource patches and particles. Why, then, tumble at all? The traditional explanation is that tumbling allows the bacteria to modify their swimming paths in a chemical gradient (8), increasing their chance of locating a food patch. In our experiments, enriching the agar spheres with attractant molecules (MB and DMSP) increased the initial colonization rate by a factor of 5 to 10 among the tumbling strains (Fig. 10A to C), whereas the nontumbling strain was not affected (Fig. 10D). Our observations therefore support the idea that tumbling is a prerequisite of the chemosensory behavior commonly found in marine pelagic bacteria (9, 15, 31).
How long does it take to find an aggregate?
The bacterial communities associated with aggregates are typically different from the bacterial communities in the surrounding water (13, 26, 32). This suggests that some bacteria specialize on colonizing aggregates and that their survival relies on how fast they encounter an aggregate. The rate at which a bacterium encounters aggregates depends on the size distribution and concentration of particles, as well as on the encounter rate kernel of individual bacteria, and can be estimated as
![]() | (9) |
![]() | (10) |
![]() | (11) |
![]() | (12) |
Can bacteria swim for long enough to encounter an aggregate? The cost of swimming for bacteria can be computed as the viscous drag (of a sphere) multiplied by the swimming speed (8): 6
ru2, where
is the dynamic viscosity (
10-2 g cm-1 s-1). Alldredge et al. (2) reported average volumes of 1 µm3 for aggregated-attached bacteria from various environments. For a 0.6-µm radius (1-µm3 volume) bacterium swimming at 100 µm s-1 the power requirement is 10-16 J s-1. A 1-µm3 bacterium contains ca. 10-13g of C (assuming 0.1 g of C cm-3), which by respiratory combustion may release 0.5 x 10-8 J (assuming 50 kJ per g of C). If the efficiency of the propulsion system is 1% and if other metabolic expenses are ignored, then a starving bacterium has energy resources enough for swimming continuously for almost 6 days. Obviously, the metabolic costs of swimming do not limit the chance that a bacterium may encounter an aggregate. This conclusion is independent of the actual swimming velocities of the bacteria because the search time decreases but the swimming cost increases with swimming velocity squared. It follows then that, in most upper-ocean habitats, aggregates are sufficiently abundant so that most bacteria can make it to an aggregate before they starve to death.
Population dynamics of attached bacteria.
Our observations suggest that bacteria that have colonized an aggregate may detach. The estimates of fractional detachment rate (
) are similar for all of the bacterial strains examined, i.e., ca. 10-4 s-1. Growth of attached bacteria may bias our
estimates low, but since conceivable magnitudes of specific growth rates are much less than our estimates of specific detachment rates (see below), this is a relatively small error. A detachment rate similar to that found here can be estimated for surface-attached Pseudomonas sp. from data in the study by Baty et al. (7). Such detachment rates imply that the bacteria have an average residence time of ca. 3 h on an aggregate. This is of the same order of magnitude as the residence time in the upper mixed layer of an aggregate sinking at 100 m day-1. Aggregates are risky environments for attached organisms because aggregates are eaten or they sink toward the seafloor. Therefore, to maintain a population in the upper ocean, the bacteria will have to leave the aggregate either by releasing progeny (6) or by detaching. The high detachment rates estimated here suggest that there is a considerable exchange of bacteria between aggregates and the surrounding water (cf. reference 16). The detachment of bacteria is probably only physically possible relatively soon after attachment or cell division because the bacteria may become embedded in the mucus film or matrix that may cover or be part of an aggregate (see reference 34).
The balance between colonization and detachment will lead to an equilibrium abundance of bacteria on an aggregate. From equations 5 and 11, the equilibrium abundance of bacteria on a sinking aggregate is calculated as:
![]() | (13) |
1.45 for aggregates with a 0.005- to 1-cm radius. However, this scaling is inconsistent with actual observations: bacterial abundances on natural aggregates scale with aggregate size to the power of 0.25 (Fig. 11; see also reference 4). Bacterial abundances normalized by ambient concentration of bacteria are also much higher on aggregates than predicted. This estimated difference is conservative, since not all ambient bacteria are potential colonizers. There are several possible explanations for these discrepancies.
![]() View larger version (16K): [in a new window] |
FIG. 11. Steady-state abundances of bacteria on aggregates as a function of aggregate size predicted from equation 13 (dashed line) compared to actually observed abundances on field-collected aggregates (symbols and solid line) (taken from the compilation of Kiørboe [22]). The different symbols refer to data from different studies. The predicted relation assumes a bacterial diffusivity of 10-5 cm2 s-1. Abundances have been normalized by ambient concentrations of bacteria. The predicted abundance increases with aggregate size raised to a power of 1.45 within the aggregate size range considered, while the log-log regression for the field-collected aggregates has a slope of 0.25.
|
Second, if bacteria grow on the aggregates at a constant rate (µ), the equilibrium abundance of bacteria would be F/(
- µ) rather than F/
, i.e., higher than predicted. Reported growth rates of attached bacteria are on the order of 1 day-1
10-5 s-1 (2, 35). Since estimated detachment rates are about 10 times higher, ca. 10-4 s-1, a growth rate of this order only increases the equilibrium abundance by 10%. Some reports suggest that attached bacteria grow much faster than free ones (see, for example, references 16 and 20). However, even a growth rate of 5 day-1 will only double the equilibrium abundance. Thus, growth alone cannot account for the different magnitudes of abundances. Moreover, growth per se does not change the scaling.
Another possible explanation is that the detachment rate declines with time as an increasing fraction of the bacteria become permanently attached. This would lead to higher equilibrium abundances of bacteria. However, it does not account for the difference in scaling, since this phenomenon will probably be more pronounced for the larger and presumably older aggregates.
Finally, the equilibrium bacterial abundances on aggregates may be governed by biological interactions other than colonization and detachment. Bacterial communities on natural aggregates consist of many species, not single-species assemblages as considered in our experiments. While the encounter rate between aggregates and free bacteria is independent of processes on the aggregate, the "decision" to attach or detach may be modified by complex inter- and intraspecific interactions among bacteria, e.g., through the release of signal molecules ("quorum sensing") and antibiotic substances. For example, swarming behavior in some pathogenic bacteria appears to be regulated through quorum sensing (see, for example, reference 14), and we have evidence that some of the strains isolated from marine particles indeed release signal molecules (17). Likewise, a large fraction of bacterial isolates from marine particles express antagonistic activity (29). Such interactions may change the population dynamics substantially. Marine aggregates also house a rich protozoan fauna that feeds on attached bacteria (10), and predator-prey interactions may render bacterial abundances on aggregates independent of aggregate size (22). Such interactions also depend on the densitynot only the abundanceof bacteria and, hence, on the fractal nature of natural aggregates. To more fully understand the population dynamics of bacteria on aggregates, processes of this kind need to be explored.
We thank Uffe H. Thygesen for help with model and software development and Andre Visser, George Jackson, and Josefin Titelman for fruitful discussions.
Present address: Max Planck Institute for Marine Microbiology, D-28359 Bremen, Germany. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»