pH Influences the Importance of Niche-Related and Neutral Processes in Lacustrine Bacterioplankton Assembly

pH is an important factor that shapes the structure of bacterial communities. However, we have very limited information about the patterns and processes by which overall bacterioplankton communities assemble across wide pH gradients in natural freshwater lakes. Here, we used pyrosequencing to analyze the bacterioplankton communities in 25 discrete freshwater lakes in Denmark with pH levels ranging from 3.8 to 8.8. We found that pH was the key factor impacting lacustrine bacterioplankton community assembly. More acidic lakes imposed stronger environmental filtering, which decreased the richness and evenness of bacterioplankton operational taxonomic units (OTUs) and largely shifted community composition. Although environmental filtering was determined to be the most important determinant of bacterioplankton community assembly, the importance of neutral assembly processes must also be considered, notably in acidic lakes, where the species (OTU) diversity was low. We observed that the strong effect of environmental filtering in more acidic lakes was weakened by the enhanced relative importance of neutral community assembly, and bacterioplankton communities tended to be less phylogenetically clustered in more acidic lakes. In summary, we propose that pH is a major environmental determinant in freshwater lakes, regulating the relative importance and interplay between niche-related and neutral processes and shaping the patterns of freshwater lake bacterioplankton biodiversity.     

Phylogenetic analysis suggests that habitat filtering is structuring marine bacterial communities across the globe

The phylogenetic structure and community composition were analysed in an existing data set of marine bacterioplankton communities to elucidate the evolutionary and ecological processes dictating the assembly. The communities were sampled from coastal waters at nine locations distributed worldwide and were examined through the use of comprehensive clone libraries of 16S ribosomal RNA genes. The analyses show that the local communities are phylogenetically different from each other and that a majority of them are phylogenetically clustered, i.e. the species (operational taxonomic units) were more related to each other than expected by chance. Accordingly, the local communities were assembled non-randomly from the global pool of available bacterioplankton. Further, the phylogenetic structures of the communities were related to the water temperature at the locations. In agreement with similar studies, including both macroorganisms and bacteria, these results suggest that marine bacterial communities are structured by “habitat filtering”, i.e. through non-random colonization and invasion determined by environmental characteristics. Different bacterial types seem to have different ecological niches that dictate their survival in different habitats. Other eco-evolutionary processes that may contribute to the observed phylogenetic patterns are discussed. The results also imply a mapping between phenotype and phylogenetic relatedness which facilitates the use of community phylogenetic structure analysis to infer ecological and evolutionary assembly processes.     

Richness and diversity of bacterioplankton species along an estuarine gradient in Moreton Bay, Australia

Bacterioplankton community diversity was investigated in the subtropical Brisbane River-Moreton Bay estuary, Australia (27 degrees 25 minutes S, 153 degrees 5 minutes E). Bacterial communities were studied using automated rRNA intergenic spacer analysis (ARISA), which amplifies 16S-23S ribosomal DNA internally transcribed spacer regions from mixed-community DNA and detects the separated products on a fragment analyzer. Samples were collected from eight sites throughout the estuary and east to the East Australian Current (Coral Sea). Bacterioplankton communities had the highest operational taxonomic unit (OTU) richness, as measured by ARISA at eastern bay stations (S [total richness] = 84 to 85 OTU) and the lowest richness in the Coral Sea (S = 39 to 59 OTU). Richness correlated positively with bacterial abundance; however, there were no strong correlations between diversity and salinity, NO(3)(-) and PO(4)(3-) concentrations, or chlorophyll a concentration. Bacterioplankton communities at the riverine stations were different from communities in the bay or Coral Sea. The main differences in OTU richness between stations were in taxa that each represented 0.1% (the detection limit) to 0.5% of the total amplified DNA, i.e., the "tail" of the distribution. We found that some bacterioplankton taxa are specific to distinct environments while others have a ubiquitous distribution from river to sea. Bacterioplankton richness and diversity patterns in the estuary are potentially a consequence of greater niche availability, mixing of local and adjacent environment communities, or intermediate disturbance. Furthermore, these results contrast with previous reports of spatially homogeneous bacterioplankton communities in other coastal waters.     

Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high‐throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and “populations” of various species in our communities, we examine the impact of intra‐ and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59–84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31–63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group‐specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.     

