Abstract:Illumina's MiSeq has become the dominant platform for gene amplicon sequencing in microbial ecology studies; however, various technical concerns, such as reproducibility, still exist. To assess reproducibility, 16S rRNA gene amplicons from 18 soil samples of a reciprocal transplantation experiment were sequenced on an Illumina MiSeq. The V4 region of 16S rRNA gene from each sample was sequenced in triplicate with each replicate having a unique barcode. The average OTU overlap, without considering sequence abundance, at a rarefaction level of 10,323 sequences was 33.4±2.1% and 20.2±1.7% between two and among three technical replicates, respectively. When OTU sequence abundance was considered, the average sequence abundance weighted OTU overlap was 85.6±1.6% and 81.2±2.1% for two and three replicates, respectively. Removing singletons significantly increased the overlap for both (~1-3%, p<0.001). Increasing the sequencing depth to 160,000 reads by deep sequencing increased OTU overlap both when sequence abundance was considered (95%) and when not (44%). However, if singletons were not removed the overlap between two technical replicates (not considering sequence abundance) plateaus at 39% with 30,000 sequences. Diversity measures were not affected by the low overlap as α-diversities were similar among technical replicates while β-diversities (Bray-Curtis) were much smaller among technical replicates than among treatment replicates (e.g., 0.269 vs. 0.374). Higher diversity coverage, but lower OTU overlap, was observed when replicates were sequenced in separate runs. Detrended correspondence analysis indicated that while there was considerable variation among technical replicates, the reproducibility was sufficient for detecting treatment effects for the samples examined. These results suggest that although there is variation among technical replicates, amplicon sequencing on MiSeq is useful for analyzing microbial community structure if used appropriately and with caution. For example, including technical replicates, removing spurious sequences and unrepresentative OTUs, using a clustering method with a high stringency for OTU generation, estimating treatment effects at higher taxonomic levels, and adapting the unique molecular identifier (UMI) and other newly developed methods to lower PCR and sequencing error and to identify true low abundance rare species all can increase reproducibility.

Confidence Intervals of Similarity Values Determined for Cloned Ssu Rrna Genes from Environmental Samples

Cloning and Sequence Analysis of the ITS and 5.8S Rdna of Eustrongylides Spp.isolates from Monopterus Albus

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

A feasible roadmap to identifying significant intercellular genomic heterogeneity in deep sequencing data

Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution

Updating the 97% identity threshold for 16S ribosomal RNA OTUs

Information about variations in multiple copies of bacterial 16S rRNA genes may aid in species identification

Evaluation of the Reproducibility of Amplicon Sequencing with Illumina MiSeq Platform

High consistency between replicate 454 pyrosequencing analyses of ectomycorrhizal plant root samples

ARISA data from the human gut microbiome can detect individual differences observed by 454 sequencing regardless of binning strategy

Reproducibility and quantitation of amplicon sequencing-based detection

Estimating sequence similarity from read sets for clustering next-generation sequencing data

Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis

Impact of 16S rRNA Gene Sequence Analysis for Identification of Bacteria on Clinical Microbiology and Infectious Diseases

In Silico and Experimental Evaluation of Primer Sets for Species-Level Resolution of the Vaginal Microbiota Using 16S Ribosomal RNA Gene Sequencing

Modified RNA-seq Method for Microbial Community and Diversity Analysis Using Rrna in Different Types of Environmental Samples

The enigmatic fungal genus provides a theoretical framework for studying intragenomic variation in ribosomal DNA sequences

RNA-clique: a method for computing genetic distances from RNA-seq data

Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences

Evaluation of Pcr-Generated Chimeras: Mutations, and Heteroduplexes with 16s Rrna Gene-Based Cloning

Impact of similarity metrics on single-cell RNA-seq data clustering