PCR biases distort bacterial and archaeal community structure in pyrosequencing datasets

As 16S rRNA gene targeted massively parallel sequencing has become a common tool for microbial diversity investigations, numerous advances have been made to minimize the influence of sequencing and chimeric PCR artifacts through rigorous quality control measures. However, there has been little effort towards understanding the effect of multi-template PCR biases on microbial community structure. In this study, we used three bacterial and three archaeal mock communities consisting of, respectively, 33 bacterial and 24 archaeal 16S rRNA gene sequences combined in different proportions to compare the influences of (1) sequencing depth, (2) sequencing artifacts (sequencing errors and chimeric PCR artifacts), and (3) biases in multi-template PCR, towards the interpretation of community structure in pyrosequencing datasets. We also assessed the influence of each of these three variables on α- and β-diversity metrics that rely on the number of OTUs alone (richness) and those that include both membership and the relative abundance of detected OTUs (diversity). As part of this study, we redesigned bacterial and archaeal primer sets that target the V3-V5 region of the 16S rRNA gene, along with multiplexing barcodes, to permit simultaneous sequencing of PCR products from the two domains. We conclude that the benefits of deeper sequencing efforts extend beyond greater OTU detection and result in higher precision in β-diversity analyses by reducing the variability between replicate libraries, despite the presence of more sequencing artifacts. Additionally, spurious OTUs resulting from sequencing errors have a significant impact on richness or shared-richness based α- and β-diversity metrics, whereas metrics that utilize community structure (including both richness and relative abundance of OTUs) are minimally affected by spurious OTUs. However, the greatest obstacle towards accurately evaluating community structure are the errors in estimated mean relative abundance of each detected OTU due to biases associated with multi-template PCR reactions.     

A diverse bacterial community in an anoxic quinoline-degrading bioreactor determined by using pyrosequencing and clone library analysis

There is a concern of whether the structure and diversity of a microbial community can be effectively revealed by short-length pyrosequencing reads. In this study, we performed a microbial community analysis on a sample from a high-efficiency denitrifying quinoline-degrading bioreactor and compared the results generated by pyrosequencing with those generated by clone library technology. By both technologies, 16S rRNA gene analysis indicated that the bacteria in the sample were closely related to, for example, Proteobacteria, Actinobacteria, and Bacteroidetes. The sequences belonging to Rhodococcus were the most predominant, and Pseudomonas, Sphingomonas, Acidovorax, and Zoogloea were also abundant. Both methods revealed a similar overall bacterial community structure. However, the 622 pyrosequencing reads of the hypervariable V3 region of the 16S rRNA gene revealed much higher bacterial diversity than the 130 sequences from the full-length 16S rRNA gene clone library. The 92 operational taxonomic unit (OTUs) detected using pyrosequencing belonged to 45 families, whereas the 37 OTUs found in the clone library belonged to 25 families. Most sequences obtained from the clone library had equivalents in the pyrosequencing reads. However, 64 OTUs detected by pyrosequencing were not represented in the clone library. Our results demonstrate that pyrosequencing of the V3 region of the 16S rRNA gene is not only a powerful tool for discovering low-abundance bacterial populations but is also reliable for dissecting the bacterial community structure in a wastewater environment.     

Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure

PCR-based surveys of microbial communities commonly use regions of the small-subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)- and taxonomy-based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.     

Amplification by PCR artificially reduces the proportion of the rare biosphere in microbial communities

The microbial world has been shown to hold an unimaginable diversity. The use of rRNA genes and PCR amplification to assess microbial community structure and diversity present biases that need to be analyzed in order to understand the risks involved in those estimates. Herein, we show that PCR amplification of specific sequence targets within a community depends on the fractions that those sequences represent to the total DNA template. Using quantitative, real-time, multiplex PCR and specific Taqman probes, the amplification of 16S rRNA genes from four bacterial species within a laboratory community were monitored. Results indicate that the relative amplification efficiency for each bacterial species is a nonlinear function of the fraction that each of those taxa represent within a community or multispecies DNA template. Consequently, the low-proportion taxa in a community are under-represented during PCR-based surveys and a large number of sequences might need to be processed to detect some of the bacterial taxa within the 'rare biosphere'. The structure of microbial communities from PCR-based surveys is clearly biased against low abundant taxa which are required to decipher the complete extent of microbial diversity in nature.     

Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies

The advent of next generation sequencing has coincided with a growth in interest in using these approaches to better understand the role of the structure and function of the microbial communities in human, animal, and environmental health. Yet, use of next generation sequencing to perform 16S rRNA gene sequence surveys has resulted in considerable controversy surrounding the effects of sequencing errors on downstream analyses. We analyzed 2.7×10(6) reads distributed among 90 identical mock community samples, which were collections of genomic DNA from 21 different species with known 16S rRNA gene sequences; we observed an average error rate of 0.0060. To improve this error rate, we evaluated numerous methods of identifying bad sequence reads, identifying regions within reads of poor quality, and correcting base calls and were able to reduce the overall error rate to 0.0002. Implementation of the PyroNoise algorithm provided the best combination of error rate, sequence length, and number of sequences. Perhaps more problematic than sequencing errors was the presence of chimeras generated during PCR. Because we knew the true sequences within the mock community and the chimeras they could form, we identified 8% of the raw sequence reads as chimeric. After quality filtering the raw sequences and using the Uchime chimera detection program, the overall chimera rate decreased to 1%. The chimeras that could not be detected were largely responsible for the identification of spurious operational taxonomic units (OTUs) and genus-level phylotypes. The number of spurious OTUs and phylotypes increased with sequencing effort indicating that comparison of communities should be made using an equal number of sequences. Finally, we applied our improved quality-filtering pipeline to several benchmarking studies and observed that even with our stringent data curation pipeline, biases in the data generation pipeline and batch effects were observed that could potentially confound the interpretation of microbial community data.     

Low sequencing efforts bias analyses of shared taxa in microbial communities

The potential for comparing microbial community population structures has been greatly enhanced by developments in next generation sequencing methods that can generate hundreds of thousands to millions of reads in a single run. Conversely, many microbial community comparisons have been published with no more than 1,000 sequences per sample. These studies have presented data on levels of shared operational taxonomic units (OTUs) between communities. Due to lack of coverage, that approach might compromise the conclusions about microbial diversity and the degree of difference between environments. In this study, we present data from recent studies that highlight this problem. Also, we analyzed datasets of 16 rRNA sequences with small and high sequence coverage from different environments to demonstrate that the level of sequencing effort used for analyzing microbial communities biases the results. We randomly sampled pyrosequencing-generated 16S rRNA gene libraries with increasing sequence effort. Sequences were used to calculate Good's coverage, the percentage of shared OTUs, and phylogenetic distance measures. Our data showed that simple counts of presence/absence of taxonomic unities do not reflect the real similarity in membership and structure of the bacterial communities and that community comparisons based on phylogenetic tests provide a way to test statistically significant differences between two or more environments without need an exhaustive sampling effort.     

Simultaneous Amplicon Sequencing to Explore Co-Occurrence Patterns of Bacterial,Archaeal and Eukaryotic Microorganisms in Rumen Microbial Communities

Ruminants rely on a complex rumen microbial community to convert dietary plant material to energy-yielding products. Here we developed a method to simultaneously analyze the community''s bacterial and archaeal 16S rRNA genes, ciliate 18S rRNA genes and anaerobic fungal internal transcribed spacer 1 genes using 12 DNA samples derived from 11 different rumen samples from three host species (Ovis aries, Bos taurus, Cervus elephas) and multiplex 454 Titanium pyrosequencing. We show that the mixing ratio of the group-specific DNA templates before emulsion PCR is crucial to compensate for differences in amplicon length. This method, in contrast to using a non-specific universal primer pair, avoids sequencing non-targeted DNA, such as plant- or endophyte-derived rRNA genes, and allows increased or decreased levels of community structure resolution for each microbial group as needed. Communities analyzed with different primers always grouped by sample origin rather than by the primers used. However, primer choice had a greater impact on apparent archaeal community structure than on bacterial community structure, and biases for certain methanogen groups were detected. Co-occurrence analysis of microbial taxa from all three domains of life suggested strong within- and between-domain correlations between different groups of microorganisms within the rumen. The approach used to simultaneously characterize bacterial, archaeal and eukaryotic components of a microbiota should be applicable to other communities occupying diverse habitats.     