Effects of long‐term differential fertilization on eukaryotic microbial communities in an arable soil: a multiple barcoding approach

To understand the fine‐scale effects of changes in nutrient availability on eukaryotic soil microorganisms communities, a multiple barcoding approach was used to analyse soil samples from four different treatments in a long‐term fertilization experiment. We performed PCR amplification on soil DNA with primer pairs specifically targeting the 18S rRNA genes of all eukaryotes and three protist groups (Cercozoa, Chrysophyceae‐Synurophyceae and Kinetoplastida) as well as the ITS gene of fungi and the 23S plastid rRNA gene of photoautotrophic microorganisms. Amplicons were pyrosequenced, and a total of 88 706 quality filtered reads were clustered into 1232 operational taxonomic units (OTU) across the six data sets. Comparisons of the taxonomic coverage achieved based on overlapping assignment of OTUs revealed that half of the eukaryotic taxa identified were missed by the universal eukaryotic barcoding marker. There were only little differences in OTU richness observed between organic‐ (farmyard manure), mineral‐ and nonfertilized soils. However, the community compositions appeared to be strongly structured by organic fertilization in all data sets other than that generated using the universal eukaryotic 18S rRNA gene primers, whereas mineral fertilization had only a minor effect. In addition, a co‐occurrence based network analysis revealed complex potential interaction patterns between OTUs from different trophic levels, for example between fungivorous flagellates and fungi. Our results demonstrate that changes in pH, moisture and organic nutrients availability caused shifts in the composition of eukaryotic microbial communities at multiple trophic levels.     

Patterns of fungal communities among and within decaying logs, revealed by 454 sequencing

Owing to previous methodological limitations, knowledge about the fine-scale distribution of fungal mycelia in decaying logs is limited. We investigated fungal communities in decaying Norway spruce logs at various spatial scales at two environmentally different locations in Sweden. On the basis of 454 pyrosequencing of the ITS2 region of rDNA, 1914 operational taxonomic units (OTUs) were detected in 353 samples. The communities differed significantly among logs, but the physical distance between logs was not found to have a significant effect on whether fungal communities had any resemblance to each other. Within a log, samples that were closer together generally had communities that showed more resemblance to each other than those that were further apart. OTUs characteristic for particular positions on the logs could be identified. In general, these OTUs did not overlap with the most abundant OTUs, and their ecological role was often unknown. Only a few OTUs were detected in the majority of logs, whereas numerous OTUs were rare and present in only one or a few logs. Wood-decaying Basidiomycetes were often represented by higher sequence reads in individual logs than Ascomycete OTUs, suggesting that Basidiomycete mycelia spread out more rapidly when established. OTU richness tended to increase with the decay stage of the sample; however, the known wood decayers were most abundant in less-decomposed samples. The fungi identified in the logs represented different ecological strategies. Our findings differ from previously published sporocarp studies, indicating that the highly abundant fruiting species may respond to environment in different ways than the rest of the fungal community.     

Locality or habitat? Exploring predictors of biodiversity in Amazonia

Amazonia is an environmentally heterogeneous and biologically megadiverse region, and its biodiversity varies considerably over space. However, existing knowledge on Amazonian biodiversity and its environmental determinants stems almost exclusively from studies of macroscopic above‐ground organisms, notably vertebrates and trees. In contrast, diversity patterns of most other organisms remain elusive, although some of them, for instance microorganisms, constitute the overwhelming majority of taxa in any given location, both in terms of diversity and abundance. Here, we use DNA metabarcoding to estimate prokaryote and eukaryote diversity in environmental soil and litter samples from 39 survey plots in a longitudinal transect across Brazilian Amazonia using 16S and 18S gene sequences, respectively. We characterize richness and community composition based on operational taxonomic units (OTUs) and test their correlation with longitude and habitat. We find that prokaryote and eukaryote OTU richness and community composition differ significantly among localities and habitats, and that prokaryotes are more strongly structured by locality and habitat type than eukaryotes. Our results 1) provide a first large‐scale mapping of Amazonian soil biodiversity, suggesting that OTU richness patterns might follow substantially different patterns from those observed for macro‐organisms; and 2) indicate that locality and habitat factors interact in determining OTU richness patterns and community composition. This study shows the potential of DNA metabarcoding in unveiling Amazonia's outstanding diversity, despite the lack of complete reference sequence databases for the organisms sequenced.     

Richness and Diversity of Bacterioplankton Species along an Estuarine Gradient in Moreton Bay, Australia

Bacterioplankton community diversity was investigated in the subtropical Brisbane River-Moreton Bay estuary, Australia (27°25′S, 153°5′E). Bacterial communities were studied using automated rRNA intergenic spacer analysis (ARISA), which amplifies 16S-23S ribosomal DNA internally transcribed spacer regions from mixed-community DNA and detects the separated products on a fragment analyzer. Samples were collected from eight sites throughout the estuary and east to the East Australian Current (Coral Sea). Bacterioplankton communities had the highest operational taxonomic unit (OTU) richness, as measured by ARISA at eastern bay stations (S [total richness] = 84 to 85 OTU) and the lowest richness in the Coral Sea (S = 39 to 59 OTU). Richness correlated positively with bacterial abundance; however, there were no strong correlations between diversity and salinity, NO₃⁻ and PO₄³⁻ concentrations, or chlorophyll a concentration. Bacterioplankton communities at the riverine stations were different from communities in the bay or Coral Sea. The main differences in OTU richness between stations were in taxa that each represented 0.1% (the detection limit) to 0.5% of the total amplified DNA, i.e., the “tail” of the distribution. We found that some bacterioplankton taxa are specific to distinct environments while others have a ubiquitous distribution from river to sea. Bacterioplankton richness and diversity patterns in the estuary are potentially a consequence of greater niche availability, mixing of local and adjacent environment communities, or intermediate disturbance. Furthermore, these results contrast with previous reports of spatially homogeneous bacterioplankton communities in other coastal waters.     

Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of bacteria, archaea, fungi, and viruses in soil

Recent studies have highlighted the surprising richness of soil bacterial communities; however, bacteria are not the only microorganisms found in soil. To our knowledge, no study has compared the diversities of the four major microbial taxa, i.e., bacteria, archaea, fungi, and viruses, from an individual soil sample. We used metagenomic and small-subunit RNA-based sequence analysis techniques to compare the estimated richness and evenness of these groups in prairie, desert, and rainforest soils. By grouping sequences at the 97% sequence similarity level (an operational taxonomic unit [OTU]), we found that the archaeal and fungal communities were consistently less even than the bacterial communities. Although total richness levels are difficult to estimate with a high degree of certainty, the estimated number of unique archaeal or fungal OTUs appears to rival or exceed the number of unique bacterial OTUs in each of the collected soils. In this first study to comprehensively survey viral communities using a metagenomic approach, we found that soil viruses are taxonomically diverse and distinct from the communities of viruses found in other environments that have been surveyed using a similar approach. Within each of the four microbial groups, we observed minimal taxonomic overlap between sites, suggesting that soil archaea, bacteria, fungi, and viruses are globally as well as locally diverse.     

Metagenomic and Small-Subunit rRNA Analyses Reveal the Genetic Diversity of Bacteria, Archaea, Fungi, and Viruses in Soil

