Next-generation Sequencing of 16S Ribosomal RNA Gene Amplicons

One of the major questions in microbial ecology is “who is there?” This question can be answered using various tools, but one of the long-lasting gold standards is to sequence 16S ribosomal RNA (rRNA) gene amplicons generated by domain-level PCR reactions amplifying from genomic DNA. Traditionally, this was performed by cloning and Sanger (capillary electrophoresis) sequencing of PCR amplicons. The advent of next-generation sequencing has tremendously simplified and increased the sequencing depth for 16S rRNA gene sequencing. The introduction of benchtop sequencers now allows small labs to perform their 16S rRNA sequencing in-house in a matter of days. Here, an approach for 16S rRNA gene amplicon sequencing using a benchtop next-generation sequencer is detailed. The environmental DNA is first amplified by PCR using primers that contain sequencing adapters and barcodes. They are then coupled to spherical particles via emulsion PCR. The particles are loaded on a disposable chip and the chip is inserted in the sequencing machine after which the sequencing is performed. The sequences are retrieved in fastq format, filtered and the barcodes are used to establish the sample membership of the reads. The filtered and binned reads are then further analyzed using publically available tools. An example analysis where the reads were classified with a taxonomy-finding algorithm within the software package Mothur is given. The method outlined here is simple, inexpensive and straightforward and should help smaller labs to take advantage from the ongoing genomic revolution.     

MLSA barcoding of Marichromatium spp. and reclassification of Marichromatium fluminis (Sucharita et al., 2010) as Phaeochromatium fluminis gen. nov. comb. nov

Thirty-one members of the genus Marichromatium were analysed based on multilocus sequence analysis (MLSA) of four concatenated protein-coding genes (fusA, pufM, dnaK, recA) along with the internal transcribed spacer (ITS; 16S-23S rRNA) region and 16S rRNA gene. The restriction patterns obtained from the in silico analysis of the concatenated sequences were good barcodes for the identification of Marichromatium spp. Distinct phenotypic, chemotaxonomic and molecular differences allowed the separation of Marichromatium fluminis JA418(T) into a new genus in the family Chromatiaceae, for which we propose the name Phaeochromatium fluminis gen. nov. comb. nov.     

Analysis of 23S rRNA genes in metagenomes - a case study from the Global Ocean Sampling Expedition

As an evolutionary marker, 23S ribosomal RNA (rRNA) offers more diagnostic sequence stretches and greater sequence variation than 16S rRNA. However, 23S rRNA is still not as widely used. Based on 80 metagenome samples from the Global Ocean Sampling (GOS) Expedition, the usefulness and taxonomic resolution of 23S rRNA were compared to those of 16S rRNA. Since 23S rRNA is approximately twice as large as 16S rRNA, twice as many 23S rRNA gene fragments were retrieved from the GOS reads than 16S rRNA gene fragments, with 23S rRNA gene fragments being generally about 100 bp longer. Datasets for 16S and 23S rRNA sequences revealed similar relative abundances for major marine bacterial and archaeal taxa. However, 16S rRNA sequences had a better taxonomic resolution due to their significantly larger reference database.Reevaluation of the specificity of previously published PCR amplification primers and group specific fluorescence in situ hybridization probes on this metagenomic set of non-amplified 23S rRNA sequences revealed that out of 16 primers investigated, only two had more than 90% target group coverage. Evaluations of two probes, BET42a and GAM42a, were in accordance with previous evaluations, with a discrepancy in the target group coverage of the GAM42a probe when evaluated against the GOS metagenomic dataset.     

Development and application of a real-time PCR approach for quantification of uncultured bacteria in the central Baltic Sea

We have developed a highly sensitive approach to assess the abundance of uncultured bacteria in water samples from the central Baltic Sea by using a noncultured member of the "Epsilonproteobacteria" related to Thiomicrospira denitrificans as an example. Environmental seawater samples and samples enriched for the target taxon provided a unique opportunity to test the approach over a broad range of abundances. The approach is based on a combination of taxon- and domain-specific real-time PCR measurements determining the relative T. denitrificans-like 16S rRNA gene and 16S rRNA abundances, as well as the determination of total cell counts and environmental RNA content. It allowed quantification of T. denitrificans-like 16S rRNA molecules or 16S rRNA genes as well as calculation of the number of ribosomes per T. denitrificans-like cell. Every real-time measurement and its specific primer system were calibrated using environmental nucleic acids obtained from the original habitat for external standardization. These standards, as well as the respective samples to be measured, were prepared from the same DNA or RNA extract. Enrichment samples could be analyzed directly, whereas environmental templates had to be preamplified with general bacterial primers before quantification. Preamplification increased the sensitivity of the assay by more than 4 orders of magnitude. Quantification of enrichments with or without a preamplification step yielded comparable results. T. denitrificans-like 16S rRNA molecules ranged from 7.1 x 10(3) to 4.4 x 10(9) copies ml(-1) or 0.002 to 49.7% relative abundance. T. denitrificans-like 16S rRNA genes ranged from 9.0 x 10(1) to 2.2 x10(6) copies ml(-1) or 0.01 to 49.7% relative abundance. Detection limits of this real-time-PCR approach were 20 16S rRNA molecules or 0.2 16S rRNA gene ml(-1). The number of ribosomes per T. denitrificans-like cell was estimated to range from 20 to 200 in seawater and reached up to 2,000 in the enrichments. The results indicate that our real-time PCR approach can be used to determine cellular and relative abundances of uncultured marine bacterial taxa and to provide information about their levels of activity in their natural environment.     

OTU Analysis Using Metagenomic Shotgun Sequencing Data

Because of technological limitations, the primer and amplification biases in targeted sequencing of 16S rRNA genes have veiled the true microbial diversity underlying environmental samples. However, the protocol of metagenomic shotgun sequencing provides 16S rRNA gene fragment data with natural immunity against the biases raised during priming and thus the potential of uncovering the true structure of microbial community by giving more accurate predictions of operational taxonomic units (OTUs). Nonetheless, the lack of statistically rigorous comparison between 16S rRNA gene fragments and other data types makes it difficult to interpret previously reported results using 16S rRNA gene fragments. Therefore, in the present work, we established a standard analysis pipeline that would help confirm if the differences in the data are true or are just due to potential technical bias. This pipeline is built by using simulated data to find optimal mapping and OTU prediction methods. The comparison between simulated datasets revealed a relationship between 16S rRNA gene fragments and full-length 16S rRNA sequences that a 16S rRNA gene fragment having a length >150 bp provides the same accuracy as a full-length 16S rRNA sequence using our proposed pipeline, which could serve as a good starting point for experimental design and making the comparison between 16S rRNA gene fragment-based and targeted 16S rRNA sequencing-based surveys possible.     

Quantification of Microthrix parvicella in activated sludge bacterial communities by real-time PCR

AIMS: This study was to develop a simple and reliable method for quantifying Microthrix parvicella 16S rRNA gene copies and its application to activated sludge samples collected from wastewater treatment plants (WWTP) with and without foaming problems. METHODS AND RESULTS: The relative frequency of M. parvicella was determined by combining real-time PCR assays for quantification of total bacterial 16S rRNA gene copies and M. parvicella 16S rRNA gene copies. The developed method was applied to analyse 32 activated sludge samples obtained from German WWTP. The level of M. parvicella 16S rRNA gene copies in the 18 nonfoaming samples was below 3% of the total number of 16S rRNA gene copies and in the range of 0-18% for the 14 foaming samples. CONCLUSIONS: The described method allows reliable monitoring of the amount of M. parvicella in activated sludge samples. SIGNIFICANCE AND IMPACT OF THE STUDY: The described method may become an important component of a warning system for forthcoming bulking and foaming episodes.     

Selective phylogenetic analysis targeted at 16S rRNA genes of thermophiles and hyperthermophiles in deep-subsurface geothermal environments

Deep-subsurface samples obtained by deep drilling are likely to be contaminated with mesophilic microorganisms in the drilling fluid, and this could affect determination of the community structure of the geothermal microflora using 16S rRNA gene clone library analysis. To eliminate possible contamination by PCR-amplified 16S rRNA genes from mesophiles, a combined thermal denaturation and enzyme digestion method, based on a strong correlation between the G+C content of the 16S rRNA gene and the optimum growth temperatures of most known prokaryotic cultures, was used prior to clone library construction. To validate this technique, hot spring fluid (76 degrees C) and river water (14 degrees C) were used to mimic a deep-subsurface sample contaminated with drilling fluid. After DNA extraction and PCR amplification of the 16S rRNA genes from individual samples separately, the amplified products from river water were observed to be denatured at 82 degrees C and completely digested by exonuclease I (Exo I), while the amplified products from hot spring fluid remained intact after denaturation at 84 degrees C and enzyme digestion with Exo I. DNAs extracted from the two samples were mixed and used as a template for amplification of the 16S rRNA genes. The amplified rRNA genes were denatured at 84 degrees C and digested with Exo I before clone library construction. The results indicated that the 16S rRNA gene sequences from the river water were almost completely eliminated, whereas those from the hot spring fluid remained.     