Developing SSU rDNA metagenomic profiles of aquatic microbial communities for environmental assessments

Five water samples from three sources, two municipal reservoirs in central North Carolina and Toolik Lake in Alaska, were processed to conduct a comparative survey of microbial small subunit rDNA sequences. Genomic DNA was extracted and amplified by PCR using universal SSU rDNA primers to generate 16S and 18S rDNA clone libraries and 50 clones from each library were sequenced and placed in operational taxonomic units (OTUs). Through this recovery and analysis of SSU rRNA genes, a metagenomic profile of the microbial community emerged for each environmental sample. Analyses of these profiles, including species diversity estimates and rank-abundance curves, revealed that approximately 64% of prokaryotic OTUs and 80% of eukaryotic OTUs were novel. Diversity estimates were consistent with predicted ecosystem characteristics: they were greater for the mesotrophic to eutrophic temperate lakes, than for the oligotrophic arctic lake. Sample comparisons showed that community similarity declined as geographic distance between sites increased. Real-time quantitative PCR results showed that OTUs which had been recovered from only one library were actually present in other samples, but at much lower frequencies, suggesting that many, if not most, microorganisms are cosmopolitan. Together, these results support the potential value of using the microbial community as an indicator of local environmental conditions. In other words, it may be realistic to monitor water quality using a single, comprehensive suite of microorganisms by analyzing patterns of relative abundance.     

Exploring Microbial Diversity and Taxonomy Using SSU rRNA Hypervariable Tag Sequencing

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA. As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe. However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process. The resulting taxonomic census provides information on both composition and diversity of the microbial community. To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent. The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match. The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the “rare biosphere” than does capillary sequencing of the full-length gene. In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation. This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.     

Archaea Appear to Dominate the Microbiome of Inflatella pellicula Deep Sea Sponges

Microbes associated with marine sponges play significant roles in host physiology. Remarkable levels of microbial diversity have been observed in sponges worldwide through both culture-dependent and culture-independent studies. Most studies have focused on the structure of the bacterial communities in sponges and have involved sponges sampled from shallow waters. Here, we used pyrosequencing of 16S rRNA genes to compare the bacterial and archaeal communities associated with two individuals of the marine sponge Inflatella pellicula from the deep-sea, sampled from a depth of 2,900 m, a depth which far exceeds any previous sequence-based report of sponge-associated microbial communities. Sponge-microbial communities were also compared to the microbial community in the surrounding seawater. Sponge-associated microbial communities were dominated by archaeal sequencing reads with a single archaeal OTU, comprising ∼60% and ∼72% of sequences, being observed from Inflatella pellicula. Archaeal sequencing reads were less abundant in seawater (∼11% of sequences). Sponge-associated microbial communities were less diverse and less even than any other sponge-microbial community investigated to date with just 210 and 273 OTUs (97% sequence identity) identified in sponges, with 4 and 6 dominant OTUs comprising ∼88% and ∼89% of sequences, respectively. Members of the candidate phyla, SAR406, NC10 and ZB3 are reported here from sponges for the first time, increasing the number of bacterial phyla or candidate divisions associated with sponges to 43. A minor cohort from both sponge samples (∼0.2% and ∼0.3% of sequences) were not classified to phylum level. A single OTU, common to both sponge individuals, dominates these unclassified reads and shares sequence homology with a sponge associated clone which itself has no known close relative and may represent a novel taxon.     

Hemolymph microbiome of Pacific oysters in response to temperature,temperature stress and infection