Recent studies have highlighted the surprising richness of soil bacterial communities; however, bacteria are not the only microorganisms found in soil. To our knowledge, no study has compared the diversities of the four major microbial taxa, i.e., bacteria, archaea, fungi, and viruses, from an individual soil sample. We used metagenomic and small-subunit RNA-based sequence analysis techniques to compare the estimated richness and evenness of these groups in prairie, desert, and rainforest soils. By grouping sequences at the 97% sequence similarity level (an operational taxonomic unit [OTU]), we found that the archaeal and fungal communities were consistently less even than the bacterial communities. Although total richness levels are difficult to estimate with a high degree of certainty, the estimated number of unique archaeal or fungal OTUs appears to rival or exceed the number of unique bacterial OTUs in each of the collected soils. In this first study to comprehensively survey viral communities using a metagenomic approach, we found that soil viruses are taxonomically diverse and distinct from the communities of viruses found in other environments that have been surveyed using a similar approach. Within each of the four microbial groups, we observed minimal taxonomic overlap between sites, suggesting that soil archaea, bacteria, fungi, and viruses are globally as well as locally diverse.     

Diversity of bacterioplankton in coastal seawaters of Fildes Peninsula,King George Island,Antarctica

The bacterioplankton not only serves critical functions in marine nutrient cycles, but can also serve as indicators of the marine environment. The compositions of bacterial communities in the surface seawater of Ardley Cove and Great Wall Cove were analyzed using a 16S rRNA multiplex 454 pyrosequencing approach. Similar patterns of bacterial composition were found between the two coves, in which Bacteroidetes, Alphaproteobacteria, and Gammaproteobacteria were the dominant members of the bacterioplankton communities. In addition, a large fraction of the bacterial sequence reads (on average 5.3 % per station) could not be assigned below the domain level. Compared with Ardley Cove, Great Wall Cove showed higher chlorophyll and particulate organic carbon concentrations and exhibited relatively lower bacterial richness and diversity. Inferred metabolisms of summer bacterioplankton in the two coves were characterized by chemoheterotrophy and photoheterotrophy. Results suggest that some cosmopolitan species (e.g., Polaribacter and Sulfitobacter) belonging to a few bacterial groups that usually dominate in marine bacterioplankton communities may have similar ecological functions in similar marine environments but at different geographic locations.     

Optimization of Terminal-Restriction Fragment Length Polymorphism Analysis for Complex Marine Bacterioplankton Communities and Comparison with Denaturing Gradient Gel Electrophoresis

The potential of terminal-restriction fragment length polymorphism (T-RFLP) and the detection of operational taxonomic units (OTUs) by capillary electrophoresis (CE) to characterize marine bacterioplankton communities was compared with that of denaturing gradient gel electrophoresis (DGGE). A protocol has been developed to optimize the separation and detection of OTUs between 20 and 1,632 bp by using CE and laser-induced fluorescence detection. Additionally, we compared T-RFLP fingerprinting to DGGE optimized for detection of less abundant OTUs. Similar results were obtained with both fingerprinting techniques, although the T-RFLP approach and CE detection of OTUs was more sensitive, as indicated by the higher number of OTUs detected. We tested the T-RFLP fingerprinting technique on complex marine bacterial communities by using the 16S rRNA gene and 16S rRNA as templates for PCR. Samples from the Northern and Middle Adriatic Sea and from the South and North Aegean Sea were compared. Distinct clusters were identifiable for different sampling sites. Thus, this technique is useful for rapid evaluation of the biogeographical distribution and relationships of bacterioplankton communities.     

Optimization of terminal-restriction fragment length polymorphism analysis for complex marine bacterioplankton communities and comparison with denaturing gradient gel electrophoresis.

The potential of terminal-restriction fragment length polymorphism (T-RFLP) and the detection of operational taxonomic units (OTUs) by capillary electrophoresis (CE) to characterize marine bacterioplankton communities was compared with that of denaturing gradient gel electrophoresis (DGGE). A protocol has been developed to optimize the separation and detection of OTUs between 20 and 1, 632 bp by using CE and laser-induced fluorescence detection. Additionally, we compared T-RFLP fingerprinting to DGGE optimized for detection of less abundant OTUs. Similar results were obtained with both fingerprinting techniques, although the T-RFLP approach and CE detection of OTUs was more sensitive, as indicated by the higher number of OTUs detected. We tested the T-RFLP fingerprinting technique on complex marine bacterial communities by using the 16S rRNA gene and 16S rRNA as templates for PCR. Samples from the Northern and Middle Adriatic Sea and from the South and North Aegean Sea were compared. Distinct clusters were identifiable for different sampling sites. Thus, this technique is useful for rapid evaluation of the biogeographical distribution and relationships of bacterioplankton communities.     