基于分子标记的宏基因组16S rRNA基因高变区选择策略

随着新一代DNA测序技术出现,人们能够同时对多个DNA样本的宏基因组进行并行分析,尤其是以16S rRNA基因高变区为分子标记的测序已经成为微生物多样性研究最为简洁有效的方法. 目前二代高通量测序的读长不能覆盖16S rRNA基因的全长,需要选择一个有效的高变区进行测序.十多年来,对于16S rRNA基因高变区的选择策略没有统一的标准.本文分析了常用的高变区选择策略,指出不同环境条件是影响高变区选择的重要因素之一.在此基础上,提出了高变区选择的参考准则,同时建议应对选择的高变区进行有效评估.     

Selective Phylogenetic Analysis Targeted at 16S rRNA Genes of Thermophiles and Hyperthermophiles in Deep-Subsurface Geothermal Environments

Deep-subsurface samples obtained by deep drilling are likely to be contaminated with mesophilic microorganisms in the drilling fluid, and this could affect determination of the community structure of the geothermal microflora using 16S rRNA gene clone library analysis. To eliminate possible contamination by PCR-amplified 16S rRNA genes from mesophiles, a combined thermal denaturation and enzyme digestion method, based on a strong correlation between the G+C content of the 16S rRNA gene and the optimum growth temperatures of most known prokaryotic cultures, was used prior to clone library construction. To validate this technique, hot spring fluid (76°C) and river water (14°C) were used to mimic a deep-subsurface sample contaminated with drilling fluid. After DNA extraction and PCR amplification of the 16S rRNA genes from individual samples separately, the amplified products from river water were observed to be denatured at 82°C and completely digested by exonuclease I (Exo I), while the amplified products from hot spring fluid remained intact after denaturation at 84°C and enzyme digestion with Exo I. DNAs extracted from the two samples were mixed and used as a template for amplification of the 16S rRNA genes. The amplified rRNA genes were denatured at 84°C and digested with Exo I before clone library construction. The results indicated that the 16S rRNA gene sequences from the river water were almost completely eliminated, whereas those from the hot spring fluid remained.     

Improved detection of bifidobacteria with optimised 16S rRNA-gene based pyrosequencing

The 16S rRNA gene is conserved across all bacteria and as such is routinely targeted in PCR surveys of bacterial diversity. PCR primer design aims to amplify as many different 16S rRNA gene sequences from as wide a range of organisms as possible, though there are no suitable 100% conserved regions of the gene, leading to bias. In the gastrointestinal tract, bifidobacteria are a key genus, but are often under-represented in 16S rRNA surveys of diversity. We have designed modified, 'bifidobacteria-optimised' universal primers, which we have demonstrated detection of bifidobacterial sequence present in DNA mixtures at 2% abundance, the lowest proportion tested. Optimisation did not compromise the detection of other organisms in infant faecal samples. Separate validation using fluorescence in situ hybridisation (FISH) shows that the proportions of bifidobacteria detected in faecal samples were in agreement with those obtained using 16S rRNA based pyrosequencing. For future studies looking at faecal microbiota, careful selection of primers will be key in order to ensure effective detection of bifidobacteria.     

PCR detection and 16S rRNA sequence-based phylogeny of a novel Propionibacterium acidipropionici applicable for enhanced fermentation of high moisture corn

AIMS: The aims of this study were to develop a sensitive and more rapid detection of Propionibacterium acidipropionici DH42 in silage and rumen fluid samples, and to explore its 16S rRNA sequence-based phylogeny. METHODS AND RESULTS: Nested polymerase chain reaction (PCR) was used with DH42-specific primers dhb1 and dhb2 for the secondary amplification of a 1267-bp fragment of 16S rRNA encoding gene. Using the established protocols for PCR amplification, as low as 10(2) and 10(3) CFU ml(-1) of strain DH42 in silage extracts and rumen fluid, respectively, were detected. To determine phylogenetic relationships between DH42 and other representatives of Propionibacterineae, a 1529-bp fragment of its 16S rRNA was amplified by PCR and sequenced. The propionibacterium DH42 formed a cluster with Eubacterium combesii, P. acidipropionici and P. microaerophilus. CONCLUSIONS: 16S rRNA-based PCR detection technique was developed for DH42 in silage and rumen fluid samples. The 16S rRNA sequence confirmed the earlier identification of strain DH42 as P. acidipropionici. However, variable nucleotide positions were revealed. SIGNIFICANCE AND IMPACT OF THE STUDY: Variability of 16S rRNA sequence within the species P. acidipropionici, determined in this study, poses the need of re-sequencing for some species of the suborder Propionibacterineae for a more reliable classification.     

Quantification of Hyphomicrobium populations in activated sludge from an industrial wastewater treatment system as determined by 16S rRNA analysis

The bacterial community structure of the activated sludge from a 25 million-gal-per-day industrial wastewater treatment plant was investigated using rRNA analysis. 16S ribosomal DNA (rDNA) libraries were created from three sludge samples taken on different dates. Partial rRNA gene sequences were obtained for 46 rDNA clones, and nearly complete 16S rRNA sequences were obtained for 18 clones. Seventeen of these clones were members of the beta subdivision, and their sequences showed high homology to sequences of known bacterial species as well as published 16S rDNA sequences from other activated sludge sources. Sixteen clones belonged to the alpha subdivision, 7 of which showed similarity to Hyphomicrobium species. This cluster was chosen for further studies due to earlier work on Hyphomicrobium sp. strain M3 isolated from this treatment plant. A nearly full-length 16S rDNA sequence was obtained from Hyphomicrobium sp. strain M3. Phylogenetic analysis revealed that Hyphomicrobium sp. strain M3 was 99% similar to Hyphomicrobium denitrificans DSM 1869(T) in Hyphomicrobium cluster II. Three of the cloned sequences from the activated sludge samples also grouped with those of Hyphomicrobium cluster II, with a 96% sequence similarity to that of Hyphomicrobium sp. strain M3. The other four cloned sequences from the activated sludge sample were more closely related to those of the Hyphomicrobium cluster I organisms (95 to 97% similarity). Whole-cell fluorescence hybridization of microorganisms in the activated sludge with genus-specific Hyphomicrobium probe S-G-Hypho-1241-a-A-19 enhanced the visualization of Hyphomicrobium and revealed that Hyphomicrobium appears to be abundant both on the outside of flocs and within the floc structure. Dot blot hybridization of activated sludge samples from 1995 with probes designed for Hyphomicrobium cluster I and Hyphomicrobium cluster II indicated that Hyphomicrobium cluster II-positive 16S rRNA dominated over Hyphomicrobium cluster I-positive 16S rRNA by 3- to 12-fold. Hyphomicrobium 16S rRNA comprised approximately 5% of the 16S rRNA in the activated sludge.     

Use of real-time PCR for group-specific quantification of aceticlastic methanogens in anaerobic processes: population dynamics and community structures

The TaqMan quantitative PCR (QPCR) method was used to detect and quantify the 16S rRNA genes of aceticlastic methanogens at different taxonomic levels. Three different sets of primers coupled with a TaqMan probe for QPCR assays to detect the 16S rRNA genes of the order Methanosarcinales, as well as the families Methanosarcinaceae and Methanosaetaceae, were separately used. Using these primer and probe sets, the 16S rRNA genes of aceticlastic methanogens in samples from various anaerobic processes (i.e., nine pure cultures, batch experiment, and three different continuous processes including a full-scale digester), were monitored and quantified by QPCR assays. A batch experiment cultivating a mixture of aceticlastic methanogens, was conducted to monitor their population dynamics. Using this group-specific quantification method, the dynamics of a competition between two aceticlastic populations, as modulated by the acetate concentration, could well be described. The target 16S rRNA genes in environmental samples, collected from three different anaerobic processes treating sludge, cheese whey, and synthetic wastewaters, were additionally quantified. The quantified 16S rRNA gene concentrations for all samples successfully represented the community structures of the target methanogens, which were correlated accurately with the operational parameters of the anaerobic processes. It was also successful to demonstrate probe nesting of aceticlastic methanogens at the levels of order and family.     

Power analysis for real-time PCR quantification of genes in activated sludge and analysis of the variability introduced by DNA extraction

The aims of this study were to determine the power of discrimination of the real-time PCR assay for monitoring fluctuations in microbial populations within activated sludge and to identify sample processing points where methodological changes are needed to minimize the variability in target quantification. DNA was extracted using a commercially available kit from mixed liquor samples taken from the aeration tank of four bench-scale activated-sludge reactors operating at 2-, 5-, 10-, and 20-day solid retention times, with mixed-liquor volatile suspended solid (MLVSS) values ranging from 260 to 2,610 mg/liter. Real-time PCR assays for bacterial and Nitrospira 16S rRNA genes were chosen because they represent, respectively, a highly abundant and a less-abundant bacterial target subject to clustering within the activated sludge matrix. The mean coefficient of variation in DNA yields (measured as microgram of DNA per milligram of MLVSS) in triplicate extractions of 12 different samples was 12.2%. Based on power analyses, the variability associated with DNA extraction had a small impact on the overall variability of the real-time PCR assay. Instead, a larger variability was associated with the PCR assay. The less-abundant target (Nitrospira 16S rRNA gene) had more variability than the highly abundant target (bacterial 16S rRNA gene), and samples from the lower-biomass reactors had more variability than samples from the higher-biomass reactors. Power analysis of real-time PCR assays indicated that three to five samples were necessary to detect a twofold increase in bacterial 16S rRNA genes, whereas three to five samples were required to detect a fivefold increase in Nitrospira 16S rRNA genes.     

Analysis of Bacterial and Fungal Communities in Japanese- and Chinese-Fermented Soybean Pastes Using Nested PCR–DGGE

The microbial diversity of Japanese- and Chinese-fermented soybean pastes was investigated using nested PCR–denaturing gradient gel electrophoresis (DGGE). Five Japanese-fermented soybean paste samples and three Chinese-fermented soybean paste samples were analyzed for bacteria and fungi. Extracted DNA was used as a template for PCR to amplify 16S rRNA and 18S rRNA genes. The nearly complete 16S rRNA and 18S rRNA genes were amplified using universal primers, and the resulting products were subsequently used as a template in a nested PCR to obtain suitable fragments for DGGE. Tetragenococcus halophilus and Staphylococcus gallinarum were found to dominate the bacterial microbiota in Japanese samples, whereas Bacillus sp. was detected as the predominant species in Chinese samples. DGGE analysis of fungi in soybean pastes determined the presence of Aspergillus oryzae and Zygosaccharomyces rouxii in most of the Chinese and Japanese samples. Some differences were observed in the bacterial diversity of Japanese- and Chinese-fermented soybean pastes.     

Confident phylogenetic identification of uncultured prokaryotes through long read amplicon sequencing of the 16S-ITS-23S rRNA operon

Amplicon sequencing of the 16S rRNA gene is the predominant method to quantify microbial compositions and to discover novel lineages. However, traditional short amplicons often do not contain enough information to confidently resolve their phylogeny. Here we present a cost-effective protocol that amplifies a large part of the rRNA operon and sequences the amplicons with PacBio technology. We tested our method on a mock community and developed a read-curation pipeline that reduces the overall read error rate to 0.18%. Applying our method on four environmental samples, we captured near full-length rRNA operon amplicons from a large diversity of prokaryotes. The method operated at moderately high-throughput (22286–37,850 raw ccs reads) and generated a large amount of putative novel archaeal 23S rRNA gene sequences compared to the archaeal SILVA database. These long amplicons allowed for higher resolution during taxonomic classification by means of long (∼1000 bp) 16S rRNA gene fragments and for substantially more confident phylogenies by means of combined near full-length 16S and 23S rRNA gene sequences, compared to shorter traditional amplicons (250 bp of the 16S rRNA gene). We recommend our method to those who wish to cost-effectively and confidently estimate the phylogenetic diversity of prokaryotes in environmental samples at high throughput.     

Bacterial diversity in the bacterioneuston (sea surface microlayer): the bacterioneuston through the looking glass

The bacterioneuston is defined as the community of bacteria present within the neuston or sea surface microlayer. Bacteria within this layer were sampled using a membrane filter technique and bacterial diversity was compared with that in the underlying pelagic coastal seawater using molecular ecological techniques. 16S rRNA gene libraries of approximately 500 clones were constructed from both bacterioneuston and the pelagic water samples and representative clones from each library were sequenced for comparison of bacterial diversity. The bacterioneuston was found to have a significantly lower bacterial diversity than the pelagic seawater, with only nine clone types (ecotaxa) as opposed to 46 ecotaxa in the pelagic seawater library. Surprisingly, the bacterioneuston clone library was dominated by 16S rRNA gene sequences affiliated to two groups of organisms, Vibrio spp. which accounted for over 68% of clones and Pseudoalteromonas spp. accounting for 21% of the library. The dominance of these two 16S rRNA gene sequence types within the bacterioneuston clone library was confirmed in a subsequent gene probing experiment. 16S rRNA gene probes specific for these groups of bacteria were designed and used to probe new libraries of 1000 clones from both the bacterioneuston and pelagic seawater DNA samples. This revealed that 57% of clones from the bacterioneuston library hybridized to a Vibrio sp.-specific 16S rRNA gene probe and 32% hybridized to a Pseudoalteromonas sp.-specific 16S rRNA gene probe. In contrast, the pelagic seawater library resulted in only 13% and 8% of 16S rRNA gene clones hybridizing to the Vibrio sp. and Pseudoalteromonas sp. probes respectively. Results from this study suggest that the bacterioneuston contains a distinct population of bacteria and warrants further detailed study at the molecular level.     

Profiling of the bacteria responsible for pyogenic liver abscess by 16S rRNA gene pyrosequencing

Pyogenic liver abscess (PLA) is a severe disease with considerable mortality and is often polymicrobial. Understanding the pathogens that cause PLA is the basis for PLA treatment. Here, we profiled the bacterial composition in PLA fluid by pyrosequencing the 16S ribosomal RNA (rRNA) gene based on next-generation sequencing (NGS) technology to identify etiological agents of PLA and to provide information of their 16S rRNA sequences for application to DNA-based techniques in the hospital. Twenty patients with PLA who underwent percutaneous catheter drainage, abscess culture, and blood culture for isolates were included. Genomic DNAs from abscess fluids were subjected to polymerase chain reaction and pyrosequencing of the 16S rRNA gene with a 454 GS Junior System. The abscess and blood cultures were positive in nine (45%) and four (20%) patients, respectively. Pyrosequencing of 16S rRNA gene showed that 90% of the PLA fluid samples contained single or multiple genera of known bacteria such as Klebsiella, Fusobacterium, Streptococcus, Bacteroides, Prevotella, Peptostreptococcus, unassigned Enterobacteriaceae, and Dialister. Klebsiella was predominantly found in the PLA fluid samples. All samples that carried unassigned bacteria had 26.8% reads on average. We demonstrated that the occurrence of PLA was associated with eight known bacterial genera as well as unassigned bacteria and that 16S rRNA gene sequencing was more useful than conventional culture methods for accurate identification of bacterial pathogens from PLA.     

Microbial Community Composition and Diversity via 16S rRNA Gene Amplicons: Evaluating the Illumina Platform

As new sequencing technologies become cheaper and older ones disappear, laboratories switch vendors and platforms. Validating the new setups is a crucial part of conducting rigorous scientific research. Here we report on the reliability and biases of performing bacterial 16S rRNA gene amplicon paired-end sequencing on the MiSeq Illumina platform. We designed a protocol using 50 barcode pairs to run samples in parallel and coded a pipeline to process the data. Sequencing the same sediment sample in 248 replicates as well as 70 samples from alkaline soda lakes, we evaluated the performance of the method with regards to estimates of alpha and beta diversity. Using different purification and DNA quantification procedures we always found up to 5-fold differences in the yield of sequences between individually barcodes samples. Using either a one-step or a two-step PCR preparation resulted in significantly different estimates in both alpha and beta diversity. Comparing with a previous method based on 454 pyrosequencing, we found that our Illumina protocol performed in a similar manner – with the exception for evenness estimates where correspondence between the methods was low. We further quantified the data loss at every processing step eventually accumulating to 50% of the raw reads. When evaluating different OTU clustering methods, we observed a stark contrast between the results of QIIME with default settings and the more recent UPARSE algorithm when it comes to the number of OTUs generated. Still, overall trends in alpha and beta diversity corresponded highly using both clustering methods. Our procedure performed well considering the precisions of alpha and beta diversity estimates, with insignificant effects of individual barcodes. Comparative analyses suggest that 454 and Illumina sequence data can be combined if the same PCR protocol and bioinformatic workflows are used for describing patterns in richness, beta-diversity and taxonomic composition.