Microbiota provide their hosts with a range of beneficial services, including defense from external pathogens. However, host-associated microbial communities themselves can act as a source of opportunistic pathogens depending on the environment. Marine poikilotherms and their microbiota are strongly influenced by temperature, but experimental studies exploring how temperature affects the interactions between both parties are rare. To assess the effects of temperature, temperature stress and infection on diversity, composition and dynamics of the hemolymph microbiota of Pacific oysters (Crassostrea gigas), we conducted an experiment in a fully-crossed, three-factorial design, in which the temperature acclimated oysters (8 or 22 °C) were exposed to temperature stress and to experimental challenge with a virulent Vibrio sp. strain. We monitored oyster survival and repeatedly collected hemolymph of dead and alive animals to determine the microbiome composition by 16s rRNA gene amplicon pyrosequencing. We found that the microbial dynamics and composition of communities in healthy animals (including infection survivors) were significantly affected by temperature and temperature stress, but not by infection. The response was mediated by changes in the incidence and abundance of operational taxonomic units (OTUs) and accompanied by little change at higher taxonomic levels, indicating dynamic stability of the hemolymph microbiome. Dead and moribund oysters, on the contrary, displayed signs of community structure disruption, characterized by very low diversity and proliferation of few OTUs. We can therefore link short-term responses of host-associated microbial communities to abiotic and biotic factors and assess the potential feedback between microbiota dynamics and host survival during disease.     

Reproducibility and quantitation of amplicon sequencing-based detection

To determine the reproducibility and quantitation of the amplicon sequencing-based detection approach for analyzing microbial community structure, a total of 24 microbial communities from a long-term global change experimental site were examined. Genomic DNA obtained from each community was used to amplify 16S rRNA genes with two or three barcode tags as technical replicates in the presence of a small quantity (0.1% wt/wt) of genomic DNA from Shewanella oneidensis MR-1 as the control. The technical reproducibility of the amplicon sequencing-based detection approach is quite low, with an average operational taxonomic unit (OTU) overlap of 17.2%±2.3% between two technical replicates, and 8.2%±2.3% among three technical replicates, which is most likely due to problems associated with random sampling processes. Such variations in technical replicates could have substantial effects on estimating β-diversity but less on α-diversity. A high variation was also observed in the control across different samples (for example, 66.7-fold for the forward primer), suggesting that the amplicon sequencing-based detection approach could not be quantitative. In addition, various strategies were examined to improve the comparability of amplicon sequencing data, such as increasing biological replicates, and removing singleton sequences and less-representative OTUs across biological replicates. Finally, as expected, various statistical analyses with preprocessed experimental data revealed clear differences in the composition and structure of microbial communities between warming and non-warming, or between clipping and non-clipping. Taken together, these results suggest that amplicon sequencing-based detection is useful in analyzing microbial community structure even though it is not reproducible and quantitative. However, great caution should be taken in experimental design and data interpretation when the amplicon sequencing-based detection approach is used for quantitative analysis of the β-diversity of microbial communities.     

Evidence from GC-TRFLP that bacterial communities in soil are lognormally distributed

The Species Abundance Distribution (SAD) is a fundamental property of ecological communities and the form and formation of SADs have been examined for a wide range of communities including those of microorganisms. Progress in understanding microbial SADs, however, has been limited by the remarkable diversity and vast size of microbial communities. As a result, few microbial systems have been sampled with sufficient depth to generate reliable estimates of the community SAD. We have used a novel approach to characterize the SAD of bacterial communities by coupling genomic DNA fractionation with analysis of terminal restriction fragment length polymorphisms (GC-TRFLP). Examination of a soil microbial community through GC-TRFLP revealed 731 bacterial operational taxonomic units (OTUs) that followed a lognormal distribution. To recover the same 731 OTUs through analysis of DNA sequence data is estimated to require analysis of 86,264 16S rRNA sequences. The approach is examined and validated through construction and analysis of simulated microbial communities in silico. Additional simulations performed to assess the potential effects of PCR bias show that biased amplification can cause a community whose distribution follows a power-law function to appear lognormally distributed. We also show that TRFLP analysis, in contrast to GC-TRFLP, is not able to effectively distinguish between competing SAD models. Our analysis supports use of the lognormal as the null distribution for studying the SAD of bacterial communities as for plant and animal communities.     

Pyrosequencing Reveals Diverse and Distinct Sponge-Specific Microbial Communities in Sponges from a Single Geographical Location in Irish Waters

Marine sponges are host to numerically vast and phylogenetically diverse bacterial communities, with 26 major phyla to date having been found in close association with sponge species worldwide. Analyses of these microbial communities have revealed many sponge-specific novel genera and species. These endosymbiotic microbes are believed to play significant roles in sponge physiology including the production of an array of bioactive secondary metabolites. Here, we report on the use of culture-based and culture-independent (pyrosequencing) techniques to elucidate the bacterial community profiles associated with the marine sponges Raspailia ramosa and Stelligera stuposa sampled from a single geographical location in Irish waters and with ambient seawater. To date, little is known about the microbial ecology of sponges of these genera. Culture isolation grossly underestimated sponge-associated bacterial diversity. Four bacterial phyla (Actinobacteria, Bacteroidetes, Firmicutes, Proteobacteria) were represented amongst ~200 isolates, compared with ten phyla found using pyrosequencing. Long average read lengths of ~430 bp (V1-V3 region of 16S rRNA gene) allowed for robust resolution of sequences to genus level. Bacterial OTUs (2,109 total), at 95% sequence similarity, from ten bacterial phyla were recovered from R. ramosa, 349 OTUs were identified in S. stuposa representing eight phyla, while 533 OTUs from six phyla were found in surrounding seawater. Bacterial communities differed significantly between sponge species and the seawater. Analysis of the data for sponge-specific taxa revealed that 2.8% of classified reads from the sponge R. ramosa can be defined as sponge-specific, while 26% of S. stuposa sequences represent sponge-specific bacteria. Novel sponge-specific clusters were identified, whereas the majority of previously reported sponge-specific clusters (e.g. Poribacteria) were absent from these sponge species. This deep and robust analysis provides further evidence that the microbial communities associated with marine sponge species are highly diverse and divergent from one another and appear to be host-selected through as yet unknown processes.     

Microbial community in a pilot-scale bioreactor promoting anaerobic digestion and sulfur-driven denitrification for domestic sewage treatment

A pilot-scale reactor treating domestic sewage was operated to promote anaerobic digestion and denitrification using endogenous electron donors. While 55 % of organic matter was removed, nitrogen and sulfur showed a different dynamics during the operation. Pyrosequencing analysis clarified this behavior revealing that specific microbial communities inhabited the anaerobic (47.05 % of OTUs) and anoxic (31.39 % of OTUs) chambers. Analysis of 16S rRNA gene partial sequences obtained through pyrosequencing revealed a total of 1727 OTUs clustered at a 3 % distance cutoff. In the anaerobic chamber, microbial community was comprised of fermentative, syntrophic and sulfate-reducing bacteria. The majority of sequences were related to Aminobacterium and Syntrophorhabdus. In the anoxic chamber, the majority of sequences were related to mixotrophic and strictly autotrophic denitrifiers Arcobacter and Sulfuricurvum, respectively, both involved in sulfur-driven denitrification. These results show that pyrosequencing was a powerful tool to investigate the microbial panorama of a complex system, providing new insights to the improvement of the system.     

16S rRNA基因在微生物生态学中的应用

16S rRNA(Small subunit ribosomal RNA)基因是对原核微生物进行系统进化分类研究时最常用的分子标志物(Biomarker),广泛应用于微生物生态学研究中。近些年来随着高通量测序技术及数据分析方法等的不断进步,大量基于16S rRNA基因的研究使得微生物生态学得到了快速发展,然而使用16S rRNA基因作为分子标志物时也存在诸多问题,比如水平基因转移、多拷贝的异质性、基因扩增效率的差异、数据分析方法的选择等,这些问题影响了微生物群落组成和多样性分析时的准确性。对当前使用16S rRNA基因分析微生物群落组成和多样性的进展情况做一总结,重点讨论当前存在的主要问题以及各种分析方法的发展,尤其是与高通量测序技术有关的实验和数据处理问题。     

Small core communities and high variability in bacteria associated with the introduced ascidian Styela plicata