Water mass-specificity of bacterial communities in the North Atlantic revealed by massively parallel sequencing

Bacterial assemblages from subsurface (100 m depth), meso- (200-1000 m depth) and bathy-pelagic (below 1000 m depth) zones at 10 stations along a North Atlantic Ocean transect from 60°N to 5°S were characterized using massively parallel pyrotag sequencing of the V6 region of the 16S rRNA gene (V6 pyrotags). In a dataset of more than 830,000 pyrotags, we identified 10,780 OTUs of which 52% were singletons. The singletons accounted for less than 2% of the OTU abundance, whereas the 100 and 1000 most abundant OTUs represented 80% and 96% respectively of all recovered OTUs. Non-metric Multi-Dimensional Scaling and Canonical Correspondence Analysis of all the OTUs excluding the singletons revealed a clear clustering of the bacterial communities according to the water masses. More than 80% of the 1000 most abundant OTUs corresponded to Proteobacteria of which 55% were Alphaproteobacteria, mostly composed of the SAR11 cluster. Gammaproteobacteria increased with depth and included a relatively large number of OTUs belonging to Alteromonadales and Oceanospirillales. The bathypelagic zone showed higher taxonomic evenness than the overlying waters, albeit bacterial diversity was remarkably variable. Both abundant and low-abundance OTUs were responsible for the distinct bacterial communities characterizing the major deep-water masses. Taken together, our results reveal that deep-water masses act as bio-oceanographic islands for bacterioplankton leading to water mass-specific bacterial communities in the deep waters of the Atlantic.     

Synthetic Statistical Approach Reveals a High Degree of Richness of Microbial Eukaryotes in an Anoxic Water Column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.     

Patterns of community composition and diversity in latent fungi of living Quercus serrata trunks across a range of oak wilt prevalence and climate variables in Japan

Given that forest dieback due to emerging pests is increasing under global warming, understanding the relationships between pests, climate, and forest biodiversity is an urgent priority. In Japan, mass attacks of an ambrosia beetle, vectoring a pathogenic fungus, cause oak wilt outbreaks in recent decades. Here, the associations of oak wilt and climate with wood-inhabiting fungal communities in apparently healthy Quercus serrata trunks were investigated using DNA metabarcoding in seven sites along a climatic gradient in Japan. Amplicon sequencing of the fungal internal transcribed spacer 1 region generated 1,339,958 sequence reads containing 879 fungal operational taxonomic units (OTUs) in 234 wood samples. OTU compositions were significantly different between sites with and without oak wilt. OTU richness increased with temperature and precipitation at sites where oak wilt was present, but this relationship was not observed at sites without oak wilt, possibly due to interaction between oak wilt and climate.     

Characterization of Lysogens in Bacterioplankton Assemblages of the Southern California Borderland

