| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |


* Department of Botany, Miami University, Oxford, Ohio 45056-3653,
University of Puerto Rico, Mayaguez, Puerto Rico, 00681 and
Molecular and Computational Biology Program, Department of Biological Sciences, University of Southern California, Los Angeles, California 90089-2910
1 Corresponding author: Department of Botany, Miami University, Oxford, OH 45056.
E-mail: hstevens{at}muohio.edu
| ABSTRACT |
|---|
|
|
|---|
Dozens of hypotheses have been generated to explain variation in diversity along resource gradients, and most of them have focused on how resource concentration, per se, can influence competitive interactions and predator–prey dynamics (WAIDE et al. 1999). In nature, resource concentration gradients inevitably confound changes in ratios of potentially limiting factors, and this variation probably contributes to the variation in the patterns of diversity we see along resource and productivity gradients. It is critical that we determine experimentally the baseline patterns of diversity along gradients that are under our complete control. Only this will allow us to observe what tends to happen along resource gradients where only concentration, and not ratios, vary.
Resource gradients create not only gradients in biodiversity, but also gradients in population size, and this variation in population size may generate patterns of both genetic and species diversity (hereafter diversity) (STEVENS and CARSON 1999; VANDERMEULEN et al. 2001; SUDING et al. 2005; BAZIN et al. 2006). Neutral theory predicts that larger populations (perhaps associated with high resource levels) will have greater heterozygosity (CROW and KIMURA 1970). Further, we know that larger populations have higher rates of adaptation (DE VISSER and ROZEN 2005), and this has the potential to increase the rate of adaptive diversification. Clonal interference in asexual populations (GERRISH and LENSKI 1998) and the Hill–Robertson effect in sexual populations (HILL and ROBERTSON 1966), in which more beneficial mutations compete with each other and slow the rate of selective sweeps, will also maintain greater diversity in larger populations. It is also widely accepted, however, that microbial populations can rapidly lose diversity through selective sweeps (ATWOOD et al. 1951) and hitchhiking (GILLESPIE 2001)—increased evolutionary rates caused by larger population size would serve only to exacerbate this loss of diversity. Experimental verification of the effect of population size on diversification rate is lacking, yet would have profound implications for global patterns of biodiversity along resource and latitudinal gradients.
Here we describe dynamics of population size and per capita genetic diversity (hereafter simply "heterozygosity") in highly dynamic long-term stationary phase (LTSP) cultures of Escherichia coli. The LTSP of microbial population dynamics follows the better characterized lag, exponential, stationary, and death phases (FINKEL 2006). LTSP is characterized by relatively low population sizes (106–107 cells · ml–1), rapid turnover, and dynamic evolutionary change (ZAMBRANO and KOLTER 1996; YEISER et al. 2002; FINKEL 2006). This population turnover and evolutionary dynamism begins to appear almost immediately, even in overnight cultures, and continues indefinitely for days, weeks, and years.
One of the most notable changes is the accumulation of beneficial mutations that confer a competitive advantage to surviving cells (ZAMBRANO et al. 1993; FINKEL and KOLTER 1999; FINKEL et al. 2000; YEISER et al. 2002). Frequently, when cells are grown under aerobic conditions in Luria–Bertani medium, the first gene to receive a detectable mutation is rpoS, the gene encoding the stationary phase-specific or stress-response sigma factor RpoS (ZAMBRANO et al. 1993; HENGGE-ARONIS 2000; FARRELL and FINKEL 2003). For example, in one experiment, individual LTSP cultures contained 6–7 rpoS mutant alleles after only 10 days of incubation (M. J. FARRELL, V. PALCHEVSKIY, T. JESSNER and S. E. FINKEL, unpublished results). Such mutations frequently result in different degrees of attenuated RpoS activity.
We used well-mixed LTSP cultures of E. coli populations to test (i) whether heterozygosity varies with resource concentration in an unstructured environment, and (ii) whether larger populations associated with high resource concentrations could account for the observed pattern of heterozygosity. Within a single culture, different rpoS mutations frequently result in different degrees of attenuated RpoS activity, and we can observe directly the degree of attenuation using a reporter gene. For our experiments, observed variation (based on subsamples of equal cell numbers) in the expression of this reporter served as our measure of the heterozygosity of LTSP cultures of E. coli. In our first experiment, we established an initial resource-concentration gradient and assayed population size and heterozygosity for 10 days. In our second experiment, we manipulated both initial population size and initial resource concentration independently and then assayed the response of both population size and heterozygosity to this manipulation for 10 days.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Identification of rpoS mutants under different growth conditions:
Derivatives of the RpoS-reporter strain ZK918 were used to observe phenotypic changes due to mutation at rpoS. These strains carry a bolA::lacZ reporter fusion (BOHANNON et al. 1991). bolA gene expression is under control of the general regulatory protein RpoS (HENGGE-ARONIS 2000); maximal lacZ expression is dependent on wild-type alleles of rpoS. To determine the RpoS phenotype of evolved strains, individual colonies are plated on MacConkey agar (Difco). Colonies are scored as wild type or a category of attenuated expression on the basis of their color: red for wild-type cells, and dark pink, light pink, or white for different attenuated mutant alleles (FARRELL and FINKEL 2003). The status of the rpoS gene can be further confirmed by performing the catalase "bubble assay" (expression of the KatE catalase is dependent on RpoS) (BOHANNON et al. 1991; ZAMBRANO and KOLTER 1993).
Use of the bolA::lacZ reporter as a phenotypic marker of genetic variation, provides a very robust and reliable indicator of rpoS activity. RpoS activity is a trait that comprises a critical component of a cell's complex fitness response to stressful environments (HENGGE-ARONIS 2000; NOTLEY-MCROBB et al. 2002; FARRELL and FINKEL 2003). As a measure of phenotypic, rather than genotypic, diversity, it reflects traits that more directly influence interactions and responses to abiotic factors. Our observations of the reporter trait, however, comprise only four categories derived from a continuum of likely types: "red," "dark pink," "light pink," and "white." This may be akin to describing the diversity of a plant community with four functional types, for instance, annual weeds, perennial weeds, early-successional trees, and late-successional trees. Such gross differences provide a sound starting point for understanding ecologically relevant phenotypic diversity in an assemblage of species or phenotypes within a species. Nonetheless, such a coarse classification may result in a poor resolution of true taxonomic diversity and composition. If we observe patterns of diversity at such a coarse resolution, they are likely to represent very robust patterns.
Assaying population size and diversity of rpoS mutants:
Each culture was sampled, serially diluted, and spread onto MacConkey agar plates to determine the color of individual colonies and population densities. One plate per culture, containing 150–300 colonies, was selected for enumeration. The viable cells were counted and scored by their colors (red for wild type and dark pink, light pink, and white colonies for different levels of attenuated activity).
Nei's heterozygosity (
, aka Simpson's diversity) measured per capita culture diversity (SIMPSON 1949; NEI 1987; STEVENS et al. 2003) using the proportions, pi, of sample colonies of each color. Because we have a finite number of observed mutant phenotypes (white, light pink, dark pink, plus the red wild type), maximum heterozygosity (
), as well as the heterozygosity expected under neutrality (CROW and KIMURA 1970) are both <1. All assays of replicate cultures were based on approximately equal numbers of cells, and a regression (not shown) of sample size and diversity revealed no pattern.
Because of our use of both a diversity index and similar sample sizes across replicates, our results reflect the relative abundance of phenotypes, i.e., diversity or population-level heterozygosity. The results do not reflect sampling artifacts, wherein larger populations contain more total mutants—rather, they reflect variation in diversity.
Resource gradient experiment:
In the first experiment, LB concentrations for each experiment included the standard concentration (1x) and four twofold serial dilutions (
). Each strain was replicated twice at each resource concentration, for a total of 20 tubes. Sterile distilled water was added to each culture periodically to maintain constant volume. Population sizes and heterozygosity were assayed each day for 10 days, using methods described above. Although we report results for heterozygosity, we obtained directly analogous results when using other measures of diversity and dominance, such as the proportion of the sample composed of the wild type (red colonies). General linear mixed models with repeated measures tested effects of resource concentration, sampling day, and their interaction on the logarithm of total population size and the arcsin square root of diversity. The most parsimonious random-effects structure was determined with likelihood ratio tests (PINHEIRO and BATES 2000). The most parsimonious models included culture tubes as independent random variables, and variances estimated independently for each nutrient concentration. All analyses were performed with the R language (R DEVELOPMENT CORE TEAM 2007) using the nlme package (PINHEIRO and BATES 2000).
Population size manipulation:
We tested directly the independent effects of population size and resource concentration on diversity by growing E. coli populations at different volumes, at different resource concentrations. The selected volumes and resource concentrations (see below) manipulated population size independently of resource concentration. In particular, they allowed the population size of the low-concentration–high-volume treatment to exceed the population size of the high-concentration–low-volume treatment. We established two replicates of each strain at each unique volume–concentration combination: 7 ml and 70 ml of
x standard LB concentration, and 0.7 ml and 7 ml of standard LB. Population sizes and Nei's heterozygosity were assayed each day for 10 days, using methods described above. Statistical methods are the same as those described above, and in addition, single degree-of-freedom contrasts used "dummy" (i.e., "treatment") coding to test differences among unique treatment combinations for total population size on day 1, and heterozygosity on day 10.
| RESULTS |
|---|
|
|
|---|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Increased population size can enhance heterozygosity in several ways, and yet this flies in the face of conventional wisdom of laboratory microbial populations. Sweeps by mutants with higher fitness are often the expectation in such populations. In contrast, our results are consistent with expectations with several other mechanisms, including neutral mutation. In addition, clonal interference (GERRISH and LENSKI 1998), mutator phenotypes (TENAILLON et al. 1999), and divergent selection (DIECKMANN and DOEBELI 1999) may all contribute to increased heterozygosity of larger populations in unstructured environments. This increased heterozygosity arose in spite of mechanisms that have the potential to reduce heterozygosity in large populations, such as increased hitchhiking (GILLESPIE 2001) and more rapid adaptive evolution that could result in more rapid sweeps. Future research must partition the relative importance of these mechanisms in unstructured environments. Regardless, however, of the details of these processes, population size does not explain all of the increased heterozygosity at higher resource concentrations.
Higher resource concentration alone, independent of its effect on initial population size, caused higher heterozygosity, and several processes might contribute to this. First, populations in higher resource concentrations may have higher mutation rates for at least two reasons aside from the effect of population size. Higher resource concentrations may cause more rapid turnover of individuals within populations thereby increasing the number of mutants appearing in a culture per unit time. Higher resource concentrations may also, through a variety of mechanisms, result in a higher concentration of free radicals (e.g.,
, ·OH) that may increase DNA damage and mutation rates (FINKEL 2003; NEMOTO and FINKEL 2004).
In addition to its effects on mutation rates, high resource concentrations may enhance the maintenance of heterozygosity in at least three ways. First, increased population size increases the absolute abundance of the rarest types, decreasing the probability of stochastic extinction (ABRAMS 1995). For this to be significant, however, there would need to be a mechanism for these types to increase in frequency enough to be observed in our samples. Second, increased resources may increase the abundance of rare resources or resource combinations used by rare phenotypes and allow these specialists to increase in frequency (ABRAMS 1995). Third, it has long been suggested that even with no change in the relative abundance or ratios of resources, increased concentration of all resources favors specialization (MACARTHUR and PIANKA 1966; DIECKMANN and DOEBELI 1999). This specialization would permit a greater variety of mutants to persist on the same variety (but greater concentration) of resources. This adaptive diversification associated with narrowing of niche breadth is most likely to occur if there are not other costs associated with this change (ACKERMANN and DOEBELI 2004). This seems quite plausible in this system with a moderately complex medium (LB broth) and in which we altered resource concentration without changing the initial relative concentrations of any of the substrates.
Recent work on adaptive diversification in both asexual and sexually reproducing populations has demonstrated clearly the plausibility of sympatric diversification (DOEBELI et al. 2005). These advances suggest that the work presented here, using a clonally reproducing organism in a highly simplified unstructured laboratory environment has implications for not only bacteria, but for sexually reproducing, charismatic megafauna and flora.
The population size mechanism we describe here may have implications for the latitudinal gradient of taxonomic richness (WILLIG et al. 2003). If this mechanism applies generally, then regions or habitats with larger numbers of individuals will contain more species (or genera or genotypes) by virtue of having more individuals. Such areas may have more individuals by virtue of greater area alone or by virtue of higher population densities. The tropics, defined by climate, have a larger land surface (ROSENZWEIG 1992) and therefore may have more individuals than any other similarly uniform terrestrial habitat. Our population size hypothesis would therefore predict, all else being equal, higher diversity. Further, as climate periodically causes bottlenecks at higher latitudes, this would dramatically lower the effective population sizes at high latitudes, reducing the potential for diversification (JABLONSKI et al. 2006).
The population size mechanism could also have implications for patterns of diversity along other gradients. A unimodal pattern of diversity along productivity gradients might arise, for instance, if such habitats contain more individuals by virtue of greater area (VANDERMEULEN et al. 2001). It might also contribute to a mid-domain effect (RANGEL and DINIZ-FILHO 2005) in which diversity peaks at the midpoint of any region with physical geographic constraints. Dispersal limitation due to physical boundaries might limit effective population size, so that effective population size would be lower near boundaries. The population size hypothesis would predict lower rates of diversification near these boundaries.
The above list of plausible implications suffers from a paradoxical chicken-and-egg problem: if larger populations exhibit more rapid diversification, this will reduce the size of the derived populations through lineage splitting. The resulting smaller populations would have correspondingly slower rates of diversification. This possible slowing of diversification, however, could well be offset by increased diversification stemming from greater diversity itself (EMERSON and KOLM 2005). This might lead us to speculate that life seems destined for diversity: large populations diverge because they are large, and the resulting divergent groups further diversify because of their differences.
This work examines the processes of diversification and its subsequent maintenance over many generations, with implications for extant patterns on the globe. Extant patterns of diversity, however, reflect not only patterns of both diversification and its in situ maintenance, but also community assembly and species sorting following dispersal of allopatrically derived forms. Nonetheless, the implications of our results for community assembly seem clear: we can hypothesize that if particular conditions maintain a greater variety of sympatrically derived phenotypes, as we have shown, then it seems plausible that a greater variety of allopatrically derived phenotypes can be maintained following dispersal-mediated assembly. Whether this is generally the case, however, is untested.
It seems difficult to reconcile the arguably small differences in population size observed in the "real" world with the amazing differences in diversity between some contrasting regions of the planet. Other factors, such as temperature, may be more important than resource availability. Nonetheless, we might consider that speciation is an exponential (albeit slow) process (EMERSON and KOLM 2005), because a greater variety of genotypes results in a wider base from which to diverge. Therefore, small differences in diversification rates may result, over time, in large differences in diversity.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| LITERATURE CITED |
|---|
|
|
|---|
ABRAMS, P. A., 1995 Monotonic or unimodal diversity-productivity gradients: What does competition theory predict? Ecology 76: 2019–2027.[CrossRef]
ACKERMANN, M., and M. DOEBELI, 2004 Evolution of niche width and adaptive diversification. Evolution 58: 2599–2612.[Medline]
ATWOOD, K. C., L. K. SCHNEIDER and F. J. RYAN, 1951 Periodic selection in Escherichia coli. Proc. Natl. Acad. Sci. USA 37: 146–155.
BAZIN, E., S. GLEMIN and N. GALTIER, 2006 Population size does not influence mitochondrial genetic diversity in animals. Science 312: 570–572.
BOHANNON, D. E., N. CONNELL, J. KEENER, A. TORMO, M. ESPINOSA-URGEL et al., 1991 Stationary-phase-inducible "gearbox" promoters: differential effects of katF mutations and role of sigma 70. J. Bacteriol. 173: 4482–4492.
CROW, J. F., and M. KIMURA, 1970 An Introduction to Population Genetics Theory. Harper & Row, New York.
DE VISSER, J. A. G. M., and D. E. ROZEN, 2005 Limits to adapation in asexual populations. J. Evol. Biol. 18: 779–788.[CrossRef][Medline]
DIECKMANN, U., and M. DOEBELI, 1999 On the origin of species by sympatric speciation. Nature 400: 354–357.[CrossRef]
DOEBELI, M., U. DIECKMANN, J. A. J. METZ and D. TAUTZ, 2005 What we have also learned: adaptive speciation is theoretically plausible. Evolution 59: 691–695.[Medline]
EMERSON, B. C., and N. KOLM, 2005 Species diversity can drive speciation. Nature 434: 1015–1017.[CrossRef][Medline]
FARRELL, M. J., and S. E. FINKEL, 2003 The growth advantage in stationary-phase phenotype conferred by rpoS mutation is dependent on the pH and nutrient environment. J. Bacteriol. 185: 7044–7052.
FINKEL, S. E., 2006 Long-term survival during stationary phase: evolution and the GASP phenotype. Nat. Rev. Microbiol. 4: 113–120.[CrossRef][Medline]
FINKEL, S. E., and R. KOLTER, 1999 Evolution of microbial diversity during prolonged starvation. Proc. Natl. Acad. Sci. USA 96: 4023–4027.
FINKEL, S. E., E. R. ZINSER and R. KOLTER, 2000 Long-term survival and evolution in the stationary phase, pp. 231–238 in Bacterial Stress Responses, edited by G. STORZ and R. HENGGE-ARONIS. ASM, Washington, DC.
FINKEL, T., 2003 Oxidant signals and oxidative stress. Curr. Opin. Cell Biol. 15: 247–254.[CrossRef][Medline]
GERRISH, P. J., and R. E. LENSKI, 1998 The fate of competing beneficial mutations in an asexual population. Genetica 103: 127–144.[CrossRef]
GILLESPIE, J. H., 2001 Is the population size of a species relevant to its evolution? Evolution 55: 2161–2169.[CrossRef][Medline]
HEDRICK, P. W., 1986 Genetic polymorphism in heterogeneous envrionments: a decade later. Annu Review of Ecology and Systematics 17: 535–566.[CrossRef]
HEDRICK, P. W., M. E. GINEVAN and E. P. EWING, 1976 Genetic polyorphism in a heterogeneous environment. Annu. Rev. Ecol. Syst. 7: 1–32.[Medline]
HENGGE-ARONIS, R., 2000 The general stress response in Escherichia coli, pp. 161–178 in Bacterial Stress Responses, edited by G. STORZ and R. HENGGE-ARONIS. ASM, Washington, DC.
HILL, W. G., and A. ROBERTSON, 1966 The effect of linkage on limits to artificial selection. Genet. Res. 8: 269–294.[Medline]
HILLEBRAND, H., 2004 On the generality of the latitudinal diversity gradient. Am. Nat. 163: 192–211.[CrossRef][Medline]
JABLONSKI, D., K. ROY and J. W. VALENTINE, 2006 Out of the tropics: evolutionary dynamics of the latitudinal gradient. Science 314: 102–106.
KASSEN, R., A. BUCKLING, G. BELL and P. B. RAINEY, 2000 Diversity peaks at intermediate productivity in a laboratory microcosm. Nature 406: 508–512.[CrossRef][Medline]
MACARTHUR, R., and E. PIANKA, 1966 On optimal use of a patchy environment. Am. Nat. 100: 603–609.[CrossRef]
MITTELBACH, G. G., C. F. STEINER, S. M. SCHEINER, K. L. GROSS, H. L. REYNOLDS et al., 2001 What is the observed relationship between species richness and productivity? Ecology 82: 2381–2396.[CrossRef]
NEI, M., 1987 Molecular Evolutionary Genetics. Columbia University Press, New York.
NEMOTO, S., and T. FINKEL, 2004 Aging and the mystery at Arles. Nature 429: 149–152.[CrossRef][Medline]
NOTLEY-MCROBB, L., T. KING and T. FERENCI, 2002 rpoS mutations and loss of general stress resistance in Escherichia coli populations as a consequence of conflict between competing stress types. J. Bacteriol. 184: 806–811.
PENNISI, E., 2005 What determines species diversity? Science 309: 90.
PINHEIRO, J., and D. BATES, 2000 Mixed-Effects Models in S and S-PLUS. Springer, New York.
R DEVELOPMENT CORE TEAM, 2007 R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna.
RANGEL, T. F. L. V. B., and J. A. F. DINIZ-FILHO, 2005 Neutral community dynamics, the mid-domain effect and spatial patterns in species richness. Ecol. Lett. 8: 783–790.[CrossRef]
ROSENZWEIG, M. L., 1992 Species diversity gradients: we know more and less than we thought. J. Mammal. 73: 715–730.[CrossRef]
SIMPSON, E. H., 1949 Measurement of diversity. Nature 163: 688.
STEVENS, M. H. H., and W. P. CARSON, 1999 The significance of assemblage level thinning for species richness. J. Ecol. 87: 490–502.[CrossRef]
STEVENS, M. H. H., O. L. PETCHEY and P. E. SMOUSE, 2003 Stochastic relations between species richness and the variability of species composition. Oikos 103: 479–488.[CrossRef]
SUDING, K. N., S. L. COLLINS, L. GOUGH, C. CLARK, E. E. CLELAND et al., 2005 Functional- and abundance-based mechanisms explain diversity loss due to n fertilization. Proc. Natl. Acad. Sci. USA 102: 4387–4392.
TENAILLON, O., B. TOUPANCE, H. LE NAGARD, F. TADDEI and B. GODELLE, 1999 Mutators, population size, adaptive landscape and the adaptation of asexual populations of bacteria. Genetics 152: 485–493.
TILMAN, D., 1982 Resource Competition and Community Structure (Monographs in Population Biology). Princeton University Press, Princeton, NJ.
VANDERMEULEN, M., A. HUDSON and S. M. SCHEINER, 2001 Three evolutionary hypotheses for the hump-shaped productivity-diversity curve. Evol. Ecol. Res. 3: 379–392.
WAIDE, R. B., M. R. WILLIG, C. F. STEINER, G. MITTELBACH, L. GOUGH et al., 1999 The relationship between productivity and species richness. Annu. Rev. Ecol. Syst. 30: 257–300.[CrossRef]
WILLIG, M. R., D. M. KAUFMAN and R. D. STEVENS, 2003 Latitudinal gradients of biodiversity: pattern, process, scale, and synthesis. Annu. Rev. Ecol. Syst. 34: 273–309.[CrossRef]
YEISER, B., E. D. PEPPER, M. F. GOODMAN and S. E. FINKEL, 2002 SOS-induced DNA polymerases enhance long-term survival and evolutionary fitness. Proc. Natl. Acad. Sci. USA 99: 8737–8741.
ZAMBRANO, M. M., and R. KOLTER, 1993 Escherichia coli mutants lacking Nadh Dehydrogenase-I have a competitive disadvantage in stationary phase. J. Bacteriol. 175: 5642–5647.
ZAMBRANO, M. M., and R. KOLTER, 1996 GASPing for life in stationary phase. Cell 86: 181–184.[CrossRef][Medline]
ZAMBRANO, M. M., D. A. SIEGELE, M. ALMIRON, A. TORMO and R. KOLTER, 1993 Microbial competition—Escherichia coli mutants that take over stationary phase cultures. Science 259: 1757–1760.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |