Nutrient Limitation of Bacterioplankton Growth in Lake Dillon, Colorado
Nutrient Limitation of Bacterioplankton Growth in Lake Dillon, Colorado
Abstract
Bacterioplankton biomass, production, and growth rate were measured over a 2-yr period in
Lake Dillon, a mesotrophic Colorado reservoir. In addition, a multivariate statistical analysis and
nutrient addition experiments were used to analyze the regulation of bacterioplankton growth in
situ. Biomass ranged between 170 (winter) and 2,200 mg C m-* (summer); production ranged from
10 (winter) to 625 mg C m-* d-* (summer); annual bacterioplankton production was 47 g C m-2
yr-l in 1987 and 67 in 1988. Population growth rates ranged between 0.001 and 0.08 1 h I. Growth
rates in the epilimnion were substantially below estimated potential rates, suggesting severe nutrient
limitation during the period of stratification. Population growth rates were highly correlated with
P concentrations but not with dissolved organic C concentrations. Bacterioplankton growth in the
summer epilimnion responded strongly to the addition of P alone or in combination with N or
labile organic C. The field data show that bacterioplankton arc frequently without suflicient nutrients
to sustain maximum growth; the experimental and statistical analysts indicate that P, rather than
organic C, is the critical nutrient for bacterioplankton growth in this lake.
that will segregate the roles of C and inor- lyzed for particulate organic C (POC) with
ganic nutrients. an elemental analyzer. Particulate P (PP)
During 1987 and 1988, we studied bac- was analyzed according to the method of
terioplankton growth and nutrient limita- Solorzano and Sharp (1980).
tion in Lake Dillon, a mesotrophic Colo- Phytoplankton were collected on GF/C
rado mountain reservoir (elevation, 2,750 filters and chlorophyll a and pheophytin a
m; surface area, 1,330 ha; Morris and Lewis were measured by the spectrophotometric
1988; Lewis et al. 1984). Our purpose was method of Marker et al. (1980) and Nusch
to quantify the role of bacterioplankton in (1980), which is based on extraction of pig-
the C cycle of the lake and to determine the ment with hot 90% ethanol (vol : vol) and
relative importance of C, P, and N in reg- subsequent acidification. Results were pe-
ulating the growth of bacteria. The study riodically validated by HPLC.
relies on data from field populations, nu- Samples collected for determination of
trient limitation experiments, and multi- dissolved organic C (DOC) were filtered
variate statistical analysis to provide insight gently (< 130 mm of Hg) through ashed Gel-
into the seasonal dynamics of bacterio- man A/E filters. The filtrate was frozen
plankton in the lake. pending analysis which was usually within
2 weeks. The sample was digested with per-
Methods sulfate and phosphoric acid in sealed am-
Analytical methods - Samples were col- pules (Wetzel and Likens 1979); CO, from
lected from Lake Dillon at the point of max- the resulting decomposition of DOC was
imum depth (60 m) at biweekly intervals quantified by gas chromatography.
from May through September 1987 and Primary production was determined in
1988 and at monthly intervals during the situ by the 14C method described by Wetzel
rest of the year. On each date, samples were and Likens (1979). Samples collected from
obtained near the surface (0.5, 1.O, 2.0, and seven depths in the photic zone were spiked
5.0 m) with a 7-liter Van Dorn bottle; ad- with 74 kBq of NaH14C0, per sample and
ditional samples were taken with a 5-m in- suspended in 130-ml Pyrex bottles at the
tegrating sampler at 5-m increments from depth from which they were collected. Sam-
7.5 m to the bottom. The samplers were ple incubation (2-4 h) was usually con-
soaked in ethanol and thoroughly rinsed with ducted near midday. Volumetric productiv-
lake water before use. ity values were converted to areal estimates
Samples collected for chemical analysis as suggested by Wetzel and Likens.
were kept cold and dark until they were an- Bacterioplankton were preserved with
alyzed. Particulate and dissolved fractions buffered formaldehyde or glutaraldehyde
were separated by vacuum filtration through (2% final concn), filtered onto 0.2~pm Nu-
ashed and tared Whatman GFK filters. Sol- clepore filters, and counted by fluorescence
uble reactive P (SRP) in the filtrate was an- microscopy at 1,250 x (Hobbie et al. 1977).
alyzed by the ascorbic acid-molybdate At least 400 bacteria were counted per sam-
method of Murphy and Riley (1962). Total ple; no fewer than 10 fields were examined
dissolved P (TDP) was determined by a per filter (expected counting error <5%). It
modification of the oxidation method de- is often difficult to distinguish bacteria and
scribed by Lagler and Hendrix (1982) and photosynthetic picoplankton with acridine
Valderrama (198 1). Dissolved organic P orange. For this reason, unstained samples
(DOP) was calculated as the difference be- were occasionally examined by epifluores-
tween TDP and SRP. NO3 was determined cent microscopy, which is capable of show-
by liquid anion chromatography. ing chlorophyll autofluorescence. Less than
Suspended particulate matter was col- 5% of cells of < l-pm diameter exhibited
lected onto ashed and tared GFK filters, chlorophyll autofluorescence.
dried to constant weight at 6O”C, and re- The mean cell volume of bacterioplank-
weighed. Subsequent to determination of ton was determined for the mixed layer on
total particulate matter, filters were ana- all dates and for the hypolimnion when the
Bacteria nutrient limitation 1181
lake was stratified. Composite samples were external standard method; blanks consisted
created by abundance-weighted mixing of of samples that were fixed with formalde-
samples from different depths. Cell volumes hyde (2% final concn) before adding the la-
were determined by fluorescence micros- bel.
copy and photomicrography (Fry 1988). Conversion factors for Tdr incorporation
Cells were collected by filtration and stained were determined in triplicate for each sam-
with acridine orange. Fields were photo- pling date from the ratio of Tdr incorpo-
graphed with Kodak T-Max (ASA 2000) or ration to growth rate in unenriched cultures
Kodak Tri-X Pan (ASA 400) high-contrast of bacterioplankton that were collected from
black-and-white film. Negatives were pro- the lake and diluted by the methods given
jected and cell dimensions (100 cells per below. New cells were assumed to have a
sample) were measured with a micrometer. volume equivalent to the mean cell volume
The measurements were routinely calibrat- of bacteria in the stratum of origin on the
ed from photographs of a stage micrometer sampling date.
or standardized microspheres that were pro- Several variables have been converted
jected along with the images of bacterio- from volumetric to depth-integrated (areal)
plankton. The cell volumes were calculated values. For these calculations, each sample
as spheres (cocci) or straight-sided rods with was weighted according to the proportion
hemispheric ends (bacilli). These model of total lake volume found in the stratum
shapes were appropriate because cells were from which it came.
predominantly cocci, bacilli, or shallow vib- Nutrient response experiments-Nutrient
rio. Few deep vibrio or spirillum shapes were limitation was assessed by the growth re-
present in the plankton. Cell volumes were sponse of dilute bacterioplankton commu-
converted to C biomass with the relation- nities to additions of various organic and
ship of Simon and Azam (1989). On the inorganic nutrients. Water collected from
basis of this relationship, the average bac- the epilimnion was split into two fractions
terioplankton in the epilimnion of the lake in’ the laboratory. The first fraction was
(mean cell volume of 0.144 pm”) contained gently filtered (< 130 mm of Hg) through
0.207 pg C pm-3. ashed Whatman GF/D and Gelman A/E
Bacterioplankton production was deter- glass-fiber filters. Although the effective pore
mined from the rate of incorporation of thy- size of the A/E filter is larger than the small-
midine (Tdr) into DNA in situ. Triplicate est bacteria, counts on the filtrate indicated
15-ml samples of lake water were placed that > 98% of the cells were removed by this
into sterile polypropylene test tubes and treatment. Analysis of the filtrate also
spiked with [3H-methyllthymidine (20 nM showed that filtration caused no statistically
Tdr; 1,850-2,590 GBq mmol-‘). The sam- significant increases in dissolved P, N, or
ples were suspended in the lake at the depth DOC:
from which they were collected and allowed The filtrate was used to dilute the natural
to incubate for 1 h (summer) to 4 h (winter). plankton community (1 : 5 to 1 : 20). Dilu-
Uptake of label was stopped by adding tion was greater during summer stratifica-
formaldehyde (2% final concn). Samples tion when bacteria, phytoplankton, and
were then placed on ice and transported to bacteriovores were most numerous at this
the laboratory where they were filtered onto time. The purpose of dilution was threefold:
0.2-pm-pore nylon membrane filters (MS1 to reduce the grazing pressure on bacteria-
Magna) and thoroughly rinsed with 20 ml plankton, to reduce the demand for soluble
of filtered lake water. Filters were frozen nutrients, and to minimize the effect of in-
until the DNA was extracted and purified creased phytoplankton C fixation and exu-
(usually <24 h). DNA was extracted and dation in response to nutrient enrichment.
purified by a modification of the method of During 1987, each enrichment experi-
Findlay et al. (1984), which is based on the ment consisted of three replicates each of a
solubility of DNA in hot and cold TCA. control and seven nutrient addition treat-
Counting efficiency was estimated by the ments. The control was diluted lake water
1182 Morris and Lewis
- 40
5w
I
h. - 30
p o!
E
m
-20 o9
f”
b
m 200 H
5 - 10
100
JFMAMJJASONDJFMAMJJASOND
0
JFMAMJJASONDJFMAMJJASOND 1987 Date 1988
1987 Date 1988 Fig. 2. Areal estimates for Chl a, pheophytin a, and
Fig. 1. Areal estimates of primary and bacterial POC for Lake Dillon, 1987-l 988. Hatching-ice cov-
production for Lake Dillon, 1987-l 988. Hatching- er.
ice cover. -i
for the epilimnion, hypolimnion, and the
78 g C mm2 yr-l in 1987 and 101 in 1988. entire water column are shown in Table 2.
In both years, rates peaked in fall as the The volumes were not normally distributed
mixed layer began to thicken but before (Kolmogorov-Smirnov goodness-of-fit test
mixing was complete (Fig. 1). Chl a and for normality, P I 0.001; Sokal and Rohlf
POC peaked about the same time (Fig. 2). 198 1). Cell volumes for the three strata were
Similar primary production trends have contrasted by use of the Kruskal-Wallis test
been noted previously for the lake (Lewis et followed by a nonparametric multiple com-
al. 1984; Morris and Lewis 1988) and can parisons analysis (Devore 1982). The bac-
be explained by the relief of nutrient limi- terioplankton of the hypolimnion were
tation that develops during summer strati- smallest (mean, 0.069 pm3); bacterioplank-
fication. Pheophytin a, which is an index of ton from the whole water column during
decomposition or of zooplankton grazing winter and during spring and fall mixing
(Carpenter and Bergquist 198 5), showed were about the same size (mean, 0.075 pm3),
only slight seasonal variation; maximum although the small difference in size was
values occurred in fall and under the ice statistically significant at a I 0.05. Bacte-
(Fig. 2). rioplankton of the epilimnion were larger
The distribution of bacterioplankton in (mean, 0.144 pm3; cy I 0.05).
the lake showed extensive seasonal and ver- In both years, mean cell volume was low-
tical variation (Fig. 3A). Numbers of bac- est in winter and increased steadily in the
teria were lowest during winter under ice epilimnion during summer stratification.
and during spring and late fall mixing (0. l- Maximum cell volumes occurred near the
0.3 x lo6 cells ml-l); bacteria were distrib- peak of water-column stability and de-
uted uniformly throughout the water col- creased again as the thermocline thickened
umn. The development of larger numbers in late summer and fall. The nonparametric
of bacterioplankton generally corresponded Mann-Whitney U-test (Sokal and Rohlf
to the increasing stability of the epilimnion; 198 1) showed that cells in the epilimnion
in both 1987 and 198 8 bacterioplankton were significantly larger than those in the
abundance peaked in the epilimnion in late hypolimnion (P 5 0.001) for all 16 dates
summer (0.6-1.8 x lo6 cells ml-l). Except when the lake was stratified.
for occasional minor peaks just below the Bacterioplankton biomass showed marked
mixed layer, bacterioplankton density de- seasonal and vertical variation in 1987 and
creased precipitously between the epilim- 1988 (Fig. 3B). A large proportion of the
nion and the hypolimnion. bacterioplankton biomass was restricted to
Pooled bacterioplankton cell volume data the epilimnion. In 1987, the biomass peak
1184 Morris and Lewis
20
50
JFMAMJJASONDJFMAMJJASOND JFMAMJJASONDJFMAMJJASOND
IMI S IM-
40
50
i
JFMAMJJASONDJFMAMJJASOND JFMAMJJASONDJFMAMJJASOND
1987 1988 1987 1988
Fig. 3. Temporal and vertical variation for (A) bacterioplankton abundance (lo6 cells ml-l), (B) bacterio-
plankton biomass (mg C m-)), (C) tritiated Tdr incorporation into DNA (pmol liter-l h-l), and (D) bacterio-
plankton production (mg C m-3 d-l) for Lake Dillon, 1987-1988. I-Ice cover; M-mixing; S-stratification.
(33 mg C me3) was lower than in 1988 (75 showed little variation. Biomass per unit
mg C md3). Biomass was often uniformly area (Fig. 4) was lowest under the ice or
distributed within the epilimnion. On some during spring mixing (170-490 mg C mW2)
dates, biomass peaked at or near the ther- and highest during summer stratification
mocline. Bacterioplankton biomass in the (1 ,OOO-2,200 mg C m-‘).
hypolimnion was substantially lower than Incorporation of tritiated Tdr also showed
in the epilimnion (l-8 mg C m-3) and marked seasonal and vertical variation (Fig.
Table 2. Pooled bacterioplankton cell volume data for Lake Dillon, 1987-1988.
Min Max
Mean volume
Sample type N Cm’, SE) Cm’)
Ol’,““““““““““‘I
JFMAMJJASONDJFMAMJJASOND”
JFMAMJJASONDJFMAMJJASOND
1987 Date 1988
1987 Date 1988
Fig. 5. Estimates of population growth rates of bac-
Fig. 4. Areal estimates of bacterioplankton bio- terioplankton in surface and bottom water of Lake Dil-
mass and the ratio ofbacterioplankton biomass to POC lon, 1987-1988. Hatching-ice cover.
in Lake Dillon, 1987-1988. Hatching-ice cover.
3C). The highest incorporation rates were 52 1% of algal production. The highest ratios
observed in the summer epilimnion (1 O-l 2 consistently appeared during the winter un-
pmol liter-’ h-l); the lowest rates occurred der ice when algal and bacterial production
in the hypolimnion or throughout the water were both low.
column in winter or during spring and fall Bacterioplankton production also ex-
mixing (0.2-2.0 pmol liter-’ h-l). Conver- ceeded primary production on several dates
sion factors relating Tdr incorporation to during summer stratification. Annual areal
growth rate were estimated for 13 dates; on bacterioplankton production was 60% of
other dates (mostly winter), growth was too phytoplankton production in 1987 and 66%
slow to allow a reliable estimate. Because in 1988.
conversion factors did not follow any sea- Cell densities and production (as cells
sonal trend, a mean value was used for all produced per unit time) were used as an
calculations (2.03 x lo9 cells nmol-1 Tdr; index of the population growth rate of bac-
SE 0.25 x 109). terioplankton in the lake. The mean pop-
Bacterioplankton production ranged be- ulation growth rate (h) for bacteria was 0.0 12
tween 0.1 and 23.0 mg C m-3 d-’ (Fig. 3D). h-’ (SE 0.0005; n = 409), which is equiv-
Rates of production were lowest in winter alent to a generation time of 2.4 d. Values
under ice and in the hypolimnion (0.1-2.5 of p ranged from 0.001 to 0.08 1 h-l (gen-
mg C rnd3 d-l); the highest rates were found eration time 29.0-0.4 d). In both 1987 and
in the epilimnion during late summer (8.0- 1988, maximum values of p (0.054 and
23.0 mg C m-3 d-l). Within the epilimnion, 0.08 1 h-l) occurred in the epilimnion dur-
production often increased with depth and ing the early stages of stratification (Fig. 5).
peaked at or near the thermocline, but de- Except for these notable peaks, p in the epi-
creased markedly in the hypolimnion. This limnion was generally low (0.005-0.0 10 h-l)
trend resulted from a combination of de- and often even less than p in the hypolim-
creased activity (Tdr incorporation) and a nion.
smaller mean cell volume of bacterioplank-
ton in the hypolimnion. Areal rates of bac- Discussion
terioplankton production (Fig. 1) were low- Biomass, production, and growth rates-
est in winter (lo-50 mg C m-2 d-l) and During 198 7, mean cell volume increased
peaked during summer stratification (270- by a factor of 2.2 from winter minima to
625 mg C m-2 d-l). Annual areal bacterio- summer maxima, while cell abundance in-
plankton production was 47 g C mm2 yr-1 creased by a factor of 2.5. Increases in mean
in 1987 and 67 in 1988. Expressed on an cell volume alone accounted for 47% of the
areal basis, daily estimates of bacterio- seasonal variation of bacteria biomass dur-
plankton production ranged between 22 and ing 1987 and 35% in 1988.
1186 Morris and Lewis
_ 0.070
‘;
G 0.060
s!
$ 0.050
fi 0.040
8 0.030
0.020
0.010
O.OW
JFMAMJJASONDJFMAMJJASOND
Table 4. Results of the nutrient enrichment experiments for Lake Dillon. Treatments with significantly higher
Tdr production rates that control (P I 0.05): L-low response, < 100% increase relative to control; M-moderate
response, lOO-500% increase relative to control; H-high response, > 500% increase relative to control. Amino
acid (AA) and yeast extract (YE) additions performed only in 1988.
Limitation
16 Mar 87 ns ns ns ns ns ns ns None
6 Apr 87 ns ns L L ns L L P
10 Jun 87 ns ns H H ns H H P
24 Jun 87 ns ns H H ns H H P -
7 Jul87 ns L H H L H H P/N -
22 Jul 87 ns ns H H ns H H P -
19Aug87 ns ns M H ns M H P C
2 Sep 87 ns ns ns ns ns ns ns None
16 Sep 87 ns ns H H ns H H P
30 Sep 87 ns ns H H L H H P -
26Oct87 ns ns ns ns ns ns ns None
23 Nov87 ns ns ns ns ns ns ns None
22 Fcb 88 ns ns ns ns ns ns ns ns ns None
14 Mar 88 ns ns ns ns ns ns ns ns ns None
6 Apr 88 ns ns ns ns ns ns ns ns None
23 May 88 ns ns M iti ns M M ns M P -
8 Jun 88 ns ns M M ns M M ns M P -
22 Jun 88 ns ns M M ns M M ns M P -
6 Jul 88 ns ns ns ns ns ns M ns M P+N+C -
17Jul88 ns ns H H ns H H ns H P C
3 Aug 88 ns ns M H ns M H ns H P C
17 Aug 88 ns ns ns ns ns M M ns H P+N -
31 Aug 88 ns ns L L ns L L ns M P N/C
15 Scp 88 ns L L M L M L L H P/N -
27 Sep 88 ns L M M L M H ns H P/N C
25 Oct88 ns ns L L ns L L ns L P -
22Nov88 ns ns ns L ns ns L ns L p+c -
nities to additions of amino acids is ambig- sist of a single nutrient or several nutrients
uous because amino acids provide both C operating concurrently or reciprocally.
and N, but can be interpreted from results The nutrient limitation that was suggest-
of separate and combined additions of C ed by low population growth rates and by
and N. Yeast extract is a multiple nutrient multivariate statistical analysis was con-
source containing a variety of C, N, and P firmed by the nutrient addition experiments
compounds as well as micronutrients. Ad- (Table 4). Bacterioplankton responded to
dition of yeast extract will result in growth nutrient enrichment during most of the
response whenever communities are limit- open-water period; the response was es-
ed by N, P, or C, but could also produce a pecially strong for bacteria in the summer
response in the absence of limitation by N, epilimnion (Fig. 8). P, either alone or in
P, or C if bacterioplankton are limited by combination with N or C, was implicated
trace elements or specific organic sub- in every case of nutrient limitation for the
stances. lake during 1987 and 1988.
Secondary limiting nutrients are those that In 1987, moderate P limitation appeared
are exhausted in cultures that have been en- in April, became more severe in summer,
riched with the primary limiting nutrient. and remained severe until fall turnover.
Their ecological significance is conditional During this period, two of eight experiments
on the availability of the primary limiting showed a strong secondary response to C
nutrient. A secondary limitation may con- or N.
1190 Morris and Lewis