A multi-amplicon 16S rRNA sequencing and analysis method for improved taxonomic profiling of bacterial communities.

Metagenomic sequencing of bacterial samples has become the gold standard for profiling microbial populations, but 16S rRNA profiling remains widely used due to advantages in sample throughput, cost, and sensitivity even though the approach is hampered by primer bias and lack of specificity. We hypothesized that a hybrid approach, that combined targeted PCR amplification with high-throughput sequencing of multiple regions of the genome, would capture many of the advantages of both approaches. We developed a method that identifies and quantifies members of bacterial communities through simultaneous analysis of multiple variable regions of the bacterial 16S rRNA gene. The method combines high-throughput microfluidics for PCR amplification, short read DNA sequencing, and a custom algorithm named MVRSION (Multiple 16S Variable Region Species-Level IdentificatiON) for optimizing taxonomic assignment. MVRSION performance was compared to single variable region analyses (V3 or V4) of five synthetic mixtures of human gut bacterial strains using existing software (QIIME), and the results of community profiling by shotgun sequencing (COPRO-Seq) of fecal DNA samples collected from gnotobiotic mice colonized with a defined, phylogenetically diverse consortium of human gut bacterial strains. Positive predictive values for MVRSION ranged from 65%-91% versus 44%-61% for single region QIIME analyses (p < .01, p < .001), while the abundance estimate r2 for MVRSION compared to COPRO-Seq was 0.77 vs. 0.46 and 0.45 for V3-QIIME and V4-QIIME, respectively. MVRSION represents a generally applicable tool for taxonomic classification that is superior to single-region 16S rRNA methods, resource efficient, highly scalable for assessing the microbial composition of up to thousands of samples concurrently, with multiple applications ranging from whole community profiling to targeted tracking of organisms of interest in diverse habitats as a function of specified variables/perturbations.

Schriefer AE ，Cliften PF ，Hibberd MC ，Sawyer C ，Brown-Kennerly V ，Burcea L ，Klotz E ，Crosby SD ，Gordon JI ，Head RD ... -  《-》

Combining 16S rRNA gene variable regions enables high-resolution microbial community profiling.

Fuks G ，Elgart M ，Amir A ，Zeisel A ，Turnbaugh PJ ，Soen Y ，Shental N ... -  《Microbiome》

Optimisation of methods for bacterial skin microbiome investigation: primer selection and comparison of the 454 versus MiSeq platform.

The composition of the skin microbiome is predicted to play a role in the development of conditions such as atopic eczema and psoriasis. 16S rRNA gene sequencing allows the investigation of bacterial microbiota. A significant challenge in this field is development of cost effective high throughput methodologies for the robust interrogation of the skin microbiota, where biomass is low. Here we describe validation of methodologies for 16S rRNA (ribosomal ribonucleic acid) gene sequencing from the skin microbiome, using the Illumina MiSeq platform, the selection of primer to amplify regions for sequencing and we compare results with the current standard protocols.. DNA was obtained from two low density mock communities of 11 diverse bacterial strains (with and without human DNA supplementation) and from swabs taken from the skin of healthy volunteers. This was amplified using primer pairs covering hypervariable regions of the 16S rRNA gene: primers 63F and 519R (V1-V3); and 347F and 803R (V3-V4). The resultant libraries were indexed for the MiSeq and Roche454 and sequenced. Both data sets were denoised, cleaned of chimeras and analysed using QIIME. There was no significant difference in the diversity indices at the phylum and the genus level observed between the platforms. The capture of diversity using the low density mock community samples demonstrated that the primer pair spanning the V3-V4 hypervariable region had better capture when compared to the primer pair for the V1-V3 region and was robust to spiking with human DNA. The pilot data generated using the V3-V4 region from the skin of healthy volunteers was consistent with these results, even at the genus level (Staphylococcus, Propionibacterium, Corynebacterium, Paracoccus, Micrococcus, Enhydrobacter and Deinococcus identified at similar abundances on both platforms). The results suggest that the bacterial community diversity captured using the V3-V4 16S rRNA hypervariable region from sequencing using the MiSeq platform is comparable to the Roche454 GS Junior platform. These findings provide evidence that the optimised method can be used in human clinical samples of low bacterial biomass such as the investigation of the skin microbiota.

Castelino M ，Eyre S ，Moat J ，Fox G ，Martin P ，Ho P ，Upton M ，Barton A ... -  《BMC MICROBIOLOGY》

Impact of DNA Sequencing and Analysis Methods on 16S rRNA Gene Bacterial Community Analysis of Dairy Products.