The solitary ascidian Styela plicata is an introduced species in harbors of temperate and tropical oceans around the world. The invasive potential of this species has been studied through reproductive biology and population genetics but no study has yet examined the microbial diversity associated with this ascidian and its potential role in host ecology and invasiveness. Here, we used 16S rRNA gene tag pyrosequencing and transmission electron microscopy to characterize the abundance, diversity and host-specificity of bacteria associated with 3 Mediterranean individuals of S. plicata. Microscopy revealed low bacterial abundance in the inner tunic and their absence from gonad tissues, while pyrosequencing revealed a high diversity of S. plicata-associated bacteria (284 OTUs from 16 microbial phyla) in the inner tunic. The core symbiont community was small and consisted of 16 OTUs present in all S. plicata hosts. This core community included a recently described ascidian symbiont (Hasllibacter halocynthiae) and several known sponge and coral symbionts, including a strictly anaerobic Chloroflexi lineage. Most recovered bacterial OTUs (79.6 %) were present in single S. plicata individuals and statistical analyses of genetic diversity and community structure confirmed high variability of bacterial communities among host individuals. These results suggest that diverse and variable bacterial communities inhabit the tunic of S. plicata, including environmental and host-associated bacterial lineages that appear to be re-established each host generation. We hypothesize that bacterial communities in S. plicata are dynamic and have the potential to aid host acclimation to new habitats by establishing relationships with beneficial, locally sourced bacteria.     

Changes in the Soil Bacterial Communities in a Cedar Plantation Invaded by Moso Bamboo

Moso bamboo is fast-growing and negatively allelopathic to neighboring plants. However, there is little information on the effects of its establishment and expansion to adjacent forest soil communities. To better understand the impacts of bamboo invasion on soil communities, the phylogenetic structure and diversity of the soil bacterial communities in moso bamboo forest, adjacent Japanese cedar plantation, and bamboo-invaded transition zone were examined using a combination of 16S rRNA gene clone libraries and bar-coded pyrosequencing techniques. Based on the number of operational taxonomic units (OTUs), Shannon diversity index, Chao1 estimator, and rarefaction analysis of both techniques, the bamboo soil bacterial community was the most diverse, followed by the transition zone, with the cedar plantation possessing the lowest diversity. The results from both techniques revealed that the Acidobacteria and Proteobacteria predominated in the three communities, though the relative abundance was different. The 250 most abundant OTUs represented about 70 % of the total sequences found by pyrosequencing. Most of these OTUs were found in all three soil communities, demonstrating the overall similarity among the bacterial communities. Nonmetric multidimensional scaling analysis showed further that the bamboo and transition soil communities were more similar with each other than the cedar soils. These results suggest that bamboo invasion to the adjacent cedar plantation gradually increased the bacterial diversity and changed the soil community. In addition, while the 10 most abundant OTUs were distributed worldwide, related sequences were not abundant in soils from outside the forest studied here. This result may be an indication of the uniqueness of this region.     

Impact of sequence processing and taxonomic classification approaches on eukaryotic community structure from environmental samples with emphasis on diatoms

Next‐generation sequencing is a common method for analysing microbial community diversity and composition. Configuring an appropriate sequence processing strategy within the variety of tools and methods is a nontrivial task and can considerably influence the resulting community characteristics. We analysed the V4 region of 18S rRNA gene sequences of marine samples by 454‐pyrosequencing. Along this process, we generated several data sets with QIIME, mothur, and a custom‐made pipeline based on DNAStar and the phylogenetic tree‐based PhyloAssigner. For all processing strategies, default parameter settings and punctual variations were used. Our results revealed strong differences in total number of operational taxonomic units (OTUs), indicating that sequence preprocessing and clustering had a major impact on protist diversity estimates. However, diversity estimates of the abundant biosphere (abundance of ≥1%) were reproducible for all conducted processing pipeline versions. A qualitative comparison of diatom genera emphasized strong differences between the pipelines in which phylogenetic placement of sequences came closest to light microscopy‐based diatom identification. We conclude that diversity studies using different sequence processing strategies are comparable if the focus is on higher taxonomic levels, and if abundance thresholds are used to filter out OTUs of the rare biosphere.