Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2004, p. 6670-6677, Vol. 70, No. 11
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.11.6670-6677.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, Michigan
Received 12 March 2004/ Accepted 23 June 2004
|
|
|---|
|
|
|---|
One readily quantifiable feature of bacteria that has been proposed as a component of bacterial life histories is the number of rRNA genes (25, 26, 39). This proposal is based on the understanding that multiple copies of rRNA genes allow for increased rates of rRNA synthesis, leading to more rapid synthesis of ribosomes and ultimately conferring the potential for a quicker response to an influx of resources and rapid growth (13). This characteristic fits the classical definition of a component of life history and may also exhibit a tradeoff, namely, a metabolic cost for retention of multiple copies of rRNA genes associated with basal levels of transcription of these genes. This proposed expense would be particularly important in environments characterized by slow and constant flux of resources, where there is little selection for shorter lag times or rapid growth.
We have investigated the potential selective forces behind the occurrence of multiple rRNA gene copies in bacteria, and we report the results of experimental tests of the hypothesis that multiple copies of rRNA genes are a central component of life history and confer an advantage in environments defined by fluctuations in resource availability. Previous work addressed this hypothesis by illustrating the connection between the ecological strategies and rRNA gene copy numbers of a phylogenetically diverse collection of bacteria (25). In the present study, the number of rRNA genes was altered experimentally in an evolved lineage of Escherichia coli so as to investigate any direct link between the life history of an organism and its rRNA gene copy number.
|
|
|---|
For all experiments, cultures were first grown overnight in LB inoculated from 80°C glycerol stocks: these cultures were then used to inoculate a DM conditioning culture, which was incubated with shaking at 225 rpm and 37°C for exactly 24 h. The inoculum for the conditioning cultures represented approximately 1/100 of the population density achieved at the specific glucose concentration of the test medium to be used. This conditioning culture was then used to inoculate fresh DM (1/100) for the experiment performed.
Batch cultures for competition experiments, conditioning cultures, and cultures for the measurement of individual growth parameters were grown in 10 ml of DM in a 50-ml Erlenmeyer flask or in 50 ml of DM in a 250-ml Erlenmeyer flask shaken at 225 rpm at 37°C. Competition experiments in chemostats were carried out in a system of vessels designed like those described previously (9, 11). The chemostat vessels were incubated in a 37°C water bath, and sterile DM + 25 was delivered to each vessel by a variable-speed peristaltic pump to achieve a dilution rate of 0.11 or 0.81 per h. A vacuum was applied to each culture vessel through a glass tube that maintained the culture volume at approximately 65 ml. The vacuum was also used to draw air through a
0.30-µm-pore-size HEPA filter (HEPA-VENT; Whatman Inc., Clifton, N.J.) attached to each medium inflow tube, providing sterile aeration and mixing.
General cloning procedures.
Restriction endonuclease digestions were carried out according to the manufacturer's protocols (New England Biolabs, Beverly, Mass.). If necessary, specific restriction fragments were isolated on an agarose gel and purified (QIAEX II gel extraction kit; QIAGEN, Valencia, Calif.). A commercially available, chemically competent strain of E. coli (One Shot TOP10; Invitrogen, Carlsbad, Calif.) was used for most cloning according to the manufacturer's protocols. When E. coli strain D308 was transformed, cells were first made chemically competent by using a modification of the protocol of Chung et al. (12, 39).
PCR conditions.
DNA fragments of <4 kb were amplified by using 0.25 U of Taq DNA polymerase (Biolase; Intermountain Scientific Corporation, Kaysville, Utah) with 1x reaction buffer, 2 mM MgCl2, 2.5 mM each deoxynucleoside triphosphate, 50 to 100 ng of template DNA, and 0.5 µM forward and reverse primers (Table 1) in a total volume of 50 µl. Reactions were carried out in a PTC100 thermal cycler (MJ Research Inc., South San Francisco, Calif.) with the following incubation conditions: 94°C for 3 min; 30 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 45 s; and 72°C for 10 min. When larger DNA fragments (>4 kb) were to be amplified, the proofreading DNA polymerase (rTth) was used (GeneAmp XL PCR kit; Perkin-Elmer Applied Biosystems, Foster City, Calif.). Reactions were set up according to the manufacturer's protocols except that colony material was sometimes used as the template. The reactions were carried out under the following conditions: 94°C for 3 min; 30 cycles of 94°C for 30 s, 60°C for 10 min, and 72°C for 12 min; and 72°C for 10 min.
|
View this table: [in a new window] |
TABLE 1. PCR primers used in this study
|
![]() View larger version (20K): [in a new window] |
FIG. 1. Outline of procedure used to delete the rRNA operon rrnA from the chromosome of E. coli. Homologous recombination events are depicted as crossed dotted lines, which encompass regions of homology; plasmid-derived regions of DNA are depicted in grey.
|
rrnB strains (BS200; BS201) by using generalized transduction with P1vir (35). Transductants (Kanr Sucs) with the marker allele in place of the original rRNA operon were referred to as "intermediate" strains (BS110 and BS111, marker replacing rrnA; BS210 and BS211, marker allele replacing rrnB; BS310 and BS311, marker allele replacing rrnA in
rrnB strains). Intermediate strains were transformed at the permissive temperature (30°C) with pKTS102 (Fig. 1) or pKTS203, both of which were derived from pMAK705 [rep(Ts) Camr] (18), and contained a truncated rrnA (UpA:DnA) or rrnB (UpB:DnB) operon. Growth of transformants (Kanr Sucs Camr) at the nonpermissive temperature (44°C) promoted the integration of the plasmid into the chromosome through recombination (Fig. 1). Subsequent growth in LB without antibiotics at 44°C favored cells that had undergone a second recombination event (excision of the plasmid from the chromosome) and lost the plasmid (8). An allelic exchange occurred if this second recombination was within the region of homology opposite from the site of integration (as shown in Fig. 1). The resulting rRNA operon deletion strains (Kans Sucr Cams), which had undergone this allelic exchange and lost the plasmid containing the marker allele, were selected by their growth on LBsuc at 30°C.
The genotypes of the parental and deletion strains were verified by size analysis of PCR products from affected genomic regions and Southern hybridization of PvuII-digested genomic DNA with a digoxigenin-dUTP-labeled DNA probe for a conserved region of the 16S gene (positions 8 to 536) (as described in references 20, 25, and 39).
Relative fitness assays.
For batch culture competition experiments, relative fitness was determined as reported in reference 20 by a method described by Lenski et al. (28), by the following equations:
![]() | (1) |
![]() | (2) |
The relative fitness of each deletion strain was measured in chemostat culture competitions against the parental strain by monitoring the change in density of each population as a function of total population doublings. The fitness of each strain under chemostat conditions was also assessed by comparison of calculated J parameters. The J parameter, a weighted value of Ks (see below), represents the subsistence or "break-even" concentration of the limiting resource for each strain under chemostat conditions at equilibrium, and the steady-state concentration of the limiting resource when each strain is grown alone. The J parameter was calculated by the following equation:
![]() | (3) |
Measurement of individual components of fitness.
The time required for each strain to reach exponential growth after inoculation (lag time) was estimated by subtraction. The cultures used for these measurements were treated identically to those in batch culture competition experiments, except that each inoculation was staggered by 15-min intervals to ensure precise timing. Immediately after inoculation and mixing (t = 0), and exactly 4 h later (t = 4), population density was estimated by dilution and plating onto LB agar plates. The time interval of 4 h was determined empirically to be within mid-exponential growth for all strains.
The theoretical time required to reach the population density at t = 4 from the initial population density, at the maximal specific growth rate determined for each strain (described below), was subtracted from the actual time required (4 h). The µmax of each strain was estimated by calculating the slope of the linear regression of ln-transformed population density values (optical densities at 600 nm) versus time during exponential growth in DM+1000. DM + 1000 was used to obtain sufficient optical densities for accurate measurement with a Lambda 3 UV-Vis spectrophotometer (Perkin-Elmer Applied Biosystems).
Ks represents the concentration of a growth-limiting substrate at which a culture grows at half the maximal specific growth rate (30) and is defined by the following equation:
![]() | (4) |
The death rate of each strain during stationary phase was represented by the rate of decrease in population density between 12 and 24 h, measured as described by Vasi et al. (42). Population densities of three replicate cultures were estimated every hour after inoculation between 12 and 24 h by dilution and plating onto LB agar. The mean CFU per liter for each culture at each time point was ln transformed and plotted versus time. This ln-transformed regression was assumed to be linear, and the inverse slope was recorded as the death rate (i.e., a positive death rate represented a decrease in population density). Population densities continued to be monitored periodically for each strain up to 197 h after inoculation.
Data analyses.
Analysis of variance was used to test for significant differences in the means (P < 0.05) for relative fitness in batch culture (n = 24), lag time (n = 8), µmax (n = 12), and Ks (n = 3), followed by planned comparisons of the sum of squares to determine which means were significantly different from each other (38).
Modeling competition experiments.
A model designed to simulate the competition of two populations for a single limited resource was programmed using object-oriented software (STELLA, version 5.1.1; High Performance Systems, Hanover, N.H.) (39). The model is a three-compartment system that defines the interaction between the state variables of a limiting resource (glucose, expressed in grams per liter) and two competing populations (population 1 and population 2, expressed in cells per liter), in which the concentration of the limited resource in the culture vessel is the central state variable of the system. The dynamics of the limited resource are determined by the initial concentration (in grams per liter) and the utilization of the resource by each population (in grams of the resource per cell). The growth (cells per liter per hour) divided by the yield (cells per gram of resource) of each population determined the rate of removal of the resource from the system due to utilization (grams of the resource per liter per hour) by each population. The model was programmed to also allow for continuous or semicontinuous addition of resources and washout of cells to simulate chemostat culture competitions.
The dynamics of each population's density are controlled by growth and death. The rate of growth of each population follows the Monod equation, which defines the relationship between the specific growth rate of a population and to the concentration of a limiting resource (30). In addition to the parameters described above, lag time is programmed to delay the growth of a population at the beginning of a simulation, for the duration of the lag time (in hours). Since no significant death rate was measured during the 24-h competition experiments, this value was set to zero. All simulations were run by using Euler's method of derivation with a step time of 0.00625 h (Stella Technical Documentation; High Performance Systems, Inc.).
|
|
|---|
![]() View larger version (80K): [in a new window] |
FIG. 2. (A) Southern hybridization of PvuII-digested genomic DNA from parental, rrnB, rrnA, and rrnAB strains. Locations of hybridized bands corresponding to the 16S rRNA gene from each rRNA operon (bands A through E, G, and H) are indicated for the parental strain. (B) PCR-amplified rrnA and rrnB regions from each strain. PCR products corresponding to intact rrnA and rrnB operons are approximately 7.4 and 6.7 kb long, whereas deleted rrnA and rrnB operons are much smaller (2.6 and 1.9 kb, respectively). Sizes of markers in lanes M are given (in kilobases) to the right of each panel.
|
rrnA or
rrnB) decreased relative to that of the parental strain by 3 to 4% (P < 0.025), while strains with two rRNA operons deleted experienced a 9% decrease in fitness (P < 0.025) (Fig. 3; Table 2). The mean relative fitness values did not differ between the
rrnA and
rrnB strains, in which a single rRNA operon was deleted (P > 0.10). In reciprocal experiments, the arabinose marker proved to be neutral, and therefore Ara+ and Ara strains with the same rRNA operon genotype were treated as identical.
![]() View larger version (13K): [in a new window] |
FIG. 3. Effects of rRNA operon deletions on relative fitness (W) in batch culture. The mean relative fitness of the parental strain (Ara+ versus Ara) and of rRNA operon deletion strains ( rrnB, rrnA, and rrnAB) versus the parental strain was measured in direct (grey) and simulated (white) batch culture competition experiments. Error bars, standard deviation of the mean (n = 8).
|
|
View this table: [in a new window] |
TABLE 2. Fitness for each strain in batch and chemostat competition experiments
|
rrnB strains had very little effect on fitness in chemostats, whereas deletion of rrnA in
rrnA and
rrnAB strains resulted in a significant decrease in fitness relative to that of the parental strain (Fig. 4). Compared to the parental strain, the relatively small differences in J for
rrnB strains (
8%) were consistent with their coexistence with the parental strain in chemostat competition experiments over 14 doublings. The substantially larger differences in J for
rrnA and
rrnAB strains (
24%) compared to the parental strain accurately predicted the outcome of those chemostat competition experiments, in which
rrnA and
rrnAB strains were undetectable after as few as 8 doublings (Table 2; Fig. 4).
![]() View larger version (14K): [in a new window] |
FIG. 4. Effects of rRNA operon deletions on relative fitness in chemostat culture (dilution rate, 0.11 h1). The ratios of population densities (rrn deletion strain/parental strain) for rrnA (squares), rrnB (circles), rrnAB (triangles), and parental strains (filled diamonds) (Ara+ versus Ara) are plotted as a function of the number of doublings. Error bars, standard errors for four replicate competition cultures.
|
rrnAB strains, which was about 20% (12 min) longer, was significantly different from that of the parental strain (P < 0.05). The parental strains, therefore, had a head start of about one-third of a doubling over cells with only five rRNA operons in competition experiments. Maximal growth rates of the rRNA deletion strains were slightly lower than that of the parental strain and also correlated with the number of deleted rRNA operons (P < 0.10 for six rrn operons; P < 0.05 for five rrn operons). |
View this table: [in a new window] |
TABLE 3. Growth parameters for each straina
|
![]() View larger version (17K): [in a new window] |
FIG. 5. The decreases in relative growth rates (rrn deletion strain/parental strain) for rrnA (squares), rrnB (circles), and rrnAB (triangles) strains are plotted as a function of the growth rate of the parental strain on the same medium (DM+Ac, DM+Glu, LB+Glu). Data points marked with an asterisk represent growth rates significantly different from that of the parental strain (P < 0.05).
|
rrnA and
rrnAB strains correlated with the deletion of the rrnA operon (P < 0.05). None of the strains had an appreciable death rate between 12 and 24 h after inoculation, with no significant change in population densities until after 48 h, when the number of viable cells dropped rapidly in all strains.
Model simulations of competition experiments.
Competition experiment simulations were conducted using the individual parameters (lag time, µmax, Ks, and death rate) measured for each strain. Empirically determined values for initial population densities (2.0 x 108 cells liter1) and cell yields (1012 cells g1) were also used to parameterize the model. Relative fitness values determined from simulations of batch and chemostat competition experiments were very close to empirical values. Simulated batch culture competitions predicted relative fitness values (
rrnA, 0.958;
rrnB, 0.961;
rrnAB, 0.902) that were within 1 to 2% of empirical values (Fig. 3; Table 2). In chemostat competition simulations, the population densities of deletion strains followed dynamics similar to those indicated by empirical results, as shown by the values for J in Table 2. Furthermore, the chemostat simulations with data corroborated the mean Ks values measured for
rrnA and
rrnAB strains (Table 3).
|
|
|---|
Correlative data exists which links the number of rRNA operons on a bacterial chromosome with an organism's capacity for response to resource availability. Strains of the genus Mycobacterium, for example, can be roughly divided between slow growers (doubling time,
16 h) and fast-growers (doubling time, 3 to 4 h), having either one or two rRNA operons per chromosome, respectively (7, 34). There are also numerous examples of bacteria for which the number of rRNA operons on the chromosome coincides with their presumed life history. Bacteria such as the bioluminescent symbiont of the Caribbean flashlight fish Kryptophanaron alfredi (43), members of the aphid endosymbiotic genus Buchnera (37, 41), parasitic mycoplasmas (1), and the marine oligotroph Sphingomonas sp. strain RB2256 (17) have one or two rRNA operons. These bacteria exist in relatively stable environments, tend to have slow growth rates, and respond slowly to fluctuations in resource availability. Bacteria with a relatively high number of rRNA operons on their chromosomes include E. coli with 7 (5), Vibrio spp. with 7 to 10 (21, 24, 33), and Bacillus subtillus with 10 rRNA operons (27); these organisms have relatively fast growth rates and would be expected to respond more quickly to fluctuations in resource availability. The hypothesis that the number of rRNA operons is reflective of an organism's ecological strategy is also supported by the correlation between rRNA operon copy number and the rate of colony appearance for soil isolates (25). Bacterial isolates that formed visible colonies in a few days had an average of five rRNA operons, whereas isolates that formed colonies up to 2 weeks later had only one or two rRNA operons on their chromosomes.
Despite the considerable predictive power associated with having a readily assessed genomic marker (rRNA operon copy number) with which to infer an organism's ecological strategy, there are few direct measurements of the influence of rRNA operon copy number on an organism's life history. Previous alterations of rRNA operons in E. coli included strains in which plasmid-borne rRNA operons were introduced (as in references 23 and 40), where rRNA genes were truncated or interrupted with antibiotic resistance genes but could still be expressed from intact promoters (2-4, 13), or in which adjacent genes on the chromosome were affected (16). Potential secondary effects in each of these manipulated strains confound attempts to measure the effect of rRNA operon copy number on physiological traits.
To assess the consequences of variation in rRNA operon copy number, we developed a strategy for deleting rRNA operons in order to create strains of E. coli with one or two rRNA operons (including promoters) deleted that were otherwise isogenic (Fig. 1). rRNA operons rrnA and rrnB were targeted for deletion in order to minimize perturbation of the tRNA gene pool. The tRNAs present in the internally transcribed spacers of rrnA and rrnB are replicated in other rRNA operons, and neither of these operons contains distally located tRNA genes (31). rrnA and rrnB are also located adjacent to each other on the chromosome, reducing differences in gene dosage that could occur in cells with multiple DNA replication loops (10, 15). The construction of these strains enabled a direct measurement of the fitness effects associated with varying the rRNA operon copy number.
Results from batch culture competition experiments support the hypothesis that multiple rRNA operons are an advantage under culture conditions defined by fluctuating resource availability (Fig. 3; Table 2). As the number of rRNA operons was decreased from 7 to either 6 or 5 copies per chromosome, the ability of these cells to compete with the parental strain for the same limiting resource (glucose) decreased. Although differences between individual components of fitness for each strain were not always statistically significant, the cumulative changes in lag time and growth rate were sufficient to account for 98 to 99% of the relative fitness of the strains in simulated competition experiments modeled on Monod kinetics (Fig. 3).
The observation that deletion of rrnA or rrnB had a relatively small (1 to 2%) impact on the growth rate suggests that the mutant strains are able to compensate, in large part, for these deletions. If there were no means for compensation, deletion of any one of the seven rRNA operons would result in a 14% decrease in growth rate. The most likely mechanism to account for that compensation is the recognized capacity for increased expression from the remaining intact operons (13). However, as expression of rRNA operons approaches maximal levels, the capacity to compensate for deletions should decrease. To test for this compensatory capacity at different expression levels, the growth rate of each of the strains was varied. As predicted, deletion strains were able to compensate completely for the deletion of either rrnA or rrnB or both at low growth rates, but as expression of the rRNA operons was increased at higher growth rates, the deleterious effect of the deletions became increasingly pronounced (Fig. 5). Although the decrease in maximal growth rate that coincided with fewer rRNA operons might appear to be relatively small, sensitivity analyses of our batch culture competition model corroborated those by Vasi et al., which demonstrate that relative fitness is extremely sensitive to even minor changes in growth rates under batch culture conditions (42).
While it may be tempting to assert that the advantage conferred by multiple rRNA genes is simply the capacity for higher growth rates, deletion of rrnA or rrnB in E. coli had a similar impact on a second component of fitnesslag time. In nature, the capacity for rapid response to resources (i.e., shorter lag time) may be more important than maximal growth rate, particularly in environments with fluctuating resource availability. Lag times measured in these experiments may be more informative than measurements of "shift-up" times in strains of E. coli with inactivated rRNA operons (13). While "shift-up" times reflect the amount of time needed to shift from one exponential-growth rate to a higher growth rate, lag times include the time necessary to initiate growth from a nongrowing physiological statea common condition for microbes in nature.
The results from competition experiments in the chemostat did not reveal a selective advantage for the knockout strains that might have resulted from faster growth due to more efficient utilization of resources (Fig. 4). There was, however, a differential effect on fitness in chemostat culture with the deletion of the rrnA versus the rrnB operon. Deletion of the rrnA operon resulted in a substantial decrease in fitness, whereas no difference in fitness was found for the deletion of the rrnB operon. Although we do not know the cause of this differential effect, it might be the first significant evidence of differential properties between rRNA operons in E. coli (14), and it directly correlates with differences in measured Ks and therefore with estimates of J in these strains (Table 3). The J parameter represents the lowest concentration of a limiting resource where the specific growth rate of an organism is equal to the dilution rate, and it can be used to accurately predict the outcome of chemostat competition experiments (19). J parameter values for the
rrnB strain were within 8% of those for the parental strains, which is consistent with the coexistence of these two strains in chemostat culture competition experiments (Fig. 4). The more substantial differences in J (>24%) between rrnA deletion strains and the parental strain were consistent with the much lower relative fitness of
rrnA strains in chemostats. Perhaps it is not surprising that strains with fewer rRNA operons were unable to make more efficient use of the limiting resources, considering that these experiments were carried out on a line of E. coli with cell machinery that had evolved for 10,000 generations to take advantage of periodic fluctuations in resource availability (42). The large effect of deleting the rrnA operon in
rrnA and
rrnAB strains, relative to the parental strain or the
rrnB strain, was unexpected and remains unresolved.
It is necessary to consider other possible explanations for the experimental results reported here. The fact that genomic DNA from an intermediate construct in E. coli K-12 was used during generalized transduction with phage P1vir raises the possibility that DNA other than the targeted rRNA operons could have been incorporated into the recipient strain and could have been responsible for the observed changes in phenotype. We view this as an unlikely source of the variation, because the E. coli recipient strain and strain K-12 are closely related within one subgroup of E. coli (22), and there is only a single known region of intraspecific sequence divergence (the phn operon) located within 100 kb of the rrnA and rrnB operons. We also note that a similar decrease in growth rate was observed when the rrnB operon was inactivated in the same parental strain of E. coli by the insertion of the sacB::neo marker allele into the 16S gene (20). Since the effect on relative fitness of the inactivation of the rrnB operon was indistinguishable from that of BS200 and BS201, used in the present study, we conclude that changes in growth rate and lag time are attributable to deletion of the rRNA operons.
It is also not possible to distinguish the effect of the deletion of rRNA genes from the deletion of tRNA genes that are found in the ITS in rrnA and rrnB operons. The rrnA ITS region contains the two tRNA genes tRNAAla and tRNAIle, whereas the rrnB ITS region contains the tRNAGlu gene (32). Although the tRNA genes found in the ITS regions of rrnA and rrnB are repeated in the other rRNA operons found elsewhere on the chromosome, the pools of the respective tRNAs may have been perturbed in these strains.
In summary, the research reported here supports the hypothesis that the rRNA operon copy number is a component of bacterial life histories. The ecological strategies of bacteria with many rRNA operons most likely include quick response and fast growth upon an influx of resources, as opposed to a strategy of more efficient utilization of limited resources. We are continuing to investigate the mechanisms behind a potential tradeoff for rRNA operon copy numbers in microbes in order to better understand the role that this trait plays in evolution and life history in the microbial world.
We thank Richard E. Lenski for providing essential guidance on the measurement of relative fitness in batch culture and the parental strain REL4547. We also thank Les Dethlefsen for many valuable conversations regarding fitness.
|
|
|---|
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»