DNA sequencing and analysis methods were compared for 16S rRNA V4 PCR amplicon and genomic DNA (gDNA) mock communities encompassing nine bacterial species commonly found in milk and dairy products. The two communities comprised strain-specific DNA that was pooled before (gDNA) or after (PCR amplicon) the PCR step. The communities were sequenced on the Illumina MiSeq and Ion Torrent PGM platforms and then analyzed using the QIIME 1 (UCLUST) and Divisive Amplicon Denoising Algorithm 2 (DADA2) analysis pipelines with taxonomic comparisons to the Greengenes and Ribosomal Database Project (RDP) databases. Examination of the PCR amplicon mock community with these methods resulted in operational taxonomic units (OTUs) and amplicon sequence variants (ASVs) that ranged from 13 to 118 and were dependent on the DNA sequencing method and read assembly steps. The additional 4 to 109 OTUs/ASVs (from 9 OTUs/ASVs) included assignments to spurious taxa and sequence variants of the 9 species included in the mock community. Comparisons between the gDNA and PCR amplicon mock communities showed that combining gDNAs from the different strains prior to PCR resulted in up to 8.9-fold greater numbers of spurious OTUs/ASVs. However, the DNA sequencing method and paired-end read assembly steps conferred the largest effects on predictions of bacterial diversity, with effect sizes of 0.88 (Bray-Curtis) and 0.32 (weighted Unifrac), independent of the mock community type. Overall, DNA sequencing performed with the Ion Torrent PGM and analyzed with DADA2 and the Greengenes database resulted in the most accurate predictions of the mock community phylogeny, taxonomy, and diversity.IMPORTANCE Validated methods are urgently needed to improve DNA sequence-based assessments of complex bacterial communities. In this study, we used 16S rRNA PCR amplicon and gDNA mock community standards, consisting of nine, dairy-associated bacterial species, to evaluate the most commonly applied 16S rRNA marker gene DNA sequencing and analysis platforms used in evaluating dairy and other bacterial habitats. Our results show that bacterial metataxonomic assessments are largely dependent on the DNA sequencing platform and read curation method used. DADA2 improved sequence annotation compared with QIIME 1, and when combined with the Ion Torrent PGM DNA sequencing platform and the Greengenes database for taxonomic assignment, the most accurate representation of the dairy mock community standards was reached. This approach will be useful for validating sample collection and DNA extraction methods and ultimately investigating bacterial population dynamics in milk- and dairy-associated environments.

Xue Z ，Kable ME ，Marco ML 《mSphere》

Primer, Pipelines, Parameters: Issues in 16S rRNA Gene Sequencing.

Short-amplicon 16S rRNA gene sequencing is currently the method of choice for studies investigating microbiomes. However, comparative studies on differences in procedures are scarce. We sequenced human stool samples and mock communities with increasing complexity using a variety of commonly used protocols. Short amplicons targeting different variable regions (V-regions) or ranges thereof (V1-V2, V1-V3, V3-V4, V4, V4-V5, V6-V8, and V7-V9) were investigated for differences in the composition outcome due to primer choices. Next, the influence of clustering (operational taxonomic units [OTUs], zero-radius OTUs [zOTUs], and amplicon sequence variants [ASVs]), different databases (GreenGenes, the Ribosomal Database Project, Silva, the genomic-based 16S rRNA Database, and The All-Species Living Tree), and bioinformatic settings on taxonomic assignment were also investigated. We present a systematic comparison across all typically used V-regions using well-established primers. While it is known that the primer choice has a significant influence on the resulting microbial composition, we show that microbial profiles generated using different primer pairs need independent validation of performance. Further, comparing data sets across V-regions using different databases might be misleading due to differences in nomenclature (e.g., Enterorhabdus versus Adlercreutzia) and varying precisions in classification down to genus level. Overall, specific but important taxa are not picked up by certain primer pairs (e.g., Bacteroidetes is missed using primers 515F-944R) or due to the database used (e.g., Acetatifactor in GreenGenes and the genomic-based 16S rRNA Database). We found that appropriate truncation of amplicons is essential and different truncated-length combinations should be tested for each study. Finally, specific mock communities of sufficient and adequate complexity are highly recommended.IMPORTANCE In 16S rRNA gene sequencing, certain bacterial genera were found to be underrepresented or even missing in taxonomic profiles when using unsuitable primer combinations, outdated reference databases, or inadequate pipeline settings. Concerning the last, quality thresholds as well as bioinformatic settings (i.e., clustering approach, analysis pipeline, and specific adjustments such as truncation) are responsible for a number of observed differences between studies. Conclusions drawn by comparing one data set to another (e.g., between publications) appear to be problematic and require independent cross-validation using matching V-regions and uniform data processing. Therefore, we highlight the importance of a thought-out study design including sufficiently complex mock standards and appropriate V-region choice for the sample of interest. The use of processing pipelines and parameters must be tested beforehand.

Abellan-Schneyder I ，Matchado MS ，Reitmeier S ，Sommer A ，Sewald Z ，Baumbach J ，List M ，Neuhaus K ... -  《mSphere》