Viruses cause significant mortality of marine microorganisms; however, their role in shaping the composition of microbial assemblages has not been fully elucidated. Because viruses may form lysogenic relationships with their hosts, temperate viruses may influence bacterial assemblage structures through direct lysis of hosts when induced by environmental stimuli or by homoimmunity (i.e., immunity to closely related viruses). We investigated the components of bacterioplankton assemblages that bore prophage using the lysogenic induction agent mitomycin C. Seawater was collected at two locations (the San Pedro Ocean Time Series Station and in the Santa Barbara Channel) in the Southern California Borderland and amended with mitomycin C. After 24-h incubation, the community structure of bacterioplankton was compared with unamended controls using automated rRNA intergenic spacer analysis. The addition of mitomycin C to seawater had effects on the community structure of bacterioplankton, stimulating detectable overall diversity and richness of fingerprints and causing the assemblages within incubations to become different to control assemblages. Most negatively impacted operational taxonomic units (OTU) in mitomycin C-amended incubations individually comprised a large fraction of total amplified DNA in initial seawater (5.3-23.3% of amplified DNA fluorescence) fingerprints, and data suggest that these include organisms putatively classified as members of the gamma-Proteobacteria, SAR11 cluster, and Synechococcus groups. The stimulation of assemblage richness by induction of lysogens, and the reduction in the contribution to total DNA of common OTU (and concomitant increase in rare OTU), suggests that temperate phage have the potential to strongly influence the diversity of bacterioplankton assemblages. Because lysogenic OTU may also be resistant to closely related lytic (i.e., free-living) viruses, the impact of lytic virioplankton on assemblages may only be pronounced transiently or when conditions causing lysogenic induction arise.     

Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness

Although copious qualitative information describes the members of the diverse microbial communities on Earth, statistical approaches for quantifying and comparing the numbers and compositions of lineages in communities are lacking. We present a method that addresses the challenge of assigning sequences to operational taxonomic units (OTUs) based on the genetic distances between sequences. We developed a computer program, DOTUR, which assigns sequences to OTUs by using either the furthest, average, or nearest neighbor algorithm for each distance level. DOTUR uses the frequency at which each OTU is observed to construct rarefaction and collector's curves for various measures of richness and diversity. We analyzed 16S rRNA gene libraries derived from Scottish and Amazonian soils and the Sargasso Sea with DOTUR, which assigned sequences to OTUs rapidly and reliably based on the genetic distances between sequences and identified previous inconsistencies and errors in assigning sequences to OTUs. An analysis of the two 16S rRNA gene libraries from soil demonstrated that they do not contain enough sequences to support a claim that they contain different numbers of bacterial lineages with statistical confidence (P > 0.05), nor do they contain enough sequences to provide a robust estimate of species richness when an OTU is defined as containing sequences that are no more than 3% different from each other. In contrast, the richness of OTUs at the 3% level in the Sargasso Sea collection began to plateau after the sampling of 690 sequences. We anticipate that an equivalent extent of sampling for soil would require sampling more than 10,000 sequences, almost 100 times the size of typical sequence collections obtained from soil.     

Introducing DOTUR, a Computer Program for Defining Operational Taxonomic Units and Estimating Species Richness

Although copious qualitative information describes the members of the diverse microbial communities on Earth, statistical approaches for quantifying and comparing the numbers and compositions of lineages in communities are lacking. We present a method that addresses the challenge of assigning sequences to operational taxonomic units (OTUs) based on the genetic distances between sequences. We developed a computer program, DOTUR, which assigns sequences to OTUs by using either the furthest, average, or nearest neighbor algorithm for each distance level. DOTUR uses the frequency at which each OTU is observed to construct rarefaction and collector's curves for various measures of richness and diversity. We analyzed 16S rRNA gene libraries derived from Scottish and Amazonian soils and the Sargasso Sea with DOTUR, which assigned sequences to OTUs rapidly and reliably based on the genetic distances between sequences and identified previous inconsistencies and errors in assigning sequences to OTUs. An analysis of the two 16S rRNA gene libraries from soil demonstrated that they do not contain enough sequences to support a claim that they contain different numbers of bacterial lineages with statistical confidence (P > 0.05), nor do they contain enough sequences to provide a robust estimate of species richness when an OTU is defined as containing sequences that are no more than 3% different from each other. In contrast, the richness of OTUs at the 3% level in the Sargasso Sea collection began to plateau after the sampling of 690 sequences. We anticipate that an equivalent extent of sampling for soil would require sampling more than 10,000 sequences, almost 100 times the size of typical sequence collections obtained from soil.     

Synthetic statistical approach reveals a high degree of richness of microbial eukaryotes in an anoxic water column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.