This Article
Right arrow Full Text (PDF)
Right arrow An erratum has been published
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hughes, J. B.
Right arrow Articles by Bohannan, B. J. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hughes, J. B.
Right arrow Articles by Bohannan, B. J. M.
Agricola
Right arrow Articles by Hughes, J. B.
Right arrow Articles by Bohannan, B. J. M.

Next Article 

Applied and Environmental Microbiology, October 2001, p. 4399-4406, Vol. 67, No. 10
0099-2240/01/$04.00+0   DOI: 10.1128/AEM.67.10.4399-4406.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.

MINIREVIEW

Counting the Uncountable: Statistical Approaches to Estimating Microbial Diversity

Jennifer B. Hughes,* Jessica J. Hellmann,dagger Taylor H. Ricketts, and Brendan J. M. Bohannan

Department of Biological Sciences, Stanford University, Stanford, California 94305-5020


    INTRODUCTION
Top
Introduction
Discussion
References

All biologists who sample natural communities are plagued with the problem of how well a sample reflects a community's "true" diversity. New genetic techniques have revealed extensive microbial diversity that was previously undetected with culture-dependent methods and morphological identification (reviewed in references 2 and 46), but exhaustive inventories of microbial communities still remain impractical. As a result, we must rely on samples to inform us about the actual diversity of microbial communities.

Ecologists studying the diversity of macroorganisms also face this estimation problem and have designed tools to deal with the problems of sampling (14, 25, 36). Sparked by the availability of microbial diversity data, interest is emerging in applying these tools to microbes. Reliable estimates of microbial diversity would offer a means to address once intractable questions, such as what processes control microbial diversity? How do microbial communities affect ecosystem functioning? How are human beings affecting microbial communities?

Several microbial studies have used diversity indices (39, 44), estimated species richness (33, 43), and compared sample diversity with rarefaction curves (19, 40). Still others have proposed new diversity statistics specific to microbial samples (69). Despite the recent interest, however, the success of these tools has not yet been evaluated for microbial communities, and other potential approaches remain to be explored.

Here we compare the utility of various statistical approaches for assessing the diversity of microbial communities. First, we show examples of communities in which macroorganisms are as diverse as some microbial communities, suggesting that diversity estimation methods developed for macroorganisms may be appropriate for microbial samples. Second, we review these methods and discuss how to evaluate the success of diversity estimators for microbial communities for which the true diversity is unknown. We argue that even without knowing the "truth," it is possible to rigorously compare relative diversity among communities. Finally, we apply some of these diversity measures to microbial data sets and examine how the confidence of the measures changes with sample size.

Throughout the paper, we use the term diversity to mean richness, or the number of types. We also use the term microbial with bacteria in mind, although much of the discussion is applicable to other microbes. For clarity, we will often refer to species as the measured unit of diversity, but our discussion can be applied to any operational taxonomic units (OTUs), such as the number of unique terminal restriction fragments (35) or number of 16S ribosomal DNA (rDNA) sequence similarity groups (41). Finally, we are concerned here with estimating richness and do not address how this diversity is related to functional diversity (1).


    ARE MICROBES TOO DIVERSE TO COUNT?

In any community, the number of types of organisms observed increases with sampling effort until all types are observed. The relationship between the number of types observed and sampling effort gives information about the total diversity of the sampled community. This pattern can be visualized by plotting an accumulation or a rank-abundance curve.

An accumulation curve is a plot of the cumulative number of types observed versus sampling effort. Figure 1 shows the accumulation curves for samples from five communities: bacteria from a human mouth (33), soil bacteria (6), tropical moths (56), tropical birds (J. B. Hughes, unpublished data), and temperate forests (26). We standardized the data sets by the number of individuals collected to compare the shapes of the curves. Differences in the richness and relative abundances of species in the sampled communities underlie the differences in the shape of the curves. Because all communities contain a finite number of species, if the surveyors continued to sample, the curves would eventually reach an asymptote at the actual community richness (number of types). Thus, the curves contain information about how well the communities have been sampled (i.e., what fraction of the species in the community have been detected). The more concave-downward the curve, the better sampled the community.


View larger version (29K):
[in this window]
[in a new window]
 
FIG. 1.   Accumulation curves for Michigan plants (cross ; n = 1,783) (26), Costa Rican birds (black-triangle; n = 5,007) (J. B. Hughes, unpublished data), human oral bacteria (open circle ; n = 264) (33), Costa Rican moths (; n = 4,538) (56), and East Amazonian soil bacteria (; n = 98) (6). Curves are averaged over 100 simulations using the computer program EstimateS and are standardized for the number of individuals and species observed.

The idea that microbial diversity cannot be estimated comes from the fact that many microbial accumulation curves are linear or close to linear because of high diversity, small sample sizes, or both. Indeed, the accumulation curve of East Amazonian soil bacteria represents the worst-case scenario (Fig. 1). Every individual identified was a different type; therefore, this sample supplies no information about how well the community has been sampled. At the other extreme, the plant and bird communities plotted in Fig. 1 are well sampled, and the samples therefore contain considerable information about total richness. The two intermediate curves provide the most telling comparison, however. Even though the moth sample is much larger than the mouth bacteria sample (4,538 versus 264 individuals), the shape of the curves is similar. In other words, the communities have been sampled with roughly equivalent intensity relative to their overall richness.

Another way to compare how well communities have been sampled is to plot their rank-abundance curves. The species are ordered from most to least abundant on the x axis, and the abundance of each type observed is plotted on the y axis. The moth and soil bacteria communities exhibit a similar pattern (Fig. 2), one that is typical of superdiverse communities such as tropical insects. A few species in the sample are abundant, but most are rare, producing the long right-hand tail on the rank-abundance curve.


View larger version (19K):
[in this window]
[in a new window]
 
FIG. 2.   Rank-abundance curves for (a) tropical moths (n = 4,538) (56) and (b) temperate soil bacteria (n = 137) (39). The two most abundant species of moths (396 and 173 individuals) are excluded from panel a to shorten the y axis.

If these organisms were sampled on the same spatial scale, there is no doubt that soil bacterial diversity would be higher than moth diversity. These comparisons suggest, however, that our ability to sample bacterial diversity in a human mouth or in a few grams of some soils may be similar to our ability to sample moth diversity in a few hundred square kilometers of tropical forest. Thus, at least for some communities, microbiologists may be able to coopt techniques that ecologists use to estimate and compare the richness of macroorganisms.

Ultimately, microbes---like tropical insects---are too diverse to count exhaustively. While it would be useful to know the actual diversity of different microbial communities, most diversity questions address how diversity changes across biotic and abiotic gradients, such as disturbance, productivity, area, latitude, and resource heterogeneity. The answers to these questions require knowing only relative diversities among sites, over time, and under different treatment regimens. Using this approach, the relationships between insect diversity and many environmental variables have been well studied (50, 57, 63, 64), even though estimates of the total number of insect species range over three orders of magnitude (22, 54).


    SOME POSSIBLE TOOLS: RAREFACTION AND RICHNESS ESTIMATORS

A variety of statistical approaches have been developed to compare and estimate species richness from samples of macroorganisms. In this section, we consider the suitability of four approaches for microbial diversity studies.

The first approach, rarefaction, has been adopted recently by a number of microbiologists (4, 19, 40). Rarefaction compares observed richness among sites, treatments, or habitats that have been unequally sampled. A rarefied curve results from averaging randomizations of the observed accumulation curve (25). The variance around the repeated randomizations allows one to compare the observed richness among samples, but it is distinct from a measure of confidence about the actual richness in the communities.

In contrast to rarefaction, richness estimators estimate the total richness of a community from a sample, and the estimates can then be compared across samples. These estimators fall into three main classes: extrapolation from accumulation curves, parametric estimators, and nonparametric estimators (14, 23, 47). To date, we have found only two studies that apply richness estimators to microbial data (33, 43).

Most curve extrapolation methods use the observed accumulation curve to fit an assumed functional form that models the process of observing new species as sampling effort increases. The asymptote of this curve, or the species richness expected at infinite effort, is then estimated. These models include the Michaelis-Menten equation (13, 51) and the negative exponential function (61). The benefit of estimating diversity with such extrapolation methods is that once a species has been counted, it does not need to be counted again. Hence, a surveyor can focus effort on identifying new, generally rarer, species. The downside is that for diverse communities in which only a small fraction of species is detected, several curves often fit equally well but predict very different asymptotes (61). This approach therefore requires data from relatively well sampled communities, so at present curve extrapolation methods do not seem promising for estimating microbial diversity in most natural environments.

Parametric estimators are another class of estimation methods. These methods estimate the number of unobserved species in the community by fitting sample data to models of relative species abundances. These models include the lognormal (49) and Poisson lognormal (7). For instance, Pielou (48) derived an estimator that assumes species abundances are distributed lognormally; that is, if species are assigned to log abundance classes, the distribution of species among these classes is normal. By fitting sample data to the lognormal distribution, the parameters of the curve can be evaluated. Pielou's estimator uses these parameter values to estimate the number of species that remain unobserved and thereby estimate the total number of species in the community.

There are three main impediments to using parametric estimators for any community. First, data on relative species abundances are needed. For macroorganisms, often only the presence or absence of a species in a sample or quadrat is recorded. In contrast, data on relative OTU abundances of microbes are often collected (see discussion below about potential biases). Second, one has to make an assumption about the true abundance distribution of a community. Although most communities of macroorganisms seem to display a lognormal pattern of species abundance (17, 36, 66), there is still controversy as to which models fit best (24, 30). In the absence of a variety of large microbial data sets, it is not clear which, if any, of the proposed distribution models describe microbial communities. Finally, even if one of these models is a good approximation of relative abundances in microbial communities, parametric estimators require large data sets to evaluate the distribution parameters. The largest microbial data sets currently available include only a few hundred individuals.

The final class of estimation methods, nonparametric estimators, is the most promising for microbial studies. These estimators are adapted from mark-release-recapture (MRR) statistics for estimating the size of animal populations (32, 59). Nonparametric estimators based on MRR methods consider the proportion of species that have been observed before ("recaptured") to those that are observed only once. In a very diverse community, the probability that a species will be observed more than once will be low, and most species will only be represented by one individual in a sample. In a depauperate community, the probability that a species will be observed more than once will be higher, and many species will be observed multiple times in a sample.

The Chao1 and abundance-based coverage estimators (ACE) use this MRR-like ratio to estimate richness by adding a correction factor to the observed number of species (9, 11). (For reviews of these and other nonparametric estimators, see Colwell and Coddington [14] and Chazdon et al. [12].) For instance, Chao1 estimates total species richness as
S<SUB><UP>Chao1</UP></SUB>=S<SUB><UP>obs</UP></SUB>+<FR><NU>n<SUP>2</SUP><SUB>1</SUB></NU><DE>2n<SUB>2</SUB></DE></FR>
where Sobs is the number of observed species, n1 is the number of singletons (species captured once), and n2 is the number of doubletons (species captured twice) (9). Chao (9) noted that this index is particularly useful for data sets skewed toward the low-abundance classes, as is likely to be the case with microbes.

The ACE (10) incorporate data from all species with fewer than 10 individuals, rather than just singletons and doubletons. ACE estimates species richness as
S<SUB><UP>ACE</UP></SUB>=S<SUB><UP>abund</UP></SUB>+<FR><NU>S<SUB><UP>rare</UP></SUB></NU><DE>C<SUB><UP>ACE</UP></SUB></DE></FR>+<FR><NU>F<SUB>1</SUB></NU><DE>C<SUB><UP>ACE</UP></SUB></DE></FR> &ggr;<SUP>2</SUP><SUB><UP>ACE</UP></SUB>
where Srare is the number of rare samples (sampled abundances <= 10) and Sabund is the number of abundant species (sampled abundances >10). Note that Srare + Sabund equals the total number of species observed. CACE = 1 - F1/Nrare estimates the sample coverage, where F1 is the number of species with i individuals and N<SUB>rare</SUB><IT>=&Sgr;</IT><SUB><IT>i=1</IT></SUB><SUP><IT>10</IT></SUP><IT>iF<SUB>i</SUB>.</IT> Finally,
&ggr;<SUP>2</SUP><SUB><UP>ACE</UP></SUB>=<UP>max</UP><FENCE><FR><NU>S<SUB><UP>rare</UP></SUB><LIM><OP>∑</OP><LL>i=1</LL><UL>10</UL></LIM>i(i−1)F<SUB>i</SUB></NU><DE>C<SUB><UP>ACE</UP></SUB> (N<SUB><UP>rare</UP></SUB>) (N<SUB><UP>rare</UP></SUB>−1)</DE></FR> −1, 0</FENCE>
which estimates the coefficient of variation of the Fi's (R. Colwell, User's Guide to EstimateS 5 [http://viceroy.eeb.uconn.edu/estimates]).

Both Chao1 and ACE underestimate true richness at low sample sizes. For example, the maximum value of SChao1 is (S2obs + 1)/2 when one species in the sample is a doubleton and all others are singletons. Thus, SChao1 will strongly correlate with sample size until Sobs reaches at least the square root of twice the total richness (14).


    EVALUATING RICHNESS ESTIMATORS

Given the variety of possible diversity estimators, how does one evaluate their utility? Clearly, the most desirable estimator is one that is both precise and unbiased. Precision describes the variation of the estimates from all possible samples that can be taken from the population. Bias describes the difference between the expected value of the estimator and the true, unknown richness of the community being sampled (in other words, whether the estimator consistently under- or overestimates the true richness).

To test for bias, one needs to know the true richness to compare against the sample estimates. As yet, this comparison is impossible for microbes, because no communities have been exhaustively sampled. The bias of richness estimators has only been tested in a few natural communities in which the exact abundance of every species in an area is known (12, 14, 15, 26, 47).

In contrast, precision is a relatively simple property to assess. With multiple samples (or one large sample) from a microbial community, the variance of microbial richness estimates can be calculated and compared. Moreover, most ecological questions require only comparisons of relative diversity. For these questions, an estimator that is consistent with repeated sampling (is precise) is often more useful than one that on average correctly predicts true richness (has the lowest bias). Thus, if we use diversity measures for relative comparisons, we avoid the problem of not being able to measure bias. (This assumes that the bias of an estimator does not differ so radically among communities that it disrupts the relative order of the estimates. In the absence of alternative evidence, this initial assumption seems appropriate.)

Chao (8) derives a closed-form solution for the variance of SChao1:
Var(S<SUB><UP>Chao1</UP></SUB>)=n<SUB>2</SUB><FENCE><FR><NU>m<SUP>4</SUP></NU><DE>4</DE></FR>+m<SUP>3</SUP>+<FR><NU>m<SUP>2</SUP></NU><DE>2</DE></FR></FENCE>,<UP> where </UP>m=<FR><NU>n<SUB><UP>1</UP></SUB></NU><DE>n<SUB><UP>2</UP></SUB></DE></FR>
This formula estimates the precision of Chao1; that is, it estimates the variance of richness estimates that one expects from multiple samples. A closed-form solution of variance for the ACE has not yet been derived.

Comparisons of relative species richness based on rarefaction may seem more reliable than comparisons using extrapolations that require a number of assumptions, but rarefaction is limited for two reasons. First, rarefaction compares samples, not communities. The error bars around a rarefaction curve describe the variation due to reordering of subsamples within the collected sample, not the precision of the observed richness. In contrast, a measure of precision would describe the variation in the number of species expected to be observed if the community were sampled repeatedly. It is possible to estimate the precision of rarefaction curves, for instance, by bootstrapping (20). Error bars derived by this method allow the detection of significant differences in observed richness between communities.

Second, the rank order of observed richness values does not necessarily correspond to relative total richness, because rarefaction analyses do not exclude the possibility that the species accumulation curves cross at a higher sample size (34). In contrast, species richness estimators take the shape of the accumulation curve into account to determine total richness. Thus, in theory these estimators can predict a crossover of the accumulation curves and thereby better predict relative total richness.


    CASE STUDIES

In terms of both underlying assumptions and their ability to be evaluated, nonparametric estimators are a promising tool for assessing microbial diversity. To further investigate their potential, we applied these techniques to four microbial data sets. In particular, we compared the use of nonparametric estimators with the rarefaction approach and investigated how the precision of their estimates changes with sample size. These four data sets were among the largest available and represented a range of habitat types and environmental gradients. We came across a number of additional data sets that would also have been appropriate for these analyses (19, 53), although others of comparable size were too diverse to be analyzed with these techniques (5, 45).

The analyses were performed with EstimateS (version 5.0.1; R. Colwell, University of Connecticut [http://viceroy.eeb.uconn.edu/estimates]). For the purposes of inputting data into the program, we treated each cloned sequence as a separate sample. We ran 100 randomizations for all tests. Further randomizations did not change the results.

Human mouth and gut. Two of the best-sampled microbial communities are from human habitats. Kroes et al. (33) sampled subgingival plaque from a human mouth. They used PCR to amplify the bacterial 16S rDNA, created clone libraries from the amplified DNA, and then sequenced 264 clones. Kroes et al. defined an OTU as a 16S rDNA sequence group in which sequences differed by <= 1%. By this definition, they found 59 distinct OTUs from their sample of 264 16S rDNA sequences. Although the accumulation curve does not reach an asymptote, it is not linear (Fig. 3). Thus, we can try to estimate total OTU richness. For these data, the Chao1 estimator levels off at 123 OTUs, suggesting that, after that point, the Chao1 estimate is relatively independent of sample size. In contrast, the ACE does not plateau as sample size increases, indicating that the estimate is not independent of sample size.


View larger version (21K):
[in this window]
[in a new window]
 
FIG. 3.   Observed and estimated OTU richness of bacteria in a human mouth (33) versus sample size. The number of OTUs observed for a given sample size, or the accumulation curve, is averaged over 50 simulations (open circle ). Estimated OTU richness is plotted for Chaol () and ACE (black-triangle) estimators.

Suau et al. (65) investigated the diversity of bacteria in a human gut. Similar to Kroes et al. (33), they amplified, cloned, and sequenced 16S rDNA fragments. Their definition of an OTU differed slightly from that in the Kroes et al. study, however; they define an OTU as a 16S rDNA sequence group in which sequences differed by <= 2%. With this definition, they identified 82 OTUs from 284 clones.

Because the two studies use slightly different definitions of an OTU, the data for the mouth and gut bacteria are not entirely comparable. Their contrast does demonstrate the application of these approaches, however. After an initial increase, the mean Chao1 estimate for both communities is relatively level as sample size increases, and therefore we can compare the estimates at the highest sample size for each community (Fig. 4). We used a log transformation to calculate the confidence intervals (CIs) because the distribution of estimates is not normal (8). Given the OTU definitions, total richness of the mouth and gut bacterial communities is not significantly different, as estimated by Chao1. Chao1 estimates that the mouth community has 123 OTUs (95% CIs, 93 and 180), and the gut community has 135 OTUs (95% CIs, 110 and 170).


View larger version (50K):
[in this window]
[in a new window]
 
FIG. 4.   Chaol estimates of human mouth (open circle ) and gut () bacterial richness as a function of sample size. Error bars are 95% CIs and were calculated with the variance formula derived by Chao (8). The dashed lines are error bars for the mouth. The solid lines are error bars for the gut.

What do the CIs say about the Chao1 estimate? The CIs estimate the precision of the richness estimates. In other words, 95% of new samples of 264 clones from the same person's mouth are predicted to yield Chao1 estimates that fall within this range. Because the CIs overlap, one cannot reject the null hypothesis at the significance level of 0.05 that there is no difference between the richness of the mouth and gut communities. The CIs do not address how close the estimates are to the true total richness (i.e., bias) or whether these samples are representative of other people's mouths or guts.

Another question is how much more sampling is needed to detect a significant difference between two estimates, which in this case differ by only 12 OTUs. The range of the CIs initially increases with sample size, peaks, and then decreases exponentially. To obtain a rough idea of how much further sampling would be needed to detect a statistically significant difference, we estimated the size of the CIs for larger samples by extrapolating from the decreasing portion of these curves. Negative exponential curves for both the mouth [f(x) = 270e-0.0046x] and gut [f(x) = 120e-0.0026x] data fit well (r2 = 0.90 and r2 = 0.87, respectively). From these curves, it appears that a sample of about 1,000 clones (four times the original number) would be needed to detect a significant difference between these communities (Fig. 5).


View larger version (22K):
[in this window]
[in a new window]
 
FIG. 5.   Average size of the 95% CIs of Chaol estimates for bacteria in the human mouth (open circle ) and gut () as sample size increases. These CIs are the same as in Fig. 4, but only the decreasing portions of the CIs are plotted. The curves are fitted negative exponential curves [mouth, f(x) = 270e-0.0046x, r2 = 0.90; gut, f(x) = 120e-0.0026x, r2 = 0.87].

Rarefaction curves yield the same pattern of relative diversity as Chao1; significantly more OTUs are observed in the gut sample than the mouth sample (Fig. 6). At the highest shared sample size (264 clones), 79 OTUs are observed in the gut versus 59 OTUs in the mouth, and the 95% CIs do not overlap. As discussed in the previous section, however, rarefaction curves do not address the precision of the observed species richness. Thus, although the rarefaction curves suggest that the gut community is more diverse than the mouth community, we cannot address the statistical significance of this evidence with rarefaction curves.


View larger version (32K):
[in this window]
[in a new window]
 
FIG. 6.   Rarefaction curves of observed OTU richness in human mouth (open circle ) and gut () bacterial samples. The error bars are 95% CIs and were calculated from the variance of the number of OTUs drawn in 100 randomizations at each sample size.

Aquatic mesocosms. Bohannan and Leibold (unpublished data) sampled bacterial diversity from three outdoor aquatic mesocosms designed to mimic small ponds. The mesocosms varied along a gradient of increasing primary productivity and decreasing eukaryotic algal diversity, and all received the same inoculum. DNA was extracted from samples from each mesocosm, and a region of 16S rDNA was PCR amplified with Bacteria-specific primers, the amplicons were cloned, and the clones were sequenced. The sequences were grouped into OTUs using a definition of 95% similarity.

Bohannan and Leibold sequenced 158, 128, and 174 clones from the low-, intermediate-, and high-productivity mesocosms, respectively. The Chao1 estimates suggest that OTU richness varies positively with productivity. The lowest productivity pond contained 54 OTUs (95% CIs, 42 and 80), the intermediate pond contained 58 OTUs (43 and 90), and the high-productivity pond contained an estimated 95 OTUs (73 and 140). The richness of the high- and low-productivity ponds is significantly different at the 0.10 level (Fig. 7). Furthermore, the Chao1 estimates for the high-productivity pond have not yet stabilized (Fig. 7), suggesting that further sampling will result in a greater difference in richness between the ponds with low and high productivity.


View larger version (38K):
[in this window]
[in a new window]
 
FIG. 7.   Chaol estimates of bacterial OTU richness in low- (), intermediate- (), and high- (black-triangle) productivity ponds. Error bars are 90% CIs and were calculated with the variance formula derived by Chao (8). The dotted, solid, and dashed bars are error bars for the low-, intermediate-, and high-productivity mesocosms, respectively.

Scottish soil. The most diverse data set that we analyzed is for terrestrial soil. McCaig et al. (39) collected soil samples from two grazed grasslands, allowing us to make a direct comparison of microbial diversity between these two habitats. One grassland was previously reseeded and fertilized (improved), and the other was not (unimproved). As in the studies described above, bacterial 16S rDNA was PCR amplified and cloned.

McCaig et al. sequenced 137 clones from the improved soil and 138 clones from the unimproved soil. By their OTU definition of <3% sequence difference, they identified 113 OTUs in the improved habitat and 117 in the unimproved habitat. The Chao1 estimates level off in both habitats at about 70 clones. Bacterial richness appears to be higher in the unimproved habitat (590 OTUs) than in the improved habitat (467 OTU), but the difference is not significant (Fig. 8). As before, we can approximate how much further sampling is needed to detect a significant difference by extrapolating the range of the CIs at larger sample sizes. Negative exponential curves fit very well for the improved [f(x) = 1,500e-0.012x, r2 = 0.96] and unimproved [f(x) = 2,000e-0.011x, r2 = 0.94] soil samples. Thus, if these estimates remain stable with more sampling, about 250 clones are needed to detect a significant difference at the 0.05 level (Fig. 9).


View larger version (44K):
[in this window]
[in a new window]
 
FIG. 8.   Chaol estimates of bacterial OTU richness in improved (open circle ) and unimproved () soil as a function of sample size. Error bars are 95% CIs and were calculated with the variance formula derived by Chao (8). The solid lines are error bars for the improved sample. The dashed lines are error bars for the unimproved sample.


View larger version (20K):
[in this window]
[in a new window]
 
FIG. 9.   Average size of the 95% CIs of Chaol estimates for the improved (open circle ) and unimproved () soil as the number of clones sampled increases. These CIs are the same as in Fig. 8, but only the decreasing portions of the CIs are plotted. The curves are fitted negative exponential curves [improved, f(x) = 1,500e-0.012x, r2 = 0.96; unimproved, f(x) = 2,000e-0.011x, R2 = 0.94].


    DISCUSSION
Top
Introduction
Discussion
References

Comparisons of accumulation curves and rank-abundance plots demonstrate that some bacterial communities have been sampled as well as some macroorganism communities. Therefore, evaluating microbial diversity with statistical approaches available for macroorganisms seems feasible. We estimated and compared microbial richness in a variety of habitats and found that although the estimators depend on sample size, most of the richness estimates stabilized with the sample sizes available. We also made rough estimates of the sample sizes needed to detect significant differences in diversity between comparable samples.

Of course, these statistical approaches have their limitations. For example, diversity comparisons require clear OTU definitions. Often microbial "species" are defined by a cutoff of percent genetic similarity, leading some authors to charge that microbial diversity studies adopt arbitrary species definitions (62). This problem is not limited to microorganisms, however. In fact, the debate over species definitions in eukaryotic organisms has persisted for decades (16, 18, 37, 38), and some suggest that even in sexual organisms, "the prevalence of the clearly defined species is a myth" (21).

Similarly, most of these approaches require data on the relative frequencies of different OTUs, and many studies have revealed that sampling biases accompany genetic surveys of microbial diversity. For example, the abundances of amplified genes in PCRs may not reflect the relative abundances of template DNA because of differences in primer binding and elongation efficiency (52, 55, 67). Larger organisms differ in their ease of detection as well, and hence samples may not be representative of the species frequencies in a community. For example, butterfly species differ in their attraction to bait traps (29), and bird species' vocalizations are unequally detectable (58).

The fact that most questions about the structure and function of communities require relative comparisons overcomes many of the problems with species definitions and sampling biases. As long as the measurement unit is defined and held constant, diversity can be compared among sites or treatments. Likewise, to minimize the effect of sampling biases, multiple techniques or genes can be employed to increase the robustness of relative comparisons (44).

Further work is needed to investigate the general applicability of these approaches for microbial diversity studies. Ideally, large data sets should be gathered to evaluate better the bias and precision of different nonparametric estimators, such as Chao1 and ACE. The performance of richness estimators should also be measured in terms of their ability to predict the true ordering of richness among samples. Large data sets are also needed to investigate how often microbial accumulation curves cross with additional sampling. If the accumulation curves cross only infrequently, then, in combination with methods such as bootstrapping (20), rarefaction curves may be a valuable way to compare the relative diversity of communities.

Even without exhaustive surveys of microbial communities, computer simulations may provide useful insights. Simulated communities have already been used to compare the bias and precision of some diversity estimators (3, 27, 31, 68, 71). These studies could be extended to examine the ability of different estimators to predict the correct order of richness among samples and the conditions under which rarefaction curves are likely to cross. Of course, simulation studies cannot be used as a substitute for real data, as they require input on realistic species abundance distributions of microbial communities.

Although our discussion has been directed towards data collected from clone libraries, genetic techniques that do not depend on cloning also offer promising opportunities for quickly analyzing community diversity. For instance, denaturing gradient gel electrophoresis (DGGE) patterns of amplified 16s rDNA have been used as estimates of microbial diversity (42, 44). Incidence-based nonparametric estimators (R. Colwell, User's Guide to EstimateS 5 [http://viceroy.eeb.uconn.edu/estimates]), such as the jackknife and bootstrap (60, 70), use presence-absence data and could be used with DGGE data to estimate total richness. Likewise, oligonucleotide probes can be used to detect the presence of a subset of microbial diversity in a sample (28). Once the specific probes have been developed, many samples can be analyzed relatively quickly, and incidence estimators could be adapted to extrapolate these patterns to the entire community.

In conclusion, while microbiologists should be cautious about sampling biases and use clear OTU definitions, our results suggest that comparisons among estimates of microbial diversity are possible. Nonparametric estimators show particular promise for microbial data and in some habitats may require sample sizes of only 200 to 1,000 clones to detect richness differences of only tens of species. While daunting less than a decade ago, sequencing this number of clones is reasonable with the development of high-throughput sequencing technology. Augmenting this new technology with statistical approaches borrowed from "macrobial" biologists offers a powerful means to study the ecology and evolution of microbial diversity in natural environments.


    ACKNOWLEDGMENTS

We thank Ian Kroes, Paul Lepp, and David Relman; Allison McCaig, Jim Prosser, and the Scottish Executive Rural Affairs Department; and Antonia Suau, Joël Doré, and coworkers for sharing unpublished data. We also thank Robert Colwell, Craig Criddle, Gregory Gilbert, and Aaron Hirsh for comments on earlier drafts and Mark Tanaka, Lauren Ancel, and Michael Lachmann for useful discussions. B.B. is especially grateful to Dan Janzen and the supporters of the NSF/CRUSA workshop on microbial biocomplexity, at which the idea for this paper originated.

This work was supported by a National Science Foundation award (DEB-9907797) to B.B.


    FOOTNOTES

* Corresponding author. Present address: Dept. of Ecology and Evolutionary Biology and Center for Environmental Studies, Brown University, Providence, RI 02912. Phone: (401) 863-9676. Fax: (401) 863-2166. E-mail: Jennifer_Hughes{at}Brown.edu.

dagger Present address: Centre for Biodiversity Research, University of British Columbia, British Columbia, Canada.


    REFERENCES
Top
Introduction
Discussion
References

1. Achenbach, L. A., and J. D. Coates. 2000. Disparity between bacterial phylogeny and physiology. ASM News 66:714-715.
2. Amann, R. I., W. Ludwig, and K.-H. Schleifer. 1995. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol. Rev. 59:143-169[Abstract/Free Full Text].
3. Baltanas, A. 1992. On the use of some methods for the estimation of species richness. Oikos 65:484-492.
4. Bills, G. F., and J. D. Polishook. 1994. Abundance and diversity of microfungi in leaf-litter of a lowland rain-forest in Costa Rica. Mycologia 86:187-198.
5. Borneman, J., P. W. Skroch, K. M. O'Sullivan, J. A. Palus, N. G. Rumjanek, J. L. Jansen, J. Nienhuis, and E. W. Triplett. 1996. Molecular microbial diversity of an agricultural soil in Wisconsin. Appl. Environ. Microbiol. 62:1935-1943[Abstract].
6. Borneman, J., and E. W. Triplett. 1997. Molecular microbial diversity in soils from Eastern Amazonia: evidence for unusual microorganisms and microbial population shifts associated with deforestation. Appl. Environ. Microbiol. 63:2647-2653[Abstract].
7. Bulmer, M. G. 1974. On fitting the Poisson lognormal distribution to species abundance data. Biometrics 30:101-110[CrossRef].
8. Chao, A. 1987. Estimating the population size for capture-recapture data with unequal catchability. Biometrics 43:783-791[CrossRef][Medline].
9. Chao, A. 1984. Non-parametric estimation of the number of classes in a population. Scand. J. Stat. 11:265-270.
10. Chao, A., and S.-M. Lee. 1992. Estimating the number of classes via sample coverage. J. Am. Stat. Assoc. 87:210-217[CrossRef].
11. Chao, A., M.-C. Ma, and M. C. K. Yang. 1993. Stopping rules and estimation for recapture debugging with unequal failure rates. Biometrics 43:783-791.
12. Chazdon, R. L., R. K. Colwell, J. S. Denslow, and M. R. Guariguata. 1998. Statistical methods for estimating species richness of woody regeneration in primary and secondary rain forests of northeastern Costa Rica, p. 285-309. In F. Dallmeier, and J. A. Comiskey (ed.), Forest biodiversity reseach, monitoring, and modeling: conceptual background and old World case studies. Parthenon, Paris, France.
13. Clench, H. 1979. How to make regional lists of butterflies: Some thoughts. J. Lepid. Soc. 33:216-231.
14. Colwell, R. K., and J. A. Coddington. 1994. Estimating terrestrial biodiversity through extrapolation. Phil. Trans. R. Soc. London B 345:101-118[Medline].
15. Condit, R., S. P. Hubbell, J. V. Lafrankie, R. Sukumar, N. Manokaran, R. B. Foster, and P. S. Asthon. 1996. Species-area and species-individual relationships for tropical trees: a comparison of three 50-ha plots. J. Ecol. 84:549-562[CrossRef].
16. Coyne, J. A., H. A. Orr, and D. J. Futuyma. 1988. Do we need a new species concept? Syst. Zool. 37:190-200[CrossRef].
17. DeVries, P. J., D. Murray, and R. Lande. 1997. Species diversity in vertical, horizontal, and temporal dimensions of a fruit-feeding butterfly community in an Ecuadorian rainforest. Biol. J. Linnean Soc. 62:343-364[CrossRef].
18. Dobzhansky, T. 1935. A critique of the species concept in biology. Phil. Sci. 2:344-355.
19. Dunbar, J., S. Takala, S. M. Barnes, J. A. Davis, and C. R. Kuske. 1999. Levels of bacterial community diversity in four arid soils compared by cultivation and 16S rRNA gene cloning. Appl. Environ. Microbiol. 65:1662-1669[Abstract/Free Full Text].
20. Efron, B., and R. Tibshirani. 1993. An introduction to the bootstrap. Chapman & Hall, New York, N.Y.
21. Ehrlich, P. R. 1961. Has the biological species concept outlived its usefulness? Syst. Zool. 10:167-176[CrossRef].
22. Erwin, T. L. 1982. Tropical forests: their richness in Coleoptera and other arthropod species. Coleopt. Bull. 36:74-82.
23. Gaston, K. J. 1996. Species richness: measure and measurement, p. 77-113. In K. J. Gaston (ed.), Biodiversity: a biology of numbers and difference. Blackwell, Cambridge, England.
24. Harte, J., A. Kinzig, and J. Green. 1999. Self-similarity in the distribution and abundance of species. Science 284:334-336[Abstract/Free Full Text].
25. Heck, K. L., G. v. Belle, and D. Simberloff. 1975. Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size. Ecology 56:1459-1461[CrossRef].
26. Hellmann, J. J., and G. W. Fowler. 1999. Bias, precision, and accuracy of four measures of species richness. Ecol. Appl. 9:824-834.
27. Heltshe, J. F., and N. E. Forrester. 1983. Estimating species richness using the jackknife procedure. Biometrics 39:1-11[CrossRef][Medline].
28. Heuer, H., K. Hartung, G. Wieland, I. Kramer, and K. Smalla. 1999. Polynicleotide probes that target a hypervariable region of 16S rRNA genes to identify bacterial isolates corresponding to bands of community fingerprints. Appl. Environ. Microbiol. 65:1045-1049[Abstract/Free Full Text].
29. Hughes, J. B., G. C. Daily, and P. R. Ehrlich. 1998. Use of fruit bait traps for monitoring of butterflies (Lepidoptera: Nymphalidae). Rev. Biol. Trop. 46:697-704.
30. Hughes, R. G. 1986. Theories and models of species abundance. Am. Nat. 128:897-899[CrossRef].
31. Keating, K. A., J. F. Quinn, M. A. Ivie, and L. L. Ivie. 1998. Estimating the effectiveness of further sampling in species inventories. Ecol. Appl. 8:1239-1249.
32. Krebs, C. J. 1989. Ecological methodology. Harper and Row, New York, N.Y.
33. Kroes, I., P. W. Lepp, and D. A. Relman. 1999. Bacterial diversity within the human subgingival crevice. Proc. Natl. Acad. Sci. USA 96:14547-14552[Abstract/Free Full Text].
34. Lande, R., P. J. DeVries, and T. R. Walla. 2000. When species accumulation curves intersect: implications for ranking diversity using small samples. Oikos 89:601-605[CrossRef].
35. Liu, W. T., T. L. Marsh, H. CHeng, and L. J. Forney. 1997. Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl. Environ. Microbiol. 63:4516-4522[Abstract].
36. Magurran, A. E. 1995. Ecological diversity and its measurement. Princeton University, Princeton, N.J.
37. Masters, J. C., and H. G. Spencer. 1989. Why we need a new genetic species concept. Syst. Zool. 30:270-279[CrossRef].
38. Mayr, E. 1940. Systematics and the origin of species. Columbia University Press, New York, N.Y.
39. McCaig, A. E., L. Glover, and J. I. Prosser. 1999. Molecular analysis of bacterial community structure and diversity in unimproved and improved upland grass pastures. Appl. Environ. Microbiol. 65:1721-1730[Abstract/Free Full Text].
40. Moyer, C. L., J. M. Tiedje, F. C. Dobbs, and D. M. Karl. 1998. Diversity of deep-sea hydrothermal vent Archaea from Loihi Seamount, Hawaii. Deep-Sea Res. II 45:303-317[CrossRef].
41. Mullins, T. D., T. B. Britschgi, R. L. Krest, and S. J. Giovannoni. 1995. Genetic comparisons reveal the same unknown bacterial lineages in Atlantic and Pacific bacterioplankton communities. Limnol. Oceanogr. 40:148-158.
42. Muyzer, G. 1999. DGGE/TGGE a method for identifying genes from natural ecosystems. Curr. Opin. Microbiol. 2:317-322[CrossRef][Medline].
43. Nübel, U., F. Garcia-Pichel, M. Kühl, and G. Muyzer. 1999. Spatial scale and the diversity of benthic cyanobacteria and diatoms in a salina. Hydrobiologia 401:199-206[CrossRef].
44. Nübel, U., F. Gracia-Pichel, M. Kuhl, and G. Muyzer. 1999. Quantifying microbial diversity: morphotypes, 16S rRNA genes, and carotenoids of oxygenic phototrophs in microbial mats. App. Environ. Microbiol. 65:422-430[Abstract/Free Full Text].
45. Nusslein, K., and J. M. Tiedje. 1998. Characterization of the dominant and rare members of a young Hawaiian soil bacterial community with small-subunit ribosomal DNA amplified from DNA fractionated on the basis of its guanine and cytosine composition. Appl. Environ. Microbiol. 64:1283-1289[Abstract/Free Full Text].
46. Pace, N. R., D. A. Stahl, D. J. Lane, and G. J. Olsen. 1986. The analysis of natural microbial populations by ribosomal RNA sequences. Adv. Microb. Ecol. 9:1-55.
47. Palmer, M. W. 1990. The estimation of species richness by extrapolation. Ecology. 71:1195-1198[CrossRef].
48. Pielou, E. C. 1975. Ecological diversity. Wiley, New York, N.Y.
49. Preston, F. W. 1948. The commonness, and rarity, of species. Ecology 29:254-283[CrossRef].
50. Price, P. W. 1975. Insect ecology. Wiley, New York, N.Y.
51. Raaijmakers, J. G. W. 1987. Statistical analysis of the Michaelis-Menten equation. In Biometrics 43:793-803.
52. Raeymaekers, L. 1995. A commentary on the practical applications of quantitative PCR. Genome Res. 5:91-94[Free Full Text].
53. Rappé, M. S., P. F. Kemp, and S. J. Giovannoni. 1997. Phylogenetic diversity of marine coastal picoplankton 16S rRNA genes cloned from the continental shelf off Cape Hatteras, North Carolina. Limnol. Oceanogr. 42:811-826.
54. Raven, P. H. 1983. The challenge of tropical biology. Bull. ESA Spring:4-12.
55. Reysenbach, A.-L., L. J. Giver, G. S. Wickham, and N. R. Pace. 1992. Differential amplification of rRNA genes by polymerase chain reaction. Appl. Environ. Microbiol. 58:3417-3418[Abstract/Free Full Text].
56. Ricketts, T. H., G. C. Daily, P. R. Ehrlich, and J. P. Fay. 2001. Countryside biogeography of moths in native and human-dominated habitats. Conserv. Biol. 15:378-388[CrossRef].
57. Samways, M. J. 1994. Insect conservation biology. Chapman & Hall, London, England.
58. Schieck, J. 1997. Biased detection of bird vocalizations affects comparisons of bird abundance among forested habitats. Condor 99:179-190.
59. Seber, G. A. F. 1973. The estimation of animal abundance and related parameters. Griffin, London, England.
60. Smith, E. P., and G. Van Belle. 1984. Nonparametric estimation of species richness. Biometrics 40:119-129[CrossRef].
61. Soberón, J., and J. Llorente. 1993. The use of species accumulation functions for the prediction of species richness. Conserv. Biol. 7:480-488[CrossRef].
62. Staley, J. T. 1997. Biodiversity: are microbial species threatened? Curr. Opin. Biotechnol. 8:340-345[CrossRef][Medline].
63. Stork, N. E., J. Adis, and R. K. Didham. 1997. Canopy arthropods. Chapman & Hall, London, England.
64. Strong, D. R., J. H. Lawton, and R. Southwood. 1984. Insects on plants: community patterns and mechanisms. Harvard University Press, Cambridge, Mass.
65. Suau, A., R. Bonnet, M. Sutren, J.-J. Godon, G. Gibson, M. D. Collins, and J. Dore. 1999. Direct analysis of genes encoding 16S rRNA from complex communities reveals many novel molecular species within the human gut. Appl. Environ. Microbiol. 65:4799-4807[Abstract/Free Full Text].
66. Sugihara, G. 1980. Minimal community structure: an explanation of species abundance patterns. Am. Nat. 116:770-787[CrossRef].
67. Suzuki, M. T., and S. J. Giovannoni. 1996. Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR. Appl. Environ. Microbiol. 62:625-630[Abstract].
68. Walther, B. A., and S. Morand. 1998. Comparative performance of species richness estimation methods. Parasitology 116:395-405.
69. Watve, M. G., and R. M. Gangal. 1996. Problems in measuring bacterial diversity and a possible solution. Appl. Environ. Microbiol. 62:4299-4301[Abstract].
70. Zahl, S. 1977. Jackknifing an index of diversity. Ecology 58:907-913[CrossRef].
71. Zelmer, D. A., and G. W. Esch. 1999. Robust estimation of parasite component community richness. J. Parasitol. 85:592-594[CrossRef][Medline].


Applied and Environmental Microbiology, October 2001, p. 4399-4406, Vol. 67, No. 10
0099-2240/01/$04.00+0   DOI: 10.1128/AEM.67.10.4399-4406.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.



This article has been cited by other articles:

  • Herrera, C. M., Canto, A., Pozo, M. I., Bazaga, P. (2009). Inhospitable sweetness: nectar filtering of pollinator-borne inocula leads to impoverished, phylogenetically clustered yeast communities. Proc R Soc B 0: rspb.2009.1485v1-rspb20091485 [Abstract] [Full Text]  
  • Abnous, K., Brooks, S. P.J., Kwan, J., Matias, F., Green-Johnson, J., Selinger, L. B., Thomas, M., Kalmokoff, M. (2009). Diets Enriched in Oat Bran or Wheat Bran Temporally and Differentially Alter the Composition of the Fecal Community of Rats. J. Nutr. 139: 2024-2031 [Abstract] [Full Text]  
  • Zinger, L., Coissac, E., Choler, P., Geremia, R. A. (2009). Assessment of Microbial Communities by Graph Partitioning in a Study of Soil Fungi in Two Alpine Meadows. Appl. Environ. Microbiol. 75: 5863-5870 [Abstract] [Full Text]  
  • Youssef, N., Sheik, C. S., Krumholz, L. R., Najar, F. Z., Roe, B. A., Elshahed, M. S. (2009). Comparison of Species Richness Estimates Obtained Using Nearly Complete Fragments and Simulated Pyrosequencing-Generated Fragments in 16S rRNA Gene-Based Environmental Surveys. Appl. Environ. Microbiol. 75: 5227-5236 [Abstract] [Full Text]  
  • Valster, R. M., Wullings, B. A., Bakker, G., Smidt, H., van der Kooij, D. (2009). Free-Living Protozoa in Two Unchlorinated Drinking Water Supplies, Identified by Phylogenic Analysis of 18S rRNA Gene Sequences. Appl. Environ. Microbiol. 75: 4736-4746 [Abstract] [Full Text]  
  • Schellenberg, J., Links, M. G., Hill, J. E., Dumonceaux, T. J., Peters, G. A., Tyler, S., Ball, T. B., Severini, A., Plummer, F. A. (2009). Pyrosequencing of the Chaperonin-60 Universal Target as a Tool for Determining Microbial Community Composition. Appl. Environ. Microbiol. 75: 2889-2898 [Abstract] [Full Text]  
  • Biers, E. J., Sun, S., Howard, E. C. (2009). Prokaryotic Genomes and Diversity in Surface Ocean Waters: Interrogating the Global Ocean Sampling Metagenome. Appl. Environ. Microbiol. 75: 2221-2229 [Abstract] [Full Text]  
  • Tanriverdi, S., Grinberg, A., Chalmers, R. M., Hunter, P. R., Petrovic, Z., Akiyoshi, D. E., London, E., Zhang, L., Tzipori, S., Tumwine, J. K., Widmer, G. (2008). Inferences about the Global Population Structures of Cryptosporidium parvum and Cryptosporidium hominis. Appl. Environ. Microbiol. 74: 7227-7234 [Abstract] [Full Text]  
  • Steinberg, L. M., Regan, J. M. (2008). Phylogenetic Comparison of the Methanogenic Communities from an Acidic, Oligotrophic Fen and an Anaerobic Digester Treating Municipal Wastewater Sludge. Appl. Environ. Microbiol. 74: 6663-6671 [Abstract] [Full Text]  
  • Oliver, K. L., Hamelin, R. C., Hintz, W. E. (2008). Effects of Transgenic Hybrid Aspen Overexpressing Polyphenol Oxidase on Rhizosphere Diversity. Appl. Environ. Microbiol. 74: 5340-5348 [Abstract] [Full Text]  
  • Bryant, J. A., Lamanna, C., Morlon, H., Kerkhoff, A. J., Enquist, B. J., Green, J. L. (2008). Colloquium Paper: Microbes on mountainsides: Contrasting elevational patterns of bacterial and plant diversity. Proc. Natl. Acad. Sci. USA 105: 11505-11511 [Abstract] [Full Text]  
  • Mendez, M. O., Neilson, J. W., Maier, R. M. (2008). Characterization of a Bacterial Community in an Abandoned Semiarid Lead-Zinc Mine Tailing Site. Appl. Environ. Microbiol. 74: 3899-3907 [Abstract] [Full Text]  
  • Lepere, C., Domaizon, I., Debroas, D. (2008). Unexpected Importance of Potential Parasites in the Composition of the Freshwater Small-Eukaryote Community. Appl. Environ. Microbiol. 74: 2940-2949 [Abstract] [Full Text]  
  • Brons, J. K., van Elsas, J. D. (2008). Analysis of Bacterial Communities in Soil by Use of Denaturing Gradient Gel Electrophoresis and Clone Libraries, as Influenced by Different Reverse Primers. Appl. Environ. Microbiol. 74: 2717-2727 [Abstract] [Full Text]  
  • Qvit-Raz, N., Jurkevitch, E., Belkin, S. (2008). Drop-Size Soda Lakes: Transient Microbial Habitats on a Salt-Secreting Desert Tree. Genetics 178: 1615-1622 [Abstract] [Full Text]  
  • Fierer, N., Breitbart, M., Nulton, J., Salamon, P., Lozupone, C., Jones, R., Robeson, M., Edwards, R. A., Felts, B., Rayhawk, S., Knight, R., Rohwer, F., Jackson, R. B. (2007). Metagenomic and Small-Subunit rRNA Analyses Reveal the Genetic Diversity of Bacteria, Archaea, Fungi, and Viruses in Soil. Appl. Environ. Microbiol. 73: 7059-7066 [Abstract] [Full Text]  
  • Duriez, P., Topp, E. (2007). Temporal Dynamics and Impact of Manure Storage on Antibiotic Resistance Patterns and Population Structure of Escherichia coli Isolates from a Commercial Swine Farm. Appl. Environ. Microbiol. 73: 5486-5493 [Abstract] [Full Text]  
  • Kim, J.-S., Crowley, D. E. (2007). Microbial Diversity in Natural Asphalts of the Rancho La Brea Tar Pits. Appl. Environ. Microbiol. 73: 4579-4591 [Abstract] [Full Text]  
  • Ashby, M. N., Rine, J., Mongodin, E. F., Nelson, K. E., Dimster-Denk, D. (2007). Serial Analysis of rRNA Genes and the Unexpected Dominance of Rare Members of Microbial Communities. Appl. Environ. Microbiol. 73: 4532-4542 [Abstract] [Full Text]  
  • Buckley, D. H., Huangyutitham, V., Hsu, S.-F., Nelson, T. A. (2007). Stable Isotope Probing with 15N2 Reveals Novel Noncultivated Diazotrophs in Soil. Appl. Environ. Microbiol. 73: 3196-3204 [Abstract] [Full Text]  
  • Mathur, J., Bizzoco, R. W., Ellis, D. G., Lipson, D. A., Poole, A. W., Levine, R., Kelley, S. T. (2007). Effects of Abiotic Factors on the Phylogenetic Diversity of Bacterial Communities in Acidic Thermal Springs. Appl. Environ. Microbiol. 73: 2612-2623 [Abstract] [Full Text]  
  • Bent, S. J., Pierson, J. D., Forney, L. J., Danovaro, R., Luna, G. M., Dell'Anno, A., Pietrangeli, B. (2007). Measuring Species Richness Based on Microbial Community Fingerprints: the Emperor Has No Clothes. Appl. Environ. Microbiol. 73: 2399-2401 [Full Text]  
  • Poulain, A. J., Ni Chadhain, S. M., Ariya, P. A., Amyot, M., Garcia, E., Campbell, P. G. C., Zylstra, G. J., Barkay, T. (2007). Potential for Mercury Reduction by Microbes in the High Arctic. Appl. Environ. Microbiol. 73: 2230-2238 [Abstract] [Full Text]  
  • Lehours, A.-C., Evans, P., Bardot, C., Joblin, K., Gerard, F. (2007). Phylogenetic Diversity of Archaea and Bacteria in the Anoxic Zone of a Meromictic Lake (Lake Pavin, France). Appl. Environ. Microbiol. 73: 2016-2019 [Abstract] [Full Text]  
  • Mills, D. K., Entry, J. A., Gillevet, P. M., Mathee, K. (2007). Assessing Microbial Community Diversity Using Amplicon Length Heterogeneity Polymerase Chain Reaction. Soil Sci. 71: 572-578 [Abstract] [Full Text]  
  • Zwolinski, M. D. (2007). DNA Sequencing: Strategies for Soil Microbiology. Soil Sci. 71: 592-600 [Abstract] [Full Text]  
  • Gao, Z., Tseng, C.-h., Pei, Z., Blaser, M. J. (2007). Molecular analysis of human forearm superficial skin bacterial biota. Proc. Natl. Acad. Sci. USA 104: 2927-2932 [Abstract] [Full Text]  
  • Hartmann, M., Widmer, F. (2006). Community Structure Analyses Are More Sensitive to Differences in Soil Bacterial Communities than Anonymous Diversity Indices. Appl. Environ. Microbiol. 72: 7804-7812 [Abstract] [Full Text]  
  • Ruan, W., Sassoon, A., An, F., Simko, J. P., Liu, B. (2006). Identification of Clinically Significant Tumor Antigens by Selecting Phage Antibody Library on Tumor Cells in Situ Using Laser Capture Microdissection. Mol. Cell. Proteomics 5: 2364-2373 [Abstract] [Full Text]  
  • Cohan, F. M (2006). Towards a conceptual and operational union of bacterial systematics, ecology, and evolution. Phil Trans R Soc B 361: 1985-1996 [Abstract] [Full Text]  
  • Lynch, M. D. J., Thorn, R. G. (2006). Diversity of Basidiomycetes in Michigan Agricultural Soils. Appl. Environ. Microbiol. 72: 7050-7056 [Abstract] [Full Text]  
  • Jeon, S.-O., Bunge, J., Stoeck, T., Barger, K. J.-A., Hong, S.-H., Epstein, S. S. (2006). Synthetic statistical approach reveals a high degree of richness of microbial eukaryotes in an anoxic water column.. Appl. Environ. Microbiol. 72: 6578-6583 [Abstract] [Full Text]  
  • Schloss, P. D., Handelsman, J. (2006). Introducing SONS, a Tool for Operational Taxonomic Unit-Based Comparisons of Microbial Community Memberships and Structures.. Appl. Environ. Microbiol. 72: 6773-6779 [Abstract] [Full Text]  
  • Buckley, D. H., Huangyutitham, V., Nelson, T. A., Rumberger, A., Thies, J. E. (2006). Diversity of planctomycetes in soil in relation to soil history and environmental heterogeneity.. Appl. Environ. Microbiol. 72: 4522-4531 [Abstract] [Full Text]  
  • DeSantis, T. Z. Jr, Hugenholtz, P., Keller, K., Brodie, E. L., Larsen, N., Piceno, Y. M., Phan, R., Andersen, G. L. (2006). NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes.. Nucleic Acids Res 34: W394-W399 [Abstract] [Full Text]  
  • Schloss, P. D., Handelsman, J. (2006). Introducing TreeClimber, a Test To Compare Microbial Community Structures. Appl. Environ. Microbiol. 72: 2379-2384 [Abstract] [Full Text]  
  • Neubert, K., Mendgen, K., Brinkmann, H., Wirsel, S. G. R. (2006). Only a Few Fungal Species Dominate Highly Diverse Mycofloras Associated with the Common Reed. Appl. Environ. Microbiol. 72: 1118-1128 [Abstract] [Full Text]  
  • Hong, S.-H., Bunge, J., Jeon, S.-O., Epstein, S. S. (2006). Predicting microbial species richness. Proc. Natl. Acad. Sci. USA 103: 117-122 [Abstract] [Full Text]  
  • Lozupone, C., Knight, R. (2005). UniFrac: a New Phylogenetic Method for Comparing Microbial Communities. Appl. Environ. Microbiol. 71: 8228-8235 [Abstract] [Full Text]  
  • Acinas, S. G., Sarma-Rupavtarm, R., Klepac-Ceraj, V., Polz, M. F. (2005). PCR-Induced Sequence Artifacts and Bias: Insights from Comparison of Two 16S rRNA Clone Libraries Constructed from the Same Sample. Appl. Environ. Microbiol. 71: 8966-8969 [Abstract] [Full Text]  
  • Lewis, L. A., Lewis, P. O. (2005). Unearthing the Molecular Phylodiversity of Desert Soil Green Algae (Chlorophyta). Syst Biol 54: 936-947 [Abstract] [Full Text]  
  • Skidmore, M., Anderson, S. P., Sharp, M., Foght, J., Lanoil, B. D. (2005). Comparison of Microbial Community Compositions of Two Subglacial Environments Reveals a Possible Role for Microbes in Chemical Weathering Processes. Appl. Environ. Microbiol. 71: 6986-6997 [Abstract] [Full Text]  
  • Neufeld, J. D., Mohn, W. W. (2005). Unexpectedly High Bacterial Diversity in Arctic Tundra Relative to Boreal Forest Soils, Revealed by Serial Analysis of Ribosomal Sequence Tags. Appl. Environ. Microbiol. 71: 5710-5718 [Abstract] [Full Text]  
  • Lu, Z., Lapen, D., Scott, A., Dang, A., Topp, E. (2005). Identifying Host Sources of Fecal Pollution: Diversity of Escherichia coli in Confined Dairy and Swine Production Systems. Appl. Environ. Microbiol. 71: 5992-5998 [Abstract] [Full Text]  
  • Timke, M., Wang-Lieu, N. Q., Altendorf, K., Lipski, A. (2005). Community Structure and Diversity of Biofilms from a Beer Bottling Plant as Revealed Using 16S rRNA Gene Clone Libraries. Appl. Environ. Microbiol. 71: 6446-6452 [Abstract] [Full Text]  
  • O'Brien, H. E., Parrent, J. L., Jackson, J. A., Moncalvo, J.-M., Vilgalys, R. (2005). Fungal Community Analysis by Large-Scale Sequencing of Environmental Samples. Appl. Environ. Microbiol. 71: 5544-5550 [Abstract] [Full Text]  
  • Papineau, D., Walker, J. J., Mojzsis, S. J., Pace, N. R. (2005). Composition and Structure of Microbial Communities from Stromatolites of Hamelin Pool in Shark Bay, Western Australia. Appl. Environ. Microbiol. 71: 4822-4832 [Abstract] [Full Text]  
  • Weitz, J. S., Hartman, H., Levin, S. A. (2005). Coevolutionary arms races between bacteria and bacteriophage. Proc. Natl. Acad. Sci. USA 102: 9535-9540 [Abstract] [Full Text]  
  • Floyd, M. M., Tang, J., Kane, M., Emerson, D. (2005). Captured Diversity in a Culture Collection: Case Study of the Geographic and Habitat Distributions of Environmental Isolates Held at the American Type Culture Collection. Appl. Environ. Microbiol. 71: 2813-2823 [Full Text]  
  • Dick, L. K., Simonich, M. T., Field, K. G. (2005). Microplate Subtractive Hybridization To Enrich for Bacteroidales Genetic Markers for Fecal Source Identification. Appl. Environ. Microbiol. 71: 3179-3183 [Abstract] [Full Text]  
  • Tringe, S. G., von Mering, C., Kobayashi, A., Salamov, A. A., Chen, K., Chang, H. W., Podar, M., Short, J. M., Mathur, E. J., Detter, J. C., Bork, P., Hugenholtz, P., Rubin, E. M. (2005). Comparative Metagenomics of Microbial Communities. Science 308: 554-557 [Abstract] [Full Text]  
  • Schloss, P. D., Handelsman, J. (2005). Introducing DOTUR, a Computer Program for Defining Operational Taxonomic Units and Estimating Species Richness. Appl. Environ. Microbiol. 71: 1501-1506 [Abstract] [Full Text]  
  • Thompson, J. R., Pacocha, S., Pharino, C., Klepac-Ceraj, V., Hunt, D. E., Benoit, J., Sarma-Rupavtarm, R., Distel, D. L., Polz, M. F. (2005). Genotypic Diversity Within a Natural Coastal Bacterioplankton Population. Science 307: 1311-1313 [Abstract] [Full Text]  
  • Schloss, P. D., Handelsman, J. (2004). Status of the Microbial Census. Microbiol. Mol. Biol. Rev. 68: 686-691 [Abstract] [Full Text]  
  • Schloss, P. D., Larget, B. R., Handelsman, J. (2004). Integration of Microbial Ecology and Statistics: a Test To Compare Gene Libraries. Appl. Environ. Microbiol. 70: 5485-5492 [Abstract] [Full Text]  
  • Hewson, I., Fuhrman, J. A. (2004). Richness and Diversity of Bacterioplankton Species along an Estuarine Gradient in Moreton Bay, Australia. Appl. Environ. Microbiol. 70: 3425-3433 [Abstract] [Full Text]  
  • Roling, W. F. M., Couto de Brito, I. R., Swannell, R. P. J., Head, I. M. (2004). Response of Archaeal Communities in Beach Sediments to Spilled Oil and Bioremediation. Appl. Environ. Microbiol. 70: 2614-2620 [Abstract] [Full Text]  
  • Acinas, S. G., Marcelino, L. A., Klepac-Ceraj, V., Polz, M. F. (2004). Divergence and Redundancy of 16S rRNA Sequences in Genomes with Multiple rrn Operons. J. Bacteriol. 186: 2629-2635 [Abstract] [Full Text]  
  • Holben, W. E., Feris, K. P., Kettunen, A., Apajalahti, J. H. A. (2004). GC Fractionation Enhances Microbial Community Diversity Assessment and Detection of Minority Populations of Bacteria by Denaturing Gradient Gel Electrophoresis. Appl. Environ. Microbiol. 70: 2263-2270 [Abstract] [Full Text]  
  • Seghers, D., Wittebolle, L., Top, E. M., Verstraete, W., Siciliano, S. D. (2004). Impact of Agricultural Practices on the Zea mays L. Endophytic Community. Appl. Environ. Microbiol. 70: 1475-1482 [Abstract] [Full Text]  
  • Yang, H.-H., Vinopal, R. T., Grasso, D., Smets, B. F. (2004). High Diversity among Environmental Escherichia coli Isolates from a Bovine Feedlot. Appl. Environ. Microbiol. 70: 1528-1536 [Abstract] [Full Text]  
  • Dorigo, U., Jacquet, S., Humbert, J.-F. (2004). Cyanophage Diversity, Inferred from g20 Gene Analyses, in the Largest Natural Lake in France, Lake Bourget. Appl. Environ. Microbiol. 70: 1017-1022 [Abstract] [Full Text]  
  • Mai, V., Morris, J. G. Jr. (2004). Colonic Bacterial Flora: Changing Understandings in the Molecular Age. J. Nutr. 134: 459-464 [Abstract] [Full Text]  
  • Stres, B., Mahne, I., Avgustin, G., Tiedje, J. M. (2004). Nitrous Oxide Reductase (nosZ) Gene Fragments Differ between Native and Cultivated Michigan Soils. Appl. Environ. Microbiol. 70: 301-309 [Abstract] [Full Text]  
  • Dilly, O., Bloem, J., Vos, A., Munch, J. C. (2004). Bacterial Diversity in Agricultural Soils during Litter Decomposition. Appl. Environ. Microbiol. 70: 468-474 [Abstract] [Full Text]  
  • Egert, M., Wagner, B., Lemke, T., Brune, A., Friedrich, M. W. (2003). Microbial Community Structure in Midgut and Hindgut of the Humus-Feeding Larva of Pachnoda ephippiata (Coleoptera: Scarabaeidae). Appl. Environ. Microbiol. 69: 6659-6668 [Abstract] [Full Text]  
  • Stach, J. E. M., Maldonado, L. A., Masson, D. G., Ward, A. C., Goodfellow, M., Bull, A. T. (2003). Statistical Approaches for Estimating Actinobacterial Diversity in Marine Sediments. Appl. Environ. Microbiol. 69: 6189-6200 [Abstract] [Full Text]  
  • Seurinck, S., Verstraete, W., Siciliano, S. D. (2003). Use of 16S-23S rRNA Intergenic Spacer Region PCR and Repetitive Extragenic Palindromic PCR Analyses of Escherichia coli Isolates To Identify Nonpoint Fecal Sources. Appl. Environ. Microbiol. 69: 4942-4950 [Abstract] [Full Text]  
  • Singh, R., Stine, O. C., Smith, D. L., Spitznagel, J. K. Jr., Labib, M. E., Williams, H. N. (2003). Microbial Diversity of Biofilms in Dental Unit Water Systems. Appl. Environ. Microbiol. 69: 3412-3420 [Abstract] [Full Text]  
  • Pynaert, K., Smets, B. F., Wyffels, S., Beheydt, D., Siciliano, S. D., Verstraete, W. (2003). Characterization of an Autotrophic Nitrogen-Removing Biofilm from a Highly Loaded Lab-Scale Rotating Biological Contactor. Appl. Environ. Microbiol. 69: 3626-3635 [Abstract] [Full Text]  
  • Bowman, J. P., McCuaig, R. D. (2003). Biodiversity, Community Structural Shifts, and Biogeography of Prokaryotes within Antarctic Continental Shelf Sediment. Appl. Environ. Microbiol. 69: 2463-2483 [Abstract] [Full Text]  
  • Polz, M. F., Bertilsson, S., Acinas, S. G., Hunt, D. (2003). A(r)Ray of Hope in Analysis of the Function and Diversity of Microbial Communities. Biol. Bull. 204: 196-199 [Abstract] [Full Text]  
  • Boon, N., Top, E. M., Verstraete, W., Siciliano, S. D. (2003). Bioaugmentation as a Tool To Protect the Structure and Function of an Activated-Sludge Microbial Community against a 3-Chloroaniline Shock Load. Appl. Environ. Microbiol. 69: 1511-1520 [Abstract] [Full Text]  
  • Morris, C. E., Bardin, M., Berge, O., Frey-Klett, P., Fromin, N., Girardin, H., Guinebretiere, M.-H., Lebaron, P., Thiery, J. M., Troussellier, M. (2002). Microbial Biodiversity: Approaches to Experimental Design and Hypothesis Testing in Primary Scientific Literature from 1975 to 1999. Microbiol. Mol. Biol. Rev. 66: 592-616 [Abstract] [Full Text]  
  • Breitbart, M., Salamon, P., Andresen, B., Mahaffy, J. M., Segall, A. M., Mead, D., Azam, F., Rohwer, F. (2002). Genomic analysis of uncultured marine viral communities. Proc. Natl. Acad. Sci. USA 99: 14250-14255 [Abstract] [Full Text]  
  • Ward, B. B. (2002). How many species of prokaryotes are there?. Proc. Natl. Acad. Sci. USA 99: 10234-10236 [Full Text]  
  • Martin, A. P. (2002). Phylogenetic Approaches for Describing and Comparing the Diversity of Microbial Communities. Appl. Environ. Microbiol. 68: 3673-3682 [Full Text]  
  • Bishop-Hurley, S. L., Mounter, S. A., Laskey, J., Morris, R. O., Elder, J., Roop, P., Rouse, C., Schmidt, F. J., English, J. T. (2002). Phage-Displayed Peptides as Developmental Agonists for Phytophthora capsici Zoospores. Appl. Environ. Microbiol. 68: 3315-3320 [Abstract] [Full Text]  
  • Dunbar, J., Barns, S. M., Ticknor, L. O., Kuske, C. R. (2002). Empirical and Theoretical Bacterial Diversity in Four Arizona Soils. Appl. Environ. Microbiol. 68: 3035-3045 [Abstract] [Full Text]  

This Article
Right arrow Full Text (PDF)
Right arrow An erratum has been published
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hughes, J. B.
Right arrow Articles by Bohannan, B. J. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hughes, J. B.
Right arrow Articles by Bohannan, B. J. M.
Agricola
Right arrow Articles by Hughes, J. B.
Right arrow Articles by Bohannan, B. J. M.