Dataset Information

Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform.

ABSTRACT: With read lengths of currently up to 2 × 300 bp, high throughput and low sequencing costs Illumina's MiSeq is becoming one of the most utilized sequencing platforms worldwide. The platform is manageable and affordable even for smaller labs. This enables quick turnaround on a broad range of applications such as targeted gene sequencing, metagenomics, small genome sequencing and clinical molecular diagnostics. However, Illumina error profiles are still poorly understood and programs are therefore not designed for the idiosyncrasies of Illumina data. A better knowledge of the error patterns is essential for sequence analysis and vital if we are to draw valid conclusions. Studying true genetic variation in a population sample is fundamental for understanding diseases, evolution and origin. We conducted a large study on the error patterns for the MiSeq based on 16S rRNA amplicon sequencing data. We tested state-of-the-art library preparation methods for amplicon sequencing and showed that the library preparation method and the choice of primers are the most significant sources of bias and cause distinct error patterns. Furthermore we tested the efficiency of various error correction strategies and identified quality trimming (Sickle) combined with error correction (BayesHammer) followed by read overlapping (PANDAseq) as the most successful approach, reducing substitution error rates on average by 93%.

SUBMITTER: Schirmer M

PROVIDER: S-EPMC4381044 | biostudies-literature | 2015 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform.

Schirmer Melanie M Ijaz Umer Z UZ D'Amore Rosalinda R Hall Neil N Sloan William T WT Quince Christopher C

Nucleic acids research 20150113 6

With read lengths of currently up to 2 × 300 bp, high throughput and low sequencing costs Illumina's MiSeq is becoming one of the most utilized sequencing platforms worldwide. The platform is manageable and affordable even for smaller labs. This enables quick turnaround on a broad range of applications such as targeted gene sequencing, metagenomics, small genome sequencing and clinical molecular diagnostics. However, Illumina error profiles are still poorly understood and programs are therefore ...[more]

PMID: 25586220

Similar Datasets

Project description: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.

Project description:BackgroundIllumina sequencing platform requires base diversity in the initial 11 cycles for efficient cluster identification and colour matrix estimation. This limitation yields low-quality data for amplicon libraries having homogeneous base composition. Spike-in of PhiX library ensures base diversity but reduces the overall number of sequencing reads for data analysis. To overcome such low diversity issues during amplicon sequencing on illumina platforms, we developed a high throughput single amplicon sequencing method by introducing 'N' (0-10) spacers in target gene amplification primers that are pooled for simple handling.ResultWe evaluated the efficiency of 'N' (0-10) spacer-linked primers by targeting bacterial 16S V3-V4 region, demonstrating heterogeneous base library construction. The addition of 'N' (0-10) spacers causes sequencing frameshift at every base that leads to base diversity and produces heterogeneous high quality reads within a single amplicon library. We have written a python based command-line software,"MetReTrim", to trim the 'N' (0-10) spacers from the raw reads ( https://github.com/Mohak91/MetReTrim ). We further demonstrated the accuracy of this method by comparative mock community analysis with standard illumina V3-V4 primer method. The ZymoBIOMICS™ microbial community DNA standard was used as a control for this study. We performed data analysisusing the DADA2 pipeline where taxonomy was assigned using SILVA database as reference. We observed no difference between the communities represented by our method and standard illumina V3-V4 primer method.ConclusionThis method eliminates the need for PhiX spike-in for single amplicon sequencing on illumina MiSeq platform. This allows for sequencing of more number of samples in a run and a reduction in the overall cost. Given that Illumina sequencing works on SBS chemistry irrespective of the platform (such as HiSeq, MiSeq, NextSeq, NovaSeq, etc.) we propose that this strategy of using 'N' (0-10) spacer-linked primer design can be adopted for generating high-quality single locus amplicon sequencing in a high throughput manner across the illumina platform subject to further validation.

Dataset Information

Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform.

Publications

Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets