Applied and Environmental Microbiology, October 2001, p. 4399-4406, Vol. 67, No. 10
Counting the Uncountable: Statistical Approaches
to Estimating Microbial Diversity
Department of Biological Sciences, Stanford
University, Stanford, California 94305-5020
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).
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.
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

![]()
INTRODUCTION
Top
Introduction
Discussion
References
![]()
ARE MICROBES TOO DIVERSE TO COUNT?

View larger version (29K):
[in a new window]
FIG. 1.
Accumulation curves for Michigan plants (
;
n = 1,783) (26), Costa Rican birds (
;
n = 5,007) (J. B. Hughes, unpublished data), human
oral bacteria (
; 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.
|
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
|
The ACE (10) incorporate data from all species with fewer
than 10 individuals, rather than just singletons and doubletons. ACE
estimates species richness as
|
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

|
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:
|
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.
|
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).
|
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).
|
|
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.
|
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).
|
|
| |
DISCUSSION |
|---|
|
|
|---|
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.
Present address: Centre for Biodiversity Research, University of
British Columbia, British Columbia, Canada.
| |
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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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 |
| 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]. |